[FLASH-USERS] flash-users Digest, Vol 138, Issue 7

Thu Apr 18 06:01:04 EDT 2019

Hi, Michel:
Thank you!
I used the tabulated EOS provided by FLASH center in my simulation. If the problem comes from EOS, should I change the EOS?
Best regards,
Wenyi

From: flash-users-bounces at flash.uchicago.edu [mailto:flash-users-bounces at flash.uchicago.edu] On Behalf Of Michel de Peretti
Sent: Wednesday, April 17, 2019 8:00 PM
To: flash-users at flash.uchicago.edu
Subject: Re: [FLASH-USERS] flash-users Digest, Vol 138, Issue 7

I am not fully familiar with Flash but very often this kind of problem comes from your equation of state....try to check it.
Hoping that helps you
Michel

Le mer. 17 avr. 2019 à 19:00, <flash-users-request at flash.uchicago.edu<mailto:flash-users-request at flash.uchicago.edu>> a écrit :
Send flash-users mailing list submissions to
        flash-users at flash.uchicago.edu<mailto:flash-users at flash.uchicago.edu>

To subscribe or unsubscribe via the World Wide Web, visit
        http://flash.uchicago.edu/mailman/listinfo/flash-users
or, via email, send a message with subject or body 'help' to
        flash-users-request at flash.uchicago.edu<mailto:flash-users-request at flash.uchicago.edu>

You can reach the person managing the list at
        flash-users-owner at flash.uchicago.edu<mailto:flash-users-owner at flash.uchicago.edu>

When replying, please edit your Subject line so it is more specific
than "Re: Contents of flash-users digest..."

Today's Topics:

   1.  The problem of Negative Internal Energy in 3T MHD simulation
      with laser generated plasmas (Huo Wenyi)

----------------------------------------------------------------------

Message: 1
Date: Wed, 17 Apr 2019 07:43:27 +0000
From: Huo Wenyi <Wenyi.Huo at eli-beams.eu<mailto:Wenyi.Huo at eli-beams.eu>>
Subject: [FLASH-USERS] The problem of Negative Internal Energy in 3T
        MHD simulation with laser generated plasmas
To: "flash-users at flash.uchicago.edu<mailto:flash-users at flash.uchicago.edu>" <flash-users at flash.uchicago.edu<mailto:flash-users at flash.uchicago.edu>>
Message-ID: <4b8b49f0356446ab9b4326bb39c14550 at eli-beams.eu<mailto:4b8b49f0356446ab9b4326bb39c14550 at eli-beams.eu>>
Content-Type: text/plain; charset="us-ascii"

Dear all:

I am using FLASH to simulate the magnetohydrodynamics in  laser generated plasmas. However, I always encounter a problem of negative internal energy both in FLASH 4.5 and FLASH4.6. Does anyone encounter such a similar problem?  The error information in the Log file is as follows:

[ 04-16-2019  14:34:35.488 ] [hy_uhd_ragelike] ERROR DETECTED
[ 04-16-2019  14:34:35.488 ]   Zero or negative advected specific internal energy detected
[ 04-16-2019  14:34:35.488 ]   DENS_OLD =  1.0000904290291E-09
[ 04-16-2019  14:34:35.488 ]   DENS_NEW =  1.0000983091349E-09
[ 04-16-2019  14:34:35.488 ]   UTOT_NEW = -3.7386802822434E-01
[ 04-16-2019  14:34:35.488 ]   UEXTRA = -1.7491864582091E+01
[ 04-16-2019  14:34:35.488 ]     UELE_NEW = -3.7135596091566E-01
[ 04-16-2019  14:34:35.488 ]     UION_NEW = -2.5120492887951E-03
[ 04-16-2019  14:34:35.488 ]     URAD_NEW = -1.8019885718467E-08
[ 04-16-2019  14:34:35.488 ]     EELE_VAR =  1.7001438966846E+10
[ 04-16-2019  14:34:35.488 ]     EION_VAR =  1.1498940630502E+08
[ 04-16-2019  14:34:35.488 ]     ERAD_VAR =  8.2498749422145E+02
[ 04-16-2019  14:34:35.488 ]     EELE     = -3.7131945682109E+08
[ 04-16-2019  14:34:35.488 ]     EION     = -2.5118023556785E+06
[ 04-16-2019  14:34:35.488 ]     ERAD     = -1.8018114373230E+01
[ 04-16-2019  14:34:35.488 ]     uele_adv  =  1.7002978536045E+01
[ 04-16-2019  14:34:35.488 ]     uion_adv  =  1.1501719275908E-01
[ 04-16-2019  14:34:35.488 ]     urad_adv  =  8.2506210304961E-07
[ 04-16-2019  14:34:35.488 ]   PELE_VAR =  1.2145270532940E+00
[ 04-16-2019  14:34:35.488 ]   PION_VAR =  7.6660475248193E-02
[ 04-16-2019  14:34:35.488 ]   PRAD_VAR =  2.7502069901319E-07
[ 04-16-2019  14:34:35.488 ]   PELE_ADV =  1.2145343502140E+00
[ 04-16-2019  14:34:35.488 ]   PION_ADV =  7.6672066071296E-02
[ 04-16-2019  14:34:35.488 ]   PRAD_ADV =  2.7502070101654E-07
[ 04-16-2019  14:34:35.488 ]   TELE_VAR =  6.0996867323016E+02
[ 04-16-2019  14:34:35.488 ]   TION_VAR =  3.6907447766733E+00
[ 04-16-2019  14:34:35.489 ]   TRAD_VAR =  1.0219015193890E+02
[ 04-16-2019  14:34:35.489 ]   CELL X =  3.3222656250000E-01
[ 04-16-2019  14:34:35.489 ]   CELL Y =  6.9890625000000E-01
[ 04-16-2019  14:34:35.489 ]   CELL Z =  0.0000000000000E+00
[ 04-16-2019  14:34:35.489 ]   You might be able to get things to keep running by trying
[ 04-16-2019  14:34:35.489 ]   more diffusive hydro options (reduce CFL, order, etc...)
[ 04-16-2019  22:00:17.496 ] [hy_uhd_ragelike] Negative internal energies in 3T update
[ 04-16-2019  22:00:17.496 ]   Falling back to using internal energy ratios instead of pressure ratios

I have tried many different RiemannSolvers. It does not help. I guess that this problem is related to the grid size, but I am not sure.

Can anyone provide a solution to this problem? Thank you very much!

Best regards,

Wenyi Huo

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://flash.uchicago.edu/pipermail/flash-users/attachments/20190417/a66bdae6/attachment-0001.html

------------------------------

_______________________________________________
flash-users mailing list
flash-users at flash.uchicago.edu<mailto:flash-users at flash.uchicago.edu>
http://flash.uchicago.edu/mailman/listinfo/flash-users

End of flash-users Digest, Vol 138, Issue 7
*******************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20190418/006d460b/attachment.htm>