Lagrangian Particle Transport Research Articles

Model of multiphase flow, superheat transport, and particle capture in a steel continuous caster

A model of turbulent multiphase flow, heat transfer, and particle entrapment during continuous casting of steel is presented. The model includes the top 7m of the vertical and curved strand, and considers the effects of argon gas injection and thermal buoyancy. RANS model flow results are compared with LES. Lagrangian particle transport through the Eulerian-Eulerian multiphase flow field from the k-ɛ RANS model is based on random walk and features advanced particle capture criteria, improved from previous work by including anisotropic turbulent velocity fluctuations near the walls. Particle capture via 3 mechanisms is included: capture by solidified hooks at the meniscus, entrapment between dendrites, and engulfment by surrounding large particles. The fluid flow and bubble/inclusion capture results are validated with plant measurements, including nail board dipping tests and ultrasonic tests of particle locations, and good agreement is seen. The superheat has negligible effect on flow in the mold region but causes complex flow in the lower strand by creating multiple recirculation zones due to the thermal buoyancy. With high (30 K) superheat, this leads to less penetration, and slightly fewer and shallower capture of particles. Lower (10 K) superheat may enable significant top surface freezing, leading to very large internal defect clusters. Lower superheat also leads to deeper meniscus hooks, and more surface capture. Capture bands occur near the transition from vertical to curved, where downward liquid flow velocity balances the particle terminal velocity.

Численное моделирование переноса твёрдых частиц в атмосферном городском пограничном слое с использованием лагранжева подхода: физические задачи и параллельная реализация

This paper presents the results of the development of a numerical model of the Lagrangian particle transport and the application of parallel computation methods to increase the efficiency of the software implementation of the model. The model is a software package allowing calculations of transport and deposition of aerosol particles taking into account the properties of particles and input data describing atmospheric conditions and the underlying surface geometry. The dynamic core, physical parameterizations, numerical implementation, and algorithm of the model are described. Initially, the model has been used for computationally low-intensive problems. In this paper, given the need to use the model in computationally intensive problems, we conduct optimization of the sequential software implementation of the model, as well as creation of software implementations of the model with the use of parallel computing technologies OpenMP, MPI, CUDA. The results of testing of different implementations of the model show that optimization of the most computationally complex blocks in the sequential version of the model can reduce the execution time by 27%, at the same time the use of parallel computing technologies allows to achieve acceleration by several orders of magnitude. The use of OpenMP in dynamic block of the model resulted in acceleration of block up to 4 times, the use of MPI – up to 8 times, the use of CUDA – up to 16 times, all other conditions being equal. Recommendations on the choice of parallel computing technology depending on the properties of the computing system are proposed.

Lagrangian conditional statistics and flow topology in edge plasma turbulence

Lagrangian statistics and particle transport in edge plasma turbulence are investigated using the Hasegawa–Wakatani model and its modified version. The latter shows the emergence of pronounced zonal flows. Different values of the adiabaticity parameter are considered. The main goal is to characterize the role of coherent structures, i.e., vortices and zonal flows, and their impact on the Lagrangian statistics of particles. Computationally intensive long time simulations following ensembles of test particles over hundreds of eddy turnover times are considered in statistically stationary turbulent flows. The flow topology is characterized using the Lagrangian Okubo–Weiss criterion in order to split the flow into topologically different domains. In elliptic and hyperbolic regions, the probability density functions (PDFs) of the residence time have self-similar algebraic decaying tails. However, in the intermediate regions, the PDFs exhibit exponentially decaying tails. Topologically conditioned PDFs of the Lagrangian velocity, and acceleration and density fluctuations are likewise computed. The differences between the classical Hasegawa–Wakatani system and its modified version are assessed, and the role of zonal flows is highlighted. The density flux spectrum, which characterizes the contributions of different length scales, is studied, and its inertial scaling is found to be in agreement with predictions based on dimensional arguments. Analyzing the angular change of particle tracers at different time scales, corresponding to coarse grained curvature, completes the study, and these multiscale geometric statistics quantify the directional properties of the particle motion in different flow regimes.

Mesoscale Eddies in the Black Sea and Their Impact on River Plumes: Numerical Modeling and Satellite Observations

