Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lagrangian Particles”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Exploring ship track spreading rates with a physics-informed Langevin particle parameterization

Abstract. The rate at which aerosols spread from a point source injection, such as from a ship or other stationary pollution source, is critical for accurately representing subgrid plume spreading in a climate model. Such climate model results will guide future decisions regarding the feasibility and application of large-scale intentional marine cloud brightening (MCB). Prior modeling studies have shown that the rate at which ship plumes spread may be strongly dependent on meteorological conditions, such as precipitating versus non-precipitating boundary layers and shear. In this study, we apply a Lagrangian particle model (PM-ABL v1.0), governed by a Langevin stochastic differential equation, to create a simplified framework for predicting the rate of spreading from a ship-injected aerosol plume in sheared, precipitating, and non-precipitating boundary layers. The velocity and position of each stochastic particle is predicted with the acceleration of each particle being driven by the turbulent kinetic energy, dissipation rate, momentum variance, and mean wind. These inputs to the stochastic particle velocity equation are derived from high-fidelity large-eddy simulations (LES) equipped with a prognostic aerosol–cloud microphysics scheme (UW-SAM) to simulate an aerosol injection from a ship into a cloud-topped marine boundary layer. The resulting spreading rate from the reduced-order stochastic model is then compared to the spreading rate in the LES. The stochastic particle velocity representation is shown to reasonably reproduce spreading rates in sheared, precipitating, and non-precipitating cases using domain-averaged turbulent statistics from the LES.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Efficient bi-directional coupling of 3D computational fluid-particle dynamics and 1D Multiple Path Particle Dosimetry lung models for multiscale modeling of aerosol dosimetry

The development of predictive aerosol dosimetry models has been a major focus of environmental toxicology and pharmaceutical health research for decades. One-dimensional (1D) models successfully predict overall deposition averages but fail to accurately predict local deposition. Computational fluid-particle dynamics (CFPD) models provide site-specific predictions but at a computational cost that prohibits whole lung predictions. Thus, there is a need for developing multiscale strategies to provide a realistic subject-specific picture of the fate of inhaled aerosol in the lungs. CT-based 3D/CFPD models of the large airways were bidirectionally coupled with individualized 1D Navier-Stokes airflow and particle transport based upon the widely used Multiple Path Particle Dosimetry Model (MPPD). Distribution of airflows among lobes was adjusted by measured lobar volume changes observed in CT images between FRC and FRC + 1.5 L. Additionally, as a test of the effectiveness of the coupling procedures, deposition modeling of previous 1 µm aerosol exposure studies was performed. The complete coupled model was run for 3 breaths, with the computation-intense portion being the 3D CFD Lagrangian particle tracking calculation. The average deposition per breath was 11% in the combined multiscale model with site-specific doses available in the CFPD portion of the model and airway- or region-specific deposition available for the MPPD portion. In conclusion, the key methods developed in this study enable predictions of ventilation heterogeneities and aerosol deposition across the lungs that are not captured by 3D or 1D models alone. Overall, these methods can be used as the foundation for multi-scale modeling of the full respiratory system.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Simultaneous 10 kHz three-dimensional CH 2 O and tomographic PIV measurements in a lifted partially-premixed jet flame

High-speed, three-dimensional (3D) scalar-velocity field measurements are demonstrated in a lifted partially-premixed dimethyl-ether/air jet flame using simultaneous laser-induced fluorescence (LIF) of formaldehyde (CH 2 O) and tomographic particle image velocimetry (TPIV). The 3D LIF measurements are performed by raster scanning the laser beam from a 100 kHz pulse-burst laser across the probe volume using an acousto-optic deflector. The volumetric reconstruction of the LIF signal from ten parallel planes provides quasi-instantaneous 3D LIF measurements that are synchronized with 10 kHz TPIV measurements. Here, the temporally resolved CH 2 O-LIF and velocity field data are used to analyze Lagrangian particle trajectories and displacement speeds at the base of the lifted flame. The particle trajectories reveal flow structures that are difficult to observe in an Eulerian reference frame. Positive and negative displacement speeds are observed at the CH 2 O-LIF surfaces at the inner and outer regions of the jet flame with a maximum displacement speed of approximately eight times the laminar flame speed of a stoichiometric dimethyl-ether/air mixture. The observations in present study support formation of lean premixed flames near the flame base whose role in flame stabilization should be investigated in more detail.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Going With the Flow: Using a Particle Dispersion Model to Reveal Links Between Ocean Ecosystems and Marine Aerosols

