Passive Tracer Research Articles

Numerical Analysis of Convective Mass Transfer during Multi-Droplet Impingement on a Structured Surface in the Presence of an Adhered Liquid Film—An Application to Spray Etching of PCBs

Multi-droplet impingement is a fundamental aspect inherent to all kinds of technical spray processes which typically aim at enhancing the convective exchange of reagents or heat at the impinged surface. In this paper, the impingement of multiple droplets onto a structured surface is investigated by a comprehensive CFD model, which resolves the dynamics of the individual droplets and the film on a micro-scale level based on the Volume of Fluid (VOF) method. The considered surface topology includes cavities and is typical for protective masks used in the spray etching of Printed Circuit Boards (PCBs). The agitation of the liquid film in terms of the convective mass transfer rates across virtual horizontal evaluation planes is studied and the influence of film height and droplet impaction velocity is elaborated. Passive tracer tracking is employed to investigate the release and re-entrainment of fluid at the surface cavities. Two modes of mass exchange between the cavities and the main flow upon droplet impingement are identified, which are central inflow accompanied by lateral outflow (1) and lateral inflow with outflow at the opposing side (2). A statistical analysis of the allocation of tracer particles shows that high impaction velocities and low film heights correlate with an enhanced decay of tracer particles within the cavities. The susceptibility to re-entrainment is also reduced by high impaction velocities, whereas increased film heights are found to promote re-entrainment.

FEOTS v0.0.0: a new offline code for the fast equilibration of tracers in the ocean

Abstract. In this paper we introduce a new software framework for the offline calculation of tracer transport in the ocean. The Fast Equilibration of Ocean Tracers Software (FEOTS) is an end-to-end set of tools to efficiently calculate tracer distributions on a global or regional sub-domain using transport operators diagnosed from a comprehensive ocean model. To the best of our knowledge, this is the first application of a transport matrix model to an eddying ocean state. While a Newton–Krylov-based equilibration capability is still under development and not presented here, we demonstrate in this paper the transient modeling capabilities of FEOTS in an application focused on the Argentine Basin, where intense eddy activity and the Zapiola Anticyclone lead to strong mixing of water masses. The demonstration illustrates progress in developing offline passive tracer simulation capabilities, while highlighting the challenges of the impulse response functions approach in capturing tracer transports by a non-linear advection scheme. Our future work will focus on improving the computational efficiency of the code to reduce time-to-solution, using different basis functions to better represent non-linear advection operators, applying FEOTS to a parent model with unstructured grids (Ocean Model for Prediction Across Scales, MPAS-Ocean), and fully implementing a Newton–Krylov steady-state solver.

Fish Injury from Movements across Hydraulic Structures: A Review

Fish migration is essential to maintain healthy aquatic ecosystems, but hydraulic structures across rivers have impeded natural fish migration worldwide. While efforts have been made to allow fish to pass some hydraulic structures, there is limited understanding of hydrodynamic effects that cause fish injury in different hydraulic systems, such as spillways and stilling basins as well as hydropower systems. This study reviewed available literature on this topic to identify the current knowledge of fish injury thresholds in laboratory- and field-based studies of hydraulic systems. Often, the hydraulic effects that lead to fish injury have been described with time-averaged simplified parameters including shear stress, pressure changes, acceleration, vortical motions, aeration, collision, and strike, while these hydrodynamic effects often occur simultaneously in the turbulent flows across hydraulic structures, making it difficult to link specific fish injuries to a particular hydrodynamic effect. Strong variations of injury may occur, depending on the type and the intensity of hydrodynamic effects, as well as the fish species and fish sizes. Modelling can provide information of stressors, but real-world tests are needed to accurately assess fish injury and mortality. Fish injury mechanisms at hydropower turbines are well understood, however, clear understanding at other sites is lacking. Future studies should aim to report holistic hydrodynamic thresholds with associated fish injury rates. Multidisciplinary systematic research is required, including laboratory and field studies, using passive tracer sensor packages and state-of-the art instrumentation in conjunction with live fish. This can quantify stressors with meaningful parameters, aiming to improve fish safety with more sustainable design of water infrastructure that reduces fish injury when passing across hydraulic structures.

Pulses of Cold Atlantic Water in the Arctic Ocean From an Ocean Model Simulation

