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

Implements the steady incompressible Navier-Stokes equations for a hybridized discretization. More...

#include <NavierStokesLHDGKernel.h>

Inheritance diagram for NavierStokesLHDGKernel:
[legend]

Public Types

enum  ResidualTagType { ResidualTagType::NonReference, ResidualTagType::Reference }
 
typedef DataFileName DataFileParameterType
 

Public Member Functions

 NavierStokesLHDGKernel (const InputParameters &parameters)
 
virtual void computeResidual () override
 
virtual void computeJacobian () override
 Compute this object's entire element interior Jacobian, both on- and off-diagonal. More...
 
virtual void computeOffDiagJacobian (unsigned int jvar) override
 Forwards to computeJacobian() the first time this is called for a given element. More...
 
virtual void computeResidualOnSide () override
 
virtual void computeJacobianOnSide () override
 
virtual void initialSetup () override
 
virtual void jacobianSetup () override
 
virtual std::set< std::string > additionalROVariables () override
 
virtual void computeResidualAndJacobianOnSide ()
 
virtual void computeResidualAndJacobian () override
 
virtual void computeOffDiagJacobianScalar (unsigned int jvar) override
 
virtual const MooseVariable & variable () const override
 
virtual void computeNonlocalJacobian ()
 
virtual void computeNonlocalOffDiagJacobian (unsigned int)
 
const SubProblem & subProblem () const
 
virtual void prepareShapes (unsigned int var_num)
 
virtual bool enabled () const
 
std::shared_ptr< MooseObject > getSharedPtr ()
 
std::shared_ptr< const MooseObject > getSharedPtr () const
 
MooseApp & getMooseApp () const
 
const std::string & type () const
 
virtual const std::string & name () const
 
std::string typeAndName () const
 
std::string errorPrefix (const std::string &error_type) const
 
void callMooseError (std::string msg, const bool with_prefix) const
 
MooseObjectParameterName uniqueParameterName (const std::string &parameter_name) const
 
const InputParameters & parameters () const
 
MooseObjectName uniqueName () const
 
const T & getParam (const std::string &name) const
 
std::vector< std::pair< T1, T2 > > getParam (const std::string &param1, const std::string &param2) const
 
const T * queryParam (const std::string &name) const
 
const T & getRenamedParam (const std::string &old_name, const std::string &new_name) const
 
T getCheckedPointerParam (const std::string &name, const std::string &error_string="") const
 
bool isParamValid (const std::string &name) const
 
bool isParamSetByUser (const std::string &nm) const
 
void paramError (const std::string &param, Args... args) const
 
void paramWarning (const std::string &param, Args... args) const
 
void paramInfo (const std::string &param, Args... args) const
 
void connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const
 
void mooseError (Args &&... args) const
 
void mooseErrorNonPrefixed (Args &&... args) const
 
void mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseInfo (Args &&... args) const
 
std::string getDataFileName (const std::string &param) const
 
std::string getDataFileNameByName (const std::string &relative_path) const
 
std::string getDataFilePath (const std::string &relative_path) const
 
virtual void timestepSetup ()
 
virtual void residualSetup ()
 
virtual void subdomainSetup ()
 
virtual void customSetup (const ExecFlagType &)
 
const ExecFlagEnum & getExecuteOnEnum () const
 
const Function & getFunction (const std::string &name) const
 
const Function & getFunctionByName (const FunctionName &name) const
 
bool hasFunction (const std::string &param_name) const
 
bool hasFunctionByName (const FunctionName &name) const
 
UserObjectName getUserObjectName (const std::string &param_name) const
 
const T & getUserObject (const std::string &param_name, bool is_dependency=true) const
 
const T & getUserObjectByName (const UserObjectName &object_name, bool is_dependency=true) const
 
const UserObject & getUserObjectBase (const std::string &param_name, bool is_dependency=true) const
 
const UserObject & getUserObjectBaseByName (const UserObjectName &object_name, bool is_dependency=true) const
 
bool isImplicit ()
 
Moose::StateArg determineState () const
 
bool isDefaultPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
bool hasPostprocessor (const std::string &param_name, const unsigned int index=0) const
 
bool hasPostprocessorByName (const PostprocessorName &name) const
 
std::size_t coupledPostprocessors (const std::string &param_name) const
 
const PostprocessorName & getPostprocessorName (const std::string &param_name, const unsigned int index=0) const
 
const VectorPostprocessorValue & getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValue & getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValue & getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValue & getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValue & getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValue & getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValue & getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValue & getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const ScatterVectorPostprocessorValue & getScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValue & getScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValue & getScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValue & getScatterVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
bool hasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const
 
bool hasVectorPostprocessor (const std::string &param_name) const
 
bool hasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
bool hasVectorPostprocessorByName (const VectorPostprocessorName &name) const
 
const VectorPostprocessorName & getVectorPostprocessorName (const std::string &param_name) const
 
void setRandomResetFrequency (ExecFlagType exec_flag)
 
unsigned long getRandomLong () const
 
Real getRandomReal () const
 
unsigned int getSeed (std::size_t id)
 
unsigned int getMasterSeed () const
 
bool isNodal () const
 
ExecFlagType getResetOnTime () const
 
void setRandomDataPointer (RandomData *random_data)
 
virtual void meshChanged ()
 
void useVectorTag (const TagName &tag_name, VectorTagsKey)
 
void useVectorTag (TagID tag_id, VectorTagsKey)
 
void useMatrixTag (const TagName &tag_name, MatrixTagsKey)
 
void useMatrixTag (TagID tag_id, MatrixTagsKey)
 
bool isVectorTagged ()
 
bool isMatrixTagged ()
 
bool hasVectorTags () const
 
const std::set< TagID > & getVectorTags (VectorTagsKey) const
 
const std::set< TagID > & getMatrixTags (MatrixTagsKey) const
 
const std::vector< SubdomainName > & blocks () const
 
unsigned int numBlocks () const
 
virtual const std::set< SubdomainID > & blockIDs () const
 
unsigned int blocksMaxDimension () const
 
bool hasBlocks (const SubdomainName &name) const
 
bool hasBlocks (const std::vector< SubdomainName > &names) const
 
bool hasBlocks (const std::set< SubdomainName > &names) const
 
bool hasBlocks (SubdomainID id) const
 
bool hasBlocks (const std::vector< SubdomainID > &ids) const
 
bool hasBlocks (const std::set< SubdomainID > &ids) const
 
bool isBlockSubset (const std::set< SubdomainID > &ids) const
 
bool isBlockSubset (const std::vector< SubdomainID > &ids) const
 
bool hasBlockMaterialProperty (const std::string &prop_name)
 
const std::set< SubdomainID > & meshBlockIDs () const
 
virtual bool blockRestricted () const
 
virtual void checkVariable (const MooseVariableFieldBase &variable) const
 
virtual const VariableValue & coupledValueByName (const std::string &var_name)
 
virtual const ArrayVariableValue & coupledArrayValueByName (const std::string &var_name)
 
const std::unordered_map< std::string, std::vector< MooseVariableFieldBase *> > & getCoupledVars () const
 
const std::vector< MooseVariableFieldBase *> & getCoupledMooseVars () const
 
const std::vector< MooseVariable *> & getCoupledStandardMooseVars () const
 
const std::vector< VectorMooseVariable *> & getCoupledVectorMooseVars () const
 
const std::vector< ArrayMooseVariable *> & getCoupledArrayMooseVars () const
 
void addFEVariableCoupleableVectorTag (TagID tag)
 
void addFEVariableCoupleableMatrixTag (TagID tag)
 
std::set< TagID > & getFEVariableCoupleableVectorTags ()
 
const std::set< TagID > & getFEVariableCoupleableVectorTags () const
 
std::set< TagID > & getFEVariableCoupleableMatrixTags ()
 
const std::set< TagID > & getFEVariableCoupleableMatrixTags () const
 
auto & getWritableCoupledVariables () const
 
bool hasWritableCoupledVariables () const
 
const ADVariableValue * getADDefaultValue (const std::string &var_name) const
 
const ADVectorVariableValue * getADDefaultVectorValue (const std::string &var_name) const
 
const ADVariableGradient & getADDefaultGradient () const
 
const ADVectorVariableGradient & getADDefaultVectorGradient () const
 
const ADVariableSecond & getADDefaultSecond () const
 
const ADVectorVariableCurl & getADDefaultCurl () const
 
const std::vector< MooseVariableScalar *> & getCoupledMooseScalarVars ()
 
const std::set< TagID > & getScalarVariableCoupleableVectorTags () const
 
const std::set< TagID > & getScalarVariableCoupleableMatrixTags () const
 
const std::set< MooseVariableFieldBase *> & getMooseVariableDependencies () const
 
std::set< MooseVariableFieldBase *> checkAllVariables (const DofObjectType &dof_object, const std::set< MooseVariableFieldBase * > &vars_to_omit={})
 
std::set< MooseVariableFieldBase *> checkVariables (const DofObjectType &dof_object, const std::set< MooseVariableFieldBase * > &vars_to_check)
 
void addMooseVariableDependency (MooseVariableFieldBase *var)
 
void addMooseVariableDependency (const std::vector< MooseVariableFieldBase * > &vars)
 
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)
 
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)
 
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)
 
virtual unsigned int getElementIDIndex (const std::string &id_parameter_name, unsigned int comp=0) const
 
virtual unsigned int getElementIDIndexByName (const std::string &id_name) const
 
virtual const dof_id_type & getElementID (const std::string &id_parameter_name, unsigned int comp=0) const
 
dof_id_type getElementID (const Elem *elem, unsigned int elem_id_index) const
 
virtual const dof_id_type & getElementIDNeighbor (const std::string &id_parameter_name, unsigned int comp=0) const
 
virtual const dof_id_type & getElementIDByName (const std::string &id_name) const
 
virtual const dof_id_type & getElementIDNeighborByName (const std::string &id_name) const
 
bool hasElementID (const std::string &id_name) const
 
dof_id_type maxElementID (unsigned int elem_id_index) const
 
dof_id_type minElementID (unsigned int elem_id_index) const
 
bool areElemIDsIdentical (const std::string &id_name1, const std::string &id_name2) const
 
std::unordered_map< dof_id_type, std::set< dof_id_type > > getElemIDMapping (const std::string &id_name1, const std::string &id_name2) const
 
std::set< dof_id_type > getAllElemIDs (unsigned int elem_id_index) const
 
std::set< dof_id_type > getElemIDsOnBlocks (unsigned int elem_id_index, const std::set< SubdomainID > &blks) const
 
MooseVariableBase * mooseVariableBase () const
 
MooseVariableField< Real > & mooseVariableField ()
 
MooseVariableFE< Real > * mooseVariable () const
 
MooseVariableFV< Real > * mooseVariableFV () const
 
MooseLinearVariableFV< Real > * mooseLinearVariableFV () const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
const PostprocessorValue & getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValue & getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValue & getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValue & getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValue & getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValue & getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
virtual const PostprocessorValue & getPostprocessorValueByName (const PostprocessorName &name) const
 
virtual const PostprocessorValue & getPostprocessorValueByName (const PostprocessorName &name) const
 
const PostprocessorValue & getPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValue & getPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValue & getPostprocessorValueOlderByName (const PostprocessorName &name) const
 
const PostprocessorValue & getPostprocessorValueOlderByName (const PostprocessorName &name) const
 
bool isVectorPostprocessorDistributed (const std::string &param_name) const
 
bool isVectorPostprocessorDistributed (const std::string &param_name) const
 
bool isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const
 
bool isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const
 
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 hasGenericMaterialProperty (const std::string &name)
 
bool hasGenericMaterialProperty (const std::string &name)
 
bool hasGenericMaterialPropertyByName (const std::string &name)
 
bool hasGenericMaterialPropertyByName (const std::string &name)
 
PenetrationLocator & getPenetrationLocator (const BoundaryName &primary, const BoundaryName &secondary, Order order)
 
PenetrationLocator & getQuadraturePenetrationLocator (const BoundaryName &primary, const BoundaryName &secondary, Order order)
 
NearestNodeLocator & getNearestNodeLocator (const BoundaryName &primary, const BoundaryName &secondary)
 
NearestNodeLocator & getQuadratureNearestNodeLocator (const BoundaryName &primary, const BoundaryName &secondary)
 
bool requiresGeometricSearch () const
 
const Parallel::Communicator & comm () const
 
processor_id_type n_processors () const
 
processor_id_type processor_id () const
 
void checkCoupling ()
 

Static Public Member Functions

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

Public Attributes

const ConsoleStream _console
 

Static Public Attributes

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

Protected Member Functions

virtual Real computeQpResidual () override
 
virtual Real computeQpJacobian ()
 
virtual Real computeQpOffDiagJacobian (unsigned int)
 
virtual Real computeQpOffDiagJacobianScalar (unsigned int)
 
virtual RealEigenVector computeQpOffDiagJacobianArray (const ArrayMooseVariable &jvar)
 
virtual void precalculateResidual ()
 
virtual void precalculateJacobian ()
 
virtual void precalculateOffDiagJacobian (unsigned int)
 
const MooseVariableFieldBase & getVariable (unsigned int jvar_num) const
 
virtual void addUserObjectDependencyHelper (const UserObject &) const
 
virtual void addPostprocessorDependencyHelper (const PostprocessorName &) const
 
virtual void addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &) const
 
T & declareRestartableData (const std::string &data_name, Args &&... args)
 
ManagedValue< T > declareManagedRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
const T & getRestartableData (const std::string &data_name) const
 
T & declareRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
T & declareRecoverableData (const std::string &data_name, Args &&... args)
 
T & declareRestartableDataWithObjectName (const std::string &data_name, const std::string &object_name, Args &&... args)
 
T & declareRestartableDataWithObjectNameWithContext (const std::string &data_name, const std::string &object_name, void *context, Args &&... args)
 
std::string restartableName (const std::string &data_name) const
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar)
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)
 
void prepareVectorTagNeighbor (Assembly &assembly, unsigned int ivar)
 
void prepareVectorTagLower (Assembly &assembly, unsigned int ivar)
 
void prepareMatrixTag (Assembly &assembly, unsigned int ivar, unsigned int jvar)
 
void prepareMatrixTag (Assembly &assembly, unsigned int ivar, unsigned int jvar, DenseMatrix< Number > &k) const
 
void prepareMatrixTagNonlocal (Assembly &assembly, unsigned int ivar, unsigned int jvar)
 
void prepareMatrixTagNeighbor (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type)
 
void prepareMatrixTagNeighbor (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type, DenseMatrix< Number > &k) const
 
void prepareMatrixTagLower (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::ConstraintJacobianType type)
 
void accumulateTaggedLocalResidual ()
 
void assignTaggedLocalResidual ()
 
void accumulateTaggedLocalMatrix ()
 
void accumulateTaggedLocalMatrix (Assembly &assembly, unsigned int ivar, unsigned int jvar, const DenseMatrix< Number > &k)
 
void accumulateTaggedLocalMatrix (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type, const DenseMatrix< Number > &k)
 
void accumulateTaggedNonlocalMatrix ()
 
void assignTaggedLocalMatrix ()
 
