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

const std::string &	name () const

std::string	typeAndName () const

MooseObjectParameterName	uniqueParameterName (const std::string &parameter_name) const

MooseObjectName	uniqueName () const

const InputParameters &	parameters () const

const hit::Node *	getHitNode () const

bool	hasBase () const

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

void	connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const

void	paramError (const std::string &param, Args... args) const

void	paramWarning (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

std::string	messagePrefix (const bool hit_prefix=true) const

std::string	errorPrefix (const std::string &) const

void	mooseError (Args &&... args) const

void	mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const

void	mooseErrorNonPrefixed (Args &&... args) const

void	mooseWarning (Args &&... args) const

void	mooseWarning (Args &&... args) const

void	mooseWarningNonPrefixed (Args &&... args) const

void	mooseWarningNonPrefixed (Args &&... args) const

void	mooseDeprecated (Args &&... args) const

void	mooseDeprecated (Args &&... args) const

void	mooseDeprecatedNoTrace (Args &&... args) const

void	mooseInfo (Args &&... args) const

void	callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr, const bool show_trace=true) 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 ()

static void	callMooseError (MooseApp const app, const InputParameters &params, std::string msg, const bool with_prefix, const hit::Node node, const bool show_trace=true)

Public Attributes
	usingCombinedWarningSolutionWarnings

const ConsoleStream	_console

Static Public Attributes
static const std::string	unique_action_name_param

static const std::string	type_param

static const std::string	name_param

static const std::string	unique_name_param

static const std::string	app_param

static const std::string	moose_base_param

static const std::string	kokkos_object_param

static constexpr auto	SYSTEM

static constexpr auto	NAME

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

virtual void	addMesh () override
	Adds the mesh.

virtual void	addInitialConditions () override
	Adds the initial conditions.

virtual void	addSolutionVariables () override
	Adds the solution variables.

virtual void	addAuxVariables () override
	Adds aux variables.

virtual void	addMaterials () override
	Adds materials.

virtual void	addUserObjects () override
	Adds user objects.

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)

virtual void	addPreconditioner () override
	Adds the preconditioner.

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

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

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.

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

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

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

virtual void	addExecutioner ()
	Adds the executioner.

virtual void	addOutput ()
	Add output.

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

void	flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const

InvalidSolutionID	registerInvalidSolutionInternal (const std::string &message, const bool warning) 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

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

const VariableName	_rhoA_name
	rho*A variable name

const VariableName	_rhouA_name
	rhouA variable name

const VariableName	_rhoEA_name
	rhoEA variable name

const bool	_add_dg_kernel
	option to add DG kernel

const bool	_add_bc
	option to BC

const UserObjectName	_numerical_flux_name
	name of numerical flux user object being tested

const UserObjectName	_boundary_flux_name
	name of boundary flux user object being tested

const MooseEnum	_ic_option
	initial conditions option

const FunctionName &	_A_fn_name
	area function name

const bool	_use_slope_reconstruction
	option to use slope reconstruction

const std::string	_reconstruction_material_name
	reconstruction material name

const MaterialPropertyName	_direction_name
	direction material property name

const UserObjectName	_fp_name
	fluid properties object name

const std::string	_snes_test_err
	Finite differencing parameter.

bool	_default_use_transient_executioner
	Default for option to use a transient executioner.

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

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

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

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

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

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

const std::string	_fe_family
	Default FE family.

const std::string	_fe_order
	Default FE order.

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

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

Factory &	_factory

ActionFactory &	_action_factory

const std::string &	_type

const std::string &	_name

const InputParameters &	_pars

MooseApp &	_pg_moose_app

const std::string	_prefix

const Parallel::Communicator &	_communicator

Private Member Functions
bool	addRelationshipManager (Moose::RelationshipManagerType input_rm_type, const InputParameters &moose_object_pars, std::string rm_name, Moose::RelationshipManagerType rm_type, Moose::RelationshipManagerInputParameterCallback rm_input_parameter_func, Moose::RMSystemType sys_type=Moose::RMSystemType::NONE)

virtual std::string	meshPropertyPrefix (const std::string &data_name) const

const RestartableDataValue &	getMeshPropertyInternal (const std::string &data_name, const std::string &prefix) const

void	mooseErrorInternal (Args &&... args) const

Static Private Member Functions
static const hit::Node *	getHitNode (const InputParameters &params)

static std::string	messagePrefix (const InputParameters &params, const bool hit_prefix)

Private Attributes
const ParallelParamObject &	_parent

MooseApp &	_meta_data_app

const MooseObject *const	_meta_data_object

const MooseBase &	_si_moose_base

const FEProblemBase *	_si_problem

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 77 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 250 of file TestAction.C.

{
  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);
}

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

◆ 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 267 of file TestAction.C.

{
  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);
}

Referenced by addSolutionVariables().

◆ addExecutioner()

void TestAction::addExecutioner ( )

protectedvirtualinherited

Adds the executioner.

Definition at line 304 of file TestAction.C.

{
  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);
  }
}

Referenced by TestAction::addObjects().

◆ 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 283 of file TestAction.C.

{
  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);
}

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

◆ 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 107 of file TestAction.C.

{
  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);
  }
}

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

◆ 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 353 of file TestAction.C.

354{

355}

Referenced by TestAction::addObjects().

◆ 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 156 of file TestAction.C.

{
  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);
    }
  }
}

Referenced by TestAction::addObjects().

◆ 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 232 of file TestAction.C.

{
  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);
}

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

◆ 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.

{
  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);
  }
}

Referenced by addSolutionVariables().

◆ 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 TestAction::addAuxVariables(), ClosureTestAction::addAuxVariables(), and addAuxVariables().

◆ _fe_order

const std::string TestAction::_fe_order

protectedinherited

Default FE order.

Definition at line 152 of file TestAction.h.

Referenced by TestAction::addAuxVariables(), ClosureTestAction::addAuxVariables(), and 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

Private Member Functions

Static Private Member Functions

Private 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