Upwind numerical flux scheme for the advection equation using a cell-centered finite volume method. More...

#include <AEFVUpwindInternalSideFlux.h>

Inheritance diagram for AEFVUpwindInternalSideFlux:

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	AEFVUpwindInternalSideFlux (const InputParameters &parameters)

virtual	~AEFVUpwindInternalSideFlux ()

virtual void	calcFlux (unsigned int iside, dof_id_type ielem, dof_id_type ineig, const std::vector< Real > &uvec1, const std::vector< Real > &uvec2, const RealVectorValue &dwave, std::vector< Real > &flux) const override
	Solve the Riemann problem. More...

virtual void	calcJacobian (unsigned int iside, dof_id_type ielem, dof_id_type ineig, const std::vector< Real > &uvec1, const std::vector< Real > &uvec2, const RealVectorValue &dwave, DenseMatrix< Real > &jac1, DenseMatrix< Real > &jac2) const override
	Compute the Jacobian matrix. More...

virtual void	execute () override

virtual void	initialize () override

virtual void	finalize () override

virtual void	threadJoin (const UserObject &) override

virtual const std::vector< Real > &	getFlux (unsigned int iside, dof_id_type ielem, dof_id_type ineig, const std::vector< Real > &uvec1, const std::vector< Real > &uvec2, const RealVectorValue &dwave) const
	Get the flux vector. More...

virtual const DenseMatrix< Real > &	getJacobian (Moose::DGResidualType type, unsigned int iside, dof_id_type ielem, dof_id_type ineig, const std::vector< Real > &uvec1, const std::vector< Real > &uvec2, const RealVectorValue &dwave) const
	Get the Jacobian matrix. More...

virtual void	subdomainSetup () override

bool	needThreadedCopy () const override final

SubProblem &	getSubProblem () 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)

UserObject *	primaryThreadCopy ()

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< 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	initialSetup ()

virtual void	timestepSetup ()

virtual void	jacobianSetup ()

virtual void	residualSetup ()

virtual void	customSetup (const ExecFlagType &)

const ExecFlagEnum &	getExecuteOnEnum () 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

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

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)

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)

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

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

T &	getSampler (const std::string &name)

Sampler &	getSampler (const std::string &name)

T &	getSamplerByName (const SamplerName &name)

Sampler &	getSamplerByName (const SamplerName &name)

virtual void	meshChanged ()

virtual void	meshDisplaced ()

PerfGraph &	perfGraph ()

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 Distribution &	getDistribution (const std::string &name) const

const T &	getDistribution (const std::string &name) const

const Distribution &	getDistribution (const std::string &name) const

const T &	getDistribution (const std::string &name) const

const Distribution &	getDistributionByName (const DistributionName &name) const

const T &	getDistributionByName (const std::string &name) const