Marine biogenic particle contributions to atmospheric aerosol concentrations are not well understood though they are important for determining cloud optical and cloud nucleating properties. Here we examine the relationship between marine aerosol measurements with satellite and model fields of ocean biology and meteorological variables during the North Atlantic Aerosols and Marine Ecosystems Study (NAAMES). NAAMES consisted of four field campaigns between November 2015 and April 2018 that aligned with the four major phases of the annual phytoplankton bloom cycle. The FLEXPART Lagrangian particle dispersion model is used to connect these variables spatiotemporally to ship-based aerosol. We find that correlations between some aerosol measurements with satellite measured and modelled variables increase with increasing trajectory length. This indicates biological and meteorological processes over the air mass history are influential to measured particle properties. For example, results show organic aerosol mass is driven by biogenic volatile organic compound (VOC) emissions from photosynthesis by phytoplankton stocks during advection into the region. In contrast, primary marine aerosol (PMA) concentrations showed better correlations to variables weighted with shorter air mass residence times, which reflects their localized origin as primary emissions. Evidence also suggest variability in PMA concentrations likely play an important role in VOC emission removal. The removal of VOCs lessen their impact on particle number and mass concentrations. The sensitivity of marine aerosols to the main ecosystem-atmosphere drivers and its dependence on air mass history are to be further discussed.

Kevin J Sanchez↗
DOE OSTI2021

A Lagrangian Study of Interfaces at the Edges of Cumulus Clouds

Abstract Interfaces at the edge of an idealized, nonprecipitating, warm cloud are studied using direct numerical simulation (DNS) complemented with a Lagrangian particle tracking routine. Once a shell has formed, four zones can be distinguished: the cloud core, visible shell, invisible shell, and the environment. The union of the visible and invisible regions is the shell commonly referred to in literature. The boundary between the invisible shell and the environment is the turbulent–nonturbulent interface (TNTI), which is typically not considered in cloud studies. Three million particles were seeded homogeneously across the domain and properties were recorded along individual trajectories. The results demonstrate that the traditional cloud boundary (separating cloudy and noncloudy regions using thresholds applied on liquid condensate or updraft velocity) are some distance away from the TNTI. Furthermore, there is no dynamic difference between the traditional liquid-condensate boundary and the region extending to the TNTI. However, particles crossing the TNTI exhibit a sharp jump in enstrophy and a smooth increase in buoyancy. The traditional cloud boundary coincides with the location of minimum buoyancy in the shell. The shell premixes the entraining and detraining air and analysis reveals a highly skewed picture of entrainment and detrainment at the traditional cloud boundary. A preferential entrainment of particles with velocity and specific humidity higher than the mean values in the shell is observed. Large-eddy simulation of a more realistic setup detects an interface with similar properties using the same thresholds as in the DNS, indicating that the DNS results extrapolate beyond their idealized conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Characterizing Reactive Transport Behavior in a Three-Dimensional Discrete Fracture Network