AbstractThe intrusion of the warm and saline Atlantic Water (AW) has significant implications and feedback on the dynamics and thermodynamics of the Arctic Ocean. The AW enters the Arctic Ocean through two gateways: Fram Strait (FS) and the Barents Sea Opening (BSO). The AW is the dominant source of oceanic heat to the Arctic Ocean. In conjunction with the measurements in key hydrographic sections, numerical ocean modeling provides us with a useful tool to characterize and corroborate the temporal and spatial variability of the AW branches. Simulations are conducted using the regional configuration Arctic and North Hemispheric Atlantic (ANHA) of the ocean/sea‐ice model Nucleus for Modeling of the Ocean running at 1/4° and 1/12° resolution. Online passive tracers from the model configuration are used to trace the pathways of the AW inflow in the Arctic Ocean. While studying the variability of the AW through FS and the BSO, a cold mode of AW (Cold Atlantic Water [CAW]) with a temperature below 0℃ is detected along the rim of the eastern Eurasian Basin. It is formed primarily due to the intense sea surface cooling of the AW in the Barents and Kara Seas. The CAW exhibits the cascading process as it propagates downstream through the St. Anna Trough, leading to a considerable mixing with the ambient water. The CAW pulses eventually result in a significant heat content reduction in the eastern Arctic. This study provides insight into the changing AW property and its potential impact on the Arctic Basin.

Geostrophic Eddies Spread Near-Inertial Wave Energy to High Frequencies

Abstract The generation of broadband wave energy frequency spectra from narrowband wave forcing in geophysical flows remains a conundrum. In contrast to the long-standing view that nonlinear wave–wave interactions drive the spreading of wave energy in frequency space, recent work suggests that Doppler-shifting by geostrophic flows may be the primary agent. We investigate this possibility by ray tracing a large number of inertia–gravity wave packets through three-dimensional, geostrophically turbulent flows generated either by a quasigeostrophic (QG) simulation or by synthetic random processes. We find that, in all cases investigated, a broadband quasi-stationary inertia–gravity wave frequency spectrum forms, irrespective of the initial frequencies and wave vectors of the packets. The frequency spectrum is well represented by a power law. A possible theoretical explanation relies on the analogy between the kinematic stretching of passive tracer gradients and the refraction of wave vectors. Consistent with this hypothesis, the spectrum of eigenvalues of the background flow velocity gradients predicts a frequency spectrum that is nearly identical to that found by integration of the ray tracing equations.

Odd Mobility of a Passive Tracer in a Chiral Active Fluid.

Chiral active fluids break both time-reversal and parity symmetry, leading to exotic transport phenomena unobservable in ordinary passive fluids. We develop a generalized Green-Kubo relation for the anomalous lift experienced by a passive tracer suspended in a two-dimensional chiral active fluid subjected to an applied force. This anomalous lift is characterized by a transport coefficient termed the odd mobility. We validate our generalized response theory using molecular dynamics simulations, and we show that the asymmetric tracer mobility may be understood mechanically in terms of asymmetric deformations of the tracer-fluid density distribution function. We show that the even and odd components of the mobility decay at different rates with tracer size, suggesting the possibility of size-based particle separation using a chiral active working fluid.

Three-Dimensionally Complex Phase Behavior and Collective Phenomena in Mixtures of Acoustically Powered Chiral Microspinners.

The process of dynamic self-organization of small building blocks is fundamental to the emergent function of living systems and is characteristic of their out-of-equilibrium homeostasis. The ability to control the interactions of synthetic particles in large groups could lead to the realization of analogous macroscopic robotic systems with microscopic complexity. Rotationally induced self-organization has been observed in biological systems and modeled theoretically, but studies of fast, autonomously moving synthetic rotors remain rare. Here, we report switchable, out-of-equilibrium hydrodynamic assembly and phase separation in suspensions of acoustically powered chiral microspinners. Semiquantitative modeling suggests that three-dimensionally (3D) complex spinners interact through viscous and weakly inertial (streaming) flows. The interactions between spinners were studied over a range of densities to construct a phase diagram, which included gaseous dimer pairing at low density, collective rotation and multiphase separation at intermediate densities, and ultimately jamming at high density. The 3D chirality of the spinners leads to self-organization in parallel planes, forming a three-dimensionally hierarchical system that goes beyond the 2D systems that have so far been modeled computationally. Dense mixtures of spinners and passive tracer particles also show active-passive phase separation. These observations are consistent with recent theoretical predictions of the hydrodynamic coupling between rotlets generated by autonomous spinners and provide an exciting experimental window to the study of colloidal active matter and microrobotic systems.

