The Zapiola Anticyclone: A Lagrangian study of its kinematics in an eddy-permitting ocean model
Not Available
54 ENVIRONMENTAL SCIENCES↗
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.
A continuum meshfree method for sandbox-style numerical modeling of accretionary and doubly vergent wedges
Here, in this study, we introduce and utilize GEOSPH, a Lagrangian particle-based continuum Smoothed Particle Hydrodynamics (SPH) code to simulate the deformation of geologic materials in sandbox-style experiments using an elastoplastic constitutive model. Due to the meshfree nature of SPH, our proposed framework is capable of handling the large deformations characteristic of sandbox-style experiments. SPH is also able to capture strain localization and discontinuities, allowing us to resolve faulting patterns emerging in the sandbox experiments. We successfully validate our SPH implementation against the benchmark experiments of Buiter and co-workers, and then apply SPH to understand the deformational processes occurring in doubly vergent orogens resultant from the S-point experimental setup and from finitely high wedge-like indenter backstops of different inclinations. We find that the topographic slope of the prowedge results from the competing processes of vertical uplift due to the indenter and forward thrusting, while the topographic slope of the retrowedge is a product of repeated shallow slope failure. Our findings show that SPH is a promising numerical method which should be added to the current palette of computational tools used by structural geologists for solving problems relating to tectonics and crustal deformation.
58 GEOSCIENCES↗
Are turbulence effects on droplet collision–coalescence a key to understanding observed rain formation in clouds?
Rain formation is a critical factor governing the lifecycle and radiative forcing of clouds and therefore it is a key element of weather and climate. Cloud microphysics–turbulence interactions occur across a wide range of scales and are challenging to represent in atmospheric models with limited resolution. Based on past experiments and idealized numerical simulations, it has been postulated that cloud turbulence accelerates rain formation by enhancing drop collision–coalescence. We provide substantial evidence for significant impacts of turbulence on the evolution of cloud droplet size distributions and rain formation by comparing high-resolution observations of cumulus congestus clouds with state-of-the-art large-eddy simulations coupled with a Lagrangian particle-based microphysics scheme. Turbulent coalescence must be included in the model to accurately represent the observed drop size distributions, especially for drizzle drop sizes at lower heights in the cloud. Turbulence causes earlier rain formation and greater rain accumulation compared to simulations with gravitational coalescence only. The observed rain size distribution tail just above cloud base follows a power law scaling that deviates from theoretical scalings considering either a purely gravitation collision kernel or a turbulent kernel neglecting droplet inertial effects, providing additional evidence for turbulent coalescence in clouds. In contrast, large aerosols acting as cloud condensation nuclei (“giant CCN”) do not significantly impact rain formation owing to their long timescale to reach equilibrium wet size relative to the lifetime of rising cumulus thermals. Overall, turbulent drop coalescence exerts a dominant influence on rain initiation in warm cumulus clouds, with limited impacts of giant CCN.
54 ENVIRONMENTAL SCIENCES↗
Interpretative 3D MHD modelling of deuterium SPI into a JET H-mode plasma
Abstract The pre-thermal quench (pre-TQ) dynamics of a pure deuterium ( D 2 ) shattered pellet injection (SPI) into a 3 MA / 7 MJ JET H-mode plasma is studied via 3D non-linear MHD modelling with the JOREK code. The interpretative modelling captures the overall evolution of the measured density and radiated power. The simulations also identify the importance of the drifts of ablation plasmoids towards the tokamak low field side (LFS) and the impurities in the background plasma in fragment penetration, assimilation, radiative cooling and MHD activity in D 2 SPI experiments. It is found that plasmoid drifts lead to an about 70% reduction of the central line-integrated density (compared to a simulation without drifts) in the JET D 2 SPI discharge considered. Impurities that pre-exist before the SPI as well as those from possible impurity influxes related to the SPI are shown to dominate the radiation in the considered discharge. With inputs from JOREK simulations, modelling with the Lagrangian particle-based pellet code PELOTON reproduces the deviation of the SPI fragments in the direction of the major radius as observed by the fast camera. This confirms the role of rocket effects and plasmoid drifts in the considered discharge and reinforces the validity of the JOREK modelling. The limited core density rise due to plasmoid drifts and the strong radiative cooling and MHD activity with impurities (depending on their species and concentration) could limit the effectiveness of LFS D 2 SPI in runaway electron avoidance and are worth considering in the design of the ITER disruption mitigation system.
Physics↗
r -process nucleosynthesis and kilonovae from hypermassive neutron star post-merger remnants
ABSTRACT We investigate r-process nucleosynthesis and kilonova emission resulting from binary neutron star (BNS) mergers based on a three-dimensional (3D) general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. The simulation includes a microphysical finite-temperature equation of state (EOS) and neutrino emission and absorption effects via a leakage scheme. We track the thermodynamic properties of the ejecta using Lagrangian tracer particles and determine its composition using the nuclear reaction network SkyNet. We investigate the impact of neutrinos on the nucleosynthetic yields by varying the neutrino luminosities during post-processing. The ejecta show a broad distribution with respect to their electron fraction Ye, peaking between ∼0.25–0.4 depending on the neutrino luminosity employed. We find that the resulting r-process abundance patterns differ from solar, with no significant production of material beyond the second r-process peak when using luminosities recorded by the tracer particles. We also map the HMNS outflows to the radiation hydrodynamics code SNEC and predict the evolution of the bolometric luminosity as well as broadband light curves of the kilonova. The bolometric light curve peaks on the timescale of a day and the brightest emission is seen in the infrared bands. This is the first direct calculation of the r-process yields and kilonova signal expected from HMNS winds based on 3D GRMHD simulations. For longer-lived remnants, these winds may be the dominant ejecta component producing the kilonova emission.
79 ASTRONOMY AND ASTROPHYSICS↗
Mixing by internal gravity waves in stars: assessing numerical simulations against theory
ABSTRACT Here we present a study of radial chemical mixing in non-rotating massive main-sequence stars driven by internal gravity waves (IGWs), based on multidimensional hydrodynamical simulations with the fully compressible code MUSIC. We examine two proposed mechanisms of material mixing in stars by IGWs that are commonly quoted, relating to thermal diffusion and sub-wavelength shearing. Thermal diffusion provides a non-restorative effect to the waves, leaving material displaced from its previous equilibrium, while shearing arising within the waves drives weak localized flows, mixing the fluid there. Using IGW spectra from the simulations, we evaluate theoretical predictions of mixing rates due to these mechanisms. We show, for $20\, \mathrm{M}_\odot$ main-sequence stars, that neither of these mechanisms are likely to create mixing sufficient to correct inaccuracies in current stellar evolution models. Furthermore, we compare these predictions to results obtained from Lagrangian tracer particles, following a method recently used for global simulations of stellar interiors to measure mixing by IGWs in their radiative zones. We demonstrate that tracer particle methods face significant numerical challenges in measuring the small diffusion coefficients predicted by the aforementioned theories, for which they are prone to yielding artificially enhanced coefficients. Diffusion coefficients based on such methods are currently used with stellar evolution codes for asteroseismic studies, but should be viewed with caution. Finally, in a case where tracer particles do not suffer from numerical artefacts, we suggest that a diffusion model is not suitable for time-scales typically considered by 2D numerical simulations.
79 ASTRONOMY AND ASTROPHYSICS↗
Fully kinetic simulations of strong steady-state collisional planar plasma shocks
We report on simulations of strong, steady-state collisional planar plasma shocks with fully kinetic ions and electrons, independently confirmed by two fully kinetic codes (an Eulerian continuum and a Lagrangian particle-in-cell). While kinetic electrons do not fundamentally change the shock structure as compared with fluid electrons, we find an appreciable rearrangement of the preheat layer, associated with nonlocal electron heat transport effects. The electron heat-flux profile qualitatively agrees between kinetic- and fluid-electron models, suggesting a certain level of “stiffness,” though substantial nonlocality is observed in the kinetic heat flux. We also find good agreement with nonlocal electron heat-flux closures proposed in the literature. Finally, in contrast to the classical hydrodynamic picture, we find a significant collapse in the “precursor” electric-field shock at the preheat layer leading edge, which correlates with the electron-temperature gradient relaxation.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A Hybrid Lagrangian/Eulerian Collocated Velocity Advection and Projection Method for Fluid Simulation
Abstract We present a hybrid particle/grid approach for simulating incompressible fluids on collocated velocity grids. Our approach supports both particle‐based Lagrangian advection in very detailed regions of the flow and efficient Eulerian grid‐based advection in other regions of the flow. A novel Backward Semi‐Lagrangian method is derived to improve accuracy of grid based advection. Our approach utilizes the implicit formula associated with solutions of the inviscid Burgers’ equation. We solve this equation using Newton's method enabled by C 1 continuous grid interpolation. We enforce incompressibility over collocated, rather than staggered grids. Our projection technique is variational and designed for B‐spline interpolation over regular grids where multiquadratic interpolation is used for velocity and multilinear interpolation for pressure. Despite our use of regular grids, we extend the variational technique to allow for cut‐cell definition of irregular flow domains for both Dirichlet and free surface boundary conditions.
Gagniere, S.↗
ORNL/grit
Grit is a software library designed specifically for conducting particle-based Lagrangian simulations with CPU/GPUperformance portability. This library allows researchers and engineers to perform a wide range of simulations, including multiphase flow, and more. Grit employs Message Passing Interface (MPI) for distributed memory parallelism and Kokkos programming model for on-node shared memory parallelism with performance portability across different architectures of GPUs and multi-core/manycore CPUs.
Ge, Wenjun [Oak Ridge National Laboratory (ORNL), ↗
Lagrangian-Eulerian Simulations of Fission Particle Transport in a Full-Scale Ventilation Duct in Nonreactor Nuclear Facilities
Explore the source record for details and available documents.
Nguyen, Thien D. [ORNL] (ORCID:0000000159722006)↗
Hydrodynamic Analysis and Optimization of Aquantis Marine Turbine: Cooperative Research and Development (Final Report)
The primary aim of this proposal is to improve the accurate prediction of hydrodynamic performance and dynamic load responses of the AQ10 floating axial-flow tidal turbine with a tri-cat mooring configuration. The validation of reduced-order modeling approaches with high-fidelity model will be implemented. Additionally, the frequency response domain, Response Amplitude Floating Wind (RAFT) toolbox plus an optimizer expanded for marine hydrokinetic turbines under the Submarine Hydrokinetic And Riverine Kilo-megawatt. Systems (SHARKS) program will be used for designing and exploring different key design parameters (platform dimension, mooring layout and its parameters) of next marine hydrokinetic (MHK) turbine generation.
16 TIDAL AND WAVE POWER↗
Going Off Grid: A Comparative Study of the Lagrangian and Eulerian Perspectives of New Particle Formation Events
New particle formation and growth (NPF&G) is the process by which ultrafine particles are formed from gas-phase precursors. NPF&G is the dominant source of global aerosol number with important influences on climate. Most observations of NPF&G events are conducted at stationary sites; however, NPF&G observed from stationary sites is influenced by gradual or rapid changes in the air masses passing over the site, complicating NPF&G analysis. In this work, we use observations and a 3D aerosol model to compare aerosol size distributions at a stationary site (Southern Great Plains [SGP] observatory, Oklahoma, USA) and along Lagrangian trajectories crossing the site. The model simulates the NPF&G events reasonably well at SGP. Using the model to compare the Lagrangian and stationary perspectives, we can explain previously unanalyzable days with some evidence of NPF&G as either non-event or analyzable NPF&G days. We find most of the unanalyzable NPF&G days are due to isolated and inhomogeneous NPF&G occurring upwind of the stationary site, often in the outflow of urban regions. Finally, we compare formation rates of 3 nm particles, growth rates, and the survival probability of 3 nm particles growing to 25 nm between the stationary and Lagrangian perspectives. Because of the much larger number of analyzable days along the Lagrangian trajectories, this perspective potentially provides more robust statistics and better characterization of NPF&G event extremes. Our method for extracting chemical/physical properties along Lagrangian trajectories from 3D models can be applied to a wide range of science questions.
O’Donnell, Samuel E. [Colorado State Univ., Fort C↗
Lagrangian conditional statistics and flow topology in edge plasma turbulence
Lagrangian statistics and particle transport in edge plasma turbulence are investigated using the Hasegawa–Wakatani model and its modified version. The latter shows the emergence of pronounced zonal flows. Different values of the adiabaticity parameter are considered. The main goal is to characterize the role of coherent structures, i.e., vortices and zonal flows, and their impact on the Lagrangian statistics of particles. Computationally intensive long time simulations following ensembles of test particles over hundreds of eddy turnover times are considered in statistically stationary turbulent flows. The flow topology is characterized using the Lagrangian Okubo–Weiss criterion in order to split the flow into topologically different domains. In elliptic and hyperbolic regions, the probability density functions (PDFs) of the residence time have self-similar algebraic decaying tails. However, in the intermediate regions, the PDFs exhibit exponentially decaying tails. Topologically conditioned PDFs of the Lagrangian velocity, and acceleration and density fluctuations are likewise computed. The differences between the classical Hasegawa–Wakatani system and its modified version are assessed, and the role of zonal flows is highlighted. Here, the density flux spectrum, which characterizes the contributions of different length scales, is studied, and its inertial scaling is found to be in agreement with predictions based on dimensional arguments. Analyzing the angular change of particle tracers at different time scales, corresponding to coarse grained curvature, completes the study, and these multiscale geometric statistics quantify the directional properties of the particle motion in different flow regimes.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Settling of two-way momentum and energy coupled particles subject to Boussinesq and non-Boussinesq heating
This work establishes a procedure to accurately compute heat transfer between an Eulerian fluid and Lagrangian point-particles. Recent work has focused on accurately computing momentum transfer between fluid and particles. The coupling term for momentum involves the undisturbed fluid velocity at the particle location which is not directly accessible in the simulation. Analogously, in the context of thermal coupling, the undisturbed fluid temperature at the particle location is not directly accessible in simulations and must be estimated. In this paper, we develop a scheme to accurately estimate the undisturbed fluid temperature of a point-particle exchanging thermal energy with a surrounding fluid. Furthermore, the temperature disturbance is correlated with the enhanced temperature curvature in the vicinity of the particle and is formally valid in the low heating, low convection limit. We conduct extensive verification of the correction procedure for a settling particle subject to radiation. This setup allows the simultaneous testing of thermal and momentum corrections. By considering equations of drag and Nusselt number extended to finite Péclet and Boussinesq numbers, we establish a large range over which the correction procedure can be applied.
42 ENGINEERING↗
Numerical Analysis of Regular Material Point Method and its Application to Multiphase Flows
The material point method (MPM) is gaining wide popularity in engineering research to model and simulate complex multiphase flow dynamics. The method relies on solving the governing equations of motion and transport in a Lagrangian framework using particles also known as material points. The fluid and kinematic properties are stored on the material points while the spatial gradient calculation and temporal integration are performed on a background grid. This Lagrangian framework allows for large deformations, easy integration of constitutive models, and direct import of complex geometries as particles. However, despite their increasing popularity, very few studies have addressed the issues of numerical resolution and stability of MPM techniques. The presence of additional factors such as the number of material points-per-cell, the location of the material points, the CFL-like condition used in time update, and the grid shape functions also increase the complexity of the error analysis when compared to other finite element methods. In this presentation, we analyze the various forms of error incurred in the application of MPM to continuum mechanics and multiphase flows. The effect of the previously mentioned factors on the error dynamics is studied. The application of these principles to canonical and industrial problems is also presented.
high pressure reverse osmosis↗
Implementing a Stochastic Soot Model in a 3D Domain
In this document a highly detailed stochastic soot model limited to 0D simulations is coupled to 3D domains through the use of a Lagrangian parcel-Eulerian fluid method. Although both the stochastic soot model used and the Lagrangian parcel-Eulerian fluid method have long histories of development, their use, as described in this document, is novel and useful. Previous works significant to the model's development include documents from both the stochastic and simpler soot model genre and the fuel, aerosol, and other particle tracking and Lagrangian parcel-Eulerian fluid genres. Developing the theory behind the model is crucial to its implementation because of the unconventional Lagrangian parcel-Eulerian fluid representation. In fact, the Lagrangian parcel-Eulerian fluid coupling is only a means of simplifying the much more difficult theoretical model. Despite the accompanying assumptions, the complete model is validated from the bottom up and then again in a holistic fashion with promising results. Then the model is used in multiple case studies, showcasing capabilities not had by other models. The optical properties of soot generated in a Spray A simulation are calculated to find measured results are likely overpredicting the mass of soot along the central axis via KL extinction measurements. A thermophoresis submodel that calculates force based on the soot particles' shape is installed to find it has a high influence on the path of the heaviest soot particles. Investigations into the HACA mechanism for soot growth and the effect of syngas on soot in 3D domains are performed with much more detail than has been available before.
Strickland, Tyler↗
Learning effective physical laws for generating cosmological hydrodynamics with Lagrangian deep learning
The goal of generative models is to learn the intricate relations between the data to create new simulated data, but current approaches fail in very high dimensions. When the true data-generating process is based on physical processes, these impose symmetries and constraints, and the generative model can be created by learning an effective description of the underlying physics, which enables scaling of the generative model to very high dimensions. In this work, we propose Lagrangian deep learning (LDL) for this purpose, applying it to learn outputs of cosmological hydrodynamical simulations. The model uses layers of Lagrangian displacements of particles describing the observables to learn the effective physical laws. The displacements are modeled as the gradient of an effective potential, which explicitly satisfies the translational and rotational invariance. The total number of learned parameters is only of order 10, and they can be viewed as effective theory parameters. Further, we combine N-body solver fast particle mesh (FastPM) with LDL and apply it to a wide range of cosmological outputs, from the dark matter to the stellar maps, gas density, and temperature. The computational cost of LDL is nearly four orders of magnitude lower than that of the full hydrodynamical simulations, yet it outperforms them at the same resolution. We achieve this with only of order 10 layers from the initial conditions to the final output, in contrast to typical cosmological simulations with thousands of time steps. This opens up the possibility of analyzing cosmological observations entirely within this framework, without the need for large dark-matter simulations.
79 ASTRONOMY AND ASTROPHYSICS↗
Equation of State Measurements on Iron Near the Melting Curve at Planetary Core Conditions by Shock and Ramp Compressions
The outer core of the Earth is composed primarily of liquid iron, and the inner core boundary is governed by the intersection of the melt line and the geotherm. While there are many studies on the thermodynamic equation of state for solid iron, the equation of state of liquid iron is relatively unexplored. In this work, we use dynamic compression to diagnose the high-pressure liquid equation of state of iron by utilizing the shock-ramp capability at Sandia National Laboratories’ Z-Machine. This technique enables measurements of material states off the Hugoniot by initially shocking samples and subsequently driving a further, shockless compression. Planetary studies benefit greatly from isentropic, off-Hugoniot experiments since they can cover pressure-temperature (P-T) conditions that are close to adiabatic profiles found in planetary interiors. We used this method to drive iron to P-T conditions similar to those of the Earth’s outer-inner core boundary, along an elevated-temperature isentrope in the liquid from 275 GPa to 400 GPa. We derive the equation of state using a hybrid backward integration – forward Lagrangian technique on particle velocity traces to determine the pressure-density history of the sample. Our results are in excellent agreement with SESAME 92141, a previously published equation of state table. With our data and previous experimental data on liquid iron we provide new information on the iron melting line and derive new parameters for a Vinet-based equation of state. The table and our parameterized equation of state are applied to provide an updated means of modeling the pressure, mass, and density of liquid iron cores in exoplanetary interiors.
58 GEOSCIENCES↗