本文選題：上呼吸道 + 氣流運(yùn)動(dòng)��； 參考：《中國人民解放軍軍事醫(yī)學(xué)科學(xué)院》2008年博士論文

【摘要】： 空氣污染物氣溶膠對(duì)人類健康產(chǎn)生嚴(yán)重的威脅,人體通過呼吸道吸入有毒氣溶膠會(huì)引發(fā)哮喘、肺氣腫和支氣管炎等呼吸道疾病,SARS和禽流感等疾病的爆發(fā)同樣是病毒以生物氣溶膠的形式通過呼吸道傳播而引發(fā)感染;針對(duì)各種呼吸道疾病,氣溶膠吸入治療在呼吸道疾病的防治中表現(xiàn)出了明顯的優(yōu)勢。人體上呼吸道內(nèi)氣流流場、氣溶膠性質(zhì)、呼吸模式及其幾何特性決定了有毒氣溶膠或吸入藥物氣溶膠的沉積位置和局部濃度,進(jìn)而決定有毒氣溶膠的危害程度或藥物氣溶膠的治療效果。因此,研究人體呼吸道內(nèi)的氣流運(yùn)動(dòng)特性,探討有毒氣溶膠或藥物氣溶膠在人體呼吸道內(nèi)的沉積規(guī)律,對(duì)于認(rèn)識(shí)有毒氣溶膠對(duì)人體的危害、進(jìn)行劑量健康效應(yīng)的評(píng)價(jià)以及繼續(xù)深入探索有毒氣溶膠的致病機(jī)理具有重要的實(shí)際意義,對(duì)于提高藥物氣溶膠的治療效果具有重要的指導(dǎo)意義。 本文在總結(jié)國內(nèi)外研究工作的基礎(chǔ)上,借鑒ARLA(Aerosol Research Laboratory of Alberta)的理想口-喉模型和Weible模型的氣管-支氣管模型,建立了具有口腔-咽-喉-氣管-前三級(jí)支氣管的完整的人體上呼吸道模型。運(yùn)用CFD數(shù)值仿真和試驗(yàn)研究相結(jié)合的方法對(duì)人體上呼吸道內(nèi)的氣流運(yùn)動(dòng)特性和氣溶膠運(yùn)動(dòng)沉積規(guī)律進(jìn)行了全面、系統(tǒng)的研究。研究結(jié)果表明: 穩(wěn)態(tài)呼吸模式下,氣流在咽部外壁、氣管外壁發(fā)生分離現(xiàn)象,氣流在氣管內(nèi)壁形成局部高速區(qū)域,容易造成較多的氣溶膠沉積;在氣管內(nèi)的三個(gè)截面內(nèi)分別形成兩個(gè)對(duì)稱的二次渦流運(yùn)動(dòng),二次渦流運(yùn)動(dòng)使得氣管內(nèi)壁所受到的流動(dòng)剪應(yīng)力增大,而外壁面剪應(yīng)力減小;同時(shí)軸向速度在氣管內(nèi)壁引起高剪應(yīng)力分布,而外壁的剪應(yīng)力較小,二次渦流容易造成氣溶膠在氣管內(nèi)壁沉積較多;進(jìn)入到支氣管內(nèi)的氣流在分叉處發(fā)生分離,并且在下游支氣管的內(nèi)壁區(qū)域形成新的邊界層,靠近支氣管內(nèi)壁速度較高,并且在支氣管邊界層的外緣速度達(dá)到最大值,靠近支氣管外壁速度較低。 循環(huán)吸氣模式下,吸氣加速階段咽部、喉部和聲門下游氣管內(nèi)壁形成局部高速區(qū),容易造成氣溶膠因慣性碰撞而沉積;在咽部外壁、聲門下游氣管上部外壁氣流逐漸發(fā)生分離,形成分離區(qū),使得氣溶膠在這些部位隨氣流環(huán)流循環(huán)運(yùn)動(dòng),滯留時(shí)間增加,容易造成一部分氣溶膠在咽部外壁和氣管外壁的沉積;支氣管平面內(nèi)形成拋物線型的速度分布,氣流主流逐漸偏離支氣管外壁,流向內(nèi)壁,容易造成氣溶膠在支氣管內(nèi)壁沉積較多;支氣管內(nèi)各個(gè)截面上形成的二次渦流運(yùn)動(dòng)的強(qiáng)度逐漸增強(qiáng),二次渦流運(yùn)動(dòng)會(huì)造成氣溶膠在內(nèi)壁的沉積幾率增加。吸氣減速階段,咽部外壁的分離區(qū)逐漸向內(nèi)壁方向擴(kuò)大,同時(shí)在咽部的內(nèi)壁也發(fā)生氣流分離現(xiàn)象,氣管外壁的分離區(qū)向氣管內(nèi)壁方向擴(kuò)大,同時(shí)向氣管下游延伸;在第一級(jí)支氣管外壁區(qū)產(chǎn)生氣流分離現(xiàn)象,形成較小的分離區(qū),在第二級(jí)和第三級(jí)支氣管內(nèi)沒有發(fā)生類似現(xiàn)象。 循環(huán)呼氣模式下,呼氣加速階段與吸氣加速階段不同,氣管平面內(nèi)氣流流速分布比較均勻;氣管內(nèi)逐漸形成四個(gè)二次渦流運(yùn)動(dòng),氣流在支氣管內(nèi)經(jīng)過多級(jí)分支的效果是使速度分布均勻化,特別是在同級(jí)兩個(gè)支氣管軸線共處的平面內(nèi)的速度分布在經(jīng)過幾級(jí)匯合后,軸線處的峰值分布現(xiàn)象會(huì)消失而變得均勻;在支氣管內(nèi)呈現(xiàn)典型的拋物線速度分布形式,支氣管內(nèi)的二次氣流運(yùn)動(dòng)經(jīng)歷了從兩個(gè)渦流到四個(gè)渦流的運(yùn)動(dòng)變化過程。呼氣減速階段,呼氣終了時(shí)刻,在第一級(jí)和第二級(jí)支氣管的內(nèi)壁均出現(xiàn)了氣流分離現(xiàn)象,在內(nèi)壁形成分離區(qū)。 循環(huán)呼吸模式下,咽部、喉部、氣管以及支氣管內(nèi)高軸向速度區(qū)和二次渦流運(yùn)動(dòng)均是在呼吸過程中逐漸出現(xiàn)的,只是間歇性地產(chǎn)生,所以由此而引起的氣道壁面的氣流剪應(yīng)力集中,形成的高剪應(yīng)力區(qū)也是間歇性的,只在整個(gè)周期的部分時(shí)間出現(xiàn)。壁面受到的剪應(yīng)力周期性地改變方向,引起壁面勞損和組織損傷的可能性增大,同時(shí)在這些部位容易造成氣溶膠的沉積,還可能會(huì)引起各種呼吸道疾病。呼氣階段支氣管內(nèi)的氣流運(yùn)動(dòng)、氣流剪應(yīng)力的分布和氣溶膠的運(yùn)動(dòng)形式比吸氣階段更為復(fù)雜。 采用激光快速成型技術(shù)(Stereo-Lithography,SL)制作了人體上呼吸道的試驗(yàn)?zāi)Ｐ?應(yīng)用粒子圖像速度儀(Particle Image Velocimetry,PIV)對(duì)人體上呼吸道內(nèi)的穩(wěn)態(tài)氣流運(yùn)動(dòng)特性進(jìn)行了試驗(yàn)研究,分析了在低強(qiáng)度呼吸條件下人體上呼吸道內(nèi)的氣流運(yùn)動(dòng)特性。研究結(jié)果表明,數(shù)值仿真結(jié)果與試驗(yàn)測量結(jié)果基本一致,證明了數(shù)值仿真方法的準(zhǔn)確性和合理性。 利用拉格朗日方法對(duì)氣溶膠在上呼吸道內(nèi)的運(yùn)動(dòng)進(jìn)行仿真計(jì)算,分析了不同呼吸模式下氣溶膠的沉積特點(diǎn)。慣性碰撞對(duì)于微尺度氣溶膠沉積而言是主要的沉積機(jī)制,慣性參數(shù)是衡量碰撞作用造成顆粒沉積的一個(gè)重要的參數(shù),人體上呼吸道內(nèi)不同部位氣溶膠沉積率隨慣性參數(shù)的增加而增加;而湍流擴(kuò)散、二次氣流運(yùn)動(dòng)和環(huán)流氣流運(yùn)動(dòng)對(duì)氣溶膠在人體上呼吸道內(nèi)的沉積同樣具有重要的影響,人體的呼吸流量和氣溶膠性質(zhì)對(duì)氣溶膠在上呼吸道內(nèi)的沉積模式影響較小;慣性碰撞和湍流擴(kuò)散的影響致使在喉部氣溶膠沉積最多,氣管中氣溶膠的沉積效率要高于支氣管中的氣溶膠沉積效率。人體循環(huán)吸氣模式下,氣溶膠在人體上呼吸道內(nèi)的沉積率要高于穩(wěn)態(tài)吸氣情況下的氣溶膠的沉積率;循環(huán)吸氣模式下的氣溶膠沉積率遠(yuǎn)大于循環(huán)呼氣模式下的氣溶膠沉積率。 建立了氣溶膠在人體上呼吸道內(nèi)沉積的試驗(yàn)臺(tái),對(duì)氣溶膠在上呼吸道內(nèi)的沉積進(jìn)行了試驗(yàn)研究。數(shù)值仿真結(jié)果與試驗(yàn)結(jié)果基本一致,平均誤差為11%,沉積變化趨勢吻合較好。人體上呼吸道內(nèi)氣溶膠沉積的數(shù)值仿真方法,能夠較好地預(yù)測氣溶膠在上呼吸道內(nèi)的沉積模式以及不同部位的沉積率,是我們獲得上呼吸道內(nèi)有毒氣溶膠或藥物氣溶膠在不同部位沉積信息的一種有效的方法。
[Abstract]:Air pollutant aerosols pose a serious threat to human health. Inhalation of aerosols through the respiratory tract causes respiratory diseases such as asthma, emphysema and bronchitis. The outbreak of SARS and avian influenza is the same as the infection of the virus in the form of biological aerosols through the breathing tract; for various respiratory diseases. Aerosol inhalation therapy has shown obvious advantages in the prevention and control of respiratory diseases. The flow field, aerosol properties, breathing patterns and geometric properties of the upper respiratory tract determine the location and local concentration of toxic aerosols or inhaled aerosol aerosols, and determine the degree of harm or drug gas solubility of the aerosol. Therefore, the characteristics of the air flow in the human respiratory tract are studied, and the depositional laws of toxic aerosol or drug aerosol in the human respiratory tract are discussed. It is important to recognize the harm of the toxic aerosol to the human body, to evaluate the health effect of the dose and to further explore the pathogenic mechanism of the toxic aerosol. The intertemporal significance has important guiding significance for improving the therapeutic effect of aerosol.
On the basis of the research work at home and abroad, the complete human upper respiratory tract model with oral pharynx larynx trachea and the first three bronchus is established using the ideal oral throat model of ARLA (Aerosol Research Laboratory of Alberta) and the tracheobronchial model of Weible model. The combination of CFD numerical simulation and experimental research is combined. The characteristics of the airflow and the rule of aerosol movement and deposition in the upper respiratory tract of the human body were studied comprehensively and systematically.
In the steady breathing mode, the air flow is separated from the outer wall of the pharynx and the outer wall of the trachea is separated. The air flow in the inner wall of the trachea is formed in a local high-speed region, which can easily cause more aerosol deposition, and two symmetrical two swirl motions are formed in the three sections of the trachea, and the flow shear stress on the inner wall of the trachea is caused by the two vortexes. In addition, the shear stress of the outer wall decreases, while the axial velocity causes the high shear stress distribution on the inner wall of the trachea, while the shear stress of the outer wall is smaller. The two swirl can easily cause more deposition of aerosols in the inner wall of the trachea, and the air flow into the bronchus is separated at the fork, and a new boundary layer is formed in the inner wall area of the downstream bronchi. The velocity near the inner wall of the bronchus is high, and the velocity at the outer edge of the bronchial boundary layer reaches the maximum, and the velocity near the bronchial outer wall is relatively low.
In the cycle suction mode, a local high-speed region is formed in the throat of the suction acceleration stage, the throat and the inner wall of the trachea at the lower reaches of the glottis, which can easily cause the aerosol to be deposited by the inertia collision, and the outer wall of the upper wall of the trachea in the lower throat of the pharynx gradually separates and forms a separation area, causing the aerosol to be detained in these parts with the circulation circulation of the air flow. When the time increases, a part of the aerosol is deposited in the outer wall of the pharynx and the outer wall of the trachea; the velocity distribution of the parabolic type is formed in the plane of the bronchus, and the mainstream of the air flow gradually deviates from the outer wall of the bronchus and flows into the inner wall. It is easy to cause more deposition of aerosols in the inner wall of the bronchus, and two swirl movements formed on each section of the trachea. The strength of the two swirl will increase the deposition probability of the aerosol on the inner wall. The separation area of the outer wall of the pharynx expands gradually to the inner wall, while the inner wall of the pharynx also occurs in the inner wall of the pharynx. The separation area of the outer wall of the trachea extends to the inner wall of the trachea and extends downstream of the trachea at the same time. There was a phenomenon of air separation in the outer wall of the bronchus, and a smaller separation area was formed. There was no similar phenomenon in grade second and grade third bronchi.
In the cycle expiration mode, the air flow velocity distribution in the trachea plane is more uniform in the exhaled exhalation stage than in the inhalation acceleration stage, and four two swirl motions are formed in the trachea, and the effect of the air flow through the multistage branch in the bronchus is to homogenization the velocity distribution, especially in the plane of the two bronchial axes cooperating at the same level. The peak distribution in the axis will disappear after the convergence of the degree distribution at several levels. The typical parabolic velocity distribution in the bronchus shows a typical parabolic velocity distribution in the bronchus. The two flow of air in the bronchus has undergone the movement from two swirl to four eddy current. The air separation phenomenon appeared in the inner wall of the two grade bronchus, and the separation area was formed in the inner wall.
In the cycle breathing mode, the throat, larynx, trachea, and the high axial velocity zone and two swirl movement in the bronchus are all gradually appearing during the breathing process, only intermittent real estate, so the air shear stress concentration of the airway wall is concentrated, and the high shear stress area is intermittent, only in the whole period. The shear stress periodically changes the direction of the wall and causes the increase in the possibility of wall strain and tissue damage. At the same time, it is easy to cause the deposition of aerosol and may cause various respiratory diseases. The airflow movement in the bronchus, the distribution of the air flow shear stress and the form ratio of the aerosol movement in the expiratory stage The inhalation stage is more complex.
The experimental model of the human upper respiratory tract was made by Stereo-Lithography (SL). The dynamic characteristics of the steady air flow in the upper respiratory tract were studied by the Particle Image Velocimetry (PIV), and the airflow in the upper respiratory tract under low intensity respiration was analyzed. The results show that the numerical simulation results are basically consistent with the experimental results, which proves the accuracy and rationality of the numerical simulation method.
The Lagrange method is used to simulate the motion of aerosols in the upper respiratory tract, and the characteristics of aerosol deposition in different breathing patterns are analyzed. Inertial collisions are the main deposition mechanism for microscale aerosol deposition. The inertial parameters are an important parameter to measure the particle deposition caused by collision, and the human body calls on the body. The aerosol deposition rate in the different parts of the suction channel increases with the increase of the inertial parameters, while the turbulent diffusion, the two air flow and the circulation flow have an important influence on the deposition of aerosols in the upper respiratory tract. The respiratory flow and aerosol properties of the human body have little influence on the deposition pattern of gas soluble glue in the upper respiratory tract. The effect of inertial and turbulent diffusion causes the most deposition in the larynx. The deposition efficiency of the aerosol in the trachea is higher than the aerosol deposition efficiency in the bronchus. The deposition rate of aerosols in the human body's upper respiratory tract is higher than that in the steady inhalation condition. The aerosol deposition rate is much larger than the aerosol deposition rate under cyclic expiratory mode.
An experimental platform for aerosol deposition in the upper respiratory tract was established, and the deposition of aerosols in the upper respiratory tract was experimentally studied. The numerical simulation results are basically consistent with the experimental results, the average error is 11%, and the trend of deposition changes well. The numerical simulation method of aerosol deposition in the upper respiratory tract can be well predicted. The deposition patterns of aerosols in the upper respiratory tract and the deposition rate in different parts of the upper respiratory tract are an effective method for obtaining the information of the aerosol or aerosol in the upper respiratory tract in different parts of the respiratory tract.
【學(xué)位授予單位】：中國人民解放軍軍事醫(yī)學(xué)科學(xué)院
【學(xué)位級(jí)別】：博士
【學(xué)位授予年份】：2008
【分類號(hào)】：R363

