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

Computes the numerical flux for FlowModelGasMix using the HLLC approximate Riemann solver. More...

#include <NumericalFluxGasMixHLLC.h>

Inheritance diagram for NumericalFluxGasMixHLLC:
[legend]

Public Types

typedef DataFileName DataFileParameterType
 

Public Member Functions

 NumericalFluxGasMixHLLC (const InputParameters &parameters)
 
virtual void calcFlux (const std::vector< ADReal > &UL, const std::vector< ADReal > &UR, const RealVectorValue &nLR, const RealVectorValue &t1, const RealVectorValue &t2, std::vector< ADReal > &FL, std::vector< ADReal > &FR) const override
 Calculates the 3D flux vectors given "left" and "right" states. More...
 
virtual unsigned int getNumberOfRegions () const override
 Returns the total possible number of regions. More...
 
virtual void execute () override
 
virtual void initialize () override
 
virtual void finalize () override
 
virtual void threadJoin (const UserObject &) override
 
virtual const std::vector< ADReal > & getFlux (const unsigned int iside, const dof_id_type ielem, bool res_side_is_left, const std::vector< ADReal > &UL_1d, const std::vector< ADReal > &UR_1d, Real nLR_dot_d) const
 Gets the 1D flux vector for an element/side combination. More...
 
virtual const std::vector< ADReal > & getFlux3D (const unsigned int iside, const dof_id_type ielem, const std::vector< ADReal > &UL_3d, const std::vector< ADReal > &UR_3d, const RealVectorValue &nLR, const RealVectorValue &t1, const RealVectorValue &t2) const
 Gets the 3D flux vector for an element/side combination on the left side. More...
 
unsigned int getLastRegionIndex () const
 Returns the index of the region last entered. More...
 
virtual void subdomainSetup () override
 
bool needThreadedCopy () const override final
 
SubProblemgetSubProblem () const
 
bool shouldDuplicateInitialExecution () const
 
virtual Real spatialValue (const Point &) const
 
virtual const std::vector< Point > spatialPoints () const
 
void gatherSum (T &value)
 
void gatherMax (T &value)
 
void gatherMin (T &value)
 
void gatherProxyValueMax (T1 &proxy, T2 &value)
 
void gatherProxyValueMin (T1 &proxy, T2 &value)
 
void setPrimaryThreadCopy (UserObject *primary)
 
UserObjectprimaryThreadCopy ()
 
std::set< UserObjectName > getDependObjects () const
 
const std::set< std::string > & getRequestedItems () override
 
const std::set< std::string > & getSuppliedItems () override
 
unsigned int systemNumber () const
 
virtual bool enabled () const
 
std::shared_ptr< MooseObjectgetSharedPtr ()
 
std::shared_ptr< const MooseObjectgetSharedPtr () const
 
MooseAppgetMooseApp () 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 InputParametersparameters () 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
 
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 initialSetup ()
 
virtual void timestepSetup ()
 
virtual void jacobianSetup ()
 
virtual void residualSetup ()
 
virtual void customSetup (const ExecFlagType &)
 
const ExecFlagEnumgetExecuteOnEnum () 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 UserObjectgetUserObjectBase (const std::string &param_name, bool is_dependency=true) const
 
const UserObjectgetUserObjectBaseByName (const UserObjectName &object_name, bool is_dependency=true) const
 
const std::vector< MooseVariableScalar *> & getCoupledMooseScalarVars ()
 
const std::set< TagID > & getScalarVariableCoupleableVectorTags () const
 
const std::set< TagID > & getScalarVariableCoupleableMatrixTags () const
 
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< SubdomainIDgetMaterialPropertyBlocks (const std::string &name)
 
std::vector< SubdomainName > getMaterialPropertyBlockNames (const std::string &name)
 
std::set< BoundaryIDgetMaterialPropertyBoundaryIDs (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)
 
bool isImplicit ()
 
Moose::StateArg determineState () 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 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)
 
MaterialBasegetMaterial (const std::string &name)
 
MaterialBasegetMaterial (const std::string &name)
 
MaterialBasegetMaterialByName (const std::string &name, bool no_warn=false)
 
MaterialBasegetMaterialByName (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)
 
const FunctiongetFunction (const std::string &name) const
 
const FunctiongetFunctionByName (const FunctionName &name) const
 
bool hasFunction (const std::string &param_name) const
 
