计算物理 ›› 1987, Vol. 4 ›› Issue (3): 329-338.

• 论文 • 上一篇    下一篇

恒星的星核区坍缩与物态方程(二)——恒星星核区坍缩的计算

王贻仁1, 黄维章1, 张锁春1, 张天树2   

  1. 1. 中国科学院应用数学研究所;
    2. 应用物理与计算数学研究所 北京
  • 收稿日期:1986-12-13 出版日期:1987-09-25 发布日期:1987-09-25

STELLAR CORE COLLAPSE AND EQUATION OF STATE (Ⅱ)——COMPUTER SIMULATION OF STELLAR COLLAPSE

Wang Yi-ren1, Huang Wei-zhang1, Zhang Suo-chun1, Zhang Tian-shu2   

  1. 1. Institute of Applied Math. Acadimia-Sinica;
    2. Institute of Applied Phys, and Computational Math
  • Received:1986-12-13 Online:1987-09-25 Published:1987-09-25

摘要: 本文用文[17]所提供的物态方程等输入物理因素,以Weaver等的质量为15M的恒星演化模型为初始模型,对恒星坍缩过程作了数值计算,计算至中心密度为1014gcm-3。在由平均重核,α粒子及自由核子组成的简化核平衡系统基础上,对核的裂解与复合、重核能级激发及电子俘获使物质中子化等过程热效应的处理作了尝试。计算结果表明,我们的结果(用EOS(Ⅰ)态方程)与Van Riper等(用LLPR物态方程)的典型模型(非相对论)结果很接近。当中心密度ρ(A)达1.0×1014gcm-3。最大坍缩速度为3.8×109cm sec-1,相似区质量Ms为0.61M较Van Riper典型模型值分别大2.7%和7%。在质量壳0.15M及0.5M,相应的中微子俘陷密度ρtrap为1.55×1011及1.8×1011g cm-3,较VanRiper值分别小3%及10%,而相应的轻子浓度为0.387及0.386较VanRiper分别大895及10%。

Abstract: We take 15 M model of Weaver et al. as initial model,the collapse process of stellar core between the onset of collapse and core bounce has been calculated by using the eouation of state of us (EOSⅡ).With the simplified nuclear statistical equilibrium system mentioned above, the treatment of the thermal effects of neutronization, nuclear disintegration, nuclear combination and nuclear excitation has been attemped.The results of us are essential agreed with Van Riper's results.calculated by using the ILPR's equation of state for the same initial model.As the center density of core reaches 1014gcm-3, the maximum velocity of fall and the Ms (mass at sonic point)in model D is greater than that of the typical model of Van Riper by 2.7% and 7% respectively. In model D the neutrino trapping density is less than that of the typical model of Van Riper by 3% (at mass layer 0.15M) but the lepton fraction in model D is greater than that of typical model of Van Riper by 8%(at 0.15M) and by 10% (at 0.5M).