Размещено на http://www.allbest.ru/

Расчет температуры воды по глубине водоема

Расчет температуры воды водоемов методом суперпозиции (наложения) предложен А.И. Пеховичем и В.М. Жидких. Этот метод изложен в работе [11] и рекомендациях по термическому расчету водохранилищ [36]. Метод предусматривает использование дифференциального уравнения теплопроводности для непроточного водоема (Лекция №7):

?t/?ф = a_т ?²t/?z², (5.1)

где a_т = л_т/(cс) - коэффициент турбулентной температуропроводности.

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

Этот принцип строго применим к системам, поведение которых описывается линейными соотношениями.

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

Разложение сложной тепловой задачи на простые должно производиться таким образом, чтобы сумма значений начальной температуры (t₀₁ + t₀₂ + …) и тепловых условий на поверхности воды (Q₁ + Q₂ + …) и на дне (Q_д1 + Q_д2 + …) для слагаемых задач была равна начальной температуре (t₀ = t₀₁ + t₀₂ + …) и тепловым условиям на поверхности (Q_п = Q₁ + Q₂ + …) и на дне (Q_д = Q_д1 + Q_д2 + …) в основной задаче. Коэффициенты температуропроводности a_т, теплопроводности л_т и теплопередачи б в решаемой и слагаемых задачах должны быть одинаковыми, за исключением случаев, в которых a_т и л_т меняются во времени.

Из изложенного следует, что для решения сложной тепловой задачи необходимо иметь набор решений простых задач. Авторы метода А.И. Пехович и В.М. Жидких разработали аналитические решения для 19 таких простых задач. Эти решения представлены в виде расчетных графиков в безразмерных координатах и сводной таблицы. Решения 19 задач позволяют рассчитать температуру воды в мелких, глубоких и очень глубоких водохранилищах, как при отсутствии ледяного покрова, так и при его наличии, а также в водохранилищах при их наполнении.

Безразмерные координаты графиков в зависимости от номера задачи (начальных и граничных условий) представлены искомой относительной избыточной температурой:

и_и1 = (t - t_п)/(t₀ - t_п); и_и2 = (t - и₂)/(t₀ - и₂);

и_и3 = (t - t₀)/(bф) и т.п., (5.2)

критерием Фурье

Fo = a_тф/h², (5.3)

критерием Био

Bi = бh/ л_т (5.4)

и относительной глубиной з = z/h, где t, t₀, t_п и и₂ - соответственно температура воды в точке, начальная и на поверхности, а также температура воздуха на высоте 2 м; b - коэффициент при линейном задании температуры поверхности воды или воздуха; a_т - коэффициент турбулентной температуропроводности; ф - время; z и h - соответственно переменная и полная глубина водохранилища; б и л_т - соответственно коэффициенты теплоотдачи и турбулентной теплопроводности.

Рассмотрим метод суперпозиции на примере решения конкретной тепловой задачи, заимствованной из рекомендаций [36].

Требуется найти распределение температуры воды по глубине на конец третьей декады июня в слабопроточном водохранилище глубиной 40 м, если в начальный момент (1 июня) температура воды по глубине одинакова и равна 4 °С. Нагрев воды происходит в результате теплообмена с атмосферой, его ход показан на рис. 5.1 (схема 1): в течение первой декады (ф₁) тепловой поток постоянен (Q₁ = 150 Вт/м²), в течение двух последующих декад он возрастает, причем во второй декаде (ф₂) со скоростью Q_о = 0,4 Вт/(м²·ч), а в третьей (ф₃) - со скоростью Q_о = 0,3 Вт/(м²·ч). Коэффициенты турбулентной тепло- и температуропроводности воды соответственно равны: л_т = 1000 Вт/(м·°С) и a_т = 1 м²/ч.

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

Рис. 5.1. Разложение теплообмена с атмосферой (1) на составляющие (2, 3, 4) [8]

1. Согласно принципу суперпозиции, раскладываем тепловой поток, приходящий на поверхность воды, на три составляющие (рис. 5.1, схемы 2, 3, 4). Первый поток Q₁ действует в течение всего расчетного периода ф = ф₁ + ф₂ + ф₃ = 30 сут = 720 ч. Второй поток действует с интенсивностью в течение периода ф₂ + ф₃ = 20 сут = 480 ч он равен Q₂ = (ф₂ + ф₃) =0,4 (ф₂ + ф₃) Вт/м². Третий поток теплоты действует в течение периода ф₃ = 10 сут = 240 ч. Так как действие второго потока интенсивностью мы распространили и на период ф₃, в то время как в этот период она равна, т.е. ниже, чем во второй декаде, поэтому третий поток следует находить по формуле Q₃ = () ф₃ = - 0,1ф₃ Вт/м² (рис. 5.1, схема 4).

Итак, решение общей задачи находим в виде суммы решений трех задач - по числу соответствующих потоков (Q₁, Q₂, Q₃).

2. Для каждой из трех задач устанавливаем начальные и граничные условия. В качестве начальных условий для первой задачи принимаем условия основной задачи: t₀ = 4 °C. Тогда во второй и третьей задачах, согласно условию разложения сложной задачи на простые, в качестве начальных условий следует принять t₀₂ = t₀₃ = 0 °С. В первой задаче в качестве граничного условия на поверхности воды принят источник Q₁ (теплообмен с атмосферой постоянный), во второй - Q₂ (теплообмен с атмосферой возрастает) и в третьей - Q₃ (теплообмен с атмосферой возрастает, но его рост ниже, чем во втором периоде).

Так как распределение температуры рассматривается в летний период (период отсутствия ледяного покрова), то для всех трех декад можно принять граничное условие на дне

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

3. Находим решение общей задачи в виде суммы решений трех задач. Для этого обращаемся к перечню решений 19 простых задач, разработанных А.И. Пеховичем и В.М. Жидких, и обнаруживаем, что первая задача совпадает с задачей №6, а вторая и третья задачи - с задачей №7 этого перечня (рис. 5.2). Причем во второй задаче в качестве Q₀ (графа 5) необходимо принять.

Расчетная формула для определения температуры воды с учетом полученных решений для трех простых задач имеет вид

(5.5)

Решения слагаемых (простых) задач

суперпозиция температуропроводность водохранилище непроточный

где относительная избыточная температура и_иi определяемая формулами (5.2), находится по графикам, построенным для каждой задачи, в зависимости от критерия Фурье и относительной глубины з = z/h.

Результаты расчета температуры воды в водоеме по его глубине для рассматриваемого примера приведены в табл. 5.1.

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

Расчет температуры воды по глубине водоема