Актуальность темы
В связи с проблемой наблюдающегося глобального изменения климата все более актуальными становятся вопросы, связанные с изучением последствий воздействия изменения метеорологических условий на гидрологические и экологические процессы в водоемах и водотоках при их круглогодичной эксплуатации. Большое значение приобретает рациональное использование водных ресурсов в зимнее время, когда очень важным становится обеспечение достоверности прогноза процесса формирования ледяного покрова и его прочности, позволяющее определить начало и сроки эксплуатации ледовых трасс, проложенных по поверхности водных объектов.
Разработкой методов расчета толщины ледяного покрова занимались многие исследователи: Я. Л. Готлиб (1983), Л. Г. Шуляковский (I960,1969, 1972), В. В. Пиотрович (1968, 1969, 1970), Р. В. Донченко (1971, 1983), Т. В. Одрова (1979), Г. Эштон (1978,1986), А. И. Пехович (1983), Г. А. Трегуб (1994, 1997), И. Н. Шаталина (1990), В. М. Мишен (1997) и многие другие. К настоящему времени предложено большое количество формул и расчетных приемов для определения толщины ледяного покрова пресноводных водоемов и водотоков, опирающихся на основные методы расчета: эмпирический и теоретический методы.
Эмпирический метод основан на отыскании эмпирических связей толщины льда и отдельных факторов, определяющих изменение толщины ледяного покрова. В этом случае расчетные эмпирические соотношения получены на основании относительно тесной корреляции между некоторыми температурными характеристиками и толщиной льда. Как правило, эти характеристики привязаны к конкретному региону.
Теоретический метод основан на интегрировании исходных дифференциальных уравнений, описывающих физическую сущность нарастания толщины льда.
Основой для получения эмпирических формул послужило уравнение для подсчета теплового баланса, так называемое условие Стефана. В 1889 г. Й. Стефан в работе о промерзании грунта поставил и решил следующую задачу (Рубинштейн Л.И., 1967).
Среда, находящаяся в двух фазовых состояниях (жидком и твердом) и проводящая теплоту исключительно посредством теплопроводности, заполняет полупространство х > 0. В начальный момент времени вся она находится при постоянной температуре Ti > 0 °C. На плоскости х = 0 поддерживается постоянная температура Т < 0 °C, под воздействием которой происходит кристаллизация, протекающая изотермически при температуре Т = 0 °C без переохлаждения и с пренебрежимо малым объемным эффектом.
Требуется определить температуры Ui (x, t) и U2(x, t) твердой и жидкой фазы и положение х = y (t) границы раздела фаз.
Подсчет теплового баланса приводит, как показал Стефан, к условию.
Здесь, А — скрытая теплота кристаллизации, отнесенная к единице массыр — плотность образующейся фазык и — коэффициенты теплопроводности соответственно твердой и жидкой фазы.
Эмпирические формулы для определения толщины ледового покрова, полученные на основе решения задачи Стефана, имеют общий вид :
1) где hice — толщина льда, см, ?(—#2) — сумма средних суточных температур воздуха (°С) на высоте 2 м над уровнем водоема за расчетный период времениа — эмпирический коэффициент, определяемый по данным непосредственных наблюдений и отражающий в среднем те условия, которые имели место в период наблюдений (температуру воды, высоту и плотность снежного покрова, глубину водоема, скорость подледного течения воды). Однако ввиду различия этих факторов даже для отдельных участков водоемов и недостаточной продолжительности наблюдений указанный параметр существенно меняется. Поэтому все подобные формулы носят локальный характер. Подробный обзор работ, описывающих получение эмпирических формул различными авторами для конкретных водных объектов представлен в монографии Д. В. Козлова (2000).
Очевидно, что в подобных эмпирических формулах не находят в полной мере своего отражения такие важные факторы, как мощность и физические свойства снежного покрова, интенсивность теплопотока из водной массы, метеорологические условия. Тем не менее, снег за счет своих теплоизоляционных свойств оказывает особое влияние на рост льда, уменьшая влияние отрицательных атмосферных температур.
В своих исследованиях норвежский ученый О. Дэвик (Винников С.Д., Проскуряков Б. В., 1988) к процессу льдообразования подошел с позиций физики, также основываясь на уравнении Стефана. В предложенной им формуле учитывается теплообмен на верхней поверхности снежного покрова, теплофизические свойства снега. Его подход стал классическим и получил широкое развитие. В дальнейших разработках многие ученые уделяли больше внимания процессу теплообмена снежно-ледовой поверхности с атмосферой (В.В. Пиотрович), теплообмену с дном водоема, водным потокам, начальной стадии формирования льда, образованию шуги.
JI.Г. Шуляковский, Р. В. Донченко, С. Д. Винников, Б. В. Проскурякрв, Г. А. Трегуб, И.Н. Шаталина).
На основе обощения этих разработок сформированы методические рекомендации, выпущенные научноисследовательскими и проектными институтами ВНИИГ, ГГИ, Ленгидропроект.
С развитием вычислительных методов математического моделирования и компьютерных технологий последнее время стали чаще обращаться к теоретическим методам, более точно описывающим физические явления процесса льдообразования. Появилась возможность численно решать сложные системы дифференциальных уравнений (Цибульский В.Р. и др., 1993: Васильев О. Ф., Бочаров О. Б., Зиновьев А. Т., 1999: Васильев В. И. и др., 1997: Дебольская Б. И., 2003: Прокофьев В. А., 2004), описывающих процессы льдообразования. Не остались без внимания важные задачи прогнозирования ледотермических режимов водоемов-охладителей тепловых и атомных электростанций, влияние их эксплуатации на температурный и ледовый режим рек (Квон В.И., Филатова Т.Н.). Над задачами моделирования ледового режима бьефов Красноярской ГЭС успешно трудились В. М. Белолипецкий, С. Н. Генова, В. Б. Туговиков, Ю.И. Шо-кин (1991, 1993).
На юге Западной Сибири имеется значительное число минерализированных озер. Как показали проведенные исследования, даже слабая минерализация оказывает существенное влияние на процесс замерзания. А именно, при формировании ледового покрова в минерализированном водоеме в результате образования практически пресного льда в воде перед фронтом кристаллизации формируется слой с повышенным содержанием химических веществ, что существенно влияет на температуру фазового перехода, а тем самым и на весь процесс формирования ледового покрова (Власов В.П., 1968, 1976), что необходимо учитывать при разработке математической модели и выполнении соответствующих расчетов.
Процесс замерзания морских вод достаточно подробно изучен, описан в учебниках по океанологии и статьях (Смирнов Г. Н., 1974: Федоров К. Н., 1976: Саркисян А. С., Демин Ю. Л., Бреховских А. Л., Шаха-нова Т.В., 1986: Багно А. В., Гаращук Р. В., Залесный В. Б., 1996: Яковлев Н. Г., 2003). Гидродинамические и гидрохимические режимы морских вод отличается от режимов соленых озер, которые по солевому составу чаще являются солоноватыми. Поэтому способы описания процесса образования морских льдов не применимы к озерным.
Во всех этих работах, основанных как на эмпирических, так и на теоретических методах, снежный покров в процессе льдообразования принято учитывать через эквивалентный слой: h-ecv = hice ~Ь hsnow—, (3).
Ksnow где hecv — эквивалентный слой снежно-ледового покрова, hmmv — высота снежного покрова на льду, kice и ksnow — коэффициенты теплопроводности льда и снега, соответственно. При этом предполагается линейное распределение температуры в слое льда и снега, что далеко не соответствует действительности.
Снежный покров по высоте имеет неоднородную структуру. Соответственно, изменяются его физические свойства, такие как плотность и теплопроводность. Плотность снега меняется со временем в зависимости от усадки, ветрового воздействия, плотности свежевыпавшего снега. Это приводит к нелинейному распределению температуры в слое снега.
В работе Л. И. Рубинштейна (1967) подробно изложена история развития «проблемы Стефана». Задачи Стефановского типа, помимо изучения ледового режима водоемов, получили широкое распространение в других областях науки. Они успешно используются при моделировании процессов фазовых переходов в бинарных смесях (производство полупроводниковых материалов, очистка методом направленной кристаллизации, металлургическое производство), подробно описанных Н. А. Авдониным (1980), Рубинштейном Л. И., 1967 A.M. Мейрмановым (1986), при описание роста кристаллов (Карслоу Г., Егер Д., 1964: Любов Б. Я., 1969, 1975: Флеминге М., 1977), при описании диффузионного переноса вещества в зоне реакции и.т.д.
Известны различные численные методы решения задачи Стефана, такие как метод сквозного счета, используемый в работе (Бондарев ЭА., Васильев В. И., 1984: Васильев В. И. и др., 1997), методы с явным выделением фронта (Будак Б.М. и др., 1966: Воеводин А. Ф., Леонтьев Н. А., Петрова А. Г., 1982: Петрова А. Г., 1983: Белолипецкий В. М. и др., 1993; Овчарова А. С., 1994, 1995, 1999: Журавлева Е. Н., 1998). В работе (Бондарев Э.А., Попов Ф. С., 1989) приводятся оценки точности наиболее употребляемых приближенных методов решения задачи Стефана.
Целыо работы является разработка математических моделей и эффективных численных методов для исследования процесса формирования ледового покрова водоемов различной минерализации в условиях Западной Сибири.
Основные задачи формулируются следующим образом:
• разработка математических моделей и их численная реализация для исследования процесса формирования ледового покрова в водохранилищах при наличии снежного покрова с учетом минерализации воды;
• разработка численного метода решения двухфазной термодиффузионной задачи Стефана с температурой кристаллизации, зависящей от концентрации примеси;
• проведение расчетов ледотермического режима реальных объектов: Новосибирское водохранилище, озера Чановской системы;
• сравнительный анализ результатов расчетов по разработанным моделям и известным инженерным методикам и сравнение с натурными измерениями.
В работе основное внимание уделено разработке математических моделей для исследования процесса динамики ледяного покрова водных объектов с учетом важных физических факторов: наличие снежного покрова и зависимость температуры фазового перехода «вода-лед» от концентрации солей в воде, построению адекватных разработанным моделям численных методов и созданию реализующих эти методы программ для ПЭВМ, созданию численных моделей для исследования ледотерми-ческих процессов на реальных объектах Сибири (Новосибирское водохранилище, озера Чановской системы).
В первой главе пункт 1.1. посвящен описанию и разработке алгоритма численного решения одномерной двухфазной классической задачи Стефана, моделирующей кристаллизацию чистого вещества (Флеминге М., 1977; Рубинштейн Л. И., 1967). Рассматриваются области Qi (t) = {О < z < f{t)} и Qs (t) = {f (t) < z < Н}, занятые соответственно жидкой и твердой фазами данного вещества с гладкой границей раздела z = /(0- Распределение температуры в областях описывается уравнениями: дт, 2д2т,.
— Ж = а (4) дт, гд2т8 аГ = a'-aJ- (5).
Здесь Ti (t, z), Ts (t, z) — температура жидкой и твердой фаз соответственноfit) — неизвестная граница раздела фазof, a2s — соответствующие жидкой и твердой фазам коэффициенты температуропроводностиtпеременная по времениz — переменная по пространствуН — размер заданной области.
Для замыкания системы задаются граничные условия. На подвижной границе z = fit) выполняются условия сопряжения:
1) условие Стефана, которое описывает тепловой баланс (скачком плотности при фазовом переходе пренебрегаем):
X Т/ L дТ1.
Л psVf = ki dz дТ<
— к' vf = ^ (б) OZ z=fl dt.
2) равенство температуры среды температуре фазового перехода Т/ данного вещества, которая считается известной постоянной величиной. Для воды Tf = Т* = 0 °C.
Здесь Vf{t) — скорость нарастания льдаТ* - температура замерзания (кристаллизации) чистого веществаki, ks — коэффициенты теплопроводности, соответственно для жидкой и твердой средыЛ — скрытая теплота плавленияpi, psсоответствующая средам плотностьls — толщина твердой фазы.
На внешних границах, как правило, считается известной либо температура, либо тепловой поток. На границе z = Н должна поддерживаться температура, при которой происходит кристаллизация. Для воды Ts (t, H)<0°C.
Кроме граничных условий и условий сопряжения формулируются начальные условия. Начальное распределение температуры по толщине можно определить функцией или константой. Для задания положения подвижной границы при малых временных параметрах t ~ At, когда скорость продвижения границы теоретически равна бесконечности, значение толщины зарождающейся твердой фазы определяется из начальной асимптотики скорости движения граница раздела, предложенной А. Н. Тихоновым и А. А. Самарским (1972).
Метод численного решения основан на отображении областей с криволинейными границами в регулярные — метод спрямления фронта (Бу-дак Б.М. и др. 1966). При аппроксимации уравнений полученной системы используется неявная схема с направленными разностями для конвективных слагаемых. Алгебраическая система уравнений решается методом прогонки. Путём вычислительных экспериментов на последовательности сеток показана устойчивость предложенного алгоритма.
Пункт 1.2. посвящен описанию математической постановки и методу численного решения термодиффузионной задачи Стефана. В такой постановке задача давно успешно используется при моделировании процессов фазовых переходов в бинарных смесях (производство полупроводниковых материалов, очистка методом направленной кристаллизации, металлургическое производство), подробно описанных Б.Я. Любо-вым (1969, 1975), Н. А. Авдониным (1980), A.M. Мейрмановым (1986).
Здесь рассматривается двухфазная термодиффузионная задача Стефана о замерзании раствора. При образовании пресного льда перед фронтом кристаллизации образуется слой с повышенным содержанием примеси, что влияет на температуру замерзания.
Область Qi (t) = {0 <? < fi (t)} содержит раствор соли (солей). Распределение температуры и концентрации примеси описывается уравнениями: основные обозначения сохранены, C (t, z) — концентрация примеси, d — коэффициент диффузии соли в воде. В образующейся твердой фазе Qs (t) — {fi (t) < z < /2^)} предполагаем отсутствие диффузии примеси, а распределение температуры аналогичное: dTs 2fQ. dt ~ ^ dz2 ' W.
В области, содержащей две фазы, имеются две подвижные границы, положение которых описывается следующими формулами: fi (t) = li = Н — kpls — граница между жидкой и твердой фазой,.
2СО — fi (t) + 18- внешняя граница, мало перемещаемая за счет разности плотностей твердой и жидкой фаз, то есть плавучести льда (Од-рова Т.В., 1979). Здесь li, ls — толщина слоя жидкой и твердой фаз, соответственнокр =.
На внешней границе z = /2СО происходит охлаждение раствора. На подвижной границе z = fi (t) выполняются условия сопряжения:
1) условие Стефана — уравнение (6);
2) равенство температур жидкой и твердой фаз неизвестной температуре замерзания Т/, определяемой уравнением:
Т, = Т* - 7С/- (10).
3) баланс массы растворенного в воде вещества (Васильев В.И. и др., 1997):
На границе z = 0 для температуры и примеси определены: условие изоляции (отсутствие притока тепла и примеси), или закон теплообмена со средой, или постоянное значение.
Tf (t) — в данном случае неизвестная температура фазового переходаCf (t) — искомое значение примеси на границе раздела фаз- 7 — равновесный коэффициент распределения примеси.
Кроме граничных условий и условий сопряжения формулируются начальные условия: распределение температуры по толщине в двух фазах, значение концентрации примеси (оно, как правило постоянно С = Со).
Начальные значения V/ и ls определялись путем численного эксперимента и согласованы с асимптотикой для однородной жидкости (п. 1.1).
Для численного решения задачи использовались метод спрямления фронта. Аппроксимация уравнений проводилась как и в п. 1.1. В каждой из областей была построена равномерная сетка. Так как температура замерзания Т/ и солености С/ на границе z = fi (t) заранее неизвестны, то здесь классический метод прогонки использовать невозможно. Кроме того, в связи с тем, что начальная стадия процесса приводит к возникновению больших градиентов температуры и солености, то применение методов ловли фронта в узел сетки или сквозного счета вызывает большие трудности. Поэтому для численной реализации задачи была разработана модификация метода (Воеводин А.Ф., Леонтьев Н. А., Петрова А. Г., 1982), основанного на методе встречной прогонки, позволяющего эффективно использовать безитерационные (точные) методы решения разностной задачи, несмотря на ее нелинейность.
Идея этой модификации состоит в том, что значения Т/, Ts, С на границе z = f (t) 1 записываются через прогоночные коэффициенты. После подстановки полученных выражений в уравнения (6), (10), (11) имеем квадратное уравнение относительно значения концентрации примеси на границе раздела фаз С/:
1АТС) + [Ас — BT) Cf + ВС = 0. (12).
Результатом решения этого уравнения является один удовлетворяющий физическим условиям корень: с = Вт-Ас ЛуЦАс — Вт)2 — 4jATBc. f 27 Ат ' ^ } где выражения Ат, Вт, Ас, Вс находятся однозначно через прогоночные коэффициенты.
Получив значение примеси на границе раздела z = fi (t), из уравнения (10) находим температуру фазового перехода, соответствующую температуре жидкой и твердой фаз на этой границе. Зная граничные условия, восстанавливаем значения температуры и примеси на новом шаге по времени.
Вторая глава содержит численные модели, описывающие динамику роста ледового покрова водоемов, построенные на основе классической и термодиффузионной задачах Стефана.
Для моделирования ледового покрова естественных водоемов недостаточно описать процесс фазового перехода воды в лед. Необходимо точно учитывать наиболее важные физические факторы, оказывающие влияние на рост льда. Такими факторами являются, например, снежный покров и минерализированность водоема.
Снег за счет своих теплоизоляционных свойств оказывает особое влияние на рост льда, уменьшая влияние отрицательных атмосферных температур. Снежный покров по высоте имеет неоднородную структуру. Соответственно, изменяются его физические свойства, такие как плотность и теплопроводность. Плотность снега меняется со временем в зависимости от усадки, ветрового воздействия, плотности свежевыпавше-го снега. Это приводит к нелинейному распределению температуры в слое снега, что и было отмечено при исследовании водного и ледового режимов на Новосибирском водохранилище в 1982 г. экспедицией dt dz2'.
9Т{се о ^ ^г'се = a.
H И С И им. В. В. Куйбышева (Козлов Д.В., 2000). Следовательно, появилась необходимость модифицировать математическую модель и ввести новую область, описывающую процессы в снежном покрове.
В моделях используется описанная в первой главе базовая система уравнений двухфазной задачи Стефана. С появлением третьего слоя <? < /з (£), содержащего снег, образуется еще одна подвижная неизвестная граница z = /з (£)d^water т д Twafer (л л.
Р water ср ^ — & water '.
ЭС (Л d-^, (15) dt ~&trade-гсе dz2 ' ^ dTsnow д / dTsnow. .
CpPsnow — ^ I tisnow) 5 V-1-') h (t) = /2 W + Isnowi hnaw (t) =? ln (l + bl*snow{t)),.
T-water 1 Tice, Tsnaw — температура в слое воды, льда, снега. psnow = - формула Абе, b = 1.255, ksnow = 2.9 • 106/92nou- + 0.043 — зависимость, предложенная В.В. Пиатро-вичем (1968) для нахождения теплопроводности снега, psnow ~ плотность снега, ро — плотность свежевыпавшего снега, lsrww ~ толщина слоя снега, l*snaw — толщина слоя свежевыпавшего снега.
На дне водоема при z = 0 принимаем для температуры воды и примеси постоянные значения или условия изоляции. На границе 2 = f (t) — условие равенства температур и условия сопряжения (6), (10), (11). На границе z = /г (^) задается условие равенства тепловых потоков со стороны слоя льда и слоя снега и равенство температур. Значение температуры между льдом и снегом рассчитывалось в ходе решения. Температура снега на внешней границе z = fs (t) задавалась равной атмосфер ной температуре, измеренной на расстоянии 2 м над землей — условие, • удовлетворяющее зимнему метеорологическому режиму Сибирского региона (Пиотрович В.В., 1963).
В связи с тем, что рассматриваемые водоемы являются слабопроточными, то конвективные перенос субстанций не учитывается.
Метод решения изложенной задачи аналогичен методу, описанному в первой главе.
Для пресного водоема задача упрощается. При С = 0 температура фазового перехода Т* = 0. Уравнения (10), (11) и (15) не рассматриваются.
Также для пресного водоема предложено решение проблемы начальной асимптотики для скорости продвижения фронта, за счет введения новой переменной S (t) = lfce (t) (Атавин А.А., Гранкина Т. Б., 2004). После такой замены начальное значение /г-се можно положить равным 0.
Третья глава содержит описание общепринятых методик расчета тол> щины ледового покрова водоемов, результаты расчетов по предложенным математическим моделям для водоемов Западной Сибири: Новосибирского водохранилища — пресный водоем, озера Чаны — система озер и плесов различной минерализации. Также проведен сравнительный анализ результатов расчетов по инженерным методикам и математическим моделям и исследование области применимости приближенных инженерных методов для решения задачи динамики ледового покрова.
В заключении описаны полученные результаты:
• разработана математическая модель для расчета процесса формирования ледового покрова водоемов различной минерализации с учетом таких важных физических факторов, как наличие снежного покрова и зависимость температуры фазового перехода «вода-лед» от концентрации солей в воде;
• разработан и реализован на ЭВМ численный метод решения многофазной термодиффузионной задачи Стефана, с температурой фазового перехода, зависящей от концентрации примеси;
• разработаны и реализованы на ЭВМ методы расчета процесса формирования ледового покрова водоемов различной минерализации;
• создана численная модель и проведены методические расчеты для исследования ледотермических процессов на реальных объектах Западной Сибири (Новосибирское водохранилище, озера Чановской системы), проведен сравнительный анализ результатов расчетов с приближенными методами и сравнение с натурными натурными данными, что показало хорошую точность и эффективность расчетов по созданной модели.
Научная новизна. Разработан численный метод для решения задач типа Стефана. Создана и апробирована математическая модель, а также комплекс программ для исследования процесса динамики ледяного покрова водоемов с различной минерализацией. Впервые в математической постановке задачи учтены пространственное изменение плотности и температуры снежного покрова, а также зависимость температуры фазового перехода от минерализации водной среды. Разработаны численные модели, позволяющие исследовать ледотермические процессы реальных объектов Западной Сибири.
Практическая значимость. Полученные результаты могут быть использованы для прогноза динамики роста ледового покрова на водоемах Западной Сибири в условиях изменения климата, уровня воды в водоемах и степени минерализации водной среды. Предлагаемые математические модели, методы численного решения и программное обеспечение могут быть использованы при проведении экологического мониторинга водных объектов и организации ледовых переправ. Также, разработанные модели и метод могут быть использованы при исследовании физических, биологических и технологических процессов, приводящих к необходимости численного решения задач фазового перехода вещества, содержащего примесь.
Объект и методы исследований. Объектом исследований являются процессы формирования ледового покрова Новосибирского водохранилища и озер Чановской системы. Выполнение исследований базируется на математическом моделировании и численных методах решения термодиффузионной задачи Стефана.
На защиту выносится:
• математическая модель для расчета процесса формирования ледового покрова водоемов различной минерализации с учетом физических факторов, таких как наличие снежного покрова и зависимость температуры фазового перехода «вода-лед» от концентрации солей в воде;
• численный метод решения многофазной термодиффузионной задачи Стефана;
• методы расчета процесса формирования ледового покрова водоемов различной минерализации;
• результаты расчетов, выполненные для исследования ледотермиче-ских процессов на ряде реальных объектах Западной Сибири.
Апробация работы. Результаты докладывались на 1-ой межвузовской научно — практической конференции студентов и молодых ученых Сахалинской области (Южно-Сахалинск, 1997 г.), XXXII научно — методической конференции преподавателей ЮСГПИ (Южно-Сахалинск, 1998 г.), IV школе-семинаре «Математические проблемы механики сплошных сред» ИГиЛ (Новосибирск, 2000), XXXIX Международной научной студенческой конференции «Студент и научно — технический прогресс» (Новосибирск, 2001 г.), Международной конференции «Вычислительные технологии и математическое моделирование в науке, технике и образовании «(Алматы, 2002 г.), Всероссийской конференции «Математические методы в механике природных сред и экологии» (Барнаул, 2002 г.), 61-й научно — практической конференции профессорско — преподавательского состава НГАСУ (Новосибирск, 2004 г.), 17-м Международном симпозиуме по льду (Санкт-Петербург, 2004 г.), Международной конференции по измерениям, моделированию и информационным системам для изучения окружающей среды «ENVIROMIS-2004» (Томск, 2004 г.), VI конференции «Динамика и термика рек, водохранилищ и прибрежной зоны морей» (Москва, 2004 г.), 3-й и 4-й конференциях молодых ученых ИВЭП СО РАН (Барнаул, 2004, 2005 гг.), научной конференции «Фундаментальные проблемы изучения и использования воды и водных ресурсов» (Иркутск, 2005 г.), научных семинарах лаборатории прикладной и вычислительной гидродинамики Института гидродинамики им. М. А. Лаврентьева и Новосибирского филиала ИВЭП СО РАН.
Автор выражает глубокую благодарность научному руководителю проф. д.ф.-м.н. А. Ф. Воеводину за постановку задачи, постоянно внимание к работе, к.т.н. А. А. Атавину и д.ф.-.м.н. В. И. Квону за плодотворное сотрудничество, ценные советы и организационное содействие. Также выражает признание сотрудникам Новосибирского филиала ИВЭП СО РАН и Лаборатории прикладной гидродинамики ИГиЛ СО РАН за помощь в работе, консультации, своевременную критику и дружескую поддержку.
Заключение
.
В диссертации получены следующие основные результаты:
• разработана математическая модель для расчета процесса формирования ледового покрова водоемов различной минерализации с учетом таких важных физических факторов, как наличие снежного покрова и зависимость температуры фазового перехода «вода-лед» от концентрации солей в воде;
• разработан и реализован на ЭВМ численный метод решения многофазной термодиффузионной задачи Стефана, с температурой фазового перехода, зависящей от концентрации примеси;
• разработаны и реализованы на ЭВМ методы расчета процесса формирования ледового покрова водоемов различной минерализации;
• создана численная модель и проведены методические расчеты для исследования ледотермических процессов на реальных объектах Западной Сибири (Новосибирское водохранилище, озера Чановской системы), проведен сравнительный анализ результатов расчетов с приближенными методами и сравнение с натурными натурными данными, что показало хорошую точность и эффективность расчетов по созданной модели. ч.
