https://mooseframework.inl.gov
Public Member Functions | Static Public Member Functions | Static Public Attributes | Protected Member Functions | Protected Attributes | Static Protected Attributes | List of all members
NavierStokesStressIPHDGAssemblyHelper 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 Navier-Stokes stress term. More...

#include <NavierStokesStressIPHDGAssemblyHelper.h>

Inheritance diagram for NavierStokesStressIPHDGAssemblyHelper:
[legend]

Public Member Functions

 NavierStokesStressIPHDGAssemblyHelper (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)
 
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)
 
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 Public Attributes

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

Protected Member Functions

virtual void scalarVolume () override
 Computes a local residual vector for the weak form: (Dq, grad(w)) - (f, w) where D is the diffusivity, w are the test functions associated with the scalar field, and f is a forcing function. More...
 
virtual void scalarFace () override
 Computes a local residual vector for the weak form: -<Dq*n, w> + < * (u - {u}) * n * n, w> More...
 
virtual void scalarDirichlet (const Moose::Functor< Real > &dirichlet_value) override
 
virtual void lmFace () override
 
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)
 

Protected Attributes

const MooseVariableFE< Real > & _pressure_var
 The pressure variable on element interiors. More...
 
const MooseVariableFE< Real > & _pressure_face_var
 The pressure variable on element faces. More...
 
const ADVariableValue & _pressure_sol
 The pressure solution on element interiors. More...
 
const ADVariableValue & _pressure_face_sol
 The pressure solution on element faces. More...
 
const Moose::CoordinateSystemType & _coord_sys
 The coordinate system. More...
 
const unsigned int _rz_radial_coord
 The radial coordinate index for RZ coordinate systems. More...
 
unsigned int _component
 The velocity component this object is adding a residual for. More...
 
const ADMaterialProperty< Real > & _diff
 
const ADMaterialProperty< Real > & _face_diff
 
const Real _alpha
 
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 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 Navier-Stokes stress term.

Definition at line 18 of file NavierStokesStressIPHDGAssemblyHelper.h.

Constructor & Destructor Documentation

◆ NavierStokesStressIPHDGAssemblyHelper()

