GRINS::HeatTransferStabilizationHelper Class Reference

#include <heat_transfer_stab_helper.h>

Inheritance diagram for GRINS::HeatTransferStabilizationHelper:

Collaboration diagram for GRINS::HeatTransferStabilizationHelper:

Public Member Functions
	HeatTransferStabilizationHelper (const std::string &helper_name, const GetPot &input)
	~HeatTransferStabilizationHelper ()
libMesh::Real	compute_res_energy_steady (AssemblyContext &context, unsigned int qp, const libMesh::Real rho, const libMesh::Real Cp, const libMesh::Real k) const
void	compute_res_energy_steady_and_derivs (AssemblyContext &context, unsigned int qp, const libMesh::Real rho, const libMesh::Real Cp, const libMesh::Real k, libMesh::Real &res, libMesh::Real &d_res_dT, libMesh::Gradient &d_res_dgradT, libMesh::Tensor &d_res_dhessT, libMesh::Gradient &d_res_dU) const
libMesh::Real	compute_res_energy_transient (AssemblyContext &context, unsigned int qp, const libMesh::Real rho, const libMesh::Real Cp) const
void	compute_res_energy_transient_and_derivs (AssemblyContext &context, unsigned int qp, const libMesh::Real rho, const libMesh::Real Cp, libMesh::Real &res, libMesh::Real &d_res_dTdot) const
libMesh::Real	compute_tau_energy (AssemblyContext &c, libMesh::RealTensor &G, libMesh::Real rho, libMesh::Real cp, libMesh::Real k, libMesh::Gradient U, bool is_steady) const
void	compute_tau_energy_and_derivs (AssemblyContext &c, libMesh::RealTensor &G, libMesh::Real rho, libMesh::Real cp, libMesh::Real k, libMesh::Gradient U, libMesh::Real &tau_E, libMesh::Real &d_tau_E_d_rho, libMesh::Gradient &d_tau_E_d_U, bool is_steady) const
Public Member Functions inherited from GRINS::StabilizationHelper
	StabilizationHelper (const std::string &helper_name)
	~StabilizationHelper ()
libMesh::RealGradient	compute_g (libMesh::FEBase *fe, AssemblyContext &c, unsigned int qp) const
libMesh::RealTensor	compute_G (libMesh::FEBase *fe, AssemblyContext &c, unsigned int qp) const
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
libMesh::Real	_C
libMesh::Real	_tau_factor
const PrimitiveTempFEVariables &	_temp_vars
const VelocityFEVariables &	_flow_vars
const PressureFEVariable &	_press_var

Additional Inherited Members
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...

Detailed Description

Definition at line 42 of file heat_transfer_stab_helper.h.

Constructor & Destructor Documentation

GRINS::HeatTransferStabilizationHelper::HeatTransferStabilizationHelper	(	const std::string &	helper_name,
		const GetPot &	input
	)

Definition at line 45 of file heat_transfer_stab_helper.C.

    : StabilizationHelper(helper_name),
      _C(1),
      _tau_factor(0.5),
      _temp_vars(GRINSPrivate::VariableWarehouse::get_variable_subclass<PrimitiveTempFEVariables>(VariablesParsing::temperature_section())),
      _flow_vars(GRINSPrivate::VariableWarehouse::get_variable_subclass<VelocityFEVariables>(VariablesParsing::velocity_section())),
      _press_var(GRINSPrivate::VariableWarehouse::get_variable_subclass<PressureFEVariable>(VariablesParsing::pressure_section()))
  {
    if (input.have_variable("Stabilization/tau_constant_T"))
      this->set_parameter
        (_C, input, "Stabilization/tau_constant_T", _C );
    else
      this->set_parameter
        (_C, input, "Stabilization/tau_constant", _C );

    if (input.have_variable("Stabilization/tau_factor_T"))
      this->set_parameter
        (_tau_factor, input, "Stabilization/tau_factor_T", _tau_factor );
    else
      this->set_parameter
        (_tau_factor, input, "Stabilization/tau_factor", _tau_factor );
  }

