眼内房水流动的数值研究

蔡建程; 张宝允; 曹月红; BRAZHENKO Volodymyr

doi:10.21656/1000-0887.410113

眼内房水流动的数值研究

doi: 10.21656/1000-0887.410113

¹浙江师范大学工学院, 浙江金华 321004;²金华市人民医院眼科, 浙江金华 321000

基金项目: 国家自然科学基金(51976201);金华市科技局重点项目（2019-3-016）

详细信息

作者简介:
蔡建程（1980—），男，副教授，博士（E-mail: cai_jiancheng@foxmail.com）;曹月红（1979—），女，硕士（通讯作者. E-mail: 25258581@qq.com）.

中图分类号: O368|Q66
计量
- 文章访问数: 1010
- HTML全文浏览量: 243
- PDF下载量: 236
- 被引次数: 0
出版历程
- 收稿日期: 2020-04-20
- 修回日期: 2020-09-16
- 刊出日期: 2021-02-01

Numerical Study of Aqueous Humor Flow in Human Eyes

¹College of Engineering, Zhejiang Normal University, Jinhua, Zhejiang 321004, P.R.China;²Department of Ophthalmology, Jinhua People’s Hospital,Jinhua, Zhejiang 321000, P.R.China

Funds: The National Natural Science Foundation of China（51976201）

摘要

摘要: 研究眼内房水流动有助于认识眼睛疾病（如青光眼）发生机理.该文利用计算流体力学计算了人眼睫状体分泌的房水由后房流经5μm和30μm虹膜-晶体间隙进入前房，再经小梁网排出这一过程，分析了平视、仰视、俯视不同体位以及不考虑浮力因素的眼内房水流场.结果表明不同情况的房水压力场类似，5μm虹膜-晶体间隙时后房眼压比前房约高30 Pa，而30μm时前后房压力基本相同，小梁网压力下降明显.前后房压差驱动房水从后房流经虹膜间隙到前房，而角膜与虹膜温差引起的自然对流为前房房水的主要流动.前房自然对流引起的压差变化在mPa量级，温度较高房水汇集区域压力稍高.平视时主要为前房虹膜侧房水上升、角膜侧下降的自然对流循环；仰视时瞳孔区房水沿前房中心轴上浮，然后沿角膜向周边下沉，从虹膜周边部回到瞳孔，形成轴对称的对流模式；对于俯视，流动模式与仰视正好相反，前房中心轴上房水从角膜运动到瞳孔处；不考虑浮力作用时，房水的平均速度比自然对流时小1~2个数量级.
- 房水流动 /
- 虹膜-晶体间隙 /
- 眼压 /
- 自然对流 /
- 数值研究
Abstract: The study of intraocular aqueous humor (AH) flow is helpful to understand the mechanism of some eye diseases such as glaucoma. A numerical study of AH flow inside human eyes with the computational fluid dynamics (CFD) were carried out to address the following flow processes: AH secreted in the posterior chamber by the ciliary body entering the anterior chamber through the iris-lens gap of 5 μm and 30 μm respectively, and discharging through the trabecular meshwork (TM). Detailed flow fields in different eye orientations were analyzed. Results show that, in general the intraocular pressure distributions are similar in all cases; the pressure in the posterior chamber is around 30 Pa higher than that in the anterior chamber with the 5 μm iris-lens gap, and they are almost equal with the 30 μm iris-lens gap. The pressure drop in the TM is noticeable. The pressure difference between the anterior and posterior chambers drives AH from the posterior chamber to the anterior chamber, while the temperature difference between cornea and iris surfaces causes the natural convection, which is the dominant flow in the anterior chamber. The pressure difference due to natural convection is in the magnitude of millipascal, and the higher-pressure region forms where warmer AH gathers. For the vertical orientation, warmer fluid rises along the iris surface and then turns downwards as it encounters the higher resistance in the upper TM regions. The flow then descends along the corneal surface toward the lower TM. For the upward-facing position, AH entering the pupil rises along the center line of the anterior chamber, and moves down along the cornea surface leading to 2 large symmetric recirculation zones. For the downward-facing position, the circulation route is opposite to that of the upward-facing position. With no buoyancy, the averaged velocity will be 1 to 2 orders of magnitude smaller than that with buoyancy.
- aqueous humor flow /
- iris-lens gap /
- intraocular pressure /
- natural convection /
- numerical study

HTML全文

参考文献(19)

