<div dir="ltr">Dear Marissa, Klaus and <span class="" id=":1fi.1" tabindex="-1" style="">Hanif</span>,<div>Many thanks for your inputs - I really appreciate them! I do understand that my problem goes beyond what FLASH currently does and therefore requires significant new coding. Coding tabulated data is just a technical problem whereas Klaus' question "How to combine "two different <span class="" id=":1fi.2" tabindex="-1" style="">EoS</span>" in cells that contain both species?" is of essential importance. I am trying to find some publications on this subject. Currently I think that I can find partial pressures of both species from separate <span class="" id=":1fi.3" tabindex="-1" style="">EoS</span>. Densities of both species are already updated by separate equations of continuity. Velocities are updated by one momentum conservation equation which is linear in pressure so I just need to add up all partial pressures to find "averaged" velocity. Energy conservation equation is linear in density and pressure so again I need to use sums of densities and pressures over all species to update "averaged" energy (total or internal). This averaged internal energy should be used in <span class="" id=":1fi.4" tabindex="-1" style="">EoS</span> for each species along with corresponding density to find partial pressures of each species. Please let me know your opinion.</div><div>Thank you very much!</div><div><span class="" id=":1fi.5" tabindex="-1" style="">Maksim</span> <span class="" id=":1fi.6" tabindex="-1" style="">Kozlov</span></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Feb 19, 2021 at 11:30 PM M. B. P. Adams <<a href="mailto:madams15@ur.rochester.edu">madams15@ur.rochester.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="overflow-wrap: break-word;">Hi all,<div><br></div><div>This is a good yet potentially controversial question as from what I understand, and perhaps Petros and others would be willing to chime in here. I may also be misunderstanding the request here.</div><div><br></div><div>Effectively FLASH does not support the self-production of “tabulated” EOS. Tables are something you’ll have to find resources for elsewhere, sadly. You can <font face="Courier">grep</font> within the contents of the <i>source/</i> (much of these things can be found in the <i>Simulation</i> modules provided) for IONMIX4 files. For instance you’ll find some in the <i>LaserSlab</i> problem module. They end with <i>*.cnr</i> or <i>*.cn4</i>. FLASH also supports PrOpacEOS* tables (ends with <i>*.prp</i>) however they are not provided within the source and that requires a patch to be applied to adjust the relevant routines. </div><div><br></div><div>I hope that helps somewhat. I’d recommend building your simulation and utilizing the provided IONMIX4 tables you can find in the <i>source/</i> (even if they are not for your desired material, such as Cu) to make sure that your initialization and everything runs smoothly while in tandem working to obtain tables.</div><div><br></div><div>I’d avoid asking users to share their tables though as that may put them in an awkward position as I think many of such things, especially for higher Z materials, are export controlled. When it comes to that I’d directly refer to the Director of the FLASH Center and stuff.</div><div><br></div><div>Best,</div><div>Marissa</div><div><br></div><div>* <a href="http://www.prism-cs.com/Software/Propaceos/overview.html" target="_blank">http://www.prism-cs.com/Software/Propaceos/overview.html</a><br><div><br><blockquote type="cite"><div>On Feb 19, 2021, at 12:19, Maksim Kozlov <<a href="mailto:maksim.kozlov@nu.edu.kz" target="_blank">maksim.kozlov@nu.edu.kz</a>> wrote:</div><br><div><div dir="ltr">Thank you Hanif, but as I explained in my question it does not work for me because FLASH solves one equation of state with mass averaged adiabatic index and I need two different EoS for each species.</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Feb 19, 2021 at 10:47 PM Hanif zargarnezhad <<a href="mailto:h.zargarnejad@gmail.com" target="_blank">h.zargarnejad@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr">Hello Maksim,<div>If I understand your question correctly, you can use different approaches to do it. The simple way is that using Simulation_initSpecies.F90 you can define constant gas properties for different species and call them to EOS  wherever is located.</div><div>regards,</div><div><br></div><div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">Hanif Zargar</div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px">Fluid Mixing at Extreme conditions Laboratory</div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;font-size:12pt;line-height:inherit;vertical-align:baseline">Texas A&M Dept. of Mechanical  Engineering</span></div><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;margin-top:0px;margin-bottom:0px"><span style="margin:0px;padding:0px;border:0px;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;font-size:12pt;line-height:inherit;vertical-align:baseline">Graduate Research Assistant </span></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Feb 19, 2021 at 6:35 AM Maksim Kozlov <<a href="mailto:maksim.kozlov@nu.edu.kz" target="_blank">maksim.kozlov@nu.edu.kz</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Dear fiends,<div> I want to simulate hydrodynamic explosions of one material (say copper) into another material (say water). As far as I understand different E<span style="background-color:rgb(250,210,207)">quations of State </span>are implemented only for 3T mode and multigamma mode 

simply calculates mass averaged adiabatic index and substitutes it into equation of state for ideal gas. I could redefine EoS for ideal gas by tabulated EoS but of course I need to do it separately for different materials. Any ideas on how I can do it separately would be greatly appreciated!</div><div>Thank you!</div><div><span id="gmail-m_-2723892437778690474gmail-m_5332077176106748713gmail-m_4973181663998564798:1nk.16">MAksim</span> <span id="gmail-m_-2723892437778690474gmail-m_5332077176106748713gmail-m_4973181663998564798:1nk.17">Kozlov</span></div></div>
</blockquote></div>
</blockquote></div>
</div></blockquote></div><br></div></div></blockquote></div>