Estimation of bottom microplastic flux in the Indonesian seas.

Marine debris, particularly microplastics, is an important environmental problem for Indonesia, impacting vast coastline and diverse marine ecosystems. However, little is known about the microplastics flux to the seafloor in the Indonesian Seas. This study employs HYCOM data and Lagrangian particle tracking model to analyze the distribution of microplastics flux to the seafloor from 68 rivers in Indonesia. Microplastics accumulation can be found in coastal waters near major islands, typically within 50-100 km of the coast. Accumulation regions within the Fisheries Management Area (FMA) in Indonesia are mostly located in Karimata Strait, Java Sea and South Java. This study highlights the importance of understanding sinking particles behavior for effective waste management strategies, as well as for mitigating environmental impacts in Indonesian waters.

Seasonal change in fate and transport of plastics from Red River to the coast of Vietnam.

A Lagrangian-particle tracking model, Delft3D-PART, combined with hydrodynamics models are used to investigate the fate and transport of buoyant plastics from Ba Lat river mouth in Red River Delta, northern Vietnam. It was found that during the dry season (Dec-Feb), 23 % (26.43 ton) of the plastics reached the shoreline while 76.1 % (68.3 ton) moved towards the coast further south of Red River Delta. During the wet season (Jun-Aug), 42 % (56.3 ton) were transported offshore away from the coast and 20 % (26.43 ton) distributed along the shore. The two bays adjacent to the river mouth are major hotspots with the intensity skewed towards the upwind side relative to the seasonal monsoon. This phenomenon is exacerbated by storm events which reverse the typical transport and lead to formation of hotspots at the upwind side of the plastic source. Guidance of model results for targeted cleanup operations is discussed.

Investigation of cardiopulmonary bypass parameters on embolus transport in a patient-specific aorta.

Neurological complexities resulting from surgery requiring cardiopulmonary bypass (CPB) remain a major concern, encompassing a spectrum of complications including thromboembolic stroke and various cognitive impairments. Surgical manipulation during CPB is considered the primary cause of these neurological complications. This study addresses the overall lack of knowledge concerning CPB hemodynamics within the aorta, employing a combined experimental-computational modeling approach, featuring computational fluid dynamics simulations validated with an in vitro CPB flow loop under steady conditions. Parametric studies were systematically performed, varying parameters associated with CPB techniques (pump flow rate and hemodiluted blood viscosity) and properties related to formed emboli (size and density). This represents the first comprehensive investigation into the individual and combined effects of these factors. Our findings reveal critical insights into the operating conditions of CPB, indicating a positive correlation between pump flow rate and emboli transport into the aortic branches, potentially increasing the risk of stroke. It was also found that larger emboli were more often transported into the aortic branches at higher pump flow rates, while smaller emboli preferred lower flow rates. Further, as blood is commonly diluted during CPB to decrease its viscosity, more emboli were found to enter the aortic branches with greater hemodilution. The combined effects of these parameters are captured using the non-dimensional Stokes number, which was found to positively correlate with emboli transport into the aortic branches. These findings contribute to our understanding of embolic stroke risk factors during CPB and shed light on the complex interplay between CPB parameters.

Managing deepsea oil spills through a systematic modeling approach.

Offshore oil exploration and production in deepwater are associated with environmental risks to marine ecosystems. This research introduces DWOSM (Deep Water Oil Spill Model), a three-dimensional Lagrangian model, which is developed to simulate the transport and fate of oil spills resulting from subsea blowouts. DWOSM comprises three interconnected modules: DWOSM-DSD, which predicts the oil droplet size distribution from a blowout release; DWOSM-NearField, simulating plume dynamics and tracking oil droplets within the plume region; and DWOSM-FarField, modeling the evolution of dispersed oil beyond the near-field. Compared to other oil spill models, this integrated approach improves the transition between near and far fields using a near-field particle tracking algorithm. It also employs the thermodynamic models to enable the prediction of oil properties under varying deep water pressure and temperature. To gauge the reliability and efficacy of DWOSM, a hypothetical case situated within a North American context is employed for model testing. The DWOSM and its each module are juxtaposed with other established oil spill models. The outcomes indicate that DWOSM yields comparable results to these models by providing reasonable predictions of a deepwater blowout. The model's verification through case scenario testing and comparison underscores its potential as a decision tool for assessing and managing the potential environmental impacts of offshore petroleum activities.

