[FLASH-USERS] energy inaccuracy in Sedov simulation

Marissa Adams madams at pas.rochester.edu
Tue Jul 26 11:25:19 EDT 2016


We understand that. The total energy is less than the input value. This is
the problem.

On Tue, Jul 26, 2016 at 11:23 AM, Tomasz Plewa <tplewa at fsu.edu> wrote:

> The relevant data is the total energy shown in the first line of the dat
> file. That should stay constant with relative accuracy better than 1e-10.
>
> The actual initial value of the total will vary as one changes problem
> parameters.
>
> Tomek
> --
>
> On 07/26/16 10:30, Slavin, Jonathan wrote:
>
> Hi all,
>
> I've been running supernova remnant evolution simulations.  The
> initialization is done with code taken from the supplied Sedov simulation
> with certain small changes.  One issue that I've run into is that the total
> energy is not very accurate.  That is, though the energy is conserved, the
> total energy, initially input as thermal energy, is below that which is
> given in the sim_expEnergy parameter by ~5%.  Has anyone else run into
> this?  I don't think that it has to do with any of the changes that I've
> made to Simulation_initBlock or Simulation_init.  I'm wondering if
> increasing the sim_nsubzones parameter might help.  Any advice would be
> appreciated.
>
> Thanks,
> Jon
>
> --
> ________________________________________________________
> Jonathan D. Slavin                 Harvard-Smithsonian CfA
> jslavin at cfa.harvard.edu       60 Garden Street, MS 83
> phone: (617) 496-7981       Cambridge, MA 02138-1516
> cell: (781) 363-0035             USA
> ________________________________________________________
>
>
>


-- 
Marissa Adams                         E-mail: madams at pas.rochester.edu
Graduate Student                      Ph: (585) 402-5779
Department of Physics & Astronomy     Website:
http://www.pas.rochester.edu/~madams/
University of Rochester
478 Bausch & Lomb Hall
P.O. Box 270171
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20160726/0152e250/attachment.htm>


More information about the flash-users mailing list