Physics class for arbitrary scalar-valued ODEs. More...

#include <scalar_ode.h>

Inheritance diagram for GRINS::ScalarODE:

Collaboration diagram for GRINS::ScalarODE:

Public Member Functions
	ScalarODE (const std::string &physics_name, const GetPot &input)

	~ScalarODE ()

virtual void	init_variables (libMesh::FEMSystem *system)
	Initialization of variables. More...

virtual void	set_time_evolving_vars (libMesh::FEMSystem *system)
	Sets scalar variable(s) to be time-evolving. More...

virtual void	read_input_options (const GetPot &input)
	Read options from GetPot input file. More...

virtual void	init_context (AssemblyContext &context)
	Prepare the context for evaluations. More...

virtual void	nonlocal_time_derivative (bool compute_jacobian, AssemblyContext &context, CachedValues &cache)
	Time dependent part(s) of physics for scalar variables. More...

virtual void	nonlocal_constraint (bool compute_jacobian, AssemblyContext &context, CachedValues &cache)
	Constraint part(s) of physics for scalar variables. More...

virtual void	nonlocal_mass_residual (bool compute_jacobian, AssemblyContext &context, CachedValues &cache)
	Mass matrix part(s) for scalar variables. More...

VariableIndex	scalar_ode_var () const

virtual bool	enabled_on_elem (const libMesh::Elem *elem)
	Find if current physics is active on supplied element. More...

void	set_is_steady (bool is_steady)
	Sets whether this physics is to be solved with a steady solver or not. More...

bool	is_steady () const
	Returns whether or not this physics is being solved with a steady solver. More...

virtual void	auxiliary_init (MultiphysicsSystem &system)
	Any auxillary initialization a Physics class may need. More...

virtual void	register_postprocessing_vars (const GetPot &input, PostProcessedQuantities< libMesh::Real > &postprocessing)
	Register name of postprocessed quantity with PostProcessedQuantities. More...

virtual void	element_time_derivative (bool compute_jacobian, AssemblyContext &context, CachedValues &cache)
	Time dependent part(s) of physics for element interiors. More...

virtual void	side_time_derivative (bool compute_jacobian, AssemblyContext &context, CachedValues &cache)
	Time dependent part(s) of physics for boundaries of elements on the domain boundary. More...

virtual void	element_constraint (bool compute_jacobian, AssemblyContext &context, CachedValues &cache)
	Constraint part(s) of physics for element interiors. More...

virtual void	side_constraint (bool compute_jacobian, AssemblyContext &context, CachedValues &cache)
	Constraint part(s) of physics for boundaries of elements on the domain boundary. More...

virtual void	mass_residual (bool compute_jacobian, AssemblyContext &context, CachedValues &cache)
	Mass matrix part(s) for element interiors. All boundary terms lie within the time_derivative part. More...

void	init_bcs (libMesh::FEMSystem *system)

void	init_ics (libMesh::FEMSystem *system, libMesh::CompositeFunction< libMesh::Number > &all_ics)

void	attach_neumann_bound_func (GRINS::NBCContainer &neumann_bcs)

void	attach_dirichlet_bound_func (const GRINS::DBCContainer &dirichlet_bc)

virtual void	compute_element_time_derivative_cache (const AssemblyContext &context, CachedValues &cache)

virtual void	compute_side_time_derivative_cache (const AssemblyContext &context, CachedValues &cache)

virtual void	compute_nonlocal_time_derivative_cache (const AssemblyContext &context, CachedValues &cache)

virtual void	compute_element_constraint_cache (const AssemblyContext &context, CachedValues &cache)

virtual void	compute_side_constraint_cache (const AssemblyContext &context, CachedValues &cache)

virtual void	compute_nonlocal_constraint_cache (const AssemblyContext &context, CachedValues &cache)

virtual void	compute_mass_residual_cache (const AssemblyContext &context, CachedValues &cache)

virtual void	compute_nonlocal_mass_residual_cache (const AssemblyContext &context, CachedValues &cache)

