Специальные разностные схемы для сингулярно возмущенных краевых задач тема автореферата и диссертации по математике, 01.01.07 ВАК РФ

ГАЕВОИ Виктор Павлович

V —---—---

СПЕЦИАЛЬНЫЕ РАЗНОСТНЫЕ СХЕМЫ ДЛЯ СИНГУЛЯРНО ЕиЗМУЩЕННЫХ КРАЕВЫХ ЗАДАЧ

01.01.07 -' вычислительная математика

АВТОРЕФЕРАТ

диссертации на соискание учёной степени кандидата физико-математических наук

Новосибирск - ■ 1995

Работа выполнена в Институте катализа им. Г.К.Борескова

Научный руководитель - доктор физико-математических

наук, профессор Зеленяк Т.И.

Официальные оппоненты - доктор фйзико-математических

Ведущая организация - Институт гидродинамики

им. М.А.Лаврентьева

в 15 часов на заседании диссертационного совета К 063.98.04 в Новосибирском государственном университете по адресу: 630090, Новосибирск -90, ул. Пирогова, 2, НГУ.

С диссертацией можно ознакомиться в библиотеке Новосибирского' государственного университета.

наук, доцент Мацокин А.М

доктор физико-математических наук, доцент Лаврентьев М.М.-мл.

Защита состоится

Учёный секретарь диссертационного совета

В.В.Шелухин

Актуальность темы. Краевые задачи для дифференциальных уравнений с малым параметром при старшей производной ( сингулярно возмущенные краевые задачи ) встречаются в различных областях науки и техники, например, уравнения движения вязкой кидкости, океанических .течений, переноса тепла и вещества, химико-технологических процессов. Главной особенностью сингулярно возмущенных краевых задач, затрудняющей их численное решение, является наличие узких зон с резким изменением решения и большими значениями производных, получивших название пограничных и внутренних переходных слоев, ширина которых пропорциональна значению малого параметра в той или иной положительной степени. Для нахождения- с приемлемой, точностью приближенного решения сингулярно возмущенных краевых задач традиционные конечно-разностные и вариационно- разностные методы требуют выбора шага разностной сетки, значительно меньшего ширины пограничного слоя. С уменьшением значения малого параметра при старшей призводной это ведет к чрезмерно мелкому шагу разностной сетки. В практических задачах параметр при старшей производной может принимать значения в довольно широких пределах. Поэтому естественно возникает потребность в разработке численных методов, обладающих равномерной сходимостью относительно значения параметра при старшей производной и позволяющих находить приближенные решения сингулярно возмущенных краевых задач с требуемой точностью на разумных сетках.

К настоящему времени в развитии равномерных численных методов сложились два основных направления.Одно из них связано с применением известных конечно-разностных схем на специально сконструированных неравномерных сетках, сгущающихся по определенному закону в области пограничного или переходного слоя. Другим направлением является построение специальных разностных схем, имеющих равномерный по малому параметру порядок точности на равномерных или произвольных неравномерных разностных сетках. Такие разностные схемы называются схемами

равномерного порядка точности или просто равномерными разностными схемами. До работы [4] в литературе были известны разностные схемы равномерного первого порядка точности на равномерных разностных сетках.

Целью данной работы является разработка методов построения специальных разностных схем, учитывающих особенности поведения решения исходных краевых задач, и построение с помощью этих методов равномерных разностных-схем для сингулярно возмущенных краевых задач, имеющих равномерный второй порядок точности на равномерных или произвольных неравномерных разностных сетках.

Научная новизна.

Для обыкновенного дифференциального уравнения второго порядка с граничными условиями произвольного вида предложен метод построения специальных разностных схем, развивающий метод интегрального товдества Марчука Г.И. и метод точной и усеченных разностных схем Самарского A.A. и позволяющий строить разностные схемы, учитывающие особенности поведения точного решения исходной краевой задачи

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

первого рода и равномерной разностной сетки.

