https://mooseframework.inl.gov
Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | Protected Member Functions | Protected Attributes | Static Protected Attributes | List of all members
MassContinuityAssemblyHelper Class Reference

Implements all the methods for assembling a hybridized interior penalty discontinuous Galerkin (IPDG-H), which is a type of HDG method, discretization of the conservation of mass equation for incompressible Navier-Stokes. More...

#include <MassContinuityAssemblyHelper.h>

Inheritance diagram for MassContinuityAssemblyHelper:
[legend]

Public Member Functions

 MassContinuityAssemblyHelper (const MooseObject *const moose_obj, MooseVariableDependencyInterface *const mvdi, const TransientInterface *const ti, const MooseMesh &mesh, SystemBase &sys, const Assembly &assembly, const THREAD_ID tid, const std::set< SubdomainID > &block_ids, const std::set< BoundaryID > &boundary_ids)
 
virtual void scalarVolume () override
 
virtual void scalarFace () override
 
virtual void lmFace () override
 
virtual void scalarDirichlet (const Moose::Functor< Real > &) override
 
void resizeResiduals ()
 
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 ()
 
const MaterialProperty< T > & getFaceMaterialProperty (const std::string &name)
 
const MaterialProperty< T > & getFaceMaterialPropertyByName (const std::string &name)
 
const ADMaterialProperty< T > & getFaceADMaterialProperty (const std::string &name)
 
const MaterialProperty< T > & getFaceMaterialPropertyOld (const std::string &name)
 
const MaterialProperty< T > & getFaceMaterialPropertyOlder (const std::string &name)
 
