Steady Inspiratory Airflow Research Articles

Computational fluid dynamics calculations in inferior turbinate surgery: a cohort study

PurposeTo investigate how the results of nasal computational fluid dynamics (CFD) simulations change due to inferior turbinate surgery and how the results correlate with patient specific subjective assessment and volumetric results in the nasal cavities.MethodsThe steady inspiratory airflow of 25 patients was studied pre- and postoperatively with heat transfer from the mucous membrane by performing CFD calculations to patient-specific nasal cone beam computed tomography images. These results were then compared to the severity of the patients’ nasal obstruction Visual Analogue Scale (VAS) and Glasgow Health Status Inventory assessments, and acoustic rhinometry measurements.ResultsTotal wall shear forces decreased statistically significantly (p < 0.01) in the operated parts of the inferior turbinates. Patients’ subjective nasal obstruction VAS assessment changes between the pre- and postoperative conditions correlated statistically significantly (p = 0.04) with the wall shear force results.ConclusionInferior turbinate surgery lead to decreased total wall shear force values postoperatively. Changes in subjective nasal obstruction VAS results against total wall shear force changes between the pre- and postoperative conditions were statistically significant. CFD data have a potential to be used for the evaluation of nasal airflow.

Computational fluid dynamics assessed changes of nasal airflow after inferior turbinate surgery

ObjectiveTo demonstrate how Computational Fluid Dynamics (CFD) simulations can reveal important airflow changes in the nasal cavities due to surgical interventions. Material and methodsThe steady inspiratory airflow of eight patients was studied pre- and postoperatively with heat transfer from the mucous membrane by performing CFD calculations to patient specific cone beam computed tomography (CBCT) images. Eight patients with the largest distance from pre- and postoperative mean changes in inferior turbinate volumetry and Visual Analogue Scale (VAS) results were selected. ResultsCalculated CFD heat transfer results from the anterior parts of the inferior turbinates, where surgical interventions were performed, decreased significantly. The heat transfer results were in line with VAS changes. ConclusionSurgical interventions reduced heat transfer in the operated parts of the inferior turbinates and were in line with VAS changes. CFD is an option in assessing patient well-being as a function of airflow parameters from mucous membrane with larger data sets. The limitations of the study were the small sample size and the preliminary nature of the study.

Impact of endoscopic craniofacial resection on simulated nasal airflow and heat transport.

Endoscopic craniofacial resections (CFR) are performed for extensive anterior skull base lesions. This surgery involves removal of multiple intranasal structures, potentially leading to empty nose syndrome (ENS). However, many patients remain asymptomatic postoperatively. Our objective was to analyze the impact of CFR on nasal physiology and airflow using computational fluid dynamics (CFD). This is the first CFD analysis of post-CFR patients. Three-dimensional sinonasal models were constructed from 3 postoperative images using MimicsTM . Hybrid computational meshes were created. Steady inspiratory airflow and heat transport were simulated at patient-specific flow rates using shear stress transport k-omega turbulent flow modeling in FluentTM . Simulated average heat flux (HF) and surface area where HF exceeded 50 W/m2 (SAHF50) were compared with laminar simulations in 9 radiographically normal adults. Three adults underwent CFR without developing ENS. Average HF (W/m2 ) were 132.70, 134.84, and 142.60 in the CFR group, ranging from 156.24 to 234.95 in the nonoperative cohort. SAHF50 (m2 ) values were 0.0087, 0.0120, and 0.0110 in the CFR group, ranging from 0.0082 to 0.0114 in the radiographically normal cohort. SAHF50 was distributed throughout the CFR cavities, with increased HF at the roof and walls compared with the nonoperative cohort. Average HF was low in the CFR group compared with the nonoperative group. However, absence of ENS in most CFR patients may be due to large stimulated mucosal surface area, commensurate with the nonoperative cohort. Diffuse distribution of stimulated area may result from turbulent mixing after CFR. To better understand heat transport post-CFR, a larger cohort is necessary.

Experimental Measurements and Computational Predictions of Regional Particle Deposition in a Sectional Nasal Model