【參考文獻(xiàn)】

相關(guān)期刊論文 前9條

1 朱科倫;朱郇憫;;新發(fā)現(xiàn)的傳染病的特點(diǎn)[J];廣州醫(yī)藥;2006年01期

2 喬偉;黃昌麗;于永軍;袁欣;;城市空氣污染控制與管理[J];黑龍江環(huán)境通報(bào);2006年03期

3 孫秀珍;于馳;劉迎曦;于申;張軍;蘇英鋒;;人體上呼吸道三維有限元重建與流場數(shù)值模擬[J];航天醫(yī)學(xué)與醫(yī)學(xué)工程;2006年02期

4 熊志明,張國強(qiáng),彭建國,周舟;大氣可吸入顆粒物對(duì)IAQ的影響研究進(jìn)展[J];建筑熱能通風(fēng)空調(diào);2004年02期

5 王姝;;空氣污染對(duì)人體健康影響及經(jīng)濟(jì)損失估算[J];遼寧城鄉(xiāng)環(huán)境科技;2005年06期

6 趙彬,張穎,李先庭;生物顆粒在相鄰房間運(yùn)動(dòng)的數(shù)值研究[J];暖通空調(diào);2003年03期

7 張楚華;聞蘇平;劉陽;;人體呼吸道的二級(jí)及三級(jí)支氣管內(nèi)吸氣流動(dòng)的數(shù)值研究[J];生物醫(yī)學(xué)工程學(xué)雜志;2006年04期

8 張鴻雁,王麗,王元;二維溫室繞流流場的PIV實(shí)驗(yàn)研究[J];西安建筑科技大學(xué)學(xué)報(bào)(自然科學(xué)版);2005年02期

9 曾敏捷;胡桂林;樊建人;;微顆粒在人體上呼吸道中運(yùn)動(dòng)沉積的數(shù)值模擬[J];浙江大學(xué)學(xué)報(bào)(工學(xué)版);2006年07期

，

本文編號(hào)：1914518