Sets up a Jacobian test for 1-phase rDG. More...

#include <JacobianTest1PhaseRDGAction.h>

Inheritance diagram for JacobianTest1PhaseRDGAction:

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	JacobianTest1PhaseRDGAction (const InputParameters &params)

virtual void	act ()

void	timedAct ()

virtual void	addRelationshipManagers (Moose::RelationshipManagerType when_type)

MooseObjectName	uniqueActionName () const

const std::string &	specificTaskName () const

const std::set< std::string > &	getAllTasks () const

void	appendTask (const std::string &task)

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

PerfGraph &	perfGraph ()

const Parallel::Communicator &	comm () const

processor_id_type	n_processors () const

processor_id_type	processor_id () const

Static Public Member Functions
static InputParameters	validParams ()

Public Attributes
const ConsoleStream	_console

Static Public Attributes
static constexpr auto	SYSTEM

static constexpr auto	NAME

Protected Member Functions
virtual void	addObjects () override
	Adds all non-mesh objects. More...

virtual void	addMesh () override
	Adds the mesh. More...

virtual void	addInitialConditions () override
	Adds the initial conditions. More...

virtual void	addSolutionVariables () override
	Adds the solution variables. More...

virtual void	addAuxVariables () override
	Adds aux variables. More...

virtual void	addMaterials () override
	Adds materials. More...

virtual void	addUserObjects () override
	Adds user objects. More...

void	addSolutionVariablesRiemannIC (const std::vector< VariableName > &variables, const std::vector< Real > &values_left, const std::vector< Real > &values_right)
	Adds solution variables with Riemann problem IC (constant left and right states) More...

virtual void	addPreconditioner () override
	Adds the preconditioner. More...

void	addMeshInternal (const unsigned int &nx)
	Adds the mesh with a provided number of elements in x direction. More...

void	addScalarVariables (const std::vector< VariableName > &names, const std::vector< FunctionName > &values)
	Adds scalar variables. More...

void	addSolutionVariable (const VariableName &var_name, const std::string &family="LAGRANGE", const std::string &order="FIRST", const Real &scaling=1.0)
	Adds a solution variable. More...

void	addAuxVariable (const VariableName &var_name, const std::string &fe_family, const std::string &fe_order)
	Adds an aux variable. More...

void	addConstantIC (const VariableName &var_name, const Real &value)
	Adds a constant initial condition. More...

void	addFunctionIC (const VariableName &var_name, const FunctionName &function_name)
	Adds a function initial condition. More...

virtual void	addExecutioner ()
	Adds the executioner. More...

virtual void	addOutput ()
	Add output. More...

bool	addRelationshipManagers (Moose::RelationshipManagerType when_type, const InputParameters &moose_object_pars)

void	associateWithParameter (const std::string &param_name, InputParameters &params) const

void	associateWithParameter (const InputParameters &from_params, const std::string &param_name, InputParameters &params) 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

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

Protected Attributes
const VariableName	_A_name
	cross-sectional area variable name, elemental average More...

const VariableName	_A_linear_name
	cross-sectional area variable name, linear Lagrange More...

const VariableName	_rhoA_name
	rho*A variable name More...

const VariableName	_rhouA_name
	rhouA variable name More...

const VariableName	_rhoEA_name
	rhoEA variable name More...

const bool	_add_dg_kernel
	option to add DG kernel More...

const bool	_add_bc
	option to BC More...

const UserObjectName	_numerical_flux_name
	name of numerical flux user object being tested More...

const UserObjectName	_boundary_flux_name
	name of boundary flux user object being tested More...

const MooseEnum	_ic_option
	initial conditions option More...

const FunctionName &	_A_fn_name
	area function name More...

const bool	_use_slope_reconstruction
	option to use slope reconstruction More...

const std::string	_reconstruction_material_name
	reconstruction material name More...

const MaterialPropertyName	_direction_name
	direction material property name More...

const UserObjectName	_fp_name
	fluid properties object name More...

const std::string	_snes_test_err
	Finite differencing parameter. More...

bool	_default_use_transient_executioner
	Default for option to use a transient executioner. More...

const std::vector< VariableName >	_scalar_variables
	List of scalar variables to add. More...

const std::vector< FunctionName >	_scalar_variable_values
	List of values for the scalar variables to add. More...

const std::vector< VariableName >	_aux_variables
	List of aux variables to add. More...

const std::vector< FunctionName >	_aux_variable_values
	List of function names for aux variables to add. More...

const std::vector< std::string >	_mat_property_names
	List of material properties to add. More...

const std::vector< FunctionName >	_mat_property_values
	List of function names for material properties to add. More...

const std::string	_fe_family
	Default FE family. More...

const std::string	_fe_order
	Default FE order. More...

const bool &	_ad
	True for setting up testing with AD, false otherwise. More...

std::string	_registered_identifier

std::string	_specific_task_name

std::set< std::string >	_all_tasks

ActionWarehouse &	_awh

const std::string &	_current_task

std::shared_ptr< MooseMesh > &	_mesh

std::shared_ptr< MooseMesh > &	_displaced_mesh