Entrainment in a Simulated Supercell Thunderstorm. Part III: The Influence of Decreased Environmental Humidity and General Effects upon Precipitation Efficiency

Abstract Entrainment is a key process that can modulate the intensity of supercells, and a better understanding of its impact could help improve forecasts of thunderstorms and the precipitation they produce. In Part III of this series, the three distinct mechanisms of entrainment identified during the mature stage of idealized supercell thunderstorms in Part I (overturning “ribbons” of horizontal vorticity, “disorganized turbulent eddies,” and the “storm-relative airstream”) are examined as the absolute humidity of the environment is decreased. The existence of these mechanisms in a more realistic simulated storm environment is also established. Entrainment is calculated as fluxes of air across the storm core surface; passive fluid tracers help determine the resulting dilution of the storm updraft. Model microphysical rates are used to examine the direct impacts of entrainment on hydrometeors within the storm updraft as well as precipitation that falls to the ground. Results show that as mixed-layer humidity decreases, the “ribbons” and turbulent eddy mechanisms decrease in intensity, but their effects on precipitation production change little. With decreasing humidity in the 3–4 km AGL layer, the storm-relative airstream entrains less humid low-level air into the storm core, decreasing the vertical mass flux, and therefore the precipitation produced by the storm. When the humidity in the mid- to upper troposphere (4–20 km AGL) is decreased, precipitation is significantly reduced, but not due to the effects of the entrained air. Rather, enhanced evaporation and sublimation of falling precipitation decreases the overall precipitation efficiency of the storm.

Lagrangian description of the atmospheric flow from Eulerian tracer advection with relaxation

AbstractThe Lagrangian representation of fluid flows offers a natural perspective to study many kinds of physical mechanisms. By contrast, the Eulerian representation is more convenient from a diagnostic point of view. This article attempts to combine elements of both worlds by proposing an Eulerian method that allows one to extract Lagrangian information about the atmospheric flow. The method is based on the offline advection of passive tracer fields and includes a relaxation term. The latter device allows one to run the integration in a continuous fashion without the need for reinitialization. As a result one obtains accumulated Lagrangian information, for example, about the recent parcel displacement or the recent parcel‐based diabatic heating, at each point of an Eulerian grid at any time step. The method is implemented with a pseudospectral algorithm suitable for gridded global atmospheric data and compared with the more traditional trajectory method. The method's utility is demonstrated on the basis of a few examples, which relate to cloud formation and the development of temperature anomalies. The examples highlight that the method provides a convenient diagnostic of parcel‐based changes, paving an intuitive way to explore the physical processes involved. Due to its gridpoint‐based nature, the proposed method can be applied to large data sets in a straightforward and computationally efficient manner, suggesting that the method is particularly useful for climatological analyses.

Enhanced Aircraft Wake Decay Under Crosswind Conditions

The starting vortex generated by air blowing from the ground destroys the aircraft wake vortex within a short period of time. Thus, air blowing might shorten flight separation and promote airport operating efficiency. A crosswind is a common phenomenon affecting airplane landing, and it can greatly affect air blowing. In this paper, a large-eddy simulation was used to study the interference exerted by different blowing zones on a wake vortex. The best-performing blowing strategy was estimated by ascertaining the relationship between the blowing zone aspect ratio and the corresponding starting vortex strength. It was found that, without the crosswind, the blowing effect remained almost constant for an aspect ratio of between one and 16. A passive tracer released from the different blowing zones proved that the mechanisms of enhanced vortex decay, with or without crosswind, were identical. However, the results showed that crosswind caused asymmetrical and reduced artificial wake decay enhancement. When the crosswind was imposed, the development of the upwind vortex was sensitive to the aspect ratio, and the blowing zone with an aspect ratio of 16 showed the best performance.

