МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ МИКРОКЛИМАТА В СОЛНЕЧНОЙ ТЕПЛИЦЕ ТРАНШЕЙНОГО ТИПА Текст научной статьи по специальности «Математика»

Статья поступила в редакцию 28.05.10. Ред. рег. № 818

The article has entered in publishing office 28.05.10. Ed. reg. No. 818

УДК 621.383; 621.472

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ МИКРОКЛИМАТА В СОЛНЕЧНОЙ ТЕПЛИЦЕ ТРАНШЕЙНОГО ТИПА

А.М. Пенджиев

Туркменский политехнический институт Туркменистан,744032, Ашхабат-32, м. Бекрова, Солнечный 4/1 Тел. + (312)37-09-50), e-mail: [email protected]

Заключение совета рецензентов: 18.06.10 Заключение совета экспертов: 28.06.10 Принято к публикации: 05.07.10

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

Ключевые слова: математическая модель, солнечная теплица, микроклимат.

MATHEMATICAL MODELLING OF THE MICROCLIMATE IN SOLAR HOT-HOUSE DEEP OF TYPE

A.M. Penjiyev

Turkmen Polytechnical Institute Solar 4/1, m. Bekrova, Ashabad-32, 744032, Turkmenistan Tel.: (312) 37-09-50, e-mail: [email protected]

Referred: 18.06.10 Expertise: 28.06.10 Accepted: 05.07.10

In article are considered mathematical modelling of a microclimate of a solar hot-house deep type for cultivation of tropical and subtropical plants.

Keywords: mathematical model, a solar hot-house, a microclimate.

Актуальность задачи

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

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

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

расчету Д.А. Куртенера, А.Ф. Чудновского, Г.Г. Се-мийкиной, О.Л. Решетина [1-5], А.Н. Ануфриева, И. А. Кожинова, Г.М. Позина [6], Л.Г. Прищепа [7], Ю.К. Росса [8], Ю.Л. Раунера [9], Х.Г. Тооминга [10], Б. Хайрутдинова [14], Р.Б. Байрамова, Л.Е. Рыбаковой, Н.Г. Гурбанова, [11-13] и многих других.

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

С углублением траншеи появляется неоднородность освещения почвы и стен. Все это, в свою очередь, усложняет задачу расчета.

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

Рассчитывая математическую модель, условно перемещаем культивационное сооружение в различные климатические зоны: в северную точку - Дашо-вузская область (велаят), Конеургенчский район (эт-рап); в восточную - Лебапская область (велаят), Сердарабатский район (этрап); в центральную -Ахалская область (велаят), Ашгабатский район (этрап); в юго-восточную - Балканская область (велаят), Етрекский район (этрап), - рассматривая условия в теплице в данных зонах в течение года.

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

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

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

Основные характеристики защищенного грунта

Микроклимат теплиц, в частности, тепловой, водный, световой режимы системы почва - воздух, определяются, с одной стороны, погодой, с другой -теплотехническими и конструктивными характеристиками защищенного грунта и свойствами агрофи-тоценоза.

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

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

Существует множество типов культивационных сооружений: наземные малогабаритные теплицы (парники), теплицы-оранжереи, блочные, арочные теплицы и теплицы траншейного типа. Они представляют собой постройки со светопрозрачным укрытием из пленки или стекла. Теплица траншейного типа - сооружение, углубленное в землю и закрытое сверху светопрозрачным ограждением. Итак, защищенный грунт реализуется разными по конструкции, размерам и площади строениями, общим признаком которых является только прозрачная крыша [1, 6, 7, 11-13].

Конструктивные особенности культивационного сооружения можно охарактеризовать с помощью коэффициентов.

Коэффициент ограждения п:

П = ^огр/^, (1)

где ^огр - площадь поверхностного ограждения; ^ -площадь культивационного сооружения [1, 2].

Коэффициент ограждения показывает, насколько площадь ограждения превышает площадь культивационного сооружения. Он зависит от того, насколько развита поверхность ограждения, изменяется от 1 до 3.

