Notebook sobre o teste_suite make_o_ne_wd do MESA. Todos os testes e resultados estarão neste notebook

Make_O_Ne_WD

Teste original



In [ ]:

Ligando a Difusão de elementos

Adicionado o parâmetro:

Utilizando $\alpha_{MLT} = 2$

o inlist original do teste define $\alpha_{MLT} =1$ atravez do parametro:

gradT_excess_lambda1 = -1

assim, será feito o teste com as seguintes modificações:

comentando o parâmetro gradT_excess_lambda1
Ligando a difusão (do_element_diffusion = .true.)
Desligando a rotação (comentando os parametros de rotação)
definindo alpha = 2 ( mixing_length_alpha = 2.0 )

Modificando a quantidade de camadas

Triplicando as camadas

modificando o parâmetro de camadas:

mesh_delta_coeff = 1.0d0 (padrão do teste = 1.5d0)

e adicionando os parâmetros:

mesh_delta_coeff_for_highT = 0.33d0
logT_max_for_standard_mesh_delta_coeff = 8.0d0
logT_min_for_highT_mesh_delta_coeff = 8.6d0

pretendo ~triplicar a quantidade de camadas na faixa de temperatura entre 8 e 8.6 $\log T_c$ ($T_c =$ temperatura central)

adicionei os parametros num_trace_history_values = 1 e trace_history_value_name(1) = 'num_zones' para printar no terminal a quantidade de camadas em cada modelo

Resultados



In [4]:

    
#Diagrama HR
'''
o pontos representam a regiao entre o modelos 155 e 2150, onde deveria ocorrer os pulsos termicos
'''
from IPython.display import Image
Image(filename='figures/alpha2_hrd_3x1.png')









    Out[4]:



In [5]:

    
# Camadas por modelo
'''
O grafico mostra o aumento da quantidade de camadas, sendo que na regiao que deveria ocorrer os pulsos
ha um aumento de 1500 para 4000 camadas
'''
from IPython.display import Image
Image(filename='figures/alpha2_3x_camadas.png')









    Out[5]:



In [6]:

    
# luminosidade de helio por modelo
'''
nao ha pulsos termicos...
'''
from IPython.display import Image
Image(filename='figures/alpha2_3x_log_L_he.png')









    Out[6]:



In [7]:

    
#o pequeno pico ocorre no modelo 440. A sua posicao esta marcada no HRD abaixo
from IPython.display import Image
Image(filename='figures/alpha2_3x_hrd_zoom.png')









    Out[7]:



In [9]:

    
# abundancia de elementos
from IPython.display import Image
Image(filename='figures/alpha2_3x_abun.png')









    Out[9]:

Multiplicando por 10 as camadas

mesmas modificações do item anterior (Triplicando as camadas) porem modifiquei o parâmetro:

mesh_delta_coeff_for_highT = 0.1d0

Resultado

Erro ao no run ao chegar no modelo 590 (estava com 7495 zonas)

*** Error in `./star': corrupted double-linked list: 0x000000000cb6f130 ***

Program received signal SIGABRT: Process abort signal.

Backtrace for this error:
#0  0x7f6a9bdfbc02
#1  0x7f6a9bdfaf30
#2  0x7f6a9aeb4caf
#3  0x7f6a9aeb4c37
#4  0x7f6a9aeb8027
#5  0x7f6a9aef12a3
#6  0x7f6a9aefd55d
#7  0x4927bb
#8  0x492fa7
#9  0x68db79
#10  0x56b85a
#11  0x41372c
#12  0x42947b
#13  0x41068c
#14  0x4106dd
#15  0x410714
#16  0x7f6a9ae9ff44
#17  0x40fbcc
#18  0xffffffffffffffff
Aborted (core dumped)

O diagrama HR mostrando até onde foi calculado a evolução esta a seguir:



In [1]:

    
from IPython.display import Image
Image(filename='figures/alpha2_10x.png')









    Out[1]:

Triplicando as camadas e modificando os valores mínimos de elementos

modificando o parâmetro de camadas:

mesh_delta_coeff = 1.0d0 (padrão do teste = 1.5d0)

e adicionando os parâmetros:

mesh_delta_coeff_for_highT = 0.33d0
logT_max_for_standard_mesh_delta_coeff = 8.0d0
logT_min_for_highT_mesh_delta_coeff = 8.6d0
dX_mix_dist_limit = 1d-30



In [ ]: