Lagrangian Particle Research Articles

Deciphering the Evolution of Inertial Migration in Serpentine Channels.

Serpentine channels coupling inertial and secondary flows enable effective particle focusing and separation, showing great potential in clinical diagnostics and drug screening. However, the nonsteady secondary flows in the serpentine channel make the evolution of inertial migration unclear, hindering the development and application of the serpentine channel. Herein, to refine the inertial migration mechanism, we established a model with varying curvature ratios to study the effect of secondary flow on particle migration in the serpentine channel. This method used direct numerical simulation (DNS) to calculate inertial lift, mapped the inertial lift to cross sections of the serpentine channel, and deciphered the inertial migration by using the Lagrangian particle tracking (LPT) method. The inertial migration of microparticles is experimentally investigated to validate the established numerical model. The results indicate that particle migration in serpentine channels follows a two-stage migration. An increase in secondary flow accelerates the second stage of the migration process while slowing the first stage process. Subsequently, we investigated the effects of different parameters, including Reynolds number, aspect ratio, and blockage ratio, on the equilibrium positions of particles, providing guidelines for the high-resolution separation of particles. Taking flow resistance into account, the dimensionless study makes the separation of arbitrary-sized particles possible. This work reveals the migration mechanism in serpentine channels, paving the way for the inertial separation of the particles.

A new ship tracing technology from oil spills based on multi-source data

Oil spill from ship can cause serious pollution to the Marine environment, but it is very difficult to find and confirm the troublemaker. In order to determine the oil spill ship, this paper proposes a new method to trace the source of ship oil spills and find the suspected ship that spills oil based on SAR imagery, AIS data and related marine environment data. First, we filter AIS data based on position of oil spill areas on remote sensing imagery and convert oil spill areas into trajectory points. Secondly, based on the Lagrangian particle motion model, a bidirectional drift model is proposed to calculate the average similarity between the forward and backward drift results. Finally, the most likely oil spill ship is determined according to the average similarity results. The results of the case study show that the method is effective and practical.

Spatial source, simulating improvement, and short-term health effect of high PM2.5 exposure during mutation event in the key urban agglomeration regions in China

Air quality in China has significantly improved owing to the effective implementation of pollution control measures. However, mutation events caused by short-term spikes in PM2.5 in urban agglomeration regions continue to occur frequently. Identifying the spatial sources and influencing factors, as well as improving the prediction accuracy of high PM2.5 during mutation events, are crucial for public health. In this study, we firstly introduced discrete wavelet transform (DWT) to identify the mutation events with high PM2.5 concentration in the four key urban agglomerations, and evaluated the spatial sources for the polluted scenario using Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. Additionally, DWT was combined with a widely used artificial neural network (ANN) to improve the prediction accuracy of PM2.5 concentration seven days in advance (seven-day forecast). Results indicated that mutation events commonly occurred in the northern regions during winter time, which were under the control of both short-range transportation of dirty airmass as well as negative meteorology conditions. Compared with the ANN model alone, the average band errors decreased by 9% when using DWT-ANN model. The average correlation coefficient (R) and root mean square error (RMSE) obtained using the DWT-ANN improved by 10% and 12% compared to those obtained using the ANN, indicating the efficiency and accuracy of simulating PM2.5, by combining the DWT and ANN. The short-term mortality during mutation events was then calculated, with the total averted all-cause, cardiovascular, and respiratory deaths in the four regions, being 4751, 2554, and 582 persons, respectively. A declining trend in prevented deaths from 2018 to 2020 demonstrated that the pollution intensity during mutation events gradually decreased owing to the implementation of the Three-Year Action Plan to Win the Blue Sky Defense War. The method proposed in this study can be used by policymakers to take preventive measures in response to a sudden increase in PM2.5, thereby ensuring public health.

The screening evaluation of environmental odors: a new dispersion modelling-based tool.

