#include <gas_solid_catalytic_wall.h>

Inheritance diagram for GRINS::GasSolidCatalyticWall< Chemistry >:

Collaboration diagram for GRINS::GasSolidCatalyticWall< Chemistry >:

Public Member Functions
	GasSolidCatalyticWall (const Chemistry &chem_mixture, CatalycityBase &gamma, const unsigned int reactant_gas_species_idx, const unsigned int reactant_solid_species_idx, const unsigned int product_species_idx)

virtual	~GasSolidCatalyticWall ()

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

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

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

libMesh::Real	compute_reactant_solid_mass_consumption (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)

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

libMesh::Real	compute_reactant_solid_mass_consumption_dYs (const libMesh::Real rho, const std::vector< libMesh::Real > Y, 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_gas_species_idx

VariableIndex	_reactant_gas_var_idx

const unsigned int	_reactant_solid_species_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::GasSolidCatalyticWall< Chemistry >

Definition at line 36 of file gas_solid_catalytic_wall.h.

Constructor & Destructor Documentation

template<typename Chemistry >

GRINS::GasSolidCatalyticWall< Chemistry >::GasSolidCatalyticWall	(	const Chemistry &	chem_mixture,
		CatalycityBase &	gamma,
		const unsigned int	reactant_gas_species_idx,
		const unsigned int	reactant_solid_species_idx,
		const unsigned int	product_species_idx
	)

Definition at line 34 of file gas_solid_catalytic_wall.C.

     : CatalyticWallBase<Chemistry>(chem_mixture,gamma,reactant_gas_species_idx),
       _reactant_gas_species_idx(reactant_gas_species_idx),
       _reactant_solid_species_idx(reactant_solid_species_idx),
       _product_species_idx(product_species_idx)
   {
     return;
   }

template<typename Chemistry >

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

virtual

Definition at line 48 of file gas_solid_catalytic_wall.C.

   {
     return;
   }

Member Function Documentation

template<typename Chemistry >

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

virtual

Implements GRINS::CatalyticWallBase< Chemistry >.

Definition at line 71 of file gas_solid_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_gas_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_gas_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_gas_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_gas_species_idx];
 
         const libMesh::Real T = cache.get_cached_values(Cache::TEMPERATURE)[qp];
 
         const libMesh::Real r_value = this->compute_reactant_gas_mass_flux(rho, Y_r, T);
 
         const libMesh::Real p_value = this->compute_product_mass_flux(rho, Y_r, T);
 
         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::GasSolidCatalyticWall< Chemistry >::compute_product_mass_flux	(	const libMesh::Real	rho,
		const libMesh::Real	Y_r,
		const libMesh::Real	T
	)

inline

Definition at line 119 of file gas_solid_catalytic_wall.h.

Referenced by test().

   {
     const libMesh::Real rho_r = rho*Y_r;
 
     const libMesh::Real M_r = this->_chemistry.M(_reactant_gas_species_idx);
 
     const libMesh::Real M_p = this->_chemistry.M(_product_species_idx);
 
     const libMesh::Real omega_dot = this->omega_dot( rho_r, T )*M_p/M_r;
 
     return omega_dot;
   }

template<typename Chemistry >

libMesh::Real GRINS::GasSolidCatalyticWall< Chemistry >::compute_reactant_gas_mass_flux	(	const libMesh::Real	rho,
		const libMesh::Real	Y_r,
		const libMesh::Real	T
	)

inline

Definition at line 89 of file gas_solid_catalytic_wall.h.

Referenced by test().

   {
     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::GasSolidCatalyticWall< Chemistry >::compute_reactant_solid_mass_consumption	(	const libMesh::Real	rho,
		const libMesh::Real	Y_r,
		const libMesh::Real	T
	)

inline

Definition at line 102 of file gas_solid_catalytic_wall.h.

Referenced by test().

   {
     const libMesh::Real rho_r = rho*Y_r;
 
     const libMesh::Real M_r = this->_chemistry.M(_reactant_gas_species_idx);
 
     const libMesh::Real M_solid = this->_chemistry.M(_reactant_solid_species_idx);
 
     const libMesh::Real omega_dot = this->omega_dot( rho_r, T )*M_solid/M_r;
 
     return -omega_dot;
   }

template<typename Chemistry >

libMesh::Real GRINS::GasSolidCatalyticWall< Chemistry >::compute_reactant_solid_mass_consumption_dT	(	const libMesh::Real	rho,
		const libMesh::Real	Y_r,
		const libMesh::Real	T
	)

inline

Definition at line 136 of file gas_solid_catalytic_wall.h.

   {
     const libMesh::Real rho_r = rho*Y_r;
 
     const libMesh::Real M_r = this->_chemistry.M(_reactant_gas_species_idx);
 
     const libMesh::Real M_solid = this->_chemistry.M(_reactant_solid_species_idx);
 
     const libMesh::Real domega_dot_dT = this->domega_dot_dT( rho_r, T )*M_solid/M_r;
 
     return -domega_dot_dT;
   }

template<typename Chemistry >

libMesh::Real GRINS::GasSolidCatalyticWall< Chemistry >::compute_reactant_solid_mass_consumption_dYs	(	const libMesh::Real	rho,
		const std::vector< libMesh::Real >	Y,
		const libMesh::Real	T
	)

inline

Definition at line 153 of file gas_solid_catalytic_wall.h.

   {
     const libMesh::Real Y_r = Y[_reactant_gas_species_idx];
 
     const libMesh::Real rho_r = rho*Y_r;
 
     const libMesh::Real M_r = this->_chemistry.M(_reactant_gas_species_idx);
 
     const libMesh::Real M_solid = this->_chemistry.M(_reactant_solid_species_idx);
 
     const libMesh::Real R = this->_chemistry.R_mix(Y);
 
     const libMesh::Real domega_dot_dYs = this->domega_dot_dws( rho_r, Y_r, T, R )*M_solid/M_r;
 
     return -domega_dot_dYs;
   }

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::GasSolidCatalyticWall< Chemistry >::init ( const libMesh::FEMSystem & system )

virtual

Reimplemented from GRINS::CatalyticWallBase< Chemistry >.

Definition at line 54 of file gas_solid_catalytic_wall.C.

   {
     const std::string r_var_name = std::string("w_"+this->_chemistry.species_name( this->_reactant_gas_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_gas_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::GasSolidCatalyticWall< Chemistry >::_product_species_idx

protected

Definition at line 82 of file gas_solid_catalytic_wall.h.

template<typename Chemistry>

VariableIndex GRINS::GasSolidCatalyticWall< Chemistry >::_product_var_idx

protected

Definition at line 83 of file gas_solid_catalytic_wall.h.

template<typename Chemistry>

const unsigned int GRINS::GasSolidCatalyticWall< Chemistry >::_reactant_gas_species_idx

protected

Definition at line 77 of file gas_solid_catalytic_wall.h.

template<typename Chemistry>

VariableIndex GRINS::GasSolidCatalyticWall< Chemistry >::_reactant_gas_var_idx

protected

Definition at line 78 of file gas_solid_catalytic_wall.h.

template<typename Chemistry>

const unsigned int GRINS::GasSolidCatalyticWall< Chemistry >::_reactant_solid_species_idx

protected

Definition at line 80 of file gas_solid_catalytic_wall.h.

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

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

Public Member Functions

Protected Attributes

Detailed Description

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

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation

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