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

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

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

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

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

T0(M)\displaystyle T_0(M) =Φ(M,0),\displaystyle {} = \Phi(M, 0),​M\displaystyle M​∈S.\displaystyle {} \in S.
(2.19)

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

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

∂Φ(M,0)∂t\displaystyle \frac{\partial \Phi(M, 0)}{\partial t} =a2⋅ΔT0(M),\displaystyle {} = a^2 \cdot \Delta T_0(M),​M\displaystyle M​∈S.\displaystyle {} \in S.
(2.20)

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

∂T0∂n→\displaystyle \frac{\partial T_0}{\partial \overrightarrow{n}} =Φ(M,0),\displaystyle {} = \Phi(M, 0),​M\displaystyle M​∈S.\displaystyle {} \in S.

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

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

a2⋅ΔT(M)\displaystyle a^2 \cdot \Delta T(M)​+f(M)\displaystyle {} + f(M)​=0.\displaystyle {} = 0.
(2.21)

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

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

∭G(U(M)⋅ΔT(M)−T(M)⋅ΔU(M)) dM\displaystyle \iiint_G \left( U(M) \cdot \Delta T(M) - T(M) \cdot \Delta U(M) \right) \,dM =∬S(U(M)⋅∂T(M)∂n→\displaystyle {} = \iint_S \biggl( U(M) \cdot \frac{\partial T(M)}{\partial \overrightarrow{n}}​−T(M)⋅∂U(M)∂n→) dS.\displaystyle {} - T(M) \cdot \frac{\partial U(M)}{\partial \overrightarrow{n}} \biggr) \,dS.

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

−∭GΔT(M) dM\displaystyle {} - \iiint_G \Delta T(M) \,dM =−∬S∂T(M)∂n→ dS.\displaystyle {} = - \iint_S \frac{\partial T(M)}{\partial \overrightarrow{n}} \,dS.

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

∭Gf(M) dM\displaystyle \iiint_G f(M) \,dM =−a2⋅∬SΦ(M) dS,\displaystyle {} = - a^2 \cdot \iint_S \Phi(M) \,dS,
(2.22)

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