#include <gas_recombination_catalytic_wall.h>

Inheritance diagram for GRINS::GasRecombinationCatalyticWall< Chemistry >:

Collaboration diagram for GRINS::GasRecombinationCatalyticWall< Chemistry >:

Public Member Functions
	GasRecombinationCatalyticWall (const Chemistry &chem_mixture, CatalycityBase &gamma, const unsigned int reactant_species_idx, const unsigned int product_species_idx)

virtual	~GasRecombinationCatalyticWall ()

virtual void	init (const libMesh::FEMSystem &system)

virtual void	apply_fluxes (AssemblyContext &context, const CachedValues &cache, const bool request_jacobian)

libMesh::Real	compute_reactant_mass_flux (const libMesh::Real rho, const libMesh::Real Y_r, const libMesh::Real T)

libMesh::Real	compute_product_mass_flux (const libMesh::Real rho, const libMesh::Real Y_r, const libMesh::Real T)

void	set_axisymmetric (bool is_axisymmetric)

libMesh::Real	omega_dot (const libMesh::Real rho_s, const libMesh::Real T) const
	$\rho_s \gamma \sqrt{ \frac{R_s T}{2\pi} }$ More...

libMesh::Real	domega_dot_dws (const libMesh::Real rho_s, const libMesh::Real w_s, const libMesh::Real T, const libMesh::Real R) const

libMesh::Real	domega_dot_dT (const libMesh::Real rho_s, const libMesh::Real T) const

void	set_catalycity_params (const std::vector< libMesh::Real > &params)

Protected Attributes
const unsigned int	_reactant_species_idx

VariableIndex	_reactant_var_idx

const unsigned int	_product_species_idx

VariableIndex	_product_var_idx

const Chemistry &	_chemistry

boost::scoped_ptr< CatalycityBase >	_gamma_s

const libMesh::Real	_C
	$\sqrt{ \frac{R_s}{2\pi} }$ More...

bool	_is_axisymmetric

Detailed Description

template<typename Chemistry>
class GRINS::GasRecombinationCatalyticWall< Chemistry >

Definition at line 36 of file gas_recombination_catalytic_wall.h.

Constructor & Destructor Documentation

template<typename Chemistry >

GRINS::GasRecombinationCatalyticWall< Chemistry >::GasRecombinationCatalyticWall	(	const Chemistry &	chem_mixture,
		CatalycityBase &	gamma,
		const unsigned int	reactant_species_idx,
		const unsigned int	product_species_idx
	)

Definition at line 34 of file gas_recombination_catalytic_wall.C.

     : CatalyticWallBase<Chemistry>(chem_mixture,gamma,reactant_species_idx),
       _reactant_species_idx(reactant_species_idx),
       _product_species_idx(product_species_idx)
   {
     return;
   }

template<typename Chemistry >

GRINS::GasRecombinationCatalyticWall< Chemistry >::~GasRecombinationCatalyticWall ( )

virtual

Definition at line 46 of file gas_recombination_catalytic_wall.C.

   {
     return;
   }

Member Function Documentation

template<typename Chemistry >

void GRINS::GasRecombinationCatalyticWall< Chemistry >::apply_fluxes	(	AssemblyContext &	context,
		const CachedValues &	cache,
		const bool	request_jacobian
	)

virtual

Implements GRINS::CatalyticWallBase< Chemistry >.

Definition at line 69 of file gas_recombination_catalytic_wall.C.

