鼻腔模拟采用的边界条件的设置以及相关的文献

在COMPUTATIONAL FLUID DYNAMICS(CFD)  APPLIED IN THE DRUGDELIVERYDESIGNPROCESS  TOTHENASALPASSAGES:AREVIEW这篇文章中提到，

1. 对于大多数情况，无需解析湍流脉动的细节，因其需要非常密集的网格，通常只研究平均流的影响。雷诺平均纳维 - 斯托克斯方程（RANS）这类时均方程是湍流建模的有效方法，通过时间平均确定湍流对平均流场的影响，这种平均忽略了湍流的随机特性并产生额外的雷诺应力。RANS 方法在实际 CFD 模拟湍流时广泛应用，如 k - ε 和 k - ω 模型。使用 RANS 方程模拟时，需结合涡相互作用模型（EIM），但对于惯性冲击参数 的情况，RANS/EIM 方法计算的颗粒沉积量比实验数据显著偏高，而较重颗粒受气流波动影响较小 。
2. 许多鼻腔通道的研究假设为层流条件，但部分研究者认为在湍流条件下研究更有必要。涉及湍流模型的研究多采用 RNG k - ε 模型，试图考虑低雷诺数区域；也有研究使用 k - ω 模型等。Liu 等人使用 SST k - ω 模型并进行大涡模拟（LES），LES 能捕捉 RANS/EIM 无法再现的近壁区域流动特征，效果最佳。此前研究表明鼻腔气道会产生湍流，意味着需进行湍流建模，但也有观点认为健康成年人静息呼吸时气流主要是层流，部分研究假设层流也能得到准确预测结果。
3. 研究发现，对于所研究的双向输送装置，是否应用湍流模型的结论高度依赖于注入的颗粒大小。对于评估的颗粒尺寸，沉积仅发生在入口鼻孔，但较小尺寸颗粒会绕过鼻中隔在另一侧鼻孔沉积，当气流转向时，启用湍流模型的影响更显著。Mylavarapu 等人验证用于人体上呼吸道流动模拟的 CFD 方法时发现，标准 k - ω 湍流模型与静压测量结果吻合最佳，但该模拟是在呼气峰值流量 200L/min 下进行的。

1. 边界条件的重要性及进出口边界设定：边界条件在实际现象的数学表达中至关重要。许多研究者通过定义进出口边界的压力来模拟现实中压力驱动的呼吸，如 Garcia 等人（2007）的做法；Ozlugediks 等人（2008）则通过在出口设置流出边界条件来产生一定的吸气气流，他们认为这种方式能更好地模拟从鼻咽部抽吸空气的生理条件，使流经每个鼻孔的流量能适当自我调节。对于药物输送相关计算，通常采用主动进口边界条件模拟气流流入。壁面边界条件常设为无滑移条件，即壁面处速度为零。
2. 粒子释放条件：在模拟中，通常假设粒子悬浮在空气中，在鼻孔入口处以气流速度均匀释放。在模拟不同药物输送方式时，常将装置插入鼻腔，粒子会从装置中带着动量射出，定义计算边界条件时需考虑这一点。此外，粒子形状、注入粒子或液滴的尺寸分布在研究中也很重要，目前多数研究将粒子注入模拟为均匀分布。
3. 粒子沉积边界条件及挑战：在沉积研究中，许多研究者定义的边界条件是当粒子接近鼻腔壁的距离小于其半径时，粒子会被壁面捕获。在计算粒子轨迹时，通常假设壁面被黏液覆盖，直至粒子沉积在壁面或离开计算域。由于鼻腔中的黏液是多种细胞分泌物的混合物，且存在黏液纤毛清除和鼻周期等生理现象，所以正确定义边界条件颇具挑战

4. Fluent 计算 操作步骤

   File-Read-Mesh

   Display —可以看网格有几个面

   scale—几何尺寸（检查一下网格有没有问题）

   solver—type—pressure based

   time base—可选稳态 也可以选非稳态（看文献）

   考虑重力——gravity z轴—（-9.8）

   model 板块 除了Viscous全部关闭

   Viscous—15L以下用laminar层流模型 15L以上用k-omega-SST模型 选择低流动修正（15L指的是两个鼻孔总流量）

   disCrete phase—释放颗粒模型

   最大步数的选择是根据最后的计算结果中escap、incomplete的颗粒来决定的，如果有逃逸和未完成的颗粒，则将步数调大

   injection创建一个释放颗粒的面，创建在两个鼻孔的表面，释放类型type-surface，材料类型可以随便选择，重点是改变**材料密度（密度选择依据文献，多组分水滴的最终尺寸会不低于纯水滴），**密度的更改就在material，修改完成之后必须点击change，Diameter distribute指的是输入颗粒分布，uniform是单一直径的颗粒，R-R是符合R分布正态分布，直径的考虑（20 to 200 µm）选择injection using face normal direction

   0.9524 for liquid water and 0.0476 for the non-volatile viral species.

   cell zone comdition空气流体域 空气

   boundary condition 边界条件是由 进口速度还是进口压力决定的，进口速度有三种类型

   pressure -inlet 从面进入，pressure -outlet 从面出去，进口压力设置为0，出口压力是-4.408268，最后后处理，需要看出口面流速，压差越大流速越大

   壁面类型，可以不用设置壁面，wall 其他部分都选择trap，嗅区部分选择escape

   methods 选simple

   算流场 直接初始化 initialization，迭代完成

   计算 run calculation 初始化10000步 calculate

   inlet和outlet 选prevent reverse flow 回流