While several studies have linked network and in-fracture scale properties to conservative transport behavior in subsurface fractured media, studies on reactive transport cases remain relatively underdeveloped. In this study, we explore the behavior of an irreversible kinetic reaction during the interaction of two solute plumes, one consisting of species A and the other species B. When the plumes converge, these species react kinetically to form a new species C via A+B→kC. This reactive system is studied using a three-dimensional discrete fracture network (DFN) model coupled with reactive Lagrangian particle tracking. We find that the interplay of network topology and chemical properties of the reactive solutes controls reactive transport processes. The network topology drives species A and B together, and the chemical properties dictate whether and how quickly a reaction occurs. Results demonstrate that reactions are most likely to occur in high-velocity fractures that make up the network backbone. The interplay between species’ chemical properties and transport is characterized by a non-dimensional Damköhler (Da) number. We show that the spatial distribution of reactions is sensitive to Da, which subsequently influences late-time tailing behavior in outlet breakthrough time distributions. The results of this study provide initial insights into how an irreversible reaction occurs during transport in a fracture network, using a methodology that can be applied to study reactive transport in a wide range of fractured media environments and contexts.

58 GEOSCIENCES↗
DOE OSTI2022

Explicit quantification of residence and flushing times in the Salish Sea using a sub-basin scale shoreline resolving model

The Salish Sea, located in the Pacific Northwest region of North America has complex currents and circulation features distributed over numerous interconnected deep basins with islands. Increased risk of exposure to oil spills and untreated wastewater from maritime emergencies and treatment plant failures have led to a need for quantifying residence and flushing characteristics at a sub-basin scale using a shoreline resolving hydrodynamic model. An unstructured grid model of the Salish Sea was developed using FVCOM with a ≈75-100m shoreline resolution. In addition to 7 tides and 23 salinity and temperature monitoring stations, an extensive currents data set from 135 stations collected over a span of three years was used for skill assessment and validation. Explicit forward computations were then conducted to define and quantify residence and flushing times in various sub-basins of interest using (a) Lagrangian particles and (b) Numerical/virtual dye experiments. The results in most basins show expected seasonal variability with longer flushing time associated with summer lower tides and lower freshwater inflows. However, contrary to expectation, flushing time is significantly longer in wintertime in large fjord-like basins such as Hood Canal (≈138 days), likely due to increased stratification and reduced mixing. The flushing time for the Puget Sound region of the Salish Sea is ≈115 days, while Georgia Basin is 240 days when analyzed as stand-alone basins with zero background concentrations. When examined as part of the flushing of the entire Puget Sound filled with virtual dye, the compounded flushing times for embedded sub-basins can be significantly longer in order of magnitudes and largely dictated by the flushing time of Puget Sound. The computed residence and flushing time scales tabulated over 36 sub-basins provide an improved understanding of water renewal in the system, informing pollution management actions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Comparison of atmospheric radionuclide dispersion models for a risk-informed consequence-driven advanced reactor licensing framework

Current nuclear facility emergency planning zones (EPZs) are based on outdated distance-based criteria, predating comprehensive dose and risk-informed frameworks. Recent advancements in simulation tools have permitted the development of site-specific, dose, and risk-based consequence-driven assessment frameworks. This study investigated the computation of advanced reactor (AR) EPZs using two atmospheric dispersion models: a straight-line Gaussian plume model (GPM) and a semi-Lagrangian Particle in Cell (PIC). Two case studies were conducted: (1) benchmarking the NRC SOARCA study for the Peach Bottom Nuclear Generating Station and (2) analyzing an advanced INL Heat Pipe Design A microreactor's end-of-cycle inventory. The dose criteria for both cases were 10 mSv at mean weather conditions and 50 mSv at 95th percentile weather conditions at 96 h post-release. Results demonstrated that GPM and PIC estimated similar mean peak dose levels for large boiling water reactors in the farfield case, placing EPZ limits beyond current regulations. For ARs with source terms remaining in the nearfield, PIC modeling without specific nearfield considerations could result in excessively high doses and inaccurate EPZ designations. PIC dispersion demonstrated an order of magnitude higher estimate of nearfield inhalation dose contribution when compared to GPM results. Furthermore, both models significantly reduced EPZ sizing within the nearfield. Thus, reductions in the AR source term may eliminate the need for a separate EPZ.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2022

SimELIT: A Novel GUI-Based Comprehensive Ion Trajectory Simulation Software for Mass Spectrometry

