GRINS::SpalartAllmarasStabilizationHelper Class Reference

#include <spalart_allmaras_stab_helper.h>

Inheritance diagram for GRINS::SpalartAllmarasStabilizationHelper:

Collaboration diagram for GRINS::SpalartAllmarasStabilizationHelper:

Public Member Functions
	SpalartAllmarasStabilizationHelper (const std::string &helper_name, const GetPot &input)
	~SpalartAllmarasStabilizationHelper ()
libMesh::Real	compute_tau_spalart (AssemblyContext &c, unsigned int qp, libMesh::RealGradient &g, libMesh::RealTensor &G, libMesh::Real rho, libMesh::Gradient U, libMesh::Real mu, bool is_steady) const
void	compute_tau_spalart_and_derivs (AssemblyContext &c, unsigned int qp, libMesh::RealGradient &g, libMesh::RealTensor &G, libMesh::Real rho, libMesh::Gradient U, libMesh::Real T, libMesh::Real &tau_M, libMesh::Real &d_tau_M_d_rho, libMesh::Gradient &d_tau_M_d_U, bool is_steady) const
libMesh::Real	compute_tau (AssemblyContext &c, unsigned int qp, libMesh::Real mat_prop_sq, libMesh::RealGradient &g, libMesh::RealTensor &G, libMesh::Real rho, libMesh::Gradient U, bool is_steady) const
void	compute_tau_and_derivs (AssemblyContext &c, unsigned int qp, libMesh::Real mat_prop_sq, libMesh::RealGradient &g, libMesh::RealTensor &G, libMesh::Real rho, libMesh::Gradient U, libMesh::Real &tau, libMesh::Real &d_tau_d_rho, libMesh::Gradient &d_tau_d_U, bool is_steady) const
libMesh::Real	compute_res_spalart_steady (AssemblyContext &context, unsigned int qp, const libMesh::Real rho, const libMesh::Real mu, const libMesh::Real distance_qp, const bool infinite_distance) const
void	compute_res_spalart_steady_and_derivs (AssemblyContext &context, unsigned int qp, const libMesh::Real rho, const libMesh::Real mu, libMesh::Gradient &res_M, libMesh::Tensor &d_res_M_dgradp, libMesh::Tensor &d_res_M_dU, libMesh::Gradient &d_res_Muvw_dgraduvw, libMesh::Tensor &d_res_Muvw_dhessuvw) const
libMesh::Real	compute_res_spalart_transient (AssemblyContext &context, unsigned int qp, const libMesh::Real rho) const
void	compute_res_spalart_transient_and_derivs (AssemblyContext &context, unsigned int qp, const libMesh::Real rho, libMesh::RealGradient &res_M, libMesh::Real &d_res_Muvw_duvw) const
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...
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...

Protected Attributes
libMesh::Real	_C
libMesh::Real	_tau_factor
const VelocityVariable &	_flow_vars
const PressureFEVariable &	_press_var
const TurbulenceFEVariables &	_turbulence_vars
SpalartAllmarasHelper	_spalart_allmaras_helper
SpalartAllmarasParameters	_sa_params

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 46 of file spalart_allmaras_stab_helper.h.

Constructor & Destructor Documentation

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

Definition at line 43 of file spalart_allmaras_stab_helper.C.

References _C, _tau_factor, and GRINS::ParameterUser::set_parameter().

    : StabilizationHelper(helper_name),
      _C( input("Stabilization/tau_constant_vel", input("Stabilization/tau_constant", 1.0 ) ) ),
      _tau_factor( input("Stabilization/tau_factor_vel", input("Stabilization/tau_factor", 0.5 ) ) ),
      _flow_vars(GRINSPrivate::VariableWarehouse::get_variable_subclass<VelocityVariable>(VariablesParsing::velocity_variable_name(input,PhysicsNaming::spalart_allmaras(),VariablesParsing::PHYSICS))),
      _press_var(GRINSPrivate::VariableWarehouse::get_variable_subclass<PressureFEVariable>(VariablesParsing::press_variable_name(input,PhysicsNaming::spalart_allmaras(),VariablesParsing::PHYSICS))),
      _turbulence_vars(GRINSPrivate::VariableWarehouse::get_variable_subclass<TurbulenceFEVariables>(VariablesParsing::turb_variable_name(input,PhysicsNaming::spalart_allmaras(),VariablesParsing::PHYSICS))),
      _spalart_allmaras_helper(input),
      _sa_params(input)
  {
    this->set_parameter(this->_C ,input, "Stabilization/tau_constant_vel" , this->_C );
    this->set_parameter(this->_tau_factor ,input, "Stabilization/tau_factor_sa", this->_tau_factor );
  }

