Virtual class used to implement properties that are characteristic of a porous layer. More...

#include <porous_layer.h>

Inheritance diagram for FuelCellShop::Layer::PorousLayer< dim >:

Collaboration diagram for FuelCellShop::Layer::PorousLayer< dim >:

Public Member Functions
Initalization
void	set_gases_and_compute (std::vector< FuelCellShop::Material::PureGas * > &gases_in, const double &pressure_in, const double &temperature_in)
	Member function used to store all the gases that are in the pore space in the gas diffusion layer as well as their temperature [`Kelvin`] and total pressure [`atm`]. More...

void	compute_gas_diffusion (FuelCellShop::Material::PureGas solute_gas, FuelCellShop::Material::PureGas solvent_gas)
	Member function used to compute bulk diffusion coefficients and derivatives w.r.t temperature for non-isothermal case and store inside the layer. More...

void	set_gases (std::vector< FuelCellShop::Material::PureGas * > &gases_in, const double &pressure_in)
	Member function used to store all the gases that are in the pore space in the porous layer. More...

void	set_gas_mixture (FuelCellShop::Material::GasMixture &rgas_mixture)
	Set `gas_mixture`. More...

void	set_porosity_permeability_tortuosity_booleans (const bool &rporosity_is_constant, const bool &rpermeability_is_constant, const bool &rtortuosity_is_constant)
	Set. More...

void	set_temperature (const SolutionVariable &T_in)
	Member function used to set the temperature [`Kelvin`] at every quadrature point inside the cell. More...

void	set_saturation (const SolutionVariable &s_in)
	Member function used to set the liquid water saturation at every quadrature point inside the cell. More...

Accessors and info
FuelCellShop::Material::PureGas *	get_gas_pointer (int index) const
	Return the FuelCellShop::Material::PureGas pointer that is stored inside the class in the ith position. More...

std::vector < FuelCellShop::Material::PureGas * >	get_gases () const
	Returns the vector of FuelCellShop::Material::PureGas pointers stored in the porous layer. More...

const FuelCellShop::Material::GasMixture *const	get_gas_mixture () const
	This function returns `gas_mixture`. More...

void	get_gas_index (FuelCellShop::Material::PureGas *gas_type, int &index) const
	Return the gas index in the GDL class. More...

void	get_T_and_p (double &T, double &p) const
	Return the constant temperature [`Kelvin`] and constant pressure [`atm`] inside the layer. More...

void	get_p (double &p) const
	Return the constant pressure [`atm`] inside the layer. More...

const bool &	get_porosity_is_constant () const
	This function returns `porosity_is_constant`. More...

const bool &	get_permeability_is_constant () const
	This function returns `permeability_is_constant`. More...

const bool &	get_tortuosity_is_constant () const
	This function returns `tortuosity_is_constant`. More...

double	get_porosity () const
	This function computes constant porosity in quadrature points of a mesh entity. More...

void	get_porosity (std::vector< double > &dst) const
	This function computes constant porosity in quadrature points of a mesh entity. More...

void	get_porosity (std::vector< double > &dst, const std::vector< Point< dim > > &points) const
	This function computes variable porosity in quadrature points of a mesh entity. More...

void	get_permeability (std::vector< SymmetricTensor< 2, dim > > &dst) const
	This function computes constant permeability in quadrature points of a mesh entity. More...

void	get_permeability (std::vector< SymmetricTensor< 2, dim > > &dst, const std::vector< Point< dim > > &points) const
	This function computes variable permeability in quadrature points of a mesh entity. More...

void	get_SQRT_permeability (std::vector< SymmetricTensor< 2, dim > > &dst) const
	This function computes square root of constant permeability in quadrature points of a mesh entity. More...

void	get_SQRT_permeability (std::vector< SymmetricTensor< 2, dim > > &dst, const std::vector< Point< dim > > &points) const
	This function computes square root of variable permeability in quadrature points of a mesh entity. More...

void	get_permeability_INV (std::vector< SymmetricTensor< 2, dim > > &dst) const
	This function computes inverse of constant permeability in quadrature points of a mesh entity. More...

void	get_permeability_INV (std::vector< SymmetricTensor< 2, dim > > &dst, const std::vector< Point< dim > > &points) const
	This function computes inverse of variable permeability in quadrature points of a mesh entity. More...

