Performance Portable Opacity and Emissivity library for simulation codes
singularity-opac provides a uniform API for all opacity models. The following functions are provided
(here, \sigma is the frequency- and angle-dependent cross section in units of cm^2):
| Function | Expression | Description | Units |
|---|---|---|---|
| AbsorptionCoefficient | n \sigma |
Absorption coefficient | cm^{-1} |
| AngleAveragedAbsorptionCoefficient | \frac{1}{4 \pi}\int n \sigma d\Omega |
Absorption coefficient averaged over solid angle | cm^{-1} |
| EmissivityPerNuOmega | j_{\nu} = \frac{dE}{d^3x dt d\Omega d\nu} |
Frequency- and angle-dependent emissivity | erg cm^{-3} s^{-1} Sr^{-1} Hz^{-1} |
| EmissivityPerNu | \int j_{\nu} d\Omega |
Frequency-dependent emissivity | erg cm^{-3} s^{-1} Hz^{-1} |
| Emissivity | \int j_{\nu} d\nu d\Omega |
Total emissivity | erg cm^{-3} s^{-1} |
| NumberEmissivity | \int \frac{1}{h \nu} j_{\nu} d\Omega d\nu |
Total number emissivity | cm^{-3} s^{-1} |
| ThermalDistributionOfTNu | B_{\nu} = \frac{dE}{dA dt d\Omega d\nu} |
Specific intensity of thermal distribution | erg cm^{-2} s^{-1} Sr^{-1} Hz^{-1} |
| DThermalDistributionOfTNuDT | dB_{\nu}/dT |
Temperature derivative of specific intensity of thermal distribution | erg cm^{-2} s^{-1} Sr^{-1} Hz^{-1} K^{-1} |
| ThermalDistributionOfT | B = \int B_{\nu} d\Omega d\nu |
Frequency- and angle-integrated intensity of thermal distribution | erg cm^{-2} s^{-1} |
| ThermalNumberDistributionOfT | B = \int \frac{1}{h \nu} B_{\nu} d\Omega d\nu |
Frequency- and angle-integrated intensity of thermal distribution | erg cm^{-2} s^{-1} |
Note that the thermal radiation energy density u = 1/c ThermalDistributionOfT and the thermal radiation number density n = 1/c ThermalNumberDistributionOfT.
Internally singularity-opac always uses CGS units, as in the above table. However, arbitrary units are supported through the units modifier.
At its most basic:
git clone --recursive [email protected]:lanl/singularity-opac.git
cd singularity-opac
mkdir bin
cd bin
cmake ..
makeThen, in the singularity-opac root directory:
mkdir -p bin
cd bin
cmake -DSINGULARITY_BUILD_TESTS=ON ..
make -j
make testA number of options are avaialable for compiling:
| Option | Default | Comment |
|---|---|---|
| SINGULARITY_BUILD_TESTS | OFF | Build test infrastructure. |
| SINGULARITY_USE_HDF5 | ON | Enables HDF5. Required for Spiner opacities. |
| SINGULARITY_KOKKOS_IN_TREE | OFF | Force cmake to use Kokkos source included in tree. |
Currently, the MeanOpacity class, defined in singularity-opac/photons/mean_opacity_photons.hpp, supports loading grey Rosseland and Planck opacity data in an ASCII format. An example of this format is provided by singularity-opac/photons/example_ascii/kap_plaw.txt. The 1st row of the header has the number of density points, NRho, then the number of temperature points, NT. The 2nd (3rd) row of the header has min and max density (temperature) bounds. These bounds are inclusive, so the opacity data in the file should have evaluations at these min and max values. The rest of the ASCII file is a two-column table, where the 1st (2nd) column is Rosseland (Planck) opacity. The number of rows in each column is NRhoxNT, where density is the slow index and temperature is the fast index (thus the row index = temperature index
- NT x (density index), indexing from 0). Each opacity is assumed to be evaluated on log-spaced density and temperature grids, where these grids are defined by NRho, NT, and the (again inclusive) min and max bounds the header.
© 2021. Triad National Security, LLC. All rights reserved. This program was produced under U.S. Government contract 89233218CNA000001 for Los Alamos National Laboratory (LANL), which is operated by Triad National Security, LLC for the U.S. Department of Energy/National Nuclear Security Administration. All rights in the program are reserved by Triad National Security, LLC, and the U.S. Department of Energy/National Nuclear Security Administration. The Government is granted for itself and others acting on its behalf a nonexclusive, paid-up, irrevocable worldwide license in this material to reproduce, prepare derivative works, distribute copies to the public, perform publicly and display publicly, and to permit others to do so.