人体上呼吸道内气流运动特性的数值模拟分析

运用计算流体动力学(CFD)方法对人体上呼吸道内的气流运动特性进行了数值模拟,通过PIV实验对数值模拟结果进行了验证。研究结果表明:气流在咽部外壁、气管外壁发生分离现象,气流在气管内壁形成局部高速区,支气管内的气流在分叉处发生分离,靠近支气管内壁速度较高,并且在支气管边界层的外缘速度达到最大值。气管和支气管内的二次涡流运动和轴向速度的分布使得气管支气管内壁受到的剪应力较大,内壁粘膜更容易受到损伤。     

循环呼吸模式口喉模型内气流运动特性数值模拟

应用计算流体动力学方法对人体在循环呼吸模式下口喉模型内的气流运动特性进行数值模拟,分析了循环呼吸模式下口喉模型内的气流组织形式以及气流运动对呼吸道壁面以及气溶胶运动沉积的影响. 吸气阶段,在咽部外壁和声门下游气管上部外壁气流发生分离,形成分离区. 呼气阶段,分别在咽部外壁和喉部外壁形成高速区. 循环呼吸模式下,咽部、喉部与气管内的高轴向速度区和二次涡流运动均是在呼吸过程中间歇性的产生,形成的高剪应力区也是间歇性的. 壁面受到的剪应力周期性地改变方向,引起壁面劳损和组织损伤的可能性增大,同时在这些部位容易造成气溶胶的沉积.      

西安地区活动断裂活动性的Fuzzy Grey模式评价

应用计算流体动力学方法对人体在循环呼吸模式下口喉模型内的气流运动特性进行数值模 拟,分析了循环呼吸模式下口喉模型内的气流组织形式以及气流运动对呼吸道壁面以及气溶 胶运动沉积的影响. 吸气阶段,在咽部外壁和声门下游气管上部外壁气流发生分离,形成分 离区. 呼气阶段,分别在咽部外壁和喉部外壁形成高速区. 循环呼吸模式下,咽部、喉部与 气管内的高轴向速度区和二次涡流运动均是在呼吸过程中间歇性的产生,形成的高剪应力区 也是间歇性的. 壁面受到的剪应力周期性地改变方向,引起壁面劳损和组织损伤的可能性增 大,同时在这些部位容易造成气溶胶的沉积.     

颗粒物在人体上呼吸道模型中沉降规律的数值研究

本文采用RNG k-ε湍流模型对简化的上呼吸道内的流场进行数值仿真,并结合随机涡相互作用模型对上呼吸道内颗粒的局部运动沉积特性进行了数值模拟研究.结果表明,惯性冲击是微米级颗粒物主要的沉降机理,惯性参量能很好地衡量颗粒的沉降规律.口腔中流动较平稳,颗粒受扰动较小.颗粒由于惯性冲击而运动靠近咽部后部并易于沉积在咽部后壁面...     

超声速平面混合层小激波的形成与演变

为了揭示超声速混合层中小激波形成机理及其与涡相互作用的演变过程,本文基于大涡模拟（LES）方法,结合五阶精度混合TCD／WENO格式,对超声速平面混合层在对流马赫数为Mc=0．65条件下的流场结构进行了数值模拟,数值结果详细描述了超声速混合层中小激波的形成过程。研究了小激波形成后,随涡运动而产生的变形、脱落及发展过程。同时,对混合层双涡合并过程中,小激波与相邻涡相互作用所产生的变形与演变过程进行了讨论。     

送风方式对飞机座舱内污染物传播影响研究

采用数值模拟方法研究了波音737-200座舱模型内的气流场以及污染物传播过程,并通过实验数据加以验证;分析了天花板加侧壁送风(混合送风)、侧壁送风以及天花板送风三种送风方式下座舱内的速度场、涡量场、污染物浓度场和空气龄,并且解释了舱内流场结构对污染物传播过程的影响。研究结果表明:不同送风方式下飞机座舱内污染物的传播过程差异明显,这是由流场结构特征的差异造成的,特别是在小尺度范围内,涡结构十分复杂,严重影响污染物的传播规律;在聚集过程中,污染物更容易在漩涡范围内聚集,而变形主导的气流运动会阻碍污染物向更大范围传播;在排除过程中,天花板送风形式下污染物的排除效率最快,混合送风次之,侧壁送风最慢。这对于防控新冠肺炎是有参考价值的。     