bool hasFunctionByName (const FunctionName &name) 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 VectorPostprocessorValuegetVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueOldByName (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
 
T & getSampler (const std::string &name)
 
SamplergetSampler (const std::string &name)
 
T & getSamplerByName (const SamplerName &name)
 
SamplergetSamplerByName (const SamplerName &name)
 
virtual void meshChanged ()
 
virtual void meshDisplaced ()
 
PerfGraphperfGraph ()
 
const PostprocessorValuegetPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
virtual const PostprocessorValuegetPostprocessorValueByName (const PostprocessorName &name) const
 
virtual const PostprocessorValuegetPostprocessorValueByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOlderByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOlderByName (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 DistributiongetDistribution (const std::string &name) const
 
const T & getDistribution (const std::string &name) const
 
const DistributiongetDistribution (const std::string &name) const
 
const T & getDistribution (const std::string &name) const
 
const DistributiongetDistributionByName (const DistributionName &name) const
 
const T & getDistributionByName (const std::string &name) const
 
const DistributiongetDistributionByName (const DistributionName &name) const
 
const T & getDistributionByName (const std::string &name) const
 
const Parallel::Communicator & comm () const
 
processor_id_type n_processors () const
 
processor_id_type processor_id () const
 

Static Public Member Functions

static InputParameters validParams ()
 
static void sort (typename std::vector< T > &vector)
 
static void sortDFS (typename std::vector< T > &vector)
 
static void cyclicDependencyError (CyclicDependencyException< T2 > &e, const std::string &header)
 

Public Attributes

const ConsoleStream _console
 

Static Public Attributes

static constexpr PropertyValue::id_type default_property_id
 
static constexpr PropertyValue::id_type zero_property_id
 
static constexpr auto SYSTEM
 
static constexpr auto NAME
 

Protected Types

enum  NaNMessage { NAN_MESSAGE_NONE = 0, NAN_MESSAGE_WARNING = 1, NAN_MESSAGE_EXCEPTION = 2, NAN_MESSAGE_ERROR = 3 }
 

Protected Member Functions

virtual ADReal computeFlowArea (const std::vector< ADReal > &UL, const std::vector< ADReal > &UR) const
 Computes the flow area that is used in the numerical flux. More...
 
virtual std::vector< ADRealconvert1DInputTo3D (const std::vector< ADReal > &U_1d) const override
 Converts a 1D flux input vector to a 3D flux input vector. More...
 
virtual std::vector< ADRealconvert3DFluxTo1D (const std::vector< ADReal > &F_3d) const override
 Converts a 3D flux vector to a 1D flux vector. More...
 
virtual void transform3DFluxDirection (std::vector< ADReal > &F_3d, Real nLR_dot_d) const override
 Applies direction transformation to a 3D flux vector. More...
 
virtual void addPostprocessorDependencyHelper (const PostprocessorName &name) const override
 
virtual void addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &name) const override
 
virtual void addUserObjectDependencyHelper (const UserObject &uo) const override
 
void addReporterDependencyHelper (const ReporterName &reporter_name) override
 
const ReporterNamegetReporterName (const std::string &param_name) 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
 
const T & getMeshProperty (const std::string &data_name, const std::string &prefix)
 
const T & getMeshProperty (const std::string &data_name)
 
bool hasMeshProperty (const std::string &data_name, const std::string &prefix) const
 
bool hasMeshProperty (const std::string &data_name, const std::string &prefix) const
 
bool hasMeshProperty (const std::string &data_name) const
 
bool hasMeshProperty (const std::string &data_name) const
 
std::string meshPropertyName (const std::string &data_name) const
 
PerfID registerTimedSection (const std::string &section_name, const unsigned int level) const
 
PerfID registerTimedSection (const std::string &section_name, const unsigned int level, const std::string &live_message, const bool print_dots=true) const
 
std::string timedSectionName (const std::string &section_name) 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 VariableValuecoupledScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledScalarValue (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 VariableValuecoupledVectorTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValuecoupledMatrixTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValuecoupledScalarValueOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarValueOlder (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDu (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
const MooseVariableScalargetScalarVar (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 ()
 
const T & getReporterValue (const std::string &param_name, const std::size_t time_index=0)
 
const T & getReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)
 
const T & getReporterValue (const std::string &param_name, const std::size_t time_index=0)
 
const T & getReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)
 
const T & getReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)
 
const T & getReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)
 
const T & getReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)
 
const T & getReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
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)
 
Real getNaN () const
 Throws an error or returns a NaN with or without a warning, with a default message. More...
 
template<typename... Args>
Real getNaN (Args &&... args) const
 Throws an error or returns a NaN with or without a warning. More...
 
std::vector< RealgetNaNVector (const unsigned int &n) const
 Throws an error or returns NaNs with or without a warning, with a default message. More...
 
template<typename... Args>
std::vector< RealgetNaNVector (const unsigned int &n, Args &&... args) const
 Throws an error or returns NaNs with or without a warning. More...
 

Static Protected Member Functions

static std::string meshPropertyName (const std::string &data_name, const std::string &prefix)
 

Protected Attributes

const VaporMixtureFluidProperties_fp
 fluid properties user object More...
 
unsigned int _cached_flux_elem_id
 element ID of the cached flux values More...
 
unsigned int _cached_flux_side_id
 side ID of the cached flux values More...
 