std::shared_ptr< FEProblemBase > &	_problem

PerfID	_act_timer

MooseApp &	_app

const std::string	_type

const std::string	_name

const InputParameters &	_pars

Factory &	_factory

ActionFactory &	_action_factory

MooseApp &	_pg_moose_app

const std::string	_prefix

const Parallel::Communicator &	_communicator

Detailed Description

Sets up a Jacobian test for 1-phase rDG.

Definition at line 18 of file JacobianTest1PhaseRDGAction.h.

Constructor & Destructor Documentation

◆ JacobianTest1PhaseRDGAction()

JacobianTest1PhaseRDGAction::JacobianTest1PhaseRDGAction ( const InputParameters & params )

Definition at line 50 of file JacobianTest1PhaseRDGAction.C.

   : JacobianTestAction(params),
     _A_name("A"),
     _A_linear_name("A_linear"),
     _rhoA_name("rhoA"),
     _rhouA_name("rhouA"),
     _rhoEA_name("rhoEA"),
     _add_dg_kernel(getParam<bool>("add_dg_kernel")),
     _add_bc(getParam<bool>("add_bc")),
     _numerical_flux_name(getParam<UserObjectName>("numerical_flux")),
     _boundary_flux_name(getParam<UserObjectName>("boundary_flux")),
     _ic_option(getParam<MooseEnum>("ic_option")),
     _A_fn_name(getParam<FunctionName>("A_function")),
     _use_slope_reconstruction(getParam<bool>("use_slope_reconstruction")),
     _reconstruction_material_name("reconstruction_material"),
     _direction_name("direction"),
     _fp_name("fluid_properties")
 {
   if (_add_dg_kernel && !_pars.isParamSetByUser("numerical_flux"))
     mooseError("The parameter 'numerical_flux' must be provided when adding a DG kernel.");
   if (_add_bc && !_pars.isParamSetByUser("boundary_flux"))
     mooseError("The parameter 'boundary_flux' must be provided when adding a BC.");
 }

Member Function Documentation

◆ act()

void TestAction::act ( )

virtualinherited

Implements Action.

Definition at line 75 of file TestAction.C.

 {
   if (_current_task == "meta_action")
   {
     if (getParam<bool>("generate_mesh"))
       addMesh();
 
     addObjects();
 
     {
       const std::string class_name = "CreateProblemAction";
       auto action_params = _action_factory.getValidParams(class_name);
       action_params.set<std::string>("type") = "FEProblem";
 
       std::shared_ptr<MooseObjectAction> action = std::static_pointer_cast<MooseObjectAction>(
           _action_factory.create(class_name, "fe_problem", action_params));
 
       action->getObjectParams().set<MooseEnum>("kernel_coverage_check") = "false";
       _awh.addActionBlock(action);
     }
   }
 }

◆ addAuxVariable()

void TestAction::addAuxVariable	(	const VariableName &	var_name,
		const std::string &	fe_family,
		const std::string &	fe_order
	)

protectedinherited

Adds an aux variable.

Parameters

[in]	var_name	name of the variable to add
[in]	fe_family	finite element family
[in]	fe_order	finite element order

Definition at line 248 of file TestAction.C.

Referenced by ClosureTestAction::addAuxVariables(), addAuxVariables(), TestAction::addAuxVariables(), ClosureTestAction::setupADOutput(), and ClosureTestAction::setupOutput().

 {
   const std::string class_name = "AddAuxVariableAction";
   InputParameters params = _action_factory.getValidParams(class_name);
 
   params.set<MooseEnum>("family") = fe_family;
   params.set<MooseEnum>("order") = fe_order;
 
   std::shared_ptr<Action> action =
       std::static_pointer_cast<Action>(_action_factory.create(class_name, var_name, params));
 
   _awh.addActionBlock(action);
 }

◆ addAuxVariables()

void JacobianTest1PhaseRDGAction::addAuxVariables ( )

overrideprotectedvirtual

Adds aux variables.

Reimplemented from TestAction.

Definition at line 239 of file JacobianTest1PhaseRDGAction.C.

 {
   addAuxVariable(_A_name, _fe_family, _fe_order);
   addFunctionIC(_A_name, _A_fn_name);
   if (getParam<bool>("use_elem_area"))
     addAuxVariable(_A_linear_name, _fe_family, _fe_order);
   else
     addAuxVariable(_A_linear_name, "LAGRANGE", "FIRST");
   addFunctionIC(_A_linear_name, _A_fn_name);
 }

◆ addConstantIC()

void TestAction::addConstantIC	(	const VariableName &	var_name,
		const Real &	value
	)

protectedinherited

Adds a constant initial condition.

Parameters

[in]	var_name	name of the variable for which to add initial condition
[in]	value	value of the initial condition

Definition at line 265 of file TestAction.C.

Referenced by addSolutionVariables().

 {
   const std::string class_name = "AddInitialConditionAction";
   InputParameters params = _action_factory.getValidParams(class_name);
   params.set<std::string>("type") = "ConstantIC";
 
   std::shared_ptr<MooseObjectAction> action = std::static_pointer_cast<MooseObjectAction>(
       _action_factory.create(class_name, var_name + "_IC", params));
 
   action->getObjectParams().set<VariableName>("variable") = var_name;
   action->getObjectParams().set<Real>("value") = value;
 
   _awh.addActionBlock(action);
 }

