Глава 1. Решение как процесс на пространственно-временной сетке

Где стоит глава

Открытие Части IX

Часть VIII решала обыкновенные дифференциальные уравнения — эволюцию величины во времени, где состояние есть одно рациональное число $y_n$. Часть IX делает следующий шаг: уравнения в частных производных, описывающие эволюцию функции, зависящей и от времени, и от пространства. Это язык всей континуальной физики — теплопроводности, волн, переноса, течений жидкости.

Классическая теория берёт решение уравнения в частных производных как завершённую функцию $u(x,t)$ на области, элемент пространства $C^2$ или Соболева. ToS применяет здесь приём, доведённый в Части VIII до гарантий сходимости. Решение для нас — процесс: но теперь состояние на шаге $n$ есть уже не число, а grid-функция $u^n$, сопоставляющая каждому узлу пространственной сетки рациональное значение (мгновенный профиль). Эволюция задаётся конечно-разностной схемой

$$u_j^{n+1}=u_j^n+\Delta t\cdot (\text{дискретный пространственный оператор}{)}_j,$$

а пространственная производная замещается конечной разностью. Где в Части VIII была лестница чисел $(1+1/N{)}^N$, в Части IX стоит лестница grid-функций; метод Эйлера во времени остаётся, добавляется дискретизация пространства.

Настоящая глава — каркас всей Части. Она вводит grid-функции и конечные разности, дискретный аналог интегрирования по частям (формулу Грина) и — что задаёт структуру глав 9.2–9.4 — три типа уравнений как три роли пространственного оператора.

Что эта глава утверждает — и что нет

Доказанное ядро — конечная дискретная математика над ℚ. Grid-функция есть отображение $N\to Q$; прямая разность — дискретный градиент, вторая разность — дискретный лапласиан; суммирование по частям (формула Абеля) и дискретная формула Грина связывают энергию с градиентом; норма лапласиана ограничена нормой градиента.¹ Всё это — конечные рациональные тождества и неравенства.

Чего глава не делает: она не строит завершённую функцию-решение $u(x,t)$ как объект, не доказывает непрерывный предел (измельчение сетки до нуля) и не вводит завершённые пространства Соболева. Эти пределы — роль-пределы; их процессные представители и конечные энергетические аналоги (дискретная энергия, $L^2$-ядра Частей VI–VII) с P4 совместимы, но как завершённые объекты они остаются границей.

Три уровня строгости, как и в Части VIII. Первый — конечная дискретная математика над ℚ: grid-функции, конечные разности, суммирование по частям, формула Грина — доказаны без аксиом. Второй — процессные чтения: решение как процесс grid-функций, пространственная производная как разность, законы сохранения как конечные суммы, устойчивость как условие Куранта. Третий — завершённое непрерывное решение, непрерывный предел и Соболевы пространства — роль-пределы и направление. Первые два глава ведёт как доказанные, третий — как честно помеченную границу.

E/R/R-каркас главы

Каркас в порождающем порядке Rules $\to$ Roles $\to$ Elements; полный разбор — в § «Разбор в порождающем порядке». Глава вводит систему решения уравнения в частных производных как процесса на пространственно-временной сетке.

Rules (Правила, $\leftrightarrow$ L5). Пространственно-временная сетка ($x_j=j\Delta x$, $t_n=n\Delta t$); конечно-разностная схема $u_j^{n+1}=u_j^n+\Delta t\cdot ({\mathcal{L}}_h{u}^n{)}_j$, где ${\mathcal{L}}_h$ — дискретный пространственный оператор; конечные законы (энергия, масса, принцип максимума) как конечные суммы по узлам; устойчивость как условие Куранта (CFL).

Roles (Роли, $\leftrightarrow$ L4). Профиль $u^n$ — роль-состояние (моментальный снимок во времени); дискретный оператор ${\mathcal{L}}_h$ — роль-пространственный-закон (как профиль меняется по пространству); число Куранта — роль-условие устойчивости; завершённое $u(x,t)$ — роль-предел.

Elements (Элементы, $\leftrightarrow$ L1$+$P4). Grid-функции $u^n:\{0,\dots ,N-1\}\to Q$, рациональные значения в узлах $u_j^n$, конечная пространственная сетка, конечный номер шага. Под P4 завершённая функция-решение и непрерывный предел — вне элементов: актуальны конечные grid-функции на конечных стадиях.

От ОДУ к УЧП: состояние становится функцией

Обыкновенное уравнение $y^{\prime }=f(t,y)$ описывает эволюцию числа: состояние $y$ — одна величина, метод Эйлера двигает её по времени, $y_{n+1}=y_n+hf(t_n,y_n)$. Уравнение в частных производных ${\partial }_{t}u=\mathcal{L}u$ описывает эволюцию функции: состояние $u(\cdot ,t)$ — целый профиль по пространству, а $\mathcal{L}$ — пространственный дифференциальный оператор (первая производная ${\partial }_x$, лапласиан ${\partial }_x^2$ и т. п.).

Процессный приём один и тот же на обоих уровнях, с одним добавлением: дискретизуется не только время, но и пространство. Разобьём пространство сеткой узлов $x_j=j\Delta x$; состояние $u^n$ становится grid-функцией (значениями в узлах); оператор $\mathcal{L}$ замещается дискретным оператором ${\mathcal{L}}_h$ (конечной разностью); и схема

$$u^{n+1}=u^n+\Delta t\cdot {\mathcal{L}}_h{u}^n$$

есть в точности метод Эйлера во времени, но на пространстве grid-функций. Весь аппарат Части VIII — шаг Эйлера, сжатие, устойчивость, Гронуолл — переносится сюда; новизна лишь в том, что величина, которую двигает время, сама есть функция пространства.

Профиль $u^n$ — роль-состояние: моментальный снимок поля во времени. Решение уравнения в частных производных — роль-предел процесса таких профилей $u^0,u^1,u^2,\dots$, а не завершённая функция $u(x,t)$. Актуальна стадия-профиль; непрерывное решение — направление.

Конечные разности: производная как разность

Пространственная производная на сетке становится конечной разностью. Прямая разность — дискретный градиент:

$$({\nabla }_{h}u{)}_j=u_{j+1}-u_j,$$

а её повторное применение — дискретный лапласиан (вторая разность):

$$({\Delta }_{h}u{)}_j=u_{j+2}-2u_{j+1}+u_j=({\nabla }_h{\nabla }_{h}u{)}_j.$$

Это конечные рациональные операции над grid-функцией, без предельного перехода.² Техническая оговорка: эта вторая разность — повторная прямая ($u_{j+2}-2u_{j+1}+u_j$); для классической схемы теплопроводности (Глава 9.2) удобнее центрированная запись $u_{j+1}-2u_j+u_{j-1}$, вводимая там локально с явной оговоркой об индексации.

Ключевой инструмент анализа — дискретный аналог интегрирования по частям. На конечной сетке роль интеграла играет конечная сумма ${\sum }_j{u}_j$, и для неё выполнено суммирование по частям (формула Абеля): спаривание разности одной функции с другой переносится на вторую с граничными членами, а на периодической сетке граничные члены сокращаются.³ Отсюда — дискретная формула Грина, энергетическое тождество, связывающее лапласиан с градиентом:

$$\sum _j{u}_{j+1}({\Delta }_{h}u{)}_j=-\text{norm}{{\nabla }_{h}u}^2\text{(периодическая сетка),}$$

и сопутствующая оценка $\text{norm}{{\Delta }_{h}u}^2\le 4\text{norm}{{\nabla }_{h}u}^2$.⁴ Эти тождества — конечные суммы над ℚ; непрерывная формула Грина $\int u\Delta u=-\int |\nabla u{|}^2$ — их роль-предел.

Конечная разность — роль-производная: место, которое производная занимает на дискретной сетке. Формула Грина — роль-инструмент: правило, переводящее пространственный оператор в энергетическую величину. Они конечны и рациональны; непрерывные производная и интеграл — их роль-пределы.

Три типа уравнений — три роли

Уравнения в частных производных, которые разбирает Часть IX, различаются тем, какую роль играет пространственный оператор в процессе. Три канонических типа задают систематический триптих — три следующие главы.

Параболическое — теплопроводность (${\partial }_{t}u={\partial }_x^{2}u$, Глава 9.2). Пространственный оператор — лапласиан; процесс сглаживает: каждый шаг заменяет значение средневзвешенным с соседями, профиль растекается. Роль оператора — диссипация: энергия убывает, действует принцип максимума (значения не выходят за начальные пределы).

Гиперболическое — волна (${\partial }_t^{2}u=c^2{\partial }_x^{2}u$, Глава 9.3). Уравнение второго порядка по времени; процесс распространяет: возмущение бежит по сетке. Роль оператора — конечная скорость: причинность, фронт продвигается не более чем на один узел за шаг.

Транспортное — перенос (${\partial }_{t}u+a{\partial }_{x}u=0$, Глава 9.4). Пространственный оператор — первая производная; процесс движет профиль без изменения формы. Роль оператора — сохранение: полная масса не меняется, профиль лишь сдвигается.

Все три — процессы grid-функций на одной пространственно-временной сетке; различие — в роли пространственного оператора (сглаживание, распространение, перенос). Кульминация — уравнения Навье–Стокса (Глава 9.5), нелинейная система, где те же энергетические методы дают диссипацию и условную регулярность.

Тип уравнения в частных производных — это роль пространственного оператора в процессе: параболическое сглаживает (диссипация), гиперболическое распространяет (конечная скорость), транспортное переносит (сохранение). Один процессный каркас, три роли — и это структура всей Части.

E/R/R-разбор: решение как Роль-процесс на сетке

Постановка

Разберём построенную систему — решение уравнения в частных производных как процесс на пространственно-временной сетке — в терминах E/R/R. Оговорка об уровне: <<система решения>> здесь — содержательная интерпретация конструкции (сетка плюс конечно-разностная схема плюс дискретные законы), а не объект, буквально объявленный в коде как System. Разбор есть онтологическое осмысление структуры Части и чтение авторских E/R/R-шапок опорных файлов,⁵ а не приписывание коду новых формальных утверждений. Как и прежде, аббревиатура задаёт эпистемический порядок (Elements $\to$ Roles $\to$ Rules); разбор ведём в порождающем.

Разбор в порождающем порядке

Rules (Правила, $\leftrightarrow$ L5). В основании — пространственно-временная сетка и конечно-разностная схема $u^{n+1}=u^n+\Delta t\cdot {\mathcal{L}}_h{u}^n$, где ${\mathcal{L}}_h$ — дискретный пространственный оператор (градиент, лапласиан). Над ней — дискретное исчисление: суммирование по частям и формула Грина, переводящие оператор в энергетические величины. Венчают конструкцию конечные законы — диссипация энергии, сохранение массы, принцип максимума — и условие устойчивости (CFL). Универсальный слой — законы L1–L5; конкретный — формулы разностей и схема.

Roles (Роли, $\leftrightarrow$ L4). Правила задают роли. Профиль $u^n$ — роль-состояние (снимок поля во времени). Дискретный оператор ${\mathcal{L}}_h$ — роль-пространственный-закон, и его конкретная форма решает тип: сглаживание, распространение или перенос. Число Куранта — роль-условие устойчивости. Дискретная энергия и масса — роль-сохраняющиеся (или роль-диссипирующие) величины. А завершённое $u(x,t)$ и непрерывный предел — роль-предел.

Elements (Элементы, $\leftrightarrow$ L1$+$P4). Носители конечны: grid-функции $u^n:\{0,\dots ,N-1\}\to Q$, значения в узлах $u_j^n\in Q$, конечные суммы ${\sum }_{j<N}$, конечная сетка и конечный номер шага. Под P4 ни завершённая функция-решение, ни непрерывный предел не суть элементы: актуальны конечные профили на конечных стадиях; актуализируется процесс измельчения, а не его предел.

Компонент	Что фиксирует	E/R/R
сетка $x_j,t_n$; схема $u^{n+1}=u^n+\Delta t{\mathcal{L}}_h{u}^n$; формула Грина; законы; CFL	КАК профиль эволюционирует и контролируется	Rules ($L5$)
профиль $u^n$; оператор ${\mathcal{L}}_h$ (тип); число Куранта; энергия/масса; $u(x,t)$ как предел	ЗАЧЕМ значимы носители: роли в процессе	Roles ($L4$)
grid-функции $u_j^n\in Q$; конечные суммы; сетка; номер шага	ЧТО есть на каждой стадии (конечно, P4)	Elements ($L1$, P4)

Проверка сформированности и что даёт разбор

{ Система сформирована корректно: каждый компонент попадает ровно в одну E/R/R-категорию, и нет самоотнесения (P1). Существенно, чего здесь нет: процесс решения не берёт своим элементом <<функцию $u(x,t)$ как завершённый объект>> или <<непрерывную область>>. Пространство дискретизовано конечной сеткой, производная есть конечная разность, а законы (энергия, масса, максимум) проверяются как конечные суммы. Именно эта конечность каждой стадии отделяет доказанное (дискретное исчисление над ℚ, без аксиом; законы и устойчивость конкретных схем — в следующих главах, тоже как конечные рациональные факты) от границы (непрерывное решение и предел — роль-предел).}

Что даёт разбор. Он переводит замысел Части на язык структуры: решение — это Роль-процесс grid-функций; тип уравнения — это Роль пространственного оператора; конечная разность — Роль-производная; а сами профили $u_j^n$ и суммы — Элементы (конечные на каждой стадии). Брать <<решение>> или непрерывную функцию как завершённый объект-Элемент значит ставить роль-предел на место Элемента — то же смешение категорий, что <<решение ОДУ как готовая функция>> (Глава 8.1) или <<точка как объект>> вместо класса (Глава 4.3). Уравнение в частных производных решается не извлечением готовой функции, а разворачиванием процесса профилей на сетке.

Следующая глава предъявляет первый и канонический процесс этого рода — уравнение теплопроводности: диффузию как процесс сглаживания, с доказанным принципом максимума и диссипацией энергии над ℚ.

---

---

Часть: Часть IX. Уравнения в частных производных · Том: «Математика»

Навигация: ← Глава 5. Экспонента, линейные уравнения, применения — Часть VIII · Глава 2. Уравнение теплопроводности: диффузия как сглаживание →

Footnotes

1. Опорная инфраструктура — GridFunction.v и FiniteDifference.v (каталог navier_stokes/ Rocq-репозитория ToS; grid-функции, суммы, нормы, fd, dd, суммирование по частям abel_summation; $\sim 45$ доказанных утверждений, 0 аксиом), плюс новый файл ProcessPDEGreen.v (process/: дискретная формула Грина и оценка $\text{norm}{\Delta u}^2\le 4\text{norm}{\nabla u}^2$; 8 утверждений, 0 аксиом, Print Assumptions: <>). ↩

2. fd (прямая разность) и dd (вторая разность) в FiniteDifference.v; тождество dd_eq_fd_fd: вторая разность есть композиция первых. Линейность обоих операторов — там же, 0 аксиом. ↩

3. abel_summation и abel_periodic в FiniteDifference.v: $\sum ({\nabla }_{h}f{)}_i{g}_i+\sum f_{i+1}({\nabla }_{h}g{)}_i=f_N{g}_N-f_0{g}_0$, и для периодических — правая часть нуль. 0 аксиом. ↩

4. green_dirichlet_periodic и lap_norm_bound в ProcessPDEGreen.v (новый файл Части IX, проверено Print Assumptions: <>), 0 аксиом — для прямого лапласиана. Энергетическую оценку теплопроводности (Глава 9.2) даёт не это тождество напрямую, а выпуклость (неравенство Йенсена) вместе с леммами о циклическом сдвиге суммы (sum_ns_shift_periodic, sum_pred_shift_periodic, там же); формула Грина — канонический дискретный инструмент сам по себе. ↩

5. Шапки размечают прямо: ProcessPDEGreen.v — Elements: grid-функции, конечные суммы; Roles: fd роль-градиент, dd роль-лапласиан, формула Грина роль-инструмент; Rules: суммирование по частям, $\text{norm}{\Delta u}^2\le 4\text{norm}{\nabla u}^2$. ↩