Menu
▾
▴
mesa-users
|
From: Ricky O. <ol...@sa...> - 2010-02-22 09:44:43
|
Hi, I am interested in using the MESA EOS and Opacity modules for use with my stellar pulsation code. As such, I have been poking around the MESA source files, and hacked the sample_eos.f source file to generate EOS tables in the format I need. However, it seems the MESA EOS routines automatically add the contribution due to radiation into the thermodynamic variables. I could easily subtract the contribution to radiation out, but is there a way of switching it off inside MESA? I treat radiation separately from the gas in the code I use. Can anyone help me with this? Regards, Ricky |
|
From: Aaron D. <aar...@gm...> - 2010-02-22 14:51:36
|
Hi Ricky, Actually, the eos module works with Pgas, see mesa/eos/public/eos_def.f for a list of variables that it works with. So, for example, the eos call that returns P(Rho,T) = eosDT_get actually returns Pgas. See also mesa/eos/public/eos_lib.f for a list of different eos subroutines. A function called Radiation_Pressure(T) is included in eos_lib as a convenience. Best wishes, Aaron On Mon, Feb 22, 2010 at 1:25 AM, Ricky Olivier <ol...@sa...> wrote: > Hi, > > I am interested in using the MESA EOS and Opacity modules for use with > my stellar pulsation code. As such, I have been poking around the MESA > source files, and hacked the sample_eos.f source file to generate EOS > tables in the format I need. However, it seems the MESA EOS routines > automatically add the contribution due to radiation into the > thermodynamic variables. I could easily subtract the contribution to > radiation out, but is there a way of switching it off inside MESA? I > treat radiation separately from the gas in the code I use. > > Can anyone help me with this? > > Regards, > Ricky > > > > ------------------------------------------------------------------------------ > Download Intel® Parallel Studio Eval > Try the new software tools for yourself. Speed compiling, find bugs > proactively, and fine-tune applications for parallel performance. > See why Intel Parallel Studio got high marks during beta. > http://p.sf.net/sfu/intel-sw-dev > _______________________________________________ > mesa-users mailing list > mes...@li... > https://lists.sourceforge.net/lists/listinfo/mesa-users > |
|
From: Ricky O. <ol...@sa...> - 2010-02-22 15:05:35
|
<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Hi Aaron,<br>
<br>
I have noticed that with regards to the pressure, but it appears not to
be the case for the internal energy, chiT, chiRho and Cv for example. I
need these quantities as evaluated for the gas only. In the code that
I use, I do not assume LTE and treat the radiation energy density and
radiation pressure separately from the gas.<br>
<br>
Thanks for your input, and the additional info. I am keen on using the
mesa routines, as it has all the top notch EOS in one place ...<br>
<br>
Ricky<br>
<br>
Aaron Dotter wrote:
<blockquote
cite="mid:95e...@ma..."
type="cite">Hi Ricky,
<div><br>
</div>
<div>Actually, the eos module works with Pgas, see
mesa/eos/public/eos_def.f for a list of variables that it works with.</div>
<div><br>
</div>
<div>So, for example, the eos call that returns P(Rho,T) = eosDT_get
actually returns Pgas. See also mesa/eos/public/eos_lib.f for a list
of different eos subroutines. A function called Radiation_Pressure(T)
is included in eos_lib as a convenience.</div>
<div><br>
</div>
<div>Best wishes,</div>
<div>Aaron</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
<br>
<div class="gmail_quote">On Mon, Feb 22, 2010 at 1:25 AM, Ricky
Olivier <span dir="ltr"><<a moz-do-not-send="true"
href="mailto:ol...@sa...">ol...@sa...</a>></span>
wrote:<br>
<blockquote class="gmail_quote"
style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi,<br>
<br>
I am interested in using the MESA EOS and Opacity modules for use with<br>
my stellar pulsation code. As such, I have been poking around the MESA<br>
source files, and hacked the sample_eos.f source file to generate EOS<br>
tables in the format I need. However, it seems the MESA EOS routines<br>
automatically add the contribution due to radiation into the<br>
thermodynamic variables. I could easily subtract the contribution to<br>
radiation out, but is there a way of switching it off inside MESA? I<br>
treat radiation separately from the gas in the code I use.<br>
<br>
Can anyone help me with this?<br>
<br>
Regards,<br>
Ricky<br>
<br>
<br>
------------------------------------------------------------------------------<br>
Download Intel&#174; Parallel Studio Eval<br>
Try the new software tools for yourself. Speed compiling, find bugs<br>
proactively, and fine-tune applications for parallel performance.<br>
See why Intel Parallel Studio got high marks during beta.<br>
<a moz-do-not-send="true" href="http://p.sf.net/sfu/intel-sw-dev"
target="_blank">http://p.sf.net/sfu/intel-sw-dev</a><br>
_______________________________________________<br>
mesa-users mailing list<br>
<a moz-do-not-send="true"
href="mailto:mes...@li...">mes...@li...</a><br>
<a moz-do-not-send="true"
href="https://lists.sourceforge.net/lists/listinfo/mesa-users"
target="_blank">https://lists.sourceforge.net/lists/listinfo/mesa-users</a><br>
</blockquote>
</div>
<br>
</div>
</blockquote>
</body>
</html>
|
|
From: Aaron D. <aar...@gm...> - 2010-02-22 15:13:49
|
Hi Ricky, Sorry I misunderstood; thanks for clarifying. At present, we do not have the capability to evaluate things like specific heats, ChiT, and ChiRho as functions of Pgas rather than total P. It is good to know that this is of interest. If you write some code to fulfill your needs, you are certainly encouraged to share it with us! Best wishes, Aaron On Mon, Feb 22, 2010 at 7:10 AM, Ricky Olivier <ol...@sa...> wrote: > Hi Aaron, > > I have noticed that with regards to the pressure, but it appears not to be > the case for the internal energy, chiT, chiRho and Cv for example. I need > these quantities as evaluated for the gas only. In the code that I use, I > do not assume LTE and treat the radiation energy density and radiation > pressure separately from the gas. > > Thanks for your input, and the additional info. I am keen on using the mesa > routines, as it has all the top notch EOS in one place ... > > Ricky > > > Aaron Dotter wrote: > > Hi Ricky, > > Actually, the eos module works with Pgas, see mesa/eos/public/eos_def.f > for a list of variables that it works with. > > So, for example, the eos call that returns P(Rho,T) = eosDT_get actually > returns Pgas. See also mesa/eos/public/eos_lib.f for a list of different > eos subroutines. A function called Radiation_Pressure(T) is included in > eos_lib as a convenience. > > Best wishes, > Aaron > > > > > On Mon, Feb 22, 2010 at 1:25 AM, Ricky Olivier <ol...@sa...> wrote: > >> Hi, >> >> I am interested in using the MESA EOS and Opacity modules for use with >> my stellar pulsation code. As such, I have been poking around the MESA >> source files, and hacked the sample_eos.f source file to generate EOS >> tables in the format I need. However, it seems the MESA EOS routines >> automatically add the contribution due to radiation into the >> thermodynamic variables. I could easily subtract the contribution to >> radiation out, but is there a way of switching it off inside MESA? I >> treat radiation separately from the gas in the code I use. >> >> Can anyone help me with this? >> >> Regards, >> Ricky >> >> >> >> ------------------------------------------------------------------------------ >> Download Intel® Parallel Studio Eval >> Try the new software tools for yourself. Speed compiling, find bugs >> proactively, and fine-tune applications for parallel performance. >> See why Intel Parallel Studio got high marks during beta. >> http://p.sf.net/sfu/intel-sw-dev >> _______________________________________________ >> mesa-users mailing list >> mes...@li... >> https://lists.sourceforge.net/lists/listinfo/mesa-users >> > > |
|
From: Ricky O. <ol...@sa...> - 2010-02-22 15:49:23
|
<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Hi,<br>
<br>
Ok, looks like this is going to be harder than I thought. I could
subtract off the effect of radiation to the various thermodynamic
variables. This should be ok in those parts of the
Rho-T plane where the radiation energy density and radiation pressure
do not dominate. In other parts of the Rho-T plane it will be dodgy to
say the least.<br>
<br>
One final question, so that I don't waste my time trying something that
won't work. The EOS tables the MESA routines interpolates on, is
radiation already included in these tables, or is it included
afterwards by the routines? <br>
<br>
Thanks,<br>
R.<br>
<br>
Aaron Dotter wrote:
<blockquote
cite="mid:95e...@ma..."
type="cite">Hi Ricky,
<div><br>
</div>
<div>Sorry I misunderstood; thanks for clarifying. At present, we do
not have the capability to evaluate things like specific heats, ChiT,
and ChiRho as functions of Pgas rather than total P. It is good to
know that this is of interest. If you write some code to fulfill your
needs, you are certainly encouraged to share it with us!</div>
<div><br>
</div>
<div>Best wishes,</div>
<div>Aaron</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
<div class="gmail_quote">On Mon, Feb 22, 2010 at 7:10 AM, Ricky
Olivier <span dir="ltr"><<a moz-do-not-send="true"
href="mailto:ol...@sa...">ol...@sa...</a>></span>
wrote:<br>
<blockquote class="gmail_quote"
style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div bgcolor="#ffffff" text="#000000">
Hi Aaron,<br>
<br>
I have noticed that with regards to the pressure, but it appears not to
be the case for the internal energy, chiT, chiRho and Cv for example. I
need these quantities as evaluated for the gas only. In the code that
I use, I do not assume LTE and treat the radiation energy density and
radiation pressure separately from the gas.<br>
<br>
Thanks for your input, and the additional info. I am keen on using the
mesa routines, as it has all the top notch EOS in one place ...<br>
<font color="#888888"><br>
Ricky</font>
<div>
<div class="h5"><br>
<br>
Aaron Dotter wrote:
<blockquote type="cite">Hi Ricky,
<div><br>
</div>
<div>Actually, the eos module works with Pgas, see
mesa/eos/public/eos_def.f for a list of variables that it works with.</div>
<div><br>
</div>
<div>So, for example, the eos call that returns P(Rho,T) =
eosDT_get
actually returns Pgas. See also mesa/eos/public/eos_lib.f for a list
of different eos subroutines. A function called Radiation_Pressure(T)
is included in eos_lib as a convenience.</div>
<div><br>
</div>
<div>Best wishes,</div>
<div>Aaron</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
<br>
<div class="gmail_quote">On Mon, Feb 22, 2010 at 1:25 AM, Ricky
Olivier <span dir="ltr"><<a moz-do-not-send="true"
href="mailto:ol...@sa..." target="_blank">ol...@sa...</a>></span>
wrote:<br>
<blockquote class="gmail_quote"
style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi,<br>
<br>
I am interested in using the MESA EOS and Opacity modules for use with<br>
my stellar pulsation code. As such, I have been poking around the MESA<br>
source files, and hacked the sample_eos.f source file to generate EOS<br>
tables in the format I need. However, it seems the MESA EOS routines<br>
automatically add the contribution due to radiation into the<br>
thermodynamic variables. I could easily subtract the contribution to<br>
radiation out, but is there a way of switching it off inside MESA? I<br>
treat radiation separately from the gas in the code I use.<br>
<br>
Can anyone help me with this?<br>
<br>
Regards,<br>
Ricky<br>
<br>
<br>
------------------------------------------------------------------------------<br>
Download Intel&#174; Parallel Studio Eval<br>
Try the new software tools for yourself. Speed compiling, find bugs<br>
proactively, and fine-tune applications for parallel performance.<br>
See why Intel Parallel Studio got high marks during beta.<br>
<a moz-do-not-send="true"
href="http://p.sf.net/sfu/intel-sw-dev" target="_blank">http://p.sf.net/sfu/intel-sw-dev</a><br>
_______________________________________________<br>
mesa-users mailing list<br>
<a moz-do-not-send="true"
href="mailto:mes...@li..." target="_blank">mes...@li...</a><br>
<a moz-do-not-send="true"
href="https://lists.sourceforge.net/lists/listinfo/mesa-users"
target="_blank">https://lists.sourceforge.net/lists/listinfo/mesa-users</a><br>
</blockquote>
</div>
<br>
</div>
</blockquote>
</div>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</blockquote>
</body>
</html>
|
|
From: Bill P. <pa...@ki...> - 2010-02-22 17:04:49
|
On Feb 22, 2010, at 7:53 AM, Ricky Olivier wrote:
> Ok, looks like this is going to be harder than I thought. I could
> subtract off the effect of radiation to the various thermodynamic
> variables. This should be ok in those parts of the Rho-T plane where
> the radiation energy density and radiation pressure do not dominate.
> In other parts of the Rho-T plane it will be dodgy to say the least.
Hi Ricky,
As you know, the MESA/EOS is a blend of SCVH/OPAL/HELM eos's.
HELM is for situations in which it is okay to assume complete
ionization,
and HELM does provide separate results with/without radiation.
So one possible approach is to by-pass the standard routine and
call the mesa/public/eos_lib routine 'eos_get_helm_results'.
In mesa/public/eos_def you'll find a long list of integer parameter
definitions for the results from HELM. There are separate items
for pressure ('p'), energy ('e'), and entropy ('s'). The partials are
wrt temperature ('d_t'), density ('d_d' ), zbar ('d_z'), and abar
('d_a').
And there are 2nd partials, 'd_dt', etc. The 'total' value is the
sum of the contributions from gas and radiation -- that what mesa/eos
reports of course. But you will also find the separate terms for
gas and radiation -- e.g., ptot, pgas, prad, etc.
For derived quantities such as specific heats and gammas, you
will find both values based on totals and on for the gas alone:
e.g., h_cp for the total and h_cp_gas for the gas alone.
So for the regions where you are willing to assume complete ionization,
you can go straight to the HELM results and get what you want.
Perhaps that will help?
Let me know.
Cheers,
Bill
|
×
Want the latest updates on software, tech news, and AI?
Get latest updates about software, tech news, and AI from SourceForge directly in your inbox once a month.
Thanks for helping keep SourceForge clean.
X