[FLASH-USERS] value of Ye for single species runs
Ernesto Zurbriggen
ernesto at mail.oac.uncor.edu
Fri Oct 2 16:45:40 EDT 2015
Jon,
looking at the end of the subroutine Eos/EosMain/Eos_getAbarZbar.F90
you can see for a single species run Ye=Zbar/Abar too. And I think
that is not a stub implementation for the single-species gamma-law EOS.
On the other hand, if you turn on the thermal conduction managed by
Diffuse/Split or Diffuse/Unplit units, I think the runtime parameter
"addThermalFlux" must be false. On the contrary, you might be taking
into account thermal conduction twice, within the Diffuse unit and
Hydro unit via hy_uhd_addThermalFluxes.F90
Best!
Quoting John ZuHone <jzuhone at space.mit.edu>:
> Just for reference, Jon and I have been discussing this off-list.
> The single-species gamma-law EOS only has a stub for this
> subroutine. I suggested copying the Eos_getAbarZbar stub to the
> simulation directory and changing it to set the quantities to
> whatever he needs, but as I am not familiar with this particular
> diffusion module I think someone with more knowledge should comment
> if this is the right way to go.
>
>> On Oct 2, 2015, at 3:32 PM, Slavin, Jonathan
>> <jslavin at cfa.harvard.edu> wrote:
>>
>> Hi,
>>
>> I have been doing runs trying to simulate the evaporation of
>> interstellar clouds via thermal conduction. The runs have been
>> working but the clouds are evaporating a bit faster than I
>> expected. This led me to look into the implementation of flux
>> limitation. Looking at diff_advanceTherm there is a call to
>> Eos_getAbarZbar, which returns the value of Ye, but it's unclear
>> what value will be returned for a single species run. I have been
>> setting eos_singleSpeciesA and that seems to work as expected. I
>> know in multispecies runs Ye = zbar/abar, so I wonder if its
>> 1/abar. Can anyone help me out on this?
>>
>> Regards,
>> Jon
>>
>> --
>> ________________________________________________________
>> Jonathan D. Slavin Harvard-Smithsonian CfA
>> jslavin at cfa.harvard.edu <mailto:jslavin at cfa.harvard.edu> 60
>> Garden Street, MS 83
>> phone: (617) 496-7981 Cambridge, MA 02138-1516
>> cell: (781) 363-0035 USA
>> ________________________________________________________
>>
>
>
----------------------------------------------------------------
This message was sent using IMP, the Internet Messaging Program.
More information about the flash-users
mailing list