#include <spalart_allmaras_stab_helper.h>

Inheritance diagram for GRINS::SpalartAllmarasStabilizationHelper:

Collaboration diagram for GRINS::SpalartAllmarasStabilizationHelper:

[legend]

Public Member Functions
	SpalartAllmarasStabilizationHelper (const std::string &helper_name, const GetPot &input)

	~SpalartAllmarasStabilizationHelper ()

void	init (libMesh::FEMSystem &system)

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
unsigned int	_dim
	Physical dimension of problem. More...

libMesh::Real	_C

libMesh::Real	_tau_factor

const VelocityFEVariables &	_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 44 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<VelocityFEVariables>(VariablesParsing::velocity_section())),
       _press_var(GRINSPrivate::VariableWarehouse::get_variable_subclass<PressureFEVariable>(VariablesParsing::pressure_section())),
       _turbulence_vars(GRINSPrivate::VariableWarehouse::get_variable_subclass<TurbulenceFEVariables>(VariablesParsing::turbulence_section())),
       _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 ( )

Definition at line 67 of file spalart_allmaras_stab_helper.C.

   {
     return;
   }

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 80 of file spalart_allmaras_stab_helper.C.

References _dim, _flow_vars, _sa_params, _spalart_allmaras_helper, _turbulence_vars, GRINS::SpalartAllmarasParameters::destruction_fn(), GRINS::SpalartAllmarasParameters::get_cb1(), GRINS::SpalartAllmarasParameters::get_cb2(), GRINS::SpalartAllmarasParameters::get_kappa(), GRINS::SpalartAllmarasParameters::get_sigma(), GRINS::TurbulenceFEVariables::nu(), GRINS::SpalartAllmarasParameters::source_fn(), GRINS::VelocityFEVariables::u(), GRINS::VelocityFEVariables::v(), GRINS::SpalartAllmarasHelper::vorticity(), and GRINS::VelocityFEVariables::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 ( context.get_system().get_mesh().mesh_dimension() == 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->_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 131 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 145 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 154 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 191 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 211 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 158 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 172 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::init ( libMesh::FEMSystem & system )

Definition at line 72 of file spalart_allmaras_stab_helper.C.

References _dim, _spalart_allmaras_helper, and GRINS::SpalartAllmarasHelper::init_variables().

   {
     // Init the variables belonging to SA helper
     _spalart_allmaras_helper.init_variables(&system);
 
     this->_dim = system.get_mesh().mesh_dimension();
   }

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

Referenced by compute_tau(), and SpalartAllmarasStabilizationHelper().

unsigned int GRINS::SpalartAllmarasStabilizationHelper::_dim

protected

Physical dimension of problem.

Todo:: Do we really need to cache this?

Definition at line 140 of file spalart_allmaras_stab_helper.h.

Referenced by compute_res_spalart_steady(), and init().

const VelocityFEVariables& GRINS::SpalartAllmarasStabilizationHelper::_flow_vars

protected

Definition at line 144 of file spalart_allmaras_stab_helper.h.

Referenced by compute_res_spalart_steady().

const PressureFEVariable& GRINS::SpalartAllmarasStabilizationHelper::_press_var

protected

Definition at line 145 of file spalart_allmaras_stab_helper.h.

SpalartAllmarasParameters GRINS::SpalartAllmarasStabilizationHelper::_sa_params

protected

Definition at line 151 of file spalart_allmaras_stab_helper.h.

Referenced by compute_res_spalart_steady().

SpalartAllmarasHelper GRINS::SpalartAllmarasStabilizationHelper::_spalart_allmaras_helper

protected

Definition at line 149 of file spalart_allmaras_stab_helper.h.

Referenced by compute_res_spalart_steady(), and init().

libMesh::Real GRINS::SpalartAllmarasStabilizationHelper::_tau_factor

protected

Definition at line 142 of file spalart_allmaras_stab_helper.h.

Referenced by compute_tau(), and SpalartAllmarasStabilizationHelper().

const TurbulenceFEVariables& GRINS::SpalartAllmarasStabilizationHelper::_turbulence_vars

protected

Definition at line 147 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

Public Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation