Python火焰锋动力学和浅水表面波浪偏微分方程

🎯要点

📜微分方程 | 本文 - 用例

📜数学模型：Python流感常微分方程房室数学模型

📜图计算和算法：Python流感传播感染康复图模型计算和算法

📜水文模型：Python流体数据统计模型和浅水渗流平流模型模拟

在这里插入图片描述

🍇Python三个维度扩散方程

偏微分方程是具有多个独立变量、依赖于这些变量的未知函数以及未知函数关于独立变量的偏导数的方程。

通常样式为：
$\frac{\partial^2 u}{\partial x^2}+ B \frac{\partial^2 u}{\partial x \partial y}+ C \frac{\partial^2 u}{\partial y^2}+ D \frac{\partial u}{\partial x}+ E \frac{\partial u}{\partial y}+ F u= G$
通过仅考虑前三个系数 A B 和 C，我们可以确定我们正在处理什么方程或我们正在解决什么问题。

如果 $B^2-4 A C<0$ ，我们有一个椭圆偏微分方程，为拉普拉斯方程：
$\frac{\partial^2 u}{\partial x^2}+\frac{\partial^2 u}{\partial y^2}=0$
该方程的两个导数是空间 $x^2$ 和 $y^2$ 的导数，没有时间导数。

如果 $B^2-4 A C>0$ ，则我们有双曲偏微分方程，为波动方程：
$\frac{\partial^2 u}{\partial t^2}=\frac{\partial^2 u}{\partial y^2}$
如果 $B^2-4 A C=0$ ，则我们有抛物线偏微分方程，为扩散方程：
$\frac{\partial u}{\partial t}=\frac{\partial^2 u}{\partial y^2}$
尽管计算机非常擅长数学，但它们不懂微分方程。为了告诉计算机求解微分方程，我们需要对方程进行离散化。即：
$U = u (x, y)$
以 $x$ 表示的部分展开式为
$u(x+\Delta x, y)=u(x, y)+\Delta x \frac{\partial u}{\partial x}+\frac{(\Delta x)^2}{2} \frac{\partial^2 u}{\partial x^2}+\ldots$
通过丢弃二阶项，我们最终将得到一个非常近似的公式：
$\frac{\partial u}{\partial x} \approx \frac{u(x+\Delta x, y)-u(x, y)}{\Delta x}$
这是一个著名的欧拉方法，在常微分方程中可以看到。在某些有限差分中，它被称为前向差分。

📜流体力学电磁学运动学动力学化学和电路中欧拉法示例：Python重力弹弓流体晃动微分方程模型和交直流电阻电容电路

看一下下面的例子：
$\begin{aligned} & h=\Delta x \\ & u_{i, j}=u\left(x_i, y_j\right) \\ & u_{i+1, j}=u\left(x_i+h, y_j\right) \end{aligned}$
正向差分是：
$\left(\frac{\partial u}{\partial x}\right)=\frac{1}{h}\left[u_{i+1, j}-u_{i, j}\right]+O(h)$
通过查看正向差分公式，我们注意到正向差分是通过当前步骤 $u _{ i , j }$ 减去下一步 $u _{ i +1, j}$ 得出的。直观上，我们可以想到反向差分。
$\left(\frac{\partial u}{\partial x}\right)=\frac{1}{h}\left[u_{i, j}-u_{i-1, j}\right]+O(h)$
中间差分：
$u(x+\Delta x, y)-u(x-\Delta x, y)=2 \Delta x \frac{\partial u}{\partial x}+\frac{(\Delta x)^3}{2} \frac{\partial^3 u}{\partial x^3}+\ldots$

$\frac{\partial u}{\partial x} \approx \frac{u(x+\Delta x, y)-u(x-\Delta x, y)}{2 \Delta x}$

最终得到公式：
$\left(\frac{\partial u}{\partial x}\right)_{i, j}=\frac{1}{2 h}\left[u_{i+1, j}-u_{i-1, j}\right]+O\left(h^2\right)$
一维扩散计算代码：

import numpy as np
from matplotlib import pyplot as plt

L = 100

x = np.linspace(0, 1, L)
u = np.zeros(L)


u[L//2] = 1  
u[0] = 0
u[L-1] = 0

for t in range(100):
    for i in range(1, L-1):
        u[i] += (u[i+1] - 2*u[i] + u[i-1])/4

fig, ax = plt.subplots(figsize=(20,10))
ax.axis('off')
ax.scatter(x,u, linewidth=15, c=u, cmap='jet')
plt.show()

二维扩散计算代码：

import numpy as np
import matplotlib.pyplot as plt
from matplotlib import animation
import seaborn as sns
sns.set()


T = 400
u = np.zeros((100,100))
u[30,20] = 1

fig, ax = plt.subplots(figsize=(20,10))
ax.axis('off')
plot = ax.contourf(u, cmap='jet')
def ans(f):
    global u, plot
    
    for j in range(100):
        for i in range(1,99):
            u[i,j] += (u[i+1,j] + u[i,j+1] + u[i-1,j] + u[i,j-1] - 4*u[i,j])/4
    for c in plot.collections:
        c.remove()
    plot = ax.contourf(u, cmap='jet')
    return plot

anim = animation.FuncAnimation(fig, ans, frames=T)
plt.show()