[Mesa-users] a problem of the relaxation

[曾耀田]

HiPablo,




Thanksforyourreply!

I'msorryIdidn'tclarifytheissuebefore. Iwanttoimporta3DSPHmodelintoMESA. Nowitseemssucceedsinrelaxingbutthenfailstoevolvefurther.
TheattachmentsaretheinlistfileIused, the1Daveraged-density, internalenergy, composition profiles, andtheoutputinformationfile (out) ontheruntimescreen.
WhatIdon'tquiteunderstand:
(1)  Fromtheoutfile, relaxationseemstobecompletedundertheconditions (forrelax_composition, stoppingwhenstep == 300; forrelax_entropy, stoppingwhenage_yr3.0075E-07 >= num_timescales_for_relax_entropy * timescale_for_relax_entropy = 300*1d-9). Afterrelaxiscompleted, theprogramseemstocontinuetoevolveusingthecompletedrelaxmodel. Infact, Iamnotsurewhethertheprocedureiscarriedoutthisway.
 (2)  Iftheprogramisrunasdescribedin (1), itseemsthatthemodelafterrelaxhasnotsuccessfullyevolved. Doesitmeanthatrelaxationhasfailedatthistime? Oritisbecauseofotherreasonsthatithasfailedtocontinuetoevolve. Further, myquestionishowtojudgewhetherrelaxationissuccessful?
(3) Itseemsthatnoprofilefileisoutputduringtherelaxingprocess. Isitpossibletooutputtheprofilefileduringtherelaxingprocesstocheckhowthestructureofthestarchangesintheprocesstojudgewhethertherelaxissuccessful?

Anyideas, suggestions, orhelpwouldbegreatlyappreciated!




Best wishes
Yaotian Zeng

 
 

-----原始邮件-----
发件人:"PabloMarchant" <pamarca at gmail.com>
发送时间:2020-06-1716:03:30 (星期三)
收件人: "曾耀田" <zengyaotian at ynao.ac.cn>
抄送: "FrancisTimmes" <fxt44 at mac.com>, mesa-users <mesa-users at lists.mesastar.org>
主题:Re: [Mesa-users] aproblemoftherelaxation


HiYaotian,


inordertohelpyouwithyourproblemweneedenoughinformationtoreproduceit, includingallinlistsandinputfilesattachedtotheemail. FromyouremailI'malsonotsureifyoursimulationsucceedsinrelaxing, butthenfailstoevolvefurther.


Cheers


OnMon, Jun8, 2020at4:34PM 曾耀田 <zengyaotian at ynao.ac.cn> wrote:


Hi Francis,




Thanks for your reply. 


Iwanttomap3DmodelstoMESA. ThereforethefirststepIalreadyhavedoneistoconvert3Dmodelsinto1Dmodels. However, these1Dmodelsareoutofhydrostaticequilibrium. Toobtain1Dmodelsthatareinthehydrostaticequilibrium, Iamtryingtorelaxthemforseveraldynamicaltimescalesuntilthevelocityandaccelerationofeachcellofastararesmallenough. ConsideringIamnotfamiliarwithMESA, whatisdoingnowistorelaxthemodelthatMESAitselfcanruntolearnhowtorelaxastar.




Best wishes
Yaotian Zeng


> -----原始邮件-----
> 发件人: "Francis Timmes" <fxt44 at mac.com>
> 发送时间: 2020-06-08 07:57:33 (星期一)
> 收件人: "曾耀田" <zengyaotian at ynao.ac.cn>
> 抄送: "Frank Timmes" <fxt44 at mac.com>, mesa-users <mesa-users at lists.mesastar.org>
> 主题: Re: [Mesa-users] a problem of the relaxation
> 
> hi yaotian,
> 
> it would be useful to know what you are trying to achieve,
> in addition to how you tried to achieve it. what is the goal?
> 
> fxt
> 
> 
> 
> 
> > On Jun 7, 2020, at 3:03 AM, 曾耀田 <zengyaotian at ynao.ac.cn> wrote:
> > 
> > 2020-06-07 17:52:11"曾耀田" <zengyaotian at ynao.ac.cn>写道：
> > 
> > 
> > Dear all,
> > 
> > 
> > 
> > I'm sorry to bother you, but for some reason the information in this letter is missing. I am new to MESA. I am trying to relax a model of a 1.5 solar mass star based on the result (final.mod) of 1.5M_with_diffusion in test_suite directory. To do this, as same as the example ofrelax_composition_j_entropy in the test_suite directory, I try to relax composition and entropy (but using get_entropy_for_relax_from_eos = ‘eosDT', see details below). It can run some steps but always ended with the error of dt < min_timestep_limit. I have tried to set different values for the comments in the inlist (shown below) but failed (e.g. changing the values of min_timestep_limit to 1d-30). I don't understand that even the input model has exactly the same entropy, density, and temperature profiles as the 1.5 solar mass star already had, the relaxation still doesn't work. So I try to understand the meaning of adjusting the values of these commands through the .defults file but I still don't have a clue. Now I am unsure of how and what parameters to adjust. Any ideas, suggestions, or help would be greatly appreciated!




祝好！

曾耀田

_______________________________________________
mesa-users at lists.mesastar.org
https://lists.mesastar.org/mailman/listinfo/mesa-users




--

PabloMarchantCampos
M.SconAstrophysics, UniversidadCatólicadeChile
PhDonAstrophysics, Argelander-InstitutfürAstronomie, UniversitätBonn



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20200623/f8773d37/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_relax_composition_j_entropy
Type: application/octet-stream
Size: 1911 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20200623/f8773d37/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: out
Type: application/octet-stream
Size: 421934 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20200623/f8773d37/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: composition.dat
Type: application/octet-stream
Size: 19094 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20200623/f8773d37/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: entropy.dat
Type: application/octet-stream
Size: 6366 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20200623/f8773d37/attachment-0003.obj>