[论文解读] Parallel variable-density particle-laden turbulence simulation
本文提出了一种基于MPI的并行C++代码，用于模拟具有辐射加热的可变密度、颗粒载流的湍流，采用双向耦合的拉格朗日-欧拉框架。主要贡献是针对低马赫数、可变密度流中产生的变系数泊松方程，提出了一种新型并行线性求解器，从而能够准确模拟辐射引起的湍流与浮力效应。
We have developed a fully parallel C++/MPI based simulation code for variable-density particle-laden turbulent flows. The fluid is represented through a uniform Eulerian staggered grid, while particles are modeled using a Lagrangian point-particle framework. Spatial discretization is second-order accurate, and time integration has a fourth-order accuracy. Two-way coupling of the particles with the background flow is considered in both momentum and energy equations. The code is fully modular and abstracted, and easily can be extended or modified. We have considered two different boundary conditions. We have also developed a novel parallel linear solver for the variable density Poisson equation that arises in the calculation.
研究动机与目标
- 开发一种可扩展、模块化且可扩展的模拟框架，用于可变密度、颗粒载流的湍流，具有辐射加热效应。
- 在低马赫数、可变密度流中，准确建模颗粒与流体之间的双向动量和能量耦合。
- 实现对复杂现象（如颗粒的优先聚集、优先扫掠以及辐射引起的湍流）的模拟。
- 为可变密度流中出现的变系数泊松方程实现一种新型并行线性求解器。
- 通过放宽Boussinesq近似，支持均匀各向同性湍流和浮力驱动流的模拟。
提出的方法
- 采用点粒子拉格朗日框架，对具有斯托克斯阻力及通过辐射和导热进行热能交换的颗粒进行建模。
- 在交错欧拉网格上使用二阶精度、交错的流体变量（密度、动量、温度）离散格式，时间积分采用四阶龙格-库塔方法。
- 应用线性插值和投影算子（I 和 P）实现欧拉场与拉格朗日场之间的双向耦合。
- 通过一种新型并行迭代求解器求解变系数泊松方程，以处理可变密度效应引起的流体压力。
- 通过改进的四阶龙格-库塔方案实现随时间演化的背景压力 $ P_0 $，以保持一致性并避免因数值误差导致的病态问题。
- 通过在交错网格上采用二阶、保守的有限差分格式，确保动能守恒。
实验结果
研究问题
- RQ1颗粒的辐射加热在颗粒载流湍流的发展中起到何种影响？
- RQ2可变流体密度对双向耦合颗粒载流中动量和能量传递的影响是什么？
- RQ3不采用Boussinesq近似对浮力驱动湍流的模拟有何影响？
- RQ4并行的、保守的数值格式能否准确求解低马赫数、可变密度流中的变系数泊松方程？
- RQ5辐射引起的热梯度如何改变颗粒的优先聚集特性及流体湍流结构？
主要发现
- 该代码成功再现了诸如颗粒优先聚集和优先扫掠等典型现象，验证了双向耦合框架的有效性。
- 模拟结果再现了先前在恒定密度、Boussinesq近似代码中观察到的辐射引起的湍流效应（Zamansky et al., 2014），现扩展至可变密度情形。
- 该代码在不使用Boussinesq近似的情况下，准确模拟了均匀浮力驱动湍流，其结果与Batchelor等人（1992）的研究一致。
- 开发了一种针对变系数泊松方程的新型并行线性求解器，并在可变密度流中验证了其数值稳定性与一致性。
- 通过改进的RK4方案计算随时间演化的背景压力 $ P_0 $，以确保与泊松方程的一致性，并避免因数值误差导致的病态问题。
- 该代码实现了二阶空间精度和四阶时间精度，通过在交错网格上的保守离散化确保了动能守恒。
更好的研究，从现在开始
从阅读论文到最终审阅，大幅缩短您的研究时间。
无需绑定信用卡
本解读由 AI 生成，并经人工编辑审核。