Его значение зависит от типа сооружения, к примеру, малогабаритная теплица с однослойным пленочным покрытием без обогрева имеет коэффициент ограждения 1,6. Такая же теплица, но при использовании теплозащитных пленочных штор имеет коэффициент 2,6; теплица траншейного типа - 1,1, самый наименьший коэффициент ограждения - у малогабаритных теплиц.

Другой важной характеристикой является £ - коэффициент объема, определяющийся по следующей формуле:

£ = (2)

где У - объем культивационного сооружения [1, 2].

Обычно коэффициент равен средней высоте культивационного сооружения кср. Действительно, коэффициент объема для сооружения кср У = ксрК Подставим значение объема, получаем £ = кср. Коэффициент объема для разных теплиц различен [1, 2]. К примеру, для углубленного парника с укрытием -0,25, малогабаритного наземного сооружения - 0,6, многоблочной теплицы с однослойным стеклянным

International Scientific Journal for Alternative Energy and Ecology № 7 (87) 2010

© Scientific Technical Centre «TATA», 2010

ограждением - 2,7, ангарной с однослойным стеклянным ограждением - 3,0, для теплицы траншейного типа - 4,8.

Конструктивная характеристика и

и = e/4,

(3)

где е - наименьший горизонтальный габаритный размер культивационного сооружения (ширина) [4].

Коэффициент и показывает, насколько ближние стены удалены друг от друга [4]. К примеру, малогабаритные наземные сооружения имеют и = 1,5-3,0; арочные теплицы и = 3-5; малоблочные теплицы и = 30-50; траншейные теплицы и = 1.

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

ю = FCT /F„.

(4)

где Fст - площадь стен траншеи; Fп - площадь почвы (в данном случае она равна ю = 0,94).

Коэффициент затененности почвы:

^п = Fп7Fп ; (5)

коэффициент затенения стены:

^ = F°/Fet , (6)

1 ' 0 7-^0

где и гст - площадь затененной почвы и стен.

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

Расчеты влияния ориентации и формы конструкции гелиотеплицы траншейного типа на величину энергии прямой солнечной радиации, проникающей на внутреннюю поверхность, рассмотрены в ряде работ [15, 16].

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

Общая схема математического моделирования микроклимата солнечной теплицы траншейного типа

В работах [1-4, 6-9, 11-14, 17, 18] рассматривали методику математического моделирования динамически сложных систем.

Она состоит из следующих этапов: 1 - постановка задачи; 2 - концептуализация; 3 - спецификация; 4 -идентификация; 5 - эксперименты; 6 - реализация модели; 7 - проверка модели; 8 - исследование модели; 9 - оптимизация; 10 - заключение. Остановимся на этих вопросах.

1. Постановка задачи

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

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

Первая - построение динамической модели системы: растение - почва - воздух, описывающей как параметры микроклимата, так и ростовые процессы в зависимости от температуры.

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

2. Концептуализация

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

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

Физическая картина процесса солнечной теплицы представлена на рис. 1.

Рис. 1. Схема преобразования поступающей солнечной радиации Fig. 1. The scheme of transformation of arriving solar radiation

Рис. 2. Одномерная расчетная модель теплового баланса

солнечной теплицы траншейного типа Fig. 2. One-dimensional settlement model of thermal balance of a solar hot-house deep type

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

Для описания математической модели многомерная задача пока недостаточно разработана. Это связано, с одной стороны, с недостаточной изученностью анизотропности тепломассопередачи в культивационном сооружении, с другой - с возможностью сравнительно простого решения поставленной задачи с помощью одномерной математической модели [1-4, 6-9, 11-15, 17, 18].

Детально разработанная одномерная модель, в основу которой положено представление об объекте как о многослойной системе с внутренними экранами, имитирующими поверхности листьев растений, боковых стен, освещенных и неосвещенных. Эта

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

Количественное обоснование одномерной модели теплицы траншейного типа можно представить в грубо описанном виде (рис. 2).

3. Спецификация

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

На второй и третьей стадиях развития агрофито-ценоза состав множества переменных состояния расширяется и в него входят также параметры микроклимата в посеве.

За основу математической модели расчета температурного режима теплицы траншейного типа принята методика нестационарного расчета наземных солнечных теплиц [1-4, 6, 7, 11-14, 17, 18].

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

1. Метеорологические - интенсивность солнечной радиации, температура наружного воздуха и его влажность, скорость ветра, облачность, атмосферные осадки, прозрачность неба, температура подстилающей почвы и др.

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

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

4. Тепломелиоративные системы сооружения, почвенного и воздушного отопления, вентиляционного устройства, полива, аккумулирования тепла и др.

5. Биометрические, тепло- и массофизические характеристики агрофитоциноза и др.

4. Идентификация

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

International Scientific Journal for Alternative Energy and Ecology № 7 (87) 2010

© Scientific Technical Centre «TATA», 2010

при котором записывается система баланса энергии и массы, связующая основные показатели микроклимата с входными переменными.

Определение температуры почвы тп, стенки тст, воздуха тв сводится к решению задачи

Ти = fj(tb, тСт, мь м2, м3, м4, м5);

тСт = f2(tb, тп, mi, м2, м3, м4, м5); тв = fn(tn, тСт, mi, м2, м3, МА, м5),

(7)

где М1 - совокупность метеорологических факторов; М2 - совокупность теплотехнических характеристик сооружений; М3 - совокупность теплофизических характеристик воздушной среды, ограждения, почвы, стенки; М4 - совокупность применяемых тепломе-лиоративных мероприятий; М5 - совокупность биометрических, тепло- и массофизических характеристик агрофитоценоза.

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

Уравнение теплового баланса воздушного пространства (рис. 1, 2) для промежутка времени ¿Г с учетом обогрева можно записать в следующем виде:

dQоб + ¿0рЕ - - - ^бст - dQп - - й'бр = 0, (8)

где dQоб - тепловыделение системы отопления;

- поток тепла солнечной радиации, поступающей в культивационное сооружение; ¿0вт - количество тепла, отдаваемое окружающей среде в результате воздухообмена и теплопередачи через ограждения; ¿0ст = ¿Ост + ¿0с0т - поток тепла в стенку; ¿0п = ¿01 + + - поток тепла в почву; (ст), - поток

тепла на освещенную и неосвещенную части стенки и почвы; ¿0а - поток тепла, аккумулируемый в воздухе за время ¿т; ¿0р - поток тепла на растительный покров (при условии, если теплица с растениями).

Уравнения теплового баланса поверхности почвы для освещенной и неосвещенной частей при отсутствии растительного покрова будем записывать отдельно:

dQn = dQn - - dgr -

dQ;0 = Q + dQ;ta + dQ°K + dQ

(9) (10)

где ¿0пЕрс - суммарный поток тепла приходящей солнечной радиации в почву; ¿0п - поток тепла приходящей рассеянной радиации на неосвещенную часть почвы; - поток тепла, излучаемого почвой; - затраты тепла на испарение; - по-

ток тепла, конденсированного поверхностью почвы;

- конвективный тепловой поток почвой за тот

же промежуток времени ¿т, здесь индексы; «п», «*», «0» - соответственно, почвы, освещенной и неосвещенной поверхностей.

Уравнение теплового баланса поверхности стенки записывается аналогичным образом, как уравнение теплового баланса для почвы, освещенной и неосвещенной ее части, оно будет иметь вид:

¿ост = ¿от - ¿о;; - ¿остк - ¿ости; (11)

¿ост = ¿0сртс + ¿от + ¿0с0тк + ¿0Гон; (12)

¿0с;и ~ 0 - не учитывали.

Уравнение теплового баланса ограждения можно записать в следующем виде:

¿Оогр = ¿0вл + ¿ок + ¿0вкон + ¿орс - ¿0нл - ¿онк - ¿они. (13)

Ясно, что при кратковременных периодах осадков учитывать ¿он нет смысла.

Здесь - поток тепла солнечной радиации,

поглощаемой стеклом. Индексы «н» и «в» обозначают наружную и внутреннюю поверхности ограждения.