◆ addExecutioner()

void TestAction::addExecutioner ( )

protectedvirtualinherited

Adds the executioner.

Definition at line 302 of file TestAction.C.

Referenced by TestAction::addObjects().

 {
   const std::string class_name = "CreateExecutionerAction";
 
   // determine whether to use a transient executioner - different derived
   // classes have different defaults: the declared default is irrelevant
   bool use_transient_executioner;
   if (isParamValid("use_transient_executioner"))
     use_transient_executioner = getParam<bool>("use_transient_executioner");
   else
     use_transient_executioner = _default_use_transient_executioner;
 
   // Due to more consistent divergence status reporting in PETSc (as of 5f3c5e7a), users should have
   // the option to abort on the first fail if desired. Otherwise Jacobian testing, for example, could
   // fail in undesired ways, even if the Jacobian test achieves a passing result.
   bool abort_on_solve_fail = getParam<bool>("abort_on_solve_fail");
 
   // if a time kernel is being tested, then use a transient executioner instead of steady
   if (use_transient_executioner)
   {
     InputParameters params = _action_factory.getValidParams(class_name);
     params.set<std::string>("type") = "Transient";
 
     std::shared_ptr<MooseObjectAction> action = std::static_pointer_cast<MooseObjectAction>(
         _action_factory.create(class_name, "executioner", params));
 
     action->getObjectParams().set<unsigned int>("num_steps") = 1;
 
     if (abort_on_solve_fail)
       action->getObjectParams().set<bool>("abort_on_solve_fail") = abort_on_solve_fail;
 
     _awh.addActionBlock(action);
   }
   else
   {
     InputParameters params = _action_factory.getValidParams(class_name);
     params.set<std::string>("type") = "Steady";
 
     std::shared_ptr<MooseObjectAction> action = std::static_pointer_cast<MooseObjectAction>(
         _action_factory.create(class_name, "executioner", params));
 
     if (abort_on_solve_fail)
       action->getObjectParams().set<bool>("abort_on_solve_fail") = abort_on_solve_fail;
 
     _awh.addActionBlock(action);
   }
 }

◆ addFunctionIC()

void TestAction::addFunctionIC	(	const VariableName &	var_name,
		const FunctionName &	function_name
	)

protectedinherited

Adds a function initial condition.

Parameters

[in]	var_name	name of the variable for which to add initial condition
[in]	function_name	names of the IC function

Definition at line 281 of file TestAction.C.

Referenced by addAuxVariables(), TestAction::addAuxVariables(), ClosureTestAction::addInitialConditions(), JacobianTestGeneralAction::addSolutionVariables(), and addSolutionVariablesRiemannIC().

 {
   const std::string class_name = "AddInitialConditionAction";
   InputParameters params = _action_factory.getValidParams(class_name);
   params.set<std::string>("type") = "FunctionIC";
 
   std::shared_ptr<MooseObjectAction> action = std::static_pointer_cast<MooseObjectAction>(
       _action_factory.create(class_name, var_name + "_IC", params));
 
   action->getObjectParams().set<VariableName>("variable") = var_name;
   action->getObjectParams().set<FunctionName>("function") = function_name;
 
   _awh.addActionBlock(action);
 }

◆ addInitialConditions()

void JacobianTest1PhaseRDGAction::addInitialConditions ( )

overrideprotectedvirtual

Adds the initial conditions.

Implements TestAction.

Definition at line 152 of file JacobianTest1PhaseRDGAction.C.

153 {

154 }

◆ addMaterials()

void JacobianTest1PhaseRDGAction::addMaterials ( )

overrideprotectedvirtual

Adds materials.

Reimplemented from TestAction.

