空气动力学数值方法：有限体积法(FVM)：多物理场耦合问题的FVM方法.pdfVIP

空气动力学数值方法：有限体积法(FVM)：多物理场耦合问

题的FVM方法

1空气动力学数值方法：有限体积法(FVM)：多物理场耦合

问题的FVM方法

1.1绪论

1.1.1有限体积法的基本概念

有限体积法(FiniteVolumeMethod,FVM)是一种广泛应用于流体力学、热传

导、电磁学等领域的数值方法。它基于守恒定律，将计算域划分为一系列控制

体积，然后在每个控制体积上应用守恒方程，从而将偏微分方程转化为代数方

程组。这种方法的优势在于它能够直接处理守恒形式的方程，确保了质量、动

量和能量的守恒性，因此在处理复杂流体问题时特别有效。

1.1.2多物理场耦合问题的引入

在实际工程应用中，许多问题涉及多个物理场的相互作用，如流体流动与

热传导、电磁场与流体动力学的耦合等。这些多物理场耦合问题的求解通常比

单一物理场问题更为复杂，因为它们需要同时满足不同物理场的守恒条件，并

且物理场之间的影响需要被准确地模拟。有限体积法在处理这类问题时，通过

在每个控制体积上同时应用多个物理场的守恒方程，能够有效地耦合不同物理

场，提供一个统一的数值框架。

1.2有限体积法在多物理场耦合问题中的应用

在多物理场耦合问题中，有限体积法的应用通常涉及以下步骤：

1.网格划分：将计算域划分为一系列控制体积，每个控制体积代表

一个离散的计算单元。

2.守恒方程的离散化：在每个控制体积上应用守恒方程，将连续的

偏微分方程转化为离散的代数方程。

3.耦合条件的处理：确定不同物理场之间的耦合条件，如热传导与

流体流动之间的能量交换。

4.迭代求解：通过迭代方法求解离散后的代数方程组，直到满足收

敛条件。

5.后处理：分析和可视化求解结果，以评估物理场的耦合效果。

1

1.2.1示例：流体流动与热传导的耦合

假设我们有一个二维的流体流动与热传导耦合问题，其中流体流动由

Navier-Stokes方程描述，热传导由能量守恒方程描述。我们将使用有限体积法

来求解这个问题。

1.2.1.1网格划分

我们首先将计算域划分为一个由正方形控制体积组成的网格。每个控制体

积的大小为×。

1.2.1.2守恒方程的离散化

对于流体流动，我们有连续性方程和动量方程。在控制体积上，连续性

方程可以表示为：

​∂​

∂+⋅=0

∂

动量方程可以表示为：

​∂+⋅=−+​​

∂∂∂

对于热传导，能量守恒方程可以表示为：

​∂​​

∂+⋅=⋅

∂∂

其中，是密度，是速度向量，是压力，是体积力，是总能量，是热

流向量。

1.2.1.3耦合条件的处理

在流体流动与热传导的耦合问题中，耦合条件主要体现在能量方程中流体

速度的影响。这意味着在求解能量方程时，需要使用由Navier-Stokes方程求

得的速度场。

1.2.1.4迭代求解

使用有限体积法，我们将上述方程转化为代数方程组。然后，通过迭代方

法求解这些方程，直到满足收敛条件。迭代求解过程可以使用如SIMPLE算法或

其变体。

1.2.1.5后处理

最后，我们分析和可视化求解结果，以评估流体流动与热传导的耦合效果。