[1]	FRIEDLAND A B. A hydrodynamic model of aqueous flow in the posterior chamber of the eye[J]. Bulletin of Mathematical Biology,1978,40(2): 223-235.
[2]	EL-SHAHED M, ABD ELMABOUD Y. On the fluid flow in the anterior chamber of a human eye with slip velocity[J]. International Communications in Heat and Mass Transfer,2005,32(8): 1104-1110.
[3]	CANNING C R, GREANEY M J, DEWYNNE J N, et al. Fluid flow in the anterior chamber of a human eye[J]. Mathematical Medicine and Biology: a Journal of the IMA,2002,19(1): 31-60.
[4]	FITT A D, GONZALEZ G. Fluid mechanics of the human eye: aqueous humour flow in the anterior chamber[J]. Bulletin of Mathematical Biology,2006,68(1): 53-71.
[5]	HEYS J J, BAROCAS V H, TARAVELLA M J. Modeling passive mechanical interaction between aqueous humor and iris[J]. Journal of Biomechanical Engineering,2001,123(6): 540-547.
[6]	HEYS J J, BAROCAS V H. Computational evaluation of the role of accommodation in pigmentary glaucoma[J]. Investigative Ophthalmology & Visual Science,2002,43(3): 700-708.
[7]	KUMAR S, ACHARYA S, BEUERMAN R, et al. Numerical solution of ocular fluid dynamics in a rabbit eye: parametric effects[J]. Annals of Biomedical Engineering,2006,34(3): 530-544.
[8]	VILLAMARIN A, ROY S, HASBALLA R, et al. 3D simulation of the aqueous flow in the human eye[J]. Medical Engineering & Physics,2012,34(10): 1462-1470.
[9]	CROWDER T R, ERVIN V J. Numerical simulations of fluid pressure in the human eye[J]. Applied Mathematics and Computation,2013,219(24): 11119-11133.
[10]	郭竞敏, 张虹, 王军明. 人眼前房三维重建与房水流场数值模拟[J]. 眼科新进展, 2015,35(4): 346-350.(GUO Jingmin, ZHANG Hong, WANG Junming. 3D model reconstruction and numerical simulation of fluid dynamics in anterior chamber[J]. Recent Advances in Ophthalmology,2015,35(4): 346-350.(in Chinese))
[11]	陈伟, 张向东, 余涵, 等. 虹膜膨隆对房水流动影响的数值模拟分析[J]. 眼科新进展, 2017,37(1): 72-76.(CHEN Wei, ZHANG Xiangdong, YU Han, et al. Numerical simulation of aqueous humor hydrodynamics in human eye with iris bombe[J]. Recent Advances in Ophthalmology,2017,37(1): 72-76.(in Chinese))
[12]	梁书秀, 郝菲菲, 孙昭晨, 等. 基于Fluent的白内障术后房水流体力学数值分析[J]. 暨南大学学报(自然科学与医学版), 2014,35(4): 350-356.(LIANG Shuxiu, HAO Feifei, SUN Zhaochen, et al. Numerical analysis of aqueous humor hydromechanics after cataract surgery based on Fluent[J]. Journal of Jinan University(Natural Science & Medicine Edition),2014,35(4): 350-356.(in Chinese))
[13]	WANG W, QIAN X, SONG H, et al. Fluid and structure coupling analysis of the interaction between aqueous humor and iris[J]. Biomedical Engineering Online,2016,15(2): 570-586.
[14]	HEYS J J, BAROCAS V H. A Boussinesq model of natural convection in the human eye and the formation of Krukenberg’s spindle[J]. Annals of Biomedical Engineering,2002,30(3): 392-401.
[15]	密斯让·J·C, 辛哈·A, 斯特·G·C, 等. 生物磁粘弹性流体的流动: 应用动脉电磁过热评估血液的流动,癌症治疗进程[J]. 应用数学和力学, 2010,31(11): 1330-1343.(MISRA J C, SINHA A, SHIT G C, et al. Flow of a biomagnetic viscoelastic fluid: application to estimation of blood flow in arteries during electromagnetic hyperthermia,a therapeutic procedure for cancer treatment[J]. Applied Mathematics and Mechanics,2010,31(11): 1330-1343.(in Chinese))
[16]	李方方, 刘静, 乐恺. 生物组织在冻结过程中的三维相变传热问题精确解[J]. 应用数学和力学, 2009,30(1): 64-74.(LI Fangfang, LIU Jing, YUE Kai. Exact analytical solution to three dimensional phase change heat transfer problems in biological tissues subject to freezing[J]. Applied Mathematics and Mechanics,2009,30(1): 64-74.(in Chinese))
[17]	BERGMAN T L, LAVINE A S, INCROPERA F P,et al. Fundamentals of Heat and Mass Transfer [M]. Hoboken: John Wiley & Sons, Inc, 2017.
[18]	AVTAR R, SRIVASTAVA R. Modelling aqueous humor outflow through trabecular meshwork[J]. Applied Mathematics and Computation,2007,189(1): 734-745.
[19]	陶文铨. 数值传热学[M]. 西安: 西安交通大学出版社, 2001.(TAO Wenquan. Numerical Heat Transfer [M]. Xi’an: Xi’an Jiaotong University Press, 2001.(in Chinese))