void	get_SQRT_permeability_INV (std::vector< SymmetricTensor< 2, dim > > &dst) const
	This function computes inverse of square root of constant permeability in quadrature points of a mesh entity. More...

void	get_SQRT_permeability_INV (std::vector< SymmetricTensor< 2, dim > > &dst, const std::vector< Point< dim > > &points) const
	This function computes inverse of square root of variable permeability in quadrature points of a mesh entity. More...

void	get_Forchheimer_permeability (std::vector< SymmetricTensor< 2, dim > > &dst) const
	This function computes constant Forchheimer permeability in quadrature points of a mesh entity. More...

void	get_Forchheimer_permeability (std::vector< SymmetricTensor< 2, dim > > &dst, const std::vector< Point< dim > > &points) const
	This function computes variable Forchheimer permeability in quadrature points of a mesh entity. More...

void	get_tortuosity (std::vector< SymmetricTensor< 2, dim > > &dst) const
	This function computes constant tortuosity in quadrature points of a mesh entity. More...

void	get_tortuosity (std::vector< SymmetricTensor< 2, dim > > &dst, const std::vector< Point< dim > > &points) const
	This function computes variable tortuosity in quadrature points of a mesh entity. More...

virtual void	print_layer_properties () const
	This member function is a virtual class that can be used to output to screen information from the layer. More...

Computing member functions
void	effective_gas_diffusion_coefficient (const FuelCellShop::Material::PureGas solute_gas, const FuelCellShop::Material::PureGas solvent_gas, const SolutionVariable &T_in, const SolutionVariable &p_in, std::vector< double > &D_b) const
	Member function used to compute diffusionfor a solute_gas, solvent_gas combination at a given temperature and pressure using the Bousinesq approximation is used in order to compute the bulk combined diffusivity such that $\frac{1}{D} = \frac{1}{D_{i,j}} + \frac{1}{D^K_i}$ where $D_{i,j}$ is the molecular diffusion coefficient obtained using effective_molecular_gas_diffusion_coefficient and is the Knudsen diffusivity computed using Knudsen_diffusion. More...

void	effective_gas_diffusion_coefficient (const FuelCellShop::Material::PureGas solute_gas, const FuelCellShop::Material::PureGas solvent_gas, const SolutionVariable &T_in, const SolutionVariable &p_in, std::vector< double > &D_b, std::vector< double > &dD_b_dT) const
	Member function used to compute diffusion and its derivative for a solute_gas, solvent_gas combination at a given temperature and pressure using the Bousinesq approximation is used in order to compute the bulk combined diffusivity such that $\frac{1}{D} = \frac{1}{D_{i,j}} + \frac{1}{D^K_i}$ where $D_{i,j}$ is the molecular diffusion coefficient obtained using effective_molecular_gas_diffusion_coefficient and is the Knudsen diffusivity computed using Knudsen_diffusion. More...

void	effective_molecular_gas_diffusion_coefficient (const FuelCellShop::Material::PureGas solute_gas, const FuelCellShop::Material::PureGas solvent_gas, const SolutionVariable &T_in, const SolutionVariable &p_in, std::vector< double > &D_m) const
	Member function used to compute molecular diffusion for a solute_gas, solvent_gas combination at a given temperature and pressure. More...

void	effective_molecular_gas_diffusion_coefficient (const FuelCellShop::Material::PureGas solute_gas, const FuelCellShop::Material::PureGas solvent_gas, const SolutionVariable &T_in, const SolutionVariable &p_in, std::vector< double > &D_m, std::vector< double > &dD_m_dT) const
	Member function used to compute molecular diffusion for a solute_gas, solvent_gas combination at a given temperature and pressure. More...

void	Knudsen_diffusion (const FuelCellShop::Material::PureGas *solute_gas, const SolutionVariable &T_in, std::vector< double > &D_k) const
	Member function used to compute the Knudsen diffusivity in the layer. More...

void	Knudsen_diffusion (const FuelCellShop::Material::PureGas *solute_gas, const SolutionVariable &T_in, std::vector< double > &D_k, std::vector< double > &dD_k_dT) const
	Member function used to compute the Knudsen diffusivity in the layer. More...

