<div dir="ltr">Hi Gabriel,<div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">By the way I test, the "SodSpherical" setup command written in the flash<br>documentation:<br>"./setup SodSpherical -auto -2d -geometry=spherical"</blockquote><div>You'll want to add "+splitHydro" since the default unsplit hydro solver does not support spherical geometry.</div><div><br></div><div>Adam </div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Feb 17, 2021 at 8:10 PM RIGON Gabriel <<a href="mailto:gabriel.rigon@nagoya-u.jp">gabriel.rigon@nagoya-u.jp</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">Hi,<br>
<br>
Thank you for your answer.<br>
So If I understand it, you can't use "periodic" boundary conditions in <br>
the angular direction (y), is that right?<br>
<br>
I tried the "SodSpherical" as suggest, using the<br>
<br>
By the way I test, the "SodSpherical" setup command written in the flash <br>
documentation:<br>
"./setup SodSpherical -auto -2d -geometry=spherical"<br>
<br>
While the simulation can be initialized normally, it does not give any <br>
results:<br>
The first checkpoint (sod_sph_hdf5_chk_0000) contains the expected data, <br>
but starting from the second one (0001) the data are all wrong, or more <br>
precisely they don't exist anymore ("dens" is -none- everywhere).<br>
<br>
I'm currently using FLASH4.6.2. I suppose something is wrong with either <br>
my installation of FLASH or with this version of FLASH.<br>
Could someone try the "SodSpherical" problem and tell me if they obtain <br>
the same results.<br>
<br>
Best regards,<br>
Gabriel<br>
<br>
On 18/02/2021 08:57, Klaus Weide wrote:<br>
> On Wed, 17 Feb 2021, RIGON Gabriel wrote:<br>
><br>
>> I'm trying to setup a 2D spherical simulation, but I have some trouble with<br>
>> it.<br>
>> Do anyone know a way to initialize a 2D spherical simulation ?<br>
> This should work:<br>
><br>
>> - ymin = 0.0 (angle in degree)<br>
>> - ymax = 180.0 (angle in degree)<br>
>> - yl/r_boundary_condition = "reflecting"<br>
> So should this:<br>
><br>
>> - ymin = 0.0 (angle in degree)<br>
>> - ymax = 90.0 (angle in degree)<br>
>> - yl/r_boundary_condition = "reflecting"<br>
> Or this:<br>
><br>
>> - ymin = 0.0 (angle in degree)<br>
>> - ymax = 90.0 (angle in degree)<br>
>> - yl_boundary_condition = "reflecting"<br>
>> - yr_boundary_condition = "outflow"<br>
><br>
> You may want to look at the separate "SodSpherical" simulation and compare<br>
> it to "Sod".<br>
><br>
> Klaus<br>
<br>
-- <br>
Gabriel Rigon<br>
Laboratory for Theoretical Astronomy & Astrophysics (TA-Lab)<br>
Graduate School of Science<br>
Division of Particle and Astrophysical Science<br>
Nagoya University, Japan<br>
<br>
</blockquote></div>