ElementAndTraceScalarHDGAssemblyHelper.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#pragma once

#include "MooseArray.h"
#include "MooseFunctorForward.h"
#include "MooseTypes.h"
#include "ThreeMaterialPropertyInterface.h"

#include <array>
#include <set>
#include <string>
#include <vector>

template <typename>
class MooseVariableFE;
class Assembly;
class MooseObject;
class MooseVariableDependencyInterface;
class SystemBase;
class TransientInterface;
struct ADResidualsPacket;

class ElementAndTraceScalarHDGAssemblyHelper : public ThreeMaterialPropertyInterface
{
public:
  static InputParameters validParams();

  ElementAndTraceScalarHDGAssemblyHelper(const MooseObject * moose_obj,
                                         MooseVariableDependencyInterface * mvdi,
                                         const TransientInterface * ti,
                                         SystemBase & sys,
                                         const Assembly & assembly,
                                         THREAD_ID tid,
                                         const std::set<SubdomainID> & blocks_ids,
                                         const std::set<BoundaryID> & boundary_ids);

  void resizeResiduals();

  virtual void scalarVolume() = 0;

  virtual void scalarFace() = 0;

  virtual void lmFace() = 0;

  virtual void scalarDirichlet(const Moose::Functor<Real> & dirichlet_value) = 0;

  void lmDirichlet(const Moose::Functor<Real> & dirichlet_value);

  void lmPrescribedFlux(const Moose::Functor<Real> & flux_value);

  std::array<ADResidualsPacket, 2> taggingData() const;

  std::set<std::string> additionalROVariables();

  virtual ~ElementAndTraceScalarHDGAssemblyHelper() = default;

protected:
  const TransientInterface & _ti;

  const MooseVariableFE<Real> & _u_var;

  const MooseVariableFE<Real> & _u_face_var;

  const std::vector<dof_id_type> & _u_dof_indices;

  const std::vector<dof_id_type> & _lm_u_dof_indices;

  const MooseArray<ADReal> & _u_sol;

  const MooseArray<ADReal> & _lm_u_sol;

  const MooseArray<std::vector<Real>> & _scalar_phi;

  const MooseArray<std::vector<RealVectorValue>> & _grad_scalar_phi;

  const MooseArray<std::vector<Real>> & _scalar_phi_face;

  const MooseArray<std::vector<Real>> & _lm_phi_face;

  const Real & _elem_volume;

  const Real & _side_area;

  const Elem * const & _current_elem;

  const unsigned int & _current_side;

  const Elem * const & _current_side_elem;

  const MooseArray<Real> & _JxW;

  const libMesh::QBase * const & _qrule;

  const MooseArray<Point> & _q_point;

  const MooseArray<Real> & _JxW_face;

  const libMesh::QBase * const & _qrule_face;

  const MooseArray<Point> & _q_point_face;

  const MooseArray<Point> & _normals;

  DenseVector<ADReal> _scalar_re;

  DenseVector<ADReal> _lm_re;
};

inline void
ElementAndTraceScalarHDGAssemblyHelper::resizeResiduals()
{
  _scalar_re.resize(_u_dof_indices.size());
  _lm_re.resize(_lm_u_dof_indices.size());
}