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> |