Апробация работы. Основные результаты работы докладывались на Втором Советско - Французском семинаре по математическому моделированию каталитических процессов и реакторов (Новосибирск, 1976), Всесоюзном семинаре по методу сплайн-функций (Новосибирск, 1978), Всесибирской школе по вычислительной математике (пос. Шушенское Красноярского края, 1979), на семинаре по методам вычислительной и прикладной математики ВЦ СО АН СССР (руководитель -'академик Г.И.Марчук), семинаре по численным методам механики сплошной среды ИТПМ СО АН СССР (руководитель - академик Н.Н.Яненко), семинаре по качественной теории дифференциальных уравнений ИМ СО РАН и НГУ (руководитель - д.ф.- м.н., проф. Т.И.Зеленяк)

Публикации. По теме диссертации опубликовано 5 работ.

Объём работы. Диссертация состоит из введения и трёх глав, содержит Ю5 страниц машинописного текста. Список литературы включает 75 наименований.

СОДЕРЖАНИЕ ДИССЕРТАЦИИ

Во введении к диссертации обоснована актуальность теш, сформулирована цель работы,приведён обзор основных результатов по обсуздаемой проблеме, приведено краткое содержание диссертации с указанием главных результатов.

В параграфе 1.1 рассматривается следующая краевая задача.

Ш - (р(х)их)х + г(х)их- ц(х)и = - Г(х) , х € (0,1) , р(х) > 0 , 10и - - а0их(О) + рои(0) = 70 > (1)

1.,и - а,их(1) +-^^(1) = 71 \а.\ + 1(^1 > 0 , 1 = 0,1.

Предполагается, что р(х), г(х), q(x), г(х) е 0°[0,1] •

Здесь и'далее используются следующие обозначения. Оп[о,1] - класс функций ограниченных и кусочно-непрерывных вместе с производными до п-го порядка включительно, имеющих конечное число точек разрыва первого рода, с®[0,1] - класс функций из сд_1[0,1] , имеющих ограниченные

¿Xй

производные п-го порядка , (иЦ^о = ||и||,, + Бир

• П • XI—1

с* [0,1] - класс функций из с°[о,1] , имеющих ограниченные кусочно-непрерывные производные п-го порядка с конечным числом точек разрыва первого рода.

Под решением задачи (1) понимается непрерывная кусочно-дифференцируемая функция, обладавшая непрерывным кусочно-дифференцируемым потоком ^(х) = р(х)их(х) .

В дальнейшем предполагается, что на отрезке [0,1] задана произвольная разностная сетка: о = х0<х1<... <хк = 1 , = , 1=1,2,...N. Ь. = тах

Вводится вспомогательный дифференциальный оператор

Ь° = — р°(х)— + г°(х) — - ч°(х) , йх с1х 4х

коэффициенты которого выбираются таким образом, чтобы при

х € (х^.х^), 1=1,2,...N оператор Ь° был определен на

решении уравнения (1) и обратил в явном виде на отрезках

[х1_1,х1] с граничными условиями первого рода.

Обозначим через , [х1_1,х;.] решения

следующих краевых задач

о о

ЬУ01(х) = ЫГП(Х) = 0 , х€ (ХЬ1,Х1) ,

= уи<х1> = У01(х1).= = 0 • .

В работе показано, что для решения задачи (1) справедливы равенства

"«о^о У01(0) + и1 у11(0) + Ф + Рои = Ъ •

Pi [Ui-1 vol<xi) + Ui + Fi] =

= Pi [Ui Voi+1 <xi> + Ui+1 vii+i<xi> + pi] •

i=1,2,...N-1 , <MUN-1 V0N<1> + UN V1N<1> + V + = Ъ '

где Uj = UCxj), pT = p (х±-0), pt = p (xt+0) x.

(3)

Fi = J ^(*!.?)((Ir-b°)in-f)dg ,

xi-1 Л .

Fi-1 = J Cx(xi-1'Ü((L_L )U+r)d? ' xi-1

G^x.C), x,? с [xl_1,xi] - функция Грина задач (2).

Выбирая специальным образом коэффициенты дифференциального оператора L0 и строя различные приближенные формулы для вычисления функционалов Р^, pt, можно получать с помощью равенств (3) специальные разностные схемы, учитывающие особенности поведения решения задачи (1)

Равенства (3), построенные на равномерной разностной сетке для уравнения самосопряжённого вида ( г(х)- о ) с граничными условиями первого рода, при р°(х) - р(х), г°(х) - 0, q°(x) ■ о совпадают с интегральным товдеством Марчука Г.И., а при р°(х) - р(х), г°(х) - о, q°(x) - q(x) - с точной разностной схемой Самарского A.A.

В § 1.2 рассматривается краевая задача следующего вида.

L°U(x) + P(L-L°)SU(x) = - Pf(x), х € (х1_1,xi), i = 1,2,...N ,

p°(xi-0)U(xi-0) = p0(x1+O)U(xi+O), 1=1,2,...N-1 , (4)

l0U(x) = 70 , l.,U(x) = 7, ,

где S - оператор интерполирования функции и(х) е G[0,1] в узлах разностной сетки su(xi) = и^), i= 0,1,2,...N, P - оператор локального приближения функций на интервалах разностной сетки.