void addResiduals (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addResiduals (Assembly &assembly, const DenseVector< T > &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addResiduals (Assembly &assembly, const ADResidualsPacket &packet)
 
void addResidualsAndJacobian (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addResidualsAndJacobian (Assembly &assembly, const ADResidualsPacket &packet)
 
void addJacobian (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addJacobian (Assembly &assembly, const ADResidualsPacket &packet)
 
void addJacobian (Assembly &assembly, DenseMatrix< Real > &local_k, const std::vector< dof_id_type > &row_indices, const std::vector< dof_id_type > &column_indices, Real scaling_factor)
 
void addResidualsWithoutConstraints (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addResidualsAndJacobianWithoutConstraints (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addJacobianWithoutConstraints (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addJacobianElement (Assembly &assembly, Real value, dof_id_type row_index, dof_id_type column_index, Real scaling_factor)
 
void setResidual (SystemBase &sys, const T &residual, MooseVariableFE< T > &var)
 
void setResidual (SystemBase &sys, Real residual, dof_id_type dof_index)
 
void setResidual (SystemBase &sys, SetResidualFunctor set_residual_functor)
 
virtual bool hasBlockMaterialPropertyHelper (const std::string &prop_name)
 
void initializeBlockRestrictable (const MooseObject *moose_object)
 
Moose::CoordinateSystemType getBlockCoordSystem ()
 
virtual void coupledCallback (const std::string &, bool) const
 
virtual bool isCoupled (const std::string &var_name, unsigned int i=0) const
 
virtual bool isCoupledConstant (const std::string &var_name) const
 
unsigned int coupledComponents (const std::string &var_name) const
 
VariableName coupledName (const std::string &var_name, unsigned int comp=0) const
 
std::vector< VariableName > coupledNames (const std::string &var_name) const
 
virtual unsigned int coupled (const std::string &var_name, unsigned int comp=0) const
 
std::vector< unsigned int > coupledIndices (const std::string &var_name) const
 
virtual const VariableValue & coupledValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledValues (const std::string &var_name) const
 
std::vector< const VectorVariableValue *> coupledVectorValues (const std::string &var_name) const
 
const GenericVariableValue< is_ad > & coupledGenericValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< false > & coupledGenericValue (const std::string &var_name, unsigned int comp) const
 
const GenericVariableValue< true > & coupledGenericValue (const std::string &var_name, unsigned int comp) const
 
const GenericVectorVariableValue< is_ad > & coupledGenericVectorValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVectorVariableValue< false > & coupledGenericVectorValue (const std::string &var_name, unsigned int comp) const
 
const GenericVectorVariableValue< true > & coupledGenericVectorValue (const std::string &var_name, unsigned int comp) const
 
std::vector< const GenericVariableValue< is_ad > *> coupledGenericValues (const std::string &var_name) const
 
std::vector< const GenericVariableValue< false > *> coupledGenericValues (const std::string &var_name) const
 
std::vector< const GenericVariableValue< true > *> coupledGenericValues (const std::string &var_name) const
 
const GenericVariableValue< is_ad > & coupledGenericDofValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< false > & coupledGenericDofValue (const std::string &var_name, unsigned int comp) const
 
const GenericVariableValue< true > & coupledGenericDofValue (const std::string &var_name, unsigned int comp) const
 
const GenericVariableValue< is_ad > & coupledGenericDot (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< false > & coupledGenericDot (const std::string &var_name, unsigned int comp) const
 
const GenericVariableValue< true > & coupledGenericDot (const std::string &var_name, unsigned int comp) const
 
const GenericVariableValue< is_ad > & coupledGenericDotDot (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< false > & coupledGenericDotDot (const std::string &var_name, unsigned int comp) const
 
const GenericVariableValue< true > & coupledGenericDotDot (const std::string &var_name, unsigned int comp) const
 
virtual const VariableValue & coupledValueLower (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValue & adCoupledValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVariableValue *> adCoupledValues (const std::string &var_name) const
 
const ADVariableValue & adCoupledLowerValue (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableValue & adCoupledVectorValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVectorVariableValue *> adCoupledVectorValues (const std::string &var_name) const
 
virtual const VariableValue & coupledVectorTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const VariableValue & coupledVectorTagValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
std::vector< const VariableValue *> coupledVectorTagValues (const std::string &var_names, TagID tag) const
 
std::vector< const VariableValue *> coupledVectorTagValues (const std::string &var_names, const std::string &tag_name) const
 
virtual const ArrayVariableValue & coupledVectorTagArrayValue (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const ArrayVariableValue & coupledVectorTagArrayValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
std::vector< const ArrayVariableValue *> coupledVectorTagArrayValues (const std::string &var_names, TagID tag) const
 
std::vector< const ArrayVariableValue *> coupledVectorTagArrayValues (const std::string &var_names, const std::string &tag_name) const
 
virtual const VariableGradient & coupledVectorTagGradient (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const VariableGradient & coupledVectorTagGradient (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
std::vector< const VariableGradient *> coupledVectorTagGradients (const std::string &var_names, TagID tag) const
 
std::vector< const VariableGradient *> coupledVectorTagGradients (const std::string &var_names, const std::string &tag_name) const
 
virtual const ArrayVariableGradient & coupledVectorTagArrayGradient (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const ArrayVariableGradient & coupledVectorTagArrayGradient (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
std::vector< const ArrayVariableGradient *> coupledVectorTagArrayGradients (const std::string &var_names, TagID tag) const
 
std::vector< const ArrayVariableGradient *> coupledVectorTagArrayGradients (const std::string &var_names, const std::string &tag_name) const
 
virtual const VariableValue & coupledVectorTagDofValue (const std::string &var_name, TagID tag, unsigned int index=0) const
 
virtual const VariableValue & coupledVectorTagDofValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
const ArrayVariableValue & coupledVectorTagArrayDofValue (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledVectorTagDofValues (const std::string &var_names, TagID tag) const
 
std::vector< const VariableValue *> coupledVectorTagDofValues (const std::string &var_names, const std::string &tag_name) const
 
virtual const VariableValue & coupledMatrixTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const VariableValue & coupledMatrixTagValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
std::vector< const VariableValue *> coupledMatrixTagValues (const std::string &var_names, TagID tag) const
 
std::vector< const VariableValue *> coupledMatrixTagValues (const std::string &var_names, const std::string &tag_name) const
 
virtual const VectorVariableValue & coupledVectorValue (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValue & coupledArrayValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ArrayVariableValue *> coupledArrayValues (const std::string &var_name) const
 
MooseWritableVariable & writableVariable (const std::string &var_name, unsigned int comp=0)
 
virtual VariableValue & writableCoupledValue (const std::string &var_name, unsigned int comp=0)
 
void checkWritableVar (MooseWritableVariable *var)
 
virtual const VariableValue & coupledValueOld (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledValuesOld (const std::string &var_name) const
 
std::vector< const VectorVariableValue *> coupledVectorValuesOld (const std::string &var_name) const
 
virtual const VariableValue & coupledValueOlder (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledValuesOlder (const std::string &var_name) const
 
virtual const VariableValue & coupledValuePreviousNL (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValue & coupledVectorValueOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValue & coupledVectorValueOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValue & coupledArrayValueOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValue & coupledArrayValueOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradient & coupledGradient (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableGradient *> coupledGradients (const std::string &var_name) const
 
const ADVariableGradient & adCoupledGradient (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableGradient & adCoupledGradientDot (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVariableGradient *> adCoupledGradients (const std::string &var_name) const
 
const GenericVariableGradient< is_ad > & coupledGenericGradient (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableGradient< false > & coupledGenericGradient (const std::string &var_name, unsigned int comp) const
 
const GenericVariableGradient< true > & coupledGenericGradient (const std::string &var_name, unsigned int comp) const
 
std::vector< const GenericVariableGradient< is_ad > *> coupledGenericGradients (const std::string &var_name) const
 
std::vector< const GenericVariableGradient< false > *> coupledGenericGradients (const std::string &var_name) const
 
std::vector< const GenericVariableGradient< true > *> coupledGenericGradients (const std::string &var_name) const
 
const ADVectorVariableGradient & adCoupledVectorGradient (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableSecond & adCoupledSecond (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableSecond & adCoupledVectorSecond (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradient & coupledGradientOld (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableGradient *> coupledGradientsOld (const std::string &var_name) const
 
virtual const VariableGradient & coupledGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradient & coupledGradientPreviousNL (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradient & coupledGradientDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradient & coupledGradientDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradient & coupledVectorGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradient & coupledVectorGradientOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradient & coupledVectorGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradient & coupledArrayGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradient & coupledArrayGradientOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradient & coupledArrayGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradient & coupledArrayGradientDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableCurl & coupledCurl (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableCurl & coupledCurlOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableCurl & coupledCurlOlder (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableCurl & adCoupledCurl (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableDivergence & coupledDiv (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableDivergence & coupledDivOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableDivergence & coupledDivOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecond & coupledSecond (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecond & coupledSecondOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecond & coupledSecondOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecond & coupledSecondPreviousNL (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValue & coupledDot (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledDots (const std::string &var_name) const
 
virtual const VariableValue & coupledDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValue & coupledDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValue & coupledDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValue & adCoupledDot (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVariableValue *> adCoupledDots (const std::string &var_name) const
 
const ADVariableValue & adCoupledDotDot (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableValue & adCoupledVectorDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValue & coupledVectorDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValue & coupledVectorDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValue & coupledVectorDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValue & coupledVectorDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValue & coupledVectorDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValue & coupledVectorDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValue & coupledArrayDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValue & coupledArrayDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValue & coupledArrayDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValue & coupledArrayDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValue & coupledDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValue & coupledDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
const VariableValue & coupledArrayDotDu (const std::string &var_name, unsigned int comp=0) const
 
const T & coupledNodalValue (const std::string &var_name, unsigned int comp=0) const
 
const Moose::ADType< T >::type & adCoupledNodalValue (const std::string &var_name, unsigned int comp=0) const
 
const T & coupledNodalValueOld (const std::string &var_name, unsigned int comp=0) const
 
const T & coupledNodalValueOlder (const std::string &var_name, unsigned int comp=0) const
 
const T & coupledNodalValuePreviousNL (const std::string &var_name, unsigned int comp=0) const
 
const T & coupledNodalDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValue & coupledNodalDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValue & coupledNodalDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValue & coupledNodalDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValue & coupledDofValues (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledAllDofValues (const std::string &var_name) const
 
virtual const VariableValue & coupledDofValuesOld (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledAllDofValuesOld (const std::string &var_name) const
 
virtual const VariableValue & coupledDofValuesOlder (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledAllDofValuesOlder (const std::string &var_name) const
 
virtual const ArrayVariableValue & coupledArrayDofValues (const std::string &var_name, unsigned int comp=0) const
 
virtual const ADVariableValue & adCoupledDofValues (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValue & adZeroValue () const
 
const ADVariableGradient & adZeroGradient () const
 
const ADVariableSecond & adZeroSecond () const
 
const GenericVariableValue< is_ad > & genericZeroValue ()
 
const GenericVariableValue< false > & genericZeroValue ()
 
const GenericVariableValue< true > & genericZeroValue ()
 
const GenericVariableGradient< is_ad > & genericZeroGradient ()
 
const GenericVariableGradient< false > & genericZeroGradient ()
 
const GenericVariableGradient< true > & genericZeroGradient ()
 
const GenericVariableSecond< is_ad > & genericZeroSecond ()
 
const GenericVariableSecond< false > & genericZeroSecond ()
 
const GenericVariableSecond< true > & genericZeroSecond ()
 
bool checkVar (const std::string &var_name, unsigned int comp=0, unsigned int comp_bound=0) const
 
const MooseVariableFieldBase * getFEVar (const std::string &var_name, unsigned int comp) const
 
const MooseVariableFieldBase * getFieldVar (const std::string &var_name, unsigned int comp) const
 
MooseVariableFieldBase * getFieldVar (const std::string &var_name, unsigned int comp)
 
const T * getVarHelper (const std::string &var_name, unsigned int comp) const
 
T * getVarHelper (const std::string &var_name, unsigned int comp)
 
MooseVariable * getVar (const std::string &var_name, unsigned int comp)
 
const MooseVariable * getVar (const std::string &var_name, unsigned int comp) const
 
VectorMooseVariable * getVectorVar (const std::string &var_name, unsigned int comp)
 
const VectorMooseVariable * getVectorVar (const std::string &var_name, unsigned int comp) const
 
ArrayMooseVariable * getArrayVar (const std::string &var_name, unsigned int comp)
 
const ArrayMooseVariable * getArrayVar (const std::string &var_name, unsigned int comp) const
 
void validateExecutionerType (const std::string &name, const std::string &fn_name) const
 
std::vector< T > coupledVectorHelper (const std::string &var_name, const Func &func) const
 
bool isCoupledScalar (const std::string &var_name, unsigned int i=0) const
 
unsigned int coupledScalarComponents (const std::string &var_name) const
 
unsigned int coupledScalar (const std::string &var_name, unsigned int comp=0) const
 
libMesh::Order coupledScalarOrder (const std::string &var_name, unsigned int comp=0) const
 
const VariableValue & coupledScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValue & adCoupledScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< is_ad > & coupledGenericScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< false > & coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const
 
const GenericVariableValue< true > & coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const
 
const VariableValue & coupledVectorTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValue & coupledMatrixTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValue & coupledScalarValueOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValue & coupledScalarValueOlder (const std::string &var_name, unsigned int comp=0) const
 
const VariableValue & coupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValue & adCoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValue & coupledScalarDotDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValue & coupledScalarDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValue & coupledScalarDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValue & coupledScalarDotDu (const std::string &var_name, unsigned int comp=0) const
 
const VariableValue & coupledScalarDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
const MooseVariableScalar * getScalarVar (const std::string &var_name, unsigned int comp) const
 
virtual void checkMaterialProperty (const std::string &name, const unsigned int state)
 
void markMatPropRequested (const std::string &)
 
MaterialPropertyName getMaterialPropertyName (const std::string &name) const
 
void checkExecutionStage ()
 
virtual const OutputTools< Real >::VariableValue & value ()
 
virtual const OutputTools< Real >::VariableValue & valueOld ()
 
virtual const OutputTools< Real >::VariableValue & valueOlder ()
 
virtual const OutputTools< Real >::VariableValue & dot ()
 
virtual const OutputTools< Real >::VariableValue & dotDot ()
 
virtual const OutputTools< Real >::VariableValue & dotOld ()
 
virtual const OutputTools< Real >::VariableValue & dotDotOld ()
 
virtual const VariableValue & dotDu ()
 
virtual const VariableValue & dotDotDu ()
 
virtual const OutputTools< Real >::VariableGradient & gradient ()
 
virtual const OutputTools< Real >::VariableGradient & gradientOld ()
 
virtual const OutputTools< Real >::VariableGradient & gradientOlder ()
 
virtual const OutputTools< Real >::VariableSecond & second ()
 
virtual const OutputTools< Real >::VariableSecond & secondOld ()
 
virtual const OutputTools< Real >::VariableSecond & secondOlder ()
 
virtual const OutputTools< Real >::VariableTestSecond & secondTest ()
 
virtual const OutputTools< Real >::VariableTestSecond & secondTestFace ()
 
virtual const OutputTools< Real >::VariablePhiSecond & secondPhi ()
 
virtual const OutputTools< Real >::VariablePhiSecond & secondPhiFace ()
 
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)
 
RealVectorValue rhoVelCrossVelResidual (const MooseArray< Number > &u_sol, const MooseArray< Number > &v_sol, const unsigned int qp, const unsigned int vel_component)
 
RealVectorValue rhoVelCrossVelJacobian (const MooseArray< Number > &u_sol, const MooseArray< Number > &v_sol, const unsigned int qp, const unsigned int vel_component, const unsigned int vel_j_component, const MooseArray< std::vector< Real >> &phi, const unsigned int j)
 
void scalarVolumeResidual (const MooseArray< Gradient > &vel_gradient, const unsigned int vel_component, const Moose::Functor< Real > &body_force, const MooseArray< Real > &JxW, const libMesh::QBase &qrule, const Elem *const current_elem, const MooseArray< Point > &q_point, DenseVector< Number > &scalar_re)
 Compute the volumetric contributions to a velocity residual for a provided velocity gradient and stress. More...
 
void scalarVolumeResidual (const MooseArray< Gradient > &vector_field, const Moose::Functor< Real > &source, const MooseArray< Real > &JxW, const libMesh::QBase &qrule, const Elem *const current_elem, const MooseArray< Point > &q_point, DenseVector< Number > &scalar_re)
 
void scalarVolumeJacobian (const unsigned int vel_component, const MooseArray< Real > &JxW, const libMesh::QBase &qrule, DenseMatrix< Number > &scalar_vector_jac, DenseMatrix< Number > &scalar_u_vel_jac, DenseMatrix< Number > &scalar_v_vel_jac, DenseMatrix< Number > &scalar_p_jac)
 Compute the volumetric contributions to a velocity Jacobian. More...
 
void scalarVolumeJacobian (const MooseArray< Real > &JxW, const libMesh::QBase &qrule, DenseMatrix< Number > &scalar_vector_jac)
 
void pressureVolumeResidual (const Moose::Functor< Real > &pressure_mms_forcing_function, const MooseArray< Real > &JxW, const libMesh::QBase &qrule, const Elem *const current_elem, const MooseArray< Point > &q_point, DenseVector< Number > &pressure_re, DenseVector< Number > &global_lm_re)
 Compute the volumetric contributions to the pressure residual, e.g. More...
 
void pressureVolumeJacobian (const MooseArray< Real > &JxW, const libMesh::QBase &qrule, DenseMatrix< Number > &p_u_vel_jac, DenseMatrix< Number > &p_v_vel_jac, DenseMatrix< Number > &p_global_lm_jac, DenseMatrix< Number > &global_lm_p_jac)
 Compute the volumetric contributions to the pressure Jacobian, e.g. More...
 
void pressureFaceResidual (const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseVector< Number > &pressure_re)
 
void pressureFaceJacobian (const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseMatrix< Number > &p_lm_u_vel_jac, DenseMatrix< Number > &p_lm_v_vel_jac)
 
void scalarFaceResidual (const MooseArray< Gradient > &vector_sol, const MooseArray< Number > &scalar_sol, const MooseArray< Number > &lm_sol, const unsigned int vel_component, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseVector< Number > &scalar_re)
 
void scalarFaceResidual (const MooseArray< Gradient > &vector_sol, const MooseArray< Number > &scalar_sol, const MooseArray< Number > &lm_sol, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseVector< Number > &scalar_re)
 
void scalarFaceJacobian (const unsigned int vel_component, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseMatrix< Number > &scalar_vector_jac, DenseMatrix< Number > &scalar_scalar_jac, DenseMatrix< Number > &scalar_lm_jac, DenseMatrix< Number > &scalar_p_jac, DenseMatrix< Number > &scalar_lm_u_vel_jac, DenseMatrix< Number > &scalar_lm_v_vel_jac)
 
void scalarFaceJacobian (const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseMatrix< Number > &scalar_vector_jac, DenseMatrix< Number > &scalar_scalar_jac, DenseMatrix< Number > &scalar_lm_jac)
 
void lmFaceResidual (const MooseArray< Gradient > &vector_sol, const MooseArray< Number > &scalar_sol, const MooseArray< Number > &lm_sol, const unsigned int vel_component, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, const Elem *const neigh, DenseVector< Number > &lm_re)
 
void lmFaceResidual (const MooseArray< Gradient > &vector_sol, const MooseArray< Number > &scalar_sol, const MooseArray< Number > &lm_sol, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseVector< Number > &lm_re)
 
void lmFaceJacobian (const unsigned int vel_component, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, const Elem *const neigh, DenseMatrix< Number > &lm_vec_jac, DenseMatrix< Number > &lm_scalar_jac, DenseMatrix< Number > &lm_lm_jac, DenseMatrix< Number > &lm_p_jac, DenseMatrix< Number > &lm_lm_u_vel_jac, DenseMatrix< Number > &lm_lm_v_vel_jac)
 
void lmFaceJacobian (const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseMatrix< Number > &lm_vec_jac, DenseMatrix< Number > &lm_scalar_jac, DenseMatrix< Number > &lm_lm_jac)
 
void pressureDirichletResidual (const std::array< const Moose::Functor< Real > *, 3 > &dirichlet_vel, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, const Elem *const current_elem, const unsigned int current_side, const MooseArray< Point > &q_point_face, DenseVector< Number > &pressure_re)
 
void scalarDirichletResidual (const MooseArray< Gradient > &vector_sol, const MooseArray< Number > &scalar_sol, const unsigned int vel_component, const std::array< const Moose::Functor< Real > *, 3 > &dirichlet_vel, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, const Elem *const current_elem, const unsigned int current_side, const MooseArray< Point > &q_point_face, DenseVector< Number > &scalar_re)
 
void scalarDirichletResidual (const MooseArray< Gradient > &vector_sol, const MooseArray< Number > &scalar_sol, const Moose::Functor< Real > &dirichlet_value, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, const Elem *const current_elem, const unsigned int current_side, const MooseArray< Point > &q_point_face, DenseVector< Number > &scalar_re)
 
void scalarDirichletJacobian (const unsigned int vel_component, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseMatrix< Number > &scalar_vector_jac, DenseMatrix< Number > &scalar_scalar_jac, DenseMatrix< Number > &scalar_pressure_jac)
 
void scalarDirichletJacobian (const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseMatrix< Number > &scalar_vector_jac, DenseMatrix< Number > &scalar_scalar_jac)
 
std::string deduceFunctorName (const std::string &name) const
 
void vectorVolumeResidual (const MooseArray< Gradient > &vector_sol, const MooseArray< Number > &scalar_sol, const MooseArray< Real > &JxW, const libMesh::QBase &qrule, DenseVector< Number > &vector_re)
 
void vectorVolumeJacobian (const MooseArray< Real > &JxW, const libMesh::QBase &qrule, DenseMatrix< Number > &vector_vector_jac, DenseMatrix< Number > &vector_scalar_jac)
 
void vectorFaceResidual (const MooseArray< Number > &lm_sol, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseVector< Number > &vector_re)
 
void vectorFaceJacobian (const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseMatrix< Number > &vector_lm_jac)
 
void vectorDirichletResidual (const Moose::Functor< Real > &dirichlet_value, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, const Elem *const current_elem, const unsigned int current_side, const MooseArray< Point > &q_point_face, DenseVector< Number > &vector_re)
 
void createIdentityResidual (const MooseArray< Real > &JxW, const libMesh::QBase &qrule, const MooseArray< std::vector< Real >> &phi, const MooseArray< Number > &sol, DenseVector< Number > &re)
 
void createIdentityJacobian (const MooseArray< Real > &JxW, const libMesh::QBase &qrule, const MooseArray< std::vector< Real >> &phi, DenseMatrix< Number > &ke)
 
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)
 

Protected Attributes

const Moose::Functor< Real > & _body_force_x
 
const Moose::Functor< Real > & _body_force_y
 
const Moose::Functor< Real > & _body_force_z
 
std::vector< const Moose::Functor< Real > * > _body_forces
 
const Moose::Functor< Real > & _pressure_mms_forcing_function
 
const QBase *const & _qrule_face
 The face quadrature rule. More...
 
const MooseArray< Point > & _q_point_face
 The physical locations of the quadrature points on the face. More...
 
const MooseArray< Real > & _JxW_face
 transformed Jacobian weights on the current element face More...
 
const MooseArray< Point > & _normals
 face normals More...
 
const unsigned int & _current_side
 The current side when doing face evaluations. More...
 
MooseVariable & _var
 
const VariableTestValue & _test
 
const VariableTestGradient & _grad_test
 
const VariablePhiValue & _phi
 
const VariablePhiGradient & _grad_phi
 
const VariableValue & _u
 
const VariableGradient & _grad_u
 
const Elem *const & _current_elem
 
const Real & _current_elem_volume
 
unsigned int _qp
 
const MooseArray< Point > & _q_point
 
const QBase *const & _qrule
 
const MooseArray< Real > & _JxW
 
const MooseArray< Real > & _coord
 
unsigned int _i
 
unsigned int _j
 
bool _has_save_in
 
std::vector< MooseVariableFEBase *> _save_in
 
std::vector< AuxVariableName > _save_in_strings
 
bool _has_diag_save_in
 
std::vector< MooseVariableFEBase *> _diag_save_in
 
std::vector< AuxVariableName > _diag_save_in_strings
 
std::vector< unsigned int > _displacements
 
const bool _use_displaced_mesh
 
SubProblem & _subproblem
 
FEProblemBase & _fe_problem
 
SystemBase & _sys
 
THREAD_ID _tid
 
Assembly & _assembly
 
MooseMesh & _mesh
 
const bool & _enabled
 
MooseApp & _app
 
const std::string _type
 
const std::string _name
 
const InputParameters & _pars
 
Factory & _factory
 
ActionFactory & _action_factory
 
const ExecFlagEnum & _execute_enum
 
const ExecFlagType & _current_execute_flag
 
const InputParameters & _ti_params
 
FEProblemBase & _ti_feproblem
 
bool _is_implicit
 
Real & _t
 
const Real & _t_old
 
int & _t_step
 
Real & _dt
 
Real & _dt_old
 
bool _is_transient
 
MooseApp & _restartable_app
 
const std::string _restartable_system_name
 
const THREAD_ID _restartable_tid
 
const bool _restartable_read_only
 
FEProblemBase & _mci_feproblem
 
DenseVector< Number > _local_re
 
DenseMatrix< Number > _local_ke
 
DenseMatrix< Number > _nonlocal_ke
 
const MaterialData * _blk_material_data
 
const InputParameters & _c_parameters
 
const std::string & _c_name
 
const std::string & _c_type
 
FEProblemBase & _c_fe_problem
 
const SystemBase *const _c_sys
 
std::unordered_map< std::string, std::vector< MooseVariableFieldBase *> > _coupled_vars
 
std::vector< MooseVariableFieldBase *> _coupled_moose_vars
 
std::vector< MooseVariable *> _coupled_standard_moose_vars
 
std::vector< VectorMooseVariable *> _coupled_vector_moose_vars
 
std::vector< ArrayMooseVariable *> _coupled_array_moose_vars
 
std::vector< MooseVariableFV< Real > *> _coupled_standard_fv_moose_vars
 
std::vector< MooseLinearVariableFV< Real > *> _coupled_standard_linear_fv_moose_vars
 
const std::unordered_map< std::string, std::string > & _new_to_deprecated_coupled_vars
 
bool _c_nodal
 
bool _c_is_implicit
 
const bool _c_allow_element_to_nodal_coupling
 
THREAD_ID _c_tid
 
std::unordered_map< std::string, std::vector< std::unique_ptr< VariableValue > > > _default_value
 
std::unordered_map< std::string, std::unique_ptr< MooseArray< ADReal > > > _ad_default_value
 
std::unordered_map< std::string, std::unique_ptr< VectorVariableValue > > _default_vector_value
 
std::unordered_map< std::string, std::unique_ptr< ArrayVariableValue > > _default_array_value
 
std::unordered_map< std::string, std::unique_ptr< MooseArray< ADRealVectorValue > > > _ad_default_vector_value
 
VariableValue _default_value_zero
 
VariableGradient _default_gradient
 
MooseArray< ADRealVectorValue > _ad_default_gradient
 
MooseArray< ADRealTensorValue > _ad_default_vector_gradient
 
VariableSecond _default_second
 
MooseArray< ADRealTensorValue > _ad_default_second
 
MooseArray< ADRealVectorValue > _ad_default_curl
 
const VariableValue & _zero
 
const VariablePhiValue & _phi_zero
 
const MooseArray< ADReal > & _ad_zero
 
const VariableGradient & _grad_zero
 
const MooseArray< ADRealVectorValue > & _ad_grad_zero
 
const VariablePhiGradient & _grad_phi_zero
 
const VariableSecond & _second_zero
 
const MooseArray< ADRealTensorValue > & _ad_second_zero
 
const VariablePhiSecond & _second_phi_zero
 
const VectorVariableValue & _vector_zero
 
const VectorVariableCurl & _vector_curl_zero
 
VectorVariableValue _default_vector_value_zero
 
VectorVariableGradient _default_vector_gradient
 
VectorVariableCurl _default_vector_curl
 
VectorVariableDivergence _default_div
 
ArrayVariableValue _default_array_value_zero
 
ArrayVariableGradient _default_array_gradient
 
bool _coupleable_neighbor
 
FEProblemBase & _sc_fe_problem
 
const THREAD_ID _sc_tid
 
const Real & _real_zero
 
const VariableValue & _scalar_zero
 
const Point & _point_zero
 
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 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
 
GeometricSearchData & _geometric_search_data
 
bool _requires_geometric_search
 
bool _nodal
 
MooseVariableFE< Real > * _variable
 
MooseVariableFV< Real > * _fv_variable
 
MooseLinearVariableFV< Real > * _linear_fv_variable
 
MooseVariableField< Real > * _field_variable
 
Assembly * _mvi_assembly
 
const Parallel::Communicator & _communicator
 
const MooseVariableFE< Real > & _v_var
 
const MooseVariableFE< Real > *const _w_var
 
const MooseVariableFE< RealVectorValue > & _grad_v_var
 
const MooseVariableFE< RealVectorValue > *const _grad_w_var
 
const MooseVariableFE< Real > & _v_face_var
 
const MooseVariableFE< Real > *const _w_face_var
 
const MooseVariableFE< Real > & _pressure_var
 
const MooseVariableScalar *const _enclosure_lm_var
 
const std::vector< dof_id_type > & _qv_dof_indices
 Containers for dof indices. More...
 
const std::vector< dof_id_type > & _v_dof_indices
 
const std::vector< dof_id_type > & _lm_v_dof_indices
 
const std::vector< dof_id_type > *const _qw_dof_indices
 
const std::vector< dof_id_type > *const _w_dof_indices
 
const std::vector< dof_id_type > *const _lm_w_dof_indices
 
const std::vector< dof_id_type > & _p_dof_indices
 
const std::vector< dof_id_type > *const _global_lm_dof_indices
 
const MooseArray< Gradient > & _qv_sol
 local solutions at quadrature points More...
 
const MooseArray< Number > & _v_sol
 
const MooseArray< Number > & _lm_v_sol
 
const MooseArray< Gradient > *const _qw_sol
 
const MooseArray< Number > *const _w_sol
 
const MooseArray< Number > *const _lm_w_sol
 
const MooseArray< Number > & _p_sol
 
const MooseArray< Number > *const _global_lm_dof_value
 
const Real _rho
 The density. More...
 
DenseVector< Number > _grad_u_vel_re
 
DenseVector< Number > _grad_v_vel_re
 
DenseVector< Number > _u_vel_re
 
DenseVector< Number > _v_vel_re
 
DenseVector< Number > _lm_u_vel_re
 
DenseVector< Number > _lm_v_vel_re
 
DenseVector< Number > _p_re
 
DenseVector< Number > _global_lm_re
 
DenseMatrix< Number > _grad_u_grad_u_jac
 
DenseMatrix< Number > _grad_u_u_jac
 
DenseMatrix< Number > _grad_v_grad_v_jac
 
DenseMatrix< Number > _grad_v_v_jac
 
DenseMatrix< Number > _u_grad_u_jac
 
DenseMatrix< Number > _v_grad_v_jac
 
DenseMatrix< Number > _u_u_jac
 
DenseMatrix< Number > _u_v_jac
 
DenseMatrix< Number > _v_u_jac
 
DenseMatrix< Number > _v_v_jac
 
DenseMatrix< Number > _u_p_jac
 
DenseMatrix< Number > _v_p_jac
 
DenseMatrix< Number > _p_u_jac
 
DenseMatrix< Number > _p_v_jac
 
DenseMatrix< Number > _p_global_lm_jac
 
DenseMatrix< Number > _global_lm_p_jac
 
DenseMatrix< Number > _grad_u_lm_u_jac
 
DenseMatrix< Number > _grad_v_lm_v_jac
 
DenseMatrix< Number > _u_lm_u_jac
 
DenseMatrix< Number > _v_lm_v_jac
 
DenseMatrix< Number > _u_lm_v_jac
 
DenseMatrix< Number > _v_lm_u_jac
 
DenseMatrix< Number > _lm_u_grad_u_jac
 
DenseMatrix< Number > _lm_v_grad_v_jac
 
DenseMatrix< Number > _lm_u_u_jac
 
DenseMatrix< Number > _lm_v_v_jac
 
DenseMatrix< Number > _lm_u_lm_u_jac
 
DenseMatrix< Number > _lm_v_lm_v_jac
 
DenseMatrix< Number > _lm_u_p_jac
 
DenseMatrix< Number > _lm_v_p_jac
 
DenseMatrix< Number > _lm_u_lm_v_jac
 
DenseMatrix< Number > _lm_v_lm_u_jac
 
DenseMatrix< Number > _p_lm_u_jac
 
DenseMatrix< Number > _p_lm_v_jac
 
const MooseVariableFE< Real > & _u_var
 
const MooseVariableFE< RealVectorValue > & _grad_u_var
 
const MooseVariableFE< Real > & _u_face_var
 
const std::vector< dof_id_type > & _qu_dof_indices
 
const std::vector< dof_id_type > & _u_dof_indices
 
const std::vector< dof_id_type > & _lm_u_dof_indices
 
const MooseArray< libMesh::Gradient > & _qu_sol
 
const MooseArray< Number > & _u_sol
 
const MooseArray< Number > & _lm_u_sol
 
const MooseArray< std::vector< RealVectorValue > > & _vector_phi
 
const MooseArray< std::vector< Real > > & _scalar_phi
 
const MooseArray< std::vector< RealVectorValue > > & _grad_scalar_phi
 
const MooseArray< std::vector< Real > > & _div_vector_phi
 
const MooseArray< std::vector< RealVectorValue > > & _vector_phi_face
 
const MooseArray< std::vector< Real > > & _scalar_phi_face
 
const MooseArray< std::vector< Real > > & _lm_phi_face
 
const MaterialProperty< Real > & _diff
 
const TransientInterface & _ti
 
const Real _tau
 
const Elem * _cached_elem
 
DenseVector< Number > _vector_re
 
DenseVector< Number > _scalar_re
 
DenseVector< Number > _lm_re
 
DenseMatrix< Number > _vector_vector_jac
 
DenseMatrix< Number > _vector_scalar_jac
 
DenseMatrix< Number > _scalar_vector_jac
 
DenseMatrix< Number > _scalar_scalar_jac
 
DenseMatrix< Number > _scalar_lm_jac
 
DenseMatrix< Number > _lm_scalar_jac
 
DenseMatrix< Number > _lm_lm_jac
 
DenseMatrix< Number > _vector_lm_jac
 
DenseMatrix< Number > _lm_vector_jac
 

Static Protected Attributes

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

Detailed Description

Implements the steady incompressible Navier-Stokes equations for a hybridized discretization.

Definition at line 20 of file NavierStokesLHDGKernel.h.

Constructor & Destructor Documentation

◆ NavierStokesLHDGKernel()

NavierStokesLHDGKernel::NavierStokesLHDGKernel ( const InputParameters &  parameters)

Definition at line 38 of file NavierStokesLHDGKernel.C.

39  : HDGKernel(parameters),
40  NavierStokesLHDGAssemblyHelper(this, this, this, this, _fe_problem, _sys, _mesh, _tid),
41  // body forces
42  _body_force_x(getFunctor<Real>("body_force_x")),
43  _body_force_y(getFunctor<Real>("body_force_y")),
44  _body_force_z(getFunctor<Real>("body_force_z")),
45  _pressure_mms_forcing_function(getFunctor<Real>("pressure_mms_forcing_function")),
46  _qrule_face(_assembly.qRuleFace()),
47  _q_point_face(_assembly.qPointsFace()),
48  _JxW_face(_assembly.JxWFace()),
49  _normals(_assembly.normals()),
50  _current_side(_assembly.side())
51 {
52  if (_mesh.dimension() > 2)
53  mooseError("This class only supports 2D simulations at this time");
54 
55  _body_forces.push_back(&_body_force_x);
56  _body_forces.push_back(&_body_force_y);
57  _body_forces.push_back(&_body_force_z);
58 }
HDGKernel::_mesh
MooseMesh & _mesh
HDGKernel::HDGKernel
HDGKernel(const InputParameters &parameters)
NavierStokesLHDGKernel::_current_side
const unsigned int & _current_side
The current side when doing face evaluations.
Definition: NavierStokesLHDGKernel.h:63
NavierStokesLHDGKernel::_body_force_x
const Moose::Functor< Real > & _body_force_x
Definition: NavierStokesLHDGKernel.h:44
Assembly::qRuleFace
const libMesh::QBase *const & qRuleFace() const
NavierStokesLHDGKernel::_q_point_face
const MooseArray< Point > & _q_point_face
The physical locations of the quadrature points on the face.
Definition: NavierStokesLHDGKernel.h:54
NavierStokesLHDGKernel::_pressure_mms_forcing_function
const Moose::Functor< Real > & _pressure_mms_forcing_function
Definition: NavierStokesLHDGKernel.h:48
HDGKernel::_tid
THREAD_ID _tid
HDGKernel::_sys
SystemBase & _sys
NavierStokesLHDGKernel::_body_forces
std::vector< const Moose::Functor< Real > * > _body_forces
Definition: NavierStokesLHDGKernel.h:47
NavierStokesLHDGKernel::_JxW_face
const MooseArray< Real > & _JxW_face
transformed Jacobian weights on the current element face
Definition: NavierStokesLHDGKernel.h:57
Assembly::normals
const MooseArray< Point > & normals() const
MooseMesh::dimension
virtual unsigned int dimension() const
NavierStokesLHDGKernel::_normals
const MooseArray< Point > & _normals
face normals
Definition: NavierStokesLHDGKernel.h:60
HDGKernel::_fe_problem
FEProblemBase & _fe_problem
HDGKernel::_assembly
Assembly & _assembly
Assembly::JxWFace
const MooseArray< Real > & JxWFace() const
Assembly::qPointsFace
const MooseArray< Point > & qPointsFace() const
NavierStokesLHDGAssemblyHelper::NavierStokesLHDGAssemblyHelper
NavierStokesLHDGAssemblyHelper(const MooseObject *moose_obj, MaterialPropertyInterface *mpi, MooseVariableDependencyInterface *mvdi, const TransientInterface *const ti, const FEProblemBase &fe_problem, SystemBase &sys, const MooseMesh &mesh, const THREAD_ID tid)
Definition: NavierStokesLHDGAssemblyHelper.C:44
HDGKernel::mooseError
void mooseError(Args &&... args) const
NavierStokesLHDGKernel::_body_force_z
const Moose::Functor< Real > & _body_force_z
Definition: NavierStokesLHDGKernel.h:46
HDGKernel::parameters
const InputParameters & parameters() const
NavierStokesLHDGKernel::_qrule_face
const QBase *const & _qrule_face
The face quadrature rule.
Definition: NavierStokesLHDGKernel.h:51
Assembly::side
const unsigned int & side() const
NavierStokesLHDGKernel::_body_force_y
const Moose::Functor< Real > & _body_force_y
Definition: NavierStokesLHDGKernel.h:45

Member Function Documentation

◆ additionalROVariables()

std::set< std::string > NavierStokesLHDGKernel::additionalROVariables ( )
overridevirtual

Reimplemented from HDGKernel.

Definition at line 342 of file NavierStokesLHDGKernel.C.

343 {
344  std::set<std::string> covered_vars = {_grad_u_var.name(),
345  _u_face_var.name(),
346  _v_var.name(),
347  _grad_v_var.name(),
348  _v_face_var.name(),
349  _pressure_var.name()};
350  if (_enclosure_lm_var)
351  covered_vars.insert(_enclosure_lm_var->name());
352  return covered_vars;
353 }
NavierStokesLHDGAssemblyHelper::_v_var
const MooseVariableFE< Real > & _v_var
Definition: NavierStokesLHDGAssemblyHelper.h:209
NavierStokesLHDGAssemblyHelper::_v_face_var
const MooseVariableFE< Real > & _v_face_var
Definition: NavierStokesLHDGAssemblyHelper.h:213
MooseVariableFE< RealVectorValue >::name
const std::string & name() const override
NavierStokesLHDGAssemblyHelper::_grad_v_var
const MooseVariableFE< RealVectorValue > & _grad_v_var
Definition: NavierStokesLHDGAssemblyHelper.h:211
NavierStokesLHDGAssemblyHelper::_pressure_var
const MooseVariableFE< Real > & _pressure_var
Definition: NavierStokesLHDGAssemblyHelper.h:215
DiffusionLHDGAssemblyHelper::_u_face_var
const MooseVariableFE< Real > & _u_face_var
DiffusionLHDGAssemblyHelper::_grad_u_var
const MooseVariableFE< RealVectorValue > & _grad_u_var
NavierStokesLHDGAssemblyHelper::_enclosure_lm_var
const MooseVariableScalar *const _enclosure_lm_var
Definition: NavierStokesLHDGAssemblyHelper.h:216

◆ computeJacobian()

void NavierStokesLHDGKernel::computeJacobian ( )
overridevirtual

Compute this object's entire element interior Jacobian, both on- and off-diagonal.

Reimplemented from HDGKernel.

Definition at line 95 of file NavierStokesLHDGKernel.C.

Referenced by computeOffDiagJacobian().

96 {
97  _grad_u_grad_u_jac.resize(_qu_dof_indices.size(), _qu_dof_indices.size());
98  _grad_v_grad_v_jac.resize(_qv_dof_indices.size(), _qv_dof_indices.size());
99  _grad_u_u_jac.resize(_qu_dof_indices.size(), _u_dof_indices.size());
100  _grad_v_v_jac.resize(_qv_dof_indices.size(), _v_dof_indices.size());
101  _u_grad_u_jac.resize(_u_dof_indices.size(), _qu_dof_indices.size());
102  _v_grad_v_jac.resize(_v_dof_indices.size(), _qv_dof_indices.size());
103  _u_u_jac.resize(_u_dof_indices.size(), _u_dof_indices.size());
104  _u_v_jac.resize(_u_dof_indices.size(), _v_dof_indices.size());
105  _v_u_jac.resize(_v_dof_indices.size(), _u_dof_indices.size());
106  _v_v_jac.resize(_v_dof_indices.size(), _v_dof_indices.size());
107  _u_p_jac.resize(_u_dof_indices.size(), _p_dof_indices.size());
108  _v_p_jac.resize(_v_dof_indices.size(), _p_dof_indices.size());
109  _p_u_jac.resize(_p_dof_indices.size(), _u_dof_indices.size());
110  _p_v_jac.resize(_p_dof_indices.size(), _v_dof_indices.size());
111  _p_global_lm_jac.resize(_p_dof_indices.size(),
112  _global_lm_dof_indices ? _global_lm_dof_indices->size() : 0);
113  _global_lm_p_jac.resize(_global_lm_dof_indices ? _global_lm_dof_indices->size() : 0,
114  _p_dof_indices.size());
115 
116  // qu and u
117  vectorVolumeJacobian(_JxW, *_qrule, _grad_u_grad_u_jac, _grad_u_u_jac);
118  scalarVolumeJacobian(0, _JxW, *_qrule, _u_grad_u_jac, _u_u_jac, _u_v_jac, _u_p_jac);
119 
120  // qv and v
121  vectorVolumeJacobian(_JxW, *_qrule, _grad_v_grad_v_jac, _grad_v_v_jac);
122  scalarVolumeJacobian(1, _JxW, *_qrule, _v_grad_v_jac, _v_u_jac, _v_v_jac, _v_p_jac);
123 
124  // p
125  pressureVolumeJacobian(_JxW, *_qrule, _p_u_jac, _p_v_jac, _p_global_lm_jac, _global_lm_p_jac);
126 
127  addJacobian(
128  _assembly, _grad_u_grad_u_jac, _qu_dof_indices, _qu_dof_indices, _grad_u_var.scalingFactor());
129  addJacobian(
130  _assembly, _grad_v_grad_v_jac, _qv_dof_indices, _qv_dof_indices, _grad_v_var.scalingFactor());
131  addJacobian(
132  _assembly, _grad_u_u_jac, _qu_dof_indices, _u_dof_indices, _grad_u_var.scalingFactor());
133  addJacobian(
134  _assembly, _grad_v_v_jac, _qv_dof_indices, _v_dof_indices, _grad_v_var.scalingFactor());
135  addJacobian(_assembly, _u_grad_u_jac, _u_dof_indices, _qu_dof_indices, _u_var.scalingFactor());
136  addJacobian(_assembly, _v_grad_v_jac, _v_dof_indices, _qv_dof_indices, _v_var.scalingFactor());
137  addJacobian(_assembly, _u_u_jac, _u_dof_indices, _u_dof_indices, _u_var.scalingFactor());
138  addJacobian(_assembly, _u_v_jac, _u_dof_indices, _v_dof_indices, _u_var.scalingFactor());
139  addJacobian(_assembly, _v_u_jac, _v_dof_indices, _u_dof_indices, _v_var.scalingFactor());
140  addJacobian(_assembly, _v_v_jac, _v_dof_indices, _v_dof_indices, _v_var.scalingFactor());
141  addJacobian(_assembly, _u_p_jac, _u_dof_indices, _p_dof_indices, _u_var.scalingFactor());
142  addJacobian(_assembly, _v_p_jac, _v_dof_indices, _p_dof_indices, _v_var.scalingFactor());
143  addJacobian(_assembly, _p_u_jac, _p_dof_indices, _u_dof_indices, _pressure_var.scalingFactor());
144  addJacobian(_assembly, _p_v_jac, _p_dof_indices, _v_dof_indices, _pressure_var.scalingFactor());
145  if (_global_lm_dof_indices)
146  {
147  addJacobian(_assembly,
148  _p_global_lm_jac,
149  _p_dof_indices,
150  *_global_lm_dof_indices,
151  _pressure_var.scalingFactor());
152  addJacobian(_assembly,
153  _global_lm_p_jac,
154  *_global_lm_dof_indices,
155  _p_dof_indices,
156  _enclosure_lm_var->scalingFactor());
157  }
158 }
NavierStokesLHDGAssemblyHelper::pressureVolumeJacobian
void pressureVolumeJacobian(const MooseArray< Real > &JxW, const libMesh::QBase &qrule, DenseMatrix< Number > &p_u_vel_jac, DenseMatrix< Number > &p_v_vel_jac, DenseMatrix< Number > &p_global_lm_jac, DenseMatrix< Number > &global_lm_p_jac)
Compute the volumetric contributions to the pressure Jacobian, e.g.
Definition: NavierStokesLHDGAssemblyHelper.C:258
NavierStokesLHDGAssemblyHelper::_global_lm_dof_indices
const std::vector< dof_id_type > *const _global_lm_dof_indices
Definition: NavierStokesLHDGAssemblyHelper.h:226
NavierStokesLHDGAssemblyHelper::_v_var
const MooseVariableFE< Real > & _v_var
Definition: NavierStokesLHDGAssemblyHelper.h:209
DiffusionLHDGAssemblyHelper::vectorVolumeJacobian
void vectorVolumeJacobian(const MooseArray< Real > &JxW, const libMesh::QBase &qrule, DenseMatrix< Number > &vector_vector_jac, DenseMatrix< Number > &vector_scalar_jac)
DiffusionLHDGAssemblyHelper::_u_dof_indices
const std::vector< dof_id_type > & _u_dof_indices
HDGKernel::_JxW
const MooseArray< Real > & _JxW
NavierStokesLHDGAssemblyHelper::_qv_dof_indices
const std::vector< dof_id_type > & _qv_dof_indices
Containers for dof indices.
Definition: NavierStokesLHDGAssemblyHelper.h:219
NavierStokesLHDGAssemblyHelper::scalarVolumeJacobian
void scalarVolumeJacobian(const unsigned int vel_component, const MooseArray< Real > &JxW, const libMesh::QBase &qrule, DenseMatrix< Number > &scalar_vector_jac, DenseMatrix< Number > &scalar_u_vel_jac, DenseMatrix< Number > &scalar_v_vel_jac, DenseMatrix< Number > &scalar_p_jac)
Compute the volumetric contributions to a velocity Jacobian.
Definition: NavierStokesLHDGAssemblyHelper.C:178
DiffusionLHDGAssemblyHelper::_qu_dof_indices
const std::vector< dof_id_type > & _qu_dof_indices
HDGKernel::addJacobian
void addJacobian(Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
HDGKernel::_qrule
const QBase *const & _qrule
NavierStokesLHDGAssemblyHelper::_grad_v_var
const MooseVariableFE< RealVectorValue > & _grad_v_var
Definition: NavierStokesLHDGAssemblyHelper.h:211
NavierStokesLHDGAssemblyHelper::_v_dof_indices
const std::vector< dof_id_type > & _v_dof_indices
Definition: NavierStokesLHDGAssemblyHelper.h:220
HDGKernel::_assembly
Assembly & _assembly
NavierStokesLHDGAssemblyHelper::_pressure_var
const MooseVariableFE< Real > & _pressure_var
Definition: NavierStokesLHDGAssemblyHelper.h:215
NavierStokesLHDGAssemblyHelper::_p_dof_indices
const std::vector< dof_id_type > & _p_dof_indices
Definition: NavierStokesLHDGAssemblyHelper.h:225
DiffusionLHDGAssemblyHelper::_u_var
const MooseVariableFE< Real > & _u_var
DenseMatrix< Number >::resize
void resize(const unsigned int new_m, const unsigned int new_n)
DiffusionLHDGAssemblyHelper::_grad_u_var
const MooseVariableFE< RealVectorValue > & _grad_u_var
MooseVariableFE< RealVectorValue >::scalingFactor
void scalingFactor(const std::vector< Real > &factor)
NavierStokesLHDGAssemblyHelper::_enclosure_lm_var
const MooseVariableScalar *const _enclosure_lm_var
Definition: NavierStokesLHDGAssemblyHelper.h:216

◆ computeJacobianOnSide()

void NavierStokesLHDGKernel::computeJacobianOnSide ( )
overridevirtual

Implements HDGKernel.

Definition at line 198 of file NavierStokesLHDGKernel.C.

199 {
200  const Elem * const neigh = _current_elem->neighbor_ptr(_current_side);
201 
202  _grad_u_lm_u_jac.resize(_qu_dof_indices.size(), _lm_u_dof_indices.size());
203  _u_grad_u_jac.resize(_u_dof_indices.size(), _qu_dof_indices.size());
204  _u_u_jac.resize(_u_dof_indices.size(), _u_dof_indices.size());
205  _u_lm_u_jac.resize(_u_dof_indices.size(), _lm_u_dof_indices.size());
206  _u_lm_v_jac.resize(_u_dof_indices.size(), _lm_v_dof_indices.size());
207  _u_p_jac.resize(_u_dof_indices.size(), _p_dof_indices.size());
208  _lm_u_grad_u_jac.resize(_lm_u_dof_indices.size(), _qu_dof_indices.size());
209  _lm_u_u_jac.resize(_lm_u_dof_indices.size(), _u_dof_indices.size());
210  _lm_u_lm_u_jac.resize(_lm_u_dof_indices.size(), _lm_u_dof_indices.size());
211  _lm_u_lm_v_jac.resize(_lm_u_dof_indices.size(), _lm_v_dof_indices.size());
212  _lm_u_p_jac.resize(_lm_u_dof_indices.size(), _p_dof_indices.size());
213  _grad_v_lm_v_jac.resize(_qv_dof_indices.size(), _lm_v_dof_indices.size());
214  _v_grad_v_jac.resize(_v_dof_indices.size(), _qv_dof_indices.size());
215  _v_v_jac.resize(_v_dof_indices.size(), _v_dof_indices.size());
216  _v_lm_u_jac.resize(_v_dof_indices.size(), _lm_u_dof_indices.size());
217  _v_lm_v_jac.resize(_v_dof_indices.size(), _lm_v_dof_indices.size());
218  _v_p_jac.resize(_v_dof_indices.size(), _p_dof_indices.size());
219  _lm_v_grad_v_jac.resize(_lm_v_dof_indices.size(), _qv_dof_indices.size());
220  _lm_v_v_jac.resize(_lm_v_dof_indices.size(), _v_dof_indices.size());
221  _lm_v_lm_u_jac.resize(_lm_v_dof_indices.size(), _lm_u_dof_indices.size());
222  _lm_v_lm_v_jac.resize(_lm_v_dof_indices.size(), _lm_v_dof_indices.size());
223  _lm_v_p_jac.resize(_lm_v_dof_indices.size(), _p_dof_indices.size());
224  _p_lm_u_jac.resize(_p_dof_indices.size(), _lm_u_dof_indices.size());
225  _p_lm_v_jac.resize(_p_dof_indices.size(), _lm_v_dof_indices.size());
226 
227  // qu, u, lm_u
228  vectorFaceJacobian(_JxW_face, *_qrule_face, _normals, _grad_u_lm_u_jac);
229  scalarFaceJacobian(0,
230  _JxW_face,
231  *_qrule_face,
232  _normals,
233  _u_grad_u_jac,
234  _u_u_jac,
235  _u_lm_u_jac,
236  _u_p_jac,
237  _u_lm_u_jac,
238  _u_lm_v_jac);
239  lmFaceJacobian(0,
240  _JxW_face,
241  *_qrule_face,
242  _normals,
243  neigh,
244  _lm_u_grad_u_jac,
245  _lm_u_u_jac,
246  _lm_u_lm_u_jac,
247  _lm_u_p_jac,
248  _lm_u_lm_u_jac,
249  _lm_u_lm_v_jac);
250 
251  // qv, v, lm_v
252  vectorFaceJacobian(_JxW_face, *_qrule_face, _normals, _grad_v_lm_v_jac);
253  scalarFaceJacobian(1,
254  _JxW_face,
255  *_qrule_face,
256  _normals,
257  _v_grad_v_jac,
258  _v_v_jac,
259  _v_lm_v_jac,
260  _v_p_jac,
261  _v_lm_u_jac,
262  _v_lm_v_jac);
263  lmFaceJacobian(1,
264  _JxW_face,
265  *_qrule_face,
266  _normals,
267  neigh,
268  _lm_v_grad_v_jac,
269  _lm_v_v_jac,
270  _lm_v_lm_v_jac,
271  _lm_v_p_jac,
272  _lm_v_lm_u_jac,
273  _lm_v_lm_v_jac);
274 
275  // p
276  pressureFaceJacobian(_JxW_face, *_qrule_face, _normals, _p_lm_u_jac, _p_lm_v_jac);
277 
278  addJacobian(
279  _assembly, _grad_u_lm_u_jac, _qu_dof_indices, _lm_u_dof_indices, _grad_u_var.scalingFactor());
280  addJacobian(_assembly, _u_grad_u_jac, _u_dof_indices, _qu_dof_indices, _u_var.scalingFactor());
281  addJacobian(_assembly, _u_u_jac, _u_dof_indices, _u_dof_indices, _u_var.scalingFactor());
282  addJacobian(_assembly, _u_lm_u_jac, _u_dof_indices, _lm_u_dof_indices, _u_var.scalingFactor());
283  addJacobian(_assembly, _u_lm_v_jac, _u_dof_indices, _lm_v_dof_indices, _u_var.scalingFactor());
284  addJacobian(_assembly, _u_p_jac, _u_dof_indices, _p_dof_indices, _u_var.scalingFactor());
285  addJacobian(
286  _assembly, _lm_u_grad_u_jac, _lm_u_dof_indices, _qu_dof_indices, _u_face_var.scalingFactor());
287  addJacobian(
288  _assembly, _lm_u_u_jac, _lm_u_dof_indices, _u_dof_indices, _u_face_var.scalingFactor());
289  addJacobian(
290  _assembly, _lm_u_lm_u_jac, _lm_u_dof_indices, _lm_u_dof_indices, _u_face_var.scalingFactor());
291  addJacobian(
292  _assembly, _lm_u_lm_v_jac, _lm_u_dof_indices, _lm_v_dof_indices, _u_face_var.scalingFactor());
293  addJacobian(
294  _assembly, _lm_u_p_jac, _lm_u_dof_indices, _p_dof_indices, _u_face_var.scalingFactor());
295  addJacobian(
296  _assembly, _grad_v_lm_v_jac, _qv_dof_indices, _lm_v_dof_indices, _grad_v_var.scalingFactor());
297  addJacobian(_assembly, _v_grad_v_jac, _v_dof_indices, _qv_dof_indices, _v_var.scalingFactor());
298  addJacobian(_assembly, _v_v_jac, _v_dof_indices, _v_dof_indices, _v_var.scalingFactor());
299  addJacobian(_assembly, _v_lm_u_jac, _v_dof_indices, _lm_u_dof_indices, _v_var.scalingFactor());
300  addJacobian(_assembly, _v_lm_v_jac, _v_dof_indices, _lm_v_dof_indices, _v_var.scalingFactor());
301  addJacobian(_assembly, _v_p_jac, _v_dof_indices, _p_dof_indices, _v_var.scalingFactor());
302  addJacobian(
303  _assembly, _lm_v_grad_v_jac, _lm_v_dof_indices, _qv_dof_indices, _v_face_var.scalingFactor());
304  addJacobian(
305  _assembly, _lm_v_v_jac, _lm_v_dof_indices, _v_dof_indices, _v_face_var.scalingFactor());
306  addJacobian(
307  _assembly, _lm_v_lm_u_jac, _lm_v_dof_indices, _lm_u_dof_indices, _v_face_var.scalingFactor());
308  addJacobian(
309  _assembly, _lm_v_lm_v_jac, _lm_v_dof_indices, _lm_v_dof_indices, _v_face_var.scalingFactor());
310  addJacobian(
311  _assembly, _lm_v_p_jac, _lm_v_dof_indices, _p_dof_indices, _v_face_var.scalingFactor());
312  addJacobian(
313  _assembly, _p_lm_u_jac, _p_dof_indices, _lm_u_dof_indices, _pressure_var.scalingFactor());
314  addJacobian(
315  _assembly, _p_lm_v_jac, _p_dof_indices, _lm_v_dof_indices, _pressure_var.scalingFactor());
316 }
DiffusionLHDGAssemblyHelper::_lm_u_dof_indices
const std::vector< dof_id_type > & _lm_u_dof_indices
NavierStokesLHDGAssemblyHelper::lmFaceJacobian
void lmFaceJacobian(const unsigned int vel_component, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, const Elem *const neigh, DenseMatrix< Number > &lm_vec_jac, DenseMatrix< Number > &lm_scalar_jac, DenseMatrix< Number > &lm_lm_jac, DenseMatrix< Number > &lm_p_jac, DenseMatrix< Number > &lm_lm_u_vel_jac, DenseMatrix< Number > &lm_lm_v_vel_jac)
Definition: NavierStokesLHDGAssemblyHelper.C:478
NavierStokesLHDGAssemblyHelper::_v_var
const MooseVariableFE< Real > & _v_var
Definition: NavierStokesLHDGAssemblyHelper.h:209
NavierStokesLHDGKernel::_current_side
const unsigned int & _current_side
The current side when doing face evaluations.
Definition: NavierStokesLHDGKernel.h:63
DiffusionLHDGAssemblyHelper::_u_dof_indices
const std::vector< dof_id_type > & _u_dof_indices
NavierStokesLHDGAssemblyHelper::_v_face_var
const MooseVariableFE< Real > & _v_face_var
Definition: NavierStokesLHDGAssemblyHelper.h:213
NavierStokesLHDGAssemblyHelper::pressureFaceJacobian
void pressureFaceJacobian(const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseMatrix< Number > &p_lm_u_vel_jac, DenseMatrix< Number > &p_lm_v_vel_jac)
Definition: NavierStokesLHDGAssemblyHelper.C:337
NavierStokesLHDGAssemblyHelper::_qv_dof_indices
const std::vector< dof_id_type > & _qv_dof_indices
Containers for dof indices.
Definition: NavierStokesLHDGAssemblyHelper.h:219
DiffusionLHDGAssemblyHelper::_qu_dof_indices
const std::vector< dof_id_type > & _qu_dof_indices
HDGKernel::addJacobian
void addJacobian(Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
NavierStokesLHDGAssemblyHelper::scalarFaceJacobian
void scalarFaceJacobian(const unsigned int vel_component, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseMatrix< Number > &scalar_vector_jac, DenseMatrix< Number > &scalar_scalar_jac, DenseMatrix< Number > &scalar_lm_jac, DenseMatrix< Number > &scalar_p_jac, DenseMatrix< Number > &scalar_lm_u_vel_jac, DenseMatrix< Number > &scalar_lm_v_vel_jac)
Definition: NavierStokesLHDGAssemblyHelper.C:393
NavierStokesLHDGKernel::_JxW_face
const MooseArray< Real > & _JxW_face
transformed Jacobian weights on the current element face
Definition: NavierStokesLHDGKernel.h:57
NavierStokesLHDGKernel::_normals
const MooseArray< Point > & _normals
face normals
Definition: NavierStokesLHDGKernel.h:60
NavierStokesLHDGAssemblyHelper::_grad_v_var
const MooseVariableFE< RealVectorValue > & _grad_v_var
Definition: NavierStokesLHDGAssemblyHelper.h:211
NavierStokesLHDGAssemblyHelper::_lm_v_dof_indices
const std::vector< dof_id_type > & _lm_v_dof_indices
Definition: NavierStokesLHDGAssemblyHelper.h:221
NavierStokesLHDGAssemblyHelper::_v_dof_indices
const std::vector< dof_id_type > & _v_dof_indices
Definition: NavierStokesLHDGAssemblyHelper.h:220
HDGKernel::_assembly
Assembly & _assembly
DiffusionLHDGAssemblyHelper::vectorFaceJacobian
void vectorFaceJacobian(const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseMatrix< Number > &vector_lm_jac)
NavierStokesLHDGAssemblyHelper::_pressure_var
const MooseVariableFE< Real > & _pressure_var
Definition: NavierStokesLHDGAssemblyHelper.h:215
DiffusionLHDGAssemblyHelper::_u_face_var
const MooseVariableFE< Real > & _u_face_var
NavierStokesLHDGAssemblyHelper::_p_dof_indices
const std::vector< dof_id_type > & _p_dof_indices
Definition: NavierStokesLHDGAssemblyHelper.h:225
DiffusionLHDGAssemblyHelper::_u_var
const MooseVariableFE< Real > & _u_var
DenseMatrix< Number >::resize
void resize(const unsigned int new_m, const unsigned int new_n)
HDGKernel::_current_elem
const Elem *const & _current_elem
NavierStokesLHDGKernel::_qrule_face
const QBase *const & _qrule_face
The face quadrature rule.
Definition: NavierStokesLHDGKernel.h:51
DiffusionLHDGAssemblyHelper::_grad_u_var
const MooseVariableFE< RealVectorValue > & _grad_u_var
MooseVariableFE< RealVectorValue >::scalingFactor
void scalingFactor(const std::vector< Real > &factor)

◆ computeOffDiagJacobian()

void NavierStokesLHDGKernel::computeOffDiagJacobian ( unsigned int  jvar)
overridevirtual

Forwards to computeJacobian() the first time this is called for a given element.

Reimplemented from HDGKernel.

Definition at line 332 of file NavierStokesLHDGKernel.C.

333 {
334  if (_cached_elem != _current_elem)
335  {
336  computeJacobian();
337  _cached_elem = _current_elem;
338  }
339 }
NavierStokesLHDGKernel::computeJacobian
virtual void computeJacobian() override
Compute this object&#39;s entire element interior Jacobian, both on- and off-diagonal.
Definition: NavierStokesLHDGKernel.C:95
HDGKernel::_current_elem
const Elem *const & _current_elem

◆ computeResidual()

void NavierStokesLHDGKernel::computeResidual ( )
overridevirtual

Reimplemented from HDGKernel.

Definition at line 61 of file NavierStokesLHDGKernel.C.

62 {
63  _grad_u_vel_re.resize(_qu_dof_indices.size());
64  _grad_v_vel_re.resize(_qv_dof_indices.size());
65  _u_vel_re.resize(_u_dof_indices.size());
66  _v_vel_re.resize(_v_dof_indices.size());
67  _p_re.resize(_p_dof_indices.size());
68  _global_lm_re.resize(_global_lm_dof_indices ? _global_lm_dof_indices->size() : 0);
69 
70  // qu and u
71  vectorVolumeResidual(_qu_sol, _u_sol, _JxW, *_qrule, _grad_u_vel_re);
72  scalarVolumeResidual(
73  _qu_sol, 0, _body_force_x, _JxW, *_qrule, _current_elem, _q_point, _u_vel_re);
74 
75  // qv and v
76  vectorVolumeResidual(_qv_sol, _v_sol, _JxW, *_qrule, _grad_v_vel_re);
77  scalarVolumeResidual(
78  _qv_sol, 1, _body_force_y, _JxW, *_qrule, _current_elem, _q_point, _v_vel_re);
79 
80  // p
81  pressureVolumeResidual(
82  _pressure_mms_forcing_function, _JxW, *_qrule, _current_elem, _q_point, _p_re, _global_lm_re);
83 
84  addResiduals(_assembly, _grad_u_vel_re, _qu_dof_indices, _grad_u_var.scalingFactor());
85  addResiduals(_assembly, _grad_v_vel_re, _qv_dof_indices, _grad_v_var.scalingFactor());
86  addResiduals(_assembly, _u_vel_re, _u_dof_indices, _u_var.scalingFactor());
87  addResiduals(_assembly, _v_vel_re, _v_dof_indices, _v_var.scalingFactor());
88  addResiduals(_assembly, _p_re, _p_dof_indices, _pressure_var.scalingFactor());
89  if (_global_lm_dof_indices)
90  addResiduals(
91  _assembly, _global_lm_re, *_global_lm_dof_indices, _enclosure_lm_var->scalingFactor());
92 }
DiffusionLHDGAssemblyHelper::_u_sol
const MooseArray< Number > & _u_sol
NavierStokesLHDGAssemblyHelper::_global_lm_dof_indices
const std::vector< dof_id_type > *const _global_lm_dof_indices
Definition: NavierStokesLHDGAssemblyHelper.h:226
NavierStokesLHDGAssemblyHelper::_v_var
const MooseVariableFE< Real > & _v_var
Definition: NavierStokesLHDGAssemblyHelper.h:209
NavierStokesLHDGKernel::_body_force_x
const Moose::Functor< Real > & _body_force_x
Definition: NavierStokesLHDGKernel.h:44
DiffusionLHDGAssemblyHelper::_u_dof_indices
const std::vector< dof_id_type > & _u_dof_indices
HDGKernel::_JxW
const MooseArray< Real > & _JxW
NavierStokesLHDGAssemblyHelper::pressureVolumeResidual
void pressureVolumeResidual(const Moose::Functor< Real > &pressure_mms_forcing_function, const MooseArray< Real > &JxW, const libMesh::QBase &qrule, const Elem *const current_elem, const MooseArray< Point > &q_point, DenseVector< Number > &pressure_re, DenseVector< Number > &global_lm_re)
Compute the volumetric contributions to the pressure residual, e.g.
Definition: NavierStokesLHDGAssemblyHelper.C:220
NavierStokesLHDGKernel::_pressure_mms_forcing_function
const Moose::Functor< Real > & _pressure_mms_forcing_function
Definition: NavierStokesLHDGKernel.h:48
DenseVector< Number >::resize
void resize(const unsigned int n)
HDGKernel::addResiduals
void addResiduals(Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
DiffusionLHDGAssemblyHelper::_qu_sol
const MooseArray< libMesh::Gradient > & _qu_sol
NavierStokesLHDGAssemblyHelper::_qv_dof_indices
const std::vector< dof_id_type > & _qv_dof_indices
Containers for dof indices.
Definition: NavierStokesLHDGAssemblyHelper.h:219
DiffusionLHDGAssemblyHelper::_qu_dof_indices
const std::vector< dof_id_type > & _qu_dof_indices
DiffusionLHDGAssemblyHelper::vectorVolumeResidual
void vectorVolumeResidual(const MooseArray< Gradient > &vector_sol, const MooseArray< Number > &scalar_sol, const MooseArray< Real > &JxW, const libMesh::QBase &qrule, DenseVector< Number > &vector_re)
HDGKernel::_qrule
const QBase *const & _qrule
NavierStokesLHDGAssemblyHelper::_grad_v_var
const MooseVariableFE< RealVectorValue > & _grad_v_var
Definition: NavierStokesLHDGAssemblyHelper.h:211
NavierStokesLHDGAssemblyHelper::_v_dof_indices
const std::vector< dof_id_type > & _v_dof_indices
Definition: NavierStokesLHDGAssemblyHelper.h:220
HDGKernel::_assembly
Assembly & _assembly
NavierStokesLHDGAssemblyHelper::_pressure_var
const MooseVariableFE< Real > & _pressure_var
Definition: NavierStokesLHDGAssemblyHelper.h:215
NavierStokesLHDGAssemblyHelper::_p_dof_indices
const std::vector< dof_id_type > & _p_dof_indices
Definition: NavierStokesLHDGAssemblyHelper.h:225
DiffusionLHDGAssemblyHelper::_u_var
const MooseVariableFE< Real > & _u_var
NavierStokesLHDGAssemblyHelper::_qv_sol
const MooseArray< Gradient > & _qv_sol
local solutions at quadrature points
Definition: NavierStokesLHDGAssemblyHelper.h:229
HDGKernel::_current_elem
const Elem *const & _current_elem
DiffusionLHDGAssemblyHelper::_grad_u_var
const MooseVariableFE< RealVectorValue > & _grad_u_var
NavierStokesLHDGAssemblyHelper::scalarVolumeResidual
void scalarVolumeResidual(const MooseArray< Gradient > &vel_gradient, const unsigned int vel_component, const Moose::Functor< Real > &body_force, const MooseArray< Real > &JxW, const libMesh::QBase &qrule, const Elem *const current_elem, const MooseArray< Point > &q_point, DenseVector< Number > &scalar_re)
Compute the volumetric contributions to a velocity residual for a provided velocity gradient and stre...
Definition: NavierStokesLHDGAssemblyHelper.C:148
HDGKernel::_q_point
const MooseArray< Point > & _q_point
MooseVariableFE< RealVectorValue >::scalingFactor
void scalingFactor(const std::vector< Real > &factor)
NavierStokesLHDGAssemblyHelper::_enclosure_lm_var
const MooseVariableScalar *const _enclosure_lm_var
Definition: NavierStokesLHDGAssemblyHelper.h:216
NavierStokesLHDGKernel::_body_force_y
const Moose::Functor< Real > & _body_force_y
Definition: NavierStokesLHDGKernel.h:45

◆ computeResidualOnSide()

void NavierStokesLHDGKernel::computeResidualOnSide ( )
overridevirtual

Implements HDGKernel.

Definition at line 161 of file NavierStokesLHDGKernel.C.

162 {
163  const Elem * const neigh = _current_elem->neighbor_ptr(_current_side);
164 
165  _grad_u_vel_re.resize(_qu_dof_indices.size());
166  _u_vel_re.resize(_u_dof_indices.size());
167  _lm_u_vel_re.resize(_lm_u_dof_indices.size());
168  _grad_v_vel_re.resize(_qv_dof_indices.size());
169  _v_vel_re.resize(_v_dof_indices.size());
170  _lm_v_vel_re.resize(_lm_v_dof_indices.size());
171  _p_re.resize(_p_dof_indices.size());
172 
173  // qu, u, lm_u
174  vectorFaceResidual(_lm_u_sol, _JxW_face, *_qrule_face, _normals, _grad_u_vel_re);
175  scalarFaceResidual(_qu_sol, _u_sol, _lm_u_sol, 0, _JxW_face, *_qrule_face, _normals, _u_vel_re);
176  lmFaceResidual(
177  _qu_sol, _u_sol, _lm_u_sol, 0, _JxW_face, *_qrule_face, _normals, neigh, _lm_u_vel_re);
178 
179  // qv, v, lm_v
180  vectorFaceResidual(_lm_v_sol, _JxW_face, *_qrule_face, _normals, _grad_v_vel_re);
181  scalarFaceResidual(_qv_sol, _v_sol, _lm_v_sol, 1, _JxW_face, *_qrule_face, _normals, _v_vel_re);
182  lmFaceResidual(
183  _qv_sol, _v_sol, _lm_v_sol, 1, _JxW_face, *_qrule_face, _normals, neigh, _lm_v_vel_re);
184 
185  // p
186  pressureFaceResidual(_JxW_face, *_qrule_face, _normals, _p_re);
187 
188  addResiduals(_assembly, _grad_u_vel_re, _qu_dof_indices, _grad_u_var.scalingFactor());
189  addResiduals(_assembly, _u_vel_re, _u_dof_indices, _u_var.scalingFactor());
190  addResiduals(_assembly, _lm_u_vel_re, _lm_u_dof_indices, _u_face_var.scalingFactor());
191  addResiduals(_assembly, _grad_v_vel_re, _qv_dof_indices, _grad_v_var.scalingFactor());
192  addResiduals(_assembly, _v_vel_re, _v_dof_indices, _v_var.scalingFactor());
193  addResiduals(_assembly, _lm_v_vel_re, _lm_v_dof_indices, _v_face_var.scalingFactor());
194  addResiduals(_assembly, _p_re, _p_dof_indices, _pressure_var.scalingFactor());
195 }
DiffusionLHDGAssemblyHelper::_lm_u_dof_indices
const std::vector< dof_id_type > & _lm_u_dof_indices
DiffusionLHDGAssemblyHelper::_u_sol
const MooseArray< Number > & _u_sol
NavierStokesLHDGAssemblyHelper::_v_var
const MooseVariableFE< Real > & _v_var
Definition: NavierStokesLHDGAssemblyHelper.h:209
NavierStokesLHDGKernel::_current_side
const unsigned int & _current_side
The current side when doing face evaluations.
Definition: NavierStokesLHDGKernel.h:63
DiffusionLHDGAssemblyHelper::_u_dof_indices
const std::vector< dof_id_type > & _u_dof_indices
DiffusionLHDGAssemblyHelper::vectorFaceResidual
void vectorFaceResidual(const MooseArray< Number > &lm_sol, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseVector< Number > &vector_re)
DenseVector< Number >::resize
void resize(const unsigned int n)
HDGKernel::addResiduals
void addResiduals(Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
NavierStokesLHDGAssemblyHelper::_v_face_var
const MooseVariableFE< Real > & _v_face_var
Definition: NavierStokesLHDGAssemblyHelper.h:213
DiffusionLHDGAssemblyHelper::_qu_sol
const MooseArray< libMesh::Gradient > & _qu_sol
NavierStokesLHDGAssemblyHelper::scalarFaceResidual
void scalarFaceResidual(const MooseArray< Gradient > &vector_sol, const MooseArray< Number > &scalar_sol, const MooseArray< Number > &lm_sol, const unsigned int vel_component, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseVector< Number > &scalar_re)
Definition: NavierStokesLHDGAssemblyHelper.C:363
NavierStokesLHDGAssemblyHelper::_qv_dof_indices
const std::vector< dof_id_type > & _qv_dof_indices
Containers for dof indices.
Definition: NavierStokesLHDGAssemblyHelper.h:219
DiffusionLHDGAssemblyHelper::_qu_dof_indices
const std::vector< dof_id_type > & _qu_dof_indices
NavierStokesLHDGKernel::_JxW_face
const MooseArray< Real > & _JxW_face
transformed Jacobian weights on the current element face
Definition: NavierStokesLHDGKernel.h:57
NavierStokesLHDGKernel::_normals
const MooseArray< Point > & _normals
face normals
Definition: NavierStokesLHDGKernel.h:60
NavierStokesLHDGAssemblyHelper::_grad_v_var
const MooseVariableFE< RealVectorValue > & _grad_v_var
Definition: NavierStokesLHDGAssemblyHelper.h:211
NavierStokesLHDGAssemblyHelper::_lm_v_dof_indices
const std::vector< dof_id_type > & _lm_v_dof_indices
Definition: NavierStokesLHDGAssemblyHelper.h:221
NavierStokesLHDGAssemblyHelper::_v_dof_indices
const std::vector< dof_id_type > & _v_dof_indices
Definition: NavierStokesLHDGAssemblyHelper.h:220
HDGKernel::_assembly
Assembly & _assembly
NavierStokesLHDGAssemblyHelper::_pressure_var
const MooseVariableFE< Real > & _pressure_var
Definition: NavierStokesLHDGAssemblyHelper.h:215
DiffusionLHDGAssemblyHelper::_u_face_var
const MooseVariableFE< Real > & _u_face_var
NavierStokesLHDGAssemblyHelper::_p_dof_indices
const std::vector< dof_id_type > & _p_dof_indices
Definition: NavierStokesLHDGAssemblyHelper.h:225
NavierStokesLHDGAssemblyHelper::pressureFaceResidual
void pressureFaceResidual(const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseVector< Number > &pressure_re)
Definition: NavierStokesLHDGAssemblyHelper.C:322
DiffusionLHDGAssemblyHelper::_u_var
const MooseVariableFE< Real > & _u_var
DiffusionLHDGAssemblyHelper::_lm_u_sol
const MooseArray< Number > & _lm_u_sol
NavierStokesLHDGAssemblyHelper::_qv_sol
const MooseArray< Gradient > & _qv_sol
local solutions at quadrature points
Definition: NavierStokesLHDGAssemblyHelper.h:229
HDGKernel::_current_elem
const Elem *const & _current_elem
NavierStokesLHDGKernel::_qrule_face
const QBase *const & _qrule_face
The face quadrature rule.
Definition: NavierStokesLHDGKernel.h:51
DiffusionLHDGAssemblyHelper::_grad_u_var
const MooseVariableFE< RealVectorValue > & _grad_u_var
NavierStokesLHDGAssemblyHelper::lmFaceResidual
void lmFaceResidual(const MooseArray< Gradient > &vector_sol, const MooseArray< Number > &scalar_sol, const MooseArray< Number > &lm_sol, const unsigned int vel_component, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, const Elem *const neigh, DenseVector< Number > &lm_re)
Definition: NavierStokesLHDGAssemblyHelper.C:444
MooseVariableFE< RealVectorValue >::scalingFactor
void scalingFactor(const std::vector< Real > &factor)

◆ initialSetup()

void NavierStokesLHDGKernel::initialSetup ( )
overridevirtual

Reimplemented from HDGKernel.

Definition at line 319 of file NavierStokesLHDGKernel.C.

320 {
321  // This check must occur after FEProblemBase::init()
322  checkCoupling();
323 }

◆ jacobianSetup()

void NavierStokesLHDGKernel::jacobianSetup ( )
overridevirtual

Reimplemented from HDGKernel.

Definition at line 326 of file NavierStokesLHDGKernel.C.

327 {
328  _cached_elem = nullptr;
329 }

◆ lmFaceJacobian()

void NavierStokesLHDGAssemblyHelper::lmFaceJacobian ( const unsigned int  vel_component,
const MooseArray< Real > &  JxW_face,
const libMesh::QBase &  qrule_face,
const MooseArray< Point > &  normals,
const Elem *const  neigh,
DenseMatrix< Number > &  lm_vec_jac,
DenseMatrix< Number > &  lm_scalar_jac,
DenseMatrix< Number > &  lm_lm_jac,
DenseMatrix< Number > &  lm_p_jac,
DenseMatrix< Number > &  lm_lm_u_vel_jac,
DenseMatrix< Number > &  lm_lm_v_vel_jac 
)
protectedinherited

Definition at line 478 of file NavierStokesLHDGAssemblyHelper.C.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

489 {
490  DiffusionLHDGAssemblyHelper::lmFaceJacobian(
491  JxW_face, qrule_face, normals, lm_vec_jac, lm_scalar_jac, lm_lm_jac);
492 
493  for (const auto i : make_range(lm_p_jac.m()))
494  for (const auto qp : make_range(qrule_face.n_points()))
495  {
496  for (const auto j : make_range(lm_p_jac.n()))
497  {
498  Gradient p_phi;
499  p_phi(vel_component) = _scalar_phi_face[j][qp];
500  lm_p_jac(i, j) += JxW_face[qp] * _lm_phi_face[i][qp] * (p_phi * normals[qp]);
501  }
502 
503  for (const auto j : make_range(lm_lm_u_vel_jac.n()))
504  if (neigh)
505  {
506  // derivatives wrt 0th component of velocity
507  {
508  const auto rho_vel_cross_vel =
509  rhoVelCrossVelJacobian(_lm_u_sol, _lm_v_sol, qp, vel_component, 0, _lm_phi_face, j);
510  lm_lm_u_vel_jac(i, j) +=
511  JxW_face[qp] * _lm_phi_face[i][qp] * rho_vel_cross_vel * normals[qp];
512  }
513  // derivatives wrt 1th component of velocity
514  {
515  const auto rho_vel_cross_vel =
516  rhoVelCrossVelJacobian(_lm_u_sol, _lm_v_sol, qp, vel_component, 1, _lm_phi_face, j);
517  lm_lm_v_vel_jac(i, j) +=
518  JxW_face[qp] * _lm_phi_face[i][qp] * rho_vel_cross_vel * normals[qp];
519  }
520  }
521  }
522 }
DiffusionLHDGAssemblyHelper::_scalar_phi_face
const MooseArray< std::vector< Real > > & _scalar_phi_face
NavierStokesLHDGAssemblyHelper::rhoVelCrossVelJacobian
RealVectorValue rhoVelCrossVelJacobian(const MooseArray< Number > &u_sol, const MooseArray< Number > &v_sol, const unsigned int qp, const unsigned int vel_component, const unsigned int vel_j_component, const MooseArray< std::vector< Real >> &phi, const unsigned int j)
Definition: NavierStokesLHDGAssemblyHelper.C:303
DiffusionLHDGAssemblyHelper::_lm_phi_face
const MooseArray< std::vector< Real > > & _lm_phi_face
libMesh::DenseMatrix::m
unsigned int m() const
DiffusionLHDGAssemblyHelper::lmFaceJacobian
void lmFaceJacobian(const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseMatrix< Number > &lm_vec_jac, DenseMatrix< Number > &lm_scalar_jac, DenseMatrix< Number > &lm_lm_jac)
libMesh::QBase::n_points
unsigned int n_points() const
DiffusionLHDGAssemblyHelper::_lm_u_sol
const MooseArray< Number > & _lm_u_sol
make_range
IntRange< T > make_range(T beg, T end)
EM::j
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
libMesh::DenseMatrix::n
unsigned int n() const

◆ lmFaceResidual()

void NavierStokesLHDGAssemblyHelper::lmFaceResidual ( const MooseArray< Gradient > &  vector_sol,
const MooseArray< Number > &  scalar_sol,
const MooseArray< Number > &  lm_sol,
const unsigned int  vel_component,
const MooseArray< Real > &  JxW_face,
const libMesh::QBase &  qrule_face,
const MooseArray< Point > &  normals,
const Elem *const  neigh,
DenseVector< Number > &  lm_re 
)
protectedinherited

Definition at line 444 of file NavierStokesLHDGAssemblyHelper.C.

Referenced by NavierStokesLHDGOutflowBC::computeResidual(), and computeResidualOnSide().

453 {
454  DiffusionLHDGAssemblyHelper::lmFaceResidual(
455  vector_sol, scalar_sol, lm_sol, JxW_face, qrule_face, normals, lm_re);
456 
457  for (const auto qp : make_range(qrule_face.n_points()))
458  {
459  Gradient qp_p;
460  qp_p(vel_component) = _p_sol[qp];
461  const auto rho_vel_cross_vel = rhoVelCrossVelResidual(_lm_u_sol, _lm_v_sol, qp, vel_component);
462 
463  for (const auto i : make_range(lm_re.size()))
464  {
465  // pressure
466  lm_re(i) += JxW_face[qp] * _lm_phi_face[i][qp] * (qp_p * normals[qp]);
467 
468  // If we are an internal face we add the convective term. On the outflow boundary we do not
469  // zero out the convection term, e.g. we are going to set q + p + tau * (u - u_hat) to zero
470  if (neigh)
471  // lm from convection term
472  lm_re(i) += JxW_face[qp] * _lm_phi_face[i][qp] * rho_vel_cross_vel * normals[qp];
473  }
474  }
475 }
DiffusionLHDGAssemblyHelper::_lm_phi_face
const MooseArray< std::vector< Real > > & _lm_phi_face
NavierStokesLHDGAssemblyHelper::rhoVelCrossVelResidual
RealVectorValue rhoVelCrossVelResidual(const MooseArray< Number > &u_sol, const MooseArray< Number > &v_sol, const unsigned int qp, const unsigned int vel_component)
Definition: NavierStokesLHDGAssemblyHelper.C:293
DiffusionLHDGAssemblyHelper::lmFaceResidual
void lmFaceResidual(const MooseArray< Gradient > &vector_sol, const MooseArray< Number > &scalar_sol, const MooseArray< Number > &lm_sol, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseVector< Number > &lm_re)
libMesh::QBase::n_points
unsigned int n_points() const
DiffusionLHDGAssemblyHelper::_lm_u_sol
const MooseArray< Number > & _lm_u_sol
make_range
IntRange< T > make_range(T beg, T end)
libMesh::DenseVector::size
virtual unsigned int size() const override final

◆ pressureDirichletResidual()

void NavierStokesLHDGAssemblyHelper::pressureDirichletResidual ( const std::array< const Moose::Functor< Real > *, 3 > &  dirichlet_vel,
const MooseArray< Real > &  JxW_face,
const libMesh::QBase &  qrule_face,
const MooseArray< Point > &  normals,
const Elem *const  current_elem,
const unsigned int  current_side,
const MooseArray< Point > &  q_point_face,
DenseVector< Number > &  pressure_re 
)
protectedinherited

Definition at line 525 of file NavierStokesLHDGAssemblyHelper.C.

Referenced by NavierStokesLHDGVelocityDirichletBC::computeResidual().

534 {
535  for (const auto qp : make_range(qrule_face.n_points()))
536  {
537  const Moose::ElemSideQpArg elem_side_qp_arg{
538  current_elem, current_side, qp, &qrule_face, q_point_face[qp]};
539  const auto time_arg = _ti.determineState();
540  const RealVectorValue dirichlet_velocity((*dirichlet_vel[0])(elem_side_qp_arg, time_arg),
541  (*dirichlet_vel[1])(elem_side_qp_arg, time_arg),
542  (*dirichlet_vel[2])(elem_side_qp_arg, time_arg));
543  const auto vdotn = dirichlet_velocity * normals[qp];
544  for (const auto i : make_range(pressure_re.size()))
545  pressure_re(i) += JxW_face[qp] * vdotn * _scalar_phi_face[i][qp];
546  }
547 }
DiffusionLHDGAssemblyHelper::_scalar_phi_face
const MooseArray< std::vector< Real > > & _scalar_phi_face
TransientInterface::determineState
Moose::StateArg determineState() const
libMesh::QBase::n_points
unsigned int n_points() const
DiffusionLHDGAssemblyHelper::_ti
const TransientInterface & _ti
make_range
IntRange< T > make_range(T beg, T end)
libMesh::DenseVector::size
virtual unsigned int size() const override final

◆ pressureFaceJacobian()

void NavierStokesLHDGAssemblyHelper::pressureFaceJacobian ( const MooseArray< Real > &  JxW_face,
const libMesh::QBase &  qrule_face,
const MooseArray< Point > &  normals,
DenseMatrix< Number > &  p_lm_u_vel_jac,
DenseMatrix< Number > &  p_lm_v_vel_jac 
)
protectedinherited

Definition at line 337 of file NavierStokesLHDGAssemblyHelper.C.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

342 {
343  mooseAssert((p_lm_u_vel_jac.m() == p_lm_v_vel_jac.m()) &&
344  (p_lm_u_vel_jac.n() == p_lm_v_vel_jac.n()),
345  "We already checked that lm finite element types are the same, so these matrices "
346  "should be the same size");
347  for (const auto i : make_range(p_lm_u_vel_jac.m()))
348  for (const auto j : make_range(p_lm_u_vel_jac.n()))
349  for (const auto qp : make_range(qrule_face.n_points()))
350  {
351  {
352  const Gradient phi(_lm_phi_face[j][qp], 0);
353  p_lm_u_vel_jac(i, j) += JxW_face[qp] * phi * normals[qp] * _scalar_phi_face[i][qp];
354  }
355  {
356  const Gradient phi(0, _lm_phi_face[j][qp]);
357  p_lm_v_vel_jac(i, j) += JxW_face[qp] * phi * normals[qp] * _scalar_phi_face[i][qp];
358  }
359  }
360 }
DiffusionLHDGAssemblyHelper::_scalar_phi_face
const MooseArray< std::vector< Real > > & _scalar_phi_face
DiffusionLHDGAssemblyHelper::_lm_phi_face
const MooseArray< std::vector< Real > > & _lm_phi_face
libMesh::DenseMatrix::m
unsigned int m() const
libMesh::QBase::n_points
unsigned int n_points() const
make_range
IntRange< T > make_range(T beg, T end)
EM::j
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
libMesh::DenseMatrix::n
unsigned int n() const

◆ pressureFaceResidual()

void NavierStokesLHDGAssemblyHelper::pressureFaceResidual ( const MooseArray< Real > &  JxW_face,
const libMesh::QBase &  qrule_face,
const MooseArray< Point > &  normals,
DenseVector< Number > &  pressure_re 
)
protectedinherited

Definition at line 322 of file NavierStokesLHDGAssemblyHelper.C.

Referenced by NavierStokesLHDGOutflowBC::computeResidual(), and computeResidualOnSide().

326 {
327  for (const auto qp : make_range(qrule_face.n_points()))
328  {
329  const Gradient vel(_lm_u_sol[qp], _lm_v_sol[qp]);
330  const auto vdotn = vel * normals[qp];
331  for (const auto i : make_range(pressure_re.size()))
332  pressure_re(i) += JxW_face[qp] * vdotn * _scalar_phi_face[i][qp];
333  }
334 }
DiffusionLHDGAssemblyHelper::_scalar_phi_face
const MooseArray< std::vector< Real > > & _scalar_phi_face
libMesh::QBase::n_points
unsigned int n_points() const
DiffusionLHDGAssemblyHelper::_lm_u_sol
const MooseArray< Number > & _lm_u_sol
make_range
IntRange< T > make_range(T beg, T end)
libMesh::DenseVector::size
virtual unsigned int size() const override final

◆ pressureVolumeJacobian()

void NavierStokesLHDGAssemblyHelper::pressureVolumeJacobian ( const MooseArray< Real > &  JxW,
const libMesh::QBase &  qrule,
DenseMatrix< Number > &  p_u_vel_jac,
DenseMatrix< Number > &  p_v_vel_jac,
DenseMatrix< Number > &  p_global_lm_jac,
DenseMatrix< Number > &  global_lm_p_jac 
)
protectedinherited

Compute the volumetric contributions to the pressure Jacobian, e.g.

the conservation of mass equation

Parameters
i_offsetThe local degree of freedom offset for the pressure
u_j_offsetThe local degree of freedom offset for the x-component of velocity
v_j_offsetThe local degree of freedom offset for the y-component of velocity
p_j_offsetThe local degree of freedom offset for the pressure
global_lm_i_offsetThe local degree of freedom offset for the global Lagrange multiplier that removes the pressure nullspace

Definition at line 258 of file NavierStokesLHDGAssemblyHelper.C.

Referenced by computeJacobian().

264 {
265  for (const auto qp : make_range(qrule.n_points()))
266  {
267  for (const auto i : make_range(p_u_vel_jac.m()))
268  {
269  for (const auto j : make_range(p_u_vel_jac.n()))
270  {
271  {
272  const Gradient phi(_scalar_phi[j][qp], 0);
273  p_u_vel_jac(i, j) -= JxW[qp] * (_grad_scalar_phi[i][qp] * phi);
274  }
275  {
276  const Gradient phi(0, _scalar_phi[j][qp]);
277  p_v_vel_jac(i, j) -= JxW[qp] * (_grad_scalar_phi[i][qp] * phi);
278  }
279  }
280  if (_enclosure_lm_var)
281  p_global_lm_jac(i, 0) -= JxW[qp] * _scalar_phi[i][qp];
282  }
283 
284  if (_enclosure_lm_var)
285  {
286  for (const auto j : make_range(global_lm_p_jac.n()))
287  global_lm_p_jac(0, j) -= JxW[qp] * _scalar_phi[j][qp];
288  }
289  }
290 }
DiffusionLHDGAssemblyHelper::_scalar_phi
const MooseArray< std::vector< Real > > & _scalar_phi
DiffusionLHDGAssemblyHelper::_grad_scalar_phi
const MooseArray< std::vector< RealVectorValue > > & _grad_scalar_phi
libMesh::DenseMatrix::m
unsigned int m() const
libMesh::QBase::n_points
unsigned int n_points() const
make_range
IntRange< T > make_range(T beg, T end)
EM::j
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
libMesh::DenseMatrix::n
unsigned int n() const
NavierStokesLHDGAssemblyHelper::_enclosure_lm_var
const MooseVariableScalar *const _enclosure_lm_var
Definition: NavierStokesLHDGAssemblyHelper.h:216

◆ pressureVolumeResidual()

void NavierStokesLHDGAssemblyHelper::pressureVolumeResidual ( const Moose::Functor< Real > &  pressure_mms_forcing_function,
const MooseArray< Real > &  JxW,
const libMesh::QBase &  qrule,
const Elem *const  current_elem,
const MooseArray< Point > &  q_point,
DenseVector< Number > &  pressure_re,
DenseVector< Number > &  global_lm_re 
)
protectedinherited

Compute the volumetric contributions to the pressure residual, e.g.

the conservation of mass equation

Parameters
i_offsetThe local degree of freedom offset for the pressure
global_lm_i_offsetThe local degree of freedom offset for the global Lagrange multiplier that removes the pressure nullspace

Definition at line 220 of file NavierStokesLHDGAssemblyHelper.C.

Referenced by computeResidual().

228 {
229  for (const auto qp : make_range(qrule.n_points()))
230  {
231  // Prepare forcing function
232  const auto f = pressure_mms_forcing_function(
233  Moose::ElemQpArg{current_elem, qp, &qrule, q_point[qp]}, _ti.determineState());
234 
235  const Gradient vel(_u_sol[qp], _v_sol[qp]);
236  for (const auto i : make_range(pressure_re.size()))
237  {
238  pressure_re(i) -= JxW[qp] * (_grad_scalar_phi[i][qp] * vel);
239 
240  // Pressure equation forcing function RHS
241  pressure_re(i) -= JxW[qp] * _scalar_phi[i][qp] * f;
242 
243  if (_enclosure_lm_var)
244  {
245  mooseAssert(
246  _global_lm_dof_value->size() == 1,
247  "There should only be one degree of freedom for removing the pressure nullspace");
248  pressure_re(i) -= JxW[qp] * _scalar_phi[i][qp] * (*_global_lm_dof_value)[0];
249  }
250  }
251 
252  if (_enclosure_lm_var)
253  global_lm_re(0) -= JxW[qp] * _p_sol[qp];
254  }
255 }
DiffusionLHDGAssemblyHelper::_u_sol
const MooseArray< Number > & _u_sol
DiffusionLHDGAssemblyHelper::_scalar_phi
const MooseArray< std::vector< Real > > & _scalar_phi
TransientInterface::determineState
Moose::StateArg determineState() const
DiffusionLHDGAssemblyHelper::_grad_scalar_phi
const MooseArray< std::vector< RealVectorValue > > & _grad_scalar_phi
MooseArray< Number >::size
unsigned int size() const
f
Real f(Real x)
Test function for Brents method.
Definition: BrentsMethodTest.C:21
libMesh::QBase::n_points
unsigned int n_points() const
DiffusionLHDGAssemblyHelper::_ti
const TransientInterface & _ti
make_range
IntRange< T > make_range(T beg, T end)
libMesh::DenseVector::size
virtual unsigned int size() const override final
NavierStokesLHDGAssemblyHelper::_enclosure_lm_var
const MooseVariableScalar *const _enclosure_lm_var
Definition: NavierStokesLHDGAssemblyHelper.h:216
NavierStokesLHDGAssemblyHelper::_global_lm_dof_value
const MooseArray< Number > *const _global_lm_dof_value
Definition: NavierStokesLHDGAssemblyHelper.h:236

◆ rhoVelCrossVelJacobian()

RealVectorValue NavierStokesLHDGAssemblyHelper::rhoVelCrossVelJacobian ( const MooseArray< Number > &  u_sol,
const MooseArray< Number > &  v_sol,
const unsigned int  qp,
const unsigned int  vel_component,
const unsigned int  vel_j_component,
const MooseArray< std::vector< Real >> &  phi,
const unsigned int  j 
)
protectedinherited
Parameters
u_solThe x-velocity solution, can correspond to either the volumetric or face velocity
v_solThe y-velocity solution, can correspond to either the volumetric or face velocity

Definition at line 303 of file NavierStokesLHDGAssemblyHelper.C.

Referenced by NavierStokesLHDGAssemblyHelper::lmFaceJacobian(), NavierStokesLHDGAssemblyHelper::scalarFaceJacobian(), and NavierStokesLHDGAssemblyHelper::scalarVolumeJacobian().

310 {
311  const RealVectorValue U(u_sol[qp], v_sol[qp]);
312  RealVectorValue vector_phi;
313  vector_phi(vel_j_component) = phi[j][qp];
314  auto ret = vector_phi * U(vel_component);
315  if (vel_component == vel_j_component)
316  ret += U * phi[j][qp];
317  ret *= _rho;
318  return ret;
319 }
EM::j
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")

◆ rhoVelCrossVelResidual()

RealVectorValue NavierStokesLHDGAssemblyHelper::rhoVelCrossVelResidual ( const MooseArray< Number > &  u_sol,
const MooseArray< Number > &  v_sol,
const unsigned int  qp,
const unsigned int  vel_component 
)
protectedinherited
Parameters
u_solThe x-velocity solution, can correspond to either the volumetric or face velocity
v_solThe y-velocity solution, can correspond to either the volumetric or face velocity

Definition at line 293 of file NavierStokesLHDGAssemblyHelper.C.

Referenced by NavierStokesLHDGAssemblyHelper::lmFaceResidual(), NavierStokesLHDGAssemblyHelper::scalarFaceResidual(), and NavierStokesLHDGAssemblyHelper::scalarVolumeResidual().

297 {
298  const RealVectorValue U(u_sol[qp], v_sol[qp]);
299  return _rho * U * U(vel_component);
300 }

◆ scalarDirichletJacobian()

void NavierStokesLHDGAssemblyHelper::scalarDirichletJacobian ( const unsigned int  vel_component,
const MooseArray< Real > &  JxW_face,
const libMesh::QBase &  qrule_face,
const MooseArray< Point > &  normals,
DenseMatrix< Number > &  scalar_vector_jac,
DenseMatrix< Number > &  scalar_scalar_jac,
DenseMatrix< Number > &  scalar_pressure_jac 
)
protectedinherited

Definition at line 600 of file NavierStokesLHDGAssemblyHelper.C.

Referenced by NavierStokesLHDGVelocityDirichletBC::computeJacobian().

607 {
608  DiffusionLHDGAssemblyHelper::scalarDirichletJacobian(
609  JxW_face, qrule_face, normals, scalar_vector_jac, scalar_scalar_jac);
610 
611  for (const auto i : make_range(scalar_pressure_jac.m()))
612  for (const auto j : make_range(scalar_pressure_jac.n()))
613  for (const auto qp : make_range(qrule_face.n_points()))
614  {
615  Gradient p_phi;
616  p_phi(vel_component) = _scalar_phi_face[j][qp];
617  // pressure
618  scalar_pressure_jac(i, j) += JxW_face[qp] * _scalar_phi_face[i][qp] * (p_phi * normals[qp]);
619  }
620 }
DiffusionLHDGAssemblyHelper::_scalar_phi_face
const MooseArray< std::vector< Real > > & _scalar_phi_face
libMesh::DenseMatrix::m
unsigned int m() const
DiffusionLHDGAssemblyHelper::scalarDirichletJacobian
void scalarDirichletJacobian(const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseMatrix< Number > &scalar_vector_jac, DenseMatrix< Number > &scalar_scalar_jac)
libMesh::QBase::n_points
unsigned int n_points() const
make_range
IntRange< T > make_range(T beg, T end)
EM::j
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
libMesh::DenseMatrix::n
unsigned int n() const

◆ scalarDirichletResidual()

void NavierStokesLHDGAssemblyHelper::scalarDirichletResidual ( const MooseArray< Gradient > &  vector_sol,
const MooseArray< Number > &  scalar_sol,
const unsigned int  vel_component,
const std::array< const Moose::Functor< Real > *, 3 > &  dirichlet_vel,
const MooseArray< Real > &  JxW_face,
const libMesh::QBase &  qrule_face,
const MooseArray< Point > &  normals,
const Elem *const  current_elem,
const unsigned int  current_side,
const MooseArray< Point > &  q_point_face,
DenseVector< Number > &  scalar_re 
)
protectedinherited

Definition at line 550 of file NavierStokesLHDGAssemblyHelper.C.

Referenced by NavierStokesLHDGVelocityDirichletBC::computeResidual().

562 {
563  DiffusionLHDGAssemblyHelper::scalarDirichletResidual(vector_sol,
564  scalar_sol,
565  *dirichlet_vel[vel_component],
566  JxW_face,
567  qrule_face,
568  normals,
569  current_elem,
570  current_side,
571  q_point_face,
572  scalar_re);
573 
574  for (const auto qp : make_range(qrule_face.n_points()))
575  {
576  Gradient qp_p;
577  qp_p(vel_component) = _p_sol[qp];
578 
579  const Moose::ElemSideQpArg elem_side_qp_arg{
580  current_elem, current_side, qp, &qrule_face, q_point_face[qp]};
581  const auto time_arg = _ti.determineState();
582  const RealVectorValue dirichlet_velocity((*dirichlet_vel[0])(elem_side_qp_arg, time_arg),
583  (*dirichlet_vel[1])(elem_side_qp_arg, time_arg),
584  (*dirichlet_vel[2])(elem_side_qp_arg, time_arg));
585  const auto scalar_value = dirichlet_velocity(vel_component);
586 
587  for (const auto i : make_range(scalar_re.size()))
588  {
589  // pressure
590  scalar_re(i) += JxW_face[qp] * _scalar_phi_face[i][qp] * (qp_p * normals[qp]);
591 
592  // dirichlet lm from advection term
593  scalar_re(i) += JxW_face[qp] * _scalar_phi_face[i][qp] *
594  (_rho * dirichlet_velocity * normals[qp]) * scalar_value;
595  }
596  }
597 }
DiffusionLHDGAssemblyHelper::_scalar_phi_face
const MooseArray< std::vector< Real > > & _scalar_phi_face
DiffusionLHDGAssemblyHelper::scalarDirichletResidual
void scalarDirichletResidual(const MooseArray< Gradient > &vector_sol, const MooseArray< Number > &scalar_sol, const Moose::Functor< Real > &dirichlet_value, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, const Elem *const current_elem, const unsigned int current_side, const MooseArray< Point > &q_point_face, DenseVector< Number > &scalar_re)
TransientInterface::determineState
Moose::StateArg determineState() const
libMesh::QBase::n_points
unsigned int n_points() const
DiffusionLHDGAssemblyHelper::_ti
const TransientInterface & _ti
make_range
IntRange< T > make_range(T beg, T end)
libMesh::DenseVector::size
virtual unsigned int size() const override final

◆ scalarFaceJacobian()

void NavierStokesLHDGAssemblyHelper::scalarFaceJacobian ( const unsigned int  vel_component,
const MooseArray< Real > &  JxW_face,
const libMesh::QBase &  qrule_face,
const MooseArray< Point > &  normals,
DenseMatrix< Number > &  scalar_vector_jac,
DenseMatrix< Number > &  scalar_scalar_jac,
DenseMatrix< Number > &  scalar_lm_jac,
DenseMatrix< Number > &  scalar_p_jac,
DenseMatrix< Number > &  scalar_lm_u_vel_jac,
DenseMatrix< Number > &  scalar_lm_v_vel_jac 
)
protectedinherited

Definition at line 393 of file NavierStokesLHDGAssemblyHelper.C.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

404 {
405  DiffusionLHDGAssemblyHelper::scalarFaceJacobian(
406  JxW_face, qrule_face, normals, scalar_vector_jac, scalar_scalar_jac, scalar_lm_jac);
407 
408  for (const auto i : make_range(scalar_lm_u_vel_jac.m()))
409  for (const auto qp : make_range(qrule_face.n_points()))
410  {
411  for (const auto j : make_range(scalar_p_jac.n()))
412  {
413  Gradient p_phi;
414  p_phi(vel_component) = _scalar_phi_face[j][qp];
415  // pressure
416  scalar_p_jac(i, j) += JxW_face[qp] * _scalar_phi_face[i][qp] * (p_phi * normals[qp]);
417  }
418 
419  for (const auto j : make_range(scalar_lm_u_vel_jac.n()))
420  {
421  //
422  // from convection term
423  //
424 
425  // derivatives wrt 0th component of velocity
426  {
427  const auto rho_vel_cross_vel =
428  rhoVelCrossVelJacobian(_lm_u_sol, _lm_v_sol, qp, vel_component, 0, _lm_phi_face, j);
429  scalar_lm_u_vel_jac(i, j) +=
430  JxW_face[qp] * _scalar_phi_face[i][qp] * rho_vel_cross_vel * normals[qp];
431  }
432  // derivatives wrt 1th component of velocity
433  {
434  const auto rho_vel_cross_vel =
435  rhoVelCrossVelJacobian(_lm_u_sol, _lm_v_sol, qp, vel_component, 1, _lm_phi_face, j);
436  scalar_lm_v_vel_jac(i, j) +=
437  JxW_face[qp] * _scalar_phi_face[i][qp] * rho_vel_cross_vel * normals[qp];
438  }
439  }
440  }
441 }
DiffusionLHDGAssemblyHelper::_scalar_phi_face
const MooseArray< std::vector< Real > > & _scalar_phi_face
NavierStokesLHDGAssemblyHelper::rhoVelCrossVelJacobian
RealVectorValue rhoVelCrossVelJacobian(const MooseArray< Number > &u_sol, const MooseArray< Number > &v_sol, const unsigned int qp, const unsigned int vel_component, const unsigned int vel_j_component, const MooseArray< std::vector< Real >> &phi, const unsigned int j)
Definition: NavierStokesLHDGAssemblyHelper.C:303
DiffusionLHDGAssemblyHelper::_lm_phi_face
const MooseArray< std::vector< Real > > & _lm_phi_face
libMesh::DenseMatrix::m
unsigned int m() const
DiffusionLHDGAssemblyHelper::scalarFaceJacobian
void scalarFaceJacobian(const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseMatrix< Number > &scalar_vector_jac, DenseMatrix< Number > &scalar_scalar_jac, DenseMatrix< Number > &scalar_lm_jac)
libMesh::QBase::n_points
unsigned int n_points() const
DiffusionLHDGAssemblyHelper::_lm_u_sol
const MooseArray< Number > & _lm_u_sol
make_range
IntRange< T > make_range(T beg, T end)
EM::j
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
libMesh::DenseMatrix::n
unsigned int n() const

◆ scalarFaceResidual()

void NavierStokesLHDGAssemblyHelper::scalarFaceResidual ( const MooseArray< Gradient > &  vector_sol,
const MooseArray< Number > &  scalar_sol,
const MooseArray< Number > &  lm_sol,
const unsigned int  vel_component,
const MooseArray< Real > &  JxW_face,
const libMesh::QBase &  qrule_face,
const MooseArray< Point > &  normals,
DenseVector< Number > &  scalar_re 
)
protectedinherited

Definition at line 363 of file NavierStokesLHDGAssemblyHelper.C.

Referenced by NavierStokesLHDGOutflowBC::computeResidual(), and computeResidualOnSide().

371 {
372  DiffusionLHDGAssemblyHelper::scalarFaceResidual(
373  vector_sol, scalar_sol, lm_sol, JxW_face, qrule_face, normals, scalar_re);
374 
375  for (const auto qp : make_range(qrule_face.n_points()))
376  {
377  Gradient qp_p;
378  qp_p(vel_component) = _p_sol[qp];
379  const auto rho_vel_cross_vel = rhoVelCrossVelResidual(_lm_u_sol, _lm_v_sol, qp, vel_component);
380 
381  for (const auto i : make_range(scalar_re.size()))
382  {
383  // pressure
384  scalar_re(i) += JxW_face[qp] * _scalar_phi_face[i][qp] * (qp_p * normals[qp]);
385 
386  // lm from convection term
387  scalar_re(i) += JxW_face[qp] * _scalar_phi_face[i][qp] * rho_vel_cross_vel * normals[qp];
388  }
389  }
390 }
DiffusionLHDGAssemblyHelper::_scalar_phi_face
const MooseArray< std::vector< Real > > & _scalar_phi_face
NavierStokesLHDGAssemblyHelper::rhoVelCrossVelResidual
RealVectorValue rhoVelCrossVelResidual(const MooseArray< Number > &u_sol, const MooseArray< Number > &v_sol, const unsigned int qp, const unsigned int vel_component)
Definition: NavierStokesLHDGAssemblyHelper.C:293
libMesh::QBase::n_points
unsigned int n_points() const
DiffusionLHDGAssemblyHelper::_lm_u_sol
const MooseArray< Number > & _lm_u_sol
make_range
IntRange< T > make_range(T beg, T end)
libMesh::DenseVector::size
virtual unsigned int size() const override final
DiffusionLHDGAssemblyHelper::scalarFaceResidual
void scalarFaceResidual(const MooseArray< Gradient > &vector_sol, const MooseArray< Number > &scalar_sol, const MooseArray< Number > &lm_sol, const MooseArray< Real > &JxW_face, const libMesh::QBase &qrule_face, const MooseArray< Point > &normals, DenseVector< Number > &scalar_re)

◆ scalarVolumeJacobian()

void NavierStokesLHDGAssemblyHelper::scalarVolumeJacobian ( const unsigned int  vel_component,
const MooseArray< Real > &  JxW,
const libMesh::QBase &  qrule,
DenseMatrix< Number > &  scalar_vector_jac,
DenseMatrix< Number > &  scalar_u_vel_jac,
DenseMatrix< Number > &  scalar_v_vel_jac,
DenseMatrix< Number > &  scalar_p_jac 
)
protectedinherited

Compute the volumetric contributions to a velocity Jacobian.

Parameters
i_offsetThe local degree of freedom offset for the velocity component
vel_gradient_j_offsetThe local degree of freedom offset for the associated velocity gradient
p_j_offsetThe local degree of freedom offset for the pressure
vel_componentThe velocity component
u_j_offsetThe local degree of freedom offset for the x-component velocity
v_j_offsetThe local degree of freedom offset for the y-component velocity

Definition at line 178 of file NavierStokesLHDGAssemblyHelper.C.

Referenced by computeJacobian().

185 {
186  DiffusionLHDGAssemblyHelper::scalarVolumeJacobian(JxW, qrule, scalar_vector_jac);
187 
188  for (const auto i : make_range(scalar_vector_jac.m()))
189  for (const auto qp : make_range(qrule.n_points()))
190  {
191  // Scalar equation dependence on pressure dofs
192  for (const auto j : make_range(scalar_p_jac.n()))
193  {
194  Gradient p_phi;
195  p_phi(vel_component) = _scalar_phi[j][qp];
196  scalar_p_jac(i, j) -= JxW[qp] * (_grad_scalar_phi[i][qp] * p_phi);
197  }
198 
199  // Scalar equation dependence on scalar dofs
200  mooseAssert(scalar_u_vel_jac.n() == scalar_v_vel_jac.n(), "These must be the same size");
201  for (const auto j : make_range(scalar_u_vel_jac.n()))
202  {
203  // derivatives wrt 0th component of velocity
204  {
205  const auto rho_vel_cross_vel =
206  rhoVelCrossVelJacobian(_u_sol, _v_sol, qp, vel_component, 0, _scalar_phi, j);
207  scalar_u_vel_jac(i, j) -= JxW[qp] * (_grad_scalar_phi[i][qp] * rho_vel_cross_vel);
208  }
209  // derivatives wrt 1th component of velocity
210  {
211  const auto rho_vel_cross_vel =
212  rhoVelCrossVelJacobian(_u_sol, _v_sol, qp, vel_component, 1, _scalar_phi, j);
213  scalar_v_vel_jac(i, j) -= JxW[qp] * (_grad_scalar_phi[i][qp] * rho_vel_cross_vel);
214  }
215  }
216  }
217 }
DiffusionLHDGAssemblyHelper::_u_sol
const MooseArray< Number > & _u_sol
DiffusionLHDGAssemblyHelper::_scalar_phi
const MooseArray< std::vector< Real > > & _scalar_phi
NavierStokesLHDGAssemblyHelper::rhoVelCrossVelJacobian
RealVectorValue rhoVelCrossVelJacobian(const MooseArray< Number > &u_sol, const MooseArray< Number > &v_sol, const unsigned int qp, const unsigned int vel_component, const unsigned int vel_j_component, const MooseArray< std::vector< Real >> &phi, const unsigned int j)
Definition: NavierStokesLHDGAssemblyHelper.C:303
DiffusionLHDGAssemblyHelper::_grad_scalar_phi
const MooseArray< std::vector< RealVectorValue > > & _grad_scalar_phi
libMesh::DenseMatrix::m
unsigned int m() const
libMesh::QBase::n_points
unsigned int n_points() const
make_range
IntRange< T > make_range(T beg, T end)
DiffusionLHDGAssemblyHelper::scalarVolumeJacobian
void scalarVolumeJacobian(const MooseArray< Real > &JxW, const libMesh::QBase &qrule, DenseMatrix< Number > &scalar_vector_jac)
EM::j
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
libMesh::DenseMatrix::n
unsigned int n() const

◆ scalarVolumeResidual()

void NavierStokesLHDGAssemblyHelper::scalarVolumeResidual ( const MooseArray< Gradient > &  vel_gradient,
const unsigned int  vel_component,
const Moose::Functor< Real > &  body_force,
const MooseArray< Real > &  JxW,
const libMesh::QBase &  qrule,
const Elem *const  current_elem,
const MooseArray< Point > &  q_point,
DenseVector< Number > &  scalar_re 
)
protectedinherited

Compute the volumetric contributions to a velocity residual for a provided velocity gradient and stress.

Parameters
i_offsetThe local degree of freedom offset for the velocity component
vel_gradientThe velocity gradient component
vel_componentThe velocity component

Definition at line 148 of file NavierStokesLHDGAssemblyHelper.C.

Referenced by computeResidual().

156 {
157  DiffusionLHDGAssemblyHelper::scalarVolumeResidual(
158  vel_gradient, body_force, JxW, qrule, current_elem, q_point, scalar_re);
159 
160  for (const auto qp : make_range(qrule.n_points()))
161  {
162  const auto rho_vel_cross_vel = rhoVelCrossVelResidual(_u_sol, _v_sol, qp, vel_component);
163  Gradient qp_p;
164  qp_p(vel_component) = _p_sol[qp];
165 
166  for (const auto i : index_range(scalar_re))
167  {
168  // Scalar equation dependence on pressure dofs
169  scalar_re(i) -= JxW[qp] * (_grad_scalar_phi[i][qp] * qp_p);
170 
171  // Scalar equation dependence on scalar dofs
172  scalar_re(i) -= JxW[qp] * (_grad_scalar_phi[i][qp] * rho_vel_cross_vel);
173  }
174  }
175 }
DiffusionLHDGAssemblyHelper::_u_sol
const MooseArray< Number > & _u_sol
DiffusionLHDGAssemblyHelper::_grad_scalar_phi
const MooseArray< std::vector< RealVectorValue > > & _grad_scalar_phi
NavierStokesLHDGAssemblyHelper::rhoVelCrossVelResidual
RealVectorValue rhoVelCrossVelResidual(const MooseArray< Number > &u_sol, const MooseArray< Number > &v_sol, const unsigned int qp, const unsigned int vel_component)
Definition: NavierStokesLHDGAssemblyHelper.C:293
libMesh::QBase::n_points
unsigned int n_points() const
DiffusionLHDGAssemblyHelper::scalarVolumeResidual
void scalarVolumeResidual(const MooseArray< Gradient > &vector_field, const Moose::Functor< Real > &source, const MooseArray< Real > &JxW, const libMesh::QBase &qrule, const Elem *const current_elem, const MooseArray< Point > &q_point, DenseVector< Number > &scalar_re)
make_range
IntRange< T > make_range(T beg, T end)
index_range
auto index_range(const T &sizable)

◆ validParams()

InputParameters NavierStokesLHDGKernel::validParams ( )
static

Definition at line 18 of file NavierStokesLHDGKernel.C.

19 {
20  auto params = HDGKernel::validParams();
21  params += NavierStokesLHDGAssemblyHelper::validParams();
22  params.addParam<MooseFunctorName>(
23  "body_force_x", 0, "Body force for the momentum equation in the x-direction");
24  params.addParam<MooseFunctorName>(
25  "body_force_y", 0, "Body force for the momentum equation in the y-direction");
26  params.addParam<MooseFunctorName>(
27  "body_force_z", 0, "Body force for the momentum equation in the z-direction");
28  params.addParam<MooseFunctorName>(
29  "pressure_mms_forcing_function",
30  0,
31  "A forcing function for the pressure (mass) equation for conducting MMS studies");
32  params.addClassDescription("Implements the steady incompressible Navier-Stokes equations for a "
33  "hybridized discretization");
34  params.renameParam("variable", "u", "The x-component of velocity");
35  return params;
36 }
HDGKernel::validParams
static InputParameters validParams()

Member Data Documentation

◆ _body_force_x

const Moose::Functor<Real>& NavierStokesLHDGKernel::_body_force_x
protected

Definition at line 44 of file NavierStokesLHDGKernel.h.

Referenced by computeResidual(), and NavierStokesLHDGKernel().

◆ _body_force_y

const Moose::Functor<Real>& NavierStokesLHDGKernel::_body_force_y
protected

Definition at line 45 of file NavierStokesLHDGKernel.h.

Referenced by computeResidual(), and NavierStokesLHDGKernel().

◆ _body_force_z

const Moose::Functor<Real>& NavierStokesLHDGKernel::_body_force_z
protected

Definition at line 46 of file NavierStokesLHDGKernel.h.

Referenced by NavierStokesLHDGKernel().

◆ _body_forces

std::vector<const Moose::Functor<Real> *> NavierStokesLHDGKernel::_body_forces
protected

Definition at line 47 of file NavierStokesLHDGKernel.h.

Referenced by NavierStokesLHDGKernel().

◆ _current_side

const unsigned int& NavierStokesLHDGKernel::_current_side
protected

The current side when doing face evaluations.

Definition at line 63 of file NavierStokesLHDGKernel.h.

Referenced by computeJacobianOnSide(), and computeResidualOnSide().

◆ _enclosure_lm_var

const MooseVariableScalar* const NavierStokesLHDGAssemblyHelper::_enclosure_lm_var
protectedinherited

Definition at line 216 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by additionalROVariables(), computeJacobian(), computeResidual(), NavierStokesLHDGAssemblyHelper::pressureVolumeJacobian(), and NavierStokesLHDGAssemblyHelper::pressureVolumeResidual().

◆ _global_lm_dof_indices

const std::vector<dof_id_type>* const NavierStokesLHDGAssemblyHelper::_global_lm_dof_indices
protectedinherited

Definition at line 226 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeJacobian(), and computeResidual().

◆ _global_lm_dof_value

const MooseArray<Number>* const NavierStokesLHDGAssemblyHelper::_global_lm_dof_value
protectedinherited

Definition at line 236 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGAssemblyHelper::pressureVolumeResidual().

◆ _global_lm_p_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_global_lm_p_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeJacobian().

◆ _global_lm_re

DenseVector<Number> NavierStokesLHDGAssemblyHelper::_global_lm_re
protectedinherited

Definition at line 241 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeResidual().

◆ _grad_u_grad_u_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_grad_u_grad_u_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeJacobian().

◆ _grad_u_lm_u_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_grad_u_lm_u_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _grad_u_u_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_grad_u_u_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeJacobian().

◆ _grad_u_vel_re

DenseVector<Number> NavierStokesLHDGAssemblyHelper::_grad_u_vel_re
protectedinherited

Definition at line 241 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeResidual(), NavierStokesLHDGVelocityDirichletBC::computeResidual(), NavierStokesLHDGOutflowBC::computeResidual(), and computeResidualOnSide().

◆ _grad_v_grad_v_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_grad_v_grad_v_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeJacobian().

◆ _grad_v_lm_v_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_grad_v_lm_v_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _grad_v_v_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_grad_v_v_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeJacobian().

◆ _grad_v_var

const MooseVariableFE<RealVectorValue>& NavierStokesLHDGAssemblyHelper::_grad_v_var
protectedinherited

Definition at line 211 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by additionalROVariables(), computeJacobian(), NavierStokesLHDGOutflowBC::computeJacobian(), computeJacobianOnSide(), computeResidual(), NavierStokesLHDGVelocityDirichletBC::computeResidual(), NavierStokesLHDGOutflowBC::computeResidual(), computeResidualOnSide(), and NavierStokesLHDGAssemblyHelper::NavierStokesLHDGAssemblyHelper().

◆ _grad_v_vel_re

DenseVector<Number> NavierStokesLHDGAssemblyHelper::_grad_v_vel_re
protectedinherited

Definition at line 241 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeResidual(), NavierStokesLHDGVelocityDirichletBC::computeResidual(), NavierStokesLHDGOutflowBC::computeResidual(), and computeResidualOnSide().

◆ _grad_w_var

const MooseVariableFE<RealVectorValue>* const NavierStokesLHDGAssemblyHelper::_grad_w_var
protectedinherited

Definition at line 212 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGAssemblyHelper::NavierStokesLHDGAssemblyHelper().

◆ _JxW_face

const MooseArray<Real>& NavierStokesLHDGKernel::_JxW_face
protected

transformed Jacobian weights on the current element face

Definition at line 57 of file NavierStokesLHDGKernel.h.

Referenced by computeJacobianOnSide(), and computeResidualOnSide().

◆ _lm_u_grad_u_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_lm_u_grad_u_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _lm_u_lm_u_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_lm_u_lm_u_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGVelocityDirichletBC::computeJacobian(), NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _lm_u_lm_v_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_lm_u_lm_v_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _lm_u_p_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_lm_u_p_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _lm_u_u_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_lm_u_u_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _lm_u_vel_re

DenseVector<Number> NavierStokesLHDGAssemblyHelper::_lm_u_vel_re
protectedinherited

Definition at line 241 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGVelocityDirichletBC::computeResidual(), NavierStokesLHDGOutflowBC::computeResidual(), and computeResidualOnSide().

◆ _lm_v_dof_indices

const std::vector<dof_id_type>& NavierStokesLHDGAssemblyHelper::_lm_v_dof_indices
protectedinherited

Definition at line 221 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGVelocityDirichletBC::computeJacobian(), NavierStokesLHDGOutflowBC::computeJacobian(), computeJacobianOnSide(), NavierStokesLHDGVelocityDirichletBC::computeResidual(), NavierStokesLHDGOutflowBC::computeResidual(), and computeResidualOnSide().

◆ _lm_v_grad_v_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_lm_v_grad_v_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _lm_v_lm_u_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_lm_v_lm_u_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _lm_v_lm_v_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_lm_v_lm_v_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGVelocityDirichletBC::computeJacobian(), NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _lm_v_p_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_lm_v_p_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _lm_v_sol

const MooseArray<Number>& NavierStokesLHDGAssemblyHelper::_lm_v_sol
protectedinherited

Definition at line 231 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGVelocityDirichletBC::computeResidual(), NavierStokesLHDGOutflowBC::computeResidual(), computeResidualOnSide(), NavierStokesLHDGAssemblyHelper::lmFaceJacobian(), NavierStokesLHDGAssemblyHelper::lmFaceResidual(), NavierStokesLHDGAssemblyHelper::pressureFaceResidual(), NavierStokesLHDGAssemblyHelper::scalarFaceJacobian(), and NavierStokesLHDGAssemblyHelper::scalarFaceResidual().

◆ _lm_v_v_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_lm_v_v_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _lm_v_vel_re

DenseVector<Number> NavierStokesLHDGAssemblyHelper::_lm_v_vel_re
protectedinherited

Definition at line 241 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGVelocityDirichletBC::computeResidual(), NavierStokesLHDGOutflowBC::computeResidual(), and computeResidualOnSide().

◆ _lm_w_dof_indices

const std::vector<dof_id_type>* const NavierStokesLHDGAssemblyHelper::_lm_w_dof_indices
protectedinherited

Definition at line 224 of file NavierStokesLHDGAssemblyHelper.h.

◆ _lm_w_sol

const MooseArray<Number>* const NavierStokesLHDGAssemblyHelper::_lm_w_sol
protectedinherited

Definition at line 234 of file NavierStokesLHDGAssemblyHelper.h.

◆ _normals

const MooseArray<Point>& NavierStokesLHDGKernel::_normals
protected

face normals

Definition at line 60 of file NavierStokesLHDGKernel.h.

Referenced by computeJacobianOnSide(), and computeResidualOnSide().

◆ _p_dof_indices

const std::vector<dof_id_type>& NavierStokesLHDGAssemblyHelper::_p_dof_indices
protectedinherited

Definition at line 225 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGVelocityDirichletBC::computeJacobian(), computeJacobian(), NavierStokesLHDGOutflowBC::computeJacobian(), computeJacobianOnSide(), computeResidual(), NavierStokesLHDGVelocityDirichletBC::computeResidual(), NavierStokesLHDGOutflowBC::computeResidual(), and computeResidualOnSide().

◆ _p_global_lm_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_p_global_lm_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeJacobian().

◆ _p_lm_u_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_p_lm_u_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _p_lm_v_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_p_lm_v_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _p_re

DenseVector<Number> NavierStokesLHDGAssemblyHelper::_p_re
protectedinherited

Definition at line 241 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeResidual(), NavierStokesLHDGVelocityDirichletBC::computeResidual(), NavierStokesLHDGOutflowBC::computeResidual(), and computeResidualOnSide().

◆ _p_sol

const MooseArray<Number>& NavierStokesLHDGAssemblyHelper::_p_sol
protectedinherited

Definition at line 235 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGAssemblyHelper::lmFaceResidual(), NavierStokesLHDGAssemblyHelper::pressureVolumeResidual(), NavierStokesLHDGAssemblyHelper::scalarDirichletResidual(), NavierStokesLHDGAssemblyHelper::scalarFaceResidual(), and NavierStokesLHDGAssemblyHelper::scalarVolumeResidual().

◆ _p_u_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_p_u_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeJacobian().

◆ _p_v_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_p_v_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeJacobian().

◆ _pressure_mms_forcing_function

const Moose::Functor<Real>& NavierStokesLHDGKernel::_pressure_mms_forcing_function
protected

Definition at line 48 of file NavierStokesLHDGKernel.h.

Referenced by computeResidual().

◆ _pressure_var

const MooseVariableFE<Real>& NavierStokesLHDGAssemblyHelper::_pressure_var
protectedinherited

Definition at line 215 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by additionalROVariables(), computeJacobian(), NavierStokesLHDGOutflowBC::computeJacobian(), computeJacobianOnSide(), computeResidual(), NavierStokesLHDGVelocityDirichletBC::computeResidual(), NavierStokesLHDGOutflowBC::computeResidual(), computeResidualOnSide(), and NavierStokesLHDGAssemblyHelper::NavierStokesLHDGAssemblyHelper().

◆ _q_point_face

const MooseArray<Point>& NavierStokesLHDGKernel::_q_point_face
protected

The physical locations of the quadrature points on the face.

Definition at line 54 of file NavierStokesLHDGKernel.h.

◆ _qrule_face

const QBase* const& NavierStokesLHDGKernel::_qrule_face
protected

The face quadrature rule.

Definition at line 51 of file NavierStokesLHDGKernel.h.

Referenced by computeJacobianOnSide(), and computeResidualOnSide().

◆ _qv_dof_indices

const std::vector<dof_id_type>& NavierStokesLHDGAssemblyHelper::_qv_dof_indices
protectedinherited

Containers for dof indices.

Definition at line 219 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGVelocityDirichletBC::computeJacobian(), computeJacobian(), NavierStokesLHDGOutflowBC::computeJacobian(), computeJacobianOnSide(), computeResidual(), NavierStokesLHDGVelocityDirichletBC::computeResidual(), NavierStokesLHDGOutflowBC::computeResidual(), and computeResidualOnSide().

◆ _qv_sol

const MooseArray<Gradient>& NavierStokesLHDGAssemblyHelper::_qv_sol
protectedinherited

local solutions at quadrature points

Definition at line 229 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeResidual(), NavierStokesLHDGVelocityDirichletBC::computeResidual(), NavierStokesLHDGOutflowBC::computeResidual(), and computeResidualOnSide().

◆ _qw_dof_indices

const std::vector<dof_id_type>* const NavierStokesLHDGAssemblyHelper::_qw_dof_indices
protectedinherited

Definition at line 222 of file NavierStokesLHDGAssemblyHelper.h.

◆ _qw_sol

const MooseArray<Gradient>* const NavierStokesLHDGAssemblyHelper::_qw_sol
protectedinherited

Definition at line 232 of file NavierStokesLHDGAssemblyHelper.h.

◆ _rho

const Real NavierStokesLHDGAssemblyHelper::_rho
protectedinherited

The density.

Definition at line 239 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGAssemblyHelper::rhoVelCrossVelJacobian(), NavierStokesLHDGAssemblyHelper::rhoVelCrossVelResidual(), and NavierStokesLHDGAssemblyHelper::scalarDirichletResidual().

◆ _u_grad_u_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_u_grad_u_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGVelocityDirichletBC::computeJacobian(), computeJacobian(), NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _u_lm_u_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_u_lm_u_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _u_lm_v_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_u_lm_v_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _u_p_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_u_p_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGVelocityDirichletBC::computeJacobian(), computeJacobian(), NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _u_u_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_u_u_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGVelocityDirichletBC::computeJacobian(), computeJacobian(), NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _u_v_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_u_v_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeJacobian().

◆ _u_vel_re

DenseVector<Number> NavierStokesLHDGAssemblyHelper::_u_vel_re
protectedinherited

Definition at line 241 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeResidual(), NavierStokesLHDGVelocityDirichletBC::computeResidual(), NavierStokesLHDGOutflowBC::computeResidual(), and computeResidualOnSide().

◆ _v_dof_indices

const std::vector<dof_id_type>& NavierStokesLHDGAssemblyHelper::_v_dof_indices
protectedinherited

Definition at line 220 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGVelocityDirichletBC::computeJacobian(), computeJacobian(), NavierStokesLHDGOutflowBC::computeJacobian(), computeJacobianOnSide(), computeResidual(), NavierStokesLHDGVelocityDirichletBC::computeResidual(), NavierStokesLHDGOutflowBC::computeResidual(), and computeResidualOnSide().

◆ _v_face_var

const MooseVariableFE<Real>& NavierStokesLHDGAssemblyHelper::_v_face_var
protectedinherited

Definition at line 213 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by additionalROVariables(), NavierStokesLHDGVelocityDirichletBC::computeJacobian(), NavierStokesLHDGOutflowBC::computeJacobian(), computeJacobianOnSide(), NavierStokesLHDGVelocityDirichletBC::computeResidual(), NavierStokesLHDGOutflowBC::computeResidual(), computeResidualOnSide(), and NavierStokesLHDGAssemblyHelper::NavierStokesLHDGAssemblyHelper().

◆ _v_grad_v_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_v_grad_v_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGVelocityDirichletBC::computeJacobian(), computeJacobian(), NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _v_lm_u_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_v_lm_u_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _v_lm_v_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_v_lm_v_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _v_p_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_v_p_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGVelocityDirichletBC::computeJacobian(), computeJacobian(), NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _v_sol

const MooseArray<Number>& NavierStokesLHDGAssemblyHelper::_v_sol
protectedinherited

Definition at line 230 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeResidual(), NavierStokesLHDGVelocityDirichletBC::computeResidual(), NavierStokesLHDGOutflowBC::computeResidual(), computeResidualOnSide(), NavierStokesLHDGAssemblyHelper::pressureVolumeResidual(), NavierStokesLHDGAssemblyHelper::scalarVolumeJacobian(), and NavierStokesLHDGAssemblyHelper::scalarVolumeResidual().

◆ _v_u_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_v_u_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeJacobian().

◆ _v_v_jac

DenseMatrix<Number> NavierStokesLHDGAssemblyHelper::_v_v_jac
protectedinherited

Definition at line 243 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGVelocityDirichletBC::computeJacobian(), computeJacobian(), NavierStokesLHDGOutflowBC::computeJacobian(), and computeJacobianOnSide().

◆ _v_var

const MooseVariableFE<Real>& NavierStokesLHDGAssemblyHelper::_v_var
protectedinherited

Definition at line 209 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by additionalROVariables(), NavierStokesLHDGVelocityDirichletBC::computeJacobian(), computeJacobian(), NavierStokesLHDGOutflowBC::computeJacobian(), computeJacobianOnSide(), computeResidual(), NavierStokesLHDGVelocityDirichletBC::computeResidual(), NavierStokesLHDGOutflowBC::computeResidual(), computeResidualOnSide(), and NavierStokesLHDGAssemblyHelper::NavierStokesLHDGAssemblyHelper().

◆ _v_vel_re

DenseVector<Number> NavierStokesLHDGAssemblyHelper::_v_vel_re
protectedinherited

Definition at line 241 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by computeResidual(), NavierStokesLHDGVelocityDirichletBC::computeResidual(), NavierStokesLHDGOutflowBC::computeResidual(), and computeResidualOnSide().

◆ _w_dof_indices

const std::vector<dof_id_type>* const NavierStokesLHDGAssemblyHelper::_w_dof_indices
protectedinherited

Definition at line 223 of file NavierStokesLHDGAssemblyHelper.h.

◆ _w_face_var

const MooseVariableFE<Real>* const NavierStokesLHDGAssemblyHelper::_w_face_var
protectedinherited

Definition at line 214 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGAssemblyHelper::NavierStokesLHDGAssemblyHelper().

◆ _w_sol

const MooseArray<Number>* const NavierStokesLHDGAssemblyHelper::_w_sol
protectedinherited

Definition at line 233 of file NavierStokesLHDGAssemblyHelper.h.

◆ _w_var

const MooseVariableFE<Real>* const NavierStokesLHDGAssemblyHelper::_w_var
protectedinherited

Definition at line 210 of file NavierStokesLHDGAssemblyHelper.h.

Referenced by NavierStokesLHDGAssemblyHelper::NavierStokesLHDGAssemblyHelper().

The documentation for this class was generated from the following files:
Generated on Thu Jul 17 2025 01:49:32 for https://mooseframework.inl.gov by   doxygen 1.8.14