水下发射航行体尾涡不稳定性分析

在水下连续发射过程中前一发航行体尾流会对后一发航行体运动姿态稳定性产生流动干扰现象. 因此, 研究尾流中涡旋结构演变机理对解决多弹体水下连续发射流动干扰难题具有重要的意义. 本文采用改进型分离涡模型与能量方程, VOF多相流模型与重叠网格技术相结合方法, 对航行体水下发射尾流演变过程开展精细化模拟研究, 其中模拟结果和实验吻合度较好, 验证了本文数值方法的有效性. 以航行体尾流区域为重点研究对象, 分析了尾流区瞬态流场分布, 讨论了横流强度和雷诺数对尾涡结构演变以及脉动压力分布特性的影响. 结果表明: 由于尾流区高速流体核心区与低速自由流相互作用导致Kelvin-Helmholtz不稳定现象出现, 可以清晰地发现涡旋结构在剪切力的作用下发生脱落. 在横流条件下, 航行体尾端脱落的涡环与涡腿形成发卡涡, 而多个发卡涡沿轴向间隔排列组成发卡涡包存在于尾流中. 随着横流强度增大, 形成多级发卡涡包结构, 而导致脉动压力二次峰值均出现的主要原因是尾流涡旋流场演变引起的. 随着雷诺数的增大, 尾流中由圆柱形涡和U型涡组成的二次涡结构逐渐明显, 不稳定性加强.      

基于三路舵机驱动的扑翼飞行模型气动特性研究

扑动而形成非定常气动现象是扑翼飞行过程中产生高升力的主要原因。本文以Ellington实验的鹰蛾翅膀为原形,设计扑翼实验及数值计算模型。通过压差传感器对翅膀模型上翼面固定位置进行测压,分析前缘涡的产生及脱落情况(考虑动压效应)。测量上下翼面固定位置处的压差,揭示扑翼飞行中产生高升力的主要原因。利用烟风洞观察扑翼模型周围流场结构及特殊涡产生变化情况。另外,根据Ellington提供的升力关系式估算了扑翼模型在一个周期内的平均升力。最后,基于三维欧拉方程对扑翼飞行气动特性进行数值模拟,计算结果与实验吻合良好。     

格栅-空腔自激振荡流场扰动形成机理研究

格栅-空腔流动会引发流场自激振荡现象,产生结构振动及噪声问题.研究此类现象的形成机理对相关设备的减振降噪设计具有重要意义.目前已知这种现象是在流场正逆向两种扰动的耦合作用下形成的,但上述两种扰动的形成机理尚未明确.针对该问题,建立了格栅-空腔流动数值模型,并进行非稳态数值模拟.通过分析压力振荡数据及流场结构演化过程,对两种扰动的形成过程进行了研究.研究结果表明,正逆向两种扰动的形成均与格栅间隔中形成的小尺度涡团有关.上游小尺度涡团从格栅间隔中依次脱离并聚集演化为大尺度涡团,形成该流动系统的正向扰动.下游小尺度涡团脱离格栅间隔激发的逆向扰流构成该流动系统的逆向扰动,并对正向扰动产生持续的正反馈作用.     

激波诱导环形SF6气柱演化的机理

基于可压缩多组分Navier-Stokes控制方程,结合5阶加权本质无振荡格式以及网格自适应加密技术和level-set方法,数值模拟了平面激波(Ma=1.23)与环形SF₆气柱(内外半径分别为8和17.5 mm)界面的相互作用过程。相比于之前的实验结果,数值模拟结果揭示了入射激波在界面内4次透射过程中的复杂波系结构,观察到透射激波在内部界面传播时形成自由前导折射结构并向自由前导冯诺依曼折射结构转换的波系演变过程;另外,界面内的复杂激波结构诱导内部下游界面上的涡量发生了3次反向;在界面演化后期,内部界面形成的“射流”结构与下游界面相互作用,诱导界面形成一对主涡、一对次级涡以及一个反向“射流”结构。定量分析了环形界面长度、宽度、位移、环量以及混合率的变化情况,结果表明,内部气柱的存在减弱了前期小涡结构合并形成大涡结构过程中对界面高度与长度的影响,同时提高了重质气体与环境气体的混合率。     

