книги / Метод крупных частиц в газовой динамике
..pdf$ 4] |
|
КОНСЕРВАТИВНЫЙ МЕТОД ЧАСТИЦ В ЯЧЕЙКАХ |
8 1 |
В |
связи |
с тем, что в (3.38) параметры на теле находятся |
с точностью до* |
0 (h a), |
ti> \ |
(это следует из способа построения разностных граничных усло |
вий), внутри граничной ячейки может быть определено положение поверхности, тела. Таким образом, размеры дробной ячейки находятся, исходя из значений газодинамических параметров в потоке и в фиктивных граничных ячейках.
Отсюда следует, что для задач со вдувом при надлежащей постановке краевых условий ячейки, через которые реализуется вдув (отсос) струи с поверхности тела, могут рассчитываться по формулам для целых ячеек. При, этом на поверхности, соответствующей границе тела, будут с необходимой точностью выполняться граничные условия (3.38).
Аналогичные рассуждения можно провести и для поверхности тела, на которой формулируются условия другого типа. Например, для условий при липания граничные условия следует ставить так, чтобы в любой точке поверх ности тела выполнялось
w = 0(/i2), v = 0(h*),
а для условий непротекания —
Wn = О (/I2), WX = WTQ+ 0 (/ia) и т. д.
Так моделируется граница тела. Разностные краевые условия для вели чин с «волнами» конструируются при этом с учетом расщепления. После по становки граничных условий производится сквозной счет по единым формулам для целых ячеек во всех точках расчетной сетки (в том числе и на границе о телом).
Таким образом, новая методика расчета граничных ячеек заключается в- надлежащей формулировке краевых условий для дробных ячеек и в расчете по следних по единому алгоритму для целых ячеек, приведенному в § 1 настоящей главы. При этом примерно вдвое сокращается объем программы и заметно* уменьшается время счета. Такой подход кажется особенно оправданным при рассмотрении тел сложной, формы и для пространственно-трехмерных задач. [213, 216, 218, 220 и др.].
§4. Консервативный метод частиц в ячейках
1.На основе метода крупных частиц можно построить, как частный случай,, численные схемы и для дискретной модели частиц [23].
Вконсервативном методе частиц [436] сплошная среда моделируется; совокупностью точечных частиц фиксированной массы, обладающих опреде ленными значениями импульса и энергии. Такое множество частиц представ ляет собой лагранжеву сетку. Параллельно с ней исследуемую область по крывают неподвижной эйлеровой сеткой — на ней и определяются поля всех
газодинамических величин. Лагранжевы частицы движутся, как это обычно принято, через неподвижную эйлерову сетку. При этом считается, что пока частица находится в какой-либо эйлеровой ячейке, ее параметры определяются соответственными газодинамическими параметрами данной ячейки. Ограни чимся здесь описанием основной идеи такого подхода.
Пусть каждая частица обладает фиксированной массой т , импульсом m W и полной энергией Е. В качестве эйлеровой сетки примем прямоугольное раз
биение области с шагами Ах, Ау* Поскольку сплошная среда заменена ансамблем частиц, обладающих
определенной массой, то уравнение неразрывности из рассмотрения исклю чается, так как закон сохранения массы выполняется автоматически. Каждый расчетный цикл, осуществляемый за время Aty разбивается, как обычно, на три этапа путем расщепления исходных дифференциальных уравнений.
S 2 |
МЕТОД КРУПНЫХ ЧАСТИЦ |
[ГЛ. III |
Э й л е р о в |
э т а п . Разностные формулы этого шага |
полностью совпа |
дают со схемами эйлерова этапа метода крупных частиц. Здесь вычисляются промежуточные параметры потока в предположении отсутствия потока массы частиц через границы эйлеровых ячеек. Поле плотности р при этом предпола гается «замороженным». Следовательно, изменение газодинамических величин происходит только за счет сил давления, поэтому в (3.1) опускаются все члены
вида div(\|)pW), где ф ={и, |
и, |
Е}. |
частиц и |
|
Л а г р а н ж е в |
э т а п . |
На этом шаге происходит движение |
||
вычисляются потоки |
массы |
АМ п через границы эйлеровых ячеек. |
За время |
At каждая частица, находящаяся в момент tn в ячейке (i, /), пройдет расстояние
Аг = (/•« ?,/+ /•» ?,/) At.
Новые координаты частицы теперь определятся так:
|
|
|
X n + 1 = X n -f- |
j At, |
0gv |
|
|
|
y n+1 = y n + Vit j A t. |
* ' |
|
.Для |
повышения запаса |
устойчивости в (3.39) могут использоваться скорости |
|||
/7? /» |
и?,/» |
полученные |
на эйлеровом этапе. Сравнивая координаты частиц |
||
<3’.39) с |
координатами |
границ эйлеровых |
ячеек |
|
x ~ i - Ах,
(3.40)
y j= i - b y ,
можно определить, остались ли частицы в прежней ячейке или перешли в другую (и в какую именно).
Поток массы AM '1 через границу эйлеровой ячейки (3.40) определяется
•числом частиц, пересекших эту границу за время At: ДМя = 2=Ьт 1*- Здесь знак
i
плюс относится к частицам, пришедшим в данную ячейку извне, знак минус —
кчастицам, ушедшим из нее.
За к л ю ч и т е л ь н ы й э т а п . Данный шаг консервативного метода
частиц полностью совпадает с соответствующим этапом метода крупных
Рис. 3.6. Задача о распаде разрыва (расчет методом крупных частиц, /=15AQ.
частиц. На основании законов сохранения вычисляются окончательные значе ния параметров потока в новый момент времени t+ A t. При этом предпола гается, что, переходя из одной ячейки в другую, частица переносит свои зна чения массы, импульса и энергии.
КОНСЕРВАТИВНЫЙ МЕТОД ЧАСТИЦ В ЯЧЕЙКАХ |
83 |
2. Рассмотрим в качестве примера модельную задачу о распаде ударной волны. Пусть в начальный момент времени t=Q газ покоится, правое полупро странство содержит газ с нормальной плотностью рп и давлением /?п, а левое полупространство заполнено газом с повышенной плотностью рл и давлением Рл (рис. 3.6, а). Параметры газа задаются следующим образом:
Рл/Рп = 5, рл/Рп = 5, Рп = 1» |
«л = Ип=0. |
|
Газ во всем пространстве считается идеальным |
с уравнением состояния р = |
|
= (х—1)р2. Пунктиром здесь и на последующих рисунках |
показано распре |
|
деление плотности в момент /= 0 . |
|
|
t^ZAt
х/Ах |
4) |
х/Ах |
|
|
1Рис. 3.7. Задача о распаде разрыва (расчет консервативным методом частиц).
Рис. 3.8. Задача о распаде разрыва (расчет консервативным методом частиц, /=100Д/).
Для решения одномерной задачи о распаде разрыва возьмем эйлерову сетку, состоящую из 50 ячеек. На обоих ее концах зададим условия свободного вытекания (втекания). Через некоторое время отчетливо наблюдается образо-
34 МЕТОД КРУПНЫХ ЧАСТИЦ [ГЛ. III
вание ударной волны и контактного разрыва (движущихся направо), а также волны разрежения, движущейся налево.
На рис. 3.6 приведены профили плотности и скорости в момент t/A t= 15, полученные методом крупных частиц. Из поведения профиля скорости на рис. 3.6, б видно, например, что при переходе через контактную поверхность скорость не изменяется, что и следовало ожидать.
Рис. 3.7—3.8 иллюстрируют расчеты, выполненные консервативным методом дискретных частиц. Здесь показана динамика распада разрыва: на
рис. 3.7(1) — 3.7(5) |
даны |
профили плотности |
р в |
последовательные моменты |
||
времени t/At= O y 2, |
4, 8, |
16, 32, 64 и 78; на рис. |
3.8 — профили плотности р |
|||
(рис. 3.8, а) и скорости и |
(рис. 3.8, б) в момент /= 100А/. Видно, что начинает |
|||||
образовываться |
аналогичная |
конфигурация: |
ударная волна — контактный |
|||
разрыв — волна |
разрежения. |
Однако, благодаря |
конечному числу частиц, |
здесь возникают флуктуации газодинамических величин. Их можно умень шить, если увеличивать число частиц в ячейке (ограничителем является объем машинной памяти). В расчетах, результаты которых приведены на рис. 3.7— 3.8, в начальный момент времени в левых 25 ячейках использовалось по 20 ча стиц в каждой ячейке, а в правых — по 4 частицы. Все частицы считались одинаковыми — масса каждой из них т = 0 ,2 5 .
Указанный подход, использующий дискретную модель сплошной среды, позволяет, таким образом, на основе метода крупных частиц получить кон сервативную версию метода частиц в ячейках.
Г Л А В А IV
ИССЛЕДОВАНИЕ ЧИСЛЕННЫХ СХЕМ МЕТОДА КРУПНЫХ ЧАСТИЦ
Займемся теперь исследованием полученных разностных схем метода крупных частиц — изучим здесь вопросы аппроксимации, образования меха низма диссипации (вязкостные эффекты схемы), устойчивости вычислений и т. п. Для этих целей весьма эффективным оказалось использование так назы ваемых дифференциальных приближений (д. п.) конечно-разностных уравяений *).
Для одномерных квазилинейных уравнений гиперболического типа ис следования в этой области были проведены Н. Н. Яненко и Ю. И. Шокиным 113, 27, 58—61, 235 и др.]. Как показал К- Хирт [25], вполне оправдано при менение метода дифференциальных приближений и для изучения свойств некоторых разностных схем, аппроксимирующих нелинейные уравнения (и, в частности, уравнения газовой динамики). Здесь строгого математического обоснования пока не имеется, и эта теория носит скорее эвристический харак тер, хотя практические расчеты показали ее эффективность [29, 62, 220].
Заметим, что с помощью использования аппарата дифференциальных приближений можно управлять свойствами консервативности [374, 350], ввести понятие инвариантности разностных схем [235, 354 и др.] и т. д.
Путем последовательного рассмотрения нулевого, первого и второго дифференциальных приближений оказалось возможным изучить все основные свойства полученных разностных; схем метода крупных частиц и выбрать наи более рациональные подходы для их конструирования [1, 22, 23, 29]. Деталь ный обзор состояния и развития метода дифференциальных приближений, а также анализ схем метода крупных частиц содержится в работах Ю. М. Да выдова [62, 213, 216—221, 301, 345—353 и др.].
§1. Аппроксимация уравнений
1.Введем здесь понятие первого дифференциального приближения, как это сделано в работах Н. Н. Яненко и Ю. И. Шокина [59—61 и др.].
Пусть Л(£, Ху т, hy Т) — разностный оператор, аппроксимирующий диф
ференциальный оператор J& (tt х, *2>). Здесь
&>i}\ Т = { Т 0у'П}у <2>0 = d/dt\ Я>х = {д!дх1% . . %д!дхж)\ Т ^ { Т Ху...у T J ;
Т 0— оператор сдвига по /, Ту— оператор сдвига по Xj\ т, h= {hlf . . . , hs\ — •сеточные параметры.
*) Для простейших уравнений алгоритм первого дифференциального приближения впер вые был предложен в 50-х годах А. И. Жуковым 127]. Работа Ф. X а р л о у (Н а г 1 о w F. Н. Stability of Difference Equations. Selected Topics.—Los Alamos Scientific Lab. Rept. № LAMS — 2452.—Los Alamos: I960) была, по-видимому, одной из первых попыток применения первого дифференциального приближения разностной схемы, аппроксимирующей дифференциальные уравнения, для исследования различных свойств этой схемы.
86 |
ИССЛЕДОВАНИЕ ЧИСЛЕННЫХ СХЕМ МЕТОДА КРУПНЫХ ЧАСТИЦ |
[ГЛ. IV |
Воператорном представлении будем иметь
**
|
|
|
|
7\, =<**•, |
Tj — e J |
dxi |
|
|
|
||||
Разложим оператор A (t, |
х, т, A, Т )—А (t, х, т, А, ет55°, ewS*) в ряд по парамет- |
||||||||||||
рам т, h. Тогда, вводя зависимость h=h(%)> получим |
|
|
|||||||||||
где |
|
A(tf, х, т, |
hy |
T )= J?(ty |
Ху т, |
hy £D) + Ry |
|||||||
JV(ty Ху |
т, |
hy |
|
|
|
х, |
2))+ т aP a (ty |
Ху ® ), |
|||||
|
|
|
|
||||||||||
|
R ( t , х, т, |
А, |
@>)= |
У, |
трР р (/, |
дг, |
® ). |
|
|
||||
Здесь т<*Ра (г, х, |
*2)) — первый за |
J2? (t, х, ^)) |
ненулевой |
член разложения. |
|||||||||
Оператор |
g |
(/, х, |
т, |
ht Щ |
будет называться первым дифференциальным |
||||||||
приближением (сокращенно п. д. п.) |
разностного оператора А(?, х, т, hy Т)> |
||||||||||||
а уравнение |
J3? (t, |
х, |
т, |
h, |
&)) и = 0 — п. д. |
п. |
разностного уравнения |
||||||
A(ty Ху т, к, |
Г )= 0 . |
|
|
|
|
|
|
\ |
|
|
|
|
|
Рассмотрим |
задачу |
Коши |
для гиперболической |
системы уравнений пер |
вого порядка с m неизвестными функциями и s независимыми пространствен
ными переменными [13, |
61, |
235] |
|
|
|
|
|
|
|
S |
|
|
|
|
|
% = |
Е |
А к (х, |
t) — -, |
и(х, |
0) = и0(х). |
(4.1) |
|
a |
|
k = , |
|
од* |
|
|
|
Здесь и = и ( Х у t) — вектор-функция с m |
компонентами, х = {х ь |
xs} — |
|||||
точка вещественного s-мерного пространства, |
ЛА(х, t) — веществен |
||||||
ные тхга-матрицы. |
|
|
|
|
|
|
|
Аппроксимируем систему (4.1) разностной схемой |
|
||||||
|
|
w'1+1 (х) = 'ЕВаиЛ(х + |
тХа), |
(4.2) |
|||
|
|
|
|
а |
|
|
|
где Ва= В а (Ху t) — вещественные m X m |
матрицы, /=/гт, ^ а = { ^ а » |
•» ^4} — |
|||||
векторы смещения, |
un(x)= u(xt пт). |
|
|
|
|||
При заданных ка назовем разностную схему (4.2) устойчивой, если вы |
|||||||
полнено следующее |
неравенство |
[27]: |
|
|
|
|и » « ( * ) |< [ 1 + 0 ( т ) ] |« » ( * ) |.
Норма понимается в смысле пространства L2.
Порядком точности разностной схемы (4.2) называется такое наибольшеечисло 5*, для которого все решения системы уравнений (4.1), имеющие непре рывные производные до (5*+1)-го порядка, удовлетворяют уравнению (4.2)' с точностью до. величин порядка 0 ( т ^ -1).
Если разностная схема (4.2) аппроксимирует уравнения (4.1) с первым порядком точности, то должны выполняться условия совместности
2 В а = Л УУкВа = А к (А = 1........ S), (4.3).
аа
где / — единичная матрица.
Условия (4.3) имеют весьма наглядную интерпретацию. Так, первоеиз них означает, что константа является решением данной системы уравнений, и т. д.
Если а принимает конечное число значений, то схема называется явнойу.
впротивном случае — неявной.
Вдальнейшем цри получении дифференциальных приближений будемиспользовать обычное в вычислительной математике предположение о том,.
$ 1] |
АППРОКСИМАЦИЯ УРАВНЕНИЙ |
87 |
|
что рассматриваемая функция имеет столько производных, сколько необхо димо по ходу рассуждений [222].
Гиперболическая форма (Т-форма) п. д. п . схемы (4.2) имеет вид
ь г и |
S |
дЧ |
|
|
|
= Е |
N. |
(4.4) |
|||
1 F |
' i k dXjdXf, |
||||
где |
*=1 |
|
|
|
|
|
|
|
|
||
Р „ - 2 * А в . , M - 1 ( - % + ± A , £ ) . |
|
||||
Г-форма п. д. п. получается из |
(4.2) при |
помощи разложения |
функций |
и п+1(х), ип(х+т%а) в ряды по параметру т в точке (х , t) и отбрасывания чле нов выше первого порядка по т.
Наряду с гиперболической рассматривается параболическая форма (17-
форма) п . д. л. |
|
|
|
i,k= i |
у * |
/= 1 |
(4.5) |
1 |
|||
-где |
|
|
|
т / сМ,- |
v^i |
d-Д/ \ |
|
DJ (х’ О = Л/ — |
|
у> |
|
<7* =у (Z |
|
*J = у |
*). |
|
Ее получают при |
рассмотрении Г-формы п. д. п. (4.4) на решении исходной |
|||
системы, т. е. используя соотношение |
|
|||
д 2и |
у |
+ |
± ( % + ± * , Щ |
) § г , - |
<эг2 ~~ |
A l A l‘ d x J d x lt |
|||
|
/.*= I |
|
/= 1 |
|
Из (4.5) видим, что разностная схема (4.2) добавляет к исходной системе |
||||
уравнений (4.1) вязкие члены |
|
|
|
|
|
Z |
c |
д 2и |
|
|
j kdxjdxk' |
|
||
|
/,Л=1 |
|
|
где С — матрица аппроксимационной вязкости, структура которой и опреде ляет диссипативные свойства схемы.
Практика показала, что можно рассматривать не все члены матриц Л, С, D. Обычно достаточно ограничиться исследованием диагональных членов матрицы аппроксимационной вязкости, которые и определяют в основном диссипативные свойства рассматриваемой схемы. При этом упрощении зна чительно сокращается объем выкладок, результаты же исследований для дан ного класса задач фактически не меняются. Такое упрощение мы и будем ис
пользовать ниже.
Можно получить дифференциальные приближения и более высокого порядка. Для этого каждый член разностной схемы представляется как глад кая функция своих аргументов и проводится разложение в ряды Тейлора по
•сеточным параметрам |
в некоторой окрестности точки t, x t у (например, ui+1= |
= и(х+ А х, t) и т. д.). |
Члены нулевого (низшего) порядка представляют при |
этом исходные дифференциальные уравнения, учет в разложении членов более высокого (&+1)-го порядка (k= 0, 1, 2, .) позволяет определить дополни тельные (к уравнениям) члены разложения, а следовательно, и структуру диф ференциального приближения k-то порядка (k-eo дифференциального прибли жения).
88 ИССЛЕДОВАНИЕ ЧИСЛЕННЫХ СХЕМ МЕТОДА КРУПНЫХ ЧАСТИЦ [ГЛ. IV
По своей сути, дифференциальное приближение занимает промежуточное место между исходным дифференциальным уравнением и аппроксимирующей его разностной схемой. Дифференциальное приближение имеет структуру дифференциального уравнения, но его коэффициенты зависят от параметров рассматриваемой схемы. Таким образом, нулевое дифференциальное прибли жение сохраняет информацию об исходных уравнениях, а дополнительные члены разложения более высоких дифференциальных приближений зависят от свойств разностной схемы [223]. Дифференциальная структура приближения облегчает аналитические исследования свойств разностной схемы (особенно вопросов устойчивости и образования диссипативного механизма), делает эти исследования более наглядными.
2. Обсудим вопрос о порядке аппроксимации исходной системы диффе ренциальных уравнений (3.1) разностной схемы метода крупных частиц (3.4) — (3.4'), (3.10)— (3.12), (3.25) — (3.25') [29,217, 219]. Будем считать, что общий порядок аппроксимации системы равен наименьшему порядку аппроксимации отдельных уравнений системы. Следует отметить при этом, что необходимо рассматривать аппроксимацию каждого уравнения на всех трех этапах сразу [386]. Суммарный порядок аппроксимации нельзя принять равным минимальному порядку аппроксимации на этих этапах, рассматриваемых от дельно, так как на каждом этапе аппроксимируются уравнения, отнюдь не совпадающие с исходными дифференциальными уравнениями системы (3.1).
Воспользуемся здесь введенным выше понятием дифференциального при ближения разностного оператора. Для определения порядка аппроксимации следует выписать нулевое дифференциальное приближение исследуемой системы конечно-разностных уравнений с оценкой порядка следующих членов. Если эта выражение совпадает с исходной системой дифференциальных уравнений с точностью до О (As"), то будем говорить, что система конечно-разностных уравнений обладает я-ым порядком точности по s [326].
Рассмотрим случай, например, когда поток газа течет слева-направо и снизу-вверх (рис. 3.2). Другие ориентации потока могут быть получены из вышеуказанной путем элементарных преобразований (поворот, зеркальное отражение) и поэтому не меняют существа дела.
Определим отдельно порядок аппроксимации для внутренних точек потока и для точек, находящихся у поверхности тела.
Выпишем вначале нулевое дифференциальное приближение уравнения неразрывности для внутренней ячейки поля в случае использования для
потока массы ДА!" формул (3.12) первого порядка точности |
|
||||
Р?.У ДХ АУ= Р" i Ах АУ+ АЩ - г/2. j + Ш 1 / - ! / . ■ - ■ i |
АЩ, /+1/2 = |
||||
= P i , Ах Ау + |
Ш |
[pU , (u'U ,■+ u l ,) - |
Pi,, (Hi j + |
й?+1, ,)] + |
|
Ддс Д / |
|
ы + < / ) — Р ? , / К / + » ? , м ) ] = |
|||
- [ P “ |
|||||
= p i ,-AxAy + ^ |
{ |
[ p - | Ax + O (Дх2)] [ 2 5 - g |
Ax + O (Дх2)] - |
||
|
— iP [ 2 “ + T x A x + 0 (Д * 2) ] } " / + |
|
|
||
+ Ц г { [P- |
% ь у + о (V )] [ & - § |
ау + ° |
(A</2)] - |
||
|
— P [2" + Щ Ay + 0 (Af/Z)] |
; = |
|
|
|
= P'/, / Ax Ay + |
1 — 2« A x — 2p ^ Ax + O (Дх2) j" f + |
||||
+ Ц М { - 2ид/ у А у - 2 р ^ А у + 0 (Ду2)}" |
.+ О (Д/2 Дх At). |
АППРОКСИМАЦИЯ УРАВНЕНИЙ |
8 9 |
Отсюда |
|
или |
|
p V H P l i - ^ [ { ^ p W } l i + 0 ( M , Ах, At/)]. |
(4.6) |
Мы видим, что выражение (4.6) совпадает с исходным дифференциальным уравнением неразрывности (3.1) с точностью до О (At, Ах, Ау). Таким образом, разностное уравнение неразрывности при использовании формул (3.12) имеет первый порядок точности. Легко проверить также, что при использовании для А М п формул второго порядка точности (3.11) порядок аппроксимации по х, у будет вторым.
Получим теперь нулевое дифференциальное приближение уравнения для горизонтальной составляющей скорости в случае использования формул (3.12)
Отсюда
Таким образом, и уравнение для и также имеет первый порядок аппроксимации при использовании формул (3.12). Легко показать, что и в случае при менения формул (3.11) порядок аппроксимации будет первым.
Аналогично можно провести рассуждения для уравнения энергии и вто рого уравнения импульса.
Итак, общий порядок аппроксимации метода крупных частиц для внут ренних точек поля — первый, хотя при использовании формул (3.11) урав нение неразрывности имеет второй порядок точности по пространственным координатам.
Следует отметить, что при расщеплении исходной системы дифференци альных уравнений на части и аппроксимации каждой из этих частей схемой /1-го порядка точности суммарный порядок аппроксимации первоначальной системы может быть отличен от п. Его величина зависит, вообще говоря, от способа комплексации расщепленных частей [379].
Перейдем к рассмотрению г р а н и ч н ы х ячеек. Определим порядок аппроксимации в ячейке а, лежащей на поверхности тела (рис. 3.2). Если на теле ставятся условия непротекания, т. е.
то
;-+1/2
90 ИССЛЕДОВАНИЕ ЧИСЛЕННЫХ СХЕМ МЕТОДА КРУПНЫХ ЧАСТИЦ [ГЛ. IV
И вновь получим, учитывая условие др/ду=0, первый порядок аппроксимации
для уравнения |
неразрывности |
|
|
|
Р/,+/ = P l i ~ At [{div PW}<. i + ° (A t>Ax>ЬУ)]- |
(4 -6') |
|||
Рассмотрим теперь нулевое дифференциальное приближение уравнения |
||||
для и с оценкой порядков следующих членов: |
|
|||
ui j P ? j AxAy — |
|
|
|
|
|
й - А , ] ^ * 4 » + ' | ( “ - |
# |
л * ) ( р — з г 4 - ) х |
|
Х_Н Г ^}" .+ |
° ( д*2. b y2) = [«?.уР?./— |
А^] AxAy— A xA yA tx |
|
|
|
х [{ 2Р“ 1 г + |
+ Р“Ц } " j + 0 (Д<- А*> А*)] • |
(4.9) |
|
Здесь было использовано следствие из граничных условий |
|
Покажем, что выражения, стоящие справа в фигурных скобках в (4.9)„ являются div (puW )+0(A /, Ах, А#).
Действительно,
dlv(puW) = 2 p u % + u * £ + u o & + pu2jj + p o ^ .
В силу граничных условий (4.8)
[ ° + f т+ °»4 л ]
Тогда из (4.7)
« Е М ? = «?. f i t i - A t [ { % - + div (Р“ W0}“ у+ О (At, Ах, At/)], (4.7')
т. е. уравнение для и и вблизи поверхности тела имеет первый порядок аппрок симации. Аналогичные выкладки, сделанные для уравнения энергии и вто рого уравнения импульса, приводят к тому же результату.
Здесь рассмотрен вопрос о порядке аппроксимации на границе с обтека емым телом в случае условий непротекания (4.8). Можно показать (выкладки здесь не будем приводить ввиду их полной аналогии), что и в случае условий прилипания разностная схема вблизи поверхности тела имеет также первый порядок точности.
Таким образом, порядок аппроксимации используемой разностной схемы метода крупных частиц равен единице во всей области течения: как во внут ренних точках, так и на границе с твердым телом. Эти результаты, как отме чалось, получены из рассмотрения нулевого дифференциального приближения с оценкой порядков следующих членов.
Сделаем здесь одно замечание. При расчете сложных задач газовой ди намики (например, для нелинейных уравнений в частных производных, при наличии сильных разрывов) использование схем формально более высокого порядка точности кажется нецелесообразным. Дело в том, что хотя в качестве одной из априорных характеристик схемы и используется погрешность