References GRINS::CachedValues::get_cached_values(), GRINS::CachedValues::get_cached_vector_values(), GRINS::Cache::MASS_FRACTIONS, GRINS::Cache::MIXTURE_DENSITY, and GRINS::Cache::TEMPERATURE.

   {
     libMesh::FEGenericBase<libMesh::Real>* side_fe = NULL; 
     context.get_side_fe( _reactant_var_idx, side_fe );
 
     // The number of local degrees of freedom in each variable.
     const unsigned int n_var_dofs = context.get_dof_indices(_reactant_var_idx).size();
 
     libmesh_assert_equal_to( n_var_dofs, context.get_dof_indices(_product_var_idx).size() );
 
     // Element Jacobian * quadrature weight for side integration.
     const std::vector<libMesh::Real> &JxW_side = side_fe->get_JxW();
 
     // The var shape functions at side quadrature points.
     const std::vector<std::vector<libMesh::Real> >& var_phi_side = side_fe->get_phi();
 
     // Physical location of the quadrature points
     const std::vector<libMesh::Point>& var_qpoint = side_fe->get_xyz();
 
     // reactant residual
     libMesh::DenseSubVector<libMesh::Number> &F_r_var = context.get_elem_residual(_reactant_var_idx); 
 
     // product residual
     libMesh::DenseSubVector<libMesh::Number> &F_p_var = context.get_elem_residual(_product_var_idx);
 
     unsigned int n_qpoints = context.get_side_qrule().n_points();
 
     for (unsigned int qp=0; qp != n_qpoints; qp++)
       {
         libMesh::Real jac = JxW_side[qp];
 
         if( this->_is_axisymmetric )
           {
             const libMesh::Number r = var_qpoint[qp](0);
             jac *= r;
           }
 
         const libMesh::Real rho = cache.get_cached_values(Cache::MIXTURE_DENSITY)[qp];
 
         const libMesh::Real Y_r = cache.get_cached_vector_values(Cache::MASS_FRACTIONS)[qp][this->_reactant_species_idx];
 
         const libMesh::Real T = cache.get_cached_values(Cache::TEMPERATURE)[qp];
 
         const libMesh::Real r_value = this->compute_reactant_mass_flux(rho, Y_r, T);
 
         const libMesh::Real p_value = -r_value;
 
         for (unsigned int i=0; i != n_var_dofs; i++)
           {
             F_r_var(i) += r_value*var_phi_side[i][qp]*jac;
 
             F_p_var(i) += p_value*var_phi_side[i][qp]*jac;
 
             if( request_jacobian )
               {
                 libmesh_not_implemented();
               }
           }
       }
 
     return;
   }

template<typename Chemistry >

libMesh::Real GRINS::GasRecombinationCatalyticWall< Chemistry >::compute_product_mass_flux	(	const libMesh::Real	rho,
		const libMesh::Real	Y_r,
		const libMesh::Real	T
	)

inline

Definition at line 86 of file gas_recombination_catalytic_wall.h.

   {
     return -(this->compute_reactant_mass_flux(rho,Y_r,T));
   }

template<typename Chemistry >

libMesh::Real GRINS::GasRecombinationCatalyticWall< Chemistry >::compute_reactant_mass_flux	(	const libMesh::Real	rho,
		const libMesh::Real	Y_r,
		const libMesh::Real	T
	)

inline

Definition at line 73 of file gas_recombination_catalytic_wall.h.

   {  
     const libMesh::Real rho_r = rho*Y_r;
     
     const libMesh::Real omega_dot = this->omega_dot( rho_r, T );
 
     return -omega_dot;
   }

template<typename Chemistry >

libMesh::Real GRINS::CatalyticWallBase< Chemistry >::domega_dot_dT	(	const libMesh::Real	rho_s,
		const libMesh::Real	T
	)		const

inlineinherited

Definition at line 108 of file catalytic_wall_base.h.

Referenced by test().

   {
     libMesh::Real sqrtT = std::sqrt(T);
 
     return rho_s*_C*( 0.5/sqrtT*(*_gamma_s)(T) + sqrtT*(*_gamma_s).dT(T) );
   }

template<typename Chemistry >

libMesh::Real GRINS::CatalyticWallBase< Chemistry >::domega_dot_dws	(	const libMesh::Real	rho_s,
		const libMesh::Real	w_s,
		const libMesh::Real	T,
		const libMesh::Real	R
	)		const

inlineinherited

Definition at line 100 of file catalytic_wall_base.h.

Referenced by test().

   {
     return (1.0/w_s - rho_s/R)*(this->omega_dot( rho_s, T ));
   }

template<typename Chemistry >