const Distribution &	getDistributionByName (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 Member Functions
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 ReporterName &	getReporterName (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 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 ()

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)

Static Protected Member Functions
static std::string	meshPropertyName (const std::string &data_name, const std::string &prefix)

Protected Attributes
const Real	_velocity
	advective velocity More...

unsigned int	_cached_flux_elem_id
	element ID of the cached flux values More...

unsigned int	_cached_flux_neig_id
	neighbor element ID of the cached flux values More...

unsigned int	_cached_jacobian_elem_id
	element ID of the cached Jacobian values More...

unsigned int	_cached_jacobian_neig_id
	neighbor element ID of the cached Jacobian values More...

std::vector< Real >	_flux
	flux vector of this side More...

DenseMatrix< Real >	_jac1
	Jacobian matrix contribution to the "left" cell. More...

DenseMatrix< Real >	_jac2
	Jacobian matrix contribution to the "right" cell. 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

Static Protected Attributes
static const std::string	_interpolated_old

static const std::string	_interpolated_older

Detailed Description

Upwind numerical flux scheme for the advection equation using a cell-centered finite volume method.

Definition at line 21 of file AEFVUpwindInternalSideFlux.h.

Constructor & Destructor Documentation

◆ AEFVUpwindInternalSideFlux()

AEFVUpwindInternalSideFlux::AEFVUpwindInternalSideFlux ( const InputParameters & parameters )

Definition at line 24 of file AEFVUpwindInternalSideFlux.C.

   : InternalSideFluxBase(parameters), _velocity(getParam<Real>("velocity"))
 {
 }

◆ ~AEFVUpwindInternalSideFlux()

AEFVUpwindInternalSideFlux::~AEFVUpwindInternalSideFlux ( )

virtual

Definition at line 29 of file AEFVUpwindInternalSideFlux.C.

29 {}

Member Function Documentation

◆ calcFlux()

void AEFVUpwindInternalSideFlux::calcFlux	(	unsigned int	iside,
		dof_id_type	ielem,
		dof_id_type	ineig,
		const std::vector< Real > &	uvec1,
		const std::vector< Real > &	uvec2,
		const RealVectorValue &	dwave,
		std::vector< Real > &	flux
	)		const

overridevirtual

Solve the Riemann problem.

Parameters

[in]	iside	local index of current side
[in]	ielem	global index of the current element
[in]	ineig	global index of the neighbor element
[in]	uvec1	vector of variables on the "left"
[in]	uvec2	vector of variables on the "right"
[in]	dwave	vector of unit normal
[out]	flux	flux vector across the side

Implements InternalSideFluxBase.

Definition at line 32 of file AEFVUpwindInternalSideFlux.C.

 {
   mooseAssert(uvec1.size() == 1, "Invalid size for uvec1. Must be single variable coupling.");
   mooseAssert(uvec2.size() == 1, "Invalid size for uvec1. Must be single variable coupling.");
 
   // assign the size of flux vector, e.g. = 1 for the advection equation
   flux.resize(1);
 
   // assume a constant velocity on the left
   RealVectorValue uadv1(_velocity, 0.0, 0.0);
 
   // assume a constant velocity on the right
   RealVectorValue uadv2(_velocity, 0.0, 0.0);
 
   // normal velocity on the left and right
   Real vdon1 = uadv1 * dwave;
   Real vdon2 = uadv2 * dwave;
 
   // calculate the so-called a^plus and a^minus
   Real aplus = 0.5 * (vdon1 + std::abs(vdon1));
   Real amins = 0.5 * (vdon2 - std::abs(vdon2));
 
   // finally calculate the flux
   flux[0] = aplus * uvec1[0] + amins * uvec2[0];
 }

◆ calcJacobian()

void AEFVUpwindInternalSideFlux::calcJacobian	(	unsigned int	iside,
		dof_id_type	ielem,
		dof_id_type	ineig,
		const std::vector< Real > &	uvec1,
		const std::vector< Real > &	uvec2,
		const RealVectorValue &	dwave,
		DenseMatrix< Real > &	jac1,
		DenseMatrix< Real > &	jac2
	)		const

overridevirtual

Compute the Jacobian matrix.

Parameters

[in]	iside	local index of current side
[in]	ielem	global index of the current element
[in]	ineig	global index of the neighbor element
[in]	uvec1	vector of variables on the "left"
[in]	uvec2	vector of variables on the "right"
[in]	dwave	vector of unit normal
[out]	jac1	Jacobian matrix contribution to the "left" cell
[out]	jac2	Jacobian matrix contribution to the "right" cell

Implements InternalSideFluxBase.

Definition at line 65 of file AEFVUpwindInternalSideFlux.C.

 {
   mooseAssert(uvec1.size() == 1, "Invalid size for uvec1. Must be single variable coupling.");
   mooseAssert(uvec2.size() == 1, "Invalid size for uvec1. Must be single variable coupling.");
 
   // assign the size of Jacobian matrix, e.g. = (1, 1) for the advection equation
   jac1.resize(1, 1);
   jac2.resize(1, 1);
 
   // assume a constant velocity on the left
   RealVectorValue uadv1(_velocity, 0.0, 0.0);
 
   // assume a constant velocity on the right
   RealVectorValue uadv2(_velocity, 0.0, 0.0);
 
   // normal velocity on the left and right
   Real vdon1 = uadv1 * dwave;
   Real vdon2 = uadv2 * dwave;
 
   // calculate the so-called a^plus and a^minus
   Real aplus = 0.5 * (vdon1 + std::abs(vdon1));
   Real amins = 0.5 * (vdon2 - std::abs(vdon2));
 
   // finally calculate the Jacobian matrix
   jac1(0, 0) = aplus;
   jac2(0, 0) = amins;
 }

◆ execute()

void InternalSideFluxBase::execute ( )

overridevirtualinherited

Implements ThreadedGeneralUserObject.

Definition at line 39 of file InternalSideFluxBase.C.

40 {

41 }

◆ finalize()

void InternalSideFluxBase::finalize ( )

overridevirtualinherited

Implements ThreadedGeneralUserObject.

Definition at line 44 of file InternalSideFluxBase.C.

45 {

46 }

◆ getFlux()

const std::vector< Real > & InternalSideFluxBase::getFlux	(	unsigned int	iside,
		dof_id_type	ielem,
		dof_id_type	ineig,
		const std::vector< Real > &	uvec1,
		const std::vector< Real > &	uvec2,
		const RealVectorValue &	dwave
	)		const

virtualinherited

Get the flux vector.

Parameters

[in]	iside	local index of current side
[in]	ielem	global index of the current element
[in]	ineig	global index of the neighbor element
[in]	uvec1	vector of variables on the "left"
[in]	uvec2	vector of variables on the "right"
[in]	dwave	vector of unit normal

Definition at line 54 of file InternalSideFluxBase.C.

Referenced by AEFVKernel::computeQpResidual().

 {
   if (_cached_flux_elem_id != ielem || _cached_flux_neig_id != ineig)
   {
     _cached_flux_elem_id = ielem;
     _cached_flux_neig_id = ineig;
 
     calcFlux(iside, ielem, ineig, uvec1, uvec2, dwave, _flux);
   }
   return _flux;
 }

◆ getJacobian()

const DenseMatrix< Real > & InternalSideFluxBase::getJacobian	(	Moose::DGResidualType	type,
		unsigned int	iside,
		dof_id_type	ielem,
		dof_id_type	ineig,
		const std::vector< Real > &	uvec1,
		const std::vector< Real > &	uvec2,
		const RealVectorValue &	dwave
	)		const

virtualinherited

Get the Jacobian matrix.

Parameters

[in]	iside	local index of current side
[in]	ielem	global index of the current element
[in]	ineig	global index of the neighbor element
[in]	uvec1	vector of variables on the "left"
[in]	uvec2	vector of variables on the "right"
[in]	dwave	vector of unit normal

Definition at line 72 of file InternalSideFluxBase.C.

Referenced by AEFVKernel::computeQpJacobian().

 {
   if (_cached_jacobian_elem_id != ielem || _cached_jacobian_neig_id != ineig)
   {
     _cached_jacobian_elem_id = ielem;
     _cached_jacobian_neig_id = ineig;
 
     calcJacobian(iside, ielem, ineig, uvec1, uvec2, dwave, _jac1, _jac2);
   }
 
   if (type == Moose::Element)
     return _jac1;
   else
     return _jac2;
 }

◆ initialize()

void InternalSideFluxBase::initialize ( )

overridevirtualinherited

Implements ThreadedGeneralUserObject.

Definition at line 30 of file InternalSideFluxBase.C.

 {
   _cached_flux_elem_id = libMesh::invalid_uint;
   _cached_flux_neig_id = libMesh::invalid_uint;
   _cached_jacobian_elem_id = libMesh::invalid_uint;
   _cached_jacobian_neig_id = libMesh::invalid_uint;
 }

◆ threadJoin()

void InternalSideFluxBase::threadJoin ( const UserObject & )

overridevirtualinherited

Reimplemented from ThreadedGeneralUserObject.

Definition at line 49 of file InternalSideFluxBase.C.

50 {

51 }

◆ validParams()

InputParameters AEFVUpwindInternalSideFlux::validParams ( )

static

Definition at line 15 of file AEFVUpwindInternalSideFlux.C.

 {
   InputParameters params = InternalSideFluxBase::validParams();
   params.addClassDescription("Upwind numerical flux scheme for the advection equation using a "
                              "cell-centered finite volume method.");
   params.addParam<Real>("velocity", 1.0, "Advective velocity");
   return params;
 }

Member Data Documentation

◆ _cached_flux_elem_id

unsigned int InternalSideFluxBase::_cached_flux_elem_id

mutableprotectedinherited

element ID of the cached flux values

Definition at line 112 of file InternalSideFluxBase.h.

Referenced by InternalSideFluxBase::getFlux(), and InternalSideFluxBase::initialize().

◆ _cached_flux_neig_id

unsigned int InternalSideFluxBase::_cached_flux_neig_id

mutableprotectedinherited

neighbor element ID of the cached flux values

Definition at line 114 of file InternalSideFluxBase.h.

Referenced by InternalSideFluxBase::getFlux(), and InternalSideFluxBase::initialize().

◆ _cached_jacobian_elem_id

unsigned int InternalSideFluxBase::_cached_jacobian_elem_id

mutableprotectedinherited

element ID of the cached Jacobian values

Definition at line 117 of file InternalSideFluxBase.h.

Referenced by InternalSideFluxBase::getJacobian(), and InternalSideFluxBase::initialize().

◆ _cached_jacobian_neig_id

unsigned int InternalSideFluxBase::_cached_jacobian_neig_id

mutableprotectedinherited

neighbor element ID of the cached Jacobian values

Definition at line 119 of file InternalSideFluxBase.h.

Referenced by InternalSideFluxBase::getJacobian(), and InternalSideFluxBase::initialize().

◆ _flux

std::vector<Real> InternalSideFluxBase::_flux

mutableprotectedinherited

flux vector of this side

Definition at line 122 of file InternalSideFluxBase.h.

Referenced by InternalSideFluxBase::getFlux().

◆ _jac1

DenseMatrix<Real> InternalSideFluxBase::_jac1

mutableprotectedinherited

Jacobian matrix contribution to the "left" cell.

Definition at line 124 of file InternalSideFluxBase.h.

Referenced by InternalSideFluxBase::getJacobian().

◆ _jac2

DenseMatrix<Real> InternalSideFluxBase::_jac2

mutableprotectedinherited

Jacobian matrix contribution to the "right" cell.

Definition at line 126 of file InternalSideFluxBase.h.

Referenced by InternalSideFluxBase::getJacobian().

◆ _velocity

const Real AEFVUpwindInternalSideFlux::_velocity

protected

advective velocity

Definition at line 48 of file AEFVUpwindInternalSideFlux.h.

Referenced by calcFlux(), and calcJacobian().

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

rdg/include/userobjects/AEFVUpwindInternalSideFlux.h
rdg/src/userobjects/AEFVUpwindInternalSideFlux.C

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Static Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ AEFVUpwindInternalSideFlux()

◆ ~AEFVUpwindInternalSideFlux()

Member Function Documentation

◆ calcFlux()

◆ calcJacobian()

◆ execute()

◆ finalize()

◆ getFlux()

◆ getJacobian()

◆ initialize()

◆ threadJoin()

◆ validParams()

Member Data Documentation

◆ _cached_flux_elem_id

◆ _cached_flux_neig_id

◆ _cached_jacobian_elem_id

◆ _cached_jacobian_neig_id

◆ _flux

◆ _jac1

◆ _jac2

◆ _velocity