std::vector< ADReal_FL_3d
 flux vector for the "left" cell for 3D More...
 
std::vector< ADReal_FR_3d
 flux vector for the "right" cell for 3D More...
 
std::vector< ADReal_FL_1d
 flux vector for the "left" cell for 1D More...
 
std::vector< ADReal_FR_1d
 flux vector for the "right" cell for 1D More...
 
unsigned int _last_region_index
 Index describing the region last entered, which is useful for testing and debugging. More...
 
SubProblem_subproblem
 
FEProblemBase_fe_problem
 
SystemBase_sys
 
const THREAD_ID _tid
 
Assembly_assembly
 
const Moose::CoordinateSystemType_coord_sys
 
const bool _duplicate_initial_execution
 
std::set< std::string > _depend_uo
 
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
 
MooseApp_restartable_app
 
const std::string _restartable_system_name
 
const THREAD_ID _restartable_tid
 
const bool _restartable_read_only
 
FEProblemBase_mci_feproblem
 
FEProblemBase_mdi_feproblem
 
MooseApp_pg_moose_app
 
const std::string _prefix
 
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
 
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
 
const Parallel::Communicator & _communicator
 
const MooseObject *const _moose_object
 
enum NaNMessage _emit_on_nan
 Raise mooseWarning or mooseError? More...
 

Static Protected Attributes

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

Detailed Description

Computes the numerical flux for FlowModelGasMix using the HLLC approximate Riemann solver.

Definition at line 21 of file NumericalFluxGasMixHLLC.h.

Member Enumeration Documentation

◆ NaNMessage

enum NaNInterface::NaNMessage
protectedinherited
Enumerator
NAN_MESSAGE_NONE 
NAN_MESSAGE_WARNING 
NAN_MESSAGE_EXCEPTION 
NAN_MESSAGE_ERROR 

Definition at line 30 of file NaNInterface.h.

Constructor & Destructor Documentation

◆ NumericalFluxGasMixHLLC()

NumericalFluxGasMixHLLC::NumericalFluxGasMixHLLC ( const InputParameters parameters)

Definition at line 29 of file NumericalFluxGasMixHLLC.C.

31  NaNInterface(this),
32  _fp(getUserObject<VaporMixtureFluidProperties>("fluid_properties"))
33 {
34 }
NaNInterface(const MooseObject *moose_object)
Definition: NaNInterface.C:35
NumericalFluxGasMixBase(const InputParameters &parameters)
const InputParameters & parameters() const
const VaporMixtureFluidProperties & _fp
fluid properties user object

Member Function Documentation

◆ calcFlux()

void NumericalFluxGasMixHLLC::calcFlux ( const std::vector< ADReal > &  UL_3d,
const std::vector< ADReal > &  UR_3d,
const RealVectorValue nLR,
const RealVectorValue t1,
const RealVectorValue t2,
std::vector< ADReal > &  FL,
std::vector< ADReal > &  FR 
) const
overridevirtual

Calculates the 3D flux vectors given "left" and "right" states.

This function is called only if the values are not already cached.

Parameters
[in]UL_3dVector of 3D flux inputs on the "left"
[in]UR_3dVector of 3D flux inputs on the "right"
[in]nLRDirection from "left" to "right"
[in]t11st tangent direction
[in]t22nd tangent direction
[out]FLFlux vector to be added to "left" side
[out]FRFlux vector to be added to "right" side

Implements NumericalFlux1D.

Definition at line 37 of file NumericalFluxGasMixHLLC.C.

