Adds Axisymmetric Boussinesq bouyancy source term. More...

#include <axisym_boussinesq_buoyancy.h>

Inheritance diagram for GRINS::AxisymmetricBoussinesqBuoyancy:

Collaboration diagram for GRINS::AxisymmetricBoussinesqBuoyancy:

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

	~AxisymmetricBoussinesqBuoyancy ()

virtual void	init_context (AssemblyContext &context)
	Initialize context for added physics variables. More...

virtual void	element_time_derivative (bool compute_jacobian, AssemblyContext &context)
	Source term contribution for AxisymmetricBoussinesqBuoyancy. More...

Public Member Functions inherited from GRINS::Physics
	Physics (const GRINS::PhysicsName &physics_name, const GetPot &input)

virtual	~Physics ()

virtual void	init_variables (libMesh::FEMSystem *)
	Initialize variables for this physics. More...

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	set_time_evolving_vars (libMesh::FEMSystem *system)
	Set which variables are time evolving. 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	preassembly (MultiphysicsSystem &)
	Perform any necessary setup before element assembly begins. More...

virtual void	reinit (MultiphysicsSystem &)
	Any reinitialization that needs to be done. More...

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

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

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

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

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

virtual void	damping_residual (bool, AssemblyContext &)
	Damping matrix part(s) for element interiors. All boundary terms lie within the time_derivative part. More...

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

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

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

virtual void	compute_element_time_derivative_cache (AssemblyContext &)

virtual void	compute_side_time_derivative_cache (AssemblyContext &)

virtual void	compute_nonlocal_time_derivative_cache (AssemblyContext &)

virtual void	compute_element_constraint_cache (AssemblyContext &)

virtual void	compute_side_constraint_cache (AssemblyContext &)

virtual void	compute_nonlocal_constraint_cache (AssemblyContext &)

virtual void	compute_damping_residual_cache (AssemblyContext &)

virtual void	compute_mass_residual_cache (AssemblyContext &)

virtual void	compute_nonlocal_mass_residual_cache (AssemblyContext &)

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

ICHandlingBase *	get_ic_handler ()

Public Member Functions inherited from GRINS::ParameterUser
	ParameterUser (const std::string &user_name)

virtual	~ParameterUser ()

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	set_parameter (libMesh::ParsedFunction< libMesh::Number, libMesh::Gradient > &func, const GetPot &input, const std::string &func_param_name, const std::string &param_default)
	Each subclass can simultaneously read a parsed function from. More...

virtual void	set_parameter (libMesh::ParsedFEMFunction< libMesh::Number > &func, const GetPot &input, const std::string &func_param_name, const std::string &param_default)
	Each subclass can simultaneously read a parsed function from. More...

virtual void	move_parameter (const libMesh::Number &old_parameter, libMesh::Number &new_parameter)
	When cloning an object, we need to update parameter pointers. More...

virtual void	move_parameter (const libMesh::ParsedFunction< libMesh::Number, libMesh::Gradient > &old_func, libMesh::ParsedFunction< libMesh::Number, libMesh::Gradient > &new_func)
	When cloning an object, we need to update parameter pointers. More...

virtual void	move_parameter (const libMesh::ParsedFEMFunction< libMesh::Number > &old_func, libMesh::ParsedFEMFunction< libMesh::Number > &new_func)
	When cloning an object, we need to update parameter pointers. More...

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

Protected Attributes
const VelocityVariable &	_flow_vars

const PressureFEVariable &	_press_var

const PrimitiveTempFEVariables &	_temp_vars

libMesh::Number	_rho
	$\rho =$ density More...

libMesh::Number	_T_ref
	reference temperature More...

libMesh::Number	_beta_T
	$\beta_T =$ coefficient of thermal expansion More...

libMesh::Point	_g
	Gravitational vector. More...

Protected Attributes inherited from GRINS::Physics
const PhysicsName	_physics_name
	Name of the physics object. Used for reading physics specific inputs. More...

GRINS::ICHandlingBase *	_ic_handler

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

Private Member Functions
	AxisymmetricBoussinesqBuoyancy ()

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

Additional Inherited Members
Static Public Member Functions inherited from GRINS::Physics
static void	set_is_axisymmetric (bool is_axisymmetric)
	Set whether we should treat the problem as axisymmetric. More...

