[FLASH-USERS] Numerical Effects in Spherical Coordinates

[Sean M. Couch]

Also, it would be helpful if you shared the .log file from this run so that all built-in units and runtime parameter values can be seen.

Cheers,
Sean


----------------------------------------------------------------------------------------------------------
Sean M. Couch
Assistant Professor
Department of Physics and Astronomy
Department of Computational Mathematics, Science, and Engineering
National Superconducting Cyclotron Laboratory/Facility for Rare Isotope Beams
Michigan State University
567 Wilson Rd, 3250 BPS
East Lansing, MI 48824
(517) 884-5035 —— couch at pa.msu.edu —— www.pa.msu.edu/~couch

On Mar 23, 2017, 12:37 PM -0400, Jason Galyardt <jason.galyardt at gmail.com>, wrote:
Hi Michiel,

If you suspect that your temperature fluctuations are part of the problem, try removing them. You mentioned that you need high resolution; in order to get full refinement, you can set the minimum refinement level equal to the maximum refinement level:

lrefine_min = 10
lrefine_max = 10

Of course, this effectively turns off the AMR feature, so you may have to add processors in order to get the resolution you need across the entire domain. You might also try reducing your domain dimensions for this test run, in order to reduce the cluster resource requirements.

Also, you mention that the gas might be bouncing off the inner boundary: are your boundary conditions reflective, then? If not, there will be no bouncing; if so, turn the boundary conditions to 'outflow' and see what happens.

Cheers,
Jason

------
Jason Galyardt
University of Georgia


On Wed, Mar 22, 2017 at 9:10 AM, Michiel Bustraan <michiel.bustraan at astro.su.se<mailto:michiel.bustraan at astro.su.se>> wrote:

Dear FLASH Users,


I've been running 2D simulations of shocks in spherical simulations.

The initial conditions involve a high density around the origin.


I had to introduce random temperature fluctuations in the inner region to ensure refinement.

I started with quite a high level of refinement, and the maximum level of refinement decreases over time (from 10 layers to 5).


I've found that ripples appear in the inner portion of the domain, as well as larger numerical artifacts.

I think that these might either be due to:


- The initial temperature fluctuations I added, though I feel these would be too small (spatially) to cause all these issues.

- Could it be that the higher pressure region in the center, combined with the lower radial velocity, cause the gas to bounce off the inner boundary?

This might create ripples and internal shocks.


Either way, I'm not certain how to improve my results from here.

If anyone has had similar problems, I would greatly appreciate some advice.


The figures attached highlight my problem.


Kind Regards,

Michiel Bustraan

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20170323/5508e617/attachment-0001.htm>