GRINS::SpalartAllmarasStabilizationHelper::~SpalartAllmarasStabilizationHelper ( )

inline

Definition at line 52 of file spalart_allmaras_stab_helper.h.

{};

Member Function Documentation

libMesh::Real GRINS::SpalartAllmarasStabilizationHelper::compute_res_spalart_steady	(	AssemblyContext &	context,
		unsigned int	qp,
		const libMesh::Real	rho,
		const libMesh::Real	mu,
		const libMesh::Real	distance_qp,
		const bool	infinite_distance
	)		const

Definition at line 66 of file spalart_allmaras_stab_helper.C.

References _flow_vars, _sa_params, _spalart_allmaras_helper, _turbulence_vars, GRINS::SpalartAllmarasParameters::destruction_fn(), GRINS::MultcomponentVectorVariable::dim(), GRINS::SpalartAllmarasParameters::get_cb1(), GRINS::SpalartAllmarasParameters::get_cb2(), GRINS::SpalartAllmarasParameters::get_kappa(), GRINS::SpalartAllmarasParameters::get_sigma(), GRINS::TurbulenceFEVariables::nu(), GRINS::SpalartAllmarasParameters::source_fn(), GRINS::MultcomponentVectorVariable::u(), GRINS::MultcomponentVectorVariable::v(), GRINS::SpalartAllmarasHelper::vorticity(), and GRINS::MultcomponentVectorVariable::w().

  {
    // The flow velocity
    libMesh::Number u,v;
    u = context.interior_value(this->_flow_vars.u(), qp);
    v = context.interior_value(this->_flow_vars.v(), qp);

    libMesh::NumberVectorValue U(u,v);
    if ( this->_flow_vars.dim() == 3 )
      U(2) = context.interior_value(this->_flow_vars.w(), qp);

    libMesh::RealGradient grad_u = context.fixed_interior_gradient(this->_flow_vars.u(), qp);
    libMesh::RealGradient grad_v = context.fixed_interior_gradient(this->_flow_vars.v(), qp);

    libMesh::Number nu_value = context.interior_value(this->_turbulence_vars.nu(), qp);

    libMesh::RealGradient grad_nu = context.fixed_interior_gradient(this->_turbulence_vars.nu(), qp);

    libMesh::RealTensor hess_nu = context.fixed_interior_hessian(this->_turbulence_vars.nu(), qp);

    // The convection term
    libMesh::Number rhoUdotGradnu = rho*(U*grad_nu);

    // The diffusion term
    libMesh::Number inv_sigmadivnuplusnuphysicalGradnu = (1./this->_sa_params.get_sigma())*(grad_nu*grad_nu + ((nu_value + mu)*(hess_nu(0,0) + hess_nu(1,1) + (this->_flow_vars.dim() == 3)?hess_nu(2,2):0)) + this->_sa_params.get_cb2()*grad_nu*grad_nu);

    // The source term
    libMesh::Real vorticity_value_qp = this->_spalart_allmaras_helper.vorticity(context, qp);
    libMesh::Real S_tilde = this->_sa_params.source_fn(nu_value, mu, distance_qp, vorticity_value_qp, infinite_distance);
    libMesh::Real source_term = this->_sa_params.get_cb1()*S_tilde*nu_value;

    libMesh::Real kappa2 = (this->_sa_params.get_kappa())*(this->_sa_params.get_kappa());
    libMesh::Real cw1 = this->_sa_params.get_cb1()/kappa2 + (1.0 + this->_sa_params.get_cb2())/this->_sa_params.get_sigma();

    // The destruction term
    libMesh::Real fw = this->_sa_params.destruction_fn(nu_value, distance_qp, S_tilde, infinite_distance);
    libMesh::Real destruction_term = 0.0;
    if(infinite_distance)
      {
        destruction_term = 0.0;
      }
    else
      {
        destruction_term =  cw1*fw*pow(nu_value/distance_qp, 2.);
      }

    return rhoUdotGradnu + source_term + inv_sigmadivnuplusnuphysicalGradnu - destruction_term;
  }

