Порядок моделирования

Моделирование каскада включает в себя два этапа. На первом этапе проводится предварительная подгонка. Следя за значением давления на фронте ударной волны, определяют примерное время окончания каскадного процесса t _ок, т.е. момент, когда давление упадет ниже величины 1/2p = 0,1592, и подбирают величину числа Куранта. При малых энерговыделениях наблюдать ударную волну на графике Р (r) не так просто, так как на малых временах она проявляет себя лишь небольшим перегибом на спадающей части графика. Здесь может оказаться более удобным следить за положением фронта по графику u (r).

Далее, выйдя из программы командой [E]/[Ent] и запустив ее снова (при этом необходимо повторно ввести начальные пара­метры), провести расчет для 10–12 последовательных времен, более-менее равномерно охватывающих интервал от 0 (обязательно включая саму точку t = 0) до t _ок. Последнее время t долж­но несколько превышать t _ок, чтобы иметь возможность интерполя­ции.

Для каждого момента времени необходимо выписать в лаборатор­ный журнал из таблиц зависимость Р (r) и величину перемещенной массы, а также определить положение фронта волны. для последнего (наиболее близкого к t _ок) момента также численные значения для зависимости r (r).

ВНИМАНИЕ! Вследствие необратимости процесса релаксации расчеты допускаются только для возрастающей последовательности времен. Если вы неправильно ввели очередное значение времени или не успели записать результаты, от­носящиеся к более ранним моментам времени, вам придется выйти из программы и повторить расчет для нужных моментов времени снова.

ОБРАБОТКА РЕЗУЛЬТАТОВ

Обработка результатов моделирования включает в себя постро­ение графиков, расчет характеристик каскадного процесса и срав­нение полученных результатов с предсказаниями теоретических моделей.

Перевод из безразмерных величин в размерные

По величине плотности начального энерговыделения рассчитать энергию Е первично выбитого атома. Для расчета положить, что на кинетической стадии энергия передается среде в объеме V = 4p L ³/3 и воспользоваться приближенной оценкой (см. [l]):

L = Е /(2p² СN)

где N = r ₀/ m _a – атомная плотность (m _a – средняя масса ато­ма), а С = 2.5·10^–17·Z^5/3 (эВ·см²).

Определив энергию ПВА, рассчитать значения масштабных мно­жителей преобразования безразмерных величин в размерные. При пе­ресчете воспользоваться характеристиками материала, приведенными в таблице I.