Recently Updated Pages
12.4 结构网格转非结构网格程序
12.4 结构网格转非结构网格程序 grid_unstr_2d 目录中的程序，将二维结构网格转换为非结构三角形网格。它连接距离较短的一对对角节点，将结构网格的每个四边形单元划分为两个三角形。所得...
12.3 二维结构网格生成器
12.3 二维结构网格生成器 grid_struct_2d 目录包含三个生成二维外流与内流结构网格的程序。第一个的源代码位于 cgrid 子目录，用于生成翼型周围 C 型网格，见11.1.1节及...
12.2 一维结构网格生成器
12.2 一维结构网格生成器 grid_1d 目录包含生成一维结构网格的程序。每个网格节点对应一个截面积，沿 x 轴的面积分布对应 Laval 喷管，按以下关系计算： $$ \begin{gat...
12.1 稳定性分析程序
第 12 章 软件应用 软件包包含若干一维、二维流动求解器及网格生成器的源代码和可执行程序，并提供多种一维、二维流动算例的输入数据与网格。此外，还包含两个用于显式、隐式时间推进格式 von Ne...
11 参考文献
参考文献 [1] Brodersen O, Hepperle M, Ronzheimer A, Rossow CC, Schöning B. The parametric grid genera...
11.2 非结构网格
11.2 非结构网格 非结构网格通常由二维三角形或三维四面体组成。但如今，用多种单元类型构造非结构网格越来越常见。例如，以六面体或棱柱体离散边界层，流动区域其余部分用四面体填充，并在六面体或棱柱...
11.1 结构网格
第 11 章 网格生成原理 在数值求解控制方程之前，必须离散所有边界表面，并在流动区域内部生成体网格。如第3章开头的3.1节所述，基本上可选择： 结构网格； 非结构网格。 图11.1...
10 参考文献
参考文献 [1] Roache PJ. Quantification of uncertainty in computational fluid dynamics. Ann Rev Fluid ...
10.3 von Neumann 稳定性分析
10.3 von Neumann 稳定性分析 在实现新的离散方法之前，至少应大致了解它对数值格式稳定性与收敛行为的影响。已有大量实践证实，von Neumann 稳定性分析能够可靠地评估求解格式...
10.2 离散格式的精度
10.2 离散格式的精度 离散格式的精度与其截断误差有关。例如，若截断误差的主项与 $ \Delta x $ 成正比，则称为空间一阶精度格式；若主项随 $ (\Delta x)^{2} $ 变化...
10.1 相容性要求
第 10 章 相容性、精度与稳定性 在有限体积方法中，需要采用适当方法，近似计算每个控制体上 Navier–Stokes 方程（2.19）的面积分。第4章和第5章介绍了多种空间离散方案。由于对穿...
9 参考文献
参考文献 [1] Kroll N, Jain RK. Solution of two-dimensional Euler equations—experience with a finite v...
9.6 并行化
9.6 并行化 如今，几乎所有计算机——台式机、笔记本电脑和平板电脑——的中央处理器（CPU）都由多个核心组成。此外，现代图形处理器（GPU）还提供数百乃至数千个可进行通用编程的处理单元。然而，...
9.5 低马赫数预处理
9.5 低马赫数预处理 在低亚声速马赫数范围，当流速远小于声速时，控制方程（2.19）的对流项出现刚性。以下示例可以说明这一点。三维情况下，五个特征值为 $$ \begin{aligned}(\...
9.4 多重网格
9.4 多重网格 多重网格是一种强有力的加速技术（参见图 3.7）。它在一系列逐渐变粗的网格上求解控制方程，然后将粗网格的解更新量合并并加到最细网格的解上。Brandt $ [19] $ 最初为...
9.3 残差光顺
9.3 残差光顺 通过优化级系数 $ [5-7] $，可在一定程度上改变显式多级时间推进格式（6.1.1、6.1.2 节）的最大 CFL 数和收敛性质。Jameson 和 Baker $ [8]...
9.2 焓阻尼
9.2 焓阻尼 某些情况下，总焓 H（式（2.12））在整个流场中保持常数，例如由 Euler 方程描述、没有热源及外力的外部流动。可以利用这一特点降低计算成本。 图 9.1 非结构网格无黏...
9.1 局部时间推进
第 9 章 加速技术 为加速定常问题中控制方程（2.19）的求解，人们发展了多种方法。这些加速技术也适用于双时间推进格式（6.3 节）的内迭代。本章讨论以下数值方法： 局部时间推进； 焓...
8 参考文献
参考文献 [1] Kroll N, Jain RK. Solution of two-dimensional Euler equations—experience with a finite v...
8.10 非结构网格边界处的流动梯度
8.10 非结构网格边界处的流动梯度 5.3.4 节已指出，中位对偶格式计算流动梯度时需要谨慎处理。如果在三角形或四面体网格上采用式（5.50）的 Green–Gauss 方法计算梯度，区域边界...