virtual void	pcapillary (std::vector< double > &) const
	Compute $p_c \quad \left[ dyne \cdot cm^{-2}\right]$ , at all quadrature points in the cell. More...

Public Member Functions inherited from FuelCellShop::Layer::BaseLayer< dim >
virtual void	set_derivative_flags (const std::vector< VariableNames > &flags)
	Set the variables for which you would like to compute the derivatives. More...

void	set_position (const std::vector< Point< dim > > &p)
	Member function used by some applications such as dummyGDL in order to know which value to return. More...

void	set_local_material_id (const unsigned int &id)
	Function for setting local material id, for unit testing purposes. More...

void	unset_local_material_id ()
	Function for unsetting local material id, so that it isn't incorrectly used later Once the key is "unset" to some invalid value, an error will be thrown if the key is requested again without being set. More...

virtual void	set_constant_solution (const double &value, const VariableNames &name)
	Set those solution variables which are constant in the particular application. More...

virtual void	set_solution (const std::vector< SolutionVariable > &)
	If the effective properties in the layer depend on the solution, the solution for a given cell should be passed to the class using this member function. More...

bool	belongs_to_material (const unsigned int material_id)
	Check if a given cell belongs to the catalyst layer and assign. More...

const std::string &	name_layer () const
	Return the name of the layer. More...

virtual const std::type_info &	get_base_type () const
	This member function return the name of the type of layer, i.e. More...

virtual bool	test_layer ()
	This virtual class should be used for any derived class to be able to test the functionality of the class. More...

std::vector< unsigned int >	get_material_ids ()
	Return the local material id of the layer. More...

unsigned int	local_material_id () const
	Return the local material id of the layer, performs a check. More...

Protected Member Functions
Constructors, destructor, and initialization
	PorousLayer (const std::string &name)
	Constructor. More...

	PorousLayer ()
	Constructor. More...

	PorousLayer (const std::string &name, FuelCellShop::Material::GasMixture &gas_mixture)
	Constructor. More...

virtual	~PorousLayer ()
	Destructor. More...

virtual void	declare_parameters (const std::string &name, ParameterHandler &param) const
	Declare parameters for a parameter file. More...

virtual void	declare_parameters (ParameterHandler &param) const
	Declare parameters for a parameter file. More...

virtual void	initialize (ParameterHandler &param)
	Member function used to read in data and initialize the necessary data to compute the coefficients. More...

Minor functions
void	print_caller_name (const std::string &caller_name) const
	This function is used to print out the name of another function that has been declared in the scope of this class, but not yet been implemented. More...

Internal accessors and info
virtual void	gas_diffusion_coefficients (Table< 2, double > &) const
	Return the molecular diffusivty all the gases assigned to the layer using set_gases_and_compute. More...

virtual void	derivative_gas_diffusion_coefficients (std::vector< Table< 2, double > > &) const
	Return the derivative of the molecular diffusion coefficient with respect to the derivative flags for all the gases assigned to the layer using set_gases_and_compute. More...

Protected Member Functions inherited from FuelCellShop::Layer::BaseLayer< dim >
	BaseLayer ()
	Constructor. More...

	BaseLayer (const std::string &name)
	Constructor. More...

virtual	~BaseLayer ()
	Destructor. More...

virtual void	set_parameters (const std::string &object_name, const std::vector< std::string > &name_dvar, const std::vector< double > &value_dvar, ParameterHandler &param)
	Member function used to change the values in the parameter file for a given list of parameters. More...

virtual void	set_parameters (const std::vector< std::string > &name_dvar, const std::vector< double > &value_dvar, ParameterHandler &param)
	Set parameters in parameter file. More...

Protected Attributes
Layer properties
FuelCellShop::Material::GasMixture *	gas_mixture
	Gas mixture. More...

std::vector < FuelCellShop::Material::PureGas * >	gases
	Gases inside a porous layer. More...

bool	porosity_is_constant
	Variable defining if the porosity is constant. More...

bool	permeability_is_constant
	Variable defining if the permeability is constant. More...

bool	tortuosity_is_constant
	Variable defining if the tortuosity is constant. More...

double	porosity
	User defined constant porosity. More...

double	Knudsen_radius
	Parameter used to define Knudsen pore radius. More...

SymmetricTensor< 2, dim >	permeability
	User defined constant permeability, m^2. More...

SymmetricTensor< 2, dim >	SQRT_permeability
	Square root of user defined constant permeability, m. More...

SymmetricTensor< 2, dim >	permeability_INV
	Inverse of user defined constant permeability, 1/m^2. More...

SymmetricTensor< 2, dim >	SQRT_permeability_INV
	Inverse of square root of user defined constant permeability, 1/m. More...

SymmetricTensor< 2, dim >	Forchheimer_permeability
	User defined constant Forchheimer permeability, 1/m. More...

SymmetricTensor< 2, dim >	tortuosity
	User defined constant tortuosity. More...

std::string	diffusion_species_name
	If GDL properties are stored inside the class (e.g DummyGDL) then, return the property stored under coefficient_name name. More...

double	temperature
	Temperature [`K`] used to compute gas diffusivity. More...

double	pressure
	Total pressure [`atm`] used to compute gas diffusivity. More...

SolutionVariable	T_vector
	Temperature at every quadrature point inside the cell. More...

SolutionVariable	s_vector
	Liquid water saturation at every quadrature point inside the cell. More...

Table< 2, double >	D_ECtheory
	Tensor of diffusion coefficients – This are computed with setting up the gas so that they do not need to be recomputed all the time. More...

std::vector< Table< 2, double > >	dD_ECtheory_dx
	Vector of tensors for the derivative of the diffusion coefficients – This are computed with setting up the gas so that they do not need to be recomputed all the time. More...

std::vector< double >	D_bulk
	Vector of bulk diffusion coefficients at every quadrature point inside the cell. More...

std::vector< double >	dD_bulk_dT
	Vector of derivative of bulk diffusion coefficients w.r.t temperature, at every quadrature point inside the cell. More...

Protected Attributes inherited from FuelCellShop::Layer::BaseLayer< dim >
const std::string	name
	Name of the layer. More...

std::vector< unsigned int >	material_ids
	List of material IDs that belong to the layer. More...

std::vector< Point< dim > >	point
	Coordinates of the point where we would like to compute the effective properties. More...

std::vector< VariableNames >	derivative_flags
	Flags for derivatives: These flags are used to request derivatives. More...

std::map< VariableNames, double >	constant_solutions
	Map storing values of solution variables constant in a particular application. More...

Detailed Description

template<int dim>
class FuelCellShop::Layer::PorousLayer< dim >

Virtual class used to implement properties that are characteristic of a porous layer.

This layer is used to setup information needed to compute molecular diffusivities and Knudsen diffusivities (pending implementation) for any gases. Therefore, it provides member functions to specify the gases that exist within the layer as well as their temperature and pressure.

It also inherits all properties of BaseLayer.

Note you cannot make an object of this class since it still lacks many member functions that are necessary for a full layer implementation.

With respect to gases: you can either use gases or gas_mixture. Whatever you prefer. However, gas_mixture offers more methods in its structure. You cannot use both.

Author: M. Secanell, 2013; M. Bhaiya, 2013; V. Zingan, 2013

Constructor & Destructor Documentation

template<int dim>

FuelCellShop::Layer::PorousLayer< dim >::PorousLayer ( const std::string & name )

inlineprotected

Constructor.

References temperature_of_REV, and total_pressure.

template<int dim>

FuelCellShop::Layer::PorousLayer< dim >::PorousLayer ( )

inlineprotected

Constructor.

Warning: For internal use only.

template<int dim>

FuelCellShop::Layer::PorousLayer< dim >::PorousLayer	(	const std::string &	name,
		FuelCellShop::Material::GasMixture &	gas_mixture
	)

inlineprotected

Constructor.

References temperature_of_REV, and total_pressure.

template<int dim>

virtual FuelCellShop::Layer::PorousLayer< dim >::~PorousLayer ( )

inlineprotectedvirtual

Destructor.

Member Function Documentation

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::compute_gas_diffusion	(	FuelCellShop::Material::PureGas *	solute_gas,
		FuelCellShop::Material::PureGas *	solvent_gas
	)

Member function used to compute bulk diffusion coefficients and derivatives w.r.t temperature for non-isothermal case and store inside the layer.