virtual void	compute_postprocessed_quantity (unsigned int quantity_index, const AssemblyContext &context, const libMesh::Point &point, libMesh::Real &value)

BCHandlingBase *	get_bc_handler ()

ICHandlingBase *	get_ic_handler ()

virtual void	set_parameter (libMesh::Number &param_variable, const GetPot &input, const std::string &param_name, libMesh::Number param_default)
	Each subclass can simultaneously read a parameter value from. More...

virtual void	register_parameter (const std::string &param_name, libMesh::ParameterMultiPointer< libMesh::Number > &param_pointer) const
	Each subclass will register its copy of an independent. More...

Protected Attributes
const PhysicsName	_physics_name
	Name of the physics object. Used for reading physics specific inputs. More...

GRINS::BCHandlingBase *	_bc_handler

GRINS::ICHandlingBase *	_ic_handler

std::set< libMesh::subdomain_id_type >	_enabled_subdomains
	Subdomains on which the current Physics class is enabled. More...

bool	_is_axisymmetric

Static Protected Attributes
static bool	_is_steady = false
	Caches whether or not the solver that's being used is steady or not. More...

Private Member Functions
	ScalarODE ()

Private Attributes
std::string	time_deriv_function_string

std::string	constraint_function_string

std::string	mass_residual_function_string

libMesh::AutoPtr< libMesh::FEMFunctionBase< libMesh::Number > >	time_deriv_function

libMesh::AutoPtr< libMesh::FEMFunctionBase< libMesh::Number > >	constraint_function

libMesh::AutoPtr< libMesh::FEMFunctionBase< libMesh::Number > >	mass_residual_function

libMesh::Number	_order

libMesh::Number	_epsilon

VariableIndex	_scalar_ode_var

std::string	_scalar_ode_var_name

Detailed Description

Physics class for arbitrary scalar-valued ODEs.

Definition at line 50 of file scalar_ode.h.

Constructor & Destructor Documentation

GRINS::ScalarODE::ScalarODE	(	const std::string &	physics_name,
		const GetPot &	input
	)

Definition at line 40 of file scalar_ode.C.

References GRINS::Physics::_ic_handler, and read_input_options().

     : Physics(physics_name, input), _order(1), _epsilon(1e-6)
   {
     this->read_input_options(input);
 
     this->_ic_handler = new GenericICHandler( physics_name, input );
 
     return;
   }

GRINS::ScalarODE::~ScalarODE ( )

Definition at line 50 of file scalar_ode.C.

   {
     return;
   }

GRINS::ScalarODE::ScalarODE ( )

private

Member Function Documentation

void GRINS::Physics::attach_dirichlet_bound_func ( const GRINS::DBCContainer & dirichlet_bc )

inherited

Definition at line 150 of file physics.C.

References GRINS::Physics::_bc_handler, and GRINS::BCHandlingBase::attach_dirichlet_bound_func().

   {
     _bc_handler->attach_dirichlet_bound_func( dirichlet_bc );
     return;
   }

void GRINS::Physics::attach_neumann_bound_func ( GRINS::NBCContainer & neumann_bcs )

inherited

Definition at line 144 of file physics.C.

References GRINS::Physics::_bc_handler, and GRINS::BCHandlingBase::attach_neumann_bound_func().

   {
     _bc_handler->attach_neumann_bound_func( neumann_bcs );
     return;
   }

void GRINS::Physics::auxiliary_init ( MultiphysicsSystem & system )

virtualinherited

Any auxillary initialization a Physics class may need.

This is called after all variables are added, so this method can safely query the MultiphysicsSystem about variable information.

Definition at line 113 of file physics.C.

   {
     return;
   }

void GRINS::Physics::compute_element_constraint_cache	(	const AssemblyContext &	context,
		CachedValues &	cache
	)

virtualinherited

Definition at line 185 of file physics.C.

Referenced by GRINS::MultiphysicsSystem::element_constraint().

   {
     return;
   }