Systematic CFD-based evaluation of physical factors influencing the spatiotemporal distribution patterns of microplastic particles in lakes.

Spatiotemporal distribution patterns of microplastic (MP) particles in lakes hinge on both the physical conditions in the lake and particle properties. Using numerical simulations, we systematically investigated the influence of lake depth and bathymetry, wind and temperature conditions, MP particle release location and timing, as well as particle diameter (10, 20, and 50 µm). Our results indicate that maximum lake depth had the greatest effect on the residence time in the water column, as it determines the settling timescale and occurrence of hydrodynamic complexity such as density-driven flows in the lake. Increasing particle size from 10 to 20 and 50 µm also significantly reduced the residence time making particle size the factor with the second strongest effect on the residence time and, in turn, on the availability of MP particles for uptake by organisms. Changing bathymetry from a uniform to a non-uniform had a less pronounced effect on particle residence time compared to maximum depth and particle size. Release location, wind conditions, and release time had comparably little effect on particle behavior but became more important as MP particle size decreased. The release of the 10 µm MP particles in the deeper lakes with uniform bathymetry during summer with stable thermal stratification, resulted in a nearly month-long turnover phase in the fall in which both settling and rising of particles occurred simultaneously. This was caused by convective heat and water transport during this period. In these scenarios about 2.6 to 5.4 % of the released MP particles were held in or returned to the water layers near the lake surface. While acknowledging the dominant role of lake depth and MP particle size on the particle residence time, this study further emphasizes that it is ultimately a particular combination of different factors and their interactions that shape MP distribution patterns in lakes.

Large-scale flow field and aerosol particle transport investigations in a classroom using 2D-Shake-The-Box Lagrangian Particle Tracking.

Infections with COVID-19 in enclosed public spaces, where virus-laden aerosol particles can accumulate over time, have significantly contributed to the rapid spread of the virus. It is therefore of great importance to understand the transport and dispersion process of aerosol particles in such spaces, especially against the background of future pandemics. In this work, we present a Lagrangian-Particle-Tracking experiment to assess the mixed convective flow in a classroom with different ventilation strategies. For this purpose, thermal plumes were created by heated dummies, and a collimated LED light-sheet with ∼0.4 m thickness was used for illumination of helium filled soap bubbles (HFSB) acting as passive tracer particles. In this way, the Lagrangian trajectories of the particles were recorded at two approximately 4.2 m × 2.8 m large fields using the novel 2D-Shake-The-Box-Method. As a result, time-resolved trajectories of over 300,000 simultaneously tracked HFSB have been reconstructed, so that both small-scale and large-scale properties of the flow are visualized quantitatively across the entire cross-section of the room. The trajectories show that the thermal plumes create lengthwise circulating vortices, which cannot be destroyed across the entire cross-section of the room by opening or tilting a window. Furthermore, the mixing in the room through the operation of an air purifier is higher compared to opening a window, which suggests that this strategy in combination with its air filtering capability is the most effective strategy to prevent infections.

Origin and pathway of a subsurface maximum of 137Cs detected in the winter of 2012 after the Fukushima accident.

This study investigates the dispersion behavior of 137Cs and evaluates its origin (atmospheric deposition or direct ocean release) from the Fukushima Dai-ichi Nuclear Power Plant (FDNPP) accident using a Lagrangian particle tracking model. The ocean circulation fields based on the Modular Ocean Model Version 5 (MOM5) were adopted for the simulation. The MOM5 results represented the formation and migration of subtropical mode water (STMW) comparable with observations and reanalysis data. Particularly, anticyclonic eddies south of the Kuroshio extension promoted surface mixing over 300 m in the cooling season. The particle tracking simulation reproduced well the maximum subsurface activity between 142 and 146°E, where STMW is deep owing to anticyclonic eddies, compared to the activity found via measurements conducted around 149°E in the winter of 2012. It also demonstrated that the 137Cs of the tropical and subtropical regions (10-35°N, 142-146°E) in the winter of 2012 almost entirely originated from atmospheric deposition.

Modelling marine microplastics accumulation zones through freshwater discharge sources: Simulated release from the Onkaparinga River and Torrens River, South Australia.