This function takes solute and solvent gas as an input argument (Fick's diffusion model). Before calling this function, pressure [atm] must be set inside the layer, and temperature [Kelvin] must be set using set_temperature method.

Warning: This will not work if we want to use layer in parallel!

template<int dim>

virtual void FuelCellShop::Layer::PorousLayer< dim >::declare_parameters	(	const std::string &	name,
		ParameterHandler &	param
	)		const

protectedvirtual

Declare parameters for a parameter file.

Reimplemented from FuelCellShop::Layer::BaseLayer< dim >.

template<int dim>

virtual void FuelCellShop::Layer::PorousLayer< dim >::declare_parameters ( ParameterHandler & param ) const

inlineprotectedvirtual

Declare parameters for a parameter file.

Note: Since there are two declare parameters, this one is hidden by the former, so it has to be redefined otherwise it cannot be called :(

Reimplemented from FuelCellShop::Layer::BaseLayer< dim >.

Reimplemented in FuelCellShop::Layer::GasDiffusionLayer< dim >, FuelCellShop::Layer::MultiScaleCL< dim >, FuelCellShop::Layer::MultiScaleCL< deal_II_dimension >, FuelCellShop::Layer::DesignMPL< dim >, FuelCellShop::Layer::SGL24BA< dim >, FuelCellShop::Layer::SGL24BC< dim >, FuelCellShop::Layer::DummyGDL< dim >, and FuelCellShop::Layer::DummyCL< dim >.

template<int dim>

virtual void FuelCellShop::Layer::PorousLayer< dim >::derivative_gas_diffusion_coefficients ( std::vector< Table< 2, double > > & ) const

protectedvirtual

Return the derivative of the molecular diffusion coefficient with respect to the derivative flags for all the gases assigned to the layer using set_gases_and_compute.

This routine uses FuelCellShop::Mixture::ChapmanEnskog to compute the binary diffusivity for each gas with respect to each other. If the layer contains three gases the vector returns D12, D13, D23. For 4 gases, it returns D12, D13, D14, D23, D24, D34 . Derivatives are provided in the vector w.r.t following parameters in the same order as below:

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::effective_gas_diffusion_coefficient	(	const FuelCellShop::Material::PureGas *	solute_gas,
		const FuelCellShop::Material::PureGas *	solvent_gas,
		const SolutionVariable &	T_in,
		const SolutionVariable &	p_in,
		std::vector< double > &	D_b
	)		const

Member function used to compute diffusionfor a solute_gas, solvent_gas combination at a given temperature and pressure using the Bousinesq approximation is used in order to compute the bulk combined diffusivity such that

$\frac{1}{D} = \frac{1}{D_{i,j}} + \frac{1}{D^K_i}$

where $D_{i,j}$ is the molecular diffusion coefficient obtained using effective_molecular_gas_diffusion_coefficient and $ D^K_i $ is the Knudsen diffusivity computed using Knudsen_diffusion.

By defaul the diffusion coefficient is computed using the pore radius in:

* subsection Fuel Cell Data
*   subsection layer_name
*     set Knudsen pore radius [um] = 1e6
*   end
* end
* 

If you would like to neglect the diffusion coefficient, simply make this value very large (the default basically results on negligible Knudsen effects – the value is: 161.0949 (SI) compared to 1e-5 (SI) for molecular diff for gases – Note that you add the inverse of these values).

Note: If you would like

References

[1] Pollard and Present, On gaseous self-diffusion in long capillary tubes, Physical Review, Vol 73, No. 7, 1948
[2] Burganos, Gas diffusion in random binary media, Journal of Chemical Physics, 109(16), 1998

Author: M. Secanell, 2014

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::effective_gas_diffusion_coefficient	(	const FuelCellShop::Material::PureGas *	solute_gas,
		const FuelCellShop::Material::PureGas *	solvent_gas,
		const SolutionVariable &	T_in,
		const SolutionVariable &	p_in,
		std::vector< double > &	D_b,
		std::vector< double > &	dD_b_dT
	)		const

Member function used to compute diffusion and its derivative for a solute_gas, solvent_gas combination at a given temperature and pressure using the Bousinesq approximation is used in order to compute the bulk combined diffusivity such that

$\frac{1}{D} = \frac{1}{D_{i,j}} + \frac{1}{D^K_i}$

where $D_{i,j}$ is the molecular diffusion coefficient obtained using effective_molecular_gas_diffusion_coefficient and $ D^K_i $ is the Knudsen diffusivity computed using Knudsen_diffusion.

References

[1] Pollard and Present, On gaseous self-diffusion in long capillary tubes, Physical Review, Vol 73, No. 7, 1948
[2] Burganos, Gas diffusion in random binary media, Journal of Chemical Physics, 109(16), 1998

Author: M. Secanell, 2014

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::effective_molecular_gas_diffusion_coefficient	(	const FuelCellShop::Material::PureGas *	solute_gas,
		const FuelCellShop::Material::PureGas *	solvent_gas,
		const SolutionVariable &	T_in,
		const SolutionVariable &	p_in,
		std::vector< double > &	D_m
	)		const

Member function used to compute molecular diffusion for a solute_gas, solvent_gas combination at a given temperature and pressure.

Author: M. Secanell, 2014

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::effective_molecular_gas_diffusion_coefficient	(	const FuelCellShop::Material::PureGas *	solute_gas,
		const FuelCellShop::Material::PureGas *	solvent_gas,
		const SolutionVariable &	T_in,
		const SolutionVariable &	p_in,
		std::vector< double > &	D_m,
		std::vector< double > &	dD_m_dT
	)		const

Member function used to compute molecular diffusion for a solute_gas, solvent_gas combination at a given temperature and pressure.

Author: M. Secanell, 2014

template<int dim>

virtual void FuelCellShop::Layer::PorousLayer< dim >::gas_diffusion_coefficients ( Table< 2, double > & ) const

protectedvirtual

Return the molecular diffusivty all the gases assigned to the layer using set_gases_and_compute.

This routine uses FuelCellShop::Mixture::ChapmanEnskog to compute the binary diffusivity for each gas with respect to each other. If the layer contains three gases the vector returns D12, D13, D23. For 4 gases, it returns D12, D13, D14, D23, D24, D34.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::get_Forchheimer_permeability ( std::vector< SymmetricTensor< 2, dim > > & dst ) const

This function computes constant Forchheimer permeability in quadrature points of a mesh entity.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::get_Forchheimer_permeability	(	std::vector< SymmetricTensor< 2, dim > > &	dst,
		const std::vector< Point< dim > > &	points
	)		const

This function computes variable Forchheimer permeability in quadrature points of a mesh entity.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::get_gas_index	(	FuelCellShop::Material::PureGas *	gas_type,
		int &	index
	)		const

Return the gas index in the GDL class.

If the gas does not exist in the GDL it will return -1 as the index.

template<int dim>

const FuelCellShop::Material::GasMixture* const FuelCellShop::Layer::PorousLayer< dim >::get_gas_mixture ( ) const

inline

This function returns gas_mixture.

template<int dim>

FuelCellShop::Material::PureGas* FuelCellShop::Layer::PorousLayer< dim >::get_gas_pointer ( int index ) const

inline

Return the FuelCellShop::Material::PureGas pointer that is stored inside the class in the ith position.

The diffusion coefficient D[i][j] would be the diffusion coefficient of gas type i with j.

template<int dim>

std::vector<FuelCellShop::Material::PureGas*> FuelCellShop::Layer::PorousLayer< dim >::get_gases ( ) const

inline

Returns the vector of FuelCellShop::Material::PureGas pointers stored in the porous layer.

This method is useful to access all gases which are being transported inside a particular porous layer object.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::get_p ( double & p ) const

inline

Return the constant pressure [atm] inside the layer.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::get_permeability ( std::vector< SymmetricTensor< 2, dim > > & dst ) const

This function computes constant permeability in quadrature points of a mesh entity.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::get_permeability	(	std::vector< SymmetricTensor< 2, dim > > &	dst,
		const std::vector< Point< dim > > &	points
	)		const

This function computes variable permeability in quadrature points of a mesh entity.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::get_permeability_INV ( std::vector< SymmetricTensor< 2, dim > > & dst ) const

This function computes inverse of constant permeability in quadrature points of a mesh entity.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::get_permeability_INV	(	std::vector< SymmetricTensor< 2, dim > > &	dst,
		const std::vector< Point< dim > > &	points
	)		const

This function computes inverse of variable permeability in quadrature points of a mesh entity.

template<int dim>

const bool& FuelCellShop::Layer::PorousLayer< dim >::get_permeability_is_constant ( ) const

inline

This function returns permeability_is_constant.

template<int dim>

double FuelCellShop::Layer::PorousLayer< dim >::get_porosity ( ) const

inline

This function computes constant porosity in quadrature points of a mesh entity.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::get_porosity ( std::vector< double > & dst ) const

inline

This function computes constant porosity in quadrature points of a mesh entity.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::get_porosity	(	std::vector< double > &	dst,
		const std::vector< Point< dim > > &	points
	)		const

This function computes variable porosity in quadrature points of a mesh entity.

template<int dim>

const bool& FuelCellShop::Layer::PorousLayer< dim >::get_porosity_is_constant ( ) const

inline

This function returns porosity_is_constant.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::get_SQRT_permeability ( std::vector< SymmetricTensor< 2, dim > > & dst ) const

This function computes square root of constant permeability in quadrature points of a mesh entity.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::get_SQRT_permeability	(	std::vector< SymmetricTensor< 2, dim > > &	dst,
		const std::vector< Point< dim > > &	points
	)		const

This function computes square root of variable permeability in quadrature points of a mesh entity.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::get_SQRT_permeability_INV ( std::vector< SymmetricTensor< 2, dim > > & dst ) const

This function computes inverse of square root of constant permeability in quadrature points of a mesh entity.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::get_SQRT_permeability_INV	(	std::vector< SymmetricTensor< 2, dim > > &	dst,
		const std::vector< Point< dim > > &	points
	)		const

This function computes inverse of square root of variable permeability in quadrature points of a mesh entity.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::get_T_and_p	(	double &	T,
		double &	p
	)		const

inline

Return the constant temperature [Kelvin] and constant pressure [atm] inside the layer.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::get_tortuosity ( std::vector< SymmetricTensor< 2, dim > > & dst ) const

This function computes constant tortuosity in quadrature points of a mesh entity.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::get_tortuosity	(	std::vector< SymmetricTensor< 2, dim > > &	dst,
		const std::vector< Point< dim > > &	points
	)		const

This function computes variable tortuosity in quadrature points of a mesh entity.

template<int dim>

const bool& FuelCellShop::Layer::PorousLayer< dim >::get_tortuosity_is_constant ( ) const

inline

This function returns tortuosity_is_constant.

template<int dim>

virtual void FuelCellShop::Layer::PorousLayer< dim >::initialize ( ParameterHandler & param )

protectedvirtual

Member function used to read in data and initialize the necessary data to compute the coefficients.

Reimplemented from FuelCellShop::Layer::BaseLayer< dim >.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::Knudsen_diffusion	(	const FuelCellShop::Material::PureGas *	solute_gas,
		const SolutionVariable &	T_in,
		std::vector< double > &	D_k
	)		const

Member function used to compute the Knudsen diffusivity in the layer.

Author: M. Secanell, 2014

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::Knudsen_diffusion	(	const FuelCellShop::Material::PureGas *	solute_gas,
		const SolutionVariable &	T_in,
		std::vector< double > &	D_k,
		std::vector< double > &	dD_k_dT
	)		const

Member function used to compute the Knudsen diffusivity in the layer.

Author: M. Secanell, 2014

template<int dim>

virtual void FuelCellShop::Layer::PorousLayer< dim >::pcapillary ( std::vector< double > & ) const

inlinevirtual

Compute $p_c \quad \left[ dyne \cdot cm^{-2}\right]$ , at all quadrature points in the cell.

Reimplemented in FuelCellShop::Layer::CatalystLayer< dim >, FuelCellShop::Layer::GasDiffusionLayer< dim >, FuelCellShop::Layer::MicroPorousLayer< dim >, FuelCellShop::Layer::DesignMPL< dim >, FuelCellShop::Layer::ConventionalCL< dim >, and FuelCellShop::Layer::DesignFibrousGDL< dim >.

References FcstUtilities::log.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::print_caller_name ( const std::string & caller_name ) const

protected

This function is used to print out the name of another function that has been declared in the scope of this class, but not yet been implemented.

template<int dim>

virtual void FuelCellShop::Layer::PorousLayer< dim >::print_layer_properties ( ) const

virtual

This member function is a virtual class that can be used to output to screen information from the layer.

This function should be re-implemented in each layer with data that is relevant in each case.

Reimplemented from FuelCellShop::Layer::BaseLayer< dim >.

Reimplemented in FuelCellShop::Layer::MultiScaleCL< dim >, FuelCellShop::Layer::MultiScaleCL< deal_II_dimension >, and FuelCellShop::Layer::ConventionalCL< dim >.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::set_gas_mixture ( FuelCellShop::Material::GasMixture & rgas_mixture )

inline

Set gas_mixture.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::set_gases	(	std::vector< FuelCellShop::Material::PureGas * > &	gases_in,
		const double &	pressure_in
	)

inline

Member function used to store all the gases that are in the pore space in the porous layer.

Besides this, it takes total pressure [atm] as input. This method can be particulary used for non-isothermal applications when temperature changes in every cell, hence this method doesn't compute the diffusion coefficients but store the relevant constant data.

Note: Please ensure that in the case of Fick's diffusion model, solvent gas should be the last object inside the input vector. On not ensuring this may not give any errors, as this function may also be used for multi-component gas transport but for fickian diffusion, it may give wrong results.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::set_gases_and_compute	(	std::vector< FuelCellShop::Material::PureGas * > &	gases_in,
		const double &	pressure_in,
		const double &	temperature_in
	)

Member function used to store all the gases that are in the pore space in the gas diffusion layer as well as their temperature [Kelvin] and total pressure [atm].

Then, when computing diffusion, the class will return each one of the diffusion coefficients for the gases. For example, if we pass a vector with for gases, say oxygen and nitrogen, effective_gas_diffusivity will return $D_{N_2,O_2}$

Deprecated:

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::set_porosity_permeability_tortuosity_booleans	(	const bool &	rporosity_is_constant,
		const bool &	rpermeability_is_constant,
		const bool &	rtortuosity_is_constant
	)

inline

Set.

porosity_is_constant,
permeability_is_constant,
tortuosity_is_constant.

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::set_saturation ( const SolutionVariable & s_in )

inline

Member function used to set the liquid water saturation at every quadrature point inside the cell.

This function should particulary be used in the case of multi-phase application.

References FuelCellShop::SolutionVariable::get_variablename(), and liquid_water_saturation.

Here is the call graph for this function:

template<int dim>

void FuelCellShop::Layer::PorousLayer< dim >::set_temperature ( const SolutionVariable & T_in )

inline

Member function used to set the temperature [Kelvin] at every quadrature point inside the cell.

This function should particulary be used in the case of non-isothermal application.

References FuelCellShop::SolutionVariable::get_variablename(), and temperature_of_REV.

Here is the call graph for this function:

Member Data Documentation

template<int dim>

std::vector<double> FuelCellShop::Layer::PorousLayer< dim >::D_bulk

protected

Vector of bulk diffusion coefficients at every quadrature point inside the cell.

This particulary stores non-isothermal diffusion coefficients computed using compute_gas_diffusion method.

template<int dim>

Table< 2, double > FuelCellShop::Layer::PorousLayer< dim >::D_ECtheory

protected

Tensor of diffusion coefficients – This are computed with setting up the gas so that they do not need to be recomputed all the time.

template<int dim>

std::vector<double> FuelCellShop::Layer::PorousLayer< dim >::dD_bulk_dT

protected

Vector of derivative of bulk diffusion coefficients w.r.t temperature, at every quadrature point inside the cell.

These are computed using compute_gas_diffusion method.

template<int dim>

std::vector< Table< 2, double > > FuelCellShop::Layer::PorousLayer< dim >::dD_ECtheory_dx

protected

Vector of tensors for the derivative of the diffusion coefficients – This are computed with setting up the gas so that they do not need to be recomputed all the time.

template<int dim>

std::string FuelCellShop::Layer::PorousLayer< dim >::diffusion_species_name

protected

If GDL properties are stored inside the class (e.g DummyGDL) then, return the property stored under coefficient_name name.

template<int dim>

SymmetricTensor<2,dim> FuelCellShop::Layer::PorousLayer< dim >::Forchheimer_permeability

protected

User defined constant Forchheimer permeability, 1/m.

template<int dim>

FuelCellShop::Material::GasMixture* FuelCellShop::Layer::PorousLayer< dim >::gas_mixture

protected

Gas mixture.