void GRINS::Physics::compute_element_time_derivative_cache	(	const AssemblyContext &	context,
		CachedValues &	cache
	)

virtualinherited

Reimplemented in GRINS::LowMachNavierStokes< Viscosity, SpecificHeat, ThermalConductivity >, and GRINS::ReactingLowMachNavierStokes< Mixture, Evaluator >.

Definition at line 167 of file physics.C.

Referenced by GRINS::MultiphysicsSystem::element_time_derivative().

   {
     return;
   }

void GRINS::Physics::compute_mass_residual_cache	(	const AssemblyContext &	context,
		CachedValues &	cache
	)

virtualinherited

Definition at line 203 of file physics.C.

Referenced by GRINS::MultiphysicsSystem::mass_residual().

   {
     return;
   }

void GRINS::Physics::compute_nonlocal_constraint_cache	(	const AssemblyContext &	context,
		CachedValues &	cache
	)

virtualinherited

Definition at line 197 of file physics.C.

Referenced by GRINS::MultiphysicsSystem::nonlocal_constraint().

   {
     return;
   }

void GRINS::Physics::compute_nonlocal_mass_residual_cache	(	const AssemblyContext &	context,
		CachedValues &	cache
	)

virtualinherited

Definition at line 209 of file physics.C.

Referenced by GRINS::MultiphysicsSystem::nonlocal_mass_residual().

   {
     return;
   }

void GRINS::Physics::compute_nonlocal_time_derivative_cache	(	const AssemblyContext &	context,
		CachedValues &	cache
	)

virtualinherited

Definition at line 179 of file physics.C.

Referenced by GRINS::MultiphysicsSystem::nonlocal_time_derivative().

   {
     return;
   }

void GRINS::Physics::compute_postprocessed_quantity	(	unsigned int	quantity_index,
		const AssemblyContext &	context,
		const libMesh::Point &	point,
		libMesh::Real &	value
	)

virtualinherited

Reimplemented in GRINS::IncompressibleNavierStokes< Viscosity >, GRINS::LowMachNavierStokes< Viscosity, SpecificHeat, ThermalConductivity >, GRINS::ReactingLowMachNavierStokes< Mixture, Evaluator >, GRINS::HeatTransfer< Conductivity >, GRINS::ParsedVelocitySource< Viscosity >, GRINS::VelocityPenalty< Viscosity >, GRINS::ElasticMembrane< StressStrainLaw >, and GRINS::ElasticCable< StressStrainLaw >.

Definition at line 271 of file physics.C.

   {
     return;
   }

void GRINS::Physics::compute_side_constraint_cache	(	const AssemblyContext &	context,
		CachedValues &	cache
	)

virtualinherited

Definition at line 191 of file physics.C.

Referenced by GRINS::MultiphysicsSystem::side_constraint().

   {
     return;
   }

void GRINS::Physics::compute_side_time_derivative_cache	(	const AssemblyContext &	context,
		CachedValues &	cache
	)

virtualinherited

Reimplemented in GRINS::ReactingLowMachNavierStokes< Mixture, Evaluator >.

Definition at line 173 of file physics.C.

Referenced by GRINS::MultiphysicsSystem::side_time_derivative().

   {
     return;
   }

void GRINS::Physics::element_constraint	(	bool	compute_jacobian,
		AssemblyContext &	context,
		CachedValues &	cache
	)

virtualinherited

Constraint part(s) of physics for element interiors.

Definition at line 236 of file physics.C.

Referenced by GRINS::MultiphysicsSystem::element_constraint().

   {
     return;
   }

void GRINS::Physics::element_time_derivative	(	bool	compute_jacobian,
		AssemblyContext &	context,
		CachedValues &	cache
	)

virtualinherited

Time dependent part(s) of physics for element interiors.

Definition at line 215 of file physics.C.

Referenced by GRINS::MultiphysicsSystem::element_time_derivative().

   {
     return;
   }

