Jiyuan Tu

Numerical simulations for detailed airflow dynamics in a human nasal cavity

Nasal physiology is dependent on the physical structure of the nose. Individual aspects of the nasal cavity such as the geometry and flow rate collectively affect nasal function such as the filtration of foreign particles by bringing inspired air into contact with mucous-coated walls, humidifying and warming the air before it enters the lungs and the sense of smell. To better understand the physiology of the nose, this study makes use of CFD methods and post-processing techniques to present flow patterns between the left and right nasal cavities and compared the results with experimental and numerical data that are available in literature. The CFD simulation adopted a laminar steady flow for flow rates of 7.5 L/min and 15 L/min. General agreement of gross flow features were found that included high velocities in the constrictive nasal valve area region, high flow close to the septum walls, and vortex formations posterior to the nasal valve and olfactory regions. The differences in the left and right cavities were explored and the effects it had on the flow field were discussed especially in the nasal valve and middle turbinate regions. Geometrical differences were also compared with available models.

A numerical study of spray particle deposition in a human nasal cavity

Particle depositional studies from nasal sprays are important for efficient drug delivery. The main influences on deposition involve the nasal cavity geometry and the nasal spray device of which its parameters are controlled by the product design. It is known that larger particle sizes (≫ 10 μm) at a flow rate of 333 ml/s impact in the anterior portion of the nose, leaving a significant portion of the nasal cavity unexposed to the drugs. Studies have found correlations for the spray cone angles and particle sizes with deposition efficiencies. This study extends these ideas to incorporate other parameters such as the insertion angle of the nasal spray and the injected particle velocity to observe its effect on deposition. A numerical method utilizing a particle tracking procedure found that the most important parameter was the particles Stokes number which affected all other parameters on the deposition efficiency.

Comparison of micron- and nanoparticle deposition patterns in a realistic human nasal cavity

Knowledge regarding particle deposition processes in the nasal cavity is important in aerosol therapy and inhalation toxicology applications. This paper presents a comparative study of the deposition of micron and submicron particles under different steady laminar flow rates using a Lagrangian approach. A computational model of a nasal cavity geometry was developed from CT scans and the simulation of the fluid and particle flow within the airway was performed using the commercial software GAMBIT and FLUENT. The air flow patterns in the nasal cavities and the detailed local deposition patterns of micron and submicron particles were presented and discussed. It was found that the majority of micron particles are deposited near the nasal valve region and some micron particles are deposited on the septum wall in the turbinate region. The deposition patterns of micron particles in the left cavity are different compared with that in the right one especially in the turbinate regions. In contrast, the deposition for nanoparticles shows a moderately even distribution of particles throughout the airway. Furthermore the particles releasing position obviously influences the local deposition patterns. The influence of the particle releasing position is mainly shown near the nasal valve region for micron particle deposition, while for submicron particles deposition, both the nasal valve and turbinate region are influenced. The results of the paper are valuable in aerosol therapy and inhalation toxicology.

Gas-liquid flows

Gas–liquid flows appear in natural and industrial processes in various forms and often feature complex inter-phase mass, momentum, and energy transfers. One example of naturally occurring gas–liquid flow is the dispersion of marine droplets. Gas–liquid flows are also found in abundance in industrial processes. One significant industrial application is venting of mixture vapors to liquid pools in chemical reactors. This chapter is concerned with gas–liquid flows. Within this flow system, the two phases that coexist simultaneously in the fluid flow often exhibit relative motion among the phases and heat and mass exchanges across the interface boundary. Owing to the complexities of interfaces and resultant discontinuities in fluid properties as well as from physical scaling issues, it is rather customary to apply a statistical, averaged approach in the form of a two-fluid model to resolve such a flow system. Separate transport equations governing the conservation of mass, momentum, and energy are solved for each phase and exchanges that take place at the interfaces between the two phases are explicitly accounted for in which the dynamics of the interaction between the two phases can be effectively described via suitable models of the inter-phase mass, momentum, and energy exchanges. Normally, the coupling between the two phases is very tight, which demands special numerical strategies and solution algorithms to be adopted. This particular flow system is also complicated considerably by the prevalence of particle–particle collisions. A number of population balance methods, along with suitable coalescence and break-up mechanisms, are discussed. In the context of computational fluid dynamics, the application of population balance models to describe the coalescence and break-up dynamics of these gas particles can be coupled with the two-fluid model to predict the wide range of particle sizes within the two-phase flow.

Effect of calcification on the mechanical stability of plaque based on a three-dimensional carotid bifurcation model

BackgroundThis study characterizes the distribution and components of plaque structure by presenting a three-dimensional blood-vessel modelling with the aim of determining mechanical properties due to the effect of lipid core and calcification within a plaque. Numerical simulation has been used to answer how cap thickness and calcium distribution in lipids influence the biomechanical stress on the plaque.MethodModelling atherosclerotic plaque based on structural analysis confirms the rationale for plaque mechanical examination and the feasibility of our simulation model. Meaningful validation of predictions from modelled atherosclerotic plaque model typically requires examination of bona fide atherosclerotic lesions. To analyze a more accurate plaque rupture, fluid-structure interaction is applied to three-dimensional blood-vessel carotid bifurcation modelling. A patient-specific pressure variation is applied onto the plaque to influence its vulnerability.ResultsModelling of the human atherosclerotic artery with varying degrees of lipid core elasticity, fibrous cap thickness and calcification gap, which is defined as the distance between the fibrous cap and calcification agglomerate, form the basis of our rupture analysis. Finite element analysis shows that the calcification gap should be conservatively smaller than its threshold to maintain plaque stability. The results add new mechanistic insights and methodologically sound data to investigate plaque rupture mechanics.ConclusionStructural analysis using a three-dimensional calcified model represents a more realistic simulation of late-stage atherosclerotic plaque. We also demonstrate that increases of calcium content that is coupled with a decrease in lipid core volume can stabilize plaque structurally.

