< Previous | Contents | Next >
Как уже отмечалось, главным методом исследования нуклеосинтеза тяже- лых ядер и в целом процессов термоядерного горения в звeздах является ком- пьютерное моделирование. Математически задача ядерных трансмутаций пред- ставляет собой систему обыкновенных дифференциальных уравнений первого порядка. Специфика рассматриваемых процессов обуславливает ряд особенно- стей этой системы, усложняющих ее решение: огромные размеры матриц, раз- личие коэффициентов уравнений на десятки порядков, неопределенность на- чальных условий. В связи с этим многие существующие численные методы ре- шения систем обыкновенных дифференциальных уравнений при расчете про- цессов нуклеосинтеза оказываются либо недостаточно точны, либо слишком трудоемки для использования в программах глобальной симуляции звездной эволюции.
В данном разделе проводится исследование различных методов численного решения систем обыкновенных дифференциальных уравнений на предмет их применимости к задаче термоядерного астрофизического горения. В частно- сти, предпринимается попытка программной реализации методов QSS [17, 18] и специальной явной схемы [19, 20] на основе системы астрофизической симу- ляции MESA [21]. Выводы, сделанные в настоящей работе, дожны помочь в поиске наиболее оптимальной расчетной схемы, на основе которой в дальней- шем можно будет с большей эффективностью моделировать r-процесс и другие механизмы звездного нуклеосинтеза.