The Northeast Caucasian Current (NCC) is the northeastern part of the cyclonic Rim Current (RC) in the Black Sea. As it sometimes approaches the narrow shelf very closely, topographically generated cyclonic eddies (TGEs) can be triggered. These eddies contribute to intense, along- and cross-shelf transport of trapped water with enhanced self-cleaning effects of the coastal zone. Despite intense studies of eddy dynamics in the Black Sea, the mechanisms of the generation of such coastal eddies, their unpredictability, and their capacity to capture and transport impurities are still poorly understood. We applied a 3-D low-dissipation model DieCAST/Die2BS coupled with a Lagrangian particle transport model supported by analysis of optical satellite images to study generation and evolution of TGEs and their effect on river plumes unevenly distributed along the northeastern Caucasian coast. Using the Furrier and wavelet analyses of kinetic energy time series, it was revealed that the occurrence of mesoscale TGEs ranges from 10 up to 50 days. We focused on one particular isolated anticyclonic TGE that emerged in late fall as a result of instability of the RC impinging on the abrupt submarine area adjoining the Pitsunda and Iskuria capes. Being shed, the eddy with a 30-km radius traveled along the coast as a coherent structure during ~1.5 months at a velocity of ~3 km/day and vertical vorticity normalized by the Coriolis parameter ~(0.1 ÷ 1.2). This eddy captured water from river plumes localized along the coast and then ejected it to the open sea, providing an intense cross-shelf transport of riverine matter.

Hydrometeorological assessments of the transport of microplastic pellets in the Eastern Mediterranean

Microplastic pellets were sampled in May and November 2018 during one-week surveys at 13 coastal beaches in Iskenderun Bay/Turkey. Pellet pollution index (PPI) was calculated for the beaches as a tool to assess beach pollution by microplastic pellets. Hydrometeorological conditions, including wind, current, wave, surface run-off, and precipitation, were examined during 2018 to reveal the effect on the transport of microplastic pellets within the study area. Sea-surface heights, including the astronomical tide and the storm surge and the wave runup heights, were also considered in the analysis to study the extent of hydrodynamic forcing on the beach. Hydrometeorological assessments indicated that the pellet concentrations in the coastal zone are mostly related to wind-induced transport. Three major river discharges are considered as the main source of microplastic pellets effluents. A Lagrangian particle transport model was conducted to reveal the possible beaching hotspots of microplastic pellets released from these river mouths. Average microplastic pellets were calculated as 126.04 ± 54.08 items/m2 for May 2018 and 70.22 ± 18.25 items/m2 for November 2018. An overall mean PPI for May 2018 was calculated as 1.13, indicating a moderate degree of pellet pollution, and 0.56 for November 2018, indicating a low degree of pellet pollution. The simulations showed that Orontes River effluents affected the inner Iskenderun Bay coasts more than the Seyhan and Ceyhan River.

Three-dimensional direct numerical simulation of flow induced by an oscillating sphere close to a plane boundary

The dynamic interaction of submerged spherical structures with ocean waves alters the local flow field and is responsible for the movement of sediments over the seabed, which would have detrimental effects on the ecosystems. The present research aims to gain insight into the effect of a plane boundary on the vorticity dynamics and shear modifications on the plane boundary induced by the flow due to an oscillating sphere. The flow field around an oscillating sphere close to a plane boundary is investigated numerically by three-dimensional direct numerical simulation of the Navier–Stokes equations. The sphere oscillates perpendicular to the plane boundary modeling the seabed. The flow structure, resulting mean and fluctuating shear stress distribution on the plane boundary, and the Lagrangian particle transport are investigated as a function of the sphere oscillation amplitude, Reynolds number, and the distance of the sphere from the plane boundary. The systematic correlations for the variation of maximum shear stress on the plane boundary are estimated from scaling parameters. The time-mean results are then linked to the Lagrangian particle transport where it is observed that the Lagrangian path lines track the mean streamlines closely. Further, the particles at the stagnation zone cease drifting and oscillate about a constant mean position. The practical significance of this result is that a ring of zero mean shear on the seabed is found toward where particles congregate. This provides a guide for the optimum location for the mooring anchors, minimizing the risk of anchor scour or undermining.

Global simulations of marine plastic transport show plastic trapping in coastal zones