Предполагается, что выражения (L-L°)SU(x) определены и ограничены внутри интервалов разностной сетки (xi_1,sl). Под решением задачи (4) понимается непрерывная кусочно-дифференцируемая функция, обладавшая непрерывным кусочно-дифференцируемым потоком g°(x) = p°(x)Ux(x) . Интегральные тождества (з), построенные для задачи (4), дают для неб точную разностную схему.

ЛЕММ А I. Пусть дифференциальный оператор L0 однозначно обратим на отрезках [х^.х.^, i= 1,2,...N, с граничными условиями первого рода, тогда из однозначной разрешимости задач (4) следует однозначная разрешимость точной для неб разностной схемы , и наоборот, из однозначной разрешимости точной для задачи (4) разностной схемы следует однозначная разрешимость задач (4).

В § 1.3. изложен метод построения специальных усечённых разностных схем, развивающий известный метод усечённых разностных схем Самарского A.A.

Пусть выбраны некоторые последовательности функций рк(х), гк(х), qk(x), Гк(х), к = 0,1,...п-1, аппроксимирующих функции р(х), г(х),q(x), £(х) соответственно, и введены в рассмотрение вспомогательные дифференциалльные операторы

k d 1г d k d lr L = — p (x)— + г (x)--q (x) , k = 0,1,...11-1 .

dx dx dx

Предполагается, что внутри интервалов разностной сетки дифференциальные операторы Lk, к = 0,1,...п-1 определены на одних и тех не функциях, и дифференциальный оператор L0 однозначно обратим на отрезках [х^ .х^ , i= 1,2,...N, с граничными условиями первого рода .

Обозначим через vj^U), vi^íx), w!?(x) , к = 0,1,...п , i = 1,2,...N функции, являющиеся последовательными решениями следующих краевых зададч

V°1(x) - О, V^íx) - О, WО (х) - О ,

0 lr+1 k 0 V

L V^t (x) = - (L - L ^..(x) , x e (xi_1,xi) , i = 0,1 ,

vw<xi-i> = = 1 ' = <xi-1 > = o •

L°W^+1(x) = - (Lk- (x)- fk(x) , x € (х±_1,х±) ,

wf1 (x.) = "¿^(x^) = O, k=0,1, —n-1 , i=1,2,...N .

Для функций u(x) ё с[о,1] определим последовательность операторов интерполирования sk, к = 1,2, п

S^x) = u(xi_1)Yki(x) + U(x£) V*.(x) + (x),

x € [xi_1,xi] , 1 = 1,2,...N , к = 0,1,...n .

Приближенное решение задачи (1) определим как решение следующей краевой задачи

L°Un(x) + (Ln-1-L°)Sn-1Un(x) = - fn_1(x),

х € (xi_1,xl), i = 1,2,...N , (5)

p°(xi-0)U^(xi-0) = p°(xi+0)U^(xi+0), i=1,2,...N-1 ,

lQUn(x) = T0 • l-,Un(x) = . Точная разностная схема для задачи (5) имеет следующий вид

-VUS А? + и? в? + F?> + Pouíl= То •

p¡ [о»., h + U? Bl + íj = (6)

где

= Pi [и? А°+1 + В°+1 + F°+1 ] , 1=1.2....N-1. а1(UN-1AN + UN % + V + Мн = Ъ '

Ai= VSix<Xi> ' Bi= ' Pi= •

4 = YSix<Xi-1> ' Bi= VUx<xi-1> ' Fi= «ix^i-l^

Используя понятие ранга усеченной разностной схемы, введённого Самарским A.A., разностную схему (6) будем называть специальной усечённой разностной схеммой ранга п , если выражения для ее коэффициентов и правых частей получены путем п - кратного обращения на интервалах разностной сетки вспомогательного дифференциального оператора L0 .

В второй главе рассматривается сингулярно возмущённая краевая задача для уравнения самосопряжённого вида.

LU - е2(р(х)их)х - q(х)и = - г(х) ,

х € <0,1) , р(х), q(x) »1, 1 i£> О,

10U - - aQp(0)ux(0) + ß0u(0) = т0 f (7)

ци - сцр(1 )их(1) + р1и(1) = 71

а.. * о , * о , + = 1, з = 1,2 .

Предполагается, что е, Oq, ol, . ß0, ß1 могут принимать сколь угодно малые значения

Через М, MQ, М1, М2,... в дальнейшем обозначаюся константы не зависящие от a, aQ, о,, ß0, ß1

В параграфе 2.1 доказаны следующие утвервдения.

Л Е М М А 2. Пусть функции р(х), q(x), f(x) е Q1[0,1], 1 *р(х), q(x) £ R ; о < g1< g2 <...< gjj < 1 - точки разрывов

функций р(х), q(x), f(x), f¡0 = о, ?N+1 = 1- Тогда для решения задачи (7) и его первой производной справедливы оценки

IU(x)'l s М1 , IUx(x)l s M[1+e_1 (e"*(x_?i-1 )/Re + e~(5i~x )/R£)],

2 € Ii), i = 1.2....N+1 ,

Обозначим через U(x) решение краевой задачи Lu- s2(p(x)Ux)x - q(x)U = - г(х) ,

10U - - 0qP(0)Ux(0) + ß0U(O) = 70 f (8)

l.,u - о,р(1 )UX(1) + ß^uc-t) = 71