const auto & getGenericFaceMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialProperty (const std::string &name, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialProperty (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialProperty (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialProperty (const std::string &name, const unsigned int state=0)
 
const ADMaterialProperty< T > & getADMaterialProperty (const std::string &name, MaterialData &material_data)
 
const ADMaterialProperty< T > & getADMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > & getADMaterialProperty (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOld (const std::string &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOld (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOld (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOlder (const std::string &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOlder (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOlder (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const MaterialPropertyName &name)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const MaterialPropertyName &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const MaterialPropertyName &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const MaterialPropertyName &name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyByName (const std::string &prop_name_in)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyOldByName (const std::string &prop_name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyOlderByName (const std::string &prop_name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialProperty (const std::string &name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyOld (const std::string &name)
 
Moose::Kokkos::MaterialProperty< T, dimension > getKokkosMaterialPropertyOlder (const std::string &name)
 
std::pair< const MaterialProperty< T > *, std::set< SubdomainID > > getBlockMaterialProperty (const MaterialPropertyName &name)
 
std::pair< Moose::Kokkos::MaterialProperty< T, dimension >, std::set< SubdomainID > > getKokkosBlockMaterialProperty (const MaterialPropertyName &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialProperty ()
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialPropertyByName (const std::string &prop_name)
 
const MaterialProperty< T > & getZeroMaterialProperty (Ts... args)
 
std::set< SubdomainID > getMaterialPropertyBlocks (const std::string &name)
 
std::vector< SubdomainName > getMaterialPropertyBlockNames (const std::string &name)
 
std::set< BoundaryID > getMaterialPropertyBoundaryIDs (const std::string &name)
 
std::vector< BoundaryName > getMaterialPropertyBoundaryNames (const std::string &name)
 
void checkBlockAndBoundaryCompatibility (std::shared_ptr< MaterialBase > discrete)
 
std::unordered_map< SubdomainID, std::vector< MaterialBase *> > buildRequiredMaterials (bool allow_stateful=true)
 
void statefulPropertiesAllowed (bool)
 
virtual bool getMaterialPropertyCalled () const
 
virtual const std::unordered_set< unsigned int > & getMatPropDependencies () const
 
virtual void resolveOptionalProperties ()
 
const GenericMaterialProperty< T, is_ad > & getPossiblyConstantGenericMaterialPropertyByName (const MaterialPropertyName &prop_name, MaterialData &material_data, const unsigned int state)
 
const GenericMaterialProperty< T, is_ad > & getGenericNeighborMaterialProperty (const std::string &name, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericNeighborMaterialProperty (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getNeighborMaterialProperty (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getNeighborMaterialProperty (const std::string &name, const unsigned int state=0)
 
const ADMaterialProperty< T > & getNeighborADMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > & getNeighborADMaterialProperty (const std::string &name)
 
const MaterialProperty< T > & getNeighborMaterialPropertyOld (const std::string &name)
 
const MaterialProperty< T > & getNeighborMaterialPropertyOld (const std::string &name)
 
const MaterialProperty< T > & getNeighborMaterialPropertyOlder (const std::string &name)
 
const MaterialProperty< T > & getNeighborMaterialPropertyOlder (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericNeighborMaterialPropertyByName (const std::string &name, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericNeighborMaterialPropertyByName (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getNeighborMaterialPropertyByName (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getNeighborMaterialPropertyByName (const std::string &name, const unsigned int state=0)
 
const ADMaterialProperty< T > & getNeighborADMaterialPropertyByName (const std::string &name)
 
const ADMaterialProperty< T > & getNeighborADMaterialPropertyByName (const std::string &name)
 
const GenericOptionalMaterialProperty< T, is_ad > & getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const GenericOptionalMaterialProperty< T, is_ad > & getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const OptionalMaterialProperty< T > & getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const OptionalMaterialProperty< T > & getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const OptionalADMaterialProperty< T > & getOptionalADMaterialProperty (const std::string &name)
 
const OptionalADMaterialProperty< T > & getOptionalADMaterialProperty (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOld (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOld (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOlder (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOlder (const std::string &name)
 
MaterialBase & getMaterial (const std::string &name)
 
MaterialBase & getMaterial (const std::string &name)
 
MaterialBase & getMaterialByName (const std::string &name, bool no_warn=false)
 
MaterialBase & getMaterialByName (const std::string &name, bool no_warn=false)
 
bool hasMaterialProperty (const std::string &name)
 
bool hasMaterialProperty (const std::string &name)
 
bool hasMaterialPropertyByName (const std::string &name)
 
bool hasMaterialPropertyByName (const std::string &name)
 
bool hasADMaterialProperty (const std::string &name)
 
bool hasADMaterialProperty (const std::string &name)
 
bool hasADMaterialPropertyByName (const std::string &name)
 
bool hasADMaterialPropertyByName (const std::string &name)
 
bool hasKokkosMaterialProperty (const std::string &name)
 
bool hasKokkosMaterialProperty (const std::string &name)
 
bool hasKokkosMaterialPropertyByName (const std::string &name)
 
bool hasKokkosMaterialPropertyByName (const std::string &name)
 
bool hasGenericMaterialProperty (const std::string &name)
 
bool hasGenericMaterialProperty (const std::string &name)
 
bool hasGenericMaterialPropertyByName (const std::string &name)
 
bool hasGenericMaterialPropertyByName (const std::string &name)
 

Static Public Member Functions

static InputParameters validParams ()
 
static std::string deduceFunctorName (const std::string &name, const InputParameters &params)
 

Public Attributes

const Moose::CoordinateSystemType _coord_sys
 The coordinate system. More...
 
const unsigned int _rz_radial_coord
 The radial coordinate index for RZ coordinate systems. More...
 
std::vector< const ADVariableValue * > _interior_vels
 The velocities at interior element quadrature points. More...
 
std::vector< const ADVariableGradient * > _interior_vel_grads
 The velocity gradients at interior element quadrature points. More...
 
std::vector< const Moose::Functor< ADReal > * > _face_vels
 The velocity functors used to evalute the velocity on element face quadrature points. More...
 

Static Public Attributes

static constexpr PropertyValue::id_type default_property_id
 
static constexpr PropertyValue::id_type zero_property_id
 

Protected Member Functions

virtual void checkMaterialProperty (const std::string &name, const unsigned int state)
 
virtual void getKokkosMaterialPropertyHook (const std::string &, const unsigned int)
 
void markMatPropRequested (const std::string &)
 
MaterialPropertyName getMaterialPropertyName (const std::string &name) const
 
void checkExecutionStage ()
 
const GenericMaterialProperty< T, is_ad > * defaultGenericMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > * defaultGenericMaterialProperty (const std::string &name)
 
const MaterialProperty< T > * defaultMaterialProperty (const std::string &name)
 
const MaterialProperty< T > * defaultMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > * defaultADMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > * defaultADMaterialProperty (const std::string &name)
 
std::string deduceFunctorName (const std::string &name) const
 
const Moose::Functor< T > & getFunctor (const std::string &name)
 
const Moose::Functor< T > & getFunctor (const std::string &name, THREAD_ID tid)
 
const Moose::Functor< T > & getFunctor (const std::string &name, SubProblem &subproblem)
 
const Moose::Functor< T > & getFunctor (const std::string &name, SubProblem &subproblem, THREAD_ID tid)
 
bool isFunctor (const std::string &name) const
 
bool isFunctor (const std::string &name, const SubProblem &subproblem) const
 
Moose::ElemArg makeElemArg (const Elem *elem, bool correct_skewnewss=false) const
 
void checkFunctorSupportsSideIntegration (const std::string &name, bool qp_integration)
 
const Moose::Functor< T > & getFunctorByName (const std::string &name)
 

Protected Attributes

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< ADRealVectorValue > & _grad_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< RealVectorValue > > & _grad_scalar_phi_face
 
const MooseArray< std::vector< Real > > & _lm_phi_face
 
const Real & _elem_volume
 
const Real & _side_area
 
const Elem *const & _ip_current_elem
 
const unsigned int & _ip_current_side
 
const Elem *const & _ip_current_side_elem
 
const MooseArray< Real > & _ip_JxW
 
const QBase *const & _ip_qrule
 
const MooseArray< Point > & _ip_q_point
 
const MooseArray< Real > & _ip_JxW_face
 
const QBase *const & _ip_qrule_face
 
const MooseArray< Point > & _ip_q_point_face
 
const MooseArray< Point > & _ip_normals
 
DenseVector< ADReal > _scalar_re
 
DenseVector< ADReal > _lm_re
 
MaterialData & _face_material_data
 
MaterialData & _neighbor_material_data
 
const InputParameters & _mi_params
 
const std::string _mi_name
 
const MooseObjectName _mi_moose_object_name
 
FEProblemBase & _mi_feproblem
 
SubProblem & _mi_subproblem
 
const THREAD_ID _mi_tid
 
const bool _is_kokkos_object
 
const Moose::MaterialDataType _material_data_type
 
MaterialData & _material_data
 
bool _stateful_allowed
 
bool _get_material_property_called
 
std::vector< std::unique_ptr< PropertyValue > > _default_properties
 
std::unordered_set< unsigned int > _material_property_dependencies
 
const MaterialPropertyName _get_suffix
 
const bool _use_interpolated_state
 

Static Protected Attributes

static const std::string _interpolated_old
 
static const std::string _interpolated_older
 

Detailed Description

Implements all the methods for assembling a hybridized interior penalty discontinuous Galerkin (IPDG-H), which is a type of HDG method, discretization of the conservation of mass equation for incompressible Navier-Stokes.

These routines may be called by both HDG kernels and integrated boundary conditions.

Definition at line 23 of file MassContinuityAssemblyHelper.h.

Constructor & Destructor Documentation

◆ MassContinuityAssemblyHelper()

MassContinuityAssemblyHelper::MassContinuityAssemblyHelper ( const MooseObject *const  moose_obj,
MooseVariableDependencyInterface *const  mvdi,
const TransientInterface *const  ti,
const MooseMesh &  mesh,
SystemBase &  sys,
const Assembly &  assembly,
const THREAD_ID  tid,
const std::set< SubdomainID > &  block_ids,
const std::set< BoundaryID > &  boundary_ids 
)

Definition at line 29 of file MassContinuityAssemblyHelper.C.

39  : IPHDGAssemblyHelper(moose_obj, mvdi, ti, sys, assembly, tid, block_ids, boundary_ids),
40  ADFunctorInterface(moose_obj),
41  _coord_sys(assembly.coordSystem()),
42  _rz_radial_coord(mesh.getAxisymmetricRadialCoord())
43 
44 {
45  for (const auto & interior_vel_var_name :
46  moose_obj->getParam<std::vector<VariableName>>("interior_velocity_vars"))
47  {
48  const auto & var = sys.getFieldVariable<Real>(tid, interior_vel_var_name);
49  _interior_vels.push_back(&var.adSln());
50  _interior_vel_grads.push_back(&var.adGradSln());
51  }
52  for (const auto & face_vel_functor_name :
53  moose_obj->getParam<std::vector<MooseFunctorName>>("face_velocity_functors"))
54  _face_vels.push_back(&getFunctorByName<ADReal>(face_vel_functor_name));
55 
56  if (mesh.dimension() != _interior_vels.size())
57  moose_obj->paramError(
58  "interior_velocity_vars",
59  "The number of interior velocity variables must be equal to the mesh dimension");
60  if (mesh.dimension() != _face_vels.size())
61  moose_obj->paramError(
62  "face_velocity_functors",
63  "The number of face velocity functors must be equal to the mesh dimension");
64 }
MooseObject::paramError
void paramError(const std::string &param, Args... args) const
MooseObject::getParam
const T & getParam(const std::string &name) const
MassContinuityAssemblyHelper::_rz_radial_coord
const unsigned int _rz_radial_coord
The radial coordinate index for RZ coordinate systems.
Definition: MassContinuityAssemblyHelper.h:52
MassContinuityAssemblyHelper::_face_vels
std::vector< const Moose::Functor< ADReal > * > _face_vels
The velocity functors used to evalute the velocity on element face quadrature points.
Definition: MassContinuityAssemblyHelper.h:61
mesh
MeshBase & mesh
MassContinuityAssemblyHelper::_interior_vel_grads
std::vector< const ADVariableGradient * > _interior_vel_grads
The velocity gradients at interior element quadrature points.
Definition: MassContinuityAssemblyHelper.h:58
Assembly::coordSystem
const Moose::CoordinateSystemType & coordSystem() const
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
MassContinuityAssemblyHelper::_coord_sys
const Moose::CoordinateSystemType _coord_sys
The coordinate system.
Definition: MassContinuityAssemblyHelper.h:49
SystemBase::getFieldVariable
MooseVariableFE< T > & getFieldVariable(THREAD_ID tid, const std::string &var_name)
ADFunctorInterface::ADFunctorInterface
ADFunctorInterface(const MooseObject *moose_object)
MassContinuityAssemblyHelper::_interior_vels
std::vector< const ADVariableValue * > _interior_vels
The velocities at interior element quadrature points.
Definition: MassContinuityAssemblyHelper.h:55
IPHDGAssemblyHelper::IPHDGAssemblyHelper
IPHDGAssemblyHelper(const MooseObject *const moose_obj, MooseVariableDependencyInterface *const mvdi, const TransientInterface *const ti, SystemBase &sys, const Assembly &assembly, const THREAD_ID tid, const std::set< SubdomainID > &blocks_ids, const std::set< BoundaryID > &boundary_ids)

Member Function Documentation

◆ lmFace()

void MassContinuityAssemblyHelper::lmFace ( )
overridevirtual

Implements IPHDGAssemblyHelper.

Definition at line 90 of file MassContinuityAssemblyHelper.C.

91 {
92  for (const auto qp : make_range(_ip_qrule_face->n_points()))
93  {
94  ADRealVectorValue interior_vel, face_vel;
95  for (const auto d : index_range(_interior_vels))
96  {
97  interior_vel(d) = (*_interior_vels[d])[qp];
98  face_vel(d) = (*_face_vels[d])(
99  Moose::ElemSideQpArg{
100  _ip_current_elem, _ip_current_side, qp, _ip_qrule_face, _ip_q_point_face[qp]},
101  _ti.determineState());
102  }
103  const auto qp_term = (interior_vel - face_vel) * _ip_normals[qp] * _ip_JxW_face[qp];
104  for (const auto i : index_range(_lm_re))
105  _lm_re(i) += _lm_phi_face[i][qp] * qp_term;
106  }
107 }
IPHDGAssemblyHelper::_lm_re
DenseVector< ADReal > _lm_re
IPHDGAssemblyHelper::_ip_current_side
const unsigned int & _ip_current_side
TransientInterface::determineState
Moose::StateArg determineState() const
MassContinuityAssemblyHelper::_face_vels
std::vector< const Moose::Functor< ADReal > * > _face_vels
The velocity functors used to evalute the velocity on element face quadrature points.
Definition: MassContinuityAssemblyHelper.h:61
IPHDGAssemblyHelper::_lm_phi_face
const MooseArray< std::vector< Real > > & _lm_phi_face
IPHDGAssemblyHelper::_ip_qrule_face
const QBase *const & _ip_qrule_face
IPHDGAssemblyHelper::_ti
const TransientInterface & _ti
IPHDGAssemblyHelper::_ip_current_elem
const Elem *const & _ip_current_elem
IPHDGAssemblyHelper::_ip_normals
const MooseArray< Point > & _ip_normals
IPHDGAssemblyHelper::_ip_q_point_face
const MooseArray< Point > & _ip_q_point_face
libMesh::QBase::n_points
unsigned int n_points() const
IPHDGAssemblyHelper::_ip_JxW_face
const MooseArray< Real > & _ip_JxW_face
make_range
IntRange< T > make_range(T beg, T end)
index_range
auto index_range(const T &sizable)
MassContinuityAssemblyHelper::_interior_vels
std::vector< const ADVariableValue * > _interior_vels
The velocities at interior element quadrature points.
Definition: MassContinuityAssemblyHelper.h:55

◆ scalarDirichlet()

virtual void MassContinuityAssemblyHelper::scalarDirichlet ( const Moose::Functor< Real > &  )
inlineoverridevirtual

Implements IPHDGAssemblyHelper.

Definition at line 41 of file MassContinuityAssemblyHelper.h.

42  {
43  mooseError(
44  "This kernel is applied for the pressure variable, but the pressure does not itself appear "
45  "in the weak form this kernel implements, so the scalarDirichlet method does not fit");
46  }
mooseError
void mooseError(Args &&... args)

◆ scalarFace()

void MassContinuityAssemblyHelper::scalarFace ( )
overridevirtual

Implements IPHDGAssemblyHelper.

Definition at line 85 of file MassContinuityAssemblyHelper.C.

86 {
87 }

◆ scalarVolume()

void MassContinuityAssemblyHelper::scalarVolume ( )
overridevirtual

Implements IPHDGAssemblyHelper.

Definition at line 67 of file MassContinuityAssemblyHelper.C.

68 {
69  for (const auto qp : make_range(_ip_qrule->n_points()))
70  {
71  ADReal divergence = 0;
72  for (const auto d : index_range(_interior_vel_grads))
73  {
74  divergence += (*_interior_vel_grads[d])[qp](d);
75  if (_coord_sys == Moose::COORD_RZ && (_rz_radial_coord == d))
76  divergence += (*_interior_vels[d])[qp] / _ip_q_point[qp](_rz_radial_coord);
77  }
78  const auto qp_term = _ip_JxW[qp] * divergence;
79  for (const auto i : index_range(_scalar_phi))
80  _scalar_re(i) -= qp_term * _scalar_phi[i][qp];
81  }
82 }
IPHDGAssemblyHelper::_scalar_re
DenseVector< ADReal > _scalar_re
NS::divergence
T divergence(const TensorValue< T > &gradient, const VectorType &value, const PointType &point, const Moose::CoordinateSystemType &coord_sys, const unsigned int rz_radial_coord)
Compute the divergence of a vector given its matrix of derivatives.
Definition: NavierStokesMethods.h:172
MassContinuityAssemblyHelper::_rz_radial_coord
const unsigned int _rz_radial_coord
The radial coordinate index for RZ coordinate systems.
Definition: MassContinuityAssemblyHelper.h:52
IPHDGAssemblyHelper::_ip_JxW
const MooseArray< Real > & _ip_JxW
IPHDGAssemblyHelper::_ip_qrule
const QBase *const & _ip_qrule
libMesh::QBase::n_points
unsigned int n_points() const
MassContinuityAssemblyHelper::_interior_vel_grads
std::vector< const ADVariableGradient * > _interior_vel_grads
The velocity gradients at interior element quadrature points.
Definition: MassContinuityAssemblyHelper.h:58
MassContinuityAssemblyHelper::_coord_sys
const Moose::CoordinateSystemType _coord_sys
The coordinate system.
Definition: MassContinuityAssemblyHelper.h:49
make_range
IntRange< T > make_range(T beg, T end)
IPHDGAssemblyHelper::_ip_q_point
const MooseArray< Point > & _ip_q_point
index_range
auto index_range(const T &sizable)
MassContinuityAssemblyHelper::_interior_vels
std::vector< const ADVariableValue * > _interior_vels
The velocities at interior element quadrature points.
Definition: MassContinuityAssemblyHelper.h:55
IPHDGAssemblyHelper::_scalar_phi
const MooseArray< std::vector< Real > > & _scalar_phi

◆ validParams()

InputParameters MassContinuityAssemblyHelper::validParams ( )
static

Definition at line 19 of file MassContinuityAssemblyHelper.C.

Referenced by MassContinuityIPHDGBC::validParams(), and MassContinuityIPHDGKernel::validParams().

20 {
21  auto params = IPHDGAssemblyHelper::validParams();
22  params.addRequiredCoupledVar("interior_velocity_vars",
23  "The velocity variables on the element interiors");
24  params.addRequiredParam<std::vector<MooseFunctorName>>(
25  "face_velocity_functors", "The velocity variables on element faces/facets");
26  return params;
27 }
IPHDGAssemblyHelper::validParams
static InputParameters validParams()

Member Data Documentation

◆ _coord_sys

const Moose::CoordinateSystemType MassContinuityAssemblyHelper::_coord_sys

The coordinate system.

Definition at line 49 of file MassContinuityAssemblyHelper.h.

Referenced by scalarVolume().

◆ _face_vels

std::vector<const Moose::Functor<ADReal> *> MassContinuityAssemblyHelper::_face_vels

The velocity functors used to evalute the velocity on element face quadrature points.

Definition at line 61 of file MassContinuityAssemblyHelper.h.

Referenced by lmFace(), and MassContinuityAssemblyHelper().

◆ _interior_vel_grads

std::vector<const ADVariableGradient *> MassContinuityAssemblyHelper::_interior_vel_grads

The velocity gradients at interior element quadrature points.

Definition at line 58 of file MassContinuityAssemblyHelper.h.

Referenced by MassContinuityAssemblyHelper(), and scalarVolume().

◆ _interior_vels

std::vector<const ADVariableValue *> MassContinuityAssemblyHelper::_interior_vels

The velocities at interior element quadrature points.

Definition at line 55 of file MassContinuityAssemblyHelper.h.

Referenced by lmFace(), MassContinuityAssemblyHelper(), and scalarVolume().

◆ _rz_radial_coord

const unsigned int MassContinuityAssemblyHelper::_rz_radial_coord

The radial coordinate index for RZ coordinate systems.

Definition at line 52 of file MassContinuityAssemblyHelper.h.

Referenced by scalarVolume().

The documentation for this class was generated from the following files:
Generated on Tue Jun 2 2026 21:52:50 for https://mooseframework.inl.gov by   doxygen 1.8.14