文章详细介绍了基于连续性方程和动量方程的流体动力学求解算法，包括SIMPLE、PISO和它们的变种PLMPL。SIMPLE算法通过离散线性化动量方程推导出压力方程，而PISO则在SIMPLE基础上增加了对速度场的修正。PIMPLE是两者的结合，旨在每个时间步长中找到稳定的解。这些算法在解决流体力学问题时，通常涉及速度场和压力场的迭代求解过程。
SIMPLE & PISO &PLMPLE
参考链接
不可压缩连续性方程和动量方程为：
∇ ⋅ u = 0 ∂ u ∂ t + ∇ ⋅ ( u u ) − ∇ ⋅ ( ν ∇ u ) = − ∇ p \begin{aligned} &\nabla \cdot \mathbf{u} = 0 \\ &\frac{\partial \mathbf{u}}{\partial t} +\nabla \cdot(\mathbf{u}\mathbf{u})-\nabla\cdot(\nu \nabla \mathbf{u}) = -\nabla p \end{aligned} ∇⋅u=0∂t∂u+∇⋅(uu)−∇⋅(ν∇u)=−∇p
下面利用连续性方程与动量方程推到压力方程。
离散线性化动量方程，同时压力梯度的形式保持不变。得到：
[ a 1 , 1 a 1 , 2 a 1 , 3 ⋯ a 1 , n a 2 , 1 a 2 , 2 a 2 , 3 ⋯ a 2 , n a 3 , 1 a 3 , 2 a 3 , 3 ⋯ a 3 , n ⋮ ⋮ ⋮ ⋱ ⋮ a n , 1 a n , 2 a n , 3 ⋯ a n , n ] [ u 1 u 2 u 3 ⋮ u n ] = [ r 1 r 2 r 3 ⋮ r n ] + [ ( ∂ p / ∂ x ) 1 ( ∂ p / ∂ x ) 2 ( ∂ p / ∂ x ) 3 ⋮ ( ∂ p / ∂ x ) n ] \left[ \begin{matrix} a_{1,1} & a_{1,2} & a_{1,3} & \cdots & a_{1,n} \\ a_{2,1} & a_{2,2} & a_{2,3} & \cdots & a_{2,n} \\ a_{3,1} & a_{3,2} & a_{3,3} & \cdots & a_{3,n} \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ a_{n,1} & a_{n,2} & a_{n,3} & \cdots & a_{n,n} \\ \end{matrix} \right] \left[\begin{matrix} u_1 \\ u_2 \\ u_3 \\ \vdots \\ u_n \end{matrix}\right]= \left[\begin{matrix} r_1 \\ r_2 \\ r_3 \\ \vdots \\ r_n \end{matrix}\right]+ \left[\begin{matrix} (\partial p/\partial x)_{1} \\ (\partial p/\partial x)_{2} \\ (\partial p/\partial x)_{3} \\ \vdots \\ (\partial p/\partial x)_{n} \end{matrix}\right] a1,1a2,1a3,1⋮an,1a1,2a2,2a3,2⋮an,2a1,3a2,3a3,3⋮an,3⋯⋯⋯⋱⋯a1,na2,na3,n⋮an,n
扫一扫
7075
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