Л E M M А 3. Пусть p(x),q(x),f(x) € Q1[0,1], p(x),q(x) г 1 , p(x), q(x), f(x) e Q1[0,1], p(z),q(x) i б > О ,

pío) = pío) , p(i) = pd) , p(x)/p(x) e c*[o,i ] .

Тогда для разности решений краевых задач (7) и (8) справедлива оценка

||U(x)-U(x)||c i М (Up-pH 0 + ||q-q|| 0 + ||f-I|| 0) .

С С G

В § 2.2. строится приближённая краевая задача и точная для неё специальная усечённая разностная схема второго ранга.

Обозначим через а°(х), ь°(х) кусочно-постоянные функции

а°(х)= а?, Ь°(х)= b9, х е (х1_1,х1), а?>0, i=1,2,...N.

Зададим вспомогательный дифференциальный оператор

L0 = б2 — р(х)— - q°(x) , q°(x) = а°(х)/ р(х) . dx • dx

В данном случае функции voi(x),v1i(x), являющиеся решениями

краевых задач (2), находятся в явном Биде

Приближённое р е шение за д ачи (7) строиться как решение следующей краевой задачи

е2(р(х)и£)х - q°(x)Uh - q1 (x)SUh = - f1 (x) ,

I0uh = To. ¡У = T1f (9)

q1(x) = q (x) - q°(x), i1 (x) =f(x) - q1(х)ф(х) SUh(x) = Uh(xi_1)Voi(x) + Uh(xi) Vu(x) ,

.ф(х)=(ь5/а^)(1-У01(х)-Уп(х)), x € [x^.x.], i=1,2,...N .

Для задачи (9) строится точная разностная схема, являющаяся специальной усечённой разностной схемой второго ранга для задачи (8), и доказываются следующие утверждения

ЛЕММА .4. Пусть функции р(х), q(x), f(x) € Q°[0,1],

р(х) > 0 ,q(x) > О , а°(х) >0, тогда краевая задача (9) и точная для неб разностная схема однозначно разрешемы.

ТЕОРЕМА I. Пусть выполнены предположения ЛЕММЫ 3 , разностная сетка и кусочно-постоянная функция а°(х) г. ö >0 выбраны таким образом, что

I(р(x)q(x) - а°(х))/а°(х)| * d < 1 .

Тогда для разности решений краевых задач (7) и (9) справедлива оценка

||U(x)-Uh(x)||c * M^IIp-pII Q+ Ilq-q|| 0+||f-î|| Q) +

С С С

+ M2||pq - a0H 0(||й - а°|| 0+||pf - Ь°И 0).

С С С

Далее функции Р(х),q(x),i(х) выбираются таким образом, что они удовлетворяют требованиям ТЕОРЕМЫ I, достаточно точно

аппроксимируют функции Р(х), д(х), 1(х) , и позволяют точно вычислять интегралы, входящие в выражения для коэффициентов и правых частей точной для задачи (9) разностной схемы. В этом случае точная разностная схема для задачи (9) имеет следующий вид

- <е2р0 + аов?)и£+ VIй? = - е\ - аор? •

А1и1-1- <в?+1+ в1>и1 + А1+1и±+1 = - 4 - р1+1- <10>

1 = 1,2,...N-1 ,

¥1иЙ-г + аХ)ик = " е2Т1 - -