44 {
45  // extract the conserved variables and area
46 
47  const ADReal xirhoAL = UL[THMGasMix3D::XIRHOA];
48  const ADReal rhoAL = UL[THMGasMix3D::RHOA];
49  const ADReal rhouAL = UL[THMGasMix3D::RHOUA];
50  const ADReal rhovAL = UL[THMGasMix3D::RHOVA];
51  const ADReal rhowAL = UL[THMGasMix3D::RHOWA];
52  const ADReal rhoEAL = UL[THMGasMix3D::RHOEA];
53  const ADReal AL = UL[THMGasMix3D::AREA];
54 
55  const ADReal xirhoAR = UR[THMGasMix3D::XIRHOA];
56  const ADReal rhoAR = UR[THMGasMix3D::RHOA];
57  const ADReal rhouAR = UR[THMGasMix3D::RHOUA];
58  const ADReal rhovAR = UR[THMGasMix3D::RHOVA];
59  const ADReal rhowAR = UR[THMGasMix3D::RHOWA];
60  const ADReal rhoEAR = UR[THMGasMix3D::RHOEA];
61  const ADReal AR = UR[THMGasMix3D::AREA];
62 
63  // compute the primitive variables
64 
65  const ADReal rhoL = rhoAL / AL;
66  const ADRealVectorValue uvecL(rhouAL / rhoAL, rhovAL / rhoAL, rhowAL / rhoAL);
67  const ADReal unL = uvecL * nLR;
68  const ADReal ut1L = uvecL * t1;
69  const ADReal ut2L = uvecL * t2;
70  const ADReal rhoEL = rhoEAL / AL;
71  const ADReal vL = 1.0 / rhoL;
72  const ADReal EL = rhoEAL / rhoAL;
73  const ADReal eL = EL - 0.5 * uvecL * uvecL;
74  const ADReal xiL = xirhoAL / rhoAL;
75  const ADReal pL = _fp.p_from_v_e(vL, eL, {xiL});
76  const ADReal TL = _fp.T_from_v_e(vL, eL, {xiL});
77  const ADReal cL = _fp.c_from_p_T(pL, TL, {xiL});
78 
79  const ADReal rhoR = rhoAR / AR;
80  const ADRealVectorValue uvecR(rhouAR / rhoAR, rhovAR / rhoAR, rhowAR / rhoAR);
81  const ADReal unR = uvecR * nLR;
82  const ADReal ut1R = uvecR * t1;
83  const ADReal ut2R = uvecR * t2;
84  const ADReal rhoER = rhoEAR / AR;
85  const ADReal vR = 1.0 / rhoR;
86  const ADReal ER = rhoEAR / rhoAR;
87  const ADReal eR = ER - 0.5 * uvecR * uvecR;
88  const ADReal xiR = xirhoAR / rhoAR;
89  const ADReal pR = _fp.p_from_v_e(vR, eR, {xiR});
90  const ADReal TR = _fp.T_from_v_e(vR, eR, {xiR});
91  const ADReal cR = _fp.c_from_p_T(pR, TR, {xiR});
92 
93  // compute wave speeds
94  const ADReal sL = std::min(unL - cL, unR - cR);
95  const ADReal sR = std::max(unL + cL, unR + cR);
96  const ADReal sm = (rhoR * unR * (sR - unR) - rhoL * unL * (sL - unL) + pL - pR) /
97  (rhoR * (sR - unR) - rhoL * (sL - unL));
98 
99  // compute Omega_L, Omega_R
100  const ADReal omegL = 1.0 / (sL - sm);
101  const ADReal omegR = 1.0 / (sR - sm);
102 
103  // compute p^*
104  const ADReal ps = rhoL * (sL - unL) * (sm - unL) + pL;
105 
106  // compute U_L^*, U_R^*
107 
108  const ADReal rhoLs = omegL * (sL - unL) * rhoL;
109  const ADReal rhounLs = omegL * ((sL - unL) * rhoL * unL + ps - pL);
110  const ADReal rhoELs = omegL * ((sL - unL) * rhoEL - pL * unL + ps * sm);
111 
112  const ADReal rhoRs = omegR * (sR - unR) * rhoR;
113  const ADReal rhounRs = omegR * ((sR - unR) * rhoR * unR + ps - pR);
114  const ADReal rhoERs = omegR * ((sR - unR) * rhoER - pR * unR + ps * sm);
115 
116  std::vector<ADReal> UL_1d(THMGasMix1D::N_FLUX_INPUTS);
122 
123  std::vector<ADReal> UR_1d(THMGasMix1D::N_FLUX_INPUTS);
129 
130  const ADReal A_flow = computeFlowArea(UL_1d, UR_1d);
131 
132  // compute the fluxes
133  FL.resize(THMGasMix3D::N_FLUX_OUTPUTS);
134  if (sL > 0.0)
135  {
136  FL[THMGasMix3D::SPECIES] = unL * xiL * rhoL * A_flow;
137  FL[THMGasMix3D::MASS] = unL * rhoL * A_flow;
138  FL[THMGasMix3D::MOM_NORM] = (unL * rhoL * unL + pL) * A_flow;
139  FL[THMGasMix3D::MOM_TAN1] = rhoL * unL * ut1L * A_flow;
140  FL[THMGasMix3D::MOM_TAN2] = rhoL * unL * ut2L * A_flow;
141  FL[THMGasMix3D::ENERGY] = unL * (rhoEL + pL) * A_flow;
142 
143  _last_region_index = 0;
144  }
145  else if (sL <= 0.0 && sm > 0.0)
146  {
147  FL[THMGasMix3D::SPECIES] = sm * xiL * rhoLs * A_flow;
148  FL[THMGasMix3D::MASS] = sm * rhoLs * A_flow;
149  FL[THMGasMix3D::MOM_NORM] = (sm * rhounLs + ps) * A_flow;
150  FL[THMGasMix3D::MOM_TAN1] = rhounLs * ut1L * A_flow;
151  FL[THMGasMix3D::MOM_TAN2] = rhounLs * ut2L * A_flow;
152  FL[THMGasMix3D::ENERGY] = sm * (rhoELs + ps) * A_flow;
153 
154  _last_region_index = 1;
155  }
156  else if (sm <= 0.0 && sR >= 0.0)
157  {
158  FL[THMGasMix3D::SPECIES] = sm * xiR * rhoRs * A_flow;
159  FL[THMGasMix3D::MASS] = sm * rhoRs * A_flow;
160  FL[THMGasMix3D::MOM_NORM] = (sm * rhounRs + ps) * A_flow;
161  FL[THMGasMix3D::MOM_TAN1] = rhounRs * ut1R * A_flow;
162  FL[THMGasMix3D::MOM_TAN2] = rhounRs * ut2R * A_flow;
163  FL[THMGasMix3D::ENERGY] = sm * (rhoERs + ps) * A_flow;
164 
165  _last_region_index = 2;
166  }
167  else if (sR < 0.0)
168  {
169  FL[THMGasMix3D::SPECIES] = unR * xiR * rhoR * A_flow;
170  FL[THMGasMix3D::MASS] = unR * rhoR * A_flow;
171  FL[THMGasMix3D::MOM_NORM] = (unR * rhoR * unR + pR) * A_flow;
172  FL[THMGasMix3D::MOM_TAN1] = rhoR * unR * ut1R * A_flow;
173  FL[THMGasMix3D::MOM_TAN2] = rhoR * unR * ut2R * A_flow;
174  FL[THMGasMix3D::ENERGY] = unR * (rhoER + pR) * A_flow;
175 
176  _last_region_index = 3;
177  }
178  else
179  std::fill(FL.begin(), FL.end(), getNaN());
180 
181  FR = FL;
182 
183  const ADReal A_wall_L = AL - A_flow;
184  FL[THMGasMix3D::MOM_NORM] += pL * A_wall_L;
185 
186  const ADReal A_wall_R = AR - A_flow;
187  FR[THMGasMix3D::MOM_NORM] += pR * A_wall_R;
188 }
Real getNaN() const
Throws an error or returns a NaN with or without a warning, with a default message.
Definition: NaNInterface.h:46
DualNumber< Real, DNDerivativeType, true > ADReal
unsigned int _last_region_index
Index describing the region last entered, which is useful for testing and debugging.
virtual ADReal computeFlowArea(const std::vector< ADReal > &UL, const std::vector< ADReal > &UR) const
Computes the flow area that is used in the numerical flux.
static const unsigned int N_FLUX_INPUTS
Number of numerical flux function inputs.
const VaporMixtureFluidProperties & _fp
fluid properties user object
static const std::string cL
Definition: NS.h:111
static const unsigned int N_FLUX_OUTPUTS
Number of numerical flux function outputs for 3D.