Odor pollution is the biggest source of complaints from citizens concerning environmental issues after noise. Often, the need for corrective actions is evaluated through simulations performed with atmospheric dispersion models. To save resources, air pollution control institutions perform a first-level odor impact assessment, for screening purposes. This is often based on Gaussian dispersion models (GDM), which does not need high computational power. However, their outputs tend to be conservative regarding the analyzed situation, rather than representative of the real in-site conditions. Hence, regulations and guidelines adopted at an institutional level for authorization/control purposes are based on Lagrangian particle dispersion models (LPDM). These models grant a more accurate simulation of the pollutants' dispersion even if they are more demanding regarding both technical skills and computing power. The present study aims to increase the accuracy of screening odor impact assessment by identifying the correlation function of the outputs derived from the two simulation models. The case study is placed in northern Italy, where a single-point source, with various stack heights, was considered. The case study is placed in northern Italy, where a single-point source, with various stack heights, was considered. The obtained correlation functions allow the practitioner to have a more accurate first-level odor impact assessment, to save time for training, and to reduce the site-specific meteorological data before proceeding with the simulation. The identified functions could allow institutions to estimate the results that would have been forecasted with the application of the more complex LPDM, applying, however, the much simpler GDM. This solution grants an accurate tool which can be used to address citizens' concerns while saving workforce and technical resources. Limitations are related to the specificity of the method regarding type sources, orography, and meteorological conditions. Comparison with other screening tools is also presented and discussed.

Source Levels of In‐Cloud Air in Shallow Cumulus: Consistency Between Paluch Diagram and Lagrangian Particle Tracking

