Model Description » Historique » Révision 4
Révision 3 (Florian Ferreira, 21/08/2014 12:26) → Révision 4/11 (Florian Ferreira, 21/08/2014 12:27)
h1. Model Description The content of YoGA_Ao has been validated at several levels. This page outlines the main results. [[Model_Description#Turbulence-generation|Turbulence generation]] * [[Model_Description#Algorithm-validation-on-the-CPU|Algorithm validation on the CPU]] * [[Model_Description#Multi-layer-validation-on-the-GPU|GPU version and raytracing through multiple layers]] [[Model_Description#Wavefront-sensor-model|Wavefront sensor model]] * [[Model_Description#General-model|General model]] * [[Model_Description#LGS-spot-model|LGS spot model]] * [[Model_Description#LGS-cone-effect|LGS cone effect]] h2. Turbulence generation In YoGA Ao, turbulence generation is done through extrusion of Kolmogorov-type ribbons to create an infinite phase screen (based on the work of Assémat et al, Opt. Express 14, 988–999 (2006) and Fried & Clark, J. Opt. Soc. Am. A 25, 463-468 (2008)) This process allows to generate infinite length phase screens on reasonably-sized supports. The generated screens have an infinite external scale. h3. Algorithm validation on the CPU The algorithm has been first validated on the CPU. The following figure demonstrates the validity of the model by comparing the phase structure function (computed as the average of the square of each phase screen realization minus its value at pixel [1,1]) obtained on phase screens with various sizes and the theoretical value of 6.88 (r/r0)5/3. !https://dev-lesia.obspm.fr/projets/attachments/download/571/extrudeCPU.png! legend: black: theoretical, red: 64x64 phase screen, blue: 80x80, green: 128x128 and yellow: 256x256. In these simulations, the total r0 has been taken as 0.16m, a typical value at 0.5µm for a professional astronomical site. A minimum of 100,000 column extrusion is required to get a good convergence. Even with such a number of iterations, the figure above shows that the larger the screen is the slower is convergence. h3. Multi-layer validation on the GPU The algorithm has then been validated on the GPU. YoGA_Ao provides the ability to simulate several turbulent layers at various altitude and with various r0 and wind speeds and directions. On the GPU, the overall model includes phase screen generation and raytracing through each layer along a specific direction. In this configuration, two separate targets were simulated, located at two different positions on the sky. To accommodate for the induced field of view, non-zero altitude screens have to be proportionately larger than the ground layer screen. The following figures demonstrates the validity of this model for different atmosphere configuration (for 1 layer, 4 layers, 12 layers). As a reference, the single-layer CPU result for the same total r0 (i.e. as seen from the ground) is plotted in red, the GPU result is in green. !https://dev-lesia.obspm.fr/projets/attachments/download/572/extrudeGPU.png! legend: left: one layer, center: 4 layers, right: 12 layers The discrepancies in the multi-layer cases and the theoretical structure function can at least partly be attributed to the presence of larger phase screens at altitudes > 0. As seen previously larger screens are slower to converge. Moreover this effect is only visible for low spatial frequencies which are as well the slowest to converge. h2. Wavefront sensor model h3. General model h3. LGS spot model h3. LGS cone effect The cone effect, due to the finite height of the sodium layer where the LGS is produced, is taking into account in the simulation thanks to a GPU parallelized raytracing algorithm. For each pixel of the phase screen seen by the WFS, we calculate the sum of the corresponding pixels on each phase screen of the atmosphere layers. Each term is weighting by a the Cn^2 coefficient of its atmosphere layer. A linear interpolation is performed between the 4-near pixels of the atmosphere layer seen by the pixel of the WFS phase screen. Basically, this is equivalent to do a "zoom" on each phase screen and summating them. !https://projets-lesia.obspm.fr/attachments/download/992/Target.png! !https://projets-lesia.obspm.fr/attachments/download/991/Cone_effect.png! *legend :* _first _left : atmosphere layer phase screen at 20km ; second right : corresponding phase screen seen by the WFS with the LGS at 90km_ The two figures above illustrate it : the LGS is at 90km and the atmosphere layer at 20km, so we can observe that the phase screen seen by the WFS (first (right figure) is a centered zoom of the atmosphere phase screen (second (left figure). Performances coming soon ...