Global coastlines potentially contain significant amounts of plastic debris, with harmful implications for marine and coastal ecosystems, fisheries and tourism. However, the global amount, distribution and origin of plastic debris on beaches and in coastal waters is currently unknown. Here we analyze beaching and resuspension scenarios using a Lagrangian particle transport model. Throughout the first 5 years after entering the ocean, the model indicates that at least 77% of positively buoyant marine plastic debris (PBMPD) released from land-based sources is either beached or floating in coastal waters, assuming no further plastic removal from beaches or the ocean surface. The highest concentrations of beached PBMPD are found in Southeast Asia, caused by high plastic inputs from land and limited offshore transport, although the absolute concentrations are generally overestimates compared to field measurements. The modeled distribution on a global scale is only weakly influenced by local variations in resuspension rates due to coastal geomorphology. Furthermore, there are striking differences regarding the origin of the beached plastic debris. In some exclusive economic zones (EEZ), such as the Indonesian Archipelago, plastic originates almost entirely from within the EEZ while in other EEZs, particularly remote islands, almost all beached plastic debris arrives from remote sources. Our results highlight coastlines and coastal waters as important reservoirs of marine plastic debris and limited transport of PBMPD between the coastal zone and the open ocean.

Rapid assessment tool for oil spill planning and contingencies

The Rapid Oil Spill Hazard Assessment is presented as a demonstration of concept for a tool providing a framework for managers and planners to assess potential impact areas of oil spills. The tool consists of precomputed oil spill scenarios derived from the analysis of twenty years of modeled current data using Self-Organizing Maps to identify 16 representative patterns. These patterns were used to provide boundary conditions for hydrodynamic and wave models to generate higher resolution current fields, used to drive a Lagrangian oil particle transport model creating the most probable oil spill dispersion patterns. To demonstrate the concept, the tool is applied to the Perdido region in the western Gulf of Mexico. A total of 896 oil spill simulations were performed, considering surface and bottom spills, and were stored in a database for easy access to map arrival probabilities and times to be used in risk and vulnerability analysis.

Resolution dependency of sinking Lagrangian particles in ocean general circulation models.

Any type of non-buoyant material in the ocean is transported horizontally by currents during its sinking journey. This lateral transport can be far from negligible for small sinking velocities. To estimate its magnitude and direction, the material is often modelled as a set of Lagrangian particles advected by current velocities that are obtained from Ocean General Circulation Models (OGCMs). State-of-the-art OGCMs are strongly eddying, similar to the real ocean, providing results with a spatial resolution on the order of 10 km on a daily frequency. While the importance of eddies in OGCMs is well-appreciated in the physical oceanographic community, other marine research communities may not. Further, many long term climate modelling simulations (e.g. in paleoclimate) rely on lower spatial resolution models that do not capture mesoscale features. To demonstrate how much the absence of mesoscale features in low-resolution models influences the Lagrangian particle transport, we simulate the transport of sinking Lagrangian particles using low- and high-resolution global OGCMs, and assess the lateral transport differences resulting from the difference in spatial and temporal model resolution. We find major differences between the transport in the non-eddying OGCM and in the eddying OGCM. Addition of stochastic noise to the particle trajectories in the non-eddying OGCM parameterises the effect of eddies well in some cases (e.g. in the North Pacific gyre). The effect of a coarser temporal resolution (once every 5 days versus monthly) is smaller compared to a coarser spatial resolution (0.1° versus 1° horizontally). We recommend to use sinking Lagrangian particles, representing e.g. marine snow, microplankton or sinking plastic, only with velocity fields from eddying Eulerian OGCMs, requiring high-resolution models in e.g. paleoceanographic studies. To increase the accessibility of our particle trace simulations, we launch planktondrift.science.uu.nl, an online tool to reconstruct the surface origin of sedimentary particles in a specific location.

A filtering strategy for the numerical convergence of radiation transport through purely absorbing particle clouds

Radiation transport through particle clouds plays a major role in many engineering applications. In this paper, we study this problem by means of solving the radiative transport equation (RTE) on a Eulerian continuous domain. The particle clouds are generated through direct numerical simulation of Navier-Stokes equations coupled with Lagrangian particle transport at different Stokes numbers. We confirm the earlier observation noted in the literature that the solution to the RTE on Eulerian mesh diverges when the Eulerian mesh size is of the order of the particle diameter. This observation is often called the homogenization error stemming from relegation of number density onto the Eulerian domain. In order to circumvent this divergence problem, we propose a filtering strategy that spreads the information of the particles to the neighboring cells in a way that the representation of the particles remains the same even when the mesh size is smaller than that of the particle size. We show through our simulations that our filtering strategy solves the issue of homogenization error for the cases where the Beer-Bouger law is valid.