Definition at line 251 of file JacobianTest1PhaseRDGAction.C.

 {
   // direction
   {
     const std::string class_name = "AddMaterialAction";
     InputParameters params = _action_factory.getValidParams(class_name);
     params.set<std::string>("type") = "DirectionMaterial";
 
     std::shared_ptr<MooseObjectAction> action = std::static_pointer_cast<MooseObjectAction>(
         _action_factory.create(class_name, "direction_material", params));
 
     _awh.addActionBlock(action);
   }
 
   // fluid properties
   {
     const std::string class_name = "AddMaterialAction";
     InputParameters params = _action_factory.getValidParams(class_name);
     if (_ad)
       params.set<std::string>("type") = "ADFluidProperties3EqnMaterial";
     else
       params.set<std::string>("type") = "FluidProperties3EqnMaterial";
 
     std::shared_ptr<MooseObjectAction> action = std::static_pointer_cast<MooseObjectAction>(
         _action_factory.create(class_name, "fluid_properties_material", params));
 
     action->getObjectParams().set<std::vector<VariableName>>("A") = {_A_name};
     action->getObjectParams().set<std::vector<VariableName>>("rhoA") = {_rhoA_name};
     action->getObjectParams().set<std::vector<VariableName>>("rhouA") = {_rhouA_name};
     action->getObjectParams().set<std::vector<VariableName>>("rhoEA") = {_rhoEA_name};
 
     action->getObjectParams().set<UserObjectName>("fp") = _fp_name;
 
     _awh.addActionBlock(action);
   }
 
   // reconstruction material
   {
     const std::string class_name = "AddMaterialAction";
     InputParameters params = _action_factory.getValidParams(class_name);
     params.set<std::string>("type") = "RDG3EqnMaterial";
 
     std::shared_ptr<MooseObjectAction> action = std::static_pointer_cast<MooseObjectAction>(
         _action_factory.create(class_name, _reconstruction_material_name, params));
 
     if (_use_slope_reconstruction)
       action->getObjectParams().set<MooseEnum>("scheme") = "Minmod";
     else
       action->getObjectParams().set<MooseEnum>("scheme") = "None";
 
     action->getObjectParams().set<std::vector<VariableName>>("A_elem") = {_A_name};
     action->getObjectParams().set<std::vector<VariableName>>("A_linear") = {_A_linear_name};
     action->getObjectParams().set<std::vector<VariableName>>("rhoA") = {_rhoA_name};
     action->getObjectParams().set<std::vector<VariableName>>("rhouA") = {_rhouA_name};
     action->getObjectParams().set<std::vector<VariableName>>("rhoEA") = {_rhoEA_name};
     action->getObjectParams().set<MaterialPropertyName>("direction") = _direction_name;
 
     action->getObjectParams().set<UserObjectName>("fluid_properties") = _fp_name;
 
     _awh.addActionBlock(action);
   }
 }

◆ addMesh()

void JacobianTest1PhaseRDGAction::addMesh ( )

overrideprotectedvirtual

Adds the mesh.

Reimplemented from TestAction.

Definition at line 142 of file JacobianTest1PhaseRDGAction.C.

 {
   // for internal side flux testing, use 2 elements (need an internal side)
   if (_add_dg_kernel)
     addMeshInternal(2);
   else
     addMeshInternal(1);
 }

◆ addMeshInternal()

void TestAction::addMeshInternal ( const unsigned int & nx )

protectedinherited

Adds the mesh with a provided number of elements in x direction.

Parameters

[in] nx number of elements in x direction

Definition at line 105 of file TestAction.C.

Referenced by addMesh(), and TestAction::addMesh().

 {
   std::vector<std::string> setup_mesh_tasks{"setup_mesh", "set_mesh_base", "init_mesh"};
   for (const std::string & task : setup_mesh_tasks)
   {
     const std::string class_name = "SetupMeshAction";
     InputParameters params = _action_factory.getValidParams(class_name);
     params.set<std::string>("type") = "GeneratedMesh";
     params.set<std::string>("task") = task;
 
     std::shared_ptr<MooseObjectAction> action = std::static_pointer_cast<MooseObjectAction>(
         _action_factory.create(class_name, task, params));
 
     action->getObjectParams().set<MooseEnum>("dim") = 1;
     action->getObjectParams().set<unsigned int>("nx") = nx;
     action->getObjectParams().set<bool>("allow_renumbering") = false;
 
     _awh.addActionBlock(action);
   }
 
   std::vector<std::string> setup_mesh_complete_tasks{"prepare_mesh", "setup_mesh_complete"};
   for (const std::string & task : setup_mesh_complete_tasks)
   {
     const std::string class_name = "SetupMeshCompleteAction";
     InputParameters params = _action_factory.getValidParams(class_name);
     params.set<std::string>("task") = task;
 
     std::shared_ptr<Action> action =
         std::static_pointer_cast<Action>(_action_factory.create(class_name, task, params));
 
     _awh.addActionBlock(action);
   }
 }

◆ addObjects()

void JacobianTest1PhaseRDGAction::addObjects ( )

overrideprotectedvirtual

Adds all non-mesh objects.

Reimplemented from TestAction.