bool GRINS::Physics::enabled_on_elem ( const libMesh::Elem * elem )

virtualinherited

Find if current physics is active on supplied element.

Definition at line 83 of file physics.C.

References GRINS::Physics::_enabled_subdomains.

   {
     // Check if enabled_subdomains flag has been set and if we're
     // looking at a real element (rather than a nonlocal evaluation)
     if( !elem || _enabled_subdomains.empty() ) 
       return true;
   
     // Check if current physics is enabled on elem
     if( _enabled_subdomains.find( elem->subdomain_id() ) == _enabled_subdomains.end() )
       return false;
 
     return true;
   }

BCHandlingBase * GRINS::Physics::get_bc_handler ( )

inlineinherited

Definition at line 282 of file physics.h.

References GRINS::Physics::_bc_handler.

   {
     return _bc_handler;
   }

ICHandlingBase * GRINS::Physics::get_ic_handler ( )

inlineinherited

Definition at line 288 of file physics.h.

References GRINS::Physics::_ic_handler.

   {
     return _ic_handler;
   }

void GRINS::Physics::init_bcs ( libMesh::FEMSystem * system )

inherited

Definition at line 118 of file physics.C.

References GRINS::Physics::_bc_handler, GRINS::BCHandlingBase::init_bc_data(), GRINS::BCHandlingBase::init_dirichlet_bc_func_objs(), GRINS::BCHandlingBase::init_dirichlet_bcs(), and GRINS::BCHandlingBase::init_periodic_bcs().

   {
     // Only need to init BC's if the physics actually created a handler
     if( _bc_handler )
       {
         _bc_handler->init_bc_data( *system );
         _bc_handler->init_dirichlet_bcs( system );
         _bc_handler->init_dirichlet_bc_func_objs( system );
         _bc_handler->init_periodic_bcs( system );
       }
 
     return;
   }

void GRINS::ScalarODE::init_context ( AssemblyContext & context )

virtual

Prepare the context for evaluations.

Reimplemented from GRINS::Physics.

Definition at line 116 of file scalar_ode.C.

References constraint_function, mass_residual_function, and time_deriv_function.

   {
     mass_residual_function->init_context(context);
     time_deriv_function->init_context(context);
     constraint_function->init_context(context);
   }

void GRINS::Physics::init_ics	(	libMesh::FEMSystem *	system,
		libMesh::CompositeFunction< libMesh::Number > &	all_ics
	)

inherited

Definition at line 133 of file physics.C.

References GRINS::Physics::_ic_handler, and GRINS::ICHandlingBase::init_ic_data().

   {
     if( _ic_handler )
       {
         _ic_handler->init_ic_data( *system, all_ics );
       }
 
     return;
   }

void GRINS::ScalarODE::init_variables ( libMesh::FEMSystem * system )

virtual

Initialization of variables.

Add scalar variable(s) to system.

Implements GRINS::Physics.

Definition at line 55 of file scalar_ode.C.

References _order, _scalar_ode_var, and _scalar_ode_var_name.

   {
     this->_scalar_ode_var = system->add_variable(_scalar_ode_var_name,
                                                  libMesh::Order(this->_order),
                                                  libMesh::SCALAR);
   }

bool GRINS::Physics::is_steady ( ) const

inherited

Returns whether or not this physics is being solved with a steady solver.

Definition at line 103 of file physics.C.

References GRINS::Physics::_is_steady.

Referenced by GRINS::Physics::set_is_steady().

   {
     return _is_steady;
   }

void GRINS::Physics::mass_residual	(	bool	compute_jacobian,
		AssemblyContext &	context,
		CachedValues &	cache
	)

virtualinherited

Mass matrix part(s) for element interiors. All boundary terms lie within the time_derivative part.