Northward movement of the tropical dinoflagellate Ornithocercus and Triposolenia genera in Korean coastal waters is strongly associated with the inflow of the Jeju Warm Current

Rising seawater temperatures are causing geographic shifts in the composition and abundance of the marine plankton that make up the microbial food web. Over the past few years, the water temperature around the Korea Strait has been increasing rapidly. Twelve cruise surveys were conducted in the Jeju Strait between 2017 and 2019 to monitor the presence of the tropical dinoflagellate genera Ornithocercus and Triposolenia, which are brought in by the Jeju Warm Current. Additionally, passive tracer experiments were performed using the Regional Ocean Modeling System to understand how the warm currents impact the Jeju Strait. The results of these simulations were then compared against the distribution ranges and abundance of Ornithocercus and Triposolenia. The study found that there were significant seasonal variations in abundance and spatial distribution of the tropical dinoflagellates during the sampling period, particularly during September and November when they were highly abundant and widely distributed. These patterns were closely related to the intensity and path of the Jeju Warm Current, particularly from southeastern offshore (EKB02; Eastern Kuroshio Branch02) rather than from southwestern offshore (EKB01; Eastern Kuroshio Branch01) of Jeju Island. The study also found no tropical dinoflagellates at any of the stations near the south coast of Korea during the sampling period. These findings suggest that the Ornithocercus and Triposolenia tropical dinoflagellate genera may serve as useful biological indicators to monitor the advection of warm currents into the Jeju Strait, Korea.

Implications of steep hilly terrain for modeling wind-turbine wakes

Wind energy is crucial to the transition to a carbon-free world. However, new wind farms are increasingly sited on complex terrain, whose influence on turbine performance is still not well understood. Here large-eddy simulations are performed on the flow around a wind turbine sited on three different terrain types: a flat ground, a 2-D hill, and a 3-D hill. We find that hilly terrain can increase the power generated via the speed-up effect, but that it can also increase the wake deficit and deflect the wake center, potentially affecting any turbines downstream. We interpret this deflection as being caused by the topography-induced vertical velocity component, which we model with a passive tracer method. By combining this passive tracer method with engineering-focused wake models, we obtain improved predictions of the velocity deficit and wake deflection. Such a hybrid strategy is fast and accurate, facilitating the design of wind farms sited on steep hilly terrain, which often experience flow separation at the backslope.

Passive tracer transport in peristaltic pumping of non-Newtonian blood flow: A mathematical model

The large time behavior of passive contaminant in non-Newtonian peristaltic blood flow in a two-dimensional (2D) channel (capillary) has been examined in this paper. The power-law model is employed in order to highlight the non-Newtonian blood characteristic. The study was conducted using the Reynolds decomposition technique, which converts a 2D transport problem into a 1D transport model in which species concentration can be decomposed into sectional average concentration and variation from its mean value. For flow velocity, the same decomposition method is used. This allows the derivation of the dispersion coefficient and convection coefficient. Using Fick’s law, the advection–diffusion equation is modified by replacing these coefficients by their corresponding average values and analytical solutions for the mean concentration are derived. In the absence of peristalsis effects ([Formula: see text]), i.e., for the straight rigid channel, the dispersion coefficient is invariant along the channel length. With increasing modulation (peristaltic wave) parameter, [Formula: see text], there is a strong elevation in advection coefficient in the initial half of the channel with a subsequent suppression in the second half of the channel, indicating that the location in the channel strongly influences advection characteristics. Advection coefficient is significantly elevated with increment in power-law rheological index (for shear-thinning fluids, [Formula: see text]) across the channel length and exhibits an oscillatory nature due to the peristaltic waves. In the shear-thickening range ([Formula: see text]), with progressive increase in n, an increment in peristaltic modulation parameter, [Formula: see text], induces a marked reduction in the axially average relative advection coefficient. Dispersion coefficient is initially boosted along the early section of the channel with increment in modulation parameter whereas further long the channel this trend is reversed. Increasing aspect ratio and Péclet number consistently boost dispersion coefficient along the entire channel length. The study provides a solid benchmark for further generalized simulations with computational fluid dynamics.

Cautionary tales from the mesoscale eddy transport tensor