где

А1 = Е2\1/БЬ2.± , = е2А,1(С^1 -(-1 )кС1ф1 (г^), к = 0,1, Р1 = ЬТф2(г.)-Ь°С1(ф1(г.)-2ф2(21))],

Р1 = ь7ф1(2.)-ь?с1(ф1(г1)-2ф2(г1))],

ф^г) = ОНлг - 1/г , ф2(г) = 1/е - 1/8112 ,

а1 = ч(х1± °>Р<Х:Г Ъ1 = °)Р(21±

а9=(а^_1+аТ)/2, ь9=(Ь^_1+ьТ)/2 , С1 = (аТ - а^ )/4а9 ,

Доказывается, что разностная схема (ю) монотонна и для еб решения справедливо следующее утверждение.

ТЕОРЕМА 2. Пусть функции р(х),а(х),Г(х) е <32[0,1], р(х),ч(х) г 1 ,и разностная сетка выбрана таим образом, что

p(x), q(x), f(x) € ^(x^.x^ , i = 1,2,...N ,

тогда для разности решений задачи (7) и разностной схемы (Ю) справедлива оценка

luup - U?>| s М max h| ,

В параграфе 2.3. приближённое решение краевой задачи (7) строится как точное решение следующей краевой задачи

е2(р(х)и£)х - q°(x)Uh = - ъ°(х)/ р(х),

(11)

^ =т0. Т1 ,

где функции p(x),q0(x)fb°(x) те же, что и в задаче (9). Точная разностная схема для задачи (11) имеет вид (Ю), где ,

А. = s2Ä../Shz., df = e2Ä..cthz., 1 Iii Ii (12)

pi = ^i'^iChZj-l)/ Shzlf k=q,l, i=1,2,.. .N.

Решение задачи (11) внутри интервалов разностной сетки выражается в явном виде через решение разностной схемы

Uh(x) = U^^^xJ+U^V^xi+ib^/a^d-Y^xJ+V^x)),

х € [Xj^.Xj] , i = 1.2....N , (13)

V°(x) = KT1Sh[^ihT1((h|-(x-x._1)2)/pi+(xi-x)2/pi_1)/2] ,

v9(x) = K~1Sh[Xih^((h|-(x-xi)2)/pl_1 + (x -x^, )2/p±)/2l ,

K. = Sh^h^l/p.^+l/p^) .

ТЕОРЕМА 3. Пусть функции p(x), q(x), i(x) € Q1[0,1], p(x),q(x) г 1 ,и разностная сетка выбрана таким образом, что

р(х), q(x), Г(х) € CU^.x^ , i = 1.2....N , тогда для погрешености приближённого решения задачи (9),

построенного с помощью разностной схемы (Ю), (12) . и равенств (13) справедлива оценка

|U(x) - Uh(x)| ä М шах hj ,

ТЕОРЕМА 4. Пусть функции p(x),q(x),г(х) € С*[0,1], p(x),q(x) i 1 .тогда для разности решений первой краевой задачи для уравнения (7) и разностной схемы (ю),(12), построенной на равномерной сетке с шагом h , справедлива оценка

|u(Xi) - * М h2 .

В третьей главе рассматривается сингулярно возмущённая краевая задача для уравнение несамосопряженного вида.

£Ont a(x)Ux- b(x)U = Г(х) , х € (0,1) е > О, A i а(х) ¿1, В i b(x) i О , 10U - - aQUx(0) + ßQU(1) = т0 f (14)

i.,u - a.,ux(i) + ß^d) = aJf ß^.i о, a..+ ß.=i, J = 1,2 , Ь(х)+р1г б > 0

Предполагается, что коэффициенты е, aQ, ct., , ßQ, ß1 могут принимать произвольно малые значения. Так как в граничных условиях допускается одновременное обращение в ноль значений ß0 и ß1, то требование b(x) + ß^ ö> о является существенным.

В § 3.1 доказана справедливость следующих утверждений.

Л Е М М А 5. Пусть а(х), b(x), f(x) € Q°[0,1], тогда для решения краевой задачи (14) справедливы оценки

IU(x)I s М (ITqI + l7ll + Suplf(x)l), IU (x)l sM (1 + (8+a0)~1exp(-x/s)).

Для различных частных случаев граничных условий оценки первой производной решения задачи (14) были получены в работах других авторов, но при этом предполагалось, что коэффициенты и правая часть уравнения ограничены в С., [0,1]. -

