<div dir="ltr"><p>Dear FLASH Users,</p><p>I am using FLASH to simulate a cylindrical implosion for the National Ignition Facility (NIF). The simulations are performed in 2D r-z geometry.</p><p>Initially, I utilized the 3D-to-2D ray tracing algorithm, incorporating the 128 beams of the experimental configuration. I observed a rapid growth of instabilities at the gas-shell interface and initially considered that this might be a (potentially physical) result of azimuthal asymmetry (mode-8 is particularly expected) due to the beam imprint pattern projected in 2D. However, the asymmetries were far more pronounced than anticipated.</p><p>To rule out this effect, I simplified the simulations by switching to a 2D irradiation model using an azimuthally averaged "gaussian1D" function. Surprisingly, the implosions remained highly unstable, indicating that the projected azimuthal asymmetries were not the cause this time. I experimented with varying the number of rays, ticks, and grid resolution, but the instability persisted across all configurations.</p><p>At this point, I am uncertain about how to resolve this issue, and would greatly appreciate your advice and insights.</p><p>I have attached the <code>flash.par</code> and <code>setup_call</code> files for the simulation, which was conducted using FLASH version 4.8. Below, I include mass density plots ("dens" in g/cm³) at several time steps (2.5 ns, 5 ns, and 7.5 ns out of the 10 ns pulse duration) to illustrate the instabilities observed:</p><div><br></div><div><img src="cid:ii_m4x2ca9f0" alt="image.png" width="180" height="542" class="gmail-CToWUd gmail-a6T" tabindex="0" style="cursor: pointer; outline: 0px;"><img src="cid:ii_m4x2e3qu1" alt="image.png" width="222" height="542" class="gmail-CToWUd gmail-a6T" tabindex="0" style="cursor: pointer; outline: 0px;"><img src="cid:ii_m4x2fto92" alt="image.png" width="298" height="542" class="gmail-CToWUd gmail-a6T" tabindex="0" style="cursor: pointer; outline: 0px;"><br></div><div><br></div><p>Additionally, I enabled the MHD capability with an axial B-field, but it is currently set to B = 0 T in this simulation.</p><p>Thank you in advance for your time and assistance. I look forward to hearing your thoughts and potential suggestions.</p>Best regards,<br>Mathieu Bailly-Grandvaux</div>