A ray tracing code in 2D - slab geometry based on the work of J.P. Bizarro, H. Hugon and R. Jorge, 2018 (https://journals.aps.org/pra/abstract/10.1103/PhysRevA.98.023847)
- r1, r2 (2D - positions [m])
- k = ||k|| (scalar wavevector [m^(-1)])
- k1, k2 (2D - wavenumbers [m^(-1)])
- k_0 (unperturbed wavenumber [m^(-1)])
- c (speed of light [m.s^(-1)])
- q (wavenumber of turbulence [m^(-1)]
- n0 (unperturbed refractive index)
- n (refractive index)
- < n_e > (ensemble averaged density [m^(-3)])
- \phi (random phase [0, 2 Pi])
- \delta n_e (r) = \delta n_e0 cos(q r1 + \phi) (density fluctuations << < ne > [m^(-3)])
- n(r)/n0 = 1 + \delta n_e (r) / (turbulence profile through the effects of density fluctuations in the refractive index)
- t (time [s])
- w(r) = c k /n(r)
-
x[y] = (k_0 / 2 Pi) r1[2] (position)
-
\kappa_x[y] = k_1[2] / k_0
-
\delta n_e (r) = \delta n_e0 cos(q x + \phi)
-
q --> q/k_0 (in this normalized system, wavenumber of turbulence becomes adimensional)
-
\tau = (c k_0 / 2 Pi n_0) t
-
< ne > === 1 ---> \delta n_e0 (r) becomes a percentage
-
\kappa_x[y]_0, x[y]_0 refer to initial conditions for normalized wavevectors and position, respectively
-
Output example:
TODO:
- Seed random phases so the randomness doesn't get stuck: DONE
- Solve constructor issue/bug (randomly reassigning 'self' vars) : DONE ± (STILL NEED TO ADRESS)
- Review math (Monte Carlo and formalism equations) : DONE: solved major bugs
- Include info and logging (especially at the input handler)
- Include timer in the run() method
- Include possibility of plotting dxdx, dydy, dkxdkx, dkydky with Monte Carlo: DONE
- Include trajectory mode, in which the ray trajectory is depicted with the rms-spreading, and all the Monte Carlo Rays are plotted in this case (see paper): DONE
- Refactor output locations and naming
- Refactor plot title, font, etc. according to the input: DONE
- Raise error when input has non-valid MC parameters: DONE
- Plot 1-mode turbulent background: DONE
- Include turbulence eq. in the trajectory figure: DONE
- Refactor visualizer code... its a bit amateur... (last priority)