Задача Стефана о фазовом переходе
Дипломная работа - Математика и статистика
Другие дипломы по предмету Математика и статистика
4) и условие (1.56) на фазовой поверхности. Первая часть утверждения следует немедленно, если учесть свойство ? - функции ?(x)=0 при х?0. Для доказательства второй части рассмотрим точку Р на поверхности фазового перехода Ф=0 (см.рис.1.1)
Рис
Проведем нормаль к поверхности и касательную в этой точке и построим цилиндр достаточно малого объема с осью совпадающей с нормалью, симметричной относительно касательной и образующей, параллельной к нормали. Пусть - боковая поверхность, - нижнее, - верхнее основание цилиндра; их площадь обозначим соответственно ||, ||=||=|?2|. Проинтегрируем уравнение (1.59) по объему цилиндра и устремим ||>0, а затем и ||>0. Тогда интеграл от С??Т??t в пределе обратится в нуль. Объемные интегралы
где
преобразуем в поверхностные по поперечному сечению цилиндра. Предположим, что для функции Ф(Т)=0 выполнено условие существования обратной функции dФ/dТ>0. Тогда ?(T-)=?(Ф) и ??(Ф)?t=?(Ф)(?Ф??t). Элемент объема dV=d?dn, где d? - элемент площади плоских участков, параллельных. Так как gradФ направлен по нормали к поверхности Ф(Т)=0, то интегрирование по нормали можно заменить интегрированием по Ф; при этом, принимая во внимание следующие соотношения
при ||>0, получим
где - среднее поперечное сечение цилиндра.
Для преобразования воспользуемся формулой Остроградского-Гаусса
при ||>0, интеграл по обратится в нуль, а интегралы по ||,|| перейдут в интеграл по
Подставим сюда выражение
и заменим единичный вектор нормали из равенства
Тогда для получим следующее выражение:
Вычислив значение подынтегрального выражения в точке Р, окончательно имеем
Из уравнения (1.59) следует=, т.е. получаем условие (1.56).
Таким образом, решение задачи типа Стефана (1.53)-(1.58) сводится к решению уравнения (1.59) с дополнительными условиями (1.55), (1.57), (1.58). Левая часть уравнения (1.59) содержит сосредоточенную теплоемкость L??(Т-) на поверхности фазового перехода Т=, т.е. она обращается в нуль при Т?. Теплофизический смысл этого члена заключается в том, что теплота фазового перехода L? выделяется на фазовом фронте. Если построить для этого уравнения разностное уравнение, то коэффициент левой части уравнения (1.59) будет вычисляться в узлах сетки, а меняющийся фронт фазового перехода не всегда совпадает узлом сетки. Отсюда следует, что разностная схема не всегда будет учитывать теплоту фазового перехода, т.е. она не будет обладать свойством консервативности, Возникновение такой ситуации приводит к необходимости сглаживания коэффициентов уравнения (1.59). Для этого дельта-функция приближенно заменяется дельтообразной, или размазанной, дельта-функцией ?(Т-,?)?0, где ? - величина полуинтервала, на котором отлична от нуля ?(Т-,?).
Таким образом, вводится сглаженная или эффективная теплоемкость
?)
удовлетворяющая следующим условиям:
Т<-?,
Т>+?
2. Изменение энтальпии на интервале (-?,+?) сохраняется
(+?)- (-?)=H(+?)-H(-?)
На этом же интервале (-?,+?) проводится сглаживание коэффициента теплопроводности ?.
В результате вместо задачи (1.59), (1.55), (1.57), (1.58) получается задача для уравнения теплопроводности со сглаженными коэффициентами
(1.60)
(1.61)
(1.62)
На практике сглаженную теплоемкость выбирают по разному.
Рассмотрим некоторые часто применяемые варианты сглаживания.
Пусть коэффициенты не зависят от Т.
1.Сглаженная теплоемкость на интервале сглаживания постоянна. Тогда из второго условия следует
.
2?+/2
и уравнение (1.60) решается с коэффициентом
. Сглаженная теплоемкость на интервале сглаживания есть линейная функция температуры, имеющая равные скачки на концах интервала сглаживания, т.е. строится в виде:
Подставив во второе условие, имеем
2b?=+(+)?
а равенство скачков дает
(-?)+= (+?)+
Таким образом, для этого варианта сглаживания уравнение (1.60) имеет следующий коэффициент:
Для решения задачи (1.60)-(1.62) теперь можно построить разностные схемы. Так как сглаженные коэффициенты зависят от температуры, получающаяся разностная задача будет нелинейной и ее решение будет найдено с использованием итерационного процесса.
2.4 О выборе параметра сглаживания
Разностная задача для дифференциальной задачи (1.60)-(1.62) будет вполне определенной, если указать способ выбора параметра сглаживания ?.
Коэффициенты разностного уравнения вычисляются в зависимости от интервала температуры по разным выражениям. Теплота фазового перехода, выделяющаяся на фронте, учитывается в выражении , определенном в интервале сглаживания. Если в указанном интервале не попадает ни одно узловое значение температуры, то температурное поле будет определено без учета выделения теплоты на фронте. Из сказанного следует условие выбора длины интервала сглаживания: параметр сглаживания ? должен быть выбран так, чтобы интервал сглаживания содержал в себе на каждом шаге по времени интервал, определяемый значениями температуры хотя бы в двух соседних узлах, между которыми находится фронт фазового перехода. Если фронт находится между узлами , то параметр ? должен удовлетворять условию
? ? ||
где значения т?/p>