Reimplemented in GRINS::AxisymmetricHeatTransfer< Conductivity >, GRINS::IncompressibleNavierStokes< Viscosity >, GRINS::SpalartAllmaras< Viscosity >, GRINS::LowMachNavierStokes< Viscosity, SpecificHeat, ThermalConductivity >, GRINS::HeatTransfer< Conductivity >, GRINS::ReactingLowMachNavierStokes< Mixture, Evaluator >, GRINS::Stokes< Viscosity >, GRINS::BoussinesqBuoyancySPGSMStabilization< Viscosity >, GRINS::ElasticMembrane< StressStrainLaw >, GRINS::ElasticCable< StressStrainLaw >, GRINS::HeatConduction< Conductivity >, GRINS::IncompressibleNavierStokesAdjointStabilization< Viscosity >, GRINS::IncompressibleNavierStokesSPGSMStabilization< Viscosity >, GRINS::SpalartAllmarasSPGSMStabilization< Viscosity >, GRINS::HeatTransferAdjointStabilization< Conductivity >, GRINS::HeatTransferSPGSMStabilization< Conductivity >, GRINS::LowMachNavierStokesBraackStabilization< Viscosity, SpecificHeat, ThermalConductivity >, GRINS::LowMachNavierStokesVMSStabilization< Viscosity, SpecificHeat, ThermalConductivity >, and GRINS::LowMachNavierStokesSPGSMStabilization< Viscosity, SpecificHeat, ThermalConductivity >.

Definition at line 257 of file physics.C.

Referenced by GRINS::MultiphysicsSystem::mass_residual().

   {
     return;
   }

void GRINS::ScalarODE::nonlocal_constraint	(	bool	compute_jacobian,
		AssemblyContext &	context,
		CachedValues &	cache
	)

virtual

Constraint part(s) of physics for scalar variables.

Reimplemented from GRINS::Physics.

Definition at line 216 of file scalar_ode.C.

References _epsilon, and _scalar_ode_var.

   {
     libMesh::DenseSubMatrix<libMesh::Number> &Kss =
             context.get_elem_jacobian(_scalar_ode_var, _scalar_ode_var); // R_{s},{s}
 
     libMesh::DenseSubVector<libMesh::Number> &Fs =
             context.get_elem_residual(_scalar_ode_var); // R_{s}
 
     const libMesh::Number constraint =
       (*constraint_function)(context, libMesh::Point(0),
                              context.get_time());
 
     Fs(0) += constraint;
 
     if (compute_jacobian)
       {
         // FIXME: we should replace this hacky FDM with a hook to the
         // AD fparser stuff
         libMesh::DenseSubVector<libMesh::Number> &Us =
           const_cast<libMesh::DenseSubVector<libMesh::Number>&>
             (context.get_elem_solution(_scalar_ode_var)); // U_{s}
 
         const libMesh::Number s = Us(0);
         Us(0) = s + this->_epsilon;
         libMesh::Number constraint_jacobian =
           (*constraint_function)(context, libMesh::Point(0),
                                  context.get_time());
 
         Us(0) = s - this->_epsilon;
         constraint_jacobian -=
           (*constraint_function)(context, libMesh::Point(0),
                                  context.get_time());
            
         Us(0) = s;
         constraint_jacobian /= (2*this->_epsilon);
 
         Kss(0,0) += constraint_jacobian *
           context.get_elem_solution_derivative();
       }
 
     return;
   }

void GRINS::ScalarODE::nonlocal_mass_residual	(	bool	compute_jacobian,
		AssemblyContext &	context,
		CachedValues &	cache
	)

virtual

Mass matrix part(s) for scalar variables.

Reimplemented from GRINS::Physics.

Definition at line 170 of file scalar_ode.C.