Подставляя составляющие выражения в уравнения (8) - (13), получаем:

Fn* Эт П

dт В n FcbуBZ дХ

-X

F0 дТ 0 1 ;Т U1 n

FCbYbZ дх

-X

F* дт *

СТ СТ

FCb Y BZ дХ П

-X

Fc0T дт Ст

FCb y bZ дх

= - С- (Ai + a н + в )т orp +3(т),

где

^ Н ТНВ+

m , а3 + В3

+с—z q + С—г п + ;

СВ Y BZ р СВ Y BZ

1 р (i - exp (-Qt) d т)

(14)

9об =

FCBYв? аобFo60еХР (-QT)dT

где доб взята из работы авторов [1-4, 6-7, 11-13, 17, 18], £ - кратность воздухообмена, с; - коэффициент

объема, св - теплоемкость воздуха, А1-34, 5ы5 - коэффициенты излучения и испарения.

Уравнение потока тепла для освещенной поверхности почвы и стен [(9), (11)]:

F * 1 п дТ п

F дХ Х=0

-(1 - s, )т П ( + а п + ва )+

>

+ (1 - S; )тв (ап + в5) = -а4тогр + 3; (т); (15);

-(; - S2 )тСт (( -аст) +

F* дТ*

Л СТ СТ

Л СТ

F дХ

+ (; - S2 )ТваСТ = А18ТОГР +32(т) .

(16)

Уравнение потока тепла для неосвещенной поверхности почвы и стен [(10), (12)]:

Л п

F0 1 п дТ О

F дХ

+ S,Tп ((-ап -В8) + STВ (ап + В7) = = -А12SjTогр + 3з (т); (17)

F0 дТО

СТ СТ

Л,

+ S2T(0t (А20 аСТ В11 ) +

д . дт дт , ч

--Л--с— = g (х, т).

дх дх дт

(19)

Итак, мы имеем систему из шести уравнений, содержащую пять неизвестных величин: искомые температуры Тв, Ть Т2, Тз, Т4.

Кроме того, поскольку уравнения системы содержат производные по х и т:

(д2Т. дТ. дТ. дТв Л . ч

I—г, —, —, — I (' = 1,2,3,4), ^ дх. дх. дх. дт ^

то для замыкания системы запишем следующие равенства:

4=0 = F (х,); (20)

TU = T; (21)

Tb L=0 = Tbo (i = 1,2,3,4). (22)

Рассмотрим аналитическое решение системы уравнений (20) - (22 ) и уравнений

^+MbлдП d т дх.

= 3(т); (23)

дх,.

- nt

+ RTb =3,. (т);

(24)

(i = 1,2,3,4).

^ дх

+ ^ («СТ + Во ) = ^24^гр + З4 (т) , (18) где 3 равно:

3, (т) = (& -1)[(1 -/1 )ф9р + АТт + атОт + а + В6 ]; З2 (т) =

= (& -1)[(1 - /2 )ф?р + а15тст + а16тП + а17тП + а19 ]; 3з (т) = -& [АюТ'ст + АпТст + аз + В9 ] ;

34 (т) = -[А21ТСТ + А22ТП + А23ТП + А25 + В12 ] .

К уравнениям (15) - (18) присоединим уравнение эффективной теплопроводности:

Для решения системы уравнений (20) - (24) используем возможность сведения ее к математической задаче, рассматриваемой в работе О.Л. Решетина [3], когда отыскивается решение уравнения:

д .( )д0 ( )д0 ( ) -—Л(х)-—с(х)—= я (х, т). (25) дх дх дт

В отличие от работы [3], в данном случае имеются четыре неизвестные величины вместо одной. Поэтому необходимо решить несколько уравнений типа (25):

<х'>!т=«м (2б)

при краевых условиях: д0.

-Л-

Эх,

- N0,L=o + R0B =3,. (т); (27)

= 0;

0. L=o = F (х);

(. = 1,2,3,4),

(28) (29)

где 0В - решение линейного дифференциального уравнения

d 0 4 д0.

£+a0bл. ^

1 I

(. = 1,2,3,4 ),(30)

удовлетворяющего начальному условию

0в 1т=о = 0. (31)

Величины 0! и 0В определяются по формуле:

0. = Т - Т; (32)

0в = Тв - Тво (, = 1,2,3,4). (33)

Решение уравнения (33) при условии (30) известно:

т

0в= /3(8)ехр [-к (т-8)] а 80

4 т д0

-У Л,. ехр[-к(т-8)а8 ( = 1,2,3,4). (34)

■^Т * дх х,=0 -1

, = 1 0 и л !

International Scientific Journal for Alternative Energy and Ecology № 7 (87) 2010

© Scientific Technical Centre «TATA», 2010

Если подставить выражение (34) в уравнение (27), то получим

где

-X

де,,

+ N е,

Л , J

' дх,

де

дх,

exp [-h (т - 8)] d 8 =

-Я,

дх,

+ Nе = kt (t)b(t) (' = 1,2,3,4) (37)

ев = Js^^ [-h (т-8)] d 8-¿

В' (т)

Я

(39)

0b = ¿U'V [t, Bt (t )] K'v (x,),

(40)

= К, (т); (35)

т

К (т) = -з,. (т) +13 (8) ехр [-И (т - 8)] d 8

0

(, = 1,2,3,4). (36)

В результате преобразований, исходная задача трансформируется в следующую: требуется найти решение уравнения (25) при граничных условиях (35) и (28) и начальном условии (29).

Далее в работах [1-3] показано, что поскольку решение уравнения теплопроводности при интегрированных граничных условиях существует и оно единственное, то задачу можно разделить на две.

Первая состоит в решении уравнения (25) при соответственно граничном условии на поверхности

U v = U'0v exp Ы; -т)+ +J{g(8) + [K (8) + Bi(8)]* Kiu[J exp (т - 8) ] d8 ;

(41)

i

u0l„ = J F (x, t) Km (x)C (x )dx; (42)

§ги=/8г (X, т)Кт (х)С (х^х . (43)

0

Функция Кш (х) - собственные нормированные функции граничной задачи Штурма - Лиувилля [2]:

d . dKl 2^

— К —= -^2K;

dx dx.

K dx,

+ N,K, = 0;

x,=0 1 1 U = 0

K, = 0.

1 Ix, =1

(44)

(45)

(46)

Решение уравнения (44) - (46) приводится в работах [1-3].

Составляющие для ограждения выглядят:

и краевых условиях (28) - (29).

Поскольку граничные условия (37) взаимосвязаны, первая задача решается итерационно, например, методом хорд, постоянных производных и т.д.

В результате решения этой задачи можно получить искомую величину как функцию В ,(т), то есть 91 [, Бг(т)]. Вторая задача сводится к нахождению функции Бг(т) из интегрального уравнения

т I д9 1

В1. (т) = Ш/к^ |ехр\_-И(т-8)]¿8 . (38)

Подставляя найденное решение интегрального уравнения (35) в решение первой задачи, получим искомый результат. Далее из (34) и (38) следует, что искомая величина 9в может быть получена по следующей формуле:

dT

огр

= (1 - m)qp + Т (A26 - Ai - aн + В - aв - Д4) +

огр dT v огр

n0

Решение первой задачи разыскивается с помощью метода конечных интегральных преобразований. Искомая функция 9г [, /, Бг(т)] находится в виде ряда:

+тв («в + вп) - а28т 0 + (1 - ф)9р + атв + А + антнв +

+в2тнв + В + Аэтт; - ах + а^ - АЪ1 - в,5; (47)

В итоге решение уравнений (14) - (18) имеет вид:

9в = }ф(8)ехр[-й(т-8Ж8-£ (48)

0 1=1 Шг

т

Т0р = ехР-(т-т 0)Т0 +

с

огр

+ — / ехр т (т - 8)d8 —— / (т - 8)Р(8^8 . (49)

огр 0 и огр 0

В следующей статье «Математическая модель теплотехнических расчетов микроклимата траншейной солнечной теплицы» мы рассмотрим разделы 5 - эксперименты; 6 - реализация модели; 7 -проверка модели; 8 - исследование модели; 9 -оптимизация; 10 - заключение.

Список литературы

1. Куртенер Д. А., Чудновский А.Ф. Расчет и регулирование теплового режима в открытом и защищенном грунте. Л.: Гидрометеоиздат, 1969.

+

x =0

x. =0

2. Куртенер Д.А., Чудновский А.Ф. Агрометеорологические основы тепловой мелиорации почв. Л.: Гидрометеоиздат, 1979.

3. Куртенер А.Д., Решетин О.Л., Семикина Г.Г., Чудновский А.Ф. Метод расчета температуры почвы с раздельным учетом временных изменений метеорологических параметров. Сб. трудов по агрофизике. 1970. Вып. 26. С. 16-27.

4. Куртенер Д.А., Решетин О.Л., Чудновский А.Ф. Решение уравнений теплопроводности при переменном коэффициенте переноса. Сб. трудов по агрофизике. 1970. Вып. 26. С. 80-99.

5. Чудновский А.Ф. Теплофизика почв. М.: Наука, 1976.

6. Ануфриев Л.Н., Кожинов И.А., Позин Г.М. Те-плофизические расчеты сельскохозяйственных зданий. М.: Стройиздат, 1974.

7. Прищеп Л.Г. Эффективная электрификация защищенного грунта. М.: Колос, 1980.

8. Бихеле З.Н., Молдау Х.А., Росс Ю.К. Математическое моделирование транспирации и фотосинтеза растений при недостатке почвенной влаги. Л.: Гидрометеоиздат, 1980.

9. Раунер Ю.Л. Тепловой баланс растительного покрова. Л.: Гидрометеоиздат, 1986.

10. Тооминг Х. Г. Солнечная радиация и формирование урожая. Л.: Гидрометеоиздат. 1977.

11. Байрамов Р.Б., Рыбакова Л.Е., Гурбанов Н. Упрощенная методика теплового расчета гелиотеп-лицы с учетом нестационарности ее работы // Гелиотехника. 1973. № 3. С. 45-49.

12. Байрамов Р.Б., Мезилов А., Гурбанов Н., Рыбакова Л.Е. Аналитические исследования нестационарного теплового режима гелиотеплицы // Изв. АН ССР, серия ФТХ и ГН. 1973. № 3. Ашхабад: Ылым. С. 29-32.

13. Курбанов Н., Курбанова Г. Количественное описание температурных режимов культивационных сооружений, обогреваемых солнечной радиацией. Ашхабад: ТуркмеНИИТИ, 1983.

14. Хайрутдинов Б. Разработка, исследование и внедрение гелиотеплицы-сушилки с подпочвенным аккумулятором тепла. Автореф. дисс. докт. тех. наук. Ашхабад, 1990.

15. Пенжиев А.М. Разработка, создание и исследование гелиотеплицы траншейного типа для выращивания кофейных деревьев. Дисс. на соиск. уч. степ. канд. тех. наук. Ашхабад, 1986.

16. Пенджиев А.М. Агротехника выращивания дынного дерева (Carica papaya L.) в условиях защищенного грунта в Туркменистане. Автореф. дис. уч. степени докт. наук. М., 2000.

17. Пенжиев А.М. Математическое моделирование теплотехнических расчетов микроклимата и агроклиматическое районирование гелиотеплицы // Гелиотехника. 2001. № 3. Ташкент: Изд-во «Фан».

18. Табуньщиков Ю.А. Расчеты температурного режима помещения и требуемой мощности для его отопления или охлаждения. М.: Стройиздат, 1981.

19. Рыбакова Л.Е., Пенжиев А.М. Рекомендации по выращиванию кофейных деревьев в условиях солнечной теплицы. Ашхабад: ТуркменНИИНТИ, 1990.

International Scientific Journal for Alternative Energy and Ecology № 7 (87) 2010

© Scientific Technical Centre «TATA», 2010