2025年数值流体力学题库及答案.docxVIP

2025年数值流体力学题库及答案

1.简述有限体积法（FVM）离散三维N-S方程的基本步骤，并说明其如何保证数值解的守恒性。

有限体积法离散三维N-S方程的核心步骤包括：首先对计算域进行网格划分，将连续区域离散为有限个互不重叠的控制体积（CV）；其次对控制体积应用积分形式的守恒定律（如质量守恒、动量守恒），将控制方程在时间和空间上积分；接着通过对控制体积界面上的通量（对流项、扩散项）进行近似（如采用中心差分、迎风格式或高阶格式），将积分方程转化为代数方程；最后联立所有控制体积的离散方程，形成封闭的线性或非线性方程组并求解。

守恒性的保证源于有限体积法直接基于积分形式的守恒定律，界面通量在相邻控制体积间共享。例如，对于两个相邻控制体积Vi和Vj，其公共界面f的通量Φf在Vi的离散方程中以+Φf出现，在Vj的方程中以-Φf出现，从而保证全局守恒（局部守恒的叠加即为全局守恒）。

2.推导二维非定常不可压缩流动的显式有限差分格式（采用交错网格，时间推进用向前欧拉法），并分析其稳定性条件。

考虑二维不可压缩流动的连续性方程和动量方程：

?u/?x+?v/?y=0

?u/?t+u?u/?x+v?u/?y=-1/ρ?p/?x+ν(?2u/?x2+?2u/?y2)

?v/?t+u?v/?x+v?v/?y=-1/ρ?p/?y+ν(?2v/?x2+?2v/?y2)

采用交错网格（StaggeredGrid），u定义在(i+1/2,j)，v定义在(i,j+1/2)，p定义在(i,j)。时间离散用向前欧拉法（显式），空间离散对对流项采用一阶迎风格式（避免奇性），扩散项用中心差分。

以u动量方程为例，离散形式为：

u^(n+1)_(i+1/2,j)=u^n_(i+1/2,j)+Δt[-u^n_(i+1/2,j)(u^n_(i+1,j)u^n_(i,j))/Δxv^n_(i+1/2,j)(u^n_(i+1/2,j+1)u^n_(i+1/2,j-1))/(2Δy)+ν((u^n_(i+3/2,j)2u^n_(i+1/2,j)+u^n_(i-1/2,j))/Δx2+(u^n_(i+1/2,j+1)2u^n_(i+1/2,j)+u^n_(i+1/2,j-1))/Δy2)(p^n_(i+1,j)p^n_(i,j))/(ρΔx)]

稳定性条件需考虑对流项和扩散项的限制。对于显式格式，CFL条件为：

Δt≤min(Δx/u_max,Δy/v_max)

同时，扩散项的稳定性由傅里叶数控制：

Δt≤(Δx2Δy2)/(2ν(Δx2+Δy2))

综合两者，最终稳定性条件为Δt≤min(CFL条件,扩散项限制)。

3.比较雷诺平均N-S方程（RANS）和大涡模拟（LES）在湍流模拟中的优缺点及适用场景。

RANS通过对N-S方程进行雷诺平均，将湍流脉动分解为平均量和脉动量，引入雷诺应力项（-ρu_iu_j），需用湍流模型（如k-ε、k-ωSST）封闭。优点：计算成本低（仅需求解平均量），适用于高雷诺数工程流动（如航空发动机内流、汽车气动设计）；缺点：无法解析湍流脉动细节，模型依赖性强（对复杂流动如分离流、转捩预测精度低）。

LES通过滤波分离大尺度涡和小尺度涡，直接求解大涡（含主要能量和流动特征），小涡用亚格子模型（如Smagorinsky模型、WALE模型）模拟。优点：能捕捉大尺度脉动，对复杂流动（如湍流边界层分离、涡脱落）的预测精度高于RANS；缺点：计算成本高（需加密网格解析大涡，时间步长更小），对近壁区处理仍需壁面模型（否则网格量激增），适用于中等雷诺数、需关注脉动特性的流动（如燃烧室湍流混合、风工程中的建筑绕流）。

4.说明二阶精度中心差分格式在计算对流项时的色散误差和耗散误差特性，并解释为何高马赫数流动中需采用迎风格式。

中心差分格式（CDS）对对流项?(uφ)/?x的离散为(u_iφ_iu_{i-1}φ_{i-1})/Δx（假设u0），其截断误差为O(Δx2)。通过傅里叶分析，将φ表示为谐波φ=e^(ikx)，离散格式的频散关系为：

kΔx=kΔx(kΔx)^3/6+...（色散误差项）

耗散误差（虚部）为0（中心差分无耗散）。因此，CDS的色散误差随波数k增大而显著（短波成分相位误差大），但无耗散。

高马赫数流动中，激波等强间断存在，短波成分丰富。中心差分的无耗散特性会导致间断附近出现数值振荡（如吉布斯现象），破坏解的稳定性。迎风格式（如一阶迎风格式、Roe格式）通过引入人工耗散（与速度方向相关的数值粘性），抑制高频振荡，同时保证激波捕捉的单调