References _epsilon, and _scalar_ode_var.

   {
     libMesh::DenseSubMatrix<libMesh::Number> &Kss =
             context.get_elem_jacobian(_scalar_ode_var, _scalar_ode_var); // R_{s},{s}
 
     libMesh::DenseSubVector<libMesh::Number> &Fs =
             context.get_elem_residual(_scalar_ode_var); // R_{s}
 
     const libMesh::Number mass_res =
       (*mass_residual_function)(context, libMesh::Point(0),
                                 context.get_time());
 
     Fs(0) -= mass_res;
 
     if (compute_jacobian)
       {
         // FIXME: we should replace this hacky FDM with a hook to the
         // AD fparser stuff
         libMesh::DenseSubVector<libMesh::Number> &Us =
           const_cast<libMesh::DenseSubVector<libMesh::Number>&>
             (context.get_elem_solution_rate(_scalar_ode_var)); // U_{s}
 
         const libMesh::Number s = Us(0);
         Us(0) = s + this->_epsilon;
         libMesh::Number mass_residual_jacobian =
           (*mass_residual_function)(context, libMesh::Point(0),
                                     context.get_time());
 
         Us(0) = s - this->_epsilon;
         mass_residual_jacobian -=
           (*mass_residual_function)(context, libMesh::Point(0),
                                     context.get_time());
            
         Us(0) = s;
         mass_residual_jacobian /= (2*this->_epsilon);
 
         Kss(0,0) -= mass_residual_jacobian *
           context.get_elem_solution_rate_derivative();
       }
 
     return;
   }

void GRINS::ScalarODE::nonlocal_time_derivative	(	bool	compute_jacobian,
		AssemblyContext &	context,
		CachedValues &	cache
	)

virtual

Time dependent part(s) of physics for scalar variables.

Reimplemented from GRINS::Physics.

Definition at line 124 of file scalar_ode.C.

References _epsilon, and _scalar_ode_var.

   {
     libMesh::DenseSubMatrix<libMesh::Number> &Kss =
             context.get_elem_jacobian(_scalar_ode_var, _scalar_ode_var); // R_{s},{s}
 
     libMesh::DenseSubVector<libMesh::Number> &Fs =
             context.get_elem_residual(_scalar_ode_var); // R_{s}
 
     const libMesh::Number time_deriv =
       (*time_deriv_function)(context, libMesh::Point(0),
                              context.get_time());
 
     Fs(0) += time_deriv;
 
     if (compute_jacobian)
       {
         // FIXME: we should replace this hacky FDM with a hook to the
         // AD fparser stuff
         libMesh::DenseSubVector<libMesh::Number> &Us =
           const_cast<libMesh::DenseSubVector<libMesh::Number>&>
             (context.get_elem_solution(_scalar_ode_var)); // U_{s}
 
         const libMesh::Number s = Us(0);
         Us(0) = s + this->_epsilon;
         libMesh::Number time_deriv_jacobian =
           (*time_deriv_function)(context, libMesh::Point(0),
                                  context.get_time());
 
         Us(0) = s - this->_epsilon;
         time_deriv_jacobian -=
           (*time_deriv_function)(context, libMesh::Point(0),
                                  context.get_time());
            
         Us(0) = s;
         time_deriv_jacobian /= (2*this->_epsilon);
 
         Kss(0,0) += time_deriv_jacobian *
           context.get_elem_solution_derivative();
       }
 
     return;
   }

void GRINS::ScalarODE::read_input_options ( const GetPot & input )

virtual

Read options from GetPot input file.

Reimplemented from GRINS::Physics.

Definition at line 84 of file scalar_ode.C.

References _epsilon, _order, _scalar_ode_var_name, constraint_function_string, mass_residual_function_string, GRINS::scalar_ode, GRINS::scalar_ode_var_name_default, and time_deriv_function_string.

Referenced by ScalarODE().

   {
     this->time_deriv_function_string =
       input("Physics/"+scalar_ode+"/time_deriv",
         std::string("0"));
 
     if (this->time_deriv_function_string == "0")
       std::cout << "Warning! Zero time_deriv function specified!" << std::endl;
 
     this->mass_residual_function_string =
       input("Physics/"+scalar_ode+"/mass_residual",
         std::string("0"));
 
     if (this->mass_residual_function_string == "0")
       std::cout << "Warning! Zero mass_residual function specified!" << std::endl;
 
     this->constraint_function_string =
       input("Physics/"+scalar_ode+"/constraint",
         std::string("0"));
 
     if (this->constraint_function_string == "0")
       std::cout << "Warning! Zero constraint function specified!" << std::endl;
 
     this->_epsilon = input("Physics/"+scalar_ode+"/epsilon", 1e-6);
 
     this->_order = input("Physics/"+scalar_ode+"/order", 1);
 
     _scalar_ode_var_name = input("Physics/VariableNames/scalar_ode",
                                  scalar_ode_var_name_default);
   }