Unstructured continuum–kinetic solver with fluid–particle interaction for in-space propulsion contaminant dispersal

A new parallel generalized unstructured grid solver for prediction of rarefied gas dynamics is presented. The method is based on the direct numerical solution of the semi-classical Boltzmann equation using a discrete velocity method (DVM) and supports accurate prediction of non-continuum flows in near-vacuum rarefied flow regimes. For in-space propulsion applications which involve mixed continuum–rarefied environments, the flow regimes are solved separately through the use of DVM and hybrid RANS/LES solution procedures, respectively, on the same computational domain while interacting with one another. The solver also supports coupled Lagrangian particle transport and coalescence/breakup modeling to enable prediction of plume flow impingement and contaminant dispersal through mixed continuum–rarefied environments around spacecraft. Such simulation capabilities are essential for the prediction of spacecraft environments concerning plume impingement pressure, heating, and contamination on spacecraft components such as solar panels and sensitive instruments. Example applications include venting and plume flow environments from spacecraft main engines and attitude control thrusters firing in near-vacuum conditions. The framework upon which the solvers are developed is introduced along with a detailed description of the direct Boltzmann solution and collision integral evaluation procedures. The accuracy of standalone kinetic and coupled continuum–kinetic predictions is validated by comparison against experimental data for shock structure using various different intermolecular interaction models. Results are presented to demonstrate the merits of the capability for prediction of high-speed gas–particle interaction through complex moving shock systems, and thruster plume flow impingement with liquid droplet tracking in near-vacuum conditions.

Ash Metrics for European and Trans‐Atlantic Air Routes During the Eyjafjallajökull Eruption 14 April to 23 May 2010

AbstractMetrics for the risk associated with the threat that airborne volcanic ash particles pose to commercial jet aircraft are presented using simulations based on a Lagrangian particle transport and dispersion model driven by satellite measurements for the Eyjafjallajökull volcanic eruption, Iceland, for the period 14 April to 23 May 2010. The study utilizes a four‐dimensional data set of simulated ash concentrations together with European and trans‐Atlantic air routes to determine metrics corresponding to the total mass intercepted (defined as the dose), the mass interception rate (the dose rate), and the concentration (the exposure) over time (the dosage) that a jet aircraft encounters along the air route. The methodology can be used as a logistical and flight planning tool in a forecast mode and also in hindcast mode to assess the extent of airline fleet exposure to ash following an eruption, thereby providing operators with information useful for flight safety.

CFD simulation of pesticide spray from air-assisted sprayers in an apple orchard: Tree deposition and off-target losses

The ultimate goal of a pesticide spraying system is to provide adequate coverage on intended canopies with a minimum amount of spray materials and off-target waste. Better spray coverage requires an understanding of the fate and transport of spray droplets carried by turbulent airflows in orchards. In this study, an integrated computational fluid dynamics (CFD) model was developed to predict displacement of pesticide spray droplets discharged from an air-assisted sprayer, depositions onto tree canopies, and off-target deposition and airborne drift in an apple orchard. Pesticide droplets discharged from a moving sprayer were tracked using the Lagrangian particle transport model, and the deposition model was applied to droplets entering porous canopy zones. Measurements of the droplet deposition and drift in the same orchard were used to validate the model simulations. Good agreement was found between the measured and simulated spray concentrations inside tree canopies and off-target losses (ground deposition and airborne drifts) with the overall relative errors of 22.1% and 40.6%, respectively, under three growth stages. The CFD model was able to estimate the mass balance of pesticide droplets in the orchard, which was practically difficult to investigate by measurements in field conditions. As the foliage of trees became denser, spray deposition inside canopies increased from 8.5% to 65.8% and airborne drift and ground deposition decreased from 25.8% to 7.0% and 47.8% to 21.2%, respectively. Higher wind speed also increased the spray airborne drift downwind of the orchard. This study demonstrates that CFD model can be used to evaluate spray application performance and design and operate sprayers with increased spray efficiencies and reduced drift potentials.

Quantifying oceanic moisture exports to mainland China in association with summer precipitation

Oceanic moisture exports (OMEs) are considered the major moisture sources for precipitation over Mainland China during the boreal summer season. In this study, a Lagrangian particle dispersion and transport model [FLEXible PARTicle dispersion model (FLEXPART)] driven with European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis (ERA)-Interim data was used to conduct 35-year modeling of the summer season (May–August) for 1980–2014. Based on the 6-h output over 35 years, a relatively sophisticated approach was adopted that considers the change in specific humidity with trajectory tracking to diagnose OME-based precipitation during the summer season in China. We specifically explored the spatiotemporal structure of OME-based precipitation over Mainland China with a focus on quantifying the relative contributions of three specific oceanic sub-regions: the Arabian Sea (AS), the Bay of Bengal (BOB), and the South China Sea (SCS). The relevance of the OME anomalies from the three sub-regions and the observed precipitation changes on an interannual scale were also explored. The main research conclusions are summarized as follows: (1) The diagnosed OME-based precipitation and gauge observations exhibit similar spatial patterns in both seasonal and sub-seasonal scales, further evidencing the robustness of the approach used in this study. (2) Climatologically, the OMEs originating from the AS, the BOB, and the SCS made roughly equivalent contributions to the entire areal-averaged precipitation over Mainland China on a seasonal scale, but the preferred regions influenced by the three oceanic sources differ strongly from each other. (3) The relative contributions of OME from three specific subsections to precipitation varied significantly on the sub-seasonal scale. During the onset of summer monsoons, the AS region ranked first as an important oceanic source, followed by the BOB and the SCS, whereas during the withdrawal of summer monsoons, this order was reversed. (4) The interannual anomalies of OME-based precipitation from the SCS and the BOB regions are negatively correlated with those outside the AS region.

Numerical simulations of flow pattern and particle trajectories in feline aorta for hypertrophic cardiomyopathy heart conditions

ABSTRACTNumerical simulations of pulsatile flow in a feline aorta for hypertrophic cardiomyopathy (HCM) heart conditions are performed to predict flow details and to evaluate possible thrombus trajectory patterns. The study demonstrates that average flow rate boundary conditions (QBC) at artery outlets act as a resistance-type boundary condition for pulsatile flow. For simulations when the exact artery outflows are not known, specification of estimated values from physiological conditions is a plausible approach. This boundary condition is further improved using an iterative method (I-QBC) to accurately satisfy outflow conditions when expected outflow is known. The approach is validated against experimental data for the prediction of iliac artery flow and wall stresses in a human abdominal aorta. The feline aorta simulations including Lagrangian particle transport are performed on grids with up to 11M cells for a generalized feline aorta. It is found that LES on larger grids performs significantly better than URANS for the prediction of vortical structures. Simulations for both healthy and HCM conditions show similar flow patterns in the upper abdominal aorta. However, the HCM condition shows the presence of large recirculation regions in the thoracic aorta resulting in 50% lower flow through the iliac arteries and increased entrapment of fluid-borne particles near the trifurcation region compared to the healthy condition.

Modeled Population Connectivity across the Hawaiian Archipelago.

We present the first comprehensive estimate of connectivity of passive pelagic particles released from coral reef habitat throughout the Hawaiian Archipelago. Potential connectivity is calculated using a Lagrangian particle transport model coupled offline with currents generated by an oceanographic circulation model, MITgcm. The connectivity matrices show a surprising degree of self-recruitment and directional dispersal towards the northwest, from the Main Hawaiian Islands (MHI) to the northwestern Hawaiian Islands (NWHI). We identify three predicted connectivity breaks in the archipelago, that is, areas in the mid and northern part of the archipelago that have limited connections with surrounding islands and reefs. Predicted regions of limited connectivity generally match observed patterns of genetic structure reported for coral reef species in the uninhabited NWHI, but multiple genetic breaks observed in the inhabited MHI are not explained by passive dispersal. The better congruence in our modeling results based on physical transport of passive particles in the low-lying atolls of the uninhabited NWHI, but not in the anthropogenically impacted high islands of the MHI begs the question: what ultimately controls connectivity in this system?

Time history of diesel particle deposition in cylindrical dielectric barrier discharge reactors