NavierStokesStressIPHDGAssemblyHelper::NavierStokesStressIPHDGAssemblyHelper ( 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 35 of file NavierStokesStressIPHDGAssemblyHelper.C.

45  : DiffusionIPHDGAssemblyHelper(moose_obj, mvdi, ti, sys, assembly, tid, block_ids, boundary_ids),
46  _pressure_var(sys.getFieldVariable<Real>(
47  tid, moose_obj->getParam<NonlinearVariableName>("pressure_variable"))),
48  _pressure_face_var(sys.getFieldVariable<Real>(
49  tid, moose_obj->getParam<NonlinearVariableName>("pressure_face_variable"))),
50  _pressure_sol(_pressure_var.adSln()),
51  _pressure_face_sol(_pressure_face_var.adSln()),
52  _coord_sys(assembly.coordSystem()),
53  _rz_radial_coord(mesh.getAxisymmetricRadialCoord()),
54  _component(moose_obj->getParam<unsigned int>("component"))
55 {
56 }
MooseObject::getParam
const T & getParam(const std::string &name) const
mesh
MeshBase & mesh
MooseVariableFE< Real >::adSln
const ADTemplateVariableValue< Real > & adSln() const override
NavierStokesStressIPHDGAssemblyHelper::_pressure_face_sol
const ADVariableValue & _pressure_face_sol
The pressure solution on element faces.
Definition: NavierStokesStressIPHDGAssemblyHelper.h:58
Assembly::coordSystem
const Moose::CoordinateSystemType & coordSystem() const
NavierStokesStressIPHDGAssemblyHelper::_rz_radial_coord
const unsigned int _rz_radial_coord
The radial coordinate index for RZ coordinate systems.
Definition: NavierStokesStressIPHDGAssemblyHelper.h:64
NavierStokesStressIPHDGAssemblyHelper::_pressure_var
const MooseVariableFE< Real > & _pressure_var
The pressure variable on element interiors.
Definition: NavierStokesStressIPHDGAssemblyHelper.h:51
NavierStokesStressIPHDGAssemblyHelper::_coord_sys
const Moose::CoordinateSystemType & _coord_sys
The coordinate system.
Definition: NavierStokesStressIPHDGAssemblyHelper.h:61
SystemBase::getFieldVariable
MooseVariableFE< T > & getFieldVariable(THREAD_ID tid, const std::string &var_name)
DiffusionIPHDGAssemblyHelper::DiffusionIPHDGAssemblyHelper
DiffusionIPHDGAssemblyHelper(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 > &block_ids, const std::set< BoundaryID > &boundary_ids)
NavierStokesStressIPHDGAssemblyHelper::_pressure_face_var
const MooseVariableFE< Real > & _pressure_face_var
The pressure variable on element faces.
Definition: NavierStokesStressIPHDGAssemblyHelper.h:53
NavierStokesStressIPHDGAssemblyHelper::_component
unsigned int _component
The velocity component this object is adding a residual for.
Definition: NavierStokesStressIPHDGAssemblyHelper.h:67
NavierStokesStressIPHDGAssemblyHelper::_pressure_sol
const ADVariableValue & _pressure_sol
The pressure solution on element interiors.
Definition: NavierStokesStressIPHDGAssemblyHelper.h:56

Member Function Documentation

◆ scalarDirichlet()

void NavierStokesStressIPHDGAssemblyHelper::scalarDirichlet ( const Moose::Functor< Real > &  dirichlet_value)
overrideprotectedvirtual

Reimplemented from DiffusionIPHDGAssemblyHelper.

Definition at line 84 of file NavierStokesStressIPHDGAssemblyHelper.C.

85 {
86  DiffusionIPHDGAssemblyHelper::scalarDirichlet(dirichlet_value);
87 
88  for (const auto i : index_range(_scalar_re))
89  for (const auto qp : make_range(_ip_qrule_face->n_points()))
90  _scalar_re(i) += _ip_JxW_face[qp] * _pressure_face_sol[qp] * _ip_normals[qp](_component) *
91  _scalar_phi_face[i][qp];
92 }
DiffusionIPHDGAssemblyHelper::_scalar_re
DenseVector< ADReal > _scalar_re
DiffusionIPHDGAssemblyHelper::_ip_qrule_face
const QBase *const & _ip_qrule_face
DiffusionIPHDGAssemblyHelper::_ip_normals
const MooseArray< Point > & _ip_normals
libMesh::QBase::n_points
unsigned int n_points() const
NavierStokesStressIPHDGAssemblyHelper::_pressure_face_sol
const ADVariableValue & _pressure_face_sol
The pressure solution on element faces.
Definition: NavierStokesStressIPHDGAssemblyHelper.h:58
DiffusionIPHDGAssemblyHelper::_ip_JxW_face
const MooseArray< Real > & _ip_JxW_face
DiffusionIPHDGAssemblyHelper::_scalar_phi_face
const MooseArray< std::vector< Real > > & _scalar_phi_face
make_range
IntRange< T > make_range(T beg, T end)
DiffusionIPHDGAssemblyHelper::scalarDirichlet
virtual void scalarDirichlet(const Moose::Functor< Real > &dirichlet_value) override
index_range
auto index_range(const T &sizable)
NavierStokesStressIPHDGAssemblyHelper::_component
unsigned int _component
The velocity component this object is adding a residual for.
Definition: NavierStokesStressIPHDGAssemblyHelper.h:67

◆ scalarFace()

void NavierStokesStressIPHDGAssemblyHelper::scalarFace ( )
overrideprotectedvirtual

Computes a local residual vector for the weak form: -<Dq*n, w> + < * (u - {u}) * n * n, w>

Reimplemented from DiffusionIPHDGAssemblyHelper.

Definition at line 73 of file NavierStokesStressIPHDGAssemblyHelper.C.

74 {
75  DiffusionIPHDGAssemblyHelper::scalarFace();
76 
77  for (const auto i : index_range(_scalar_re))
78  for (const auto qp : make_range(_ip_qrule_face->n_points()))
79  _scalar_re(i) += _ip_JxW_face[qp] * _pressure_face_sol[qp] * _ip_normals[qp](_component) *
80  _scalar_phi_face[i][qp];
81 }
DiffusionIPHDGAssemblyHelper::_scalar_re
DenseVector< ADReal > _scalar_re
DiffusionIPHDGAssemblyHelper::_ip_qrule_face
const QBase *const & _ip_qrule_face
DiffusionIPHDGAssemblyHelper::_ip_normals
const MooseArray< Point > & _ip_normals
libMesh::QBase::n_points
unsigned int n_points() const
NavierStokesStressIPHDGAssemblyHelper::_pressure_face_sol
const ADVariableValue & _pressure_face_sol
The pressure solution on element faces.
Definition: NavierStokesStressIPHDGAssemblyHelper.h:58
DiffusionIPHDGAssemblyHelper::_ip_JxW_face
const MooseArray< Real > & _ip_JxW_face
DiffusionIPHDGAssemblyHelper::_scalar_phi_face
const MooseArray< std::vector< Real > > & _scalar_phi_face
make_range
IntRange< T > make_range(T beg, T end)
DiffusionIPHDGAssemblyHelper::scalarFace
virtual void scalarFace() override
index_range
auto index_range(const T &sizable)
NavierStokesStressIPHDGAssemblyHelper::_component
unsigned int _component
The velocity component this object is adding a residual for.
Definition: NavierStokesStressIPHDGAssemblyHelper.h:67

◆ scalarVolume()

void NavierStokesStressIPHDGAssemblyHelper::scalarVolume ( )
overrideprotectedvirtual

Computes a local residual vector for the weak form: (Dq, grad(w)) - (f, w) where D is the diffusivity, w are the test functions associated with the scalar field, and f is a forcing function.

Reimplemented from DiffusionIPHDGAssemblyHelper.

Definition at line 59 of file NavierStokesStressIPHDGAssemblyHelper.C.

60 {
61  DiffusionIPHDGAssemblyHelper::scalarVolume();
62 
63  for (const auto qp : make_range(_ip_qrule->n_points()))
64  for (const auto i : index_range(_scalar_re))
65  {
66  _scalar_re(i) -= _ip_JxW[qp] * (_grad_scalar_phi[i][qp](_component) * _pressure_sol[qp]);
67  if (_coord_sys == Moose::COORD_RZ && (_rz_radial_coord == _component))
68  _scalar_re(i) -= _pressure_sol[qp] / _ip_q_point[qp](_rz_radial_coord) * _scalar_phi[i][qp];
69  }
70 }
DiffusionIPHDGAssemblyHelper::_scalar_re
DenseVector< ADReal > _scalar_re
DiffusionIPHDGAssemblyHelper::scalarVolume
virtual void scalarVolume() override
DiffusionIPHDGAssemblyHelper::_grad_scalar_phi
const MooseArray< std::vector< RealVectorValue > > & _grad_scalar_phi
DiffusionIPHDGAssemblyHelper::_ip_JxW
const MooseArray< Real > & _ip_JxW
DiffusionIPHDGAssemblyHelper::_ip_qrule
const QBase *const & _ip_qrule
libMesh::QBase::n_points
unsigned int n_points() const
NavierStokesStressIPHDGAssemblyHelper::_rz_radial_coord
const unsigned int _rz_radial_coord
The radial coordinate index for RZ coordinate systems.
Definition: NavierStokesStressIPHDGAssemblyHelper.h:64
NavierStokesStressIPHDGAssemblyHelper::_coord_sys
const Moose::CoordinateSystemType & _coord_sys
The coordinate system.
Definition: NavierStokesStressIPHDGAssemblyHelper.h:61
make_range
IntRange< T > make_range(T beg, T end)
DiffusionIPHDGAssemblyHelper::_ip_q_point
const MooseArray< Point > & _ip_q_point
index_range
auto index_range(const T &sizable)
DiffusionIPHDGAssemblyHelper::_scalar_phi
const MooseArray< std::vector< Real > > & _scalar_phi
NavierStokesStressIPHDGAssemblyHelper::_component
unsigned int _component
The velocity component this object is adding a residual for.
Definition: NavierStokesStressIPHDGAssemblyHelper.h:67
NavierStokesStressIPHDGAssemblyHelper::_pressure_sol
const ADVariableValue & _pressure_sol
The pressure solution on element interiors.
Definition: NavierStokesStressIPHDGAssemblyHelper.h:56

◆ validParams()

InputParameters NavierStokesStressIPHDGAssemblyHelper::validParams ( )
static

Definition at line 24 of file NavierStokesStressIPHDGAssemblyHelper.C.

Referenced by NavierStokesStressIPHDGDirichletBC::validParams(), and NavierStokesStressIPHDGKernel::validParams().

25 {
26  auto params = DiffusionIPHDGAssemblyHelper::validParams();
27  params.addRequiredParam<NonlinearVariableName>(
28  "pressure_variable", "The pressure variable that lives on element interiors.");
29  params.addRequiredParam<NonlinearVariableName>(
30  "pressure_face_variable", "The pressure variable that lives on element faces.");
31  params.addRequiredParam<unsigned int>("component", "number of component (0 = x, 1 = y, 2 = z)");
32  return params;
33 }
DiffusionIPHDGAssemblyHelper::validParams
static InputParameters validParams()

Member Data Documentation

◆ _component

unsigned int NavierStokesStressIPHDGAssemblyHelper::_component
protected

The velocity component this object is adding a residual for.

Definition at line 67 of file NavierStokesStressIPHDGAssemblyHelper.h.

Referenced by scalarDirichlet(), scalarFace(), and scalarVolume().

◆ _coord_sys

const Moose::CoordinateSystemType& NavierStokesStressIPHDGAssemblyHelper::_coord_sys
protected

The coordinate system.

Definition at line 61 of file NavierStokesStressIPHDGAssemblyHelper.h.

Referenced by scalarVolume().

◆ _pressure_face_sol

const ADVariableValue& NavierStokesStressIPHDGAssemblyHelper::_pressure_face_sol
protected

The pressure solution on element faces.

Definition at line 58 of file NavierStokesStressIPHDGAssemblyHelper.h.

Referenced by scalarDirichlet(), and scalarFace().

◆ _pressure_face_var

const MooseVariableFE<Real>& NavierStokesStressIPHDGAssemblyHelper::_pressure_face_var
protected

The pressure variable on element faces.

Definition at line 53 of file NavierStokesStressIPHDGAssemblyHelper.h.

◆ _pressure_sol

const ADVariableValue& NavierStokesStressIPHDGAssemblyHelper::_pressure_sol
protected

The pressure solution on element interiors.

Definition at line 56 of file NavierStokesStressIPHDGAssemblyHelper.h.

Referenced by scalarVolume().

◆ _pressure_var

const MooseVariableFE<Real>& NavierStokesStressIPHDGAssemblyHelper::_pressure_var
protected

The pressure variable on element interiors.

Definition at line 51 of file NavierStokesStressIPHDGAssemblyHelper.h.

◆ _rz_radial_coord

const unsigned int NavierStokesStressIPHDGAssemblyHelper::_rz_radial_coord
protected

The radial coordinate index for RZ coordinate systems.

Definition at line 64 of file NavierStokesStressIPHDGAssemblyHelper.h.

Referenced by scalarVolume().

The documentation for this class was generated from the following files:
Generated on Mon Nov 3 2025 20:28:06 for https://mooseframework.inl.gov by   doxygen 1.8.14