[FLASH-USERS] 2D Spherical and Cartesian Geometry

[Michiel Bustraan]

Klaus,

Those changes didn't seem to help, but switching to split hydro did seem to make the simulation work.
However, when I tried to increase the resolution, either through more AMR levels or a higher base resolution, I get the following error message.

Missing neighbors on PE    1, lb=  243,surrblks=( 232,   1, 2)( 237,   1, 1)( 249,   1, 1)( 238,   1, 2)( 243,   1, 2)( 251,   1, 1)( 254,   1, 1)( 255,   1, 1)(  -1,  -1,-1)
Missing neighbors on PE    1, lb=  243,surrblks=( 232,   1, 2)( 237,   1, 1)( 249,   1, 1)( 238,   1, 2)( 243,   1, 2)( 251,   1, 1)( 254,   1, 1)( 255,   1, 1)(  -1,  -1,-1)

Do you know what might cause this?
Also, I was unable to run the simulation with the 'dens_temp' EOS. Is split hydro only able to run with 'dens_ie', or is the temperature EOS supposed to work?

- Michiel
________________________________________
From: Klaus Weide <klaus at flash.uchicago.edu>
Sent: Friday, September 23, 2016 9:17 PM
To: Michiel Bustraan
Cc: flash-users at flash.uchicago.edu
Subject: Re: [FLASH-USERS] 2D Spherical and Cartesian Geometry

Michiel,

Your message makes it difficult to quote individual lines (it had some
strange line-end encoding), so I limit myself to the following comments:

The following looks unusual - try without it:
 EOSforRiemann   = .true.

What Eos are you using; does Eos/EosMain/Gamma work?

Also try without the following;
  hy_fPresInMomFlux = 0.0

Re
  order           = 5          # Interpolation order (first/second/third/fifth order)

Try 2 or 3.


The problems seem to occur around 180 deg, not around 0 deg;
so try a Y range of 0 - 90 deg (corresponding to a half sphere)
with reflecting Y b.c. at the equator.

I am also curious whether split Hydro (PPM, no +uhd) behaves better.

Klaus