[FLASH-USERS] Numerical Effects in Spherical Coordinates

[Jason Galyardt]

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> 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/cc092e49/attachment.htm>