<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Helvetica;
panose-1:0 0 0 0 0 0 0 0 0 0;}
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
span.apple-converted-space
{mso-style-name:apple-converted-space;}
span.EmailStyle21
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-US" link="blue" vlink="purple" style="word-wrap:break-word;-webkit-nbsp-mode: space;line-break:after-white-space">
<div class="WordSection1">
<p class="MsoNormal">Hi Adam/Eddie/All – <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">After spending some time with this today, I think I was overthinking this before and have the correct EOS set up for a 2T problem as is… However, to run the benchmark problem I am interested in, I have one additional challenge that needs
addressing. <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I need to set an analytic/user defined function for the specific heat, Cv. I have come to realize that this is retrieved by the EOS unit – e.g. in `path/to/Diffuse/DiffuseMain/Unsplit/diff_advanceTherm.F90`. I have been able to create custom
opacity functions, so I am somewhat familiar with creating additional functionalities within the FLASH source code, but it is not clear to me how one would define something new within the EOS unit.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Is there a way to create a user defined function for Cv?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thanks in advance,<o:p></o:p></p>
<p class="MsoNormal">Tony<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal" style="margin-bottom:12.0pt"><b><span style="font-size:12.0pt;color:black">From:
</span></b><span style="font-size:12.0pt;color:black">Alberti, Anthony L <anthony.alberti@oregonstate.edu><br>
<b>Date: </b>Friday, August 27, 2021 at 11:25 AM<br>
<b>To: </b>Reyes, Adam <adam.reyes@rochester.edu><br>
<b>Cc: </b>Hansen, Eddie <ehansen@pas.rochester.edu>, flash-users@flash.uchicago.edu <flash-users@flash.uchicago.edu><br>
<b>Subject: </b>Re: [FLASH-USERS] [EXT] Electron Energy Density, EELE_VAR, Initialization?<o:p></o:p></span></p>
</div>
<p class="MsoNormal">Hi Adam,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">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.
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">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.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">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?
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thanks,<o:p></o:p></p>
<p class="MsoNormal">Tony <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal" style="margin-bottom:12.0pt"><b><span style="font-size:12.0pt;color:black">From:
</span></b><span style="font-size:12.0pt;color:black">Reyes, Adam <adam.reyes@rochester.edu><br>
<b>Date: </b>Friday, August 27, 2021 at 11:01 AM<br>
<b>To: </b>Alberti, Anthony L <anthony.alberti@oregonstate.edu><br>
<b>Cc: </b>Hansen, Eddie <ehansen@pas.rochester.edu>, flash-users@flash.uchicago.edu <flash-users@flash.uchicago.edu><br>
<b>Subject: </b>Re: [FLASH-USERS] [EXT] Electron Energy Density, EELE_VAR, Initialization?<o:p></o:p></span></p>
</div>
<p><span style="color:#D73F09">[This email originated from outside of OSU. Use caution with links and attachments.]</span><o:p></o:p></p>
<div>
<p class="MsoNormal">Hi Tony, <o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">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 <a href="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">Table
16.1</a> from the users guide for the different options.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<div>
<p class="MsoNormal"> eosModeInit [STRING] ["dens_ie"]<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> Valid Values: "dens_ie", "dens_pres", "dens_temp", "dens_ie_all",<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> "dens_ie_scatter", "dens_ie_gather", "dens_temp_equi", "dens_temp_all",<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> "dens_temp_gather", "eos_nop"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> determines how to operate on thermodynamic quantities for the initial<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> conditions. A call to Simulation_initBlock sets initial conditions on<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> each block, and an eos call follows which ensures the initial values are<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> thermodynamically consistent. The value of eosModeInit determines the<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> mode in which these eos calls operate. Possible values are "dens_ie",<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> where density and internal energy are supplied and pressure and<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> temperature are returned; "dens_pres", where density and pressure are<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> given and energy and temperature are computed; and "dens_temp", where<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> density and temperature are given and pressure and energy are computed. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> Other values may also be available as, depending on whether a multiTemp<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> Eos is used and on implementation of additional physics.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Adam<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><o:p> </o:p></p>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<p class="MsoNormal">On Aug 27, 2021, at 1:47 PM, Alberti, Anthony L <<a href="mailto:anthony.alberti@oregonstate.edu">anthony.alberti@oregonstate.edu</a>> wrote:<o:p></o:p></p>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">Ah, the dreaded EOS unit. I figured that might be the case.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">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?<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Thanks, Eddie.<span class="apple-converted-space"> </span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Tony<span class="apple-converted-space"> </span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal" style="margin-bottom:12.0pt"><b><span style="font-size:12.0pt">From:<span class="apple-converted-space"> </span></span></b><span style="font-size:12.0pt">Eddie Hansen <<a href="mailto:ehansen@pas.rochester.edu">ehansen@pas.rochester.edu</a>><br>
<b>Date:<span class="apple-converted-space"> </span></b>Friday, August 27, 2021 at 9:15 AM<br>
<b>To:<span class="apple-converted-space"> </span></b>Alberti, Anthony L <<a href="mailto:anthony.alberti@oregonstate.edu">anthony.alberti@oregonstate.edu</a>><br>
<b>Cc:<span class="apple-converted-space"> </span></b><a href="mailto:flash-users@flash.uchicago.edu">flash-users@flash.uchicago.edu</a><span class="apple-converted-space"> </span><<a href="mailto:flash-users@flash.uchicago.edu">flash-users@flash.uchicago.edu</a>><br>
<b>Subject:<span class="apple-converted-space"> </span></b>Re: [EXT] [FLASH-USERS] Electron Energy Density, EELE_VAR, Initialization?</span><o:p></o:p></p>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:9.0pt;font-family:Helvetica;color:#D73F09">[This email originated from outside of OSU. Use caution with links and attachments.]</span><span style="font-size:9.0pt;font-family:Helvetica"><o:p></o:p></span></p>
<div>
<div>
<div>
<p class="MsoNormal">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.<o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<div>
<div>
<p class="MsoNormal">On Fri, Aug 27, 2021 at 11:05 AM Alberti, Anthony L <<a href="mailto:anthony.alberti@oregonstate.edu">anthony.alberti@oregonstate.edu</a>> wrote:<o:p></o:p></p>
</div>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0in;margin-bottom:5.0pt">
<div>
<div>
<div>
<p class="MsoNormal">Hi folks,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">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?<span class="apple-converted-space"> </span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">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.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Thanks all.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Tony Alberti<o:p></o:p></p>
</div>
</div>
</div>
</blockquote>
</div>
<div>
<p class="MsoNormal"><br clear="all">
<o:p></o:p></p>
</div>
<div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal">--<span class="apple-converted-space"> </span><o:p></o:p></p>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal">Eddie Hansen, PhD<br>
Postdoctoral Associate<o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">University of Rochester<o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">607-341-6126 | Flash Center<o:p></o:p></p>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
</div>
</body>
</html>