void GRINS::ParameterUser::register_parameter	(	const std::string &	param_name,
		libMesh::ParameterMultiPointer< libMesh::Number > &	param_pointer
	)		const

virtualinherited

Each subclass will register its copy of an independent.

Reimplemented in GRINS::AxisymmetricHeatTransfer< Conductivity >, GRINS::LowMachNavierStokesBase< Viscosity, SpecificHeat, ThermalConductivity >, GRINS::IncompressibleNavierStokesBase< Viscosity >, GRINS::BoussinesqBuoyancySPGSMStabilization< Viscosity >, GRINS::HeatConduction< Conductivity >, GRINS::HeatTransferBase< Conductivity >, and GRINS::BoussinesqBuoyancyAdjointStabilization< Viscosity >.

Definition at line 50 of file parameter_user.C.

Referenced by GRINS::BoussinesqBuoyancyAdjointStabilization< Viscosity >::register_parameter(), GRINS::HeatTransferBase< Conductivity >::register_parameter(), GRINS::HeatConduction< Conductivity >::register_parameter(), GRINS::BoussinesqBuoyancySPGSMStabilization< Viscosity >::register_parameter(), GRINS::IncompressibleNavierStokesBase< Viscosity >::register_parameter(), GRINS::LowMachNavierStokesBase< Viscosity, SpecificHeat, ThermalConductivity >::register_parameter(), and GRINS::AxisymmetricHeatTransfer< Conductivity >::register_parameter().

   {
     std::map<std::string, libMesh::Number*>::const_iterator it =
       _my_parameters.find(param_name);
 
     if (it != _my_parameters.end())
       {
         std::cout << _my_name << " uses parameter " << param_name
                   << std::endl;
         param_pointer.push_back(it->second);
       }
   }

void GRINS::Physics::register_postprocessing_vars	(	const GetPot &	input,
		PostProcessedQuantities< libMesh::Real > &	postprocessing
	)

virtualinherited

Register name of postprocessed quantity with PostProcessedQuantities.

Each Physics class will need to cache an unsigned int corresponding to each postprocessed quantity. This will be used in computing the values and putting them in the CachedVariables object.

Reimplemented in GRINS::ParsedVelocitySource< Viscosity >, GRINS::VelocityPenalty< Viscosity >, GRINS::IncompressibleNavierStokes< Viscosity >, GRINS::LowMachNavierStokes< Viscosity, SpecificHeat, ThermalConductivity >, GRINS::HeatTransfer< Conductivity >, GRINS::ElasticCable< StressStrainLaw >, GRINS::ElasticMembrane< StressStrainLaw >, and GRINS::ReactingLowMachNavierStokes< Mixture, Evaluator >.

Definition at line 161 of file physics.C.

   {
     return;
   }

VariableIndex GRINS::ScalarODE::scalar_ode_var ( ) const

inline

Definition at line 90 of file scalar_ode.h.

References _scalar_ode_var.

Referenced by set_time_evolving_vars().

90 { return _scalar_ode_var; }

GRINS::ScalarODE::_scalar_ode_var

VariableIndex _scalar_ode_var

Definition: scalar_ode.h:114

void GRINS::Physics::set_is_steady ( bool is_steady )

inherited

Sets whether this physics is to be solved with a steady solver or not.

Since the member variable is static, only needs to be called on a single physics.

Definition at line 97 of file physics.C.

References GRINS::Physics::_is_steady, and GRINS::Physics::is_steady().

   {
     _is_steady = is_steady;
     return;
   }

void GRINS::ParameterUser::set_parameter	(	libMesh::Number &	param_variable,
		const GetPot &	input,
		const std::string &	param_name,
		libMesh::Number	param_default
	)

virtualinherited

Each subclass can simultaneously read a parameter value from.

Definition at line 35 of file parameter_user.C.

References GRINS::ParameterUser::_my_name, and GRINS::ParameterUser::_my_parameters.

Referenced by GRINS::AveragedFanAdjointStabilization< Viscosity >::AveragedFanAdjointStabilization(), GRINS::AveragedTurbineAdjointStabilization< Viscosity >::AveragedTurbineAdjointStabilization(), GRINS::BoussinesqBuoyancyAdjointStabilization< Viscosity >::BoussinesqBuoyancyAdjointStabilization(), GRINS::BoussinesqBuoyancyBase::BoussinesqBuoyancyBase(), GRINS::BoussinesqBuoyancySPGSMStabilization< Viscosity >::BoussinesqBuoyancySPGSMStabilization(), GRINS::ConstantConductivity::ConstantConductivity(), GRINS::ConstantPrandtlConductivity::ConstantPrandtlConductivity(), GRINS::ConstantSourceFunction::ConstantSourceFunction(), GRINS::ConstantSourceTerm::ConstantSourceTerm(), GRINS::ConstantSpecificHeat::ConstantSpecificHeat(), GRINS::ConstantViscosity::ConstantViscosity(), GRINS::ElasticCable< StressStrainLaw >::ElasticCable(), GRINS::ElasticCableConstantGravity::ElasticCableConstantGravity(), GRINS::ElasticMembrane< StressStrainLaw >::ElasticMembrane(), GRINS::ElasticMembraneConstantPressure::ElasticMembraneConstantPressure(), GRINS::HeatConduction< Conductivity >::HeatConduction(), GRINS::HeatTransferBase< Conductivity >::HeatTransferBase(), GRINS::IncompressibleNavierStokesBase< Viscosity >::IncompressibleNavierStokesBase(), GRINS::AverageNusseltNumber::init(), GRINS::MooneyRivlin::MooneyRivlin(), GRINS::ReactingLowMachNavierStokesBase< Mixture, Evaluator >::ReactingLowMachNavierStokesBase(), GRINS::HookesLaw1D::read_input_options(), GRINS::HookesLaw::read_input_options(), GRINS::AxisymmetricBoussinesqBuoyancy::read_input_options(), and GRINS::VelocityDragAdjointStabilization< Viscosity >::VelocityDragAdjointStabilization().

   {
     param_variable = input(param_name, param_default);
 
     libmesh_assert_msg(!_my_parameters.count(param_name),
       "ERROR: " << _my_name << " double-registered parameter " <<
       param_name);
 
     _my_parameters[param_name] = &param_variable;
   }

void GRINS::ScalarODE::set_time_evolving_vars ( libMesh::FEMSystem * system )

virtual

Sets scalar variable(s) to be time-evolving.

Reimplemented from GRINS::Physics.

Definition at line 62 of file scalar_ode.C.

References constraint_function, constraint_function_string, mass_residual_function, mass_residual_function_string, scalar_ode_var(), time_deriv_function, and time_deriv_function_string.

   {
     system->time_evolving(this->scalar_ode_var());
 
     // FIXME: this doesn't fit here at all, but it's the only time
     // we've clearly got a System to grab hold of with all it's
     // variables initialized.
 
     this->time_deriv_function.reset
       (new libMesh::ParsedFEMFunction<libMesh::Number>
        (*system, this->time_deriv_function_string));
 
     this->mass_residual_function.reset
       (new libMesh::ParsedFEMFunction<libMesh::Number>
        (*system, this->mass_residual_function_string));
 
     this->constraint_function.reset
       (new libMesh::ParsedFEMFunction<libMesh::Number>
        (*system, this->constraint_function_string));
   }