Л Е М М А 6. Пусть а(х) € С*[0,1], Ь(х), Г(х) € С*[0,1], тогда решение задачи (14) представимо в следующем виде и(х) = и0(х) + и.,(х)г(х) ,

где

г(х)= ехр(- е-1/\(з)с1з) , ||и3(х)1|с0[0>13 л М , 3=1,2.

В 3.2. для задачи (14) строится на произизволльной неравномерной разностной сетке специальная усечённая разностная схема первого ранга. Значения а?, ъ9, 1 = 1,2,...к , задются одним из следующих двух способов:

а9 = а(х1-Ь1/2) , ъ9 = Ь(х±-111/2) , = £ (х^^/2) ,

а9 = (а(х1_1)+а(х±))/2, ъ9 = (Ъ(х1'_1 )+Ь(х1))/2,

Х° = (£(х1и)+Г(х1))/2, 1 = 1,2,...И .

Приближённое решение задачи (14) находится как решение следушей краевой задачи

е и^.(х) +-а°(х)1^(х)г- Ь°(х)и11(х) = Г°(х) , (15)

10^70.1^ = 7,

Точная разностная схема для краевой задачи (15) имеет следующий вид.

- е<Р0+а0в0)и£ + шосо<= -£Т0+ »

е[А1и1-Г (в?+в1)и1+ с1и1+1^ = <ф°+1 (16)

1 = 1 ,2,...11-1 ,

ewti- Е^1+а1вК = - sv -

где

2z,e~zi~yi 2z.e~zi+yi А- ■*-__П. = _

1 -e-2zi) ' ^О-е"22!) '

n Zj-y^tyy^e-22! « z1+y1+(2i-y.)e-2zi

hi(1-e-2zi) ' i"1 = hi(1-e-2zi) '

= hi[®(zi+yi) - е^Г^Ф^-у^/О-е-22!) , Qi-1 = 1Ч[Ф(21-У1) - е_21+у1Ф(z^Jl/d-e"22!) , Ф(г) =(1-e~2)/z , yt= (i^ , Zj= V^ .

ц^ a?/2e . Xi = (ц2 + b°/e)1/2

Решение задачи (15) внутри интервалов разностной сетки выражается в явном виде через решение разностной схемы (16).

Л Е М М А 7. Пусть Ъ° i о, 1а9|+ь?> о для всех i=1,2,...N,

(V. i О, i О, O-j+Pj^ О , 3=0,1, max Ъ°+ PQ+ pi > О,

тогда разностная схема (16) монотонна на произвольной разностной сетке и для неб справедлив принцип максимума.

ТЕОРЕМА 5. Пусть а(х), Ь(х), Г (х) е Q1 [0,1 ], разностная сетка выбрана таким образом, что внутри её интервалов функции а(х), ь(х), f(x) непрерывны, тогда для разности решений задач (14),(15) справедливы оценки

IU(x) - Uh(x)l sMh ,

.Ux(x)-^(x)U И h(1 + (8+a0)_1e_x/28).

В §3.3 доказано следующее утверждение.

Т Е О Р Е М А 6. Пусть функции а(х) ,Ь(х) .КхКС^ОЛ ], тогда для погрешности приближённого решения первой краевой задачи для уравнения (14), полученного с помощью разностной схемы (16) на равномерной сетке , справедлива оценка

тах|и(х1)-и^| * М И2.

Список работ автора по теме диссертации.

1. Гаевой В.П. Полуаналитический метод решения обыкновенных дифференциальных уравнений и уравнений в частных производных параболического типа // Управляемые системы.-Новосибирск, 1972.-Вып. 10.- с.62-70.

2. Гаевой В.П. Сходимость полуаналитического метода для квазилинейного уравнения парабоолического типа //

. Управляемые системы.- Новосибирск, 1972.-Вып.10.- с.70-78.

3. Гаевой В.П. Схемы высокого порядка точности для уравнения парабоолического типа // Тр. Второго Совет с ко- Француз -ского семинара по математическому моделированию каталитических процессов и реакторов.- Новосибирск, 1976.-с.260-264.

4. Гаевой В.П. Об одном методе построения разностных уравнений для двухточечных краевых задач // Вычисл. системы.- Новосибирск, 1978.- Вып. 75 : Метод сплайн-функций.- с. 96-109.

5. Гаевой В.П. Численные методы решения краевых задач с большими градиентами // Математическое моделирование химических реакторов./ Под ред. Г.И.Марчука. - Новосибирск, 1984.- с.144-162.