This archive contains data that describes a dummy medium to run htrdr for a combustion application, in the general case of an heterogeneous medium. It also proposes a program to validate htrdr-combustion.
The provided combustion medium is an axis aligned cube of 0.2 m side length. The (x=0, y=0) position is the center of the bottom face of the cube, meaning that x and y coordinates vary in the [-0.1, 0.1] m range, while the z coordinate varies in the [0, 0.2] m range.
The laser sheet is 1 mm thick, is horizontal, and splits the medium at half height: the altitude of the laser sheet therefore varies in the [0.0995, 0.105] m range, while x and y coordinates vary in the [-0.1, 0.1] m range. The laser sheet is emitted from the y=0 plane, and propagates in the direction of positive y. Its wavelength is 532 nm, and its surface power density is 1000 W/m².
The sensor on which the flux will be computed is located in the center of the bottom limit of the combustion medium. It is a square of 5 cm side length: x and y coordinates over the sensor vary in the [-0.025, 0.025] m range, and z=0.
Boundaries of the combustion chamber are black, and the sensor is also a blackbody: whenever radiation reaches a boundary, it is fully absorbed. Only the flux that reaches the sensor is accounted for.
Soot aggregates absorb and scatter radiation within the limits of the combustion chamber. Their optical cross-sections are computed according to the RDG-FA theory, using the following data:
- refraction index: 1.60+0.75.i
- fractal dimension: 1.80
- fractal prefactor: 1.30
- primary particles diameter: 20 nm
- number of primary particles per aggregate: 100
The soot volume fraction follows the following inhomogeneous axisymmetric profile:
fv_max the maximum soot volumic fraction, reached at x=-L/2, y=0, z=L/2,
L the length of the cubic combustion chamber, x and y in [-L/2, L/2] and z in
fv_max is set at 10^-6 m³ of soot per m³.
dummy_medium directory contains all the data necessary to describe the
aforementioned medium within htrdr-combustion. The available files are:
tetra_mesh.sth: it stores the tetrahedral mesh that represents the heterogeneous medium (12 millions of tetrahedra).
refraction_index.atrri: this file defines the spectrally varying refractive indices of the medium.
thermodynamic_conditions.atrtp: it saves the per node physical properties of the tetrahedral mesh.
This directory also contains two Bash scripts that run the
program to compute a flux density map on this combustion medium
compute_flux_density_map.sh) or to draw an image of it (
These scripts can be a good starting point to study how to run htrdr for a
combustion application. Assuming that
htrdr-combustion is correctly installed
and registered into the current Bash shell, simply run the desired script. You
can finally use
htpp to convert the computed images into a PPM image that
can be then displayed with a regular image viewer.
~ $ bash htrdr-Combustion-Starter-Pack-0.0.0/dummy_medium/draw_image.sh ~ $ htpp -vm default -o dummy_image.ppm dummy_image_800x600x128.txt
sw_flux directory contains a simple program, written in Fortran, that was
designed to solve the problem of interest only (fixed geometric configuration
and analytical thermodynamic properties fields). However, two separate
Monte-Carlo algorithms have been implemented in order to compute the required
- a direct Monte-Carlo algorithm: optical trajectories are followed from the laser emission surface, until they reach the sensor (or are lost).
- a reverse Monte-Carlo algorithm: this is actually the one that was
htrdr-combustion. Optical paths start from the sensor, and a first-scattering contribution increases the weight of the realization at each scattering position in the medium.
To compile the executable simply invoke
make into the
Note that the build procedure assumes that the GNU Fortran compiler
gfortran) is installed. Once built, run the generated executable named
sw_flux. Both algorithms (direct and reverse) must agree, within the limits
of the statistical uncertainty. Furthermore, these results should also agree
with the flux computed by
htrdr-combustion via the script
~ $ cd htrdr-Combustion-Starter-Pack-0.0.0/sw_flux sw_flux $ make sw_flux $ ./sw_flux sw_flux $ cd ../dummy_medium/ dummy_medium $ bash compute_flux_density_map.sh
The result should be, in terms of surface power density:
1.028e-2 +/- 1e-5 W/m²
Copyright © 2021 |Méso|Star> (firstname.lastname@example.org)
Copyright © 2021 CNRS/RAPSODEE
htrdr: Combustion Starter Pack is released under the GPLv3+ license: GNU GPL
version 3 or later. You can freely study, modify or extend it. You are also
welcome to redistribute it under certain conditions; refer to the
license for details.