Ion trajectory simulation in mass spectrometry systems from injection to detection is technically challenging but very important for better understanding the ion dynamics in instrument development. Here, in this work, we present SimELIT (Simulator of Eulerian and Lagrangian Ion Trajectories), a novel ion trajectory simulation platform. SimELIT is built upon a suite of multiphysics solvers compiled into OpenFOAM (an open-source numerical solver library particularly used for computational mechanics), with a simple web-based graphical user interface (GUI) allowing users to define the details of OpenFOAM cases and run simulations. SimELIT is a modular program and can provide extensions of physics (e.g., gas flows, electrodynamic fields) and thus enable ion trajectory simulations from the ion source to detector. The current version (SimELIT) provides two numerical solvers for ion trajectory simulations–(1) a Lagrangian particle tracker in vacuum and (2) a Eulerian ion density solver in background gas in the presence of electric fields. Here, we describe the architecture of SimELIT, including its use of Docker and the React Framework, and demonstrate the computation of ion trajectories of multiple m/z values in a static/linear voltage drop in vacuum (across a 1 m long flight tube). Further, the drift motion of ions under 1 Torr pressure conditions in a static background (N 2 ) gas through a 20 V/cm static electric field is shown. The results produced from SimELIT were compared with SIMION and theoretical estimates. In addition, we report the computation of ion trajectories in electrodynamic fields within a planar FAIMS device operating at atmospheric pressure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The Influence of Ocean Topography on the Upwelling of Carbon in the Southern Ocean

The physical circulation of the Southern Ocean sets the surface concentration and thus air-sea exchange of CO 2 . However, we have a limited understanding of the three-dimensional circulation that brings deep carbon-rich waters to the surface. Here, we introduce and analyze a novel high-resolution ocean model simulation with active biogeochemistry and online Lagrangian particle tracking. We focus our attention on a subset of particles with high dissolved inorganic carbon (DIC) that originate below 1,000 m and eventually upwell into the near-surface layer (upper 200 m). We find that 71% of the DIC-enriched water upwelling across 1,000 m is concentrated near topographic features, which occupy just 33% of the Antarctic Circumpolar Current. Once particles upwell to the near-surface layer, they exhibit relatively uniform pCO 2 levels and DIC decorrelation timescales, regardless of their origin. Furthermore, our results show that Southern Ocean bathymetry plays a key role in delivering carbon-rich waters to the surface.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Air‐Sea CO 2 Fluxes Localized by Topography in a Southern Ocean Channel

Abstract Air‐sea exchange of carbon dioxide (CO 2 ) in the Southern Ocean plays an important role in the global carbon budget. Previous studies have suggested that flow around topographic features of the Southern Ocean enhances the upward supply of carbon from the deep to the surface, influencing air‐sea CO 2 exchange. Here, we investigate the role of seafloor topography on the transport of carbon and associated air‐sea CO 2 flux in an idealized channel model. We find elevated CO 2 outgassing upstream of a seafloor ridge, driven by anomalous advection of dissolved inorganic carbon. Argo‐like Lagrangian particles in our channel model sample heterogeneously in the vicinity of the seafloor ridge, which could impact float‐based estimates of CO 2 flux.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Dust supply to close binary systems