The anisotropic mesoscale eddy transport tensor is diagnosed using passive tracers advected online in both an idealized 101-member mesoscale-resolving quasi-geostrophic (QG) double-gyre ensemble, and a realistic 24-member eddying (1/12∘) ensemble of the North Atlantic. We assert that the Reynold’s decomposition along the ensemble dimension, rather than the spatial or temporal dimension, allows us to capture the intrinsic spatiotemporal variability of the mean flow and eddies. The tensor exhibits good performance in reconstructing the eddy fluxes of passive tracers, here defined as fluctuations about the ensemble thickness-weighted averaged (TWA) mean. However, the inability of the tensor to reconstruct eddy fluxes of QG potential vorticity, which encapsulates the eddy-mean flow interaction, and other active tracers raises the question: To what extent can the diagnosed tensor be applied to inform the parametrization of mesoscale dynamics?

Can We Accurately Quantify a Lateral Diffusivity from a Single Tracer Release?

Abstract Mixing along sloping isopycnals plays a key role in the transport and uptake of heat and carbon by the ocean. This mixing is quantified by a lateral diffusivity, which can be measured by tracking the lateral spreading of point release tracer patches. We present a definition for the area of a tracer patch, the time derivative of which provides the lateral diffusivity. To accurately estimate the diffusivity, an ensemble mean concentration field of many tracer release experiments is required. We use numerical experiments to quantify how accurately the “true” lateral diffusivity (obtained from the ensemble mean concentration field) can be estimated from a single tracer release experiment (one ensemble member). To simulate observational campaigns, we also estimate the diffusivity from a single tracer release that is spatially and/or temporally subsampled, quantifying how the error between the estimated diffusivity and the true diffusivity grows as this sampling resolution worsens. We perform these numerical experiments in a two-layer quasigeostrophic model of turbulent flow on a β plane, using an ensemble of 50 passive tracer release experiments, each initialized as a 2D Gaussian but with differing realizations of the turbulent flow. We find that the diffusivity estimates from the single tracer releases have a relative root-mean-square error (RMSE) of 1.43% from the true diffusivity. Subsampling a single tracer release experiment every 956 km increases the relative RMSE from the true diffusivity to 3.1%; also subsampling every 277 days raises this figure to 6.5%.

Cross-Slope Transport by a Mesoscale Anticyclone in the Northern South China Sea

Cross-slope eddies play an important role in the exchange of water, salt, heat, nutrients, chlorophyll, phytoplankton and other biogeochemical elements between basin and shelf in the South China Sea. The cross-slope transport process by a mesoscale anticyclonic eddy is studied by ROMS model and satellite data. The 1000 m isobath was considered as a proxy for the slope. The anticyclone shows different features at different places on the slope: (a) the volume transport at the northeast of the slope was off-slope, while at the southwest was on-slope; (b) both on and off-slope transports were greatly enhanced during the cross-slope process, and gradually weakened after crossing the slope. The total cross-slope water transport was 5.97 Sv, which was higher than the along-slope component with −0.58 Sv. The Eulerian results also showed that enhanced cross-slope transport was related to the distance between the eddy and slope, the eddy radius, and the eccentricity of the eddy. The offline passive tracer experiment showed that particles were floating up during and after the crossing process, mainly due to the strong Ekman pumping.

Role of Kidney Stones in Renal Pelvis Flow.

Ureteroscopy is a commonly performed medical procedure to treat stones in the kidney and ureter using a ureteroscope. Throughout the procedure, saline is irrigated through the scope to aid visibility and wash-out debris from stone fragmentation. The key challenge that this research addresses is to build a fundamental understanding of the interaction between the kidney stones/stone fragments and the flow dynamics in the renal pelvis flow. We examine the time-dependent flow dynamics inside an idealized renal pelvis in the context of a surgical procedure for kidney stone removal. Here, we examine the time-dependent evolution of these vortical flow structures in three dimensions, and incorporate the presence of rigid kidney stones. We perform direct numerical simulations, solving the transient Navier-Stokes equations in a spherical domain. Our numerical predictions for the flow dynamics in the absence of stones are validated with available experimental and numerical data, and the governing parameters and flow regimes are chosen carefully in order to satisfy several clinical constraints. The results shed light on the crucial role of flow circulation in the renal cavity and its effect on the trajectories of rigid stones. We demonstrate that stones can either be washed out of the cavity along with the fluid, or be trapped in the cavity via their interaction with vortical flow structures. Additionally, we study the effect of multiple stones in the flow field within the cavity in terms of the kinetic energy, entrapped fluid volume, and the clearance rate of a passive tracer modeled via an advection-diffusion equation. We demonstrate that the flow in the presence of stones features a higher vorticity production within the cavity compared with the stone-free cases.