GRINS::HeatTransferStabilizationHelper::~HeatTransferStabilizationHelper

(

)

Definition at line 69 of file heat_transfer_stab_helper.C.

  {
    return;
  }

Member Function Documentation

libMesh::Real GRINS::HeatTransferStabilizationHelper::compute_res_energy_steady	(	AssemblyContext &	context,
		unsigned int	qp,
		const libMesh::Real	rho,
		const libMesh::Real	Cp,
		const libMesh::Real	k
	)		const

Definition at line 74 of file heat_transfer_stab_helper.C.

References _flow_vars, _temp_vars, GRINS::PrimitiveTempFEVariables::T(), and GRINS::VelocityFEVariables::w().

  {
    libMesh::Gradient grad_T = context.fixed_interior_gradient(this->_temp_vars.T(), qp);
    libMesh::Tensor hess_T = context.fixed_interior_hessian(this->_temp_vars.T(), qp);

    libMesh::RealGradient rhocpU( rho*Cp*context.fixed_interior_value(this->_flow_vars.u(), qp), 
                                  rho*Cp*context.fixed_interior_value(this->_flow_vars.v(), qp) );
    if(context.get_system().get_mesh().mesh_dimension() == 3)
      rhocpU(2) = rho*Cp*context.fixed_interior_value(this->_flow_vars.w(), qp);

    return rhocpU*grad_T - k*(hess_T(0,0) + hess_T(1,1) + hess_T(2,2));
  }

void GRINS::HeatTransferStabilizationHelper::compute_res_energy_steady_and_derivs	(	AssemblyContext &	context,
		unsigned int	qp,
		const libMesh::Real	rho,
		const libMesh::Real	Cp,
		const libMesh::Real	k,
		libMesh::Real &	res,
		libMesh::Real &	d_res_dT,
		libMesh::Gradient &	d_res_dgradT,
		libMesh::Tensor &	d_res_dhessT,
		libMesh::Gradient &	d_res_dU
	)		const

Definition at line 92 of file heat_transfer_stab_helper.C.

  {
    libMesh::Gradient grad_T = context.fixed_interior_gradient(this->_temp_vars.T(), qp);
    libMesh::Tensor hess_T = context.fixed_interior_hessian(this->_temp_vars.T(), qp);

    libMesh::RealGradient rhocpU( rho*Cp*context.fixed_interior_value(this->_flow_vars.u(), qp), 
                                  rho*Cp*context.fixed_interior_value(this->_flow_vars.v(), qp) );
    if(context.get_system().get_mesh().mesh_dimension() == 3)
      rhocpU(2) = rho*Cp*context.fixed_interior_value(this->_flow_vars.w(), qp);

    res = rhocpU*grad_T - k*(hess_T(0,0) + hess_T(1,1) + hess_T(2,2));
    d_res_dT = 0;
    d_res_dgradT = rhocpU;
    d_res_dhessT = 0;
    d_res_dhessT(0,0) = -k;
    d_res_dhessT(1,1) = -k;
    d_res_dhessT(2,2) = -k;
    d_res_dU = rho * Cp * grad_T;
  }

libMesh::Real GRINS::HeatTransferStabilizationHelper::compute_res_energy_transient	(	AssemblyContext &	context,
		unsigned int	qp,
		const libMesh::Real	rho,
		const libMesh::Real	Cp
	)		const

Definition at line 123 of file heat_transfer_stab_helper.C.

References _temp_vars, and GRINS::PrimitiveTempFEVariables::T().

  {
    libMesh::Real T_dot;
    context.interior_rate(this->_temp_vars.T(), qp, T_dot);

    return rho*Cp*T_dot;
  }

void GRINS::HeatTransferStabilizationHelper::compute_res_energy_transient_and_derivs	(	AssemblyContext &	context,
		unsigned int	qp,
		const libMesh::Real	rho,
		const libMesh::Real	Cp,
		libMesh::Real &	res,
		libMesh::Real &	d_res_dTdot
	)		const

