Particles::MonteCarlo::DiffusionMonteCarloParameters Struct Reference

Precomputed quantities useful for the diffusion approximation in high-scattering opacity regions. We follow [81] Note that r_d in that manuscript should be (distance diffused)/(time elapsed) i.e. the paper is missing a normalization of r_d by delta_t. The upper bound of the integral in Eq (30) should just be sqrt(3*tau/4) All quantities in this struct are fixed; we could also just hard-code them to save us the on-the-fly calculation, or have a single instance of this struct in the global cache. More...

#include <Scattering.hpp>

Public Attributes
const double	MinimumOpacityForDiffusion = 3.0
const double	MaximumOpacityForCorrection = 11.0
const std::array< double, 21 >	BvsRhoForScattering
	Definition of the vector B_i, equation (35) More...
std::array< double, 1001 >	ScatteringRofP
	Storage for the function r(P) implicitly defined by Eq (29) ScatteringRofP[i] = r(0.001*i) Calculation performed in constructor.
std::array< double, 101 >	OpacityDependentCorrectionToRofP
	Storage for the opacity dependent correction on the right-hand side of Eq (30) We use 101 points between the min and max opacities defined above. Calculation performed in constructor.

Detailed Description

Precomputed quantities useful for the diffusion approximation in high-scattering opacity regions. We follow [81] Note that r_d in that manuscript should be (distance diffused)/(time elapsed) i.e. the paper is missing a normalization of r_d by delta_t. The upper bound of the integral in Eq (30) should just be sqrt(3*tau/4) All quantities in this struct are fixed; we could also just hard-code them to save us the on-the-fly calculation, or have a single instance of this struct in the global cache.

Member Data Documentation

BvsRhoForScattering

const std::array<double, 21> Particles::MonteCarlo::DiffusionMonteCarloParameters::BvsRhoForScattering

Initial value:

{
      1000., 18.74, 7.52, 4.75, 3.51, 2.78, 2.32, 2.00,  1.77,   1.60,     1.47,
      1.36,  1.28,  1.21, 1.15, 1.10, 1.07, 1.04, 1.019, 1.0027, 1.0000001}

Definition of the vector B_i, equation (35)

---

The documentation for this struct was generated from the following file:

• src/Evolution/Particles/MonteCarlo/Scattering.hpp