Plastic pollution is fast becoming one of the most prominent contamination issues facing the marine environment. Microplastics are a major subset of plastic waste now present in all global oceans, with their numbers standing only to increase. This study applies a coupled hydrodynamic model and Lagrangian particle-tracking model to predict and simulate microplastic transport in South Australian waters. Virtual particles representing microplastics were released daily for 365 days from two major freshwater input sources along the coastline of Adelaide, Australia. These particles entered the Gulf St Vincent and were tracked over two model years using LTRANS software. The model identified general gulf circulation as circular, clockwise, with net southward particle transport from particle release sites. A potential accumulation zone associated with a local eddy was identified. Concentrations of particles that passed through local marine parks were also calculated in response to the potential concerns they pose in vulnerable protected areas.

Particle deposition and fluid flow characteristics in turbulent corrugated pipe flow using Eulerian-Lagrangian approach.

A numerical simulation of aerosol particle deposition in a horizontal circular pipe with a corrugated wall under turbulent flow has been carried out in this research. This paper uses the RNG k-ε turbulence model with Enhanced Wall Treatment to simulate fluid flow. Furthermore, the Lagrangian particle tracking model simulates particle deposition in the corrugated pipe. Air-particle interaction is influenced by Stokes number, surface roughness, flow velocity, particle diameter, and pipe diameter. For the parametric simulation, particle diameter varies from 1 to 30 µm, whereas the Reynolds number ranges from 5000 to 10,000. The effect of corrugation height and pipe diameter on deposition efficiency is also investigated. This study shows that corrugation height significantly increases particle deposition compared to the smooth wall pipe. As the pipe diameter decreases, keeping the corrugation ratio constant, deposition efficiency also increases. Moreover, high flow velocity enhances deposition efficiency for particle diameters lower than 5 µm.

Development of a web-based modelling framework for harmful algal blooms transport simulation using open-source technologies.

Desktop-based modelling packages presented typical limitations in interactive simulation. This study presents a web-based modelling framework that fully consolidated the simulation work-flow into a WebGIS application, providing a one-step solution for HABs transport simulation within an intuitive and interactive modelling environment. An improved Lagrangian particle-tracking scheme was proposed using fractional Brownian motion technique. The presented model was devoted to quickly forecast the transport pathways in both temporal and spatial dimensions, and evaluate the approximate trends and qualitative understanding of HABs development in data-poor situations. The web modelling platform was developed using multiple open-source JavaScript libraries. The developed WebGIS application provides user-friendly interfaces to prepare inputs, configure simulation settings, visualize, analyse, and validate simulation results within the same framework. The feasibility, capacity, and advantage of the proposed framework were tested and evaluated in a real-world application of red tide transport simulation in the Qinhuangdao coastal waters. The model results showed qualitative agreement with the red tide observed from remote sensing. Our experimental results demonstrated that the developed web-based modelling prototype would present a useful performance for study cases related to HABs transport simulation.

CFD Simulation of Particle-Laden Flow in a 3D Differentially Heated Cavity Using Coarse Large Eddy Simulation.

Particulate flow in closed space is involved in many engineering applications. In this paper, the prediction of particle removal is investigated in a thermally driven 3D cavity at turbulent Rayleigh number Ra = 109 using Coarse Large Eddy Simulation (CLES). The depletion dynamics of SiO2 aerosol with aerodynamic diameters between 1.4 and 14 µm is reported in an Euler/Lagrange framework. The main focus of this work is therefore to assess the effect of the subgrid-scale motions on the prediction of the particulate flow in a buoyancy driven 3D cavity flow when the mesh resolution is coarse and below optimal LES standards. The research is motivated by the feasibility of modeling more complex particulate flows with reduced CPU cost. The cubical cavity of 0.7 m side-length is set to have a temperature difference of 39 K between the two facing cold and hot vertical walls. As a first step, the carrier fluid flow was validated by comparing the first and second-moment statistics against both previous well-resolved LES and experimental databases [Kalilainen (J. Aero Sci. 100:73-87, 2016); Dehbi (J. Aero. Sci. 103:67-82, 2017)]. First moment Eulerian statistics show a very good match with the reference data both qualitatively and quantitatively, whereas higher moments show underprediction due to the lesser spatial resolution. In a second step, six particle swarms spanning a wide range of particle Stokes numbers were computed to predict particle depletion. In particular, predictions of 1.4 and 3.5 µm particles were compared to LES and available experimental data. Particles of low inertia i.e. dp < 3.5 µm are more affected by the SGS effects, while bigger ones i.e. dp = 3.5-14 µm exhibit much less grid-dependency. Lagrangian statistics reported in both qualitative and quantitative fashions show globally a very good agreement with reference LES and experimental databases at a fraction of the CPU power needed for optimal LES.