◆ computeFlowArea()

ADReal NumericalFluxGasMixHLLC::computeFlowArea ( const std::vector< ADReal > &  UL,
const std::vector< ADReal > &  UR 
) const
protectedvirtual

Computes the flow area that is used in the numerical flux.

Definition at line 191 of file NumericalFluxGasMixHLLC.C.

Referenced by calcFlux().

193 {
194  return std::min(UL[THMGasMix1D::AREA], UR[THMGasMix1D::AREA]);
195 }

◆ convert1DInputTo3D()

std::vector< ADReal > NumericalFluxGasMixBase::convert1DInputTo3D ( const std::vector< ADReal > &  U_1d) const
overrideprotectedvirtualinherited

Converts a 1D flux input vector to a 3D flux input vector.

Parameters
[in]U_1d1D flux input vector

Implements NumericalFlux1D.

Definition at line 26 of file NumericalFluxGasMixBase.C.

◆ convert3DFluxTo1D()

std::vector< ADReal > NumericalFluxGasMixBase::convert3DFluxTo1D ( const std::vector< ADReal > &  F_3d) const
overrideprotectedvirtualinherited

Converts a 3D flux vector to a 1D flux vector.

Parameters
[in]F_3d3D flux vector

Implements NumericalFlux1D.

Definition at line 41 of file NumericalFluxGasMixBase.C.

◆ execute()

void NumericalFlux1D::execute ( )
overridevirtualinherited

Implements ThreadedGeneralUserObject.

Definition at line 38 of file NumericalFlux1D.C.

39 {
40 }

◆ finalize()

void NumericalFlux1D::finalize ( )
overridevirtualinherited

Implements ThreadedGeneralUserObject.

Definition at line 43 of file NumericalFlux1D.C.

44 {
45 }

◆ getFlux()

const std::vector< ADReal > & NumericalFlux1D::getFlux ( const unsigned int  iside,
const dof_id_type  ielem,
bool  res_side_is_left,
const std::vector< ADReal > &  UL_1d,
const std::vector< ADReal > &  UR_1d,
Real  nLR_dot_d 
) const
virtualinherited

Gets the 1D flux vector for an element/side combination.

If the element/side is cached, the cached values are used. Otherwise the values are computed, cached, and returned.