Binary systems can be born surrounded by circumbinary discs. The gaseous disc around either of the two stellar companions can have its life extended by the supply of mass arriving from the circumbinary disc. The objective of this study is to investigate the gravitational interactions exerted by a compact and eccentric binary system on the circumbinary and circumprimary discs, and the resulting transport of gas and solids between the disc components. We assume that the gas in the system behaves as a fluid, and we model its evolution by means of high-resolution hydrodynamical simulations. Dust grains are modelled as Lagrangian particles that interact with the gas and the stars. Our models indicate that significant fluxes of gas and dust proceed from the circumbinary disc towards the circumprimary disc. For the applied system parameters, grains of certain sizes are segregated outside the tidal gap generated by the stars. Consequently, the size distribution of the transported dust is not continuous, but presents a gap in the millimetre size range. In close binaries, the lifetime of an isolated circumprimary disc is found to be short, ~10 5 years, because of its low mass. However, because of the influx of gas from beyond the tidal gap, the disc around the primary star can survive much longer, ~10 6 years, as long as gas accretion from the circumbinary disc continues. The supply of solids and the extended lifetime of a circumbinary disc also aids in the possible formation of giant planets. Compared to close binary systems without a circumbinary disc, we expect a higher frequency of singleplanet or multiple-planet systems. Additionally, a planetesimal or debris belt can form in the proximity of the truncation radius of the circumprimary disc and/or around the location of the exterior edge of the tidal gap.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Necessary and sufficient conditions for quasisymmetry

A necessary and sufficient set of conditions for a quasisymmetric magnetic field in the form of constraint equations is derived from first principles. Without any assumption regarding the magnetohydrodynamic (MHD) equilibrium of the plasma, conditions for quasisymmetry are constructed starting from the single-particle Lagrangian to the leading order. The conditions presented in the paper are less restrictive than the set recently obtained by Burby et al. [“Some mathematics for quasi-symmetry,” arXiv:1912.06468 (2019)], and could facilitate ongoing efforts toward investigating the existence of global quasisymmetric MHD equilibria. Finally, it is also shown that quasisymmetry implies the existence of flux surfaces, regardless of whether the field corresponds to an MHD equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Measurement and simulation of small cryogenic neon pellet Ne-I 640 nm photon efficiency during ablation in DIII-D plasma

Small (~1 mm) neon pellet fragments are fired into DIII-D H-mode plasmas, and resulting trajectory-averaged photon efficiencies (neutral neon ionization events for every photon emitted) of S/XB ≈ 85 are estimated for Ne-I 640 nm by dividing the estimated initial pellet fragment mass by the measured number of emitted Ne-I photons. The experiments are modeled by running the Lagrangian particle (LP) fluid/magneto-hydrodynamic pellet code to estimate axial ablation plume neon density profiles and temperature profiles at each pellet position. These solutions are then fed into the PrismSPECT collisional-radiative code, which calculates resulting neon charge states and photon emission rates, giving a profile-average of S/XB ≈ 109. The burnthrough plasma minor radius predicted by LP (ρ ≈ 0.63) is reasonably close to the experimental observation ρ ≈ 0.6. Furthermore, the modeling indicates that local S/XB is not constant along the pellet trajectory but tends to increase with increasing ablation rate. Non-equilibrium kinetics are predicted to be very important, while line trapping is predicted to be relatively unimportant (for Ne-I 640 nm S/XB).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Role of electrostatic perturbation on kinetic resistive wall mode with application to spherical tokamak

Abstract A more complete non-perturbative magnetohydrodynamic (MHD)-kinetic hybrid formulation is developed by including the perturbed electrostatic potential δφ in the particle Lagrangian. The fluid-like counter-parts of the hybrid equations, in the Chew-Goldberger-Low high-frequency limit, are also derived and utilized to test the new toroidal implementation in the MARS-K code. Application of the updated non-perturbative hybrid model for a high- β spherical tokamak plasma in MAST finds that the perturbed electrostatic potential generally plays a minor role in the n = 1 ( n is the toroidal mode number) resistive wall mode instability. The effect of δφ is largely destabilizing, with the growth rate of the instability increased by several (up to 20) percent as compared to the case without including δφ . A similar relative change is also obtained for the kinetic-induced resonant field amplification effect at high- β in the MAST plasma considered. The updated capability of the MARS-K code allows quantitative exploration of drift kinetic effects on various MHD instabilities and the antenna-driven plasma response where the electrostatic perturbation, coupled to magnetic perturbations, may play important roles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Estimated impacts of forest restoration scenarios on smoke exposures among outdoor agricultural workers in California

