Ошибки постановки краевой задачи

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

Согласование начальных и граничных условий

Первая типичная ошибка — задать начальное и граничные условия, которые не согласуются друг с другом на стыке. В нестационарной задаче решение ищется при t0t \ge 0 и MGM \in G; на границе SS при t=0t = 0 — начальное и граничное условия описывают одну и ту же точку. Если они дают там разные значения, в этой точке возникает разрыв.

Рассмотрим задачу Дирихле с начальным условием T(M,0)=T0(M)T(M, 0) = T_0(M) и граничным T(M,t)=Φ(M,t)T(M, t) = \Phi(M, t), MSM \in S. Условие согласования нулевого порядка требует совпадения данных

T0(M)=Φ(M,0),MS.T_0(M) = \Phi(M, 0), \quad M \in S.
(2.19)

Если условие (2.19) нарушено, ряд Фурье решения подходит к границе с разными пределами по tt и по SS: сходимость неравномерна, появляется явление Гиббса, и классического (гладкого) решения не существует.

Но и выполнения (2.19) недостаточно. Само уравнение Tt=a2ΔT\frac{\partial T}{\partial t} = a^2 \cdot \Delta T связывает производные данных: продифференцировав граничное условие по времени и взяв t=0t = 0, получаем условие согласования первого порядка

Φ(M,0)t=a2ΔT0(M),MS.\frac{\partial \Phi(M, 0)}{\partial t} = a^2 \cdot \Delta T_0(M), \quad M \in S.
(2.20)

Аналогично строятся условия более высоких порядков; для решения класса CkC^k нужно выполнение условий согласования примерно до порядка kk. Их нарушение не делает задачу неразрешимой, но ограничивает гладкость решения — об этом легко забыть, формально выписав «любые» T0(M)T_0(M) и Φ(M,t)\Phi(M, t). Для задачи Неймана условия согласования те же, но записываются для нормальной производной:

T0n=Φ(M,0),MS.\frac{\partial T_0}{\partial \vec{n}} = \Phi(M, 0), \quad M \in S.

Условие разрешимости стационарной задачи Неймана

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

a2ΔT(M)+f(M)=0,a^2 \cdot \Delta T(M) + f(M) = 0,
(2.21)

Краевые условия Неймана для него задаются граничным условием (1.14).

Известна вторая формула Грина

G(U(M)ΔT(M)T(M)ΔU(M))dM=S(U(M)T(M)nT(M)U(M)n)dS.\iiint_G \left( U(M) \cdot \Delta T(M) - T(M) \cdot \Delta U(M) \right) \,dM = \iint_S \left( U(M) \cdot \frac{\partial T(M)}{\partial \vec{n}} - T(M) \cdot \frac{\partial U(M)}{\partial \vec{n}} \right) \,dS.

Если принять U(M)=1U(M) = -1, то получим

GΔT(M)dM=ST(M)ndS.- \iiint_G \Delta T(M) \,dM = - \iint_S \frac{\partial T(M)}{\partial \vec{n}} \,dS.

Теперь подставим сюда уравнение (2.21) и граничные условия (1.14)

Gf(M)dM=a2SΦ(M)dS,\iiint_G f(M) \,dM = - a^2 \cdot \iint_S \Phi(M) \,dS,
(2.22)

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