Knowledge of the regional deposition of inhaled particles in the nose is important for drug delivery and assessment of the toxicity of inhaled materials. In this study, computational fluid dynamics (CFD) predictions and experimental measurements in a nasal replica cast were used to study regional deposition of inhaled microparticles. The replica cast was sectioned into six regions of interest based on nasal anatomy: the nasal vestibule, nasal valve, anterior turbinates, olfactory region, turbinates, and nasopharynx. Monodisperse fluorescein particles with aerodynamic diameters of 2.6-14.3 μm were passed through the assembled cast in the presence of steady inspiratory airflow at 15 L/min. After each experiment, the cast was disassembled and the deposited fluorescein in each region was washed out and quantified with fluorescence spectrometry. A nasal CFD model was developed from the same magnetic resonance imaging scans that were used to construct the replica cast. Steady-state inspiratory airflow and particle deposition calculations were conducted in the CFD model using Fluent(™) at flow rates producing Stokes numbers comparable to experimental conditions. Total and regional particle deposition predictions from the CFD model were compared with experimental measurements from the replica cast. Overall, good agreement was observed between CFD predictions and experimental measurements with similar deposition trends in each region of interest. CFD predictions in central nasal regions demonstrated well-defined maximum values of 15%, 7%, and 12% in the anterior turbinates, olfactory, and turbinates regions, respectively, at particle sizes of 10-11 μm. These results demonstrate the use of a sectioned nasal CFD model based on anatomical regions of interest for nasal drug delivery to elucidate patterns of regional deposition within a human nasal cavity.

Simulation of Enhanced Deposition Due to Magnetic Field Alignment of Ellipsoidal Particles in a Lung Bifurcation

Aerosolized chemotherapy has been recognized as a potential treatment for lung cancer. The challenge of providing sufficient therapeutic effects without reaching dose-limiting toxicity levels hinders the development of aerosolized chemotherapy. This could be mitigated by increasing drug-delivery efficiency with a noninvasive drug-targeting delivery method. The purpose of this study is to use direct numerical simulations to study the resulting local enhancement of deposition due to magnetic field alignment of high aspect ratio particles. High aspect ratio particles were approximated by a rigid ellipsoid with a minor diameter of 0.5 μm and fluid particle density ratio of 1,000. Particle trajectories were calculated by solving the coupled fluid particle equations using an in-house micro-macro grid finite element algorithm based on a previously developed fictitious domain approach. Particle trajectories were simulated in a morphologically realistic geometry modeling a symmetrical terminal bronchiole bifurcation. Flow conditions were steady inspiratory air flow due to typical breathing at 18 L/min. Deposition efficiency was estimated for two different cases: [1] particles aligned with the streamlines and [2] particles with fixed angular orientation simulating the magnetic field alignment of our previous in vitro study. The local enhancement factor defined as the ratio between deposition efficiency of Case [1] and Case [2] was found to be 1.43 and 3.46 for particles with an aspect ratio of 6 and 20, respectively. Results indicate that externally forcing local alignment of high aspect ratio particles can increase local deposition considerably.

Laminar flow and wall shear stress in a collapsed tube

Detailed description of the flow field in human airways is highly important to better understand human breathing and provide a patient's customized diagnosis. An integrated numerical simulation platform is presently proposed in order to incorporate medical images into a numerical software to calculate flow field and to analyze it in terms of fluid dynamics. The platform was set up to compute steady inspiratory airflow in realistic human airways reconstructed from tomodensitometric medical images at resting breathing conditions. This morpho-functional simulation platform has been tested retrospectively with two CT-scanned patient airway morphological models: (i) a normal airway model (subject A) with no evidence of morphological alteration and (ii) a highly altered airway model (subject B) exhibiting a severe stenosis in the right main bronchus. First, various morphological aspects proper to each airway model are provided to show the performance and interest of the reconstruction method. Second, we describe the three-dimensional flow patterns associated to the global morphological features, which are mainly shared by the present realistic models and previous idealistic airway models. Finally, the flow characteristics associated to local morphological features specific to realistic airway models are discussed. The results demonstrate that the morpho-functional simulation platform is able to capture the main features of airway velocity patterns but also more specific airflow patterns which are related to customized patient morphological features such as laminar vortex formation. The present results suggest that the proposed airway functional imaging platform is adequate to provide most of functional information related to airflow and enable a patient to patient diagnosis.

Airflow modeling of steady inspiration in two realistic proximal airway trees reconstructed from human thoracic tomodensitometric images

Detailed description of the flow field in human airways is highly important to better understand human breathing and provide a patient's customized diagnosis. An integrated numerical simulation platform is presently proposed in order to incorporate medical images into a numerical software to calculate flow field and to analyze it in terms of fluid dynamics. The platform was set up to compute steady inspiratory airflow in realistic human airways reconstructed from tomodensitometric medical images at resting breathing conditions. This morpho-functional simulation platform has been tested retrospectively with two CT-scanned patient airway morphological models: (i) a normal airway model (subject A) with no evidence of morphological alteration and (ii) a highly altered airway model (subject B) exhibiting a severe stenosis in the right main bronchus. First, various morphological aspects proper to each airway model are provided to show the performance and interest of the reconstruction method. Second, we describe the three-dimensional flow patterns associated to the global morphological features, which are mainly shared by the present realistic models and previous idealistic airway models. Finally, the flow characteristics associated to local morphological features specific to realistic airway models are discussed. The results demonstrate that the morpho-functional simulation platform is able to capture the main features of airway velocity patterns but also more specific airflow patterns which are related to customized patient morphological features such as laminar vortex formation. The present results suggest that the proposed airway functional imaging platform is adequate to provide most of functional information related to airflow and enable a patient to patient diagnosis.

Finite element analysis of particle motion in steady inspiratory airflow

To accurately model the inhaled particle motion, equations governing particle trajectories in carrier flow are solved together with the Navier–Stokes equations. Under the relatively dilute particle condition in the mixture, equations for two phases are coupled through the interface drag shown in the solid-phase momentum equations. The present study investigates bifurcation flow in the human central airway using the finite element method. In the gas phase, we employ the biquadratic streamline upwind Petrov–Galerkin finite element model to simulate the incompressible air flow. To solve the equations of motion for the inhaled particles, we apply another biquadratic streamline upwind finite element model. A feature common to two models applied to each phase of equations is that both of them provide nodally exact solutions to the convection–diffusion and the convection–reaction equations, which are prototype equations for the gas-phase and the solid-phase equations, respectively. In two dimensions, both models have ability to introduce physically meaningful artificial damping terms solely in the streamline direction. With these terms added to the formulation, the discrete system is enhanced without compromising the numerical diffusion error. Tests on inspiratory problem were conducted, and the results are presented, with an emphasis on the discussion of particle motion.

Deposition of monodisperse particles in hollow models representing adult and child-size tracheobronchial airways.

A series of experiments was performed to determine deposition efficiencies of four sizes of radiolabeled monodisperse particles in custom-made hollow tracheobronchial models. The particles had geometric diameters of about 1, 5, 10, and 15 microm. The tracheobronchial models, consisting of a trachea and two or more additional generations, had dimensions representative of a typical adult, a 7-y-old child, and a 4-y-old child; the child-size models were appropriately scaled-down replicas of the adult-size model. Each deposition experiment was conducted using a steady inspiratory airflow representative of low physical activity for the appropriate age of individual: 20 L min(-1) for the adult; 9 L min(-1) for the 7-y-old; 5.5 L min(-1) for the 4-y-old. The results indicate that deposition efficiency of the particles increased substantially (up to 35 times) in all three models as particle diameter increased from 1-15 microm, undoubtedly as a result of particle impaction and sedimentation-related phenomena. An analysis of variance demonstrated the occurrence of statistically-significant (p < 0.05) main effects of hollow model size and particle size; the interaction between the two parameters was also significant. In general, deposition efficiencies of the various sizes of particles were greater in the child-size models than in the adult-size model; this effect may have risk assessment implications. In addition, the results obtained experimentally agreed more closely with those predicted using a radiation-protection mathematical particle deposition formulation as the particle size increased for each of the sizes of models.

Numerical modeling of steady inspiratory airflow through a three-generation model of the human central airways.

Two-dimensional steady inspiratory airflow through a three-generation model of the human central airways is numerically investigated, with dimensions corresponding to those encountered in the fifth to seventh generations of the Weibel's model. Wall curvatures are added at the outer walls of the junctions for physiological purposes. Computations are carried out for Reynolds numbers in the mother branch ranging from 200 to 1200, which correspond to mouth air breathing at flow rates ranging from 0.27 to 1.63 liters per second. The difficulty of generating grids in a so complex configuration is overcome using a nonoverlapping multiblock technique. Simulations demonstrate the existence of separation regions whose number, location, and size strongly depend on the Reynolds number. Consequently, four different flow configurations are detected. Velocity profiles downstream of the bifurcations are shown to be highly skewed, thus leading to an important unbalance in the flow distribution between the medial and lateral branches of the model. These results confirm the observations of Snyder et al. and Tsuda et al. and suggest that a resistance model of flow partitioning based on Kirchhoff's laws is inadequate to simulate the flow behavior accurately within the airways. When plotted in a Moody diagram, airway resistance throughout the model is shown to fit with a linear relation of slope -0.61. This is qualitatively in good agreement with the experimental investigations of Pedley et al, and Slutsky et al.

ICIAM/GAMM 95 Applied sciences, especially Mechanics Contributed Lecture and Poster Presentations

The present paper deals with 2D numerical modelling of steady inspiratory airflow through a three-generations model of the human central airways, whose dimensions correspond to those encountered in the 5th to 7th generations of the Weibel's model. Computations have been carried out for Reynolds numbers in the mother branch ranging from 200 to 1200 in the case of forced ventilation. Simulations demonstrate the existence of separation regions whose number, location and size strongly depend on the Reynolds number. Airway resistance through the model is shown to be qualitatively in good agreement with available experimental investigations.