Parameters
[in]isidelocal index of current side
[in]ielemglobal index of the current element
[in]res_side_is_leftgetting flux on the left ("elem") side?
[in]UL_1dvector of 1D flux inputs on the "left"
[in]UR_1dvector of 1D flux inputs on the "right"
[in]nLR_dot_dDot product of direction from "left" to "right" with the flow channel direction
Returns
flux vector for an element/side combination

Definition at line 53 of file NumericalFlux1D.C.

Referenced by ADBoundaryFlux3EqnGhostBase::calcFlux(), BoundaryFluxGasMixGhostBase::calcFlux(), NumericalFluxGasMixDGKernel::computeQpResidual(), ADNumericalFlux3EqnDGKernel::computeQpResidual(), ADGateValve1PhaseUserObject::finalize(), and ADJunctionOneToOne1PhaseUserObject::finalize().

59 {
60  if (_cached_flux_elem_id != ielem || _cached_flux_side_id != iside)
61  {
62  _cached_flux_elem_id = ielem;
63  _cached_flux_side_id = iside;
64 
65  const auto UL_3d = convert1DInputTo3D(UL_1d);
66  const auto UR_3d = convert1DInputTo3D(UR_1d);
67 
68  const RealVectorValue nLR(nLR_dot_d, 0, 0);
69  RealVectorValue t1, t2;
71 
72  calcFlux(UL_3d, UR_3d, nLR, t1, t2, _FL_3d, _FR_3d);
73 
74  transform3DFluxDirection(_FL_3d, nLR_dot_d);
75  transform3DFluxDirection(_FR_3d, nLR_dot_d);
76 
79  }
80 
81  if (res_side_is_left)
82  return _FL_1d;
83  else
84  return _FR_1d;
85 }
std::vector< ADReal > _FR_3d
flux vector for the "right" cell for 3D
std::vector< ADReal > _FL_3d
flux vector for the "left" cell for 3D
unsigned int _cached_flux_elem_id
element ID of the cached flux values
virtual std::vector< ADReal > convert1DInputTo3D(const std::vector< ADReal > &U_1d) const =0
Converts a 1D flux input vector to a 3D flux input vector.
std::vector< ADReal > _FR_1d
flux vector for the "right" cell for 1D
virtual void transform3DFluxDirection(std::vector< ADReal > &F_3d, Real nLR_dot_d) const =0
Applies direction transformation to a 3D flux vector.
virtual std::vector< ADReal > convert3DFluxTo1D(const std::vector< ADReal > &F_3d) const =0
Converts a 3D flux vector to a 1D flux vector.
std::vector< ADReal > _FL_1d
flux vector for the "left" cell for 1D
void computeOrthogonalDirections(const RealVectorValue &n_unnormalized, RealVectorValue &t1, RealVectorValue &t2)
Computes two unit vectors orthogonal to the given vector.
Definition: THMUtils.C:22
virtual void calcFlux(const std::vector< ADReal > &UL_3d, const std::vector< ADReal > &UR_3d, const RealVectorValue &nLR, const RealVectorValue &t1, const RealVectorValue &t2, std::vector< ADReal > &FL, std::vector< ADReal > &FR) const =0
Calculates the 3D flux vectors given "left" and "right" states.
unsigned int _cached_flux_side_id
side ID of the cached flux values

◆ getFlux3D()

const std::vector< ADReal > & NumericalFlux1D::getFlux3D ( const unsigned int  iside,
const dof_id_type  ielem,
const std::vector< ADReal > &  UL_3d,
const std::vector< ADReal > &  UR_3d,
const RealVectorValue nLR,
const RealVectorValue t1,
const RealVectorValue t2 
) const
virtualinherited

Gets the 3D flux vector for an element/side combination on the left side.

If the element/side is cached, the cached values are used. Otherwise the values are computed, cached, and returned.

Parameters
[in]isidelocal index of current side
[in]ielemglobal index of the current element
[in]UL_3dvector of 3D flux inputs on the "left"
[in]UR_3dvector of 3D flux inputs on the "right"
[in]nLRDirection from "left" to "right"
[in]t11st tangent direction
[in]t22nd tangent direction
Returns
flux vector for an element/side combination

Definition at line 88 of file NumericalFlux1D.C.

95 {
96  if (_cached_flux_elem_id != ielem || _cached_flux_side_id != iside)
97  {
98  _cached_flux_elem_id = ielem;
99  _cached_flux_side_id = iside;
100 
101  calcFlux(UL_3d, UR_3d, nLR, t1, t2, _FL_3d, _FR_3d);
102  }
103 
104  return _FL_3d;
105 }
std::vector< ADReal > _FR_3d
flux vector for the "right" cell for 3D
std::vector< ADReal > _FL_3d
flux vector for the "left" cell for 3D
unsigned int _cached_flux_elem_id
element ID of the cached flux values
virtual void calcFlux(const std::vector< ADReal > &UL_3d, const std::vector< ADReal > &UR_3d, const RealVectorValue &nLR, const RealVectorValue &t1, const RealVectorValue &t2, std::vector< ADReal > &FL, std::vector< ADReal > &FR) const =0
Calculates the 3D flux vectors given "left" and "right" states.
unsigned int _cached_flux_side_id
side ID of the cached flux values

