<div dir="ltr"><div>Hi Mathieu,</div><div><br></div><div>To my untrained eye, I don't notice any obvious grid artefacts, so it's not clear to me the instability you observe is not unphysical. If you're trying to do a simulation without magnetic fields, you might also try setting ForceHydroLimit = .true. in your flash.par. Visualizing the pressure might be helpful as well.</div><div><br></div><div>Best wishes,</div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr">--------<div>Ryan</div></div></div></div></div></div></div><br></div><br><div class="gmail_quote gmail_quote_container"><div dir="ltr" class="gmail_attr">On Fri, Dec 20, 2024 at 11:13 AM Mathieu Bailly-Grandvaux <<a href="mailto:mbaillygrandvaux@ucsd.edu">mbaillygrandvaux@ucsd.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><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" style="outline: 0px;"><img src="cid:ii_m4x2e3qu1" alt="image.png" width="222" height="542" style="outline: 0px;"><img src="cid:ii_m4x2fto92" alt="image.png" width="298" height="542" style="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>
_______________________________________________<br>
flash-users mailing list<br>
<a href="mailto:flash-users@flash.rochester.edu" target="_blank">flash-users@flash.rochester.edu</a><br>
<br>
For list info, including unsubscribe:<br>
<a href="https://flash.rochester.edu/mailman/listinfo/flash-users" rel="noreferrer" target="_blank">https://flash.rochester.edu/mailman/listinfo/flash-users</a><br>
</blockquote></div>