#include <spalart_allmaras_viscosity.h>

Inheritance diagram for GRINS::SpalartAllmarasViscosity< Viscosity >:

Collaboration diagram for GRINS::SpalartAllmarasViscosity< Viscosity >:

Public Member Functions
	SpalartAllmarasViscosity (const GetPot &input)

	~SpalartAllmarasViscosity ()

libMesh::Real	operator() (AssemblyContext &context, unsigned int qp) const

libMesh::Real	operator() (const libMesh::Point &p, const libMesh::Real time=0)

void	init (libMesh::FEMSystem *system)

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
Viscosity	_mu
	Viscosity object (so we have access to the physical viscosity) More...

TurbulenceFEVariables	_turbulence_vars

SpalartAllmarasParameters	_sa_params

Private Member Functions
	SpalartAllmarasViscosity ()

Detailed Description

template<class Viscosity>
class GRINS::SpalartAllmarasViscosity< Viscosity >

Definition at line 52 of file spalart_allmaras_viscosity.h.

Constructor & Destructor Documentation

template<class Mu >

GRINS::SpalartAllmarasViscosity< Mu >::SpalartAllmarasViscosity ( const GetPot & input )

Definition at line 32 of file spalart_allmaras_viscosity.C.

                                                                              :
     ParameterUser("SpalartAllmarasViscosity"),
     _mu(input),
     _turbulence_vars(input, spalart_allmaras),
     _sa_params(input)
   {
     if( !input.have_variable("Materials/Viscosity/mu") )
       {
         std::cerr<<"No viscosity has been specified."<<std::endl;
 
         libmesh_error();
       }
     return;
   }

template<class Viscosity >

GRINS::SpalartAllmarasViscosity< Viscosity >::~SpalartAllmarasViscosity ( )

inline

Definition at line 57 of file spalart_allmaras_viscosity.h.

57 {};

template<class Viscosity >

GRINS::SpalartAllmarasViscosity< Viscosity >::SpalartAllmarasViscosity ( )

private

Member Function Documentation

template<class Mu >

void GRINS::SpalartAllmarasViscosity< Mu >::init ( libMesh::FEMSystem * system )

Definition at line 48 of file spalart_allmaras_viscosity.C.

   {
     this->_turbulence_vars.init(system);
   }

template<class Mu >

libMesh::Real GRINS::SpalartAllmarasViscosity< Mu >::operator()	(	AssemblyContext &	context,
		unsigned int	qp
	)		const

Definition at line 54 of file spalart_allmaras_viscosity.C.

   {
     // The physical viscosity
     libMesh::Real mu_physical = this->_mu(context, qp);
 
     // The unscaled turbulent viscosity (the nu the SA physics solves for)
     libMesh::Real nu = context.interior_value(this->_turbulence_vars.nu_var(),qp);
 
     // Assert that mu_value is greater than 0
     if(nu < 0.0)
       {
         libmesh_warning("Negative turbulent viscosity encountered !");
 
         // We are using a negative S-A model, so will set eddy viscosity to zero
         // if the turbulent viscosity nu < 0.0
         nu = 0.0;
       }
 
     // Step 1
     libMesh::Real chi = nu/mu_physical;
 
     // Step 2
     libMesh::Real fv1 = _sa_params.fv1(chi);
 
     // Step 3
     libMesh::Real mu_turbulent = nu*fv1;
 
     // Compute the value of the total viscosity and return it
     libMesh::Number mu_value = mu_turbulent + mu_physical; // Turbulent viscosity + physical viscosity
 
     return mu_value;
   }

template<class Viscosity >

libMesh::Real GRINS::SpalartAllmarasViscosity< Viscosity >::operator()	(	const libMesh::Point &	p,
		const libMesh::Real	time = `0`
	)

inline

Definition at line 61 of file spalart_allmaras_viscosity.h.

References GRINS::SpalartAllmarasViscosity< Viscosity >::_mu.

62 { return _mu(p,time); }

GRINS::SpalartAllmarasViscosity::_mu

Viscosity _mu

Viscosity object (so we have access to the physical viscosity)

Definition: spalart_allmaras_viscosity.h:69

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::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;
   }