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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Amplitude/operator basis in chiral perturbation theory

We establish a systematic construction of the on-shell amplitude/operator basis for Chiral Perturbation Theory (ChPT) in D = 4 spacetime dimensions and with an arbitrary number of flavors N f . For kinematic factors, we employ spinor-helicity variables to construct the soft blocks, which are local amplitudes satisfying the Adler’s zero condition, as well as to take into account the reduction in the kinematic basis due to the Gram determinant, which arises at O(p 10 ) when the number of multiplicity N in an amplitude becomes large: N > D. For flavor factors, we include group-theoretic relations at small N f , N f ≤ N, which decreases the flavor basis. The result is obtained by adapting the Young tensor method of constructing the operator basis for generic effective field theories to the case of non-linearly realized symmetries. Working in the massless quark limit, we present purely mesonic operators for both even- and odd-parity at O(p 6 ) and O(p 8 ) for N = 6 and arbitrary N f , and establish a direct correspondence between the amplitude basis and the operator basis. Furthermore, the redundancy due to the Gram determinant is studied at O(p 10 ) for N = 6, 8 and 10.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Radiative semileptonic $ \overline{B} $ decays

We consider the form factors for the radiative semileptonic decays $ \overline{B} $(v) → D(*)(v')ℓν¯ ℓγ in the kinematic region where the photon momentum, k, is small enough that heavy quark symmetry (HQS) can be applied without the radiated photon changing the heavy quark velocity (i.e., v(') ∙ k < m (b,c) ). We find that HQS is remarkably powerful, leaving only four new undetermined form factors at leading order in 1/m (b,c) . In addition, one of them is fixed in terms of the leading order Isgur-Wise function in the kinematic region, v(') ∙ k < Λ QCD .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

The three-pion K -matrix at NLO in ChPT

The three-particle K-matrix, Κ df,3 , is a scheme-dependent quantity that parametrizes short-range three-particle interactions in the relativistic-field-theory three-particle finite-volume formalism. In this work, we compute its value for systems of three pions in all isospin channels through next-to-leading order in Chiral Perturbation Theory, generalizing previous work done at maximum isospin. We obtain analytic expressions through quadratic order (or cubic order, in the case of zero isospin) in the expansion about the three-pion threshold.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

The isospin-3 three-particle K -matrix at NLO in ChPT

The three-particle K-matrix, K df,3 , is a scheme-dependent quantity that parametrizes short-range three-particle interactions in the relativistic-field-theory three-particle finite-volume formalism. In this work, we compute its value for systems of three pions at maximal isospin through next-to-leading order (NLO) in Chiral Perturbation Theory (ChPT). We compare the values to existing lattice QCD results and find that the agreement between lattice QCD data and ChPT in the first two coefficients of the threshold expansion of K df,3 is significantly improved with respect to leading order once NLO effects are incorporated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A lattice study of ππ scattering at large N c

We present the first lattice study of pion-pion scattering with varying number of colors, N c . We use lattice simulations with four degenerate quark flavors, N f = 4, and N c = 3 – 6. We focus on two scattering channels that do not involve vacuum diagrams. These correspond to two irreducible representations of the SU(4) flavor group: the fully symmetric one, SS, and the fully antisymmetric one, AA. The former is a repulsive channel equivalent to the isospin-2 channel of SU(2). By contrast, the latter is attractive and only exists for N f ≥ 4. A representative state is (|D$^{+}_{s}$π + > – | D + K + >)/√2. Using Lüscher’s formalism, we extract the near-threshold scattering amplitude and we match our results to Chiral Perturbation Theory (ChPT) at large N c . For this, we compute the analytical U(N f ) ChPT prediction for two-pion scattering, and use the lattice results to constrain the N c scaling of the relevant low-energy couplings.

1/N Expansion↗
DOE OSTI2023

Isospin-breaking effects in the three-pion contribution to hadronic vacuum polarization

Isospin-breaking (IB) effects are required for an evaluation of hadronic vacuum polarization at subpercent precision. While the dominant contributions arise from the e + e - → π + π - channel, also IB in the subleading channels can become relevant for a detailed understanding, e.g., of the comparison to lattice QCD. Here, we provide such an analysis for e + e - → 3π by extending our dispersive description of the process, including estimates of final-state radiation (FSR) and ρ–ω mixing. In particular, we develop a formalism to capture the leading infrared-enhanced effects in terms of a correction factor η3π that generalizes the analog treatment of virtual and final-state photons in the 2π case. The global fit to the e+e- → 3π data base, subject to constraints from analyticity, unitarity, and the chiral anomaly, gives for the total 3π contribution to the anomalous magnetic moment of the muon, of which ${a}^{FSR}_\mu [3\pi]$ = 0.51(1) x 10 -10 and ${a}^{\rho - \omega}_\mu [3\pi]$ = -2.68(70) x 10 -10 can be ascribed to IB. We argue that the resulting cancellation with ρ–ω mixing in e + e - → 2π can be understood from a narrow-resonance picture, and provide updated values for the vacuum-polarization-subtracted vector-meson parameters M ω = 782.70(3) MeV, M Φ = 1019.21(2) MeV, Γ ω = 8.71(3) MeV, and Γ Φ = 4.27(1) MeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Accurate hydrodynamic force and torque correlations for prolate spheroids from Stokes regime to high Reynolds numbers

Detailed simulations of flow around various prolate spheroids are performed using multi-relaxation time lattice Boltzmann method (LBM). The simulations are performed in the Reynolds numbers range 0.1 ≤ Re ≤ 2000 at different incident angles 0°≤ φ ≤ 90° for various prolate spheroids of aspect ratios 1 ≤ λ ≤ 8. The Re is based on the volume equivalent sphere diameter. From the simulations, accurate correlations for average drag, lift and torque coefficients (C D , C L , and C T respectively) are proposed. The mean deviations between the correlations and the simulation data for C D , C L , and C T are 2%, 6.5%, and 4.7% respectively. Furthermore, the correlations are fitted with the physics in mind such that they can be extrapolated beyond the regimes simulated. The fits are designed such that correlations reduce to analytical (C D and C L ) solutions available for prolate spheroids in the limit Re ≈ 0. At high Re, the correlations mimic the asymptotic flattening of C D as for spheres and therefore, we expect the correlations to be valid until the critical Re (≈ 10 5 ). Furthermore, extensive comparison of our correlations with other high-fidelity simulations from literature demonstrate the accuracy of our work compared to other correlations from literature. This work enables researchers to perform accurate unresolved Euler-Lagrangian simulations for a wide range of elongated particles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1985

Decay of a coherent scalar disturbance in a turbulent flow

The time evolution of an initially coherent, sinusoidal passive-scalar disturbance is considered when the wavelength q is less than the length scale of the surrounding isotropic turbulent flow. In 64 sup 3 direct numerical simulations a Gaussian prescription for the average scalar amplitude breaks down after a timescale associated with the wavenumber of the disturbance and there is a transition to a new characteristic decay. The Gaussian prescription is given by exp(-(1/2) q-squared w(t)), where a form for w(t), the Lagrangian mean square displacement of a single fluid particle, is proposed. After the transition the decay is given by exp(-t/tau), where tau is the new characteristic timescale. If q k(sub K), then 1/tau = 1/tau(sub D) + 1/tau(sub K), where k(sub K) is the Kolmogorov wavenumber, tau(sub D) is the diffusive timescale and tau(sub K) is the Kolmogorov timescale. An experiment originally proposed by de Gennesis considered in which the evolution of a coherent laser-induced pattern is read by a diffracting laser. The theory of this experiment involves the dispersion of particle pairs, but it is shown that in a certain limit it reduces to the single Fourier-mode problem and can be described in terms of single particle diffusion. The decay of a single mode after the transition in the simulation best describes the experiment.

Kerr, R. M.↗
NASA NTRS1991

Planetary engineering

Assuming commercial fusion power, heavy lift vehicles and major advances in genetic engineering, the authors survey possible late-21st century methods of working major transformations in planetary environments. Much more Earthlike climates may be produced on Mars by generating low freezing point greenhouse gases from indigenous materials; on Venus by biological conversion of CO2 to graphite, by canceling the greenhouse effect with high-altitude absorbing fine particles, or by a sunshield at the first Lagrangian point; and on Titan by greenhouses and/or fusion warming. However, in our present state of ignorance we cannot guarantee a stable endstate or exclude unanticipated climatic feedbacks or other unintended consequences. Moreover, as the authors illustrate by several examples, many conceivable modes of planetary engineering are so wasteful of scarce solar system resources and so destructive of important scientific information as to raise profound ethical issues, even if they were economically feasible, which they are not. Global warming on Earth may lead to calls for mitigation by planetary engineering, e.g., emplacement and replenishment of anti-greenhouse layers at high altitudes, or sunshields in space. But here especially we must be concerned about precision, stability, and inadvertent side-effects. The safest and most cost-effective means of countering global warming - beyond, e.g., improved energy efficiency, CFC bans and alternative energy sources - is the continuing reforestation of approximately 2 times 107 sq km of the Earth's surface. This can be accomplished with present technology and probably at the least cost.

Pollack, James B.↗
NASA NTRS1992

Numerical modeling for dilute and dense sprays

We have successfully implemented a numerical model for spray-combustion calculations. In this model, the governing gas-phase equations in Eulerian coordinate are solved by a time-marching multiple pressure correction procedure based on the operator-splitting technique. The droplet-phase equations in Lagrangian coordinate are solved by a stochastic discrete particle technique. In order to simplify the calculation procedure for the circulating droplets, the effective conductivity model is utilized. The k-epsilon models are utilized to characterize the time and length scales of the gas phase in conjunction with turbulent modulation by droplets and droplet dispersion by turbulence. This method entails random sampling of instantaneous gas flow properties and the stochastic process requires a large number of computational parcels to produce the satisfactory dispersion distributions even for rather dilute sprays. Two major improvements in spray combustion modelings were made. Firstly, we have developed a probability density function approach in multidimensional space to represent a specific computational particle. Secondly, we incorporate the Taylor Analogy Breakup (TAB) model for handling the dense spray effects. This breakup model is based on the reasonable assumption that atomization and drop breakup are indistinguishable processes within a dense spray near the nozzle exit. Accordingly, atomization is prescribed by injecting drops which have a characteristic size equal to the nozzle exit diameter. Example problems include the nearly homogeneous and inhomogeneous turbulent particle dispersion, and the non-evaporating, evaporating, and burning dense sprays. Comparison with experimental data will be discussed in detail.

Chen, C. P.↗
NASA NTRS1990

On the use of Lagrangian variables in descriptions of unsteady boundary-layer separation

The Lagrangian description of unsteady boundary layer separation is reviewed from both analytical and numerical perspectives. It is explained in simple terms how particle distortion gives rise to unsteady separation, and why a theory centered on Lagrangian coordinates provides the clearest description of this phenomenon. Some of the more recent results for unsteady three dimensional compressible separation are included. The different forms of separation that can arise from symmetries are emphasized. A possible description of separation is also included when the detaching vorticity layer exits the classical boundary layer region, but still remains much closer to the surface than a typical body-lengthscale.

Cowley, Stephen J.↗
DOE OSTI2020

An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing

The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Hybrid Particle Method for Computational Shock Physics

A long-standing area of research for Eulerian shock wave physics codes has been the treatment of strength and damage for materials. Here we present a method that will aid in the analysis of strength and failure in shock physics applications where excessive diffusion of critical variables can occur and control the solution outcome. Eulerian methods excel for large deformation simulations in general but are inaccurate in capturing structural behavior. Lagrangian methods provide better structural response, but finite element meshes can become tangled. Therefore, a technique for merging Lagrangian and Eulerian treatments of material response, within a single numerical framework, was implemented in the Multiple Component computational shock physics hydrocode. The capability is a Lagrangian/Eulerian Particle Method (LEPM) that uses particles to interface a Lagrangian treatment of material strength with a more traditional Eulerian treatment of the Equation of State (EOS). Lagrangian numerical methods avoid the advection diffusion found in Eulerian methods, which typically strongly affects strength constitutive law internal variables, such as equivalent plastic strain, porosity and/or damage. The Lagrangian capability enhances existing capabilities and permits accurate predictions of high rate, large deformation and/or shock of mechanical structures.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2004

Investigation of Polar Stratospheric Cloud Solid Particle Formation Mechanisms Using ILAS and AVHRR Observations in the Arctic

Satellite observations of denitrification and ice clouds in the Arctic lower stratosphere in February 1997 are used with Lagrangian microphysical box model calculations to evaluate nucleation mechanisms of solid polar stratospheric cloud (PSC) particles. The occurrences of ice clouds are not correlated in time and space with the locations of back trajectories of denitrified air masses, indicating that ice particle surfaces are not always a prerequisite for the formation of solid PSCs that lead to denitrification. In contrast, the model calculations incorporating a pseudoheterogeneous freezing process occurring at the vapor-liquid interface can quantitatively explain most of the observed denitrification when the nucleation activation free energy for nitric acid dihydrate formation is raised by only approx.10% relative to the current published values. Once nucleated, the conversion of nitric acid dihydrate to the stable trihydrate phase brings the computed levels of denitrification closer to the measurements. INDEX TERMS: 0305 Atmospheric Composition and Structure: Aerosols and particles (0345, 4801); 0320 Atmospheric Composition and Sb~lctureC: loud physics and chemistry; 0340 Atmospheric Composition and Structure: Middle atmosphere-composition and chemistry

Irie, H.↗
NASA NTRS1993

Supersonic gasdispersional jets and jet noise

This paper examines the potential for controlling jet noise radiation using methods developed for modifying jet infra-red thermal radiation. The control of jet noise may be possible by properly adding different solid and liquid particles into the jet flow and by using special nozzle shapes to change the jet exhaust flow structure. The numerical methods used to achieve these objectives are outlined in this paper. A combined Lagrangian-Eulerian approach is used to numerically simulate a Jet flow with particle addition. The unsteady behavior of jet impingement is examined. The techniques for grid and boundary condition definition are discussed as related to the accuracy of the calculations. Preliminary comparisons to experimental data are presented.

Gilinsky, M. M.↗
NASA NTRS1997

Solar energetic particle events and coronal mass ejections: New insights from SOHO

The scientific payload of SOHO, launched in December 1995, enables comprehensive studies of the sun from its interior, to the outer corona and solar wind. In its halo orbit around the Lagrangian point of the sun-earth system, the comprehensive suprathermal and energetic particle analyzer (COSTEP) measures in situ energetic partiles in the energy range 44 keV/particle to greater than 53 MeV/n. Although solar activity was at minimum, COSTEP detected from mid December 1995 until the end of July 1997, 30 solar energetic particle (SEP) events, including both gradual and implusive type SEPs. These minimum phase SEP events are unique in the sense that their associated solar source phenomena can be investigated in detail without interference by other simultaneous solar events as is usually the case at times around solar activity maximum. Simultaneous observations of the solar corona are provided by the large angle spectroscopic coronagraph (LASCO) and the extreme ultraviolet imaging telescope (EIT). From the correlated SOHO observations, a one to one correspondence of SEP events with coronal mass ejections (CMEs) was found. Most of the SEP events were associated with west-limb CMEs, some with halo CMEs that later passed the SOHO spacecraft and with Moreton-like disturbances in the lower solar atmosphere as observed by the EIT. Many SEP events were detected at sector boundaries of the interplanetary magnetic field (IMF) suggesting a magnetic connection to coronal streamers at the sun as supported by LASCO observations of mass ejections at the base of helmet streamers. Energetic particle and LASCO white-light observations yield evidence that CMEs often lead to large-scale disturbances of the sun's corona, probably affecting at times areas all around the sun.

Bothmer, V.↗
NASA NTRSDate not supplied

Demonstration of Initial GlennICE Relative Frame Capability: Axial-Flow Propeller

Modifications have been implemented in the GlennICE software to accommodate a non-inertial reference frame. GlennICE accepts a flow solution from an external flow solver. It then introduces particles and tracks them through the flow field in a Lagrangian manner. Centrifugal and Coriolis terms were added to the GlennICE software to account for relative frame simulations. The objective of the present paper is twofold. First, to check that the new terms are implemented correctly and that the code still behaves as expected with respect to convergence. And second, to provide some initial insight into an upcoming propeller experiment in the NASA Icing Research Tunnel. The paper presents a description of the code modifications. In addition, results are presented for two operating conditions, and three particle sizes. Each case was simulated with four different grid densities to assess grid dependence.

Aircraft Icing↗