Numerical simulation of micro-particle deposition in a realistic human upper respiratory tract model during transient breathing cycle

An more reliable human upper respiratory tract model that consisted of an oropharynx and four generations of asymmetric tracheo-bronchial (TB) airways has been constructed to investigate the micro-particle deposition pattern and mass distribution in five lobes under steady inspiratory condition in former work by Huang and Zhang (2011). In the present work, transient airflow patterns and particle deposition during both inspiratory and expiratory processes were numerically simulated in the realistic human upper respiratory tract model with 14 cartilaginous rings (CRs) in the tracheal tube. The present model was validated under steady inspiratory flow rates by comparing current results with the theoretical models and published experimental data. The transient deposition fraction was found to strongly depend on breathing flow rate and particle diameter but slightly on turbulence intensity. Particles were mainly distributed in the high axial speed zones and traveled basically following the secondary flow. “Hot spots” of deposition were found in the lower portion of mouth cavity and posterior wall of pharynx/larynx during inspiration, but transferred to upper portion of mouth and interior wall of pharynx/larynx during expiration. The deposition fraction in the trachea during expiration was found to be much higher than that during inspiration because of the stronger secondary flow.     

Effects of cartilaginous rings on airflow and particle transport through simplified and realistic models of human upper respiratory tracts

In the present study, computational fluid dynamics （CFD） is used to investigate inspiratory and expiratory airflow characteristics in the human upper respiratory tract for the purpose of identifying the probable locations of particle deposition and the wall injury. Computed tomography （CT） scan data was used to reconstruct a three dimensional respiratory tract from trachea to first generation bronchi. To compare, a simplified model of respiratory tract based on Weibel was also used in the study. The steady state results are obtained for an airflow rate of 45 L/min, corresponding to the heavy breathing condition. The velocity distribution, wall shear stress, static pressure and particle deposition are compared for inspiratory flows in simplified and realistic models and expiratory flows in realistic model only. The results show that the location of cartilaginous rings is susceptible to wall injury and local particle deposition.     

Deposition of non-spherical microparticles in the human upper respiratory tract

We investigated the deposition pattern of microparticles with different particle diameters, shape factors, and initial flow conditions in a realistic human upper respiratory tract model. We identified a close relationship between the deposition fraction and the particle shape factor. The deposition fraction of the particles decreased sharply with increasing particle shape factor because of the decreasing drag force. We also found that the deposition varied at different positions in the upper respiratory tract. At low shape factors, the highest fraction of particles deposited at the mouth and pharynx. However, with increasing shape factor, the deposition fraction in the trachea and lungs increased. Moreover, for a given shape factor, larger particles deposited at the mouth and pharynx, which indicates that the deposition fraction of microparticles in the human upper respiratory tract is affected first and foremost by particle inertia as well as by the drag force.     

Numerical analysis of respiratory flow patterns within human upper airway

A computational fluid dynamics （CFD） approach is used to study the respiratory airflow dynamics within a human upper airway. The airway model which consists of the airway from nasal cavity, pharynx, larynx and trachea to triple bifurcation is built based on the CT images of a healthy volunteer and the Weibel model. The flow character- istics of the whole upper airway are quantitatively described at any time level of respiratory cycle. Simulation results of respiratory flow show good agreement with the clinical mea- sures, experimental and computational results in the litera- ture. The air mainly passes through the floor of the nasal cavity in the common, middle and inferior nasal meatus. The higher airway resistance and wall shear stresses are distrib- uted on the posterior nasal valve. Although the airways of pharynx, larynx and bronchi experience low shear stresses, it is notable that relatively high shear stresses are distrib- uted on the wall of epiglottis and bronchial bifurcations. Besides, two-dimensional fluid-structure interaction models of normal and abnormal airways are built to discuss the flow-induced deformation in various anatomy models. The result shows that the wall deformation in normal airway is relatively small.     