AbstractThe Paluch diagram is a widely used tool for interpreting aircraft measurements of shallow cumulus clouds. A prior study conducted by Heus et al. (2008, https://doi.org/10.1175/2008jas2572.1) concluded that the source levels of in‐cloud air inferred from the Paluch diagram exhibit biases, sometimes of several hundred meters, in comparison to those derived from Lagrangian particle tracking. In this short study we revisit this comparison. The results indicate that the upper source levels of in‐cloud air determined from the Lagrangian Particle Tracking and the Paluch diagram are consistent, and the choice of statistical methods is crucial. The significance of this research lies in confirming the reliability of the Paluch analysis, enabling its confident application to aircraft data.

Lagrangian Versus Eulerian Spectral Estimates of Surface Kinetic Energy Over the Global Ocean

AbstractIn this study, we conducted a novel massive Lagrangian simulation experiment based on a global 1/48° tide‐resolving numerical simulation of the ocean circulation. This first‐time twin experiment enables a comparison between Eulerian (fixed‐point) and Lagrangian (along‐flow) estimates of kinetic energy (KE) across the global ocean, and the quantification of systematic differences between both types of estimations. This comparison represents an important step forward for the mapping of upper ocean high‐frequency variability from Lagrangian observations of the Global Drifter Program. Eulerian KE rotary frequency spectra and band‐integrated energy levels (e.g., tidal and near‐inertial) serve as references and are compared to Lagrangian estimates. Our analysis reveals that, except for the near‐inertial band, Lagrangian velocity spectra are systematically smoother, for example, with wider and lower spectral peaks compared to Eulerian counterparts. On average, Lagrangian KE levels derived from spectral band integrations tend to underestimate Eulerian KE levels at low‐frequency and tidal bands, especially in regions with strong low‐frequency KE. Better agreement between Lagrangian and Eulerian low‐frequency and tidal KE levels is generally found in regions with weak low‐frequency KE and/or convergent surface circulation, where Lagrangian particles tend to accumulate. Conversely, Lagrangian and Eulerian near‐inertial spectra and energy levels are comparable. Our results demonstrate that Lagrangian estimates may provide a distorted view of low‐frequency and tidal variance. To accurately map near‐surface velocity climatology at these frequencies from drifter database, conversion methods accounting for the Lagrangian bias need to be developed.

Updated Lagrangian particle hydrodynamics (ULPH) simulations of underwater bubble motions in three-dimensional space

Updated Lagrangian particle hydrodynamics (ULPH) simulations of underwater bubble motions in three-dimensional space

Numerical Study of Air–Vapor–Particle Flow in Gas Cyclone

AbstractThis paper presents a numerical study of the complex multiphase flow of air, vapor, and particles in an innovative gas cyclone called CAP cyclone. The air–vapor flow is modeled as a mixture by the Reynolds‐averaged Navier‐Stokes equations with the mixture species transport and the Eulerian wall film model using FLUENT. Particle flow is modeled by the Lagrangian particle tracking method and the condensational growth of particle droplets is modeled via a user‐defined function. The model is validated by reaching good agreement with experimental results. Flow field analysis shows that the added vapor does not change the major vortex characteristics in the cyclone, but the vapor distribution is not uniform. The vapor concentration is much higher in the upper part than in the lower part of the cyclone, leading to insufficient condensational growth in the lower part. A secondary vapor injection is proposed to improve the vapor concentration in the lower part, which is shown to be effective in improving the collection efficiency. The model and the results are helpful to the understanding and optimization of the CAP cyclone technology, and also the vapor and particle droplet flow in turbulent flows.

Effect of secondary flow and wall collisions on particle-laden flows in 90° pipe bends

Fully developed, statistically stationary particle-laden turbulent flows in a 90° pipe bend at a moderate Reynolds number are studied using direct numerical simulation coupled to a Lagrangian particle tracking technique. Three populations of particles are studied which are two-way coupled with the carrier phase. The focus is the investigation of the effect of strong curvature on particle transport dynamics across a range of particle inertias. A validation of the carrier phase predictions is performed, with good agreement found with available experimental data. It is demonstrated that the presence of particles notably affects the turbulence statistics of the carrier phase, decreasing the mean fluid velocity in the outer bend while slightly increasing it in the inner bend. Two intriguing phenomena emerge: the presence of a particle void near the inner bend and reflection layers near the outer bend. Centrifugal forces are responsible for driving particles to the outer wall, with some rebounding back into the main body of the flow, and others, influenced by the secondary flow in the plane of the pipe cross-section, converging towards the inner wall before re-entering the central pipe region. These effects dictate the size and shape of the particle void region and the formation of reflection layers due to particle-wall collisions. These two phenomena have a significant influence on the particle distribution in the pipe bend, as well as on the first- and second-order moments of the velocity and acceleration of the particulate phase. Finally, a heat map of particle-wall collision statistics indicates that the effect of the secondary flow on particle distribution and particle-wall collisions is not negligible.

Disentangling marine plastic impacts in Life Cycle Assessment: Spatially explicit Characterization Factors for ecosystem quality

Inputs of persistent plastic items to marine environments continue to pose a serious and long-term threat to marine fauna and ecosystem health, justifying further interventions on local and global scales. While Life Cycle Assessment (LCA) is frequently used for sustainability evaluations by industries and policymakers, plastic leakage to the environment and its subsequent impacts remains absent from the framework. Incorporating plastic pollution in the assessments requires development of both inventories and impact assessment methods. Here, we propose spatially explicit Characterization Factors (CF) for quantifying the impacts of plastic entanglement on marine megafauna (mammals, birds and reptiles) on a global scale. We utilize Lagrangian particle tracking and a Species Sensitivity Distribution (SSD) model along with species susceptibility records to estimate potential entanglement impacts stemming from lost plastic-based fishing gear. By simulating plastic losses from fishing hotspots within all Exclusive Economic Zones (EEZs) we provide country-specific impact estimates for use in LCA. The impacts were found to be similar across regions, although the median CF associated with Oceania was higher compared to Europe, Africa and Asia. Our findings underscore the presence of susceptible species across the world and the transboundary issue of plastic pollution. We discuss the application of the factors and identify areas of further refinement that can contribute towards a comprehensive assessment of macroplastic pollution in sustainability assessments. Degradation and beaching rates for different types of fishing gear remain a research gap, along with population-level effects on marine taxa beyond surface breathing megafauna. Increasing the coverage of impacts specific to the marine realm in LCA alongside other stressors can facilitate informed decision-making towards more sustainable marine resource management.

Embedment of WENO-Z reconstruction in Lagrangian WLS scheme implemented on GPU for strongly-compressible multi-phase flows

It is a notoriously challenge problem for Lagrangian particle methods to solve compressible flows with strong discontinuities and large volume variations. To overcome the problems of numerical oscillations, an accurate and stable Lagrangian Weighted-Least-Square (LWLS) method embedded with WENO-Z scheme (LWLS-WENO-Z) is proposed. To further enhance the numerical accuracy and stability, the particle shifting technique, the multi-level local-refinement technique, and the interface repulsive force are appropriately employed in the numerical scheme. Moreover, the smooth and stable coupled particle method between the LWLS and Roe's Riemann solver (LWLS-RR) is considered for the comparisons, further demonstrating the merit of the proposed LWLS-WENO-Z scheme. To reduce the simulation time, parallel computing with a CUDA-based program for the proposed scheme is implemented on a single GPU. In the numerical experiments, several 1D/2D benchmarks are solved to test the accuracy of the proposed method. Particularly, the challenge 2D underwater explosion and the cavitation bubble collapse and jetting are numerically investigated. Compared to other reference solutions, the present coupled particle method demonstrates robustness and superior accuracy in simulating these strongly-compressible multi-phase flows.

Stokes drift in crossing windsea and swell, and its effect on near-shore particle transport in Lofoten, Northern Norway

Transport assessments near the coast are often related to particles drifting near the surface. Such “particles” may be salmon lice, cod eggs, macro plastics or ship debris. Their drift depends on the Eulerian currents and the Stokes drift associated with the wind-generated surface wave field. The Stokes drift must be parameterized, and in doing so, one will inevitably make bulk estimates of the direction and speed of the drift for the wave spectrum. This paper implements a recently proposed parameterization of the Stokes drift vertical profile, which accounts for the effect of swell and windsea that propagate in different directions, to the open-source Lagrangian particle tracking model OpenDrift. We investigate how particles drifting in the Lofoten archipelago in Northern Norway depend on the Stokes drift and how we parameterize it. The parameterization accounting for crossing windsea and swell leads to lower residence time near the coastline than other popular parameterizations in the domain we studied.

Eulerian-Lagrangian Fluid Simulation on Particle Flow Maps

We propose a novel Particle Flow Map (PFM) method to enable accurate long-range advection for incompressible fluid simulation. The foundation of our method is the observation that a particle trajectory generated in a forward simulation naturally embodies a perfect flow map. Centered on this concept, we have developed an Eulerian-Lagrangian framework comprising four essential components: Lagrangian particles for a natural and precise representation of bidirectional flow maps; a dual-scale map representation to accommodate the mapping of various flow quantities; a particle-to-grid interpolation scheme for accurate quantity transfer from particles to grid nodes; and a hybrid impulse-based solver to enforce incompressibility on the grid. The efficacy of PFM has been demonstrated through various simulation scenarios, highlighting the evolution of complex vortical structures and the details of turbulent flows. Notably, compared to NFM, PFM reduces computing time by up to 49 times and memory consumption by up to 41%, while enhancing vorticity preservation as evidenced in various tests like leapfrog, vortex tube, and turbulent flow.

The Effect of Model Input Uncertainty on the Simulation of Typical Pollutant Transport in the Coastal Waters of China

Route planning to evade potential pollution holds critical importance for aquaculture vessels. This study establishes a fish-feed pollutant drift model based on the Lagrangian particle tracking algorithm and designs four sets of sensitivity experiments in the East China Sea. The research investigates the impact of model input uncertainties on the drift trajectory, centroid position, and sweeping area of the fish-feed pollutants. Numerical results indicate that the uncertainty in the background flow field significantly affects the uncertainty in the centroid position and sweeping area in the numerical simulations. Specifically, when a 35% random error is added to the background flow field, the centroid shift distance reaches its maximum, and the sweeping area also attains its largest value. The uncertainty in the background wind field affects the centroid position of particles but to a much lesser extent compared to the background flow field. When considering only the uncertainty of the background wind field, the sweeping area does not significantly differ from the control experiment as the uncertainty of the background wind field increases. The initial release position has little effect on the drift direction of the fish-feed pollutants but does affect the drift distance; it has minimal impact on the trajectory but significantly affects the final position of the pollutant centroid. By analyzing the model uncertainties, this study reveals the key factors influencing the drift of fish-feed pollutants. This information is crucial for aquaculture vessels in planning routes, considering environmental factors, and reducing potential pollution risks.

CFD Modeling of Aerosol Transport and Deposition Using a Drift-Flux Model

Aerosol transport and deposition are important processes in modeling of accident scenarios for a small modular reactor. An aerosol drift-flux model is attractive because it is computationally less expensive than Lagrangian particle tracking. It must be determined, however, how well it performs when implemented in a commercial computational fluid dynamics (CFD) code. This work presents results of modeling aerosol transport and deposition using a full Eulerian three-dimensional drift-flux model implemented in the commercial CFD code STAR-CCM+. The forces due to gravity and thermophoresis are included in the present drift-flux model along with Brownian motion and turbulent diffusion. The forces are added as a source term to a passive scalar transport equation. In addition, a drift velocity representing the forces is used in a built-in electrochemical species transport equation. The results of these two approaches are compared. An appropriate deposition velocity is used to calculate the aerosol concentration deposited on surfaces. The semiempirical relation proposed by Lai and Nazaroff (2000) is used to compute the deposition velocity due to gravitational settling, and the present results are compared with the experimental and numerical data obtained from the work of Chen et al. (2006). It was found that the concentration profile obtained from the present drift-flux model showed reasonable agreement with the literature data. A thermophoresis model showed good agreement when compared with the analytical solution of Nazaroff and Cass (1987). In addition to the particle concentration results, this work presents details of the drift-flux model implementation and the bulk flows. These extra details will enable comparisons by others developing similar models.

Transport Pathways for Iron Supply to the Australian Antarctic Ridge Phytoplankton Bloom

AbstractBiological productivity in the Southern Ocean is modulated by iron availability. Every summer, a large phytoplankton bloom forms northwest of the Ross Sea, above the Antarctic Australian Ridge (AAR), due to a plume of iron‐rich waters. Here, we investigate the origin and trajectories of these iron‐rich waters by analyzing water mass observations and Lagrangian experiments. Output from the Southern Ocean State Estimate (SOSE) and in situ measurements reveal that iron‐rich AAR bloom waters share properties with Modified Circumpolar Deep Water (MCDW), which forms on the Antarctic shelf‐slope. The Lagrangian experiments are conducted using SOSE velocities. Bloom waters tracked with virtual Lagrangian particles highlight an along isopycnal pathway of MCDW from Antarctica's shelf‐slope to the AAR bloom site, illustrating advection of these waters by the Balleny Gyre. These results are supported by temperature‐salinity analyses, which show a correlation between waters advected northwards; MCDW properties; and high iron concentrations.

The relation between aortic morphology and transcatheter aortic heart valve thrombosis: Particle tracing and platelet activation in larger aortic roots with and without neo-sinus

Transcatheter aortic heart valve thrombosis (THVT) affects long-term valve durability, transvalvular pressure gradient and leaflet mobility. In this study, we conduct high-fidelity fluid–structure interaction simulations to perform Lagrangian particle tracing in a generic model with larger aortic diameters (THVT model) with and without neo-sinus which is compared to a model of unaffected TAVI patients (control model). Platelet activation indices are computed for each particle to assess the risk of thrombus formation induced by high shear stresses followed by flow stagnation. Particle tracing indicates that fewer particles contribute to sinus washout of the THVT model with and without neo-sinus compared to the control model (−34.9%/−34.1%). Stagnating particles in the native sinus of the THVT model show higher platelet activation indices than for the control model (+39.6% without neo-sinus, +45.3% with neo-sinus). Highest activation indices are present for particles stagnating in the neo-sinus of the larger aorta representing THVT patients (+80.2% compared to control).This fluid–structure interaction (FSI) study suggests that larger aortas lead to less efficient sinus washout in combination with higher risk of platelet activation among stagnating particles, especially within the neo-sinus. This could explain (a) a higher occurrence of thrombus formation in transcatheter valves compared to surgical valves without neo-sinus and (b) the neo-sinus as the prevalent region for thrombi in TAV. Pre-procedural identification of larger aortic roots could contribute to better risk assessment of patients and improved selection of a patient-specific anti-coagulation therapy.

Near-Wall Flow Features In ZPG-TBL At Various Reynolds Numbers Using Dense 3D Lagrangian Particle Tracking

The 3D Lagrangian particle tracking method Shake-The-Box (STB) and, subsequently, the data assimilation method FlowFit3 have been applied to measure the flow field in a narrow wall-parallel volume inside the nearwall region of a zero-pressure-gradient (ZPG) turbulent boundary layer (TBL) flow at various Reynolds numbers. The STB experiments have been conducted in the 1m- wind tunnel facility of DLR in Göttingen at relatively high particle image densities (~0.07 ppp) using µm-sized DEHS particles and high-repetition pulse laser at up to 38.1 kHz. The results show that the chosen methodology is able to volumetrically resolve almost all features of the flow inside the viscous sub-, buffer- and lower logarithmic- layer in space and time at Reynolds numbers up to Re θ = 5,455. The gained data consists of time-series of spatially well resolved 3D velocity, acceleration and pressure fields, as well as a huge amount of long Lagrangian particle trajectories. Selecting the viscous sublayer part of the measurement volume provides time-series of the 2D distribution of both components of the instantaneous wall-shear stresses w,x/z . The resulting data are used for conditional and time-resolved 3D analyses of coherent structures, rare (e.g. back flow) events and two-point space-timecorrelations.

Scanning Lagrangian Particle Tracking To Measure 3D Large Scale Aerodynamics Of Quadcopter Flight

We present large-scale time-resolved and volumetric measurements of the flow surrounding a quadcopter drone. Cases with hovering flight, forward flight - with and without ground-effect - as well as lift-off scenarios have been realized. Five high-speed cameras were used to capture images of hundreds of thousands heliumfilled soap bubbles being illuminated by pulsed LED arrays. The arrays were operated in a dual-scanning mode, where the front and the back part of the volume are illuminated intermittently, each with a frequency of 3 kHz. The cameras were recording at 6 kHz, thereby capturing both sub-time-series in a single measurement run. Compared to a full illumination of the volume, we can either reduce the number of particles imaged on the cameras given a fixed tracer concentration - facilitating the reconstruction - or operate at a higher tracer concentration while keeping the number of imaged particles fixed. The dual time-series were evaluated by Lagrangian Particle Tracking, using the Shake-The-Box method, yielding dense fields of particle tracks. For the current results, the time series for both sub-volumes were evaluated independently. For the final contribution, an integrated evaluation scheme, alternatingly operating on both time-series is foreseen. After the tracking step, flow structures were identified using the data assimilation method FlowFit. The temporally and spatially highly resolved results document the strong interaction of the wakes induced by several rotors, as well as the interaction of wing tip vortices stemming from a single rotor. Brown-out effects are documented in forward flight. Averaged results of hovering flight enable an evaluation of the azimuthal distribution of the downwash- and outwash-velocities resulting from the geometric layout of the four rotors. By having highly resolved data encompassing the whole flight vehicle, the forces conveyed by the drone to the liquid can be determined via surface- or volume-integration.

Validating Volumetric Measurements Using Helium-Filled Soap Bubbles in a Turbulent Boundary Layer And Wake Of An Axisymmetric Body Of Revolution

This study focuses on the validation of a volumetric measurement technique to characterize the turbulent boundary layers (TBL) and wake of a slender axisymmetric body. Specifically, Lagrangian Particle Tracking (LPT) was employed to analyse a TBL subjected to an adverse pressure gradient (APG) and transverse curvature at the tail of an axisymmetric body, with a Reynolds number based on model length of 3.95×10 6 . For the LPT method, a seeding rake with 100 nozzles was used to disperse helium-filled soap bubbles (HFSB), which served as tracers. These tracers were illuminated by high-powered LED arrays, allowing their motion to be captured by three high-speed cameras at acquisition rates up to 20 kHz. A key aspect of this research was using shadowgraphy to measure the diameter of the bubbles exiting the nozzles and obtaining their statistical variations across the various nozzles in the seeding rake. The study identified that the diameter of bubbles averaged 0.45 mm, indicating a state of near-neutral buoyancy. The LPT results were then obtained using the shake-the-box (STB) tracking alogirthm, while ensemble bin-averaging of the scattered particle tracks was used to obtain turbulence statistics on a regular grid. To validate the LPT results, comparisons were made with time-resolved 2D Particle Image Velocimetry (TR-PIV) and highly resolved hot-wire anemometry. This comparative analysis revealed good agreement of the mean profiles across the methodologies, highlighting the robustness and accuracy of the LPT technique. While the turbulence levels in the outer layer were well captured by the LPT, some discrepancies were observed near the wall in both PIV and LPT data, indicating limitations in capturing small-scale phenomena. The implications of these findings, particularly in understanding the complex 3D flow effects of such TBLs and the resulting wake flow field are discussed.