（CFD）投影法求解二维不可压缩N-S方程.pdfVIP

（CFD）投影法求解⼆维不可压缩N-S⽅程 说明：程序旨在熟悉VOF求解流程，投影法求解Matlab编程主要参考Mark Owkes的Guide To CFD资料，对于VOF模型的对流项需要采 取特殊的离散⽅式，如：施主-受主格式，但本⽂中仅采取中⼼差分格式进⾏离散，会带来较⼤误差，其次对于每个时间步都对体积分数值 做矫正，限制其在0-1之间。后续如果使⽤了更精确的格式，再来更新。 摘要 本⽂介绍基于matlab平台的VOF编程实例，采⽤投影法求解⼆维不可压缩N-S⽅程，交错⽹格结合中⼼差分法离散控制⽅程。并对 Re=100⽅腔流进⾏计算验证，证明投影法求解程序的正确性，随后对⼆维溃坝算例进⾏数值计算。 关键词：数值模拟，VOF，Matlab，溃坝 第⼀章 数值离散 1.1投影法求解原理 投影法是求解不可压缩N-S⽅程的常⽤⽅法之⼀，对于⼆维不可压缩N-S⽅程以及连续性⽅程如公式(1-1)展开形式如公式(1-2) 2 ∂t + u⋅ ∇u= − ρ∇p+ υ∇ u { ∇ ⋅ u= 0 (1 - 1) 展开形式 + u + v = − ∇p+ υ( + ) 2 ∂v 1 ( 2 ∂ v) ∂v ∂v ∂ v + u + v = − ∇p+ υ ∂x + 2 ∂t ∂x ∂y ρ 2 ∂y ∂u ∂v ∂x + ∂y = 0 (1 - 2) 投影法原理可以分为以下三步，将时间推进分成三个⼦步，在中间步解出压⼒，⾸先需要确定时间离散格式，这⾥时间离散化采⽤显式欧拉 格式，则离散⽅程如（1- 3）所⽰。 n+ 1 n 1 u −u n+1 n n 2 n Δt = − ρ∇p − u ⋅ ∇u + υ∇ u (1 - 3) 将⽅程（1-3）分解为个⽅程（1-4）与 （1-5） 预算步 Step1： ∗ n u −u n n 2 n