[FLASH-USERS] [EXT] Electron Energy Density, EELE_VAR, Initialization?

[Alberti, Anthony L]

Hi Adam,

Ok, I think this makes sense, but I’d like to run something by you (and the group) to make sure that I understand things clearly.

I’m interested in running a 2T calculation, where the radiation and electron energy densities and temperatures are updated on each time step (the electron and ion temperatures are assumed to be the same). From Table 16.1, it seems as though “eosModeInit = dens_ie_gather” would be the best option for the initial condition and the runtime variable “eosMode” should match and be set to “dens_ie_gather” as well.

Am I on the right track here? If this is indeed right, Simulation_initBlock will need to change as well, since I would be supplying energy densities as runtime variables rather than temperatures. Correct?

Thanks,
Tony

From: Reyes, Adam <adam.reyes at rochester.edu>
Date: Friday, August 27, 2021 at 11:01 AM
To: Alberti, Anthony L <anthony.alberti at oregonstate.edu>
Cc: Hansen, Eddie <ehansen at pas.rochester.edu>, flash-users at flash.uchicago.edu <flash-users at flash.uchicago.edu>
Subject: Re: [FLASH-USERS] [EXT] Electron Energy Density, EELE_VAR, Initialization?

[This email originated from outside of OSU. Use caution with links and attachments.]
Hi Tony,

The EOS gets called during the grid initialization after Simulation_initBlock, you can look for instances of “Eos_wrapped”. The purpose is to ensure thermodynamic consistency between temperature pressure & internal energy in the initial conditions. How that is done depends on the choice of the Eos mode, controlled by the runtime parameter “eosModeInit”. You can look at Table 16.1<https://nam04.safelinks.protection.outlook.com/?url=http%3A%2F%2Fflash.uchicago.edu%2Fsite%2Fflashcode%2Fuser_support%2Fflash4_ug_4p62%2Fnode110.html%23Tab%3A3TEos_MODES&data=04%7C01%7Canthony.alberti%40oregonstate.edu%7C5272240cf9d84ca5d0d608d96984989e%7Cce6d05e13c5e4d6287a84c4a2713c113%7C0%7C0%7C637656840935112747%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&sdata=DQU5Fe9wEZOZCNEjfjbRX277E3tUasVFhuoPATjns80%3D&reserved=0> from the users guide for the different options.

    eosModeInit [STRING] ["dens_ie"]
        Valid Values: "dens_ie", "dens_pres", "dens_temp", "dens_ie_all",
        "dens_ie_scatter", "dens_ie_gather", "dens_temp_equi", "dens_temp_all",
        "dens_temp_gather", "eos_nop"
        determines how to operate on thermodynamic quantities for the initial
        conditions.  A call to Simulation_initBlock sets initial conditions on
        each block, and an eos call follows which ensures the initial values are
        thermodynamically consistent.  The value of eosModeInit determines the
        mode in which these eos calls operate.  Possible values are "dens_ie",
        where density and internal energy are supplied and pressure and
        temperature are returned; "dens_pres", where density and pressure are
        given and energy and temperature are computed; and "dens_temp", where
        density and temperature are given and pressure and energy are computed.
        Other values may also be available as, depending on whether a multiTemp
        Eos is used and on implementation of additional physics.

Adam


On Aug 27, 2021, at 1:47 PM, Alberti, Anthony L <anthony.alberti at oregonstate.edu<mailto:anthony.alberti at oregonstate.edu>> wrote:

Ah, the dreaded EOS unit. I figured that might be the case.

However, I’m not sure where the EOS unit gets called in my Simulation_initBlock.F90 file (see attached). This is based off of the “SimulationMain/MGDStep” problem that comes with the FLASH distribution. Would you be able to point me to where that unit gets called?

Thanks, Eddie.

Tony

From: Eddie Hansen <ehansen at pas.rochester.edu<mailto:ehansen at pas.rochester.edu>>
Date: Friday, August 27, 2021 at 9:15 AM
To: Alberti, Anthony L <anthony.alberti at oregonstate.edu<mailto:anthony.alberti at oregonstate.edu>>
Cc: 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>>
Subject: Re: [EXT] [FLASH-USERS] Electron Energy Density, EELE_VAR, Initialization?
[This email originated from outside of OSU. Use caution with links and attachments.]
After Simulation_initBlock, there should be a call to EOS to initialize EELE_VAR and all other relevant variables that haven't yet been set. The EOS unit source code can get a little confusing, but I hope at least pointing you in the right direction is helpful.

On Fri, Aug 27, 2021 at 11:05 AM Alberti, Anthony L <anthony.alberti at oregonstate.edu<mailto:anthony.alberti at oregonstate.edu>> wrote:
Hi folks,

I’m currently working on running a radiation transport benchmark problem in FLASH and am interested to know, where in the FLASH source code does the electron energy density, EELE_VAR, get initialized?

I can see where the radiation energy density gets initialized, Simulation_blockInit.F90, but it’s not clear where the electron energy density gets initialized. I’ve grepped around and can’t find anything obvious.

Thanks all.
Tony Alberti


--
Eddie Hansen, PhD
Postdoctoral Associate
University of Rochester
607-341-6126 | Flash Center

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20210827/69428f85/attachment-0001.htm>