[论文解读] A Numerical Approach to Solving Nonlinear Differential Equations on a Grid with Potential Applicability to Computational Fluid Dynamics

Jesper Tveit|arXiv (Cornell University)|Sep 3, 2014

Advanced Numerical Methods in Computational Mathematics参考文献 5被引用 4

一句话总结

本文提出了一种高阶有限元方法，通过在网格上同时引入函数值及其至指定阶数的导数，求解非线性微分方程，从而在减少网格分辨率的前提下实现高精度解。对于二维顶盖驱动腔流基准问题，该方法在多项式阶数最高达9时，实现了稳定且高雷诺数（最高达 Re = 4×10⁴）的稳态解，优于传统 h-型超细化方法，展现出更高的效率与精度。

ABSTRACT

A finite element method for solving nonlinear differential equations on a grid, with potential applicability to computational fluid dynamics (CFD), is developed and tested. The current method facilitates the computation of solutions of a high polynomial degree on a grid. A high polynomial degree is achieved by interpolating both the value, and the value of the derivatives up to a given order, of continuously distributed unknown variables. The two-dimensional lid-driven cavity, a common benchmark problem for CFD methods, is used as a test case. It is shown that increasing the polynomial degree has some advantages, compared to increasing the number of grid-points, when solving the given benchmark problem using the current method. The current method yields results which agree well with previously published results for this test case.

研究动机与目标

• 开发一种在结构化网格上求解非线性微分方程的数值方法，以实现高空间精度。
• 研究在求解具有挑战性的CFD基准问题时，p-型超细化（提高多项式阶数）相较于h-型超细化（提高网格分辨率）的优势。
• 计算二维顶盖驱动腔流中不可压Navier-Stokes方程的高雷诺数稳态解。
• 证明该方法能够产生稳定且精确的解，而不会出现伪时间方法中常见的虚假周期性伪影。
• 探索该方法推广至高维、非定常流或多相流问题的潜力，且结构改动最小。

提出的方法

• 该方法采用有限元方法，将每个网格点上的函数值及其至(Ω−1)阶导数显式存储，从而实现高阶连续性。
• 基于网格数据建立控制非线性PDE的弱形式，得到依赖于插值函数和导数值的方程组。
• 通过残差最小化方案进行迭代求解，每次迭代包含线性和非线性两个阶段以改进近似解。
• 采用基于矩阵的框架进行多维数据的索引与处理，包括使用四维矩阵表示网格上函数和导数值。
• 该离散化方法确保在整个区域具有C^(Ω−1)阶连续性，从而在每个单元内实现高阶多项式逼近。
• 该方法可通过提高连续性阶数（p-富集）实现局部超细化，而无需增加网格点数量。

实验结果

研究问题

• RQ1 在求解二维顶盖驱动腔问题时，提高多项式阶数（p-型超细化）是否比提高网格分辨率（h-型超细化）更具精度与效率优势？
• RQ2 该方法能否在雷诺数超过20,000时，实现不可压Navier-Stokes方程的稳定高阶稳态解？
• RQ3 在顶盖驱动腔基准问题中，该方法与现有高阶方法相比，在精度与计算成本方面表现如何？
• RQ4 数值不准确性在多大程度上会模拟出如湍流或周期性流动等物理不稳定性？
• RQ5 该方法能否在结构改动最小的前提下，推广至复杂几何形状、多相流或非定常问题？

主要发现

• 该方法成功计算了二维顶盖驱动腔流在雷诺数Re = 4×10⁴以内的稳态解，该雷诺数为该基准问题中已报道的最高值之一。
• 实现了最高达9次多项式（Ω = 5）的解，表明在高雷诺数流动中，p-型超细化在精度与效率方面可优于h-型超细化。
• 该方法结果与先前发表的基准结果高度一致，验证了其精度与可靠性。
• 使用高阶连续性（Ω = 5）在高雷诺数模拟中最为高效，尽管数值精度限制了进一步提高阶数。
• 该方法避免了伪时间积分方案中常见的虚假周期性解，表明所观察到的不稳定性为数值现象而非物理现象。
• 数值积分的计算成本随网格大小线性增长，而系统求解成本增长更迅速，提示在大规模应用中存在优化潜力。