Содержание к диссертации
Введение
Глава 1. Априорные методы исследования свойств разностных схем 18
1. Механизмы диссипации энергии в методах расчета ударных волн 18
2. Разностные схемы в дифференциальном представлении 23
3. Акустическое приближение... 27
4. Метод исследования диссипативных свойств разностной схемы 28
5. Метод исследования дистракции сильного разрыва 30
6. Метод исследования немонотонности 35
Глава 2. Анализ свойств разностных схем 39
1. Разностная схема Д.Неймана-Р.Рихтмайера 39
2. Разностная схема П. Лакса 51
3. Разностная схема С.К.Годунова 58
4. Недивергентная разностная схема В.Ф.Куропатенко 64
5. Разностная схема П.Лакса, Б.Вендрофа 71
Глава 3. Новая разностная схема 75
1. Выбор сетки. Типы интервалов 75
2. Разностные уравнения для ударной волны 77
3. Погрешности аппроксимации на ударной волне 79
4. Анализ устойчивости разностной схемы для ударной волны 83
5. Анализ монотонности и дистракции разностной схемы на ударной волне 85
6. Разностные уравнения для волны разрежения 87
7. Погрешности аппроксимации на волне разрежения 89
8. Анализ устойчивости разностной схемы на волне разрежения 90
9. Анализ монотонности разностной схемы на волне разрежения 91
10. Повышение порядка аппроксимации 92
11. Уменьшение немонотонности на слабых разрывах 96
12. Краткое описание программы КАМА 99
13. Верификация разностной схемы 104
Глава 4. Исследование влияния свойств разностных схем на моделирование разрушения веществ 111
1. Характерные погрешности за фронтом ударной волны. Дистракция и осцилляции 11
2. Выход стационарной ударной волны на свободную поверхность. Аналитическое решение и результаты расчетов 112
3. Взаимодействие двух волн разрежения с образованием откола. Аналитическое решение 114
3.1. Область стационарного течения за фронтом ударной волны 116
3.2. Область центрированной волны разрежения 117
3.3. Область взаимодействия двух волн разрежения 120
3.4. Точка смены краевого условия 124
3.5. Течение в области у свободной границы 124
3.6. Масса отколовшегося слоя 127
4. Зависимость положения трещины от дистракции и осцилляции разностной схемы 130
Выводы 134
Список литературы 136
- Метод исследования диссипативных свойств разностной схемы
- Недивергентная разностная схема В.Ф.Куропатенко
- Анализ устойчивости разностной схемы для ударной волны
- Взаимодействие двух волн разрежения с образованием откола. Аналитическое решение
Введение к работе
В основе моделей механики сплошной среды, описывающих поведение вещества под действием динамических нагрузок, лежат законы сохранения массы, количества движения и энергии. Для разных задач эти законы сохранения записываются в разных формах. Ниже мы будем рассматривать их, следуя [1]. В случае непрерывных течений законы сохранения образуют систему дифференциальных уравнений в частных производных. В прямоугольных декартовых координатах эта система имеет вид dV dt VdivU = 0. (0.1) dU dt VgradP f ЛП asx asy asz dx dy 9z —+—y-+ z = 0, (0.2) div((pS + P)U-q) dps dT У +— z dz _ч fasvu 5svu asul dy dx = 0, (0.3) где U - скорость, p - плотность, V = 1/p - удельный объем, є = Е + 0.5U , Е -удельная внутренняя энергия, Р - давление, q - тепловой поток, t - время. Векторы (0.4) Sx Sxxi +Хху J +ТХ2 gy=TyxT + SyJ + Ty Sx=xzxi+xzJ + Szzk. образуют девиатор тензора напряжений, а давление Р является шаровой частью тензора напряжений. Векторы і, j, k - единичные векторы ортогонального базиса. Система уравнений (0.1)-(0.3) замыкается уравнениями состояния F(P,V,E)=0, Ф(Т,У,Е)=0, (0.5) определяющими уравнениями для q, S , ц и уравнением траектории частиц § = 0. (0.6)
При интенсивных динамических воздействиях на материалы, таких, как электровзрыв фольги, детонация взрывчатого вещества, удар пластиной, поток электронов или частиц, мощное рентгеновское, лазерное или другие виды излучений, как правило, возникают ударные волны. Распространяясь по веществу и взаимодействуя со свободными или контактными поверхностями, они могут отражаться волнами разрежения. При взаимодействии встречных волн разрежения возникают растягивающие напряжения. Под действием растягивающих нагрузок сначала рвутся наиболее слабые связи между кристаллами и зернами, затем эти микроразрушения сливаются, и возникает трещина.
Разрушение вещества в различных условиях зависит от характера динамического воздействия. В одних экспериментах определяющей причиной разрушения является большая дисторсия и связанная с ней работа девиатора тензора напряжений. При этом дилатация, а значит, работа давления PdV может быть малой. В других экспериментах причиной разрушения являются большие отрицательные давления, а работа напряжений сдвига может быть малой. В большом количестве экспериментальных работ по изучению откольного разрушения опыты ставятся так, что нормаль к фронту ударной волны ортогональна свободной поверхности вещества, и до наступления разрушения движение можно рассматривать как плоское. В таких экспериментах разрушение происходит, в основном, из-за больших растягивающих давлений.
Основным методом исследования откольного разрушения является физический эксперимент, основным методом прогнозирования разрушения является математическое моделирование. Существует много различных моделей, описывающих отклик вещества на динамическое воздействие, в том числе моделей прочности и разрушения. Следует заметить, что в каждой модели можно выделить адиабатическое ядро, в котором присутствует только шаровая часть тензора напряжений, т.е. давление. Под адиабатическим ядром понимается система законов сохранения в дифференциальной форме без учета теплопроводности, девиаторов тензора напряжений и энерговыделения. Общая погрешность математической модели может быть представлена в виде суммы погрешностей физической модели, погрешности аппроксимации адиабатического ядра и погрешности аппроксимации уравнений, порождаемых девиатором тензора напряжений. А = АФ +АМ +АМ. - мод ая "д •
Следствием уравнений адиабатического ядра является постоянство энтропии вдоль линии тока. При расчетах реальных задач сильные и слабые разрывы взаимодействуют друг с другом и с контактными разрывами, в результате чего может произойти необратимое накопление погрешностей из-за осцилляции и дистракции разрывов, что в итоге дает существенное различие между характеристиками реального физического процесса и его математического образа. Для адекватного описания поведения материалов необходима высокая точность как кинетических моделей разрушения, так и численных методов с оптимальными дистракцией и немонотонностью.
В случае идеальной сжимаемой жидкости без учета теплопроводности и девиатора тензора напряжений (вязкости, упругости, пластичности) система законов сохранения (0.1)-(0.3) имеет вид:
dV dt
dU dt
VdivU = 0,
fVgradP = 0,
(0.7)
ёГЕ+и dt Л
J
+ VdivPU = 0.
В случае плоских, сферически-симметричных или цилиндрически-симметричных движений, когда все величины зависят только от времени и одной лагранжевой координаты, уравнения (0.7) принимают вид
дМ
at
5V «Ф(«) 7(г"- и) = 0 (0.8)
аи
аф(ос)і
_, ap_
дМ
= 0,
(0.9)
U
д_ dt
Ґтл2 \
+ Е
V
дМ
+ a(P(a) (ra"lpU) = 0 (0.10)
Здесь М - массовая лагранжева координата, dM = 9(a)pdra, г - расстояние до
центра или оси симметрии, a =1,2,3, ф(1)=1, ф(2)=7і, ф(3)=4тс/3. Для замыкания к уравнениям (0.7) или (0.8)-(0.10) добавляется уравнение состояния
F(P,V,E) = 0. (0.11)
Положение каждой частицы в пространстве определяется ее эйлеровой координатой г = r(M,t), удовлетворяющей уравнению
(-)
I ft л u = o.
(0.12)
В случае одномерных симметричных движений справедливо также уравнение
ф(а)
ґдгаЛ
KdMj
v = o.
(0.13)
которое также может использоваться для нахождения г. Во многих численных методах применяются следствия из уравнений (0.8)-(0.10)
дЕ аф(а)ага" и_
at дм " (0.14)
5Е + Р = 0.
a
at
(0.15)
Одним из важнейших следствий системы законов сохранения (0.8)-(0.10) и уравнения состояния (0.11) является постоянство энтропии вдоль линии тока в адиабатических течениях. Уравнения (0.8)-(0.15) содержат три термодинамические функции Р, V и Е. Пусть V и Е независимы. В этом случае все остальные термодинамические функции, в том числе и на изэнтропе, зависят от V и Е. Уравнение скорости изменения энтропии S(V,E) вдоль линии тока имеет вид
as at
rдs av fas aE
av
aE;v at
—+ at
УЕ
(0.16)
Поскольку
fdS\
J\
УдЕ
T
as "І
av
/E
T
из (0.16) следует
a at at
(0.17)
Из (0.15), (0.17) следует уравнение сохранения энтропии вдоль линии тока (траектории частицы) идеальной среды
Т = 0.
at
(0.18)
Продифференцировав по t уравнение состояния (0.11), получим скорость изменения давления
ар
at
fap v 3VyE
av
dt
vdEyv
дЕ dt Подставив сюда (0.15), получим еще одно следствие из уравнения состояния и законов сохранения
+ а2 = 0, (0-19)
at at
где а
ар удУуЕ
+ р
( 1
кдЕ,
- квадрат скорости распространения звуковых
возмущений в лагранжевых (массовых) координатах. Из (0.8) и (0.19) следует
ар , аф(а)ага_1и
ь a v =
dt дМ (0.20)
Ha поверхности сильного разрыва система законов сохранения принимает вид условий Гюгонио-Ренкина
(0.21) (0.22) (0.23)
w(v+-v) = -(u+-u_),
W(U+-U_) = P+-P_, E+-E_+0.5(P++P_)(V+-V_ ) = 0,
„. dm
где W = - скорость распространения разрыва в лагранжевых (массовых)
dt
координатах, величины с индексом "-" характеризуют состояние вещества
перед линией разрыва, а с индексом "+" - за линией разрыва, U+, U. нормальные к поверхности разрыва компоненты вектора скорости U.
Условия на контактных границах (КГ) получаются из (0.21)-(0.23) в
частном случае, когда поток массы через поверхность разрыва отсутствует,
т.е. при W = 0. Тогда из (0.21)-(0.23) следует непрерывность интенсивных
величин - скорости и давления - на КГ:
U+=U_, Р+=Р_. (0.24)
При этом экстенсивные величины - удельный объем и удельная внутренняя энергия - могут оставаться разрывными.
В силу нелинейности уравнений газовой динамики их решение в общем случае можно найти лишь численно. Наиболее разработанным численным методом решения задач газодинамики является метод конечных разностей. В разностных методах непрерывные функции заменяются дискретными, определенными в узлах разностной сетки. Вообще говоря, для каждой функции может быть выбрана своя сетка, однако, во избежание дополнительных интерполяций при вычислении давления из уравнения состояния по V и Е, эти величины задаются на одной и той же сетке. После этого остается две возможности: 1) значения скорости определяются на той же сетке, что и давление; 2) значения скорости определяются на сетке, отличной от той, на которой определяется давление.
Разностная схема, вообще говоря, должна отражать основные свойства сплошной среды. Поэтому естественно требовать, чтобы в разностной схеме прежде всего выполнялись разностные аналоги законов сохранения. Разностные схемы, в которых изменение массы, количества движения и полной энергии в области интегрирования определяются только потоками массы через границы, импульсом и работой сил, действующих на границах, называются консервативными. На важность требования консервативности схемы обращали внимание многие исследователи. Так, например, в начале 50-х годов А.Н. Тихонов и А.А. Самарский [2] для обоснования интегро-интерполяционного подхода к конструированию разностных схем построили пример, когда неконсервативная разностная схема, обеспечивающая второй порядок точности в классе достаточно гладких коэффициентов, расходится в классе разрывных коэффициентов [2]. Однако требование консервативности не исчерпывает всех требований к разностной схеме. Дело в том, что в так называемом дивергентном уравнении энергии «сохраняется» только полная
энергия є = Е + 0.5U2. Поэтому погрешности в определении скорости, т.е. кинетической энергии, влияют на точность вычисления внутренней энергии, которая является суммарной величиной, состоящей из упругой (холодной) энергии, тепловой энергии ядер, тепловой энергии электронов, свободной энергии и т.д. Упомянутые выше требования консервативности оставляют без контроля переходы энергии из одной формы в другую, а это может исказить температуру, давление, энтропию, энтальпию и другие термодинамические величины. В [3] приведен пример, когда погрешности аппроксимации приводят к заметному искажению внутренней энергии.
Для изучения свойств разностных схем разностные уравнения чаще всего рассматриваются в дифференциальной форме. Вопросы получения разностных уравнений газовой динамики в дифференциальной форме и исследование свойств их дифференциальных приближений подробно изучены в [4]. В [5] показано, что для того, чтобы определить диссипативные свойства разностной схемы, нужно построить для нее уравнение производства энтропии и уравнение производства массы и исследовать остаточные члены для этих уравнений. Очевидно, что изменение энтропии из-за погрешностей аппроксимации не должно превосходить ее изменений в характерных физических процессах.
Конечноразностные методы расчета нестационарных течений сжимаемых сред основываются на системе законов сохранения либо в форме Эйлера, либо в форме Лагранжа. И лагранжевы, и эйлеровы методы имеют свои достоинства и недостатки. Выбор системы координат для расчета течения газа определяется постановкой задачи. Если важны параметры потока в заданной пространственной области (например, течение газа в газопроводе, задачи обтекания жесткой поверхности и т.д.), то естественно выбрать эйлеровы координаты. В связи с тем, что в этом случае сетка является неподвижной в пространстве, не возникают проблемы, связанные с сеткой. Однако, при расчете задач, связанных с течением определенной массы вещества, применение эйлеровых координат может привести либо к неоправданному уменьшению, либо к увеличению количества точек сетки и, следовательно, к потере точности численного решения. Например, при сильном сжатии вся рассчитываемая масса вещества может попасть в один счетный интервал эйлеровой сетки, что приводит к полной потере точности.
Чтобы обеспечить необходимую точность расчета центрированных волн разрежения в самом начале их существования, когда градиенты велики, С.К. Годунов предложил использовать подвижные сетки [6], В этом случае точки сетки, связанные с контактными границами или со слабыми разрывами, движутся вместе с ними. Промежуточные точки сетки получаются по произвольному закону с сохранением определенного минимума или максимума точек.
При использовании лагранжевых координат в задачах с большими деформациями в двумерной или трехмерной постановке возникают проблемы, связанные с перестройкой сетки. Но при необходимости детально исследовать газодинамические процессы, происходящие в некоторой фиксированной массе вещества, применение лагранжевых методов является наиболее целесообразным. В этом случае легко проследить историю деформирования частицы вещества, не возникает проблем с отслеживанием контактных границ, местами зарождения микроповреждений, зародышей новой фазы, что особенно важно для описания сложных процессов, связанных с деформациями и фазовыми переходами.
Область, в которой рассматривается движение вещества, разбивается сильными и слабыми разрывами на области гладкого течения, в которых выполняются законы сохранения в дифференциальной форме, тогда как на разрывах удовлетворяются условия совместности. Наиболее общим является решение, представляющее собой совокупность гладких решений, примыкающих друг к другу через линии сильных, слабых или контактных разрывов с соблюдением условий совместности. Такая структура решения естественно описывается методом характеристик [7], учитывающим, в принципе, все особенности решения. От других разностных методов его отличает аппроксимация не законов сохранения, а характеристических уравнений и многократное использование операторов интерполирования. Сглаживание профилей, характерное для разностных схем с фиксированной сеткой, является минимальным в методе характеристик, так как применяемая в нем сетка строится с учетом области зависимости решения. Альтернативой методу характеристик являются разностные методы с нулевой дистракцией, выделяющие особенности в решении, например, ударные волны, слабые и контактные разрывы. В этом случае для расчета параметров течения используются алгебраические и дифференциальные уравнения - законы сохранения и их следствия, записанные для каждого типа разрыва. Примером является неоднородный метод В.Ф. Куропатенко, реализованный в программе "Волна" [8], где для решения системы разностных уравнений используется регулярная сетка для областей интегрирования с гладкими решениями и "размазанными" особенностями и сетка особенностей, которая накладывается на регулярную сетку.
Наиболее простыми для реализации на ЭВМ являются методы с ненулевой дистракцией, получившие название "однородных". Для описания течения в ударном слое, заменяющем сильный разрыв, в этом случае в уравнения гидродинамики вводятся диссипативные члены. Первой схемой такого рода была схема Неймана-Рихтмайера [9], в которой по аналогии с физической вязкостью в уравнения газодинамики вводится псевдовязкость. При этом были использованы результаты Беккера [10], который в 1922 году показал, что введение физической вязкости в уравнения газодинамики приводит к появлению переходного слоя, толщина которого в газах сравнима с длиной пробега молекулы. В методе Неймана-Рихтмайера после введения псевдовязкости сильный разрыв заменяется переходным слоем конечной ширины. Другим вариантом однородных методов с псевдовязкостью, является метод, предложенный в США Лаксом и Вендрофом [11]. Лаке [12] предложил однородный метод с аппроксимационной «вязкостью». В СССР были созданы разностные схемы, в которых для диссипации энергии использовались уравнения физических процессов. Первой схемой такого рода была схема, предложенная С.К.Годуновым [13], в которой для расчета диссипации энергии применяются соотношения для распада произвольного разрыва. В отличие от схем Неймана-Рихтмайера и Лакса-Вендроффа, этот метод не содержит эмпирических констант. В методе, предложенном В.Ф. Куропатенко [14], [15] применяется разнородная аппроксимация на ударных волнах (УВ) и волнах разрежения (ВР). Такие схемы в области гладкого течения аппроксимируют уравнения газодинамики, а на разрывах - условия Гюгонио. Диссипация энергии учитывается при расчете вспомогательного давления. Применение специальных разностных уравнений в ячейках сетки, содержащих сильный разрыв, приводит к конечной дистракции, т.е к замене сильного разрыва ударным слоем, шириной в несколько счетных интервалов. Таким образом, рассчитываемая ударная волна заменяется конечным числом ударных волн меньшей интенсивности. Этот метод, также не имеющий эмпирических констант, был реализован в разностных схемах, в которых давление, плотность и энергия определялись, как правило, в серединах интервалов, а скорость в узлах разностной сетки (недивергентные схемы). Предпринимались также попытки создания дивергентных разностных схем расчета ударных волн, реализующих этот метод, например [15]. Однако во всех реализациях немонотонность зависела от числа Куранта и при произвольном соотношении шагов по времени и по пространству наблюдались осцилляции за фронтом УВ.
Введение диссипативных членов в разностные уравнения приводит с одной стороны к дистракции разрывов, а с другой стороны, к осцилляциям в решении. Как правило, рост дистракции сопровождается уменьшением амплитуды осцилляции и наоборот. Все разностные методы обладают различными дистракцией и осцилляционными свойствами. Так, например, метод Неймана-Рихтмайера дает сильные осцилляции как за фронтом УВ, так и в окрестности слабого разрыва на волне разрежения. Метод С.К. Годунова, являющийся монотонным, обладает сильной дистракцией слабых разрывов. При сложной картине течения, взаимодействии сильных и слабых разрывов друг с другом, а также с контактными разрывами влияние этих свойств разностного метода может накапливаться, что может быть причиной значительного отличия параметров численного решения от характеристик реального физического процесса. Особенно существенным это влияние становится при описании таких тонких эффектов, как откольное разрушение и зарождение фазовых переходов. При применении подавляющего большинства существующих разностных методов для расчета откольного разрушения масса и начальный импульс отколов определяются с низкой точностью. Поэтому представляется актуальной задача построения разностного метода, сочетающего в себе малые амплитуды осцилляции и малую дистракцию разрывов.
В диссертации излагается разностный метод расчета неустановившихся течений сжимаемой жидкости, который на волне сжатия в акустическом приближении удовлетворяет условиям теоремы С.К. Годунова и обеспечивает монотонность профилей, обладает небольшой дистракцией разрывов, имеет нулевую диссипацию в области, где справедливы законы сохранения в дифференциальной форме, и практически не дает энтропийных следов при выходе УВ на свободную поверхность.
Созданный метод излагается в одномерной постановке в случае плоской симметрии. Уравнения в случае сферической и цилиндрической симметрии не обсуждаются в данной работе, так как для описания подавляющего большинства экспериментов по откольному разрушению достаточно рассмотреть плоский случай. В [16] показана принципиальная возможность обобщения метода на двух- или трехмерный случай, однако это выходит за рамки данной работы.
На защиту выносятся:
1. Дивергентная разностная схема расчета ударных волн, реализующая метод В.Ф. Куропатенко, основанный на применении условий Гюгонио для определения характеристик среды в области ударных слоев.
2. Уравнения для определения вспомогательных величин - давления и скорости - в области непрерывных течений и их обоснование.
3. Исследование диссипации энергии, дистракции разрывов, устойчивости и монотонности предложенной разностной схемы. Обоснование отсутствия диссипации энергии на волнах разрежения и малости амплитуды осцилляции за разрывом.
4. Аналитическое решение задачи о выходе ударной волны на свободную поверхность конденсированной среды с последующим взаимодействием двух волн разрежения и откольным разрушением среды.
5. Сравнительный анализ точности разностных схем при решении задачи об откольном разрушении конденсированной среды.
Диссертация состоит из введения и четырех глав. В первой главе излагаются априорные методы исследования свойств разностных схем и выбираются критерии их сравнения. Во второй главе проводится теоретический анализ свойств методов расчета ударных волн и волн разрежения. Поскольку известно всего четыре принципиально различных механизма описания диссипации энергии в разностных методах расчета ударных волн, то исследуются только разностные схемы Неймана-Рихтмайера, Лакса, Годунова и Куропатенко. В третьей главе излагается новая разностная схема, исследуются ее свойства и проводится верификация численного метода. Четвертая глава посвящена аналитическому решению задачи с откольным разрушением и исследованию влияния свойств разностных схем и характерных погрешностей в окрестности сильных и слабых разрывов на моделирование откольного разрушения.
Основные результаты работы были доложены на 7 Международных, 1 Российской конференциях и 8м Всероссийском съезде по теоретической и прикладной механике ("New Models and Numerical Codes for Shock Wave Processes in Condensed Media", Oxford, UK, 1997r. [17], "Зимняя школа по механике сплошных сред", Пермь, Россия, 1997г. [18], "Ударные волны в конденсированных средах", Санкт-Петербург, Россия, 1998г. [19], V Забабахинские научные чтения, Снежинск, Россия, 1998г. [20], "Конечно-разностные методы: теория и приложения", Минск, Беларусь, 1998г. [21], "Ударные волны в конденсированных средах", Санкт-Петербург, Россия, 2000г. [22], VIII Всероссийский съезд по теоретической и прикладной механике, Пермь, 2001 г. [23], VI Забабахинские научные чтения, Снежинск, Россия, 2001 г. [24, 25]) и опубликованы в [16], [26], в трудах конференций [17], [19], [21] и в Препринте РФЯЦ-ВНИИТФ [27]. По разностной схеме «Кама» в РФЯЦ-ВНИИТФ были проведены серии методических расчетов и производственных расчетов по прониканию полидисперсной среды через пористую преграду (номер государственной регистрации темы 21292, шифр "Хранилище"). В Томском государственном университете были проведены расчеты разрушения и фазовых переходов в керамиках и результаты этих работ вошли в несколько защищенных кандидатских диссертаций (например [28]).
Я выражаю глубокую благодарность моему научному руководителю профессору В.Ф.Куропатенко за постоянное внимание, моральную поддержку и помощь при получении результатов диссертации. Я благодарю сотрудников сектора 332 РФЯЦ-ВНИИТФ И.А.Доровских и Ю.В.Кайгородцеву за помощь и поддержку при создании программы «КАМА», в которой реализована описанная ниже разностная схема, при проведении и обработке результатов расчетов.
Метод исследования диссипативных свойств разностной схемы
Область, в которой рассматривается движение вещества, разбивается сильными и слабыми разрывами на области гладкого течения, в которых выполняются законы сохранения в дифференциальной форме, тогда как на разрывах удовлетворяются условия совместности. Наиболее общим является решение, представляющее собой совокупность гладких решений, примыкающих друг к другу через линии сильных, слабых или контактных разрывов с соблюдением условий совместности. Такая структура решения естественно описывается методом характеристик [7], учитывающим, в принципе, все особенности решения. От других разностных методов его отличает аппроксимация не законов сохранения, а характеристических уравнений и многократное использование операторов интерполирования. Сглаживание профилей, характерное для разностных схем с фиксированной сеткой, является минимальным в методе характеристик, так как применяемая в нем сетка строится с учетом области зависимости решения. Альтернативой методу характеристик являются разностные методы с нулевой дистракцией, выделяющие особенности в решении, например, ударные волны, слабые и контактные разрывы. В этом случае для расчета параметров течения используются алгебраические и дифференциальные уравнения - законы сохранения и их следствия, записанные для каждого типа разрыва. Примером является неоднородный метод В.Ф. Куропатенко, реализованный в программе "Волна" [8], где для решения системы разностных уравнений используется регулярная сетка для областей интегрирования с гладкими решениями и "размазанными" особенностями и сетка особенностей, которая накладывается на регулярную сетку.
Наиболее простыми для реализации на ЭВМ являются методы с ненулевой дистракцией, получившие название "однородных". Для описания течения в ударном слое, заменяющем сильный разрыв, в этом случае в уравнения гидродинамики вводятся диссипативные члены. Первой схемой такого рода была схема Неймана-Рихтмайера [9], в которой по аналогии с физической вязкостью в уравнения газодинамики вводится псевдовязкость. При этом были использованы результаты Беккера [10], который в 1922 году показал, что введение физической вязкости в уравнения газодинамики приводит к появлению переходного слоя, толщина которого в газах сравнима с длиной пробега молекулы. В методе Неймана-Рихтмайера после введения псевдовязкости сильный разрыв заменяется переходным слоем конечной ширины. Другим вариантом однородных методов с псевдовязкостью, является метод, предложенный в США Лаксом и Вендрофом [11]. Лаке [12] предложил однородный метод с аппроксимационной «вязкостью». В СССР были созданы разностные схемы, в которых для диссипации энергии использовались уравнения физических процессов. Первой схемой такого рода была схема, предложенная С.К.Годуновым [13], в которой для расчета диссипации энергии применяются соотношения для распада произвольного разрыва. В отличие от схем Неймана-Рихтмайера и Лакса-Вендроффа, этот метод не содержит эмпирических констант. В методе, предложенном В.Ф. Куропатенко [14], [15] применяется разнородная аппроксимация на ударных волнах (УВ) и волнах разрежения (ВР). Такие схемы в области гладкого течения аппроксимируют уравнения газодинамики, а на разрывах - условия Гюгонио. Диссипация энергии учитывается при расчете вспомогательного давления. Применение специальных разностных уравнений в ячейках сетки, содержащих сильный разрыв, приводит к конечной дистракции, т.е к замене сильного разрыва ударным слоем, шириной в несколько счетных интервалов. Таким образом, рассчитываемая ударная волна заменяется конечным числом ударных волн меньшей интенсивности. Этот метод, также не имеющий эмпирических констант, был реализован в разностных схемах, в которых давление, плотность и энергия определялись, как правило, в серединах интервалов, а скорость в узлах разностной сетки (недивергентные схемы). Предпринимались также попытки создания дивергентных разностных схем расчета ударных волн, реализующих этот метод, например [15]. Однако во всех реализациях немонотонность зависела от числа Куранта и при произвольном соотношении шагов по времени и по пространству наблюдались осцилляции за фронтом УВ.
Введение диссипативных членов в разностные уравнения приводит с одной стороны к дистракции разрывов, а с другой стороны, к осцилляциям в решении. Как правило, рост дистракции сопровождается уменьшением амплитуды осцилляции и наоборот. Все разностные методы обладают различными дистракцией и осцилляционными свойствами. Так, например, метод Неймана-Рихтмайера дает сильные осцилляции как за фронтом УВ, так и в окрестности слабого разрыва на волне разрежения. Метод С.К. Годунова, являющийся монотонным, обладает сильной дистракцией слабых разрывов. При сложной картине течения, взаимодействии сильных и слабых разрывов друг с другом, а также с контактными разрывами влияние этих свойств разностного метода может накапливаться, что может быть причиной значительного отличия параметров численного решения от характеристик реального физического процесса. Особенно существенным это влияние становится при описании таких тонких эффектов, как откольное разрушение и зарождение фазовых переходов. При применении подавляющего большинства существующих разностных методов для расчета откольного разрушения масса и начальный импульс отколов определяются с низкой точностью. Поэтому представляется актуальной задача построения разностного метода, сочетающего в себе малые амплитуды осцилляции и малую дистракцию разрывов. В диссертации излагается разностный метод расчета неустановившихся течений сжимаемой жидкости, который на волне сжатия в акустическом приближении удовлетворяет условиям теоремы С.К. Годунова и обеспечивает монотонность профилей, обладает небольшой дистракцией разрывов, имеет нулевую диссипацию в области, где справедливы законы сохранения в дифференциальной форме, и практически не дает энтропийных следов при выходе УВ на свободную поверхность.
Созданный метод излагается в одномерной постановке в случае плоской симметрии. Уравнения в случае сферической и цилиндрической симметрии не обсуждаются в данной работе, так как для описания подавляющего большинства экспериментов по откольному разрушению достаточно рассмотреть плоский случай. В [16] показана принципиальная возможность обобщения метода на двух- или трехмерный случай, однако это выходит за рамки данной работы.
Недивергентная разностная схема В.Ф.Куропатенко
Известно большое количество разностных схем, использующих для описания диссипации энергии искусственную вязкость (линейную [29], [30], линейно-квадратичную [31], [32], в неявной форме [33]).
В методе П. Лакса [12] диссипация энергии в ударном слое определяется специальной формой аппроксимации дифференциальных законов сохранения разностными. В результате погрешности аппроксимации образуют "аппроксимационную вязкость", обеспечивающую диссипацию энергии в ударном слое. Этот метод, как и метод Неймана-Рихтмайера, допускает большое количество форм [4], [12], [34], [35], [36].
В методе С.К. Годунова [13] в основе механизма диссипации энергии лежит оригинальная физическая идея. Все сеточные функции предполагаются кусочно-постоянными. В таком случае на гранях сеточных ячеек возникают разрывы, которые являются произвольными и распадаются с образованием различных конфигураций устойчивых разрывов (сильных, слабых, контактных). Полученные при этом значения давления и скорости на контактном разрыве используются в качестве вспомогательных величин на гранях сеточных ячеек, соответствующих постоянным значениям т. Все разностные уравнения записываются в дивергентной форме. Уравнения распада произвольного разрыва при заданных исходных состояниях имеют единственное решение. В этом смысле метод С.К. Годунова не допускает различных версий. Тем не менее, на практике для повышения экономичности применяются приближенные уравнения для решения задачи о распаде произвольного разрыва. Это порождает многочисленные версии метода Годунова, например [37].
В методе В.Ф. Куропатенко [14] для расчета величин в ударном слое привлекаются уравнения на сильном разрыве (0.21)-(0.23). Величины за сильным разрывом (со знаком «+» в (0.21)-(0.23)) используются в качестве вспомогательных сеточных величин. Метод допускает применение на сетках разных типов и с использованием разных уравнений ударной адиабаты [15], [38].
Перечисленные методы исчерпывают механизмы диссипации энергии в ударном слое. Каждый из перечисленных методов расчета ударных волн допускает различные версии в зависимости от сетки, определения основных и вспомогательных величин, выбора величин-параметров. Методов описания диссипации энергии всего четыре, разностных схем, реализующих эти методы - много. При расчетах по таким схемам сильные разрывы выделяются автоматически, как слои быстрого изменения величин, при переходе через которые параметры решения претерпевают изменения, близкие к тем, которые происходят при переходе через фронт разрыва, причем переходные слои движутся со скоростями, близкими к скоростям поверхностей разрывов.
При построении разностной схемы получается некоторая дискретная модель непрерывной среды, которая, вообще говоря, не адекватна исходной сплошной среде и обладает собственными диссипацией и дисперсией. Наличие собственных внутренних свойств разностной схемы может привести к ощутимым искажениям разностного решения, проявляющимся в виде не имеющих физического смысла осцилляции, дистракции разрывов, диссипации энергии и так далее. Очевидно, что чем меньше ширина переходных слоев, которыми заменяются сильные разрывы, тем точнее описывается амплитуда затухающей ударной волны. Немонотонность и дистракция разностного метода влияют на точность определения мест зарождения трещин или фазовых переходов.
При построении разностной схемы, аппроксимирующей некоторую дифференциальную задачу, бесконечномерное пространство функций непрерывного аргумента заменяется конечномерным пространством сеточных функций, а дифференциальное уравнение - системой алгебраических соотношений. То, что решение дифференциальной задачи и сеточное решение принадлежат разным функциональным пространствам, затрудняет исследование свойств разностных схем. Поэтому часто рассматривают разностные операторы в том же функциональном пространстве, что и аппроксимируемые дифференциальные операторы, то есть считают, что разностные уравнения удовлетворяются функциями непрерывного аргумента в каждой точке рассматриваемой области, а не только в узлах сетки. При таких предположениях получаются дифференциальные приближения разностных схем. Для этого выполняются те же самые операции, которые совершают при определении погрешности аппроксимации. Левая часть полученного таким образом уравнения совпадает с левой частью дифференциального уравнения, а в правой части » стоит погрешность аппроксимации. Вопросы получения разностных уравнений газовой динамики в дифференциальной форме и исследования их свойств подробно изучены в [4] и [5].
Система уравнений газовой динамики (0.8)-(0.10), (0.12)-(0.15), (0.18)-(0.20), содержащая 10 уравнений и 7 функций (г, U и термодинамические функции, из которых только две независимы), является переопределенной. Независимых искомых функций в ней только четыре. Конструирование разностной схемы, заключается в выборе четырех независимых разностных уравнений системы, что приводит к определению четырех конкретных независимых погрешностей аппроксимации.
Анализ устойчивости разностной схемы для ударной волны
Задача была рассчитана по методике КАМА и методикам [9],[44],[13],[14]. На Рис. 4.7, Рис. 4.8, Рис. 4.9, Рис. 4.10, Рис. 4.11 приведены зависимости Pmin от массовой координаты, соответствующие расчетам по этим методикам. Зависимость Pmjn(m), полученная при расчете по методу КАМА (Рис. 4.10) удовлетворительно согласуется с аналитическим решением. Видно, что при расчете по монотонному методу [13], минимальное давление не достигает предела прочности, а при расчетах по методикам [9], [44], [14] разрушение наступает раньше, чем в аналитическом решении. В таблице 2 приведены массы отколовшегося вещества, полученные в расчетах по этим разностным схемам, а также абсолютные и относительные погрешности по сравнению с аналитическим решением, полученным в 3. Погрешность определения массы откола в расчете по разностной схеме «Кама» значительно меньше, чем при расчете по разностным схемам [9], [13], [14], [44].
Предложена разностная схема расчета неустановившихся течений сплошных сред, относящаяся к классу схем, не содержащих эмпирических констант при описании диссипации энергии в ударном слое, заменяющем сильный разрыв. Разностная схема реализует метод В.Ф.Куропатенко, основанный на применении условий Гюгонио для определения характеристик среды в области ударных слоев и на интегрировании с наперед заданной точностью уравнений изэнтроп в случае адиабатических течений. Исследована устойчивость, аппроксимация, диссипация энергии, дистракция и монотонность новой разностной схемы. Для обоснования ее достоинств проведены аналогичные исследования разностных схем Неймана-Рихтмайера, Лакса, Годунова и недивергентной разностной схемы Куропатенко.
Построено аналитическое решение задачи с откольным разрушением конденсированной среды после выхода ударной волны на свободную поверхность. Специально подобранное граничное условие позволяет реализовать режим откола только в одной точке среды.
Создана методическая программа «Кама», в которой реализована предложенная разностная схема. Расчеты по программе показывают, что разностная схема «Кама» сочетает в себе малую амплитуду осцилляции и малую дистракцию разрывов, что позволяет наиболее точно описывать взаимодействие разрывов различных типов друг с другом и, следовательно, поведение веществ при откольном разрушении. Применение условий Гюгонио для описания диссипации энергии в области ударного слоя и выбранный вид разностных уравнений позволяют избежать появления энтропийных следов при выходе ударной волны на свободную поверхность вещества. 4. Проведен сравнительный анализ точности разностных схем при решении задачи об откольном разрушении конденсированной среды. Показано, что погрешность определения массы откола в расчете по предложенной в диссертации разностной схеме «Кама» значительно меньше, чем при расчете по разностным схемам [9], [13], [14], [44].
Взаимодействие двух волн разрежения с образованием откола. Аналитическое решение
При сложной картине течения, взаимодействии сильных и слабых разрывов друг с другом, а также с контактными разрывами влияние этих свойств разностного метода может накапливаться, что может быть причиной значительного отличия параметров численного решения от характеристик реального физического процесса. Особенно существенным это влияние становится при описании таких тонких эффектов, как откольное разрушение и зарождение фазовых переходов. При применении подавляющего большинства существующих разностных методов для расчета откольного разрушения масса и начальный импульс отколов определяются с низкой точностью. Поэтому представляется актуальной задача построения разностного метода, сочетающего в себе малые амплитуды осцилляции и малую дистракцию разрывов. В диссертации излагается разностный метод расчета неустановившихся течений сжимаемой жидкости, который на волне сжатия в акустическом приближении удовлетворяет условиям теоремы С.К. Годунова и обеспечивает монотонность профилей, обладает небольшой дистракцией разрывов, имеет нулевую диссипацию в области, где справедливы законы сохранения в дифференциальной форме, и практически не дает энтропийных следов при выходе УВ на свободную поверхность.
Созданный метод излагается в одномерной постановке в случае плоской симметрии. Уравнения в случае сферической и цилиндрической симметрии не обсуждаются в данной работе, так как для описания подавляющего большинства экспериментов по откольному разрушению достаточно рассмотреть плоский случай. В [16] показана принципиальная возможность обобщения метода на двух- или трехмерный случай, однако это выходит за рамки данной работы.
На защиту выносятся: 1. Дивергентная разностная схема расчета ударных волн, реализующая метод В.Ф. Куропатенко, основанный на применении условий Гюгонио для определения характеристик среды в области ударных слоев. 2. Уравнения для определения вспомогательных величин - давления и скорости - в области непрерывных течений и их обоснование. 3. Исследование диссипации энергии, дистракции разрывов, устойчивости и монотонности предложенной разностной схемы. Обоснование отсутствия диссипации энергии на волнах разрежения и малости амплитуды осцилляции за разрывом. 4. Аналитическое решение задачи о выходе ударной волны на свободную поверхность конденсированной среды с последующим взаимодействием двух волн разрежения и откольным разрушением среды. 5. Сравнительный анализ точности разностных схем при решении задачи об откольном разрушении конденсированной среды. Диссертация состоит из введения и четырех глав. В первой главе излагаются априорные методы исследования свойств разностных схем и выбираются критерии их сравнения. Во второй главе проводится теоретический анализ свойств методов расчета ударных волн и волн разрежения. Поскольку известно всего четыре принципиально различных механизма описания диссипации энергии в разностных методах расчета ударных волн, то исследуются только разностные схемы Неймана-Рихтмайера, Лакса, Годунова и Куропатенко. В третьей главе излагается новая разностная схема, исследуются ее свойства и проводится верификация численного метода. Четвертая глава посвящена аналитическому решению задачи с откольным разрушением и исследованию влияния свойств разностных схем и характерных погрешностей в окрестности сильных и слабых разрывов на моделирование откольного разрушения.
Основные результаты работы были доложены на 7 Международных, 1 Российской конференциях и 8м Всероссийском съезде по теоретической и прикладной механике ("New Models and Numerical Codes for Shock Wave Processes in Condensed Media", Oxford, UK, 1997r. [17], "Зимняя школа по механике сплошных сред", Пермь, Россия, 1997г. [18], "Ударные волны в конденсированных средах", Санкт-Петербург, Россия, 1998г. [19], V Забабахинские научные чтения, Снежинск, Россия, 1998г. [20], "Конечно-разностные методы: теория и приложения", Минск, Беларусь, 1998г. [21], "Ударные волны в конденсированных средах", Санкт-Петербург, Россия, 2000г. [22], VIII Всероссийский съезд по теоретической и прикладной механике, Пермь, 2001 г. [23], VI Забабахинские научные чтения, Снежинск, Россия, 2001 г. [24, 25]) и опубликованы в [16], [26], в трудах конференций [17], [19], [21] и в Препринте РФЯЦ-ВНИИТФ [27]. По разностной схеме «Кама» в РФЯЦ-ВНИИТФ были проведены серии методических расчетов и производственных расчетов по прониканию полидисперсной среды через пористую преграду (номер государственной регистрации темы 21292, шифр "Хранилище"). В Томском государственном университете были проведены расчеты разрушения и фазовых переходов в керамиках и результаты этих работ вошли в несколько защищенных кандидатских диссертаций (например [28]).
Я выражаю глубокую благодарность моему научному руководителю профессору В.Ф.Куропатенко за постоянное внимание, моральную поддержку и помощь при получении результатов диссертации. Я благодарю сотрудников сектора 332 РФЯЦ-ВНИИТФ И.А.Доровских и Ю.В.Кайгородцеву за помощь и поддержку при создании программы «КАМА», в которой реализована описанная ниже разностная схема, при проведении и обработке результатов расчетов.
