[Mesa-users] Modeling a sub-Saturn planet using MESA
Evan Bauer
evan.bauer.astro at gmail.com
Mon Mar 18 21:24:50 UTC 2024
Hi Jiapeng,
Could you also send the inlists that worked for you in r11532? It would be helpful to make sure we can start with a baseline of what worked in a previous MESA version. Thanks!
Cheers,
Evan
> On Mar 18, 2024, at 2:59 PM, Frank Timmes via Mesa-users <mesa-users at lists.mesastar.org> wrote:
>
> hi jiapeng,
>
> i'll confirm that that the make_planets test suite case in r24.03.1 has difficulties below about 0.28 m_jup. its being looked at.
>
> fxt
>
>
>
>> On Mar 12, 2024, at 11:47 AM, Gao, Jiapeng via Mesa-users <mesa-users at lists.mesastar.org> wrote:
>>
>> Hello,
>>
>> My name is Jiapeng Gao, and I'm a first year graduate student. I'm
>> trying to simulate planet models using MESA, and I'm having the
>> following problems:
>>
>>
>> 1. I am trying to use the "make_planets" test case in the MESA r23.05.1
>> test suite to model a planet of 0.2 Jupiter mass. However, I'm having
>> trouble creating the initial model. I get the following error message:
>> "failed in kap_get get_kap_from_rhoT, temp too low in integration". It
>> seems to work with MESA r11532, where I can create a 0.2 M_J planet with
>> 1 Jupiter radius. I'm not sure how I could modify r23.05.1 to fix the
>> problem.
>>
>>
>> 2. I noticed that for my simulation using MESA r11532, the planet radius
>> can only shrink to around 0.8 R_J at ~10Gyr, then the program stops and
>> reports: “dt 1.9D-23, min_timestep_limit 9.9 D-21, stopping because of
>> problems dt < min_timestep_limit, terminated evolution: cannot find
>> acceptable model, termination code: min_timestep_limit”. I'm not sure
>> what would be the best way to troubleshoot it.
>>
>> I included my inlists below. Thank you very much for your guidance!
>>
>> Sincerely,
>>
>>
>> Jiapeng Gao
>> ———————————MESA r23.05.1 inlist_create———————————————————
>>
>> ! inlist_create
>>
>> &star_job
>>
>>
>> show_log_description_at_start = .false.
>>
>> save_model_when_terminate = .true.
>> save_model_filename = "planet_create_0.9685_MJ_2.0_RJ.mod"
>> required_termination_code_string = 'max_age'
>>
>> ! new adiabatic, contracting initial model
>> create_initial_model = .true.
>> radius_in_cm_for_create_initial_model = 6991100000.d0 ! in cm
>> mass_in_gm_for_create_initial_model = 3.6778d29 ! in grams
>> initial_model_relax_num_steps = 0
>> set_initial_age = .true.
>> initial_age = 0.d0 ! in years
>>
>> !pgstar_flag = .true.
>>
>> / !end of star_job
>>
>>
>> &eos
>> / ! end of eos namelist
>>
>>
>> &kap
>> Zbase = 0.02d0
>>
>> kap_file_prefix = 'gs98'
>> kap_lowT_prefix = 'lowT_Freedman11'
>>
>> / ! end of kap namelist
>>
>> &controls
>>
>> use_gold2_tolerances = .false. ! wait a few steps
>> steps_before_use_gold2_tolerances = 31
>>
>> warning_limit_for_max_residual = 1d99 ! off until evolve
>> warn_when_large_virial_thm_rel_err = 1d99 ! off until evolve
>>
>> limit_for_rel_error_in_energy_conservation = 1d-5
>> hard_limit_for_rel_error_in_energy_conservation = 1d-3
>>
>> convergence_ignore_equL_residuals = .true.
>>
>> ! limit max_model_number as part of test_suite
>> max_model_number = 200
>>
>> ! for create_initial_model
>> initial_Z = 0.02d0
>> initial_Y = 0.24d0
>>
>>
>> ! when to stop
>> max_age = 50
>>
>> ! knobs
>> T_mix_limit = 0d0
>>
>> max_resid_jump_limit = 1d18
>> max_corr_jump_limit = 1d18
>>
>> ! output to files and terminal
>> max_num_profile_models=1000
>>
>> photo_interval = 50
>> profile_interval = 50
>> history_interval = 1
>> terminal_interval = 10
>> write_header_frequency = 10
>>
>> /
>>
>>
>> &pgstar
>>
>> ! top level controls
>>
>> !pause = .true.
>> ! if true, the code waits for user to enter a RETURN on the command line
>>
>>
>>
>> / ! end of pgstar namelist
>>
>> --------------------------MESA r11532 inlist_evolve -------------------------------
>>
>> ! inlist_evolve
>>
>> &star_job
>>
>> mesa_dir = '/home/jiapeng/Downloads/mesa-r11532'
>>
>> show_log_description_at_start = .false.
>>
>> load_saved_model = .true.
>> saved_model_name = "planet_core_1.0_MJ_10.0_ME_2.0_RJ.mod"
>>
>> save_model_when_terminate = .true.
>> save_model_filename = "planet_evolve_1.0_MJ_10.0_ME_2.0_RJ.mod"
>> required_termination_code_string = 'max_age'
>>
>> set_initial_cumulative_energy_error = .true.
>> new_cumulative_energy_error = 0d0
>>
>> !relax_initial_L_center = .true.
>> !new_L_center = 4.46863d28 ! ergs/sec
>> !dlgL_per_step = 5d-3
>> !relax_L_center_dt = 8.1558149984d1 ! secyer*1d-6
>>
>> eos_file_prefix = 'mesa'
>> kappa_file_prefix = 'gs98'
>> kappa_lowT_prefix = 'lowT_Freedman11'
>>
>> !set_tau_factor = .true.
>> !set_to_this_tau_factor=100
>>
>> change_lnPgas_flag = .true.
>> new_lnPgas_flag = .true.
>>
>> !pgstar_flag = .true.
>>
>> / !end of star_job
>>
>> &controls
>> use_gold_tolerances = .true.
>> use_eosDT2 = .true.
>> use_eosELM = .true.
>> logQ_limit = 1d99
>>
>> min_timestep_limit = 1d-20
>>
>> warn_when_large_virial_thm_rel_err = 5d-2
>>
>> use_dedt_form_of_energy_eqn = .true.
>> min_cell_energy_fraction_for_dedt_form = 1d-8
>>
>> num_trace_history_values = 2
>> trace_history_value_name(1) = 'rel_E_err'
>> trace_history_value_name(2) = 'log_rel_run_E_err'
>>
>> ! check for retries and backups as part of test_suite
>> max_number_backups = 200
>> max_number_retries = 200
>> max_model_number = 3500
>>
>>
>> ! surface heating
>> ! column_depth_for_irradiation = 300.d0 ! 3.d2 cm^2/g
>> ! irradiation_flux = 6.2d8 ! 1.d9 erg/cm^2/s ! day side flux!!! Bill puts in 1/4 by hand
>>
>> ! interior heating
>> !inject_uniform_extra_heat = 0.00052
>> !min_q_for_uniform_extra_heat = 0.0705882
>> !max_q_for_uniform_extra_heat = 0.1
>> !inject_extra_ergs_sec=4.46863d29
>> !base_of_inject_extra_ergs_sec=0.0
>> !total_mass_for_inject_extra_ergs_sec=0.00009
>>
>> ! when to stop
>> max_age = 100000000000.d0 ! 1.d10
>> Teff_lower_limit = 10.d0
>>
>> ! knobs
>> !use_lnE_for_eps_grav = .false.
>> T_mix_limit = 0
>>
>>
>> ! output to files and terminal
>> photo_interval = 50
>> profile_interval = 50
>> max_num_profile_models=1000
>> history_interval = 2
>> terminal_interval = 10
>> write_header_frequency = 10
>>
>>
>> !photo_interval = 10
>> !profile_interval = 10
>> !history_interval = 1
>> !terminal_interval = 1
>>
>>
>> ! FOR DEBUGGING
>>
>> !report_hydro_solver_progress = .true. ! set true to see info about newton iterations
>> !report_ierr = .true. ! if true, produce terminal output when have some internal error
>> !hydro_show_correction_info = .true.
>>
>> ! hydro debugging
>> !hydro_check_everything = .true.
>> !hydro_inspectB_flag = .true.
>> !hydro_sizequ_flag = .true.
>>
>> !hydro_get_a_numerical_partial = 1d-4
>> !hydro_test_partials_k = 1
>> !hydro_numerical_jacobian = .true.
>> !hydro_save_numjac_plot_data = .true.
>> !hydro_dump_call_number = 195
>> !hydro_dump_iter_number = 5
>> !hydro_epsder_struct = 1d-6
>> !hydro_epsder_chem = 1d-6
>>
>> !fill_arrays_with_NaNs = .true.
>> !stop_for_NaNs = .true.
>>
>> !max_years_for_timestep = 3.67628942044319d-05
>>
>> !report_why_dt_limits = .true.
>> !report_all_dt_limits = .true.
>> !report_hydro_dt_info = .true.
>>
>> !show_mesh_changes = .true.
>> !mesh_dump_call_number = 5189
>> !okay_to_remesh = .false.
>>
>> !trace_evolve = .true.
>>
>> !trace_newton_bcyclic_solve_input = .true. ! input is "B" j k iter B(j,k)
>> !trace_newton_bcyclic_solve_output = .true. ! output is "X" j k iter X(j,k)
>>
>> !trace_newton_bcyclic_matrix_input = .true.
>> !trace_newton_bcyclic_matrix_output = .true.
>>
>> !trace_newton_bcyclic_steplo = 1 ! 1st model number to trace
>> !trace_newton_bcyclic_stephi = 1 ! last model number to trace
>>
>> !trace_newton_bcyclic_iterlo = 2 ! 1st newton iter to trace
>> !trace_newton_bcyclic_iterhi = 2 ! last newton iter to trace
>>
>> !trace_newton_bcyclic_nzlo = 1 ! 1st cell to trace
>> !trace_newton_bcyclic_nzhi = 10000 ! last cell to trace; if < 0, then use nz as nzhi
>>
>> !trace_newton_bcyclic_jlo = 1 ! 1st var to trace
>> !trace_newton_bcyclic_jhi = 100 ! last var to trace; if < 0, then use nvar as jhi
>>
>> !trace_k = 0
>>
>> /
>>
>>
>> &pgstar
>>
>> ! top level controls
>>
>> !pause = .true.
>> ! if true, the code waits for user to enter a RETURN on the command line
>>
>>
>>
>> / ! end of pgstar namelist
>>
>>
>> _______________________________________________
>> mesa-users at lists.mesastar.org <mailto:mesa-users at lists.mesastar.org>
>> https://lists.mesastar.org/mailman/listinfo/mesa-users <https://lists.mesastar.org/mailman/listinfo/mesa-users>
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240318/d389f4e0/attachment.htm>
More information about the Mesa-users
mailing list