◆ getLastRegionIndex()

unsigned int NumericalFlux1D::getLastRegionIndex ( ) const
inlineinherited

Returns the index of the region last entered.

Here "region" refers to a code path taken. For some fluxes, such as centered fluxes, there is just a single code path, but for others, such as those using an approximate Riemann solver, there are multiple. Riemann solvers have "regions" defined by the characteristic waves.

Definition at line 106 of file NumericalFlux1D.h.

106 { return _last_region_index; }
unsigned int _last_region_index
Index describing the region last entered, which is useful for testing and debugging.

◆ getNaN() [1/2]

Real NaNInterface::getNaN ( ) const
inlineprotectedinherited

◆ getNaN() [2/2]

template<typename... Args>
Real NaNInterface::getNaN ( Args &&...  args) const
inlineprotectedinherited

Throws an error or returns a NaN with or without a warning.

Definition at line 62 of file NaNInterface.h.

63  {
64  switch (_emit_on_nan)
65  {
66  case (NAN_MESSAGE_WARNING):
67  mooseWarning(_moose_object->name(), ": ", std::forward<Args>(args)...);
68  break;
69  case (NAN_MESSAGE_EXCEPTION):
70  mooseException(_moose_object->name(), ": ", std::forward<Args>(args)...);
71  break;
72  case (NAN_MESSAGE_ERROR):
73  mooseError(_moose_object->name(), ": ", std::forward<Args>(args)...);
74  break;
75  default:
76  break;
77  }
78  // return a quiet NaN
79  return std::nan("");
80  }
void mooseError(Args &&... args)
void mooseWarning(Args &&... args)
virtual const std::string & name() const
enum NaNMessage _emit_on_nan
Raise mooseWarning or mooseError?
Definition: NaNInterface.h:41
const MooseObject *const _moose_object
Definition: NaNInterface.h:38

◆ getNaNVector() [1/2]

std::vector<Real> NaNInterface::getNaNVector ( const unsigned int n) const
inlineprotectedinherited

Throws an error or returns NaNs with or without a warning, with a default message.

Parameters
[in]nVector size

Definition at line 53 of file NaNInterface.h.

Referenced by NaNInterfaceTestFluidProperties::returnNaNVector().

54  {
55  return getNaNVector(n, "A NaN was produced.");
56  }
std::vector< Real > getNaNVector(const unsigned int &n) const
Throws an error or returns NaNs with or without a warning, with a default message.
Definition: NaNInterface.h:53

◆ getNaNVector() [2/2]

template<typename... Args>
std::vector<Real> NaNInterface::getNaNVector ( const unsigned int n,
Args &&...  args 
) const
inlineprotectedinherited

Throws an error or returns NaNs with or without a warning.

Parameters
[in]nVector size

Definition at line 88 of file NaNInterface.h.

89  {
90  switch (_emit_on_nan)
91  {
92  case (NAN_MESSAGE_WARNING):
93  mooseWarning(_moose_object->name(), ": ", std::forward<Args>(args)...);
94  break;
95  case (NAN_MESSAGE_EXCEPTION):
96  mooseException(_moose_object->name(), ": ", std::forward<Args>(args)...);
97  break;
98  case (NAN_MESSAGE_ERROR):
99  mooseError(_moose_object->name(), ": ", std::forward<Args>(args)...);
100  break;
101  default:
102  break;
103  }
104  // return quiet NaNs
105  return std::vector<Real>(n, std::nan(""));
106  }
void mooseError(Args &&... args)
void mooseWarning(Args &&... args)
virtual const std::string & name() const
enum NaNMessage _emit_on_nan
Raise mooseWarning or mooseError?
Definition: NaNInterface.h:41
const MooseObject *const _moose_object
Definition: NaNInterface.h:38

◆ getNumberOfRegions()

virtual unsigned int NumericalFluxGasMixHLLC::getNumberOfRegions ( ) const
inlineoverridevirtual

Returns the total possible number of regions.

Here "region" refers to a code path taken. For some fluxes, such as centered fluxes, there is just a single code path, but for others, such as those using an approximate Riemann solver, there are multiple. Riemann solvers have "regions" defined by the characteristic waves.

Implements NumericalFlux1D.

Definition at line 36 of file NumericalFluxGasMixHLLC.h.

36 { return 4; }

◆ initialize()

void NumericalFlux1D::initialize ( )
overridevirtualinherited

Implements ThreadedGeneralUserObject.

Definition at line 31 of file NumericalFlux1D.C.

32 {
35 }
const unsigned int invalid_uint
unsigned int _cached_flux_elem_id
element ID of the cached flux values
unsigned int _cached_flux_side_id
side ID of the cached flux values