static bool	is_axisymmetric ()

Static Public Attributes inherited from GRINS::ParameterUser
static std::string	zero_vector_function = std::string("{0}")
	A parseable function string with LIBMESH_DIM components, all 0. More...

Protected Member Functions inherited from GRINS::Physics
libMesh::UniquePtr< libMesh::FEGenericBase< libMesh::Real > >	build_new_fe (const libMesh::Elem elem, const libMesh::FEGenericBase< libMesh::Real > fe, const libMesh::Point p)

void	parse_enabled_subdomains (const GetPot &input, const std::string &physics_name)

void	check_var_subdomain_consistency (const FEVariablesBase &var) const
	Check that var is enabled on at least the subdomains this Physics is. More...

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

static bool	_is_axisymmetric = false
	Caches whether we are solving an axisymmetric problem or not. More...

Detailed Description

Adds Axisymmetric Boussinesq bouyancy source term.

This class implements the Axisymmetric Boussinesq approximation for thermal buoyancy. Namely: $\mathbf{F} = -\rho_0 \beta_T \left( T - T_0 \right) \mathbf{g}$ where $\rho =$ density, $ T_0 = $ reference temperature, $\beta_T =$ coefficient of thermal expansion, and $\mathbf{g} =$ the gravitional vector. This source term is added to the governing flow equations through the element_time_derivative routine. This class requires an axisymmetric flow physics enabled and the AxisymmetricHeatTransfer physics class enabled.

Definition at line 57 of file axisym_boussinesq_buoyancy.h.

Constructor & Destructor Documentation

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

Definition at line 45 of file axisym_boussinesq_buoyancy.C.

References read_input_options().

     : Physics(physics_name, input),
       _flow_vars(GRINSPrivate::VariableWarehouse::get_variable_subclass<VelocityVariable>(VariablesParsing::velocity_variable_name(input,physics_name,VariablesParsing::PHYSICS))),
       _press_var(GRINSPrivate::VariableWarehouse::get_variable_subclass<PressureFEVariable>(VariablesParsing::press_variable_name(input,physics_name,VariablesParsing::PHYSICS))),
       _temp_vars(GRINSPrivate::VariableWarehouse::get_variable_subclass<PrimitiveTempFEVariables>(VariablesParsing::temp_variable_name(input,physics_name,VariablesParsing::PHYSICS)))
   {
     this->read_input_options(input);
   }

GRINS::AxisymmetricBoussinesqBuoyancy::~AxisymmetricBoussinesqBuoyancy ( )

inline

Definition at line 63 of file axisym_boussinesq_buoyancy.h.

63 {};

GRINS::AxisymmetricBoussinesqBuoyancy::AxisymmetricBoussinesqBuoyancy ( )

private

Member Function Documentation

void GRINS::AxisymmetricBoussinesqBuoyancy::element_time_derivative	(	bool	compute_jacobian,
		AssemblyContext &	context
	)

virtual

Source term contribution for AxisymmetricBoussinesqBuoyancy.

This is the main part of the class. This will add the source term to the AxisymmetricIncompNavierStokes class.

Reimplemented from GRINS::Physics.