Abstract As wildfires continue to worsen across western United States, forest managers are increasingly employing prescribed burns as a way to reduce excess fuels and future wildfire risk. While the ecological benefits of these fuel treatments are clear, little is known about the smoke exposure tradeoffs of using prescribed burns to mitigate wildfires, particularly among at-risk populations. Outdoor agricultural workers are a population at increased risk of smoke exposure because of their time spent outside and the physical demands of their work. Here, we assess the smoke exposure impacts among outdoor agricultural workers resulting from the implementation of six forest management scenarios proposed for a landscape in the Central Sierra, California. We leverage emissions estimates from LANDIS-II to model daily PM 2.5 concentrations with the Hybrid Single Particle Lagrangian Integrated Trajectory Model (HYSPLIT) and link those to agricultural employment data from the Bureau of Labor Statistics. We find a u-shaped result, in that moderate amounts of prescribed burning result in the greatest reduction in total smoke exposure among outdoor agricultural workers, particularly during months of peak agricultural activity due to wildfire-specific smoke reductions. The reduction in total smoke exposure, relative to scenarios with the least amount of management, decreases as more prescribed burning is applied to the landscape due to the contributions of the fuel treatments themselves to overall smoke burden. The results of this analysis may contribute to preparedness efforts aimed at reducing smoke exposures among outdoor agricultural workers, while also informing forest management planning for this specific landscape.

Environmental Sciences & Ecology↗
DOE OSTI2021

Second-generation dust in planetary systems: the case of HD 163296

Observations indicate that large, dust-laden protoplanetary discs are common. Some features, like gaps, rings, and spirals, suggest they may host young planets, which can excite the orbits of nearby leftover planetesimals. Energetic collisions among these bodies can lead to the production of second-generation dust. Grains produced by collisions may have a dynamical behaviour different from that of first-generation, primordial dust out of which planetesimals and planets formed. Here, we aim to study these differences for the HD 163296 system and determine whether dynamical signatures in the mixture of the two dust populations can help separate their contributions. We use three-dimensional (3D) hydrodynamic models to describe the gaseous disc with three, Saturn- to Jupiter-mass, embedded planets. Dust grains, of sizes $1\, \mu \mathrm{m}$–$1\, \mathrm{mm}$, are treated as Lagrangian particles with resolved thermodynamics and mass-loss. Initial disc and planet configurations are derived from observation-based work, which indicates low gas viscosity. The 3D approach also allows us to detect the formation of vortices induced by Rossby waves, where dust becomes concentrated and may contribute to planetesimal formation. We find that the main differences in the dynamical behaviour of first- and second-generation dust occur in the vertical distribution. The two populations have similar distributions around the disc mid-plane, although second-generation dust shows longer residence times close to the radial locations of the planets’ gas gaps. Sedimentation rates of $\mu$m-sized grains are comparable to or lower than the production rates by planetesimals’ collisions, making this population potentially observable. These outcomes can be extended to similar systems harbouring giant planets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Discrete ion stochastic continuum overdamped solvent algorithm for modeling electrolytes

In this paper we develop a methodology for the mesoscale simulation of strong electrolytes. The methodology is an extension of the fluctuating immersed-boundary approach that treats a solute as discrete Lagrangian particles that interact with Eulerian hydrodynamic and electrostatic fields. In both algorithms the immersed-boundary method of Peskin is used for particle-field coupling. Hydrodynamic interactions are taken to be overdamped, with thermal noise incorporated using the fluctuating Stokes equation, including a "dry diffusion" Brownian motion to account for scales not resolved by the coarse-grained model of the solvent. Long-range electrostatic interactions are computed by solving the Poisson equation, with short-range corrections included using an immersed-boundary variant of the classical particle-particle particle-mesh technique. Also included is a short-range repulsive force based on the Weeks-Chandler-Andersen potential. This methodology is validated by comparison to Debye-Hückel theory for ion-ion pair correlation functions, and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In each case, good agreement is observed, provided that hydrodynamic interactions at the typical ion-ion separation are resolved by the fluid grid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