Definition at line 75 of file JacobianTest1PhaseRDGAction.C.

 {
   TestAction::addObjects();
 
   const std::vector<VariableName> variables = {_rhoA_name, _rhouA_name, _rhoEA_name};
 
   if (_add_dg_kernel)
   {
     // add the DG kernel to use the internal flux
     for (unsigned int i = 0; i < variables.size(); i++)
     {
       const std::string class_name = "AddDGKernelAction";
       InputParameters params = _action_factory.getValidParams(class_name);
       params.set<std::string>("type") = "NumericalFlux3EqnDGKernel";
       params.set<NonlinearVariableName>("variable") = variables[i];
       params.set<std::vector<SubdomainName>>("block") = {"0"};
 
       std::shared_ptr<MooseObjectAction> action = std::static_pointer_cast<MooseObjectAction>(
           _action_factory.create(class_name, "internal_side_flux_" + variables[i], params));
 
       action->getObjectParams().set<NonlinearVariableName>("variable") = variables[i];
       action->getObjectParams().set<std::vector<SubdomainName>>("block") = {"0"};
 
       action->getObjectParams().set<std::vector<VariableName>>("A_linear") = {_A_linear_name};
       action->getObjectParams().set<std::vector<VariableName>>("rhoA") = {_rhoA_name};
       action->getObjectParams().set<std::vector<VariableName>>("rhouA") = {_rhouA_name};
       action->getObjectParams().set<std::vector<VariableName>>("rhoEA") = {_rhoEA_name};
 
       action->getObjectParams().set<UserObjectName>("numerical_flux") = _numerical_flux_name;
 
       _awh.addActionBlock(action);
     }
   }
 
   if (_add_bc)
   {
     // add the BC to use the boundary flux
     for (unsigned int i = 0; i < variables.size(); i++)
     {
       const std::string class_name = "AddBCAction";
       InputParameters params = _action_factory.getValidParams(class_name);
       params.set<std::string>("type") = "BoundaryFlux3EqnBC";
       params.set<NonlinearVariableName>("variable") = variables[i];
       params.set<Real>("normal") = 1.;
       params.set<std::vector<BoundaryName>>("boundary") = {"0"};
 
       std::shared_ptr<MooseObjectAction> action = std::static_pointer_cast<MooseObjectAction>(
           _action_factory.create(class_name, "bc_" + variables[i], params));
 
       action->getObjectParams().set<NonlinearVariableName>("variable") = variables[i];
       action->getObjectParams().set<Real>("normal") = 1.;
       action->getObjectParams().set<std::vector<BoundaryName>>("boundary") = {"0"};
 
       action->getObjectParams().set<std::vector<VariableName>>("A_elem") = {_A_name};
       action->getObjectParams().set<std::vector<VariableName>>("A_linear") = {_A_linear_name};
       action->getObjectParams().set<std::vector<VariableName>>("rhoA") = {_rhoA_name};
       action->getObjectParams().set<std::vector<VariableName>>("rhouA") = {_rhouA_name};
       action->getObjectParams().set<std::vector<VariableName>>("rhoEA") = {_rhoEA_name};
 
       action->getObjectParams().set<UserObjectName>("boundary_flux") = _boundary_flux_name;
 
       _awh.addActionBlock(action);
     }
   }
 }

◆ addOutput()

void TestAction::addOutput ( )

protectedvirtualinherited

Add output.

Reimplemented in ClosureTestAction.

Definition at line 351 of file TestAction.C.

Referenced by TestAction::addObjects().

352 {

353 }

◆ addPreconditioner()

void JacobianTestAction::addPreconditioner ( )

overrideprotectedvirtualinherited

Adds the preconditioner.

Reimplemented from TestAction.

Definition at line 34 of file JacobianTestAction.C.

 {
   const std::string class_name = "SetupPreconditionerAction";
   InputParameters params = _action_factory.getValidParams(class_name);
   params.set<std::string>("type") = "SMP";
 
   std::shared_ptr<MooseObjectAction> action = std::static_pointer_cast<MooseObjectAction>(
       _action_factory.create(class_name, "smp", params));
 
   action->getObjectParams().set<bool>("full") = true;
   action->getObjectParams().set<MooseEnum>("solve_type") = "newton";
 #if PETSC_VERSION_LESS_THAN(3, 9, 0)
   action->getObjectParams().set<MultiMooseEnum>("petsc_options_iname") =
       "-snes_type -snes_test_err";
   action->getObjectParams().set<std::vector<std::string>>("petsc_options_value") = {"test",
                                                                                     _snes_test_err};
 #else
   action->getObjectParams().set<MultiMooseEnum>("petsc_options") = "-snes_test_jacobian";
   action->getObjectParams().set<MultiMooseEnum>("petsc_options_iname") = "-snes_test_err";
   action->getObjectParams().set<std::vector<std::string>>("petsc_options_value") = {_snes_test_err};
 #endif
 
   _awh.addActionBlock(action);
 }

◆ addScalarVariables()

void TestAction::addScalarVariables	(	const std::vector< VariableName > &	names,
		const std::vector< FunctionName > &	values
	)

protectedinherited

Adds scalar variables.

Parameters

[in]	names	names of the variables to add
[in]	values	values of the variables to add

Definition at line 154 of file TestAction.C.

Referenced by TestAction::addObjects().

 {
   for (unsigned int i = 0; i < names.size(); ++i)
   {
     // add the scalar variable
     addSolutionVariable(names[i], "SCALAR", "FIRST");
 
     // add its IC
     {
       const std::string class_name = "AddInitialConditionAction";
       InputParameters params = _action_factory.getValidParams(class_name);
       params.set<std::string>("type") = "FunctionScalarIC";
 
       std::shared_ptr<MooseObjectAction> action = std::static_pointer_cast<MooseObjectAction>(
           _action_factory.create(class_name, names[i] + "_IC", params));
 
       action->getObjectParams().set<VariableName>("variable") = names[i];
       action->getObjectParams().set<std::vector<FunctionName>>("function") = {values[i]};
 
       _awh.addActionBlock(action);
     }
   }
 }

◆ addSolutionVariable()