Modeling micro-particle deposition in human upper respiratory tract under steady inhalation

A representative human upper respiratory tract (URT) with idealized oral region and asymmetric tracheo-bronchial (TB) airway has been modeled, and laminar-to-turbulent airflow for typical inhalation modes as well as micro-particle transport and deposition has been simulated using CFX10.0 software from Ansys Inc. on a personal computer. The asymmetric TB airway could not be replaced by an extended straight tube as outlet of the oral region while investigating the tracheal airflow field and particle depositio...     

Numerical simulation and control of horseshoe vortex around an appendage-body junction

The horseshoe vortex generated around the appendage-body junction of submarines strongly influences the non-uniformity of submarine wakes at the propeller discs. The flow characteristics around the appended submarine body are numerically simulated and analyzed, and a new method on the vortex control baffle is presented. Then, the influence of the vortex control baffle on the horseshoe vortex generated at the sail-body junction is numerically studied, and the flow phenomena caused by the vortex control baffle with different transverse positions is investigated further. Results show that the vortex control baffle can induce a kind of attached vortex in a rotational direction opposite to the horseshoe vortex; these two kinds of vortices undermine each other. Furthermore, when the transverse position of the vortex control baffle is close to the horseshoe vortex, the state of the horseshoe vortex is directly affected, and the flow structure becomes even more complex. We adapt the vortex control baffle for the horseshoe vortex generated at the stern foil-body junction. Results from the numerical simulation of the flow around the fully appended submarine model indicate that the effect of the vortex control baffle greatly improves the performance of the submarine wake. The circumferential non-uniformity of the axial, tangential, and radial velocity components are decreased markedly. The engineering applicability of the vortex control baffle has been well presented.     

应用小干扰柱体控制角区马蹄涡结构的实验研究

本文提出一种简单的抑制和控制角区马蹄涡的被动控制方法.即在角区平板上游放置一个远小于主柱体的小干扰柱体,用其产生的弱马蹄涡来抑制和控制角区的马蹄涡结构.目的是使角区原来的马蹄涡结构由强变为弱、由大变为小、由多变为少、由非定常变为定常,以获得减小冲刷、抑制湍流、降低噪声、避免振动的工程效果.作者在风洞中采用烟线法和激光片光流动显示的方法开展研究,实验表明,在平板上游适当位置放置小干扰柱体的确可以有效抑制和控制角区马蹄涡结构.实验发现,当小干扰柱体放置在原角区马蹄涡生成区时,其抑制和控制效果最佳;当小干扰柱体放置在上游区或下游区时,控制效果不好.本文讨论了小干扰柱体控制角区马蹄涡的机理.此外,实验还研究了小干扰柱体相对尺度和截面形状对角区马蹄涡结构抑制和控制的影响.     

Some features of a turbulent wing–body junction vortical flow

Laser-Doppler velocimeter measurements of a wing/body junction flow field made within a plane to the side of the wing/wall junction and perpendicular both to a 3:2 elliptical nose—NACA 0020 tail wing, and a flat wall are presented. Reynolds number of the approach boundary layer was, Re_θ = 5940, and free-stream air velocity was, U_ref = 27.5 m/s. A large vortical structure residing in the outer region redirects the low-turbulence free-stream flow to the vicinity of the wing/wall junction, resulting in thin boundary layers with velocity magnitudes higher than free-stream flow. Lateral pressure gradients result in a three-dimensional separation on the uplifting side of the vortex. Additionally, a high vorticity vortical structure with opposite sense to the outer-layer vortex forms beneath the outer-layer vortex. Normal and shear stresses increase to attain values an order of magnitude larger compared to values measured in a three-dimensional boundary layer just outside the junction vortex. Bimodal histograms of the w fluctuating velocity occur under the outer-layer vortex near the wall due to the time-dependent nature of the horseshoe vortex. In such a flow the shear-stress angle (SSA) highly lags the flow-gradient angle (FGA), and the turbulence diffusion is highly altered due to presence of vortical structures.