Storm-induced coastward expansion of Margalefidinium polykrikoides bloom in Chesapeake Bay.

An unusual coastward expansion of the toxic dinoflagellate species Margalefidinium polykrikoides was observed in 2020 summer after a tropical storm passing Chesapeake Bay. Such coastward expansion was only recorded in 2007. A newly developed coupled Lagrangian particle tracking and harmful algal bloom model driven by environmental variables was used to investigate the underlying mechanisms and successfully reproduced the expansion patterns. Persistent pre-storm southerly winds favored the delivery of bloom source water originated inside the bay to the coast. Storm-induced strong upwelling of denser subsurface water interacted with the after-storm outflow plume (steered southward as the storm's impacts waned), forming a transport barrier to accumulate algae and delineate the coastwide bloom extent. Algal diel vertical migrations and transport barrier enable algae to stay in the nearshore regions. The storm-induced coastward expansion of M. polykrikoides might increase future bloom possibility in the coastal area.

A cell-resolved, Lagrangian solver for modeling red blood cell dynamics in macroscale flows.

When red blood cells (RBCs) experience non-physiologically high stresses, e.g., in medical devices, they can rupture in a process called hemolysis. Directly simulating this process is computationally unaffordable given that the length scales of a medical device are several orders of magnitude larger than that of a RBC. To overcome this separation of scales, the present work introduces an affordable computational framework that accurately resolves the stress and deformation of a RBC in a spatially and temporally varying macroscale flow field such as those found in a typical medical device. The underlying idea of the present framework is to treat RBCs as one-way coupled tracers in the macroscale flow by capturing the effect of the flow on their dynamics but neglecting their effect on the flow at the macroscale. As a result, the RBC dynamics are simulated after those of the flow in a postprocessing step by receiving the fluid velocity gradient tensor measured along the RBC trajectory as the input. To resolve the fluid velocity in the immediate vicinity of the RBC as well as the motion of the membrane, we employ the boundary integral method coupled to a structural solver. The governing equations are discretized in space using spherical harmonics, yielding spectral integration accuracy. The predictions produced by this formulation are in good agreement with those obtained from simulations of spherical capsules in shear flows and optical tweezers experiments. The accuracy of the present method is evaluated using unbounded shear flow as a benchmark. Its computational cost grows proportional to p 5, where p is the degree of the spherical harmonic. It also exhibits a fast convergence rate that is approximately O ( p 6 ) for p ⪅ 20.

Flow residence time in intracranial aneurysms evaluated by in vitro 4D flow MRI.

The process of an intracranial aneurysm development, growth, and rupture is multifaceted and complex. In addition, clinical observations have identified the potential of thrombus formation within such aneurysms. While the underlying mechanism is not fully understood, the thrombi represent a potential risk factor for ischemic stroke. Emerging studies indicate that blood residence time (RT) is a promising hemodynamic metric associated with the aneurysm rupture and formation of intra-aneurysmal thrombi. Here, we present a methodology to experimentally evaluate both trajectory-wise and local RT based on magnetic resonance imaging (MRI) velocimetry, and apply it to in vitro flow measurements in scaled-up replicas of 9 patient-specific intracranial aneurysms. Lagrangian tracks of massless tracers are integrated from the velocity fields and averaged to return the mean RT in the aneurysm sac. This is found to be closely approximated by a simple time scale based on the sac diameter and space-time average of the aneurysmal fluid velocity. The mean RT is also correlated with the inflow time scale at the parent artery. These results also provide a basis for the estimation of RT when high-resolution hemodynamic maps are not available. With the continuous increase in accuracy and resolution enabled by progress in MRI technology, the methodology described here may in the future be applicable to in vivo data.

Comparison between fully resolved and time-averaged simulations of particle cloud dispersion produced by a violent expiratory event.

In this work we compare the DNS results (Fabregat et al. 2021, Fabregat et al. 2021) for a mild cough already reported in the literarure with those obtained with a compressible URANS equations with a k-ϵ turbulence model. In both cases, the dispersed phase has been modelled as spherical Lagrangian particles using the one-way coupling assumption. Overall, the URANS model is capable of reproducing the observed tendency of light particles under 64 µm in diameter to rise due to the action of the drag exerted by the buoyant puff generated by the cough. Both DNS and URANS found that particles above 64 µm will tend to describe parabolic trajectories under the action of gravitational forces. Grid independence analysis allows to qualify the impact of increasing mesh resolution on the particle cloud statistics as flow evolves. Results suggest that the k-ϵ model overpredicts the horizontal displacement of the particles smaller than 64 µm while the opposite occurs for the particles larger than 64 µm.

Fate of river-derived microplastics from the South China Sea: Sources to surrounding seas, shores, and abysses.

Microplastics (MPs) in the ocean have been widely recognized as causing global marine environmental problems. To gain a quantitative and comprehensive understanding of oceanic MP contamination, detailed numerical Lagrangian particle tracking experiments were conducted to evaluate the regional oceanic transport and dispersal of MPs in the South China Sea (SCS) derived from three major rivers, Pearl (China), Mekong (Vietnam), and Pasig (the Philippines), which are known to discharge large amounts of plastic waste into the SCS. As previous field surveys have suggested, MP contamination spreads from the surface to the deeper ocean in the water column, we thus considered three types of MPs: (1) positively buoyant (light) MPs, (2) positively buoyant (light) MPs with random walk diffusion, and (3) full 3-D tracking of non-buoyant MPs that are passively transported by ambient currents. Transport patterns of these MPs from the three rivers clearly showed the intra-annual variability associated with seasonally varying circulations driven by the Asian monsoons in the SCS. Many MPs floating during the prevailing southwest monsoon are transported to the northwest Pacific Ocean and the East China Sea through the Luzon Strait and the Taiwan Strait to form MP hotspots. Non-buoyant MPs are broadly transported from the surface layer to depths of approximately 100 m or deeper, where in situ observations are rare. In addition, the buoyant MPs drifting on the continental shelf originating from southern China tend to be pushed toward the shore and beached by northward wind-induced currents more pronouncedly than the non-buoyant MPs. Therefore, the river-derived MPs to the SCS were found to serve as sources to adjacent basins and oceans, to be distributed not only in the upper layer but also in the abyssal ocean (non-buoyant MPs), and to be transported to the shores (buoyant MPs).

Dyssynchronous Left Ventricular Activation is Insufficient for the Breakdown of Wringing Rotation.

Cardiac resynchronization therapy is a valuable tool to restore left ventricular function in patients experiencing dyssynchronous ventricular activation. However, the non-responder rate is still as high as 40%. Recent studies suggest that left ventricular torsion or specifically the lack thereof might be a good predictor for the response of cardiac resynchronization therapy. Since left ventricular torsion is governed by the muscle fiber orientation and the heterogeneous electromechanical activation of the myocardium, understanding the relation between these components and the ability to measure them is vital. To analyze if locally altered electromechanical activation in heart failure patients affects left ventricular torsion, we conducted a simulation study on 27 personalized left ventricular models. Electroanatomical maps and late gadolinium enhanced magnetic resonance imaging data informed our in-silico model cohort. The angle of rotation was evaluated in every material point of the model and averaged values were used to classify the rotation as clockwise or counterclockwise in each segment and sector of the left ventricle. 88% of the patient models (n = 24) were classified as a wringing rotation and 12% (n = 3) as a rigid-body-type rotation. Comparison to classification based on in vivo rotational NOGA XP maps showed no correlation. Thus, isolated changes of the electromechanical activation sequence in the left ventricle are not sufficient to reproduce the rotation pattern changes observed in vivo and suggest that further patho-mechanisms are involved.

The role of wind in controlling the connectivity of blue mussels (Mytilus edulis L.) populations.

BACKGROUND: Larval connectivity between distinct benthic populations is essential for their persistence. Although connectivity is difficult to measure in situ, it can be predicted via models that simulate biophysical interactions between larval behaviour and ocean currents. The blue mussel (Mytilus Edulis L.) is widespread throughout the northern hemisphere and extensively commercialised worldwide. In the Irish Sea, this industry represents ~ 50% of Welsh shellfisheries, where cultivation is mainly based on wild spat. However, the main sources and amount of spat varied interannually (1100 tonnes harvest in 2014 against zero in 2018). The aim of this study is to characterise the structure and dynamics of the blue mussel metapopulation within the northern part of the Irish Sea. METHODS: We develop a Lagrangian particle tracking model, driven by a high-resolution (from 30 to 5000 m) validated unstructured coastal hydrodynamic model of the Irish Sea, to simulate spatial and temporal variability of larval dispersal and connectivity between distinct mussel populations and potential settlement areas. RESULTS: Our results showed that: (1) larvae positioned near the surface were strongly influenced by wind-driven currents suggesting that connectivity networks had the potential to span hundreds of kilometres; (2) in contrast, larvae positioned deeper in the water column were driven by tidal currents, producing intricate spatial patterns of connectivity between mussel beds over tens of kilometres that were consistent over time. CONCLUSIONS: Dispersal of mussel larvae in the tidally energetic Irish Sea during the April-May spawning season is potentially driven by wind-driven surface currents, as confirmed by fisherman observations of inter-annual variability in wild spat collection. These results have important implications for metapopulation dynamics within the context of climate change and sustainable shellfisheries management (i.e. gain and loss of populations and harvest areas according to wind conditions).

In-Silico Conceptualisation of Continuous Millifluidic Separators for Magnetic Nanoparticles.

Magnetic nanoparticles are researched intensively not only for biomedical applications, but also for industrial applications including wastewater treatment and catalytic processes. Although these particles have been shown to have interesting surface properties in their bare form, their magnetisation remains a key feature, as it allows for magnetic separation. This makes them a promising carrier for precious materials and enables recovery via magnetic fields that can be turned on and off on demand, rather than using complex (nano)filtration strategies. However, designing a magnetic separator is by no means trivial, as the magnetic field and its gradient, the separator dimensions, the particle properties (such as size and susceptibility), and the throughput must be coordinated. This is showcased here for a simple continuous electromagnetic separator design requiring no expensive materials or equipment and facilitating continuous operation. The continuous electromagnetic separator chosen was based on a current-carrying wire in the centre of a capillary, which generated a radially symmetric magnetic field that could be described using cylindrical coordinates. The electromagnetic separator design was tested in-silico using a Lagrangian particle-tracking model accounting for hydrodynamics, magnetophoresis, as well as particle diffusion. This computational approach enabled the determination of separation efficiencies for varying particle sizes, magnetic field strengths, separator geometries, and flow rates, which provided insights into the complex interplay between these design parameters. In addition, the model identified the separator design allowing for the highest separation efficiency and determined the retention potential in both single and multiple separators in series. The work demonstrated that throughputs of ~1/4 L/h could be achieved for 250-500 nm iron oxide nanoparticle solutions, using less than 10 separator units in series.

Investigating Pseudo-nitzschia australis introduction to the Gulf of Maine with observations and models.

In 2016, an unprecedented Pseudo-nitzschia australis bloom in the Gulf of Maine led to the first shellfishery closures due to domoic acid in the region's history. In this paper, potential introduction routes of P. australis are explored through observations, a hydrodynamic model, and a Lagrangian particle tracking model. Based on particle tracking experiments, the most likely source of P. australis to the Gulf of Maine was the Scotian Shelf. However, in 2016, connectivity between the Scotian Shelf and the bloom region was not significantly different from the other years between 2012 and 2019, nor were temperature conditions more favorable for P. australis growth. Observations indicated changes on the Scotian Shelf in 2016 preceded the introduction of P. australis: increased bottom salinity and decreased surface salinity. The increased bottom salinity on the shelf may be linked to anomalously saline water observed near the coast of Maine in 2016 via transport through Northeast Channel. The changes in upstream water mass properties may be related to the introduction of P. australis, and could be the result of either increased influence of the Labrador Current or increased outflow from the Gulf of St. Lawrence. The ultimate source of P. australis remains unknown, although the species has previously been observed in the eastern North Atlantic, and connectivity across the ocean is possible via a subpolar route. Continued and increased monitoring is warranted to track interannual Pseudo-nitzschia persistence in the Gulf of Maine, and sampling on the Scotian Shelf should be conducted to map upstream P. australis populations.