void TestAction::addSolutionVariable	(	const VariableName &	var_name,
		const std::string &	family = `"LAGRANGE"`,
		const std::string &	order = `"FIRST"`,
		const Real &	scaling = `1.0`
	)

protectedinherited

Adds a solution variable.

Parameters

[in]	var_name	name of the variable to add
[in]	family	variable family
[in]	order	variable order
[in]	scaling	scaling factor to apply to variable

Definition at line 230 of file TestAction.C.

Referenced by TestAction::addScalarVariables(), ClosureTestAction::addSolutionVariables(), addSolutionVariables(), JacobianTestGeneralAction::addSolutionVariables(), and addSolutionVariablesRiemannIC().

 {
   const std::string class_name = "AddVariableAction";
   InputParameters params = _action_factory.getValidParams(class_name);
   params.set<MooseEnum>("family") = family;
   params.set<MooseEnum>("order") = order;
   params.set<std::vector<Real>>("scaling") = {scaling};
 
   std::shared_ptr<Action> action =
       std::static_pointer_cast<Action>(_action_factory.create(class_name, var_name, params));
 
   _awh.addActionBlock(action);
 }

◆ addSolutionVariables()

void JacobianTest1PhaseRDGAction::addSolutionVariables ( )

overrideprotectedvirtual

Adds the solution variables.

Implements TestAction.

Definition at line 157 of file JacobianTest1PhaseRDGAction.C.

 {
   const std::vector<VariableName> variables = {_rhoA_name, _rhouA_name, _rhoEA_name};
 
   if (_ic_option == "constant")
   {
     const std::vector<Real> values = {2.0, 4.0, 30.0};
 
     for (std::size_t i = 0; i < variables.size(); ++i)
     {
       addSolutionVariable(variables[i]);
       addConstantIC(variables[i], values[i]);
     }
   }
   else if (_ic_option == "riemann_L") // left region
   {
     const std::vector<Real> values_left = {2.2, 4.4, 11.2};
     const std::vector<Real> values_right = {2.6, 6.4, 12.2};
 
     addSolutionVariablesRiemannIC(variables, values_left, values_right);
   }
   else if (_ic_option == "riemann_LM") // left star region
   {
     const std::vector<Real> values_left = {2.2, -4.4, 11.2};
     const std::vector<Real> values_right = {3.2, 8.4, 12.2};
 
     addSolutionVariablesRiemannIC(variables, values_left, values_right);
   }
   else if (_ic_option == "riemann_RM") // right star region
   {
     const std::vector<Real> values_left = {2.2, 4.4, 11.2};
     const std::vector<Real> values_right = {1.8, -5.0, 12.2};
 
     addSolutionVariablesRiemannIC(variables, values_left, values_right);
   }
   else if (_ic_option == "riemann_R") // right region
   {
     const std::vector<Real> values_left = {2.2, 4.4, 11.2};
     const std::vector<Real> values_right = {2.4, -16.0, 53.4};
 
     addSolutionVariablesRiemannIC(variables, values_left, values_right);
   }
   else
   {
     mooseError("Invalid IC option");
   }
 }

◆ addSolutionVariablesRiemannIC()

void JacobianTest1PhaseRDGAction::addSolutionVariablesRiemannIC	(	const std::vector< VariableName > &	variables,
		const std::vector< Real > &	values_left,
		const std::vector< Real > &	values_right
	)

protected

Adds solution variables with Riemann problem IC (constant left and right states)

Parameters

[in]	variables	names of the solution variables
[in]	values_left	values of the solution variables on the left half
[in]	values_right	values of the solution variables on the right half

Definition at line 206 of file JacobianTest1PhaseRDGAction.C.

Referenced by addSolutionVariables().

 {
   for (unsigned int i = 0; i < variables.size(); ++i)
   {
     // add the variable
     addSolutionVariable(variables[i], "MONOMIAL", "CONSTANT");
 
     // add the IC function
     FunctionName ic_fn_name;
     {
       const std::string class_name = "AddFunctionAction";
       InputParameters params = _action_factory.getValidParams(class_name);
       params.set<std::string>("type") = "PiecewiseConstant";
 
       ic_fn_name = variables[i] + "_IC_fn";
       std::shared_ptr<MooseObjectAction> action = std::static_pointer_cast<MooseObjectAction>(
           _action_factory.create(class_name, ic_fn_name, params));
 
       action->getObjectParams().set<MooseEnum>("axis") = "x";
       action->getObjectParams().set<std::vector<Real>>("x") = {0.0, 0.5};
       action->getObjectParams().set<std::vector<Real>>("y") = {values_left[i], values_right[i]};
 
       _awh.addActionBlock(action);
     }
 
     addFunctionIC(variables[i], ic_fn_name);
   }
 }

◆ addUserObjects()

void JacobianTest1PhaseRDGAction::addUserObjects ( )

overrideprotectedvirtual

Adds user objects.

Implements TestAction.

