[FLASH-USERS] setup 2D spherical simulation

acreyes acreyes acreyes at ucsc.edu
Thu Feb 18 09:12:05 EST 2021


Hi Gabriel,

By the way I test, the "SodSpherical" setup command written in the flash
> documentation:
> "./setup SodSpherical -auto -2d -geometry=spherical"

You'll want to add "+splitHydro" since the default unsplit hydro solver
does not support spherical geometry.

Adam

On Wed, Feb 17, 2021 at 8:10 PM RIGON Gabriel <gabriel.rigon at nagoya-u.jp>
wrote:

> Hi,
>
> Thank you for your answer.
> So If I understand it, you can't use "periodic" boundary conditions in
> the angular direction (y), is that right?
>
> I tried the "SodSpherical" as suggest, using the
>
> By the way I test, the "SodSpherical" setup command written in the flash
> documentation:
> "./setup SodSpherical -auto -2d -geometry=spherical"
>
> While the simulation can be initialized normally, it does not give any
> results:
> The first checkpoint (sod_sph_hdf5_chk_0000) contains the expected data,
> but starting from the second one (0001) the data are all wrong, or more
> precisely they don't exist anymore ("dens" is -none- everywhere).
>
> I'm currently using FLASH4.6.2. I suppose something is wrong with either
> my installation of FLASH or with this version of FLASH.
> Could someone try the "SodSpherical" problem and tell me if they obtain
> the same results.
>
> Best regards,
> Gabriel
>
> On 18/02/2021 08:57, Klaus Weide wrote:
> > On Wed, 17 Feb 2021, RIGON Gabriel wrote:
> >
> >> I'm trying to setup a 2D spherical simulation, but I have some trouble
> with
> >> it.
> >> Do anyone know a way to initialize a 2D spherical simulation ?
> > This should work:
> >
> >> - ymin = 0.0 (angle in degree)
> >> - ymax = 180.0 (angle in degree)
> >> - yl/r_boundary_condition = "reflecting"
> > So should this:
> >
> >> - ymin = 0.0 (angle in degree)
> >> - ymax = 90.0 (angle in degree)
> >> - yl/r_boundary_condition = "reflecting"
> > Or this:
> >
> >> - ymin = 0.0 (angle in degree)
> >> - ymax = 90.0 (angle in degree)
> >> - yl_boundary_condition = "reflecting"
> >> - yr_boundary_condition = "outflow"
> >
> > You may want to look at the separate "SodSpherical" simulation and
> compare
> > it to "Sod".
> >
> > Klaus
>
> --
> Gabriel Rigon
> Laboratory for Theoretical Astronomy & Astrophysics (TA-Lab)
> Graduate School of Science
> Division of Particle and Astrophysical Science
> Nagoya University, Japan
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20210218/6ed7c204/attachment.htm>


More information about the flash-users mailing list