Experimental validation of convection-diffusion discretisation scheme employed for computational modelling of biological mass transport
Tóm tắt
The finite volume solver Fluent (Lebanon, NH, USA) is a computational fluid dynamics software employed to analyse biological mass-transport in the vasculature. A principal consideration for computational modelling of blood-side mass-transport is convection-diffusion discretisation scheme selection. Due to numerous discretisation schemes available when developing a mass-transport numerical model, the results obtained should either be validated against benchmark theoretical solutions or experimentally obtained results. An idealised aneurysm model was selected for the experimental and computational mass-transport analysis of species concentration due to its well-defined recirculation region within the aneurysmal sac, allowing species concentration to vary slowly with time. The experimental results were obtained from fluid samples extracted from a glass aneurysm model, using the direct spectrophometric concentration measurement technique. The computational analysis was conducted using the four convection-diffusion discretisation schemes available to the Fluent user, including the First-Order Upwind, the Power Law, the Second-Order Upwind and the Quadratic Upstream Interpolation for Convective Kinetics (QUICK) schemes. The fluid has a diffusivity of 3.125 × 10-10 m2/s in water, resulting in a Peclet number of 2,560,000, indicating strongly convection-dominated flow. The discretisation scheme applied to the solution of the convection-diffusion equation, for blood-side mass-transport within the vasculature, has a significant influence on the resultant species concentration field. The First-Order Upwind and the Power Law schemes produce similar results. The Second-Order Upwind and QUICK schemes also correlate well but differ considerably from the concentration contour plots of the First-Order Upwind and Power Law schemes. The computational results were then compared to the experimental findings. An average error of 140% and 116% was demonstrated between the experimental results and those obtained from the First-Order Upwind and Power Law schemes, respectively. However, both the Second-Order upwind and QUICK schemes accurately predict species concentration under high Peclet number, convection-dominated flow conditions. Convection-diffusion discretisation scheme selection has a strong influence on resultant species concentration fields, as determined by CFD. Furthermore, either the Second-Order or QUICK discretisation schemes should be implemented when numerically modelling convection-dominated mass-transport conditions. Finally, care should be taken not to utilize computationally inexpensive discretisation schemes at the cost of accuracy in resultant species concentration.
Tài liệu tham khảo
Devereux PD, O'Callaghan SM, Walsh MT, McGloughlin TM: Mass transport disturbances in the distal graft/artery junction of a peripheral bypass graft. Proc IMechE Part H: J Engineering in Medicine 2005, 219(6):465–476. 10.1243/095441105X34446
Walsh MT, Kavanagh EG, O'Brien T, Grace PA, McGloughlin T: On the existence of an optimum end-to-side junctional geometry in peripheral bypass surgery - a computer generated study. European Journal of Vascular and Endovascular Surgery 2003, 26: 649–656. 10.1016/j.ejvs.2003.08.004
Ma P, Li X, Ku DN: Convective mass transport at the carotid bifurcation. Journal of Biomechanics 1997, 30: 565–571. 10.1016/S0021-9290(97)84506-X
Giannoglou GD, Soulis JV, Farmakis TM, Farmakis DM, Louridas GE: Haemodynamic factors and the important role of local static pressure in coronary wall thickening. International Journal of Cardiology 2002, 86: 27–40. 10.1016/S0167-5273(02)00188-2
Morris L, Delassus P, Walsh M, McGloughlin T: A mathematical model to predict the in vivo pulsatile drag forces acting on bifurcated stent grafts used in Endovascular treatment of abdominal aortic aneurysms (AAA). Journal of Biomechanics 2004, 37: 1087–1095. 10.1016/j.jbiomech.2003.11.014
Buchanan JR, Kleinstreuer C, Hyun S, Truskey GA: Hemodynamics simulation and identification of susceptible sites of atherosclerotic lesion formation in a model abdominal aorta. Journal of Biomechanics 2003, 36: 1185–1196. 10.1016/S0021-9290(03)00088-5
Lee KW, Xu XY: Modelling of flow and wall behaviour in a mildly stenosed tube. Medical Engineering and Physics 2002, 24: 575–586. 10.1016/S1350-4533(02)00048-6
Lei M, Kleinstreuer C, Truskey GA: A focal stress gradient-dependent mass transfer mechanism for atherogenesis in branching arteries. Medical Engineering and Physics 1996, 18: 326–332. 10.1016/1350-4533(95)00045-3
Prosi M, Perktold K, Ding Z, Friedman MH: Influence of curvature dynamics on pulsatile coronary artery flow in a realistic bifurcation model. Journal of Biomechanics 2004, 37: 1767–1775. 10.1016/j.jbiomech.2004.01.021
Perktold K, Leuprecht A, Prosi M, Berk T, Czerny M, Trubel W, Schima H: Fluid dynamics, wall mechanics, and oxygen transfer in peripheral bypass anastomoses. Annals of Biomedical Engineering 2002, 30: 447–460. 10.1114/1.1477445
Rappitsch G, Perktold K: Pulsatile albumin transport in large arteries: a numerical simulation study. Journal of Biomechanical Engineering 1996, 118: 511–19. 10.1115/1.2796038
Rappitsch G, Perktold K: Computer simulation of convective diffusion processes in large arteries. Journal of Biomechanics 1996, 29: 207–215. 10.1016/0021-9290(95)00045-3
Kaazempur-Mofrad MR, Ethier CR: Mass transport in an anatomically realistic human right coronary artery. Annals of Biomedical Engineering 2001, 29: 121–127. 10.1114/1.1349704
Stewart SFC, Lyman DJ: Effects of artery/vascular graft compliance mismatch on protein transport: A numerical study. Annals of Biomedical Engineering 2004, 32: 991–1006. 10.1023/B:ABME.0000032462.56207.65
Lutostansky EM: The role of convective mass transfer in atherosclerosis. In Ph.D. Thesis. Georgia Institute of Technology, Atlanta, GA, USA; 1996.
Qiu Y, Tarbell JM: Numerical simulation of oxygen mass transfer in a compliant curved tube model of a coronary artery. Annals of Biomedical Engineering 2000, 28: 26–38. 10.1114/1.251
Markou CP, Lutostansky EM, Ku DN, Hanson SR: A novel method for efficient drug delivery. Annals of Biomedical Engineering 1998, 26: 502–511. 10.1114/1.97
O'Brien T, Walsh W, McGloughlin T: On reducing abnormal hemodynamics in the femoral end-to-side anastomosis: The influence of mechanical factors. Annals of Biomedical Engineering 2005, 33: 309–321.
Fluent 6.1 users guide (2005): Spatial Discretization [online]. [http://jullio.pe.kr/fluent6.1/help/html/ug/node814.htm]
Versteeg HK, Malalasekera : An Introduction to Computational Fluid Dynamics. Longman Scientific and Technical; 1995.
Wada S, Karino T: Theoretical prediction of low-density lipoprotein concentration at the luminal surface of an artery with a multiple bend. Annals of Biomedical Engineering 2002, 30: 778–791. 10.1114/1.1495868
Lutostansky EM, Karner G, Rappitsch G, Ku DN, Perktold K: Analysis of hemodynamic fluid phase mass transport in a separated flow region. Journal of Biomechanical Engineering 2003, 125: 189–196. 10.1115/1.1543547
Xue SC, Phan-Thien N, Tanner RI: Upwinding with deferred correction (UPDC): an effective implementation of higher-order convection schemes for implicit finite volume methods. Journal of Non-Newtonian Fluid Mechanics 2002, 108: 1–24. 10.1016/S0377-0257(02)00122-2