Definition at line 315 of file JacobianTest1PhaseRDGAction.C.

 {
   // fluid properties
   {
     const std::string class_name = "AddFluidPropertiesAction";
     InputParameters params = _action_factory.getValidParams(class_name);
     params.set<std::string>("type") = "IdealGasFluidProperties";
 
     std::shared_ptr<MooseObjectAction> action = std::static_pointer_cast<MooseObjectAction>(
         _action_factory.create(class_name, _fp_name, params));
 
     action->getObjectParams().set<Real>("gamma") = 1.5;
     action->getObjectParams().set<Real>("molar_mass") = 0.83144598;
 
     _awh.addActionBlock(action);
   }
 }

◆ validParams()

InputParameters JacobianTest1PhaseRDGAction::validParams ( )

static

Definition at line 19 of file JacobianTest1PhaseRDGAction.C.

 {
   InputParameters params = emptyInputParameters();
   params += JacobianTestAction::validParams();
 
   params.addClassDescription("Sets up a Jacobian test for rDG");
 
   params.addParam<bool>("add_dg_kernel", false, "Option to add DG kernel");
   params.addParam<bool>("add_bc", false, "Option to add BC");
 
   params.addParam<UserObjectName>(
       "numerical_flux", "", "Name of the numerical flux user object to test");
   params.addParam<UserObjectName>(
       "boundary_flux", "", "Name of the boundary flux user object to test");
 
   MooseEnum ic_option("constant riemann_L riemann_LM riemann_RM riemann_R");
   params.addRequiredParam<MooseEnum>("ic_option", ic_option, "IC option");
 
   params.addParam<FunctionName>("A_function", "2.0", "Area function");
 
   params.addParam<bool>("use_slope_reconstruction", true, "Use slope reconstruction?");
 
   params.addParam<bool>(
       "use_elem_area", false, "Use the elemental area variable instead of linear area");
 
   params.set<std::string>("fe_family") = "MONOMIAL";
   params.set<std::string>("fe_order") = "CONSTANT";
 
   return params;
 }

Member Data Documentation

◆ _A_fn_name

const FunctionName& JacobianTest1PhaseRDGAction::_A_fn_name

protected

area function name

Definition at line 67 of file JacobianTest1PhaseRDGAction.h.

Referenced by addAuxVariables().

◆ _A_linear_name

const VariableName JacobianTest1PhaseRDGAction::_A_linear_name

protected

cross-sectional area variable name, linear Lagrange

Definition at line 46 of file JacobianTest1PhaseRDGAction.h.

Referenced by addAuxVariables(), addMaterials(), and addObjects().

◆ _A_name

const VariableName JacobianTest1PhaseRDGAction::_A_name

protected

cross-sectional area variable name, elemental average

Definition at line 44 of file JacobianTest1PhaseRDGAction.h.

Referenced by addAuxVariables(), addMaterials(), and addObjects().

◆ _ad

const bool& TestAction::_ad

protectedinherited

True for setting up testing with AD, false otherwise.

Definition at line 154 of file TestAction.h.

Referenced by ClosureTestAction::addMaterials(), and addMaterials().

◆ _add_bc

const bool JacobianTest1PhaseRDGAction::_add_bc

protected

option to BC

Definition at line 57 of file JacobianTest1PhaseRDGAction.h.

Referenced by addObjects(), and JacobianTest1PhaseRDGAction().

◆ _add_dg_kernel

const bool JacobianTest1PhaseRDGAction::_add_dg_kernel

protected

option to add DG kernel

Definition at line 55 of file JacobianTest1PhaseRDGAction.h.

Referenced by addMesh(), addObjects(), and JacobianTest1PhaseRDGAction().

◆ _aux_variable_values

const std::vector<FunctionName> TestAction::_aux_variable_values

protectedinherited

List of function names for aux variables to add.

Definition at line 142 of file TestAction.h.

Referenced by TestAction::addAuxVariables(), and TestAction::TestAction().

◆ _aux_variables

const std::vector<VariableName> TestAction::_aux_variables

protectedinherited

List of aux variables to add.

Definition at line 140 of file TestAction.h.

Referenced by TestAction::addAuxVariables(), and TestAction::TestAction().

◆ _boundary_flux_name

const UserObjectName JacobianTest1PhaseRDGAction::_boundary_flux_name

protected

name of boundary flux user object being tested

Definition at line 62 of file JacobianTest1PhaseRDGAction.h.

Referenced by addObjects().

◆ _default_use_transient_executioner

bool TestAction::_default_use_transient_executioner

protectedinherited

Default for option to use a transient executioner.

Definition at line 132 of file TestAction.h.

Referenced by TestAction::addExecutioner(), and ClosureTestAction::ClosureTestAction().

◆ _direction_name

const MaterialPropertyName JacobianTest1PhaseRDGAction::_direction_name

protected

direction material property name

Definition at line 76 of file JacobianTest1PhaseRDGAction.h.

Referenced by addMaterials().

◆ _fe_family

const std::string TestAction::_fe_family

protectedinherited

Default FE family.

Definition at line 150 of file TestAction.h.

Referenced by ClosureTestAction::addAuxVariables(), addAuxVariables(), and TestAction::addAuxVariables().

◆ _fe_order

const std::string TestAction::_fe_order

protectedinherited