void GRINS::SpalartAllmarasStabilizationHelper::compute_res_spalart_steady_and_derivs	(	AssemblyContext &	context,
		unsigned int	qp,
		const libMesh::Real	rho,
		const libMesh::Real	mu,
		libMesh::Gradient &	res_M,
		libMesh::Tensor &	d_res_M_dgradp,
		libMesh::Tensor &	d_res_M_dU,
		libMesh::Gradient &	d_res_Muvw_dgraduvw,
		libMesh::Tensor &	d_res_Muvw_dhessuvw
	)		const

Definition at line 117 of file spalart_allmaras_stab_helper.C.

  {
    // To be filled when we start using analytic jacobians with SA
    libmesh_not_implemented();
  }

libMesh::Real GRINS::SpalartAllmarasStabilizationHelper::compute_res_spalart_transient	(	AssemblyContext &	context,
		unsigned int	qp,
		const libMesh::Real	rho
	)		const

Definition at line 131 of file spalart_allmaras_stab_helper.C.

References _turbulence_vars, and GRINS::TurbulenceFEVariables::nu().

  {
    libMesh::Number nu_dot = context.interior_value(this->_turbulence_vars.nu(), qp);

    return rho*nu_dot;
  }

void GRINS::SpalartAllmarasStabilizationHelper::compute_res_spalart_transient_and_derivs	(	AssemblyContext &	context,
		unsigned int	qp,
		const libMesh::Real	rho,
		libMesh::RealGradient &	res_M,
		libMesh::Real &	d_res_Muvw_duvw
	)		const

Definition at line 140 of file spalart_allmaras_stab_helper.C.

  {
    libmesh_not_implemented();
  }

libMesh::Real GRINS::SpalartAllmarasStabilizationHelper::compute_tau ( AssemblyContext &  c, unsigned int  qp, libMesh::Real  mat_prop_sq, libMesh::RealGradient &  g, libMesh::RealTensor &  G, libMesh::Real  rho, libMesh::Gradient  U, bool  is_steady  ) const

inline

Definition at line 184 of file spalart_allmaras_stab_helper.h.

References _C, and _tau_factor.

