论文部分内容阅读
采用微分求积法数值求解流函数-涡度方程来模拟二维驱动方腔流体时会遇到流函数超约束的问题：虽然流函数方程为二阶偏微分方程，但在每个边界上都存在两个约束条件，一个Dirichlet 条件和一个Neumann 条件。针对该问题进行深入分析，并提出一种新的处理方法：在边界涡的计算中考虑Neumann 条件，而仅将Dirichlet 条件施加于流函数方程。数值结果显示该方法是可行且高效的。同时给出由前人提出的单层法和双层法进行比较。试算表明单层法对于网格数的奇偶性很敏感，不适于处理本文问题。对比本文方法与双层法发现：本文方法的计算精度较高，且由于回避了超约束问题而更方便于使用。