exojax
exojax copied to clipboard
HajimeKawahara
radiative transfer rethinking
In Paper I, we used 2 E_3 as a pure absorption transmission T, but not tried a numerical integral of mu e^-dtau/mu, for instance, using Gaussian quadrature etc, which might reduce a computational cost than the direct computation of 2 E_3. Even using mu=2/3, e^-dtau/mu, is not a bad approximation of 2 E3:
I think this is a first step to generalize RT to a scattering one.
Also, we would include scattering.
We can use jax.lax.linalg.tridiagonal_solve to solve the two-stream equations.
See an example.
in twostream branch, I am now developing the two stream solver (Toom/SH type).
This is the first example of the comparison with the pure absorption (using LART). It looks fine.