Non-thermal plasma (NTP) treatment reactors have recently been developed for elimination of diesel particulate matter for reducing both the mass and number concentration of particles. The role of the plasma itself is obscured by the phenomenon of particle deposition on the reactor surface. Therefore, in this study, the Lagrangian particle transport model is used to simulate the dispersion and deposition of nano-particles in the range of 5 to 500 nm in a NTP reactor in the absence of an electric field. A conventional cylindrical dielectric barrier discharge reactor is selected for the analysis. Brownian diffusion, gravity and Saffman lift forces were included in the simulations, and the deposition efficiencies of different sized diesel particles were studied. The results show that for the studied particle diameters, the effect of Saffman lift is negligible and gravity only affects the motion of particles with a diameter of 500 nm or larger. Time histories of particle transport and deposition were evaluated for one-time injection and a continuous (multiple-time) injection. The results show that the number of deposited particles for one-time injection is identical to the number of deposited particles for multiple-time injections when adjusted with the shift in time. Furthermore, the maximum number of escaped particles occurs at 0.045 s after the injection for all particle diameters. The presented results show that some particle reduction previously ascribed to plasma treatment has ignored contributions from the surface deposition.

Modeling watershed-scale solute transport using an integrated, process-based hydrologic model with applications to bacterial fate and transport

Summary Distributed hydrologic models that simulate fate and transport processes at sub-daily timescales are useful tools for estimating pollutant loads exported from watersheds to lakes and oceans downstream. There has been considerable interest in the application of integrated process-based hydrologic models in recent years. While the models have been applied to address questions of water quantity and to better understand linkages between hydrology and land surface processes, routine applications of these models to address water quality issues are currently limited. In this paper, we first describe a general process-based watershed-scale solute transport modeling framework, based on an operator splitting strategy and a Lagrangian particle transport method combined with dispersion and reactions. The transport and the hydrologic modules are tightly coupled and the interactions among different hydrologic components are explicitly modeled. We test transport modules using data from plot-scale experiments and available analytical solutions for different hydrologic domains. The numerical solutions are also compared with an analytical solution for groundwater transit times with interactions between surface and subsurface flows. Finally, we demonstrate the application of the model to simulate bacterial fate and transport in the Red Cedar River watershed in Michigan and test hypotheses about sources and transport pathways. The watershed bacterial fate and transport model is expected to be useful for making near real-time predictions at marine and freshwater beaches.

A Lagrangian Particle Transport and Diffusion Model for Non-passive Scalars

We present a model of the transport and diffusion of materials whose density differs from the ambient air density or where there is excess momentum above the ambient flow. The model describes the evolution of the material concentration and can be applied to both positively and negatively buoyant pollutants. In order to calculate the transport and diffusion, we use a Monte-Carlo technique to follow the Lagrangian trajectories of many particles. This makes the model suitable for the calculation of transport and diffusion in non-homogeneous areas such as complex terrain or within an urban canopy. Results are compared with wind-tunnel observations of neutral jet releases under different conditions. Comparisons with wind-tunnel observations are made also for releases of negatively and positively buoyant materials. The agreement with observations is good.

Aerosolized drug delivery in patient-specific lung model during invasive high frequency oscillatory ventilation

A combination of large eddy simulation along with the Lagrangian particle transport model was used to investigate a gas–solid two-phase flow in realistic human lung model under high frequency oscillatory ventilation (HFOV) conditions. The resistance and compliance of the airways were used to implement a physiological time-dependent pressure boundary condition at lung peripheral outlets and results are compared with simulations using traditional boundary conditions that prescribe an outlet mass flow rate fraction. Global and local particle deposition were quantified under different HFOV and lung conditions. The coupled resistance–compliance boundary condition demonstrated a significant effect on the local deposition and pattern inside the tracheobronchial tree, due to the change in the airflow distribution and existence of the pendelluft flow. The inspiratory flow produced five-fold greater deposition than expiratory flow, and the main bronchi and lower lung lobes received the majority of the drug aerosol. The lung compliance variation was observed to have a significant effect, particularly on the local deposition during expiration phase. Increasing the operating frequency resulted in a significant change in the global and local deposition indicating strong dependency on the frequency, which could be beneficial for the targeted drug delivery. The global deposition as a fraction of the total injected particles at the endotracheal tube inlet was equivalent to the cases of normal breathing and conventional mechanical ventilation signifying a potential for efficient drug delivery during HFOV. In addition, HFOV had a unique characterization of the local particle deposition due to the rapid ventilation process and a strong influence of the endotracheal tube jet.