Referenced by compute_tau_spalart().

  {
    libMesh::Real tau = (rho*U)*(G*(rho*U)) + this->_C*mat_prop_sq*G.contract(G);

    if(!is_steady)
      tau += (2.0*rho/c.get_deltat_value())*(2.0*rho/c.get_deltat_value());

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

void GRINS::SpalartAllmarasStabilizationHelper::compute_tau_and_derivs ( AssemblyContext &  c, unsigned int  qp, libMesh::Real  mat_prop_sq, libMesh::RealGradient &  g, libMesh::RealTensor &  G, libMesh::Real  rho, libMesh::Gradient  U, libMesh::Real &  tau, libMesh::Real &  d_tau_d_rho, libMesh::Gradient &  d_tau_d_U, bool  is_steady  ) const

inline

Definition at line 204 of file spalart_allmaras_stab_helper.h.

  {
    libMesh::Gradient rhoU = rho*U;
    libMesh::Gradient GrhoU = G*rhoU;
    libMesh::Real rhoUGrhoU = rhoU * GrhoU;
    libMesh::Real GG = G.contract(G);
    tau = rhoUGrhoU + this->_C*mat_prop_sq*GG;
    d_tau_d_rho = rhoUGrhoU*2/rho;
    d_tau_d_U = 2*rho*GrhoU;

    if(!is_steady)
      {
        libMesh::Real two_rho_over_dt = 2*rho/c.get_deltat_value();
        tau += two_rho_over_dt * two_rho_over_dt;
        d_tau_d_rho += 4*two_rho_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);
    libMesh::Real d_tau_d_oldtau = -this->_tau_factor / (tau*root_oldtau) / 2;

    d_tau_d_rho         = d_tau_d_oldtau * d_tau_d_rho;
    d_tau_d_U           = d_tau_d_oldtau * d_tau_d_U;

    tau = this->_tau_factor / root_oldtau;
  }

libMesh::Real GRINS::SpalartAllmarasStabilizationHelper::compute_tau_spalart ( AssemblyContext &  c, unsigned int  qp, libMesh::RealGradient &  g, libMesh::RealTensor &  G, libMesh::Real  rho, libMesh::Gradient  U, libMesh::Real  mu, bool  is_steady  ) const

inline

Definition at line 151 of file spalart_allmaras_stab_helper.h.

References compute_tau().

  {
    return this->compute_tau( c, qp, mu*mu, g, G, rho, U, is_steady );
  }

void GRINS::SpalartAllmarasStabilizationHelper::compute_tau_spalart_and_derivs ( AssemblyContext &  c, unsigned int  qp, libMesh::RealGradient &  g, libMesh::RealTensor &  G, libMesh::Real  rho, libMesh::Gradient  U, libMesh::Real  T, libMesh::Real &  tau_M, libMesh::Real &  d_tau_M_d_rho, libMesh::Gradient &  d_tau_M_d_U, bool  is_steady  ) const

inline

Definition at line 165 of file spalart_allmaras_stab_helper.h.

  {
    this->compute_tau_and_derivs( c, qp, mu*mu, g, G, rho, U, tau_M,
                                  d_tau_M_d_rho, d_tau_M_d_U,
                                  is_steady );
  }

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

virtual

Each subclass will register its copy of an independent.

Reimplemented from GRINS::ParameterUser.

Definition at line 59 of file spalart_allmaras_stab_helper.C.

References GRINS::ParameterUser::register_parameter().

  {
    ParameterUser::register_parameter(param_name, param_pointer);
    this->_sa_params.register_parameter(param_name, param_pointer);
  }

Member Data Documentation

libMesh::Real GRINS::SpalartAllmarasStabilizationHelper::_C

protected

Definition at line 135 of file spalart_allmaras_stab_helper.h.

Referenced by compute_tau(), and SpalartAllmarasStabilizationHelper().

const VelocityVariable& GRINS::SpalartAllmarasStabilizationHelper::_flow_vars

protected

Definition at line 137 of file spalart_allmaras_stab_helper.h.

Referenced by compute_res_spalart_steady().

const PressureFEVariable& GRINS::SpalartAllmarasStabilizationHelper::_press_var

protected

Definition at line 138 of file spalart_allmaras_stab_helper.h.

SpalartAllmarasParameters GRINS::SpalartAllmarasStabilizationHelper::_sa_params

protected

Definition at line 144 of file spalart_allmaras_stab_helper.h.

Referenced by compute_res_spalart_steady().

SpalartAllmarasHelper GRINS::SpalartAllmarasStabilizationHelper::_spalart_allmaras_helper

protected

Definition at line 142 of file spalart_allmaras_stab_helper.h.

Referenced by compute_res_spalart_steady().

libMesh::Real GRINS::SpalartAllmarasStabilizationHelper::_tau_factor

protected

Definition at line 135 of file spalart_allmaras_stab_helper.h.

Referenced by compute_tau(), and SpalartAllmarasStabilizationHelper().

const TurbulenceFEVariables& GRINS::SpalartAllmarasStabilizationHelper::_turbulence_vars

protected

Definition at line 140 of file spalart_allmaras_stab_helper.h.

Referenced by compute_res_spalart_steady(), and compute_res_spalart_transient().

---

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

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