◆ threadJoin()

void NumericalFlux1D::threadJoin ( const UserObject )
overridevirtualinherited

Reimplemented from ThreadedGeneralUserObject.

Definition at line 48 of file NumericalFlux1D.C.

49 {
50 }

◆ transform3DFluxDirection()

void NumericalFluxGasMixBase::transform3DFluxDirection ( std::vector< ADReal > &  F_3d,
Real  nLR_dot_d 
) const
overrideprotectedvirtualinherited

Applies direction transformation to a 3D flux vector.

Parameters
[in,out]F_3d3D flux vector
[in]nLR_dot_dDot product of direction from "left" to "right" with the flow channel direction

Implements NumericalFlux1D.

Definition at line 53 of file NumericalFluxGasMixBase.C.

54 {
55  F_3d[THMGasMix3D::SPECIES] *= nLR_dot_d;
56  F_3d[THMGasMix3D::MASS] *= nLR_dot_d;
57  F_3d[THMGasMix3D::ENERGY] *= nLR_dot_d;
58 }

◆ validParams()

InputParameters NumericalFluxGasMixHLLC::validParams ( )
static

Definition at line 18 of file NumericalFluxGasMixHLLC.C.

19 {
21  params += NaNInterface::validParams();
22  params.addRequiredParam<UserObjectName>("fluid_properties",
23  "Name for fluid properties user object");
24  params.addClassDescription("Computes internal side flux for the 1-D, 1-phase, variable-area "
25  "Euler equations using the HLLC approximate Riemann solver.");
26  return params;
27 }
void addRequiredParam(const std::string &name, const std::string &doc_string)
static InputParameters validParams()
void addClassDescription(const std::string &doc_string)
static InputParameters validParams()
Definition: NaNInterface.C:15

Member Data Documentation

◆ _cached_flux_elem_id

unsigned int NumericalFlux1D::_cached_flux_elem_id
mutableprotectedinherited

element ID of the cached flux values

Definition at line 143 of file NumericalFlux1D.h.

Referenced by NumericalFlux1D::getFlux(), NumericalFlux1D::getFlux3D(), and NumericalFlux1D::initialize().

◆ _cached_flux_side_id

unsigned int NumericalFlux1D::_cached_flux_side_id
mutableprotectedinherited

side ID of the cached flux values

Definition at line 145 of file NumericalFlux1D.h.

Referenced by NumericalFlux1D::getFlux(), NumericalFlux1D::getFlux3D(), and NumericalFlux1D::initialize().

◆ _emit_on_nan

enum NaNMessage NaNInterface::_emit_on_nan
protectedinherited

Raise mooseWarning or mooseError?

Definition at line 41 of file NaNInterface.h.

Referenced by NaNInterface::getNaN(), and NaNInterface::getNaNVector().

◆ _FL_1d

std::vector<ADReal> NumericalFlux1D::_FL_1d
mutableprotectedinherited

flux vector for the "left" cell for 1D

Definition at line 152 of file NumericalFlux1D.h.

Referenced by NumericalFlux1D::getFlux().

◆ _FL_3d

std::vector<ADReal> NumericalFlux1D::_FL_3d
mutableprotectedinherited

flux vector for the "left" cell for 3D

Definition at line 148 of file NumericalFlux1D.h.

Referenced by NumericalFlux1D::getFlux(), and NumericalFlux1D::getFlux3D().

◆ _fp

const VaporMixtureFluidProperties& NumericalFluxGasMixHLLC::_fp
protected

fluid properties user object

Definition at line 46 of file NumericalFluxGasMixHLLC.h.

Referenced by calcFlux().

◆ _FR_1d

std::vector<ADReal> NumericalFlux1D::_FR_1d
mutableprotectedinherited

flux vector for the "right" cell for 1D

Definition at line 154 of file NumericalFlux1D.h.

Referenced by NumericalFlux1D::getFlux().

◆ _FR_3d

std::vector<ADReal> NumericalFlux1D::_FR_3d
mutableprotectedinherited

flux vector for the "right" cell for 3D

Definition at line 150 of file NumericalFlux1D.h.

Referenced by NumericalFlux1D::getFlux(), and NumericalFlux1D::getFlux3D().

◆ _last_region_index

unsigned int NumericalFlux1D::_last_region_index
mutableprotectedinherited

Index describing the region last entered, which is useful for testing and debugging.

Definition at line 157 of file NumericalFlux1D.h.

Referenced by calcFlux(), ADNumericalFlux3EqnHLLC::calcFlux(), and NumericalFlux1D::getLastRegionIndex().

◆ _moose_object

const MooseObject* const NaNInterface::_moose_object
protectedinherited

Definition at line 38 of file NaNInterface.h.

Referenced by NaNInterface::getNaN(), and NaNInterface::getNaNVector().


The documentation for this class was generated from the following files: