parallel_planes.cc File Reference

1-D radiative solver, ray tracing. More...

#include <cmath>
#include <iostream>
#include "../../src/c++/rad.h"

Include dependency graph for parallel_planes.cc:

Go to the source code of this file.

Functions
void	I_IT (const vector< double > &x, const double theta, const vector< double > &T, const vector< vector< double > > &kabs, const vector< vector< double > > &awts, const vector< double > &Ilo, const vector< double > &Ihi, vector< vector< double > > &I)
void	parallel_planes (rad *RAD, const double L, const int ntheta, const vector< double > &T, const double P, const vector< double > &fvsoot, const vector< double > &xH2O, const vector< double > &xCO2, const vector< double > &xCO, const vector< double > &xCH4, vector< double > &q, vector< double > &Q, vector< double > &x, vector< double > &xQ, const bool LzeroIbc=false)

Detailed Description

1-D radiative solver, ray tracing.

Definition in file parallel_planes.cc.

Function Documentation

I_IT()

void I_IT	(	const vector< double > &	x,
		const double	theta,
		const vector< double > &	T,
		const vector< vector< double > > &	kabs,
		const vector< vector< double > > &	awts,
		const vector< double > &	Ilo,
		const vector< double > &	Ihi,
		vector< vector< double > > &	I
	)

Compute intensity at each grid point for a given angle. Solution using an implicit trapazoid method.

Parameters

x	input: vector of grid points
theta	input: angle measured: grid between plates aligned with theta=0
T	input: vector of temperatures (K)
kabs	input: 2D vector of absorption coefficients for each gas [igrid, igas]
awts	input: 2D vector of weights for each gas [igrid, igas]
Ilo	input: vector of gas intensities on the lower plane.
Ihi	input: vector of gas intensities on the upper plane.
I	\ouput vector of vector of gas intensities in direction theta [igrid, igas]

Definition at line 26 of file parallel_planes.cc.

Here is the caller graph for this function:

parallel_planes()

void parallel_planes	(	rad *	RAD,
		const double	L,
		const int	ntheta,
		const vector< double > &	T,
		const double	P,
		const vector< double > &	fvsoot,
		const vector< double > &	xH2O,
		const vector< double > &	xCO2,
		const vector< double > &	xCO,
		const vector< double > &	xCH4,
		vector< double > &	q,
		vector< double > &	Q,
		vector< double > &	x,
		vector< double > &	xQ,
		const bool	LzeroIbc = false
	)

Compute radiative heat flux (q) and volumetric heat source (Q) profiles between parallel planes. Solution using simple ray tracing. Pass in T, P, composition profiles on equally spaced points including boundaries.

Parameters

RAD	input: rad object
L	input: distance (m) between planes
ntheta	input: number of ray angles
T	input: vector of temperatures at each point (K)
P	input: system pressure (Pa)
fvsoot	input: vector of soot volume fractions at each point
xH2O	input: vector of H2O mole fractions at each point
xCO2	input: vector of CO2 mole fractions at each point
xCO	input: vector of CO mole fractions at each point
xCH4	input: vector of CH4 mole fractions at each point
q	output: vector of heat fluxes at each point (W/m2)
Q	output: vector of heat sources at each point (W/m3)
x	output: vector of grid locations (m)
xQ	output: vector of Q grid locations (m)
LzeroIbc	input: bool true to make Ibc=0, false otherwise (default = false)

Definition at line 92 of file parallel_planes.cc.

Here is the call graph for this function:

Here is the caller graph for this function: