Shinjiro Miyawaki

Numerical evaluation of the effects of planform geometry and inflow conditions on flow, turbulence structure, and bed shear velocity at a stream confluence with a concordant bed

This study numerically investigates the effects of variations in inflow conditions and planform geometry on large-scale coherent flow structures and bed friction velocities at a stream confluence with natural bathymetry and concordant bed morphology. Several numerical experiments are conducted in which either the Kelvin-Helmholtz mode or the wake mode dominates within the mixing interface (MI) between the two confluent streams as the junction angle and alignments of the tributaries are altered. In the Kelvin-Helmholtz mode, the MI contains mostly corotating vortices driven by the mean transverse shear across the MI, while in the wake mode the MI contains counterrotating vortices forming by the interaction of the separated shear layers on the two sides of a zone of stagnant fluid near the junction corner. A large angle between the two incoming streams is not necessary for the development of strongly coherent streamwise-oriented vortical (SOV) cells in the immediate vicinity of the MI. Results show that such SOV cells can develop and produce high bed friction velocities even for cases with a low angle between the two tributaries and for cases where the downstream channel is approximately aligned with the axes of the two tributaries (low-curvature cases). SOV cells tend not to develop only when the incoming streams are parallel and aligned with the downstream channel (junction angle of zero), and the incoming flows produce a strong Kelvin-Helmholtz mode. Under such conditions, quasi 2-D MI vortices play the primary role in mixing and the production of high bed shear velocities. Simulations with and without natural bed morphology/local bank line irregularities indicate that planform geometry and inflow conditions primarily govern the development of coherent flow structures, but that bathymetric and bank line effects can locally modify details of these structures.

Effect of static vs. dynamic imaging on particle transport in CT-based numerical models of human central airways.

Advances in quantitative computed tomography (CT) has provided methods to assess the detailed structure of the pulmonary airways and parenchyma, providing the means of applying computational fluid dynamics-based modeling to better understand subject-specific differences in structure-to-function relationships. Most of the previous numerical studies, seeking to predict patterns of inhaled particle deposition, have considered airway geometry and regional ventilation derived from static images. Because geometric alterations of the airway and parenchyma associated with regional ventilation may greatly affect particle transport, we have sought to investigate the effect of rigid vs. deforming airways, linear vs. nonlinear airway deformations, and step-wise static vs. dynamic imaging on particle deposition with varying numbers of intermediate lung volume increments. Airway geometry and regional ventilation at different time points were defined by four-dimensional (space and time) dynamic or static CT images. Laminar, transitional, and turbulent air flows were reproduced with a three-dimensional eddy-resolving computational fluid dynamics model. Finally, trajectories of particles were computed with the Lagrangian tracking algorithm. The results demonstrated that static-imaging-based models can contribute 7% uncertainty to overall particle distribution and deposition primarily due to regional flow rate (ventilation) differences as opposed to geometric alterations. The effect of rigid vs. deforming airways on serial distribution of particles over generations was significantly smaller than reported in a previous study that used the symmetric Weibel geometric model with smaller flow rate. Rigid vs. deforming airways were also shown to affect parallel particle distribution over lobes by 8% and the differences associated with use of static vs. dynamic imaging was 18%. These differences demonstrate that estimates derived from static vs. dynamic imaging can significantly affect the assessment of particle distribution heterogeneity. The effect of linear vs. nonlinear airway deformations was within the uncertainty due to mesh size.

A 4DCT imaging-based breathing lung model with relative hysteresis

To reproduce realistic airway motion and airflow, the authors developed a deforming lung computational fluid dynamics (CFD) model based on four-dimensional (4D, space and time) dynamic computed tomography (CT) images. A total of 13 time points within controlled tidal volume respiration were used to account for realistic and irregular lung motion in human volunteers. Because of the irregular motion of 4DCT-based airways, we identified an optimal interpolation method for airway surface deformation during respiration, and implemented a computational solid mechanics-based moving mesh algorithm to produce smooth deforming airway mesh. In addition, we developed physiologically realistic airflow boundary conditions for both models based on multiple images and a single image. Furthermore, we examined simplified models based on one or two dynamic or static images. By comparing these simplified models with the model based on 13 dynamic images, we investigated the effects of relative hysteresis of lung structure with respect to lung volume, lung deformation, and imaging methods, i.e., dynamic vs. static scans, on CFD-predicted pressure drop. The effect of imaging method on pressure drop was 24 percentage points due to the differences in airflow distribution and airway geometry.

Early Airway Structural Changes in Cystic Fibrosis Pigs as a Determinant of Particle Distribution and Deposition

The pathogenesis of cystic fibrosis (CF) airway disease is not well understood. A porcine CF model was recently generated, and these animals develop lung disease similar to humans with CF. At birth, before infection and inflammation, CF pigs have airways that are irregularly shaped and have a reduced caliber compared to non-CF pigs. We hypothesized that these airway structural abnormalities affect airflow patterns and particle distribution. To test this hypothesis we used computational fluid dynamics (CFD) on airway geometries obtained by computed tomography of newborn non-CF and CF pigs. For the same flow rate, newborn CF pig airways exhibited higher air velocity and resistance compared to non-CF. Moreover we found that, at the carina bifurcation, particles greater than 5-μm preferably distributed to the right CF lung despite almost equal airflow ventilation in non-CF and CF. CFD modeling also predicted that deposition efficiency was greater in CF compared to non-CF for 5- and 10-μm particles. These differences were most significant in the airways included in the geometry supplying the right caudal, right accessory, left caudal, and left cranial lobes. The irregular particle distribution and increased deposition in newborn CF pig airways suggest that early airway structural abnormalities might contribute to CF disease pathogenesis.

A Numerical Study of Water Loss Rate Distributions in MDCT-Based Human Airway Models

Both three-dimensional (3D) and one-dimensional (1D) computational fluid dynamics methods are applied to study regional water loss in three multi-detector row computed-tomography-based human airway models at the minute ventilations of 6, 15 and 30 L/min. The overall water losses predicted by both 3D and 1D models in the entire respiratory tract agree with available experimental measurements. However, 3D and 1D models reveal different regional water loss rate distributions due to the 3D secondary flows formed at bifurcations. The secondary flows cause local skewed temperature and humidity distributions on inspiration acting to elevate the local water loss rate; and the secondary flow at the carina tends to distribute more cold air to the lower lobes. As a result, the 3D model predicts that the water loss rate first increases with increasing airway generation, and then decreases as the air approaches saturation, while the 1D model predicts a monotonic decrease of water loss rate with increasing airway generation. Moreover, the 3D (or 1D) model predicts relatively higher water loss rates in lower (or upper) lobes. The regional water loss rate can be related to the non-dimensional wall shear stress (τ*) by the non-dimensional mass transfer coefficient (h0*) as

A Feasible Computational Fluid Dynamics Study for Relationships of Structural and Functional Alterations with Particle Depositions in Severe Asthmatic Lungs

{h_0}^{*}= 1.15{\tau ^{*0.272}}, R = 0.842

Structure of turbulent flow at a river confluence with momentum and velocity ratios close to 1: Insight provided by an eddy-resolving numerical simulation

h0∗=1.15τ∗0.272,R=0.842.

Numerical analysis of the effect of momentum ratio on the dynamics and sediment‐entrainment capacity of coherent flow structures at a stream confluence

The effects of planform geometry and momentum flux ratio on thermal mixing at a stream confluence with concordant bed morphology are investigated based on numerical simulations that can capture the dynamics of large-scale turbulence. In two simulations, the bathymetry and asymmetrical planform geometry are obtained from field experiments and the momentum flux ratio is set at values of one and four. These two conditions provide the basis for studying differences in thermal mixing processes at this confluence when the wake mode and the Kelvin–Helmholtz mode dominate the development of coherent structures within the mixing interface (MI). The effects of channel curvature and angle between the two incoming streams on thermal mixing processes are investigated based on simulations conducted with modified planform geometries. Two additional simulations are conducted for the case where the upstream channels are parallel but not aligned with the downstream channel and for the zero-curvature case where the upstream channels are parallel and aligned with the downstream channel. The simulations highlight the influence of large-scale coherent structures within the MI and of streamwise-oriented vortical (SOV) cells on thermal mixing processes within the confluence hydrodynamics zone. Simulation results demonstrate the critical role played by the SOV cells in promoting large-scale thermal mixing for cases when such cells form in the immediate vicinity of the MI and in modifying the shape of the thermal MI within cross sections of the downstream channel—predictions consistent with empirical measurements of thermal mixing at the confluence. The set of numerical simulations reveal that the degree of thermal mixing occurring within the confluence hydrodynamic zone varies dramatically with planform geometry and incoming flow conditions. In some cases thermal mixing at the downstream end of the confluence hydrodynamic zone is limited to the MI and its immediate vicinity, whereas in others substantial thermal mixing has occurred over most of the cross-sectional area of the flow. Overall, the simulations highlight the flow conditions and the controls of these conditions that influence mixing within the immediate vicinity of a confluence.