Default FE order.

Definition at line 152 of file TestAction.h.

Referenced by ClosureTestAction::addAuxVariables(), addAuxVariables(), and TestAction::addAuxVariables().

◆ _fp_name

const UserObjectName JacobianTest1PhaseRDGAction::_fp_name

protected

fluid properties object name

Definition at line 79 of file JacobianTest1PhaseRDGAction.h.

Referenced by addMaterials(), and addUserObjects().

◆ _ic_option

const MooseEnum JacobianTest1PhaseRDGAction::_ic_option

protected

initial conditions option

Definition at line 65 of file JacobianTest1PhaseRDGAction.h.

Referenced by addSolutionVariables().

◆ _mat_property_names

const std::vector<std::string> TestAction::_mat_property_names

protectedinherited

List of material properties to add.

Definition at line 145 of file TestAction.h.

Referenced by TestAction::addMaterials(), and TestAction::TestAction().

◆ _mat_property_values

const std::vector<FunctionName> TestAction::_mat_property_values

protectedinherited

List of function names for material properties to add.

Definition at line 147 of file TestAction.h.

Referenced by TestAction::addMaterials(), and TestAction::TestAction().

◆ _numerical_flux_name

const UserObjectName JacobianTest1PhaseRDGAction::_numerical_flux_name

protected

name of numerical flux user object being tested

Definition at line 60 of file JacobianTest1PhaseRDGAction.h.

Referenced by addObjects().

◆ _reconstruction_material_name

const std::string JacobianTest1PhaseRDGAction::_reconstruction_material_name

protected

reconstruction material name

Definition at line 73 of file JacobianTest1PhaseRDGAction.h.

Referenced by addMaterials().

◆ _rhoA_name

const VariableName JacobianTest1PhaseRDGAction::_rhoA_name

protected

rho*A variable name

Definition at line 48 of file JacobianTest1PhaseRDGAction.h.

Referenced by addMaterials(), addObjects(), and addSolutionVariables().

◆ _rhoEA_name

const VariableName JacobianTest1PhaseRDGAction::_rhoEA_name

protected

rho*E*A variable name

Definition at line 52 of file JacobianTest1PhaseRDGAction.h.

Referenced by addMaterials(), addObjects(), and addSolutionVariables().

◆ _rhouA_name

const VariableName JacobianTest1PhaseRDGAction::_rhouA_name

protected

rho*u*A variable name

Definition at line 50 of file JacobianTest1PhaseRDGAction.h.

Referenced by addMaterials(), addObjects(), and addSolutionVariables().

◆ _scalar_variable_values

const std::vector<FunctionName> TestAction::_scalar_variable_values

protectedinherited

List of values for the scalar variables to add.

Definition at line 137 of file TestAction.h.

Referenced by TestAction::addObjects(), and TestAction::TestAction().

◆ _scalar_variables

const std::vector<VariableName> TestAction::_scalar_variables

protectedinherited

List of scalar variables to add.

Definition at line 135 of file TestAction.h.

Referenced by TestAction::addObjects(), and TestAction::TestAction().

◆ _snes_test_err

const std::string JacobianTestAction::_snes_test_err

protectedinherited

Finite differencing parameter.

Definition at line 26 of file JacobianTestAction.h.

Referenced by JacobianTestAction::addPreconditioner().

◆ _use_slope_reconstruction

const bool JacobianTest1PhaseRDGAction::_use_slope_reconstruction

protected

option to use slope reconstruction

Definition at line 70 of file JacobianTest1PhaseRDGAction.h.

Referenced by addMaterials().

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

thermal_hydraulics/test/include/actions/JacobianTest1PhaseRDGAction.h
thermal_hydraulics/test/src/actions/JacobianTest1PhaseRDGAction.C

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ JacobianTest1PhaseRDGAction()

Member Function Documentation

◆ act()

◆ addAuxVariable()

◆ addAuxVariables()

◆ addConstantIC()

◆ addExecutioner()

◆ addFunctionIC()

◆ addInitialConditions()

◆ addMaterials()

◆ addMesh()

◆ addMeshInternal()

◆ addObjects()

◆ addOutput()

◆ addPreconditioner()

◆ addScalarVariables()

◆ addSolutionVariable()

◆ addSolutionVariables()

◆ addSolutionVariablesRiemannIC()

◆ addUserObjects()

◆ validParams()

Member Data Documentation

◆ _A_fn_name

◆ _A_linear_name

◆ _A_name

◆ _ad

◆ _add_bc

◆ _add_dg_kernel

◆ _aux_variable_values

◆ _aux_variables

◆ _boundary_flux_name

◆ _default_use_transient_executioner

◆ _direction_name

◆ _fe_family

◆ _fe_order

◆ _fp_name

◆ _ic_option

◆ _mat_property_names

◆ _mat_property_values

◆ _numerical_flux_name

◆ _reconstruction_material_name

◆ _rhoA_name

◆ _rhoEA_name

◆ _rhouA_name

◆ _scalar_variable_values

◆ _scalar_variables

◆ _snes_test_err

◆ _use_slope_reconstruction