Computational fluid and particle dynamics in the human respiratory system

From the Contents: Computational Fluid Particle Dynamics (CFPD) - An Introduction: What is CFPD.- The Human Respiratory System: Introduction.- Anatomy of the respiratory system.- Reconstruction of the Human Airways: Introduction.- Medical image acquisition.- Generation of Computational Mesh for CFPD Simulation: Introduction.- Mesh types.- Fundamentals of Fluid Dynamics: Introduction.- Fluid dynamics and governing equations.- Fundamentals of Particle Dynamics: Particle dynamics and mathematical models.- Particle trajectory models.- Continuum approach.- Modelling of further particle physics.- Basic Computational Methods: Introduction.- Case studies in the human airways: Introduction.- Modelling inhalation and heat transfer in the nasal cavity.- Inhalation of toxic particles and the effects of particle morphology.- Optimisation of nasal drug delivery.- Advanced Topics and Future Trends: Moving and Deforming Mesh.- Fluid-Structure Interaction.

Numerical analysis of micro- and nano-particle deposition in a realistic human upper airway

A computational model was developed for studying the flow field and particle deposition in a human upper airway system, including: nasal cavity, nasopharynx, oropharynx, larynx and trachea. A series of coronal CT scan images of a 24 year old woman was used to construct the 3D model. The Lagrangian and Eulerian approaches were used, respectively, to find the trajectories of micro-particles and concentration of nano-particles. The total and regional deposition fractions of micro/nanoparticles were evaluated and the major hot spots for the deposition of inhaled particles were found.

Inhalability of micron particles through the nose and mouth

Aspiration efficiencies from nose and mouth inhalations are investigated at low and high inhalation rates by using the commercial Computational Fluid Dynamics (CFD) software CFX 11. A realistic human head with detailed facial features was constructed. Facial features were matched to represent the 50th percentile of a human male, aged between 20 and 65 years old, based on anthropometric data. The constant freestream velocity was 0.2 ms−1, normal to the face, and inhalation rates through the mouth and nose were 15 liters per minute (LPM) for light breathing and 40 LPM for heavy breathing. It was found that the flow field in the near breathing region exhibited vertical direction caused by the presence of the torso where the airstream diverges as it flows around and over the body. The critical area concept was used as a tool to determine the aspiration efficiency of particles. Comparisons between critical areas for the nose and mouth inhalations show similar geometric properties such as the area’s shape, and its vertical distance location on the x-z plane located at y = 80 cm upstream. The critical area sizes were found to be slightly larger for the mouth inhalation mainly due to the larger mouth area and also the aligned orientation of the mouth to the upstream flow, whereas the nose is perpendicular to the upstream flow. This study was undertaken to establish the flow field in the near breathing region that will help to characterize the flow and particle field for initial boundary conditions leading to a more holistic modeling approach of respiration through the internal nasal cavity and mouth.

Optimising nasal spray parameters for efficient drug delivery using computational fluid dynamics

Experimental images from particle/droplet image analyser (PDIA) and particle image velocimetry (PIV) imaging techniques of particle formation from a nasal spray device were taken to determine critical parameters for the study and design of effective nasal drug delivery devices. The critical parameters found were particle size, diameter of spray cone at a break-up length and a spray cone angle. A range of values for each of the parameters were ascertained through imaging analysis which were then transposed into initial particle boundary conditions for particle flow simulation within the nasal cavity by using Computational Fluid Dynamics software. An Eulerian-Lagrangian scheme was utilised to track mono-dispersed particles (10 and 20 microm) at a breathing rate of 10 L/min. The results from this qualitative study aim to assist the pharmaceutical industry to improve and help guide the design of nasal spray devices.

Biomechanical investigation of pulsatile flow in a three-dimensional atherosclerotic carotid bifurcation model

It is a well-established fact that atherosclerosis in carotid bifurcation depends on flow parameters such as wall shear stress, flow pulsatility, and blood pressure. However, it is still not clearly verified how atherosclerosis can become aggravated when plaque experiences a high level of shear stress during advance stages of this disease. In this paper, fluid and structural properties in idealistic geometries are analyzed by using fluid-structure interaction (FSI). From our results, the relationship among blood pressure, stenotic compression, and deformation was established. We show that a high level of compression occurs at the stenotic apex, and can potentially be responsible for plaque progression. Moreover, wall shear stress and deformation are significantly affected by the degree of stenosis. Finally, through analysis of the FSI-based simulation results, we can better understand the parameters that influence flow through a stenotic artery and plaque aggravation, and apply the knowledge for the enhancement of clinical research and prediction of treatment outcomes.