#include <rad_wsgg.h>

Inheritance diagram for rad_wsgg:

Collaboration diagram for rad_wsgg:

Public Member Functions
void	get_k_a (std::vector< double > &kabs, std::vector< double > &awts, const double T_dmb, const double P, const double fvsoot, const double xH2O, const double xCO2, const double xCO_not_used, const double xCH4_not_used)
	‍absorption coefficients for pure h2o (size 5 = nGGa)

void	get_k_a_oneband (double &kabs, double &awts, const int iband, const double T, const double P, const double fvsoot, const double xH2O, const double xCO2, const double xCO_not_used, const double xCH4_not_used)

	rad_wsgg ()

virtual	~rad_wsgg ()

Public Member Functions inherited from rad
int	get_nGG ()

int	get_nGGa ()

	rad (const int p_nGG, const int p_nGGa)

virtual	~rad ()

Static Protected Attributes
static const double	cCoefs []
	map to [i,j,k] of size ni,nj,nk = 4,5,5

static const double	dCoefs []

static const double	bco2 []
	‍map to [i,k] of size ni,nk = 4,5

static const double	bh2o []
	‍map to [i,j] of size ni,nj = 4,5

static const double	kco2 [] ={0.000000e+000, 3.388079e-002, 4.544269e-001, 4.680226e+000, 1.038439e+002}
	‍map to [i,j] of size ni,nj = 4,5

static const double	kh2o [] ={0.000000e+000, 7.703541e-002, 8.242941e-001, 6.854761e+000, 6.593653e+001}
	‍absorption coefficients for pure co2 (size 5 = nGGa)

Additional Inherited Members
Static Public Attributes inherited from rad
static constexpr double	sigma = 5.670367E-8
	Stephan-Boltzmann constant.

Protected Attributes inherited from rad
int	nGG
	number of gray gases, not including the clear gas

int	nGGa
	number of grey gases including the clear gas

Detailed Description

Class implementing rad_wsgg object. Radiation properties for the WSGG model (Bordbar 2020). Four gray gases and one clear gas are assumed. This is somewhat less accurate than the RCSLW model, but less computationally expensive. get_k_a, and get_k_a_oneband are the primary interfaces.

Definition at line 19 of file rad_wsgg.h.

Constructor & Destructor Documentation

◆ rad_wsgg()

rad_wsgg::rad_wsgg ( )

inline

Definition at line 61 of file rad_wsgg.h.

◆ ~rad_wsgg()

virtual rad_wsgg::~rad_wsgg ( )

inlinevirtual

Definition at line 64 of file rad_wsgg.h.

Member Function Documentation

◆ get_k_a()

void rad_wsgg::get_k_a	(	std::vector< double > &	kabs,
		std::vector< double > &	awts,
		const double	T_dmb,
		const double	P,
		const double	fvsoot,
		const double	xH2O,
		const double	xCO2,
		const double	xCO_not_used,
		const double	xCH4_not_used
	)

virtual

‍absorption coefficients for pure h2o (size 5 = nGGa)

This is the class interface function Given the gas state, set the k and a vectors. These can then be accessed by the user. return through arg list the local gray gas coefficients (kabs) and the local weights (awts).

Parameters

kabs	output: absorption coefficients (1/m) for nGG+1 (nGG gray gases + clear gas)
awts	output: weights (unitless; sums to 1) for nGG+1 (nGG gray gases + clear gas)
T_dmb	input: gas temperature (K)
P	input: pressure (Pa)
fvsoot	input: soot volume fraction = rho*Ysoot/rhosoot
xH2O	input: mole fraction H2O
xCO2	input: mole fraction CO2
xCO_not_used	input: mole fraction CO HERE FOR THE INTERFACE, NOT USED (... pass in 0.0)
xCH4_not_used	input: mole fraction CH4 HERE FOR THE INTERFACE, NOT USED (... pass in 0.0)

See documentation for rad_rcslw::F_albdf_soot for details about the soot absorption coefficient.

\(k_{soot} = F_s fv_{soot}T\), where \(F_s = 3.72 c_{soot}/C_2\), where \(C_2 = 0.014388\, m\cdot K\) and \(c_{soot} = 36\pi n k/[(n^2 - k^2 +2)^2 + 4 (n k)^2],\) where \(k\) is the real part of the complex refractive index, and \(n\) is the imaginary part. Using Shaddix's model for \(k\), \(n\): \(k=1.03\), \(n = 1.75\), giving \(F_s = 1817\, K^{-1}m^{-1}\). Reference: Williams, Shaddix, et al. Int. J. Heat and Mass Transfer 50:1616-1630 (2007),

Note, combining the WSGG model with four gray gases and one clear gas with a single gray soot component. The soot absorption coefficient is added to that for each gas (including the clear gas, since the soot absorption spectrum is "full"). The weights don't need to be changed. In the limit of no soot, we recover gas only, and in the limit of only soot, we recover the expected behavior.

\(dI_j/ds = (kg_j + ks)I_j + (kg_j+ks)a_jIb\)

This soot addition was suggested by H. Bordbar.

Implements rad.

Definition at line 270 of file rad_wsgg.cc.

Here is the call graph for this function:

◆ get_k_a_oneband()

void rad_wsgg::get_k_a_oneband	(	double &	kabs,
		double &	awts,
		const int	iband,
		const double	T_dmb,
		const double	P,
		const double	fvsoot,
		const double	xH2O,
		const double	xCO2,
		const double	xCO_not_used,
		const double	xCH4_not_used
	)

virtual

This is the class interface function Given the gas state, set the k and a vectors. These can then be accessed by the user. return through arg list the local gray gas coefficients (kabs) and the local weights (awts).

Parameters

kabs	output: absorption coefficient (1/m) for band/gas iband: ranges from 0 to nGG inclusive
awts	output: weight (unitless; total sums to 1) for band/gas iband: ranges from 0 to nGG inclusive
iband	input: which band to compute
T_dmb	input: gas temperature (K)
P	input: pressure (Pa)
fvsoot	input: soot volume fraction = rho*Ysoot/rhosoot
xH2O	input: mole fraction H2O
xCO2	input: mole fraction CO2
xCO_not_used	input: mole fraction CO HERE FOR THE INTERFACE, NOT USED (... pass in 0.0)
xCH4_not_used	input: mole fraction CH4 HERE FOR THE INTERFACE, NOT USED (... pass in 0.0)

See documentation for rad_rcslw::F_albdf_soot for details about the soot absorption coefficient.

\(k_{soot} = F_s fv_{soot}T\), where \(F_s = 3.72 c_{soot}/C_2\), where \(C_2 = 0.014388\, m\cdot K\) and \(c_{soot} = 36\pi n k/[(n^2 - k^2 +2)^2 + 4 (n k)^2],\) where \(k\) is the real part of the complex refractive index, and \(n\) is the imaginary part. Using Shaddix's model for \(k\), \(n\): \(k=1.03\), \(n = 1.75\), giving \(F_s = 1817\, K^{-1}m^{-1}\). Reference: Williams, Shaddix, et al. Int. J. Heat and Mass Transfer 50:1616-1630 (2007),

Note, combining the WSGG model with four gray gases and one clear gas with a single gray soot component. The soot absorption coefficient is added to that for each gas (including the clear gas, since the soot absorption spectrum is "full"). The weights don't need to be changed. In the limit of no soot, we recover gas only, and in the limit of only soot, we recover the expected behavior.

\(dI_j/ds = (kg_j + ks)I_j + (kg_j+ks)a_jIb\)

This soot addition was suggested by H. Bordbar.

Implements rad.

Definition at line 93 of file rad_wsgg.cc.

Here is the caller graph for this function:

Member Data Documentation

◆ cCoefs

const double rad_wsgg::cCoefs

staticprotected

Initial value:

={7.412956e-001, -5.244441e-001,  5.822860e-001, -2.096994e-001,  2.420312e-002,
                                   -9.412652e-001,  2.799577e-001, -7.672319e-001,  3.204027e-001, -3.910174e-002,
                                    8.531866e-001,  8.230754e-002,  5.289430e-001, -2.468463e-001,  3.109396e-002,
                                   -3.342806e-001,  1.474987e-001, -4.160689e-001,  1.697627e-001, -2.040660e-002,
                                    4.314362e-002, -6.886217e-002,  1.109773e-001, -4.208608e-002,  4.918817e-003,
                                    1.552073e-001, -4.862117e-001,  3.668088e-001, -1.055508e-001,  1.058568e-002,
                                    6.755648e-001,  1.409271e+000, -1.383449e+000,  4.575210e-001, -5.019760e-002,
                                   -1.125394e+000, -5.913199e-001,  9.085441e-001, -3.334201e-001,  3.842361e-002,
                                    6.040543e-001, -5.533854e-002, -1.733014e-001,  7.916083e-002, -9.893357e-003,
                                   -1.105453e-001,  4.646634e-002, -1.612982e-003, -3.539835e-003,  6.121277e-004,
                                    2.550242e-001,  3.805403e-001, -4.249709e-001,  1.429446e-001, -1.574075e-002,
                                   -6.065428e-001,  3.494024e-001,  1.853509e-001, -1.013694e-001,  1.302441e-002,
                                    8.123855e-001, -1.102009e+000,  4.046178e-001, -8.118223e-002,  6.298101e-003,
                                   -4.532290e-001,  6.784475e-001, -3.432603e-001,  8.830883e-002, -8.415221e-003,
                                    8.693093e-002, -1.306996e-001,  7.414464e-002, -2.029294e-002,  2.010969e-003,
                                   -3.451994e-002,  2.656726e-001, -1.225365e-001,  3.001508e-002, -2.820525e-003,
                                    4.112046e-001, -5.728350e-001,  2.924490e-001, -7.980766e-002,  7.996603e-003,
                                   -5.055995e-001,  4.579559e-001, -2.616436e-001,  7.648413e-002, -7.908356e-003,
                                    2.317509e-001, -1.656759e-001,  1.052608e-001, -3.219347e-002,  3.386965e-003,
                                   -3.754908e-002,  2.295193e-002, -1.600472e-002,  5.046318e-003, -5.364326e-004}

map to [i,j,k] of size ni,nj,nk = 4,5,5

Definition at line 14 of file rad_wsgg.h.

◆ dCoefs

const double rad_wsgg::dCoefs

staticprotected

Initial value:

={3.404288e-002,  6.523048e-002, -4.636852e-002,  1.386835e-002, -1.444993e-003,
                                   3.509457e-001,  7.465138e-001, -5.293090e-001,  1.594423e-001, -1.663261e-002,
                                   4.570740e+000,  2.168067e+000, -1.498901e+000,  4.917165e-001, -5.429990e-002,
                                   1.098169e+002, -5.092359e+001,  2.343236e+001, -5.163892e+000,  4.393889e-001}

Definition at line 35 of file rad_wsgg.h.

◆ bco2

const double rad_wsgg::bco2

staticprotected

Initial value:

={ 8.425766e-001, -1.442229e+000,  1.286974e+000, -5.202712e-001,  7.581559e-002,
                                 -3.023864e-002,  5.264245e-001, -6.209696e-001,  2.704755e-001, -4.090690e-002,
                                  1.070243e-001, -1.989596e-001,  3.101602e-001, -1.737230e-001,  3.081180e-002,
                                  3.108972e-002,  1.981489e-001, -2.543676e-001,  1.061331e-001, -1.498231e-002 }

‍map to [i,k] of size ni,nk = 4,5

Definition at line 42 of file rad_wsgg.h.

◆ bh2o

const double rad_wsgg::bh2o

staticprotected

Initial value:

={ 7.129509e-001, -1.378353e+000,  1.555028e+000, -6.636291e-001,  9.773674e-002,
                                  1.589917e-001,  5.635578e-002,  2.666874e-001, -2.040335e-001,  3.742408e-002,
                                 -1.196373e-001,  1.349665e+000, -1.544797e+000,  6.397595e-001, -9.153650e-002,
                                  3.078250e-001, -6.003555e-001,  4.441261e-001, -1.468813e-001,  1.824702e-002 }

‍map to [i,j] of size ni,nj = 4,5

Definition at line 47 of file rad_wsgg.h.

◆ kco2

const double rad_wsgg::kco2 ={0.000000e+000, 3.388079e-002, 4.544269e-001, 4.680226e+000, 1.038439e+002}

staticprotected

‍map to [i,j] of size ni,nj = 4,5

Definition at line 52 of file rad_wsgg.h.

◆ kh2o

const double rad_wsgg::kh2o ={0.000000e+000, 7.703541e-002, 8.242941e-001, 6.854761e+000, 6.593653e+001}

staticprotected

‍absorption coefficients for pure co2 (size 5 = nGGa)

Definition at line 54 of file rad_wsgg.h.

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

/Users/dol4/codes/radlib/src/c++/rad_wsgg.h
/Users/dol4/codes/radlib/src/c++/rad_wsgg.cc

Public Member Functions

Static Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ rad_wsgg()

◆ ~rad_wsgg()

Member Function Documentation

◆ get_k_a()

◆ get_k_a_oneband()

Member Data Documentation

◆ cCoefs

◆ dCoefs

◆ bco2

◆ bh2o

◆ kco2

◆ kh2o