void GRINS::GasRecombinationCatalyticWall< Chemistry >::init ( const libMesh::FEMSystem & system )

virtual

Reimplemented from GRINS::CatalyticWallBase< Chemistry >.

Definition at line 52 of file gas_recombination_catalytic_wall.C.

   {
     const std::string r_var_name = std::string("w_"+this->_chemistry.species_name( this->_reactant_species_idx ) );
 
     const std::string p_var_name = std::string("w_"+this->_chemistry.species_name( this->_product_species_idx ) );
 
     libmesh_assert( system.has_variable( r_var_name ) );
     libmesh_assert( system.has_variable( p_var_name ) );
 
     this->_reactant_var_idx = system.variable_number( r_var_name );
 
     this->_product_var_idx = system.variable_number( p_var_name );
 
     return;
   }

template<typename Chemistry >

libMesh::Real GRINS::CatalyticWallBase< Chemistry >::omega_dot	(	const libMesh::Real	rho_s,
		const libMesh::Real	T
	)		const

inlineinherited

$\rho_s \gamma \sqrt{ \frac{R_s T}{2\pi} }$

Definition at line 93 of file catalytic_wall_base.h.

Referenced by test().

   {
     return rho_s*(*_gamma_s)(T)*_C*std::sqrt(T);
   }

template<typename Chemistry >

void GRINS::CatalyticWallBase< Chemistry >::set_axisymmetric ( bool is_axisymmetric )

inherited

Definition at line 64 of file catalytic_wall_base.C.

   {
     this->_is_axisymmetric = is_axisymmetric;
 
     return;
   }

template<typename Chemistry >

void GRINS::CatalyticWallBase< Chemistry >::set_catalycity_params ( const std::vector< libMesh::Real > & params )

inherited

Definition at line 73 of file catalytic_wall_base.C.

   {
     _gamma_s->set_params( params );
     return;
   }

Member Data Documentation

template<typename Chemistry>

const libMesh::Real GRINS::CatalyticWallBase< Chemistry >::_C

protectedinherited

$\sqrt{ \frac{R_s}{2\pi} }$

Definition at line 84 of file catalytic_wall_base.h.

template<typename Chemistry>

const Chemistry& GRINS::CatalyticWallBase< Chemistry >::_chemistry

protectedinherited

Definition at line 79 of file catalytic_wall_base.h.

template<typename Chemistry>

boost::scoped_ptr<CatalycityBase> GRINS::CatalyticWallBase< Chemistry >::_gamma_s

protectedinherited

Definition at line 81 of file catalytic_wall_base.h.

template<typename Chemistry>

bool GRINS::CatalyticWallBase< Chemistry >::_is_axisymmetric

protectedinherited

Definition at line 86 of file catalytic_wall_base.h.

template<typename Chemistry>

const unsigned int GRINS::GasRecombinationCatalyticWall< Chemistry >::_product_species_idx

protected

Definition at line 66 of file gas_recombination_catalytic_wall.h.

template<typename Chemistry>

VariableIndex GRINS::GasRecombinationCatalyticWall< Chemistry >::_product_var_idx

protected

Definition at line 67 of file gas_recombination_catalytic_wall.h.

template<typename Chemistry>

const unsigned int GRINS::GasRecombinationCatalyticWall< Chemistry >::_reactant_species_idx

protected

Definition at line 63 of file gas_recombination_catalytic_wall.h.

template<typename Chemistry>

VariableIndex GRINS::GasRecombinationCatalyticWall< Chemistry >::_reactant_var_idx

protected

Definition at line 64 of file gas_recombination_catalytic_wall.h.

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

src/boundary_conditions/include/grins/gas_recombination_catalytic_wall.h
src/boundary_conditions/src/gas_recombination_catalytic_wall.C

Public Member Functions

Protected Attributes

Detailed Description

template<typename Chemistry> class GRINS::GasRecombinationCatalyticWall< Chemistry >

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation

template<typename Chemistry>
class GRINS::GasRecombinationCatalyticWall< Chemistry >