Definition at line 136 of file heat_transfer_stab_helper.C.

  {
    libMesh::Real T_dot;
    context.interior_rate(this->_temp_vars.T(), qp, T_dot);

    res = rho*Cp*T_dot;
    d_res_dTdot = rho*Cp;
  }

libMesh::Real GRINS::HeatTransferStabilizationHelper::compute_tau_energy ( AssemblyContext &  c, libMesh::RealTensor &  G, libMesh::Real  rho, libMesh::Real  cp, libMesh::Real  k, libMesh::Gradient  U, bool  is_steady  ) const

inline

Definition at line 117 of file heat_transfer_stab_helper.h.

References _C, and _tau_factor.

  {
    libMesh::Real tau = (rho*cp*U)*(G*(rho*cp*U)) + this->_C*k*k*G.contract(G);
    
    if(!is_steady)
      tau += (2.0*rho*cp/c.get_deltat_value())*(2.0*rho*cp/c.get_deltat_value());

    return this->_tau_factor/std::sqrt(tau);
  }

void GRINS::HeatTransferStabilizationHelper::compute_tau_energy_and_derivs ( AssemblyContext &  c, libMesh::RealTensor &  G, libMesh::Real  rho, libMesh::Real  cp, libMesh::Real  k, libMesh::Gradient  U, libMesh::Real &  tau_E, libMesh::Real &  d_tau_E_d_rho, libMesh::Gradient &  d_tau_E_d_U, bool  is_steady  ) const

inline

Definition at line 136 of file heat_transfer_stab_helper.h.

  {
    libMesh::Gradient rhocpU = rho*cp*U;
    libMesh::Gradient GrhocpU = G*rhocpU;
    libMesh::Real rhocpUGrhocpU = rhocpU * GrhocpU;
    libMesh::Real GG = G.contract(G);
    tau_E = (rhocpU)*(GrhocpU) + this->_C*k*k*GG;
    d_tau_E_d_rho = rhocpUGrhocpU*2/rho/cp;
    d_tau_E_d_U = 2*rho*cp*GrhocpU;
    
    if(!is_steady)
      {
        libMesh::Real two_rhocp_over_dt = 2*rho*cp/c.get_deltat_value();
        tau_E += two_rhocp_over_dt * two_rhocp_over_dt;
        d_tau_E_d_rho += 4*two_rhocp_over_dt/c.get_deltat_value();
      }

    // But what we've computed so far isn't tau; we need
    // tau = _tau_factor/ sqrt(our_tau)

    libMesh::Real root_oldtau = std::sqrt(tau_E);
    libMesh::Real d_tau_d_oldtau = -this->_tau_factor / (tau_E*root_oldtau) / 2;

    d_tau_E_d_rho         = d_tau_d_oldtau * d_tau_E_d_rho;
    d_tau_E_d_U           = d_tau_d_oldtau * d_tau_E_d_U;

    tau_E = this->_tau_factor / root_oldtau;
  }

Member Data Documentation

libMesh::Real GRINS::HeatTransferStabilizationHelper::_C

protected

Definition at line 106 of file heat_transfer_stab_helper.h.

Referenced by compute_tau_energy().

const VelocityFEVariables& GRINS::HeatTransferStabilizationHelper::_flow_vars

protected

Definition at line 110 of file heat_transfer_stab_helper.h.

Referenced by compute_res_energy_steady().

const PressureFEVariable& GRINS::HeatTransferStabilizationHelper::_press_var

protected

Definition at line 111 of file heat_transfer_stab_helper.h.

libMesh::Real GRINS::HeatTransferStabilizationHelper::_tau_factor

protected

Definition at line 106 of file heat_transfer_stab_helper.h.

Referenced by compute_tau_energy().

const PrimitiveTempFEVariables& GRINS::HeatTransferStabilizationHelper::_temp_vars

protected

Definition at line 108 of file heat_transfer_stab_helper.h.

Referenced by compute_res_energy_steady(), and compute_res_energy_transient().

---

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

• src/physics/include/grins/heat_transfer_stab_helper.h
• src/physics/src/heat_transfer_stab_helper.C