sandbox/ecipriano/test/radiation.c
Optically Thin Radiation Model
We test the calculation of the radiation using the optically-thin model at different temperautres.
We define variables that are necessary for the radiation model to compile.
#define NGS 2
const char* gas_species[2] = {"CO2", "H2O"};
#include "common.h"
#include "utils.h"
#include "radiation.h"
int main (void) {
The properties for the radiation model are gathered in the following struct.
OpticallyThinProperties otp;
otp.T = 1400.;
otp.P = 101325.;
otp.xH2O = 0.21;
otp.xCO2 = 0.1;
We set the radiation model, which is set to no_radiation by default.
divq_rad = optically_thin;
We loop over different temperatures and compute the \nabla\cdot\mathbf{q}_{rad} contribution.