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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Testing the conservative character of particle simulations: I. Canonical and noncanonical guiding center model in Boozer coordinates

The guiding center (GC) Lagrangian in Boozer coordinates for toroidally confined plasmas can be cast into canonical form by eliminating terms containing the covariant component B Ψ P of the magnetic field vector with respect to the poloidal flux function Ψ P . In an unperturbed plasma, B Ψ P can be eliminated via exact coordinate transformations, but, in general, one relies on approximations, assuming that the effect of B Ψ P is small. In this work, we are interested in the question whether Hamiltonian conservation laws are still satisfied when B Ψ P is retained in the presence of fluctuations. Considering fast ions in the presence of a shear Alfvén wave field with fixed amplitude, fixed frequency, and a single toroidal mode number n, we show that simulations using the code ORBIT with and without B Ψ P yield practically the same resonant and nonresonant GC orbits. The numerical results are consistent with theoretical analyses (presented in the appendix), which show that the unabridged GC Lagrangian with B Ψ P retained yields equations of motion that possess two key properties of Hamiltonian flows: (i) phase space conservation and (ii) energy conservation. As counter-examples, we also show cases where energy conservation (ii) or both conservation laws (i) and (ii) are broken by omitting certain small terms. When testing the conservative character of the simulation code, it is found to be beneficial to apply perturbations that do not resemble normal (eigen) modes of the plasma. The deviations are enhanced and, thus, more easily spotted when one inspects wave-particle interactions using nonnormal modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A magnetic analog of pressure–strain interaction

We study the evolution equation for magnetic energy density for a non-relativistic magnetized plasma in the (Lagrangian) reference frame comoving with the electron bulk velocity. Analyzing the terms that arise due to the ideal electric field, namely, perpendicular electron compression and magnetic field line bending, we recast them to reveal a quantity with a functional form analogous to the often-studied pressure–strain interaction term that describes one piece of internal energy density evolution of the species in a plasma, except with the species pressure tensor replaced by the magnetic stress tensor. We dub it the “magnetic stress–strain interaction.” We discuss decompositions of the magnetic stress–strain interaction analogous to those used for pressure–strain interaction. These analogies facilitate the interpretation of the evolution of the various forms of energy in magnetized plasmas and should be useful for a wide array of applications, including magnetic reconnection, turbulence, collisionless shocks, and wave–particle interactions. We display and analyze all the terms that can change magnetic energy density in the Lagrangian reference frame of the electrons using a particle-in-cell simulation of magnetic reconnection.

Barbhuiya, M. Hasan (ORCID:0000000163301650)↗
DOE OSTI2020

The MFiX Particle-in-Cell Method (MFiX-PIC) Theory Guide

MFiX (Multiphase Flow with Interphase eXchanges) is an open-source multiphase flow solver developed at the National Energy Technology Laboratory. Within the code, users have access to a single phase or interpenetrating continua-based multiphase two-fluid model (TFM), a discrete element model (DEM), and a particle-in-cell model (PIC). TFM, DEM, and PIC can all be used to create multiphase simulations that include hydrodynamics, chemical reactions, and heat transfer.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Slow manifolds of classical Pauli particle enable structure-preserving geometric algorithms for guiding center dynamics

Since variational symplectic integrators for the guiding center was proposed, structure-preserving geometric algorithms have become an active research field in plasma physics. We found that the slow manifolds of the classical Pauli particle enable a family of structure-preserving geometric algorithms for guiding center dynamics with long-term stability and accuracy. This discovery overcomes the difficulty associated with the unstable parasitic modes for variational symplectic integrators when applied to the degenerate guiding center Lagrangian. It is a pleasant surprise that Pauli's Hamiltonian for electrons, which predated the Dirac equation and marks the beginning of particle physics, reappears in classical physics as an effective algorithm for solving an important plasma physics problem. This technique is applicable to other degenerate Lagrangians reduced from regular Lagrangians.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Dataset for "FLEXPART-WRF simulations to assess dust transport in the convective boundary layer at the Southern Great Plains site"

Soil dusts are an important source of aerosol in agricultural regions and can affect the Earth’s radiation budget through the modification of cloud properties, and in particular, through their ability to act as ice nucleating particles. In order to impact cloud properties, agricultural soil dusts need to be transported from the point of emission to cloud-relevant altitudes. Vertical transport within the planetary boundary layer is strongly controlled by turbulence and is challenging to represent accurately in regional and global models. Large-eddy simulations (LES) are run at resolutions capable of resolving most of the turbulent energy directly and can thus better simulate vertical transport. In this study, we leverage the LES ARM Symbiotic Simulation and Observation (LASSO) large-eddy simulations to simulate vertical transport of agricultural dust within the turbulent boundary layer using a modified version of the stochastic dispersion model FLEXPART-WRF. We find that the modified model is better capable of simulating particle transport due to turbulence, and that particle size was the greatest factor in determining particle lifetime. Individual meteorology and particle density had intermediate effects upon particle transport, while release height had little effect upon simulation results. Finally, we utilize a quasi-single column model (QSCM) approach to determine how our results compare to a parameterized treatment of turbulence. The QSCM simulations led to greater tracer transport out of the boundary layer, with ramifications for any studies utilizing a Lagrangian stochastic model to understand tracer dispersion. These results highlight the importance of accurately simulating turbulence for understanding particle transport.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A High Resolution Simulation of a Single Shock-Accelerated Particle

We report particle drag models, which capture macroviscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a two-dimensional axisymmetric Navier–Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted using the hydrocode FLAG. FLAGs capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated—a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAGs Navier–Stokes and heat diffusion solutions and to provide a rationale for recent experimental particle drag measurements.

42 ENGINEERING↗
DOE OSTI2020

A Numerical Simulation of a Single Shock-Accelerated Particle

Particle drag models, which capture macro viscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a 2D axisymmetric Navier-Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted in the hydrocode FLAG. FLAG’s capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated – a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAG’s Navier-Stokes and heat diffusion solutions, and to provide a rationale for recent experimental particle drag measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds

Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Super‐Droplet Method to Simulate Lagrangian Microphysics of Nuclear Fallout in a Homogeneous Cloud

Abstract Nuclear detonations produce hazardous local and global particles or fallout. Predicting fallout size, chemical components, and location is necessary to inform officials and determine immediate guidance for the public. However, existing nuclear detonation fallout models prescribe the particle size distributions based on limited observations. In this work, we apply the super‐droplet method, which is a numerical modeling technique developed for cloud microphysics, to simulate size distributions of particles in a mushroom cloud formed post‐detonation of a nuclear device. We model fallout formation and evolution with homogeneous nucleation and condensation of a single species and a Monte Carlo coagulation algorithm. We verify the numerical methods representing coagulation and condensation processes against analytical test problems. Additionally, we explore several scenarios for the integral system mass and yield in equivalent kilotons (kt) of TNT (trinitrotoluene). The fallout size distribution median diameter d pg follows a scaling law based on the integral system mass m v 0 kg and yield Y kt: nm. We test the effect of cloud turbulence, enhanced nucleation and growth, and vapor volatility with a sensitivity study. The range in median diameter predictions for simulations of historical tests performed over the Pacific encompass the measurements of particles sampled from the cloud caps. Predicted median particle size ranges up to 217, 123, 86, and 35 nm for historical tests with yields of 0.2, 0.7, 2, and 10 Mt, respectively. This work can be expanded in many different directions to build a more predictive model for fallout formation.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Light quark loops in ${K}^{\pm}\to {\pi}^{\pm}\nu \overline{\nu}$ from vector meson dominance and update on the Kaon Unitarity Triangle

We use vector meson dominance to calculate non-perturbative contributions to the branching ratio of the rare decay $K$ ± → π ± $v\overline{v}$ stemming from matrix elements involving up-quark loops. The importance of this observable as well as of K 0 → π 0 l + l - and of the direct CP violation parameter ϵ' K is then discussed in the context of a Unitarity Triangle sqtudy based on Kaon sector observables only.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Chiral algebra, localization, modularity, surface defects, and all that

We study the 2D vertex operator algebra (VOA) construction in 4D N = 2 superconformal field theories on S 3 × S 1 , focusing on both old puzzles and new observations. The VOA lives on a two-torus T 2 ⊂ S 3 × S 1 , it is 1 2 Z -graded, and this torus is equipped with the natural choice of spin structure (1,0) for the Z + 1 2 -graded operators, corresponding to the NS sector vacuum character. By analyzing the possible refinements of the Schur index that preserves the VOA, we find that it admits discrete deformations, which allows access to the remaining spin structures (1,1), (0,1), and (0,0), of which the latter two involve the inclusion of a particular surface defect. For Lagrangian theories, we perform the detailed analysis: we describe the natural supersymmetric background, perform localization, and derive the gauged symplectic boson action on a torus in any spin structure. In the absence of flavor fugacities, the 2D and 4D path integrals precisely match, including the Casimir factors. We further analyze the 2D theory: we identify its integration cycle and the two-point functions and interpret flavor holonomies as screening charges in the VOA. Next, we make some observations about modularity; the T-transformation acts on our four partition functions and lifts to a large diffeomorphism on S 3 × S 1 . More interestingly, we generalize the four partition functions on the torus to an infinite family labeled by both the spin structure and the integration cycle inside the complexified maximal torus of the gauge group. Members of this family transform into one another under the full modular group, and we confirm the recent observation that the S-transform of the Schur index in Lagrangian theories exhibits logarithmic behavior. Finally, we comment on how locally our background reproduces the Ω-background.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Simulated Dust Transport in the Convective Boundary Layer

Soil dusts are an important source of aerosol in agricultural regions and can affect the Earth's radiation budget through the modification of cloud properties, and in particular, through their ability to act as ice nucleating particles. In order to impact cloud properties, agricultural soil dusts need to be transported from the point of emission to cloud–relevant altitudes. Vertical transport within the planetary boundary layer is strongly controlled by turbulence and is challenging to represent accurately in regional and global models. Large–eddy simulations (LES) are run at resolutions capable of resolving most of the turbulent energy directly and can thus better simulate vertical transport. In this study, we leverage the LES ARM Symbiotic Simulation and Observation (LASSO) large–eddy simulations to simulate vertical transport of agricultural dust within the turbulent boundary layer using a modified version of the stochastic dispersion model FLEXPART–WRF. We find that the modified model is better capable of simulating particle transport due to turbulence, and that particle size was the greatest factor in determining particle lifetime. Individual meteorology and particle density had intermediate effects upon particle transport, while release height had little effect upon simulation results. Finally, we utilize a quasi–single column model (QSCM) approach to determine how our results compare to a parameterized treatment of turbulence. The QSCM simulations led to greater tracer transport out of the boundary layer, with ramifications for any studies utilizing a Lagrangian stochastic model to understand tracer dispersion. Furthermore, these results highlight the importance of accurately simulating turbulence for understanding particle transport.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Baker-Coon-Romans N -point amplitude and an exact field theory limit of the Coon amplitude

We study the N-point Coon amplitude discovered first by Baker and Coon in the 1970s and then again independently by Romans in the 1980s. This Baker-Coon-Romans (BCR) amplitude retains several properties of tree-level string amplitudes, namely duality and factorization, with a q-deformed version of the string spectrum. Although the formula for the N-point BCR amplitude is only valid for q > 1, the four-point case admits a straightforward extension to all q ≥ 0 which reproduces the usual expression for the four-point Coon amplitude. At five points, there are inconsistencies with factorization when pushing q < 1. Despite these issues, we find a new relation between the five-point BCR amplitude and Cheung and Remmen’s four-point basic hypergeometric amplitude, placing the latter within the broader family of Coon amplitudes. Finally, we compute the q → ∞ limit of the N-point BCR amplitudes and discover an exact correspondence between these amplitudes and the field theory amplitudes of a scalar transforming in the adjoint representation of a global symmetry group with an infinite set of non-derivative single-trace interaction terms. This correspondence at q = ∞ is the first definitive realization of the Coon amplitude (in any limit) from a field theory described by an explicit Lagrangian.

1/N Expansion↗
DOE OSTI2020

Heating neutron stars with GeV dark matter

An old neutron star (NS) may capture halo dark matter (DM) and get heated up by the deposited kinetic energy, thus behaving like a thermal DM detector with sensitivity to a wide range of DM masses and a variety of DM-quark interactions. Near future infrared telescopes will measure NS temperatures down to a few thousand Kelvin and probe NS heating by DM capture. We focus on GeV-mass Dirac fermion DM (which is beyond the reach of current DM direct detection experiments) in scenarios in which the DM capture rate can saturate the geometric limit. For concreteness, we study (1) a model that invokes dark decays of the neutron to explain the neutron lifetime anomaly, and (2) a framework of DM coupled to quarks through a vector current portal. In the neutron dark decay model a NS can have a substantial DM population, so that the DM capture rate can reach the geometric limit through DM self-interactions even if the DM-neutron scattering cross section is tiny. We find NS heating to have greater sensitivity than multipion signatures in large underground detectors for the neutron dark decay model, and sub-GeV gamma-ray signatures for the quark vector portal model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Three exceptions to the Grossman-Nir bound

A bstract We show that the Grossman-Nir (GN) bound, Br( K L → $$ {\pi}^0\nu \overline{\nu} $$ π 0 ν ν ¯ ) ≤ 4 . 3 Br( K + → $$ {\pi}^{+}\nu \overline{\nu} $$ π + ν ν ¯ ), can be violated in the presence of light new physics with flavor violating couplings. We construct three sample models in which the GN bound can be violated by orders of magnitude, while satisfying all other experimental bounds. In the three models the enhanced branching ratio Br( K L → π 0 + inv) is due to K L → π 0 Φ 1, K L → π 0 Φ 1 Φ 1 , K L → $$ {\pi}^0{\psi}_1{\overline{\psi}}_1 $$ π 0 ψ 1 ψ ¯ 1 transitions, respectively, where Φ 1 ( ψ 1 ) is a light scalar (fermion) that escapes the detector. In the three models Br( K + → π + + inv) remains very close to the SM value, while Br( K L → π 0 + inv) can saturate the present KOTO bound. Besides invisible particles in the final state (which may account for dark matter) the models require additional light mediators around the GeV-scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Flash-X: A multiphysics simulation software instrument

Flash-X is a highly composable multiphysics software system that can be used to simulate physical phenomena in several scientific domains. It derives some of its solvers from FLASH, which was first released in 2000. Flash-X has a new framework that relies on abstractions and asynchronous communications for performance portability across a range of increasingly heterogeneous hardware platforms. Flash-X is meant primarily for solving Eulerian formulations of applications with compressible and/or incompressible reactive flows. It also has a built-in, versatile Lagrangian framework that can be used in many different ways, including implementing tracers, particle-in-cell simulations, and immersed boundary methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Lagrangian Characterization of Surface Transport From the Equatorial Atlantic to the Caribbean Sea Using Climatological Lagrangian Coherent Structures and Self‐Organizing Maps

Abstract This study presents an assessment of the transport of suspended material by surface ocean currents, which have a critical role in determining the connectivity and distribution of living and non‐living material. Lagrangian experiments reveal pathways from the Equatorial Atlantic to 10 strategic regions within the Caribbean Sea, determined by considering the space‐time variability of climatological Lagrangian Coherent Structures, which act as recurrent attracting pathways and transport barriers. Due to windage or Stokes drift, wind forcing is a significant factor in determining the spatial locations where particles cluster and the time needed to reach the Caribbean from the Equatorial Atlantic. Pathways shift westward within the Caribbean and take less time to arrive with increasing wind influence. Depending on the wind effect, the particles show higher confluence in different areas of the Caribbean. A case study is presented for the Mexican Caribbean nearshore area, isolated from ocean‐current trajectories. Here, wind weakens the transport barrier responsible for this isolation and causes particle confluence toward that region. Spatial patterns of the Eulerian velocity identified through Self‐Organizing Maps, with time dependence given by their best matching units, can reproduce the characteristic Lagrangian patterns of surface current climate variability. Our study demonstrates the application of tools from dynamical systems and unsupervised neural networks to understand Lagrangian patterns and identify the processes that drive them. These findings improve our understanding of transport mechanisms of suspended material by surface ocean currents in the Western Atlantic and the Caribbean Sea, which is essential for managing and conserving marine ecosystems.

Allende‐Arandía, Ma. Eugenia↗