https://mooseframework.inl.gov
Loading...
Searching...
No Matches
Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | Protected Member Functions | Protected Attributes | Static Protected Attributes | Private Member Functions | Private 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)
 
const Moose::Functor< Real > * defaultFunctor (const std::string &name)
 
const Moose::Functor< ADReal > * defaultFunctor (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.
 
const unsigned int _rz_radial_coord
 The radial coordinate index for RZ coordinate systems.
 
std::vector< const ADVariableValue * > _interior_vels
 The velocities at interior element quadrature points.
 
std::vector< const ADVariableGradient * > _interior_vel_grads
 The velocity gradients at interior element quadrature points.
 
std::vector< const Moose::Functor< ADReal > * > _face_vels
 The velocity functors used to evalute the velocity on element face quadrature points.
 

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< 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
 
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
 

Private Member Functions

Moose::MaterialDataType getMaterialDataType (const std::set< BoundaryID > &boundary_ids) const
 
unsigned int getMaxQps () const
 
void addConsumedPropertyName (const MooseObjectName &obj_name, const std::string &prop_name)
 
const Moose::Functor< T > * defaultFunctor (const std::string &name)
 
const Moose::Functor< Real > * defaultFunctor (const std::string &name)
 
const Moose::Functor< ADReal > * defaultFunctor (const std::string &name)
 
const Moose::Functor< T > & getFunctorByName (const std::string &name, SubProblem &subproblem, THREAD_ID tid)
 
virtual bool isADObject () const override
 

Private Attributes

const MooseObject & _mi_moose_object
 
const bool _mi_boundary_restricted
 
const std::set< SubdomainID > & _mi_block_ids
 
const std::set< BoundaryID > & _mi_boundary_ids
 
std::vector< std::unique_ptr< OptionalMaterialPropertyProxyBase< MaterialPropertyInterface > > > _optional_property_proxies
 
const InputParameters & _fi_params
 
const std::string _fi_name
 
SubProblem *const _fi_subproblem
 
const THREAD_ID _fi_tid
 
std::vector< std::unique_ptr< Moose::Functor< Real > > > _default_real_functors
 
std::vector< std::unique_ptr< Moose::Functor< ADReal > > > _default_ad_real_functors
 

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 32 of file MassContinuityAssemblyHelper.C.

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

Member Function Documentation

◆ lmFace()

void MassContinuityAssemblyHelper::lmFace ( )
overridevirtual

Implements ElementAndTraceScalarHDGAssemblyHelper.

Definition at line 94 of file MassContinuityAssemblyHelper.C.

95{
96 for (const auto qp : make_range(_qrule_face->n_points()))
97 {
98 ADRealVectorValue interior_vel, face_vel;
99 for (const auto d : index_range(_interior_vels))
100 {
101 interior_vel(d) = (*_interior_vels[d])[qp];
102 face_vel(d) = (*_face_vels[d])(
103 Moose::ElemSideQpArg{_current_elem, _current_side, qp, _qrule_face, _q_point_face[qp]},
104 _ti.determineState());
105 }
106 const auto qp_term = (interior_vel - face_vel) * _normals[qp] * _JxW_face[qp];
107 for (const auto i : index_range(_lm_re))
108 _lm_re(i) += _lm_phi_face[i][qp] * qp_term;
109 }
110}
ElementAndTraceScalarHDGAssemblyHelper::_q_point_face
const MooseArray< Point > & _q_point_face
ElementAndTraceScalarHDGAssemblyHelper::_ti
const TransientInterface & _ti
ElementAndTraceScalarHDGAssemblyHelper::_normals
const MooseArray< Point > & _normals
ElementAndTraceScalarHDGAssemblyHelper::_lm_phi_face
const MooseArray< std::vector< Real > > & _lm_phi_face
ElementAndTraceScalarHDGAssemblyHelper::_qrule_face
const libMesh::QBase *const & _qrule_face
ElementAndTraceScalarHDGAssemblyHelper::_JxW_face
const MooseArray< Real > & _JxW_face
TransientInterface::determineState
Moose::StateArg determineState() const
libMesh::index_range
auto index_range(const T &sizable)
libMesh::make_range
IntRange< T > make_range(T beg, T end)

◆ scalarDirichlet()

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

Implements ElementAndTraceScalarHDGAssemblyHelper.

Definition at line 42 of file MassContinuityAssemblyHelper.h.

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

◆ scalarFace()

void MassContinuityAssemblyHelper::scalarFace ( )
overridevirtual

Implements ElementAndTraceScalarHDGAssemblyHelper.

Definition at line 89 of file MassContinuityAssemblyHelper.C.

90{
91}

◆ scalarVolume()

void MassContinuityAssemblyHelper::scalarVolume ( )
overridevirtual

Implements ElementAndTraceScalarHDGAssemblyHelper.

Definition at line 71 of file MassContinuityAssemblyHelper.C.

72{
73 for (const auto qp : make_range(_qrule->n_points()))
74 {
75 ADReal divergence = 0;
76 for (const auto d : index_range(_interior_vel_grads))
77 {
78 divergence += (*_interior_vel_grads[d])[qp](d);
79 if (_coord_sys == Moose::COORD_RZ && (_rz_radial_coord == d))
80 divergence += (*_interior_vels[d])[qp] / _q_point[qp](_rz_radial_coord);
81 }
82 const auto qp_term = _JxW[qp] * divergence;
83 for (const auto i : index_range(_scalar_phi))
84 _scalar_re(i) -= qp_term * _scalar_phi[i][qp];
85 }
86}
ElementAndTraceScalarHDGAssemblyHelper::_JxW
const MooseArray< Real > & _JxW
ElementAndTraceScalarHDGAssemblyHelper::_qrule
const libMesh::QBase *const & _qrule
ElementAndTraceScalarHDGAssemblyHelper::_q_point
const MooseArray< Point > & _q_point
ElementAndTraceScalarHDGAssemblyHelper::_scalar_phi
const MooseArray< std::vector< Real > > & _scalar_phi
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

◆ validParams()

InputParameters MassContinuityAssemblyHelper::validParams ( )
static

Definition at line 22 of file MassContinuityAssemblyHelper.C.

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

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

Member Data Documentation

◆ _coord_sys

const Moose::CoordinateSystemType MassContinuityAssemblyHelper::_coord_sys

The coordinate system.

Definition at line 50 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 62 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 59 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 56 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 53 of file MassContinuityAssemblyHelper.h.

Referenced by scalarVolume().

The documentation for this class was generated from the following files:
Generated on Fri Aug 21 2026 13:43:14 for https://mooseframework.inl.gov by doxygen 1.9.8