Dynamics of a spherical self-propelled tracer in a polymeric medium: interplay of self-propulsion, stickiness, and crowding.

We employ computer simulations to study the dynamics of a self-propelled spherical tracer particle in a viscoelastic medium, made of a long polymer chain. Here, the interplay between viscoelasticity, stickiness, and activity (self-propulsion) brings additional complexity to the tracer dynamics. Our simulations show that on increasing the stickiness of the tracer particle to the polymer beads, the dynamics of the tracer particle slows down as it gets stuck to the polymer chain and moves along with it. But with increasing self-propulsion velocity, the dynamics gets enhanced. In the case of increasing stickiness as well as activity, the non-Gaussian parameter (NGP) exhibits non-monotonic behavior, which also shows up in the re-scaled self part of the van-Hove function. Non-Gaussianity results owing to the enhanced binding events and the sticky motion of the tracer along with the chain with increasing stickiness. On the other hand, with increasing activity, initially non-Gaussianity increases as the tracer moves through the heterogeneous polymeric environment but for higher activity, the tracer escapes resulting in a negative NGP. For higher values of stickiness, the trapping time distributions of the passive tracer particle broaden and have long tails. On the other hand, for a given stickiness with increasing self-propulsion force, the trapping time distributions become narrower and have short tails. We believe that our current simulation study will be helpful in elucidating the complex motion of activity-driven probes in viscoelastic media.

Global Morphology Distortion of the 2021 October 9 Coronal Mass Ejection from an Ellipsoid to a Concave Shape

This paper presents a study of a 2021 October 9 coronal mass ejection (CME) with multipoint imaging and in situ observations. We also simulate this CME from the Sun to Earth with a passive tracer to tag the CME’s motion. The coronagraphic images show that the CME is observed as a full halo by SOHO and as a partial halo by STEREO-A. The heliospheric images reveal that the propagation speed of the CME approaches about 1° hr−1, suggesting a slow CME. With simulated results matching these observation results, the simulation discloses that as the CME ejects from the Sun out to interplanetary space, its global morphology is distorted from an ellipsoid to a concave shape owing to interactions with the bimodal solar wind. The cross section of the CME’s flux rope structure transforms from a circular shape into a flat one. As a result of the deflection, the propagation direction of the CME is far away from the Sun–Earth line. This means that the CME flank (or the ICME leg) likely arrives at both Solar Orbiter and the L1 point. From the CME’s eruption to 1 au, its volume and mass increase by about two orders and one order of magnitude, respectively. Its kinetic energy is about 100 times larger than its magnetic energy at 1 au. These results have important implications for our understanding of CMEs’ morphology, as well as their space weather impacts.

Life cycle of cosmic-ray electrons in the intracluster medium

We simulate the evolution of relativistic eletrons injected into the medium of a small galaxy cluster by a central radio galaxy, studying how the initial jet power affects the dispersal and the emission properties of radio plasma. By coupling passive tracer particles to adaptive-mesh cosmological magnetohydrodynamic (MHD) simulations, we study how cosmic-ray electrons are dispersed as a function of the input jet power. We also investigate how the latter affects the thermal and non-thermal properties of the intracluster medium, with differences discernible up to about one Gyr after the start of the jet. We evolved the energy spectra of cosmic-ray electrons, subject to energy losses that are dominated by synchrotron and inverse Compton emission as well as energy gains via re-acceleration by shock waves and turbulence. We find that in the absence of major mergers, the amount of re-acceleration experienced by cosmic-ray electrons is not enough to produce long-lived detectable radio emissions. However, for all simulations, the role of re-acceleration processes is crucial to maintaining a significant and volume-filling reservoir of fossil electrons (γ ∼ 103) for at least one Gyr after the first injection by jets. This is important in attempting to establish plausible explanations of recent discoveries of cluster-wide emission and other radio phenomena in galaxy clusters.