Method of calculating the dimensions of greenhouse-type single slope watermaker by taking into account the accumulation of solar energy parnik tipidagi bir nishabli suv chuchutgichi o
Download 2.25 Mb.
|
Документ
- Bu sahifa navigatsiya:
- Результаты и их обсуждение
Объект и методы исследованияПрименим видоизмененный метод квазистационарного приближения Лейбензона. Сущность метода заключается в том, что подвижная граница «замораживается», т.е. полагается l t( ) l s( ) const и решается обычная задача сопряжения с вертикальной границей раздела x l s( ). Затем, подставляя вместо l s( ) функцию l t( ) и пользуясь условием на границе x l t( ), получают уравнение для определения x l t( ). Итак, пологая l t( ) l s( ) 0 рассмотрим следующую двухслойную задачу на полупрямой x 0, решение которой зависит от l s( ): ut( )s 2 ( )x 2ux2( )s , ( , )x t t0 \{ :x x l s( ) (11) a u( )s x0 ( )t ; u( )s t t 0 ( )x ; u( )s x 0, (12) u( )sx l s ( ) 0, h1 u( )s x l s ( ) 0 h2 u( )s x l s ( ) 0 , (13) x x где и a определяется (7), (8) с учетом l t( ) l s( ); , и f обладает приведенными выше свойствами (9), символ означает, что f x( ) x a f x a 0 f x a 0. В работе [2] с помощью интегрального преобразования Лапласа-Карсона построена функция Грина двухслойной первой краевой задачи на полупрямой x 0 для уравнения теплопроводности с кусочно-постоянным коэффициентом в случае l s( ) 1. В случае l s( ) 1, l s( ) 0 она имеет вид: g x t1( , , , , ( )), l s 0 x, l s( ), g2(x t, , , , ( )), l s 0 l s( ) x , g x ti ( , , , , l s( )) g x t3( , , , , ( )), l s 0 x, l s( ) , g4(x t, , , , ( )), l s l s( ) x, , где g x ti ( , , , , ( )),l s i 1,2,3,4 являются суммами сходящихся рядов. Отметим, что из-за громоздкости выражения для g x ti ( , , , , l s( )) эти выражения здесь не приведены. Пользуясь интегральным представлением [3] решение задачи (11), (12), (13) нетрудно получить в виде решения следующего интегрального уравнения: u( )s (x t, ) A u( ) ( )s s B x t l s( , ; ( )). (13) Здесь 1 t l s( ) A u( )s ( )s d g x t( , , , , l s( )) (f u( )s 0 )d t0 0 t 2 g x t( , , , , l s( )) ( )d B x t l s( , ; ( )) 0 g x t( , , , ,t l s0 ( )) ( ) d a1 t0 0 Принимаем за приближенное решение задачи (3)-(6) функции u x t( , ) и l t( ) определяемые как решение следующей системы функциональных уравнений : u x t( , ) Au B x t l t( , ; ( )), (14) 0 [Au B x t l t( , ; ( ))] x l t ( ). Здесь Au B x t l t( , ; ( )) [A u( )s ( )s B x t l s( , ; ( ))] l s( )l t( ). Результаты численного решения задач типа Стефана незначительно отличается от точного решения, полученного в квазистационарном приближении Лейбензона для достаточно малых значений времени [1]. Однако вопрос об оценке погрешности метода остается открытым. Систему (14) запишем в виде матричного операторного уравнения: u P u l2( , ) 0 P l P u l2( , ) 0 , (15) где P u l1( , ) u Au B x t l t( , ; ( )), P u l2( , ) 0 [Au B x t l t( , ; ( ))] x l t ( ) . uСчитаем, что областями действия и значений оператора Pl является соответственно пространства C и C , где C C1 C2; C1 u x t( , ):uC[0,l] [ t T0, ], C2 l t( ): ( )l t m t t( )( t0 )1/2, m t( )C[t T0, ],1/ 2 1, b max ( )l t l T( ). t0t T П од знаком понимается прямое произведение . Определение нормы в пространствах C*, C1и C2 : u C1 maxx t, u x t( , ) ; l t( ) C2 max (t t0)1/2l t( ) ; u u l * u C1 l C2 ; P l P u l1( , ) C P u l2 ( , ) C . C Применим к уравнению (15) метод Нъютона-Канторовича [3]. В начале рассмотрим случай 0 x l t( ). Положим l t0( ) 2a t1( t0)1/2, 1/ 2 1, u x t0( , ) B x t l t( , ; 0( )). Тогда первое приближение процесса Нъютона находится решением следующего матричного уравнения относительно поправок u x t( , ) u x t1( , )u x t0( , ) и l l t1( )l t0( ) : P1u(u l0, )0 u P1l(u l0, )0 l P u l1( 0, )0 P2u(u l0, )0 u P2l(u l0, )0 l P u l2 ( 0, )0 . (16)Здесь Piu,Pil (i 1,2) производные Фреше [3]. Приступим к получению оценок для норм u0 u0 1u Pl0 Pl0 и Pl , которые участвуют в неравенстве, обеспечивающем сходимость процесса Нъютона. По определению нормы u0 Pl0 maxx t, P u l1( 0, )0 maxt P u l2 ( 0, ) .0 В силу предположений о гладкостях функции , и f , после некоторых преобразований слагаемых левой части (15) получим u0 , 1/ 2 1, (17) Pl0 M p (t t0 ) где M p положительная постоянная, зависящая от данных задачи. Если мы покажем, что матричное уравнение (16) разрешимо, то этот факт будет эквивалентен существованию и ограниченности по норме обратного оператора 1 u0 Pl0 . Решение матричного уравнения (16) равносильно системе линейных интегральных уравнений Вольтерра второго рода, которая разрешима при значениях t, достаточно близких к t0. Выполняя определенные вычисления доказывается , что операторное уравнение (15) имеет единственное в шаре u u0 1 12h0 l l0 r r1, r1 h0 0 u* Решение * , где h0,0 постоянные, связанные с данными задачи. l Замечание. Решение задачи (3)-(6) в области l t( ) x явно выписывается с помощью функции Грина для указанной области, если только известно решение указанной задачи в области 0 x l t( ). Результаты и их обсуждениеВыделен и исследован класс задач с неизвестной границей, отличающих от известных задач Стефана, Веригина, Флорина и опысывающих процесс фильтрации вблизи новых каналов и водохранилищ с учетом испарения, которое является нелинейной финитной функцией времени и уровня грунтовых вод. Результаты данной работы могут быть использованы специалистами научноисследовательских и проектных институтов, ведущих теоритические и проектноизыскательные работы по предотвращению засоления и заболачивания а также подтопления земель в районе гидротехнических сооружений (крупные каналы, водохранилища и т.п.). Заключение Проведен и применен способ раздельной линеаризации, дающий значительно меньшую погрешность по сравнению обычном способом линеаризации. Решена задача с неизвестной границей, отличающаяся от известной задачи Стефана, которая описывает процесс фильтрации вблизи новых каналов и водохранилищ с учетом испарения. При этом разработанная методика может быть применена к решению подобных задач математической физики. Download 2.25 Mb. Do'stlaringiz bilan baham: |
Ma'lumotlar bazasi mualliflik huquqi bilan himoyalangan ©fayllar.org 2024
ma'muriyatiga murojaat qiling
ma'muriyatiga murojaat qiling