基于直接力浸入边界法的颗粒沉降特性：形状效应的深度数值剖析.docxVIP

基于直接力浸入边界法的颗粒沉降特性：形状效应的深度数值剖析

一、引言

1.1研究背景与意义

颗粒沉降现象广泛存在于自然界和众多工业领域，对其深入研究具有至关重要的意义。在自然界中，大气中的尘埃沉降影响空气质量与气候变化，海洋中颗粒物的沉降对海洋生态系统和生物地球化学循环起着关键作用，如南海颗粒沉降对微生物群落结构及其垂向连通性有塑造作用。在工业领域，颗粒沉降更是涉及多个方面。在石油开采中，了解颗粒沉降有助于提高原油分离效率，降低生产成本；在污水处理中，通过掌握颗粒沉降规律，能够优化沉淀工艺，有效去除污水中的悬浮物，提升水质；在矿物加工里，颗粒沉降对于矿石的分选和提纯十分关键，直接影响矿产资源的利用率。

在模拟颗粒沉降时，直接力浸入边界法（DirectForcingImmersedBoundaryMethod）展现出独特优势。传统的数值模拟方法在处理复杂边界和流固耦合问题时存在诸多困难，如网格生成复杂、计算效率低等。而直接力浸入边界法通过将颗粒边界视为流场中的力源项，巧妙地避免了复杂的网格划分过程。它采用笛卡尔网格，无需根据颗粒形状生成贴体网格，极大地提高了计算效率，尤其适用于处理颗粒形状复杂多变的情况。此外，该方法能够精确捕捉颗粒与流体之间的相互作用力，为准确模拟颗粒沉降过程提供了有力手段。

不同形状颗粒的沉降特性存在显著差异，研究这些差异对于深入理解颗粒沉降机理至关重要。颗粒形状会影响其在流体中的受力情况，进而影响沉降速度和轨迹。例如，圆形颗粒在沉降时受力相对均匀，沉降轨迹较为稳定；而不规则形状的颗粒，由于其形状的复杂性，受力不均匀，可能会出现翻滚、摆动等复杂运动，沉降轨迹也更加复杂。研究不同形状颗粒的沉降，有助于建立更准确的颗粒沉降模型，为相关工程应用提供更可靠的理论依据。同时，这对于优化工业生产过程、提高产品质量和生产效率具有重要的指导意义，还能为解决自然界中的相关问题提供新的思路和方法。

1.2国内外研究现状

直接力浸入边界法的发展经历了多个重要阶段。1972年，CharlesPeskin为模拟血液在可收缩心脏瓣膜中的流动首次提出浸入边界法思想，将复杂结构边界模化成Navier-Stokes动量方程的力源项。此后，该方法不断发展完善。Peskin精心设计的δ函数对方法效率提升至关重要，在后续研究中，如参考文献[12,4,9]以及相关研究[13,14]，奇异力通过反馈机制获得，被称为“虚拟边界法”，但此方法引入额外自由参数，且需解决弹簧阻尼系统振荡的特征时间尺度问题，可能导致时间步长受限。为克服弊端，Fadlun等人引入力项的直接公式，修改离散动量方程隐式矩阵条目以获边界点所需速度，该方案不受虚拟边界方法时间步长限制。金等人后来提出直接强迫方法的显式变体，保持标准有限差分方法原始简单矩阵结构。目前，直接力浸入边界法已广泛应用于多相流模拟、生物流体动力学研究等领域，在模拟颗粒沉降方面也逐渐受到关注。

在不同形状颗粒沉降模拟研究方面，取得了一系列进展。对于二维圆形颗粒沉降，诸多学者采用有限元方法、格子Boltzmann方法（LBM）等进行研究，如文献[1]采用有限元方法深入研究单颗圆形颗粒在牛顿流体中的沉降规律（雷诺数≤600），文献[2]通过LBM对双圆形颗粒在偏离通道中心位置的沉降形态进行数值研究。在椭圆形颗粒沉降研究中，也有不少成果。对于矩形颗粒、三维球形颗粒以及椭球颗粒等的沉降模拟也有涉及。在研究柱状催化剂颗粒相对于球形颗粒的不同运动特性时，选择不同长度直径为2mm的5种柱状颗粒，采用CFD-DEM数值模拟仿真，研究柱状颗粒在管状容器中沉降的运动行为，并建立柱状颗粒沉降试验台，采用高速摄像拍摄的方法进行实验研究。然而，当前研究仍存在不足。一方面，对于复杂形状颗粒，如珊瑚砂颗粒（具有比重大、磨圆度低、形状各异、棱角分明且空隙较多等特点），现有沉降预测模型适用度较低；另一方面，在颗粒与流体相互作用的精细化模拟上还有待加强，不同形状颗粒沉降过程中多物理场耦合效应的研究也相对较少。

1.3研究目标与内容

本研究旨在运用直接力浸入边界法，对不同形状颗粒在流体中的沉降过程进行高精度的直接数值模拟，深入探究颗粒沉降特性以及颗粒与流体之间的相互作用机制。具体研究内容如下：

不同形状颗粒沉降特性研究：系统研究圆形、椭圆形、矩形、柱状等多种典型形状颗粒在不同工况下的沉降特性。通过数值模拟，获取颗粒的沉降速度、加速度、轨迹等关键参数，分析颗粒形状对这些参数的影响规律。例如，对比圆形与椭圆形颗粒在相同流体条件下的沉降速度变化，研究形状的扁率如何影响沉降过程。

颗粒与流体相互作用研究：深入分析不同形状颗粒在沉降过程中与流体之间的相互作用力，包括曳力、升力等