«Асимптотическое разложение решения одного параболического уравнений второго рода»

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

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

Задача, поставленная в выпускной квалификационной работе, связана с задачей о диффузии около сферической частицы при обтекании ее потоком жидкости и при наличии объемной химической реакции. Процессы массо- и теплообмена в дисперсионных средах является предметом многочисленных экспериментальных и теоретических исследований в связи с их большим значением для химической технологии и других областей техники.

Рассмотрим краевые задачи для уравнений эллиптического типа.

Задача о стационарном распределении температуры u(x,y,z) внутри тела Т, ограниченной поверхностью S, формулируется следующим образом:

Найти функцию u(x,y,z), удовлетворяющую внутри Т уравнению

 u=-f(x, y, z) и граничному условию, которое может быть взято в одном из следующих видов:

1.u=f1 на S(первая краевая задача),

2. на S (вторая краевая задача)

3. на S (третья краевая задача)

где f1,,f2,,f3 ,h заданные функции, - производная по внешней нормали к поверхности S

.

Первую краевую задачу для уравнения Лапласа  u=0 часто называют задачей Дирихле, а вторую – задачей Неймана.

Если нужно найти решение в области Т0, внутренней (или внешней) по отношению к поверхности S, то соответствующую задачу называют внутренней (или внешней краевой задачей).

Наилучшие результаты по разрешению задачи Дирихле в ограниченных областях D были доказаны Шаудером и частично Каччопполи. Они формулируются в терминах пространства Гельдера Cl+ (D). Элементы С (D)являются функции u(x), непрерывные в D в смысле Гельдера с

показателем ((,1)), т.е. u(x) непрерывна в D и для нее конечна построенная Гельдера.

Элементами Сl+ (D), l  1 , является непрерывные в D функции, имеющие всевозможные производные до порядка l, причем производные порядка l суть элементы С (D). Норма в С+l (D) определяется равенством:

В ней выражение (k) означает суммирование по всем производным порядка k. Пространства Cl+(D) , l  1 являются банаховыми. Таковыми же являются и пространства Cl(D), j=0,1… , состоящие из непрерывных в D функций, имеющих непрерывные в D производные до порядка l. Норма в Cl (D) определяется следующим образом:

Вместо С0(D) принято писать просто C(D). Будем говорить, что граница S в области D  Rn принадлежит классу Cl+,l1,  (0,1), если существует число  > 0, такое, что пересечение S с шаром В радиуса  с центром в произвольной точке x0  S есть связанная поверхность. Уравнение, которое в местной декартовой системе координат (y1, …,yn) с началом в точке x0 имеет вид

yn = (y1,…,yn-1), причем (y1…yn-1) есть функция класса Сl+ в замкнутой области В(x0,), которая является ортогональной проекцией пересечения на плоскость yn = 0.

Термин «местная декартовая система координат с началом в точке

x0  S» означает, что ось yn направлена по нормали к S в точке x0, а остальные, ортогональные друг другу yi лежит в гиперплоскости, касательной к S в точке x0 . Для функции , заданной на поверхности S класса Сl+ , принадлежность к Сk+(S), k+ l+ , будет означать, что она как Функция местных координат (y1,…yn-1) принадлежит Сk+(B(x0,)) для всех x0  S.

c0m:=0; c1m:=1;

m1:

X: =10; h: =-0.0005; b: =0.0005;

C [0] :=( c0m+c1m)/2;

C [1]:=m*f0(c [0]);

C [2]:=m*p(c [0])*c [1]/2;

C [3]:=m*(p(c [0])*c [2]-pp(c [0])*sqr(c [1])/2)/3;

C [4]:=m*(p(c [0])*c [3]-pp(c [0])*c [1]*c [2] +ppp(c [0])*exp (3*ln(c [1]))/6)/4-1*2*c [1];

c[5]:=m*(p(c[0])*c[4]+pp(c[0])/2*(c[2]*c[2]+2*c[1]*c[3])+

Ppp(c [0])*c [1]*c [1]*c [2]/2+pppp(c [0])*exp (4*ln(c [1])/24))/5-2*3*c [2];

Y: =0;

For I: =0 to 5 do

Y: =y + c [I]*exp (-I*ln(x));

y1:=y; z: =0;

For i: =0 to 5 do

z:=z-I*c[I]*exp(-(i+1)*ln(x));

Repeat

k1:=h*f(x, y, z);

p1:=h*v(x,y,z);

k2:=h*f(x+h/2, y+k1/2, z+p1/2);

p2:=h*v(x+h/2, y+k1/2, z+p1/2);

k3:=h*f(x+h/2, y+k2/2, z+p2/2);

p3:=h*v(x+h/2, y+k2/2, z+p2/2);

k4:=h*f(x+h, y+k3,z+p3);

p4:=h*v(x+h,y+k3,z+p3);

d:=(k1+2*k2+2*k3+k4)/6;

q:=(p1+2*p2+2*p3+p4)/6;

y:=y+d;

z:=z+q;

x:=x+h;

until xy :=(y+1)/2;

If (y<1) then c0m:=c[0]

Else c1m:=c[0];

if abs(y-1)>0.00001 then goto m1;

writeln ('y[0]=',y:5:10);

writeln ('z[0]=',z:5:10);

writeln ('c[0]=',c[0]:5:10);