Definition at line 80 of file axisym_boussinesq_buoyancy.C.

   {
     // The number of local degrees of freedom in each variable.
     const unsigned int n_u_dofs = context.get_dof_indices(_flow_vars.u()).size();
     const unsigned int n_T_dofs = context.get_dof_indices(_temp_vars.T()).size();
 
     // Element Jacobian * quadrature weights for interior integration.
     const std::vector<libMesh::Real> &JxW =
       context.get_element_fe(_flow_vars.u())->get_JxW();
 
     // The velocity shape functions at interior quadrature points.
     const std::vector<std::vector<libMesh::Real> >& vel_phi =
       context.get_element_fe(_flow_vars.u())->get_phi();
 
     // The temperature shape functions at interior quadrature points.
     const std::vector<std::vector<libMesh::Real> >& T_phi =
       context.get_element_fe(_temp_vars.T())->get_phi();
 
     // Physical location of the quadrature points
     const std::vector<libMesh::Point>& u_qpoint =
       context.get_element_fe(_flow_vars.u())->get_xyz();
 
     // Get residuals
     libMesh::DenseSubVector<libMesh::Number> &Fr = context.get_elem_residual(_flow_vars.u()); // R_{r}
     libMesh::DenseSubVector<libMesh::Number> &Fz = context.get_elem_residual(_flow_vars.v()); // R_{z}
 
     // Get Jacobians
     libMesh::DenseSubMatrix<libMesh::Number> &KrT = context.get_elem_jacobian(_flow_vars.u(), _temp_vars.T()); // R_{r},{T}
     libMesh::DenseSubMatrix<libMesh::Number> &KzT = context.get_elem_jacobian(_flow_vars.v(), _temp_vars.T()); // R_{z},{T}
 
     // Now we will build the element Jacobian and residual.
     // Constructing the residual requires the solution and its
     // gradient from the previous timestep.  This must be
     // calculated at each quadrature point by summing the
     // solution degree-of-freedom values by the appropriate
     // weight functions.
     unsigned int n_qpoints = context.get_element_qrule().n_points();
 
     for (unsigned int qp=0; qp != n_qpoints; qp++)
       {
         const libMesh::Number r = u_qpoint[qp](0);
 
         // Compute the solution & its gradient at the old Newton iterate.
         libMesh::Number T;
         T = context.interior_value(_temp_vars.T(), qp);
 
         // First, an i-loop over the velocity degrees of freedom.
         // We know that n_u_dofs == n_v_dofs so we can compute contributions
         // for both at the same time.
         for (unsigned int i=0; i != n_u_dofs; i++)
           {
             Fr(i) += -_rho*_beta_T*(T - _T_ref)*_g(0)*vel_phi[i][qp]*r*JxW[qp];
             Fz(i) += -_rho*_beta_T*(T - _T_ref)*_g(1)*vel_phi[i][qp]*r*JxW[qp];
 
             if (compute_jacobian && context.get_elem_solution_derivative())
               {
                 for (unsigned int j=0; j != n_T_dofs; j++)
                   {
                     const libMesh::Number val =
                       -_rho*_beta_T*vel_phi[i][qp]*T_phi[j][qp]*r*JxW[qp]
                       * context.get_elem_solution_derivative();
                     KrT(i,j) += val*_g(0);
                     KzT(i,j) += val*_g(1);
                   } // End j dof loop
               } // End compute_jacobian check
 
           } // End i dof loop
       } // End quadrature loop
   }

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

virtual

Initialize context for added physics variables.

Reimplemented from GRINS::Physics.

Definition at line 71 of file axisym_boussinesq_buoyancy.C.

References _flow_vars, _temp_vars, GRINS::PrimitiveTempFEVariables::T(), and GRINS::MultcomponentVectorVariable::u().

   {
     context.get_element_fe(_flow_vars.u())->get_JxW();
     context.get_element_fe(_flow_vars.u())->get_phi();
     context.get_element_fe(_flow_vars.u())->get_xyz();
     context.get_element_fe(_temp_vars.T())->get_phi();
   }

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

private

Read options from GetPot input file.

Definition at line 55 of file axisym_boussinesq_buoyancy.C.

References _beta_T, _g, _rho, _T_ref, GRINS::PhysicsNaming::axisymmetric_boussinesq_buoyancy(), and GRINS::ParameterUser::set_parameter().

Referenced by AxisymmetricBoussinesqBuoyancy().

   {
     this->set_parameter
       (_rho, input,
        "Physics/"+PhysicsNaming::axisymmetric_boussinesq_buoyancy()+"/rho_ref", 1.0);
 
     this->set_parameter
       (_T_ref, input, "Physics/"+PhysicsNaming::axisymmetric_boussinesq_buoyancy()+"/T_ref", 1.0);
 
     this->set_parameter
       (_beta_T, input, "Physics/"+PhysicsNaming::axisymmetric_boussinesq_buoyancy()+"/beta_T", 1.0);
 
     _g(0) = input("Physics/"+PhysicsNaming::axisymmetric_boussinesq_buoyancy()+"/g", 0.0, 0 );
     _g(1) = input("Physics/"+PhysicsNaming::axisymmetric_boussinesq_buoyancy()+"/g", 0.0, 1 );
   }