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 289 records · Page 16

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↗

Simulations and symmetries

ABSTRACT We investigate the range of applicability of a model for the real-space power spectrum based on N-body dynamics and a (quadratic) Lagrangian bias expansion. This combination uses the highly accurate particle displacements that can be efficiently achieved by modern N-body methods with a symmetries-based bias expansion which describes the clustering of any tracer on large scales. We show that at low redshifts, and for moderately biased tracers, the substitution of N-body-determined dynamics improves over an equivalent model using perturbation theory by more than a factor of two in scale, while at high redshifts and for highly biased tracers the gains are more modest. This hybrid approach lends itself well to emulation. By removing the need to identify haloes and subhaloes, and by not requiring any galaxy-formation-related parameters to be included, the emulation task is significantly simplified at the cost of modelling a more limited range in scale.

79 ASTRONOMY AND ASTROPHYSICS↗

Dividing the Concentrator Target From the Genesis Mission

The Genesis spacecraft, launched in 2001, traveled to a Lagrangian point between the Earth and Sun to collect particles from the solar wind and return them to Earth. However, during the return of the spacecraft in 2004, the parachute failed to open during descent, and the Genesis spacecraft crashed into the Utah desert. Many of the solar wind collectors were broken into smaller pieces, and the field team rapidly collected the capsule and collector pieces for later assessment. On each of the next few days, the team discovered that various collectors had survived intact, including three of four concentrator targets. Within a month, the team had imaged more than 10,000 fragments and packed them for transport to the Astromaterials Acquisition and Curation Office within the ARES Directorate at JSC. Currently, the Genesis samples are curated along with the other extraterrestrial sample collections within ARES. Although they were broken and dirty, the Genesis solar wind collectors still offered the science community the opportunity to better understand our Sun and the solar system as a whole. One of the more highly prized concentrator collectors survived the crash almost completely intact. The Genesis Concentrator was designed to concentrate the solar wind by a factor of at least 20 so that solar oxygen and nitrogen isotopes could be measured. One of these materials was the Diamond-on-Silicon (DoS) concentrator target. Unfortunately, the DoS concentrator broke on impact. Nevertheless, the scientific value of the DoS concentrator target was high. The Genesis Allocation Committee received a request for approximately 1 cm(sup 2) of the DoS specimen taken near the focal point of the concentrator for the analysis of solar wind nitrogen isotopes. The largest fragment, Genesis sample 60000, was designated for this allocation and needed to be precisely cut. The requirement was to subdivide the designated sample in a manner that prevented contamination of the sample and minimized the risk of losing or breaking the precious requested sample fragment. The Genesis curator determined that the use of laser scribing techniques to "cut" a precise line and subsequently cleave the sample (in a controlled break of the sample along that line) was the best method for accomplishing the sample subdivision. However, there were risks, including excess heating of the sample, that could cause some of the implanted solar wind to be lost via thermal diffusion. Accidentally breaking the sample during the handling and cleaving process was an additional risk. Early in fiscal year 2013, to address this delicate, complicated task, the ARES Directorate assembled its top scientists to develop a cutting plan that would ensure success when applied to the actual concentrator target wafer; i.e., to produce an approximately 1 cm(sup 2) piece from the requested area of the wafer. The team, subsequently referred to as the JSC Genesis Tiger Team, spent months researching and testing parameters and techniques related to scribing, cleaving, transporting, handling, and holding (i.e., mounting) the specimen. The investigation required considerable "thinking outside the box," and many, many trials using nonflight wafer analogs. After all preliminary testing, the following method was adopted as the final cutting plan. It was used in two final end-to-end practice runs before being used on the actual flight target wafer. The wafer was oriented on the laser cutting stage with the 100 and 010 directions of the wafer parallel to the corresponding X and Y directions of the cutting stage. The laser was programed to scribe 31 lines of the appropriate length along the Y stage direction. The programed scribe lines were separated by 5 micron in the X direction. The laser parameters were set as follows: (1) The laser power was 0.5 watts; (2) each line consisted of 50 passes, with the Z position being advanced 5 micron per pass; and (3) 30 s would elapse before the next line was scribed to allow for wafer cool down from any possible heating via the laser. The ablated material that "stuck" in the "scribe-cut" was removed from the "cut" using an ultrasonic micro-tool. After all the ablated silicon was removed from the wafer, the wafer was repositioned in exactly the same orientation on the laser stage. The laser was focused using the bottom of the wafer channel, and the 31-line scribing pattern described above was reprogrammed using the Z position of the groove bottom as the starting Z value instead of the top wafer surface, which was used previously. Upon completion of the second set of scribes, the ultrasonic micro-tool was again used to clean out the cut. The wafer was remounted on the stage in exactly the same orientation as before. The laser was again focused on the bottom of the groove. This time, however, the laser was.programed to scribe only one line down the exact center of the channel. The final scribe line consisted of 100 passes with a Z advance of 5 micron per pass and with the laser power set at 0.5 watts. As mentioned above, the final cutting plan was practiced in two end-to-end trials using non-flight, triangular-shaped silicon wafers similar in size and orientation to the actual DOS 60000 target sample. The actual scribing of the triangular-shaped wafers required scribing two lines and cleaving (i.e. scribe-cleave, then scribe-cleave) to obtain the piece requested for allocation. Early in December 2012, after many months of experiments and practicing and perfecting the techniques and procedures, the team successfully subdivided the Genesis DoS 60000 target sample, one of the most scientifically important samples from the Genesis mission (figure 2). On December 17, 2012, the allocated piece of concentrator target sample was delivered to the requesting principal investigator.The cutting plan developed for the subdivision of this sample will be used as the model for subdividing future requested Genesis flight wafers (appropriately modified for different wafer types).

Lauer, H. V., Jr.↗

Statistical properties of homogeneous and isotropic turbulence in He II measured via particle tracking velocimetry

Despite being a quantum two-fluid system, superfluid helium-4 (He II) is observed to behave similarly to classical fluids when a flow is generated by mechanical forcing. This similarity has brought up the feasibility of utilizing He II for high Reynolds number classical turbulence research, considering the small kinematic viscosity of He II. However, it has been suggested that the nonclassical dissipation mechanism in He II at small scales may alter its turbulent statistics and intermittency. In this work, we report our study of a nearly homogeneous and isotropic turbulence (HIT) generated by a towed grid in He II. We measure the velocity field using particle tracking velocimetry with solidified deuterium particles as the tracers. By correlating the velocities measured simultaneously on different particle trajectories or at different times along the same particle trajectory, we are able to conduct both Eulerian and Lagrangian flow analyses. Spatial velocity structure functions obtained through the Eulerian analysis show scaling behaviors in the inertial subrange similar to that for classical HIT but with enhanced intermittency. The Lagrangian analysis allows us to examine the flow statistics down to below the dissipation length scale. Interestingly, abnormal deviations from the classical scaling behaviors are observed in this regime. Lastly, we discuss how these deviations may relate to the motion of quantized vortices in the superfluid component in He II.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Gas-Particle Interaction Model Development in Plume Surface Interaction Erosion and Cratering

The development of a predictive simulation capability for extra-terrestrial Plume Surface Interaction (PSI) environments is undertaken by the Fluid Dynamics Branch at NASA/MSFC under the Game Changing Development (GCD) Program funded by NASA’s Space Technology Mission Directorate (STMD). Predictive simulation capability for propulsive Martian and Lunar landing PSI requires accurate modeling of the complex high-speed plume impingement and resulting gas/particle-cloud and ejecta stream formation. Physics modeling gap analysis during formulation of the PSI project scope identified two particular phenomena of first order importance in gas-particle interactions that lacked existing modeling formulations. The first arises from the lack of models capturing the unsteady drag forces imparted on particles by the rapidly expanding supersonic compressible flow. The second addresses the fluctuating forces and dispersions imparted on both the particle and the fluid resulting from the interference and wake turbulence generated by close proximity particle clouds, dubbed the Particle Turbulent Kinetic Energy(PTKE). Their first order significance has been identified in experiments, but simulation models currently do not exist for either effect. The development of models and the eventual implementation into the Eulerian Gas-Granular Flow Solver (GGFS) simulation tools was constructed as one of four major tasks of the PSI project. In this process, data on particle kinetics and gas-particle interactions are collected from carefully designed experiments of particles embedded in jets. The effects observed in the experiments are then replicated in high-fidelity particle-resolved CFD simulations to inform the formulation of improvements to particle phase drag models for implementations in the more efficient Eulerian-Lagrangian CFD simulations. The resulting models are ultimately ported to the Eulerian-Eulerian models applied for most efficient simulations in PSI production application tools. This paper describes the activities and accomplishments of the past year in the gas-particle interaction modeling task of the PSI project.

Jeff West↗

Observations and Parameterizations of Particle Size Distributions in Deep Tropical Cirrus and Stratiform Precipitation Clouds: Results from In-Situ Observations in TRMM Field Campaigns

In this study, we report on the evolution of particle size distributions (PSDs) and habits as measured during slow, Lagrangian-type spiral descents through deep subtropical and tropical cloud layers in Florida, Brazil, and Kwajalein, Marshall Islands, most of which were precipitating. The objective of the flight patterns was to learn more about how the PSDs evolved in the vertical and to obtain information of the vertical structure of microphysical properties. New instrumentation yielding better information on the concentrations of particles in the size (D) range between 0.2 and 2 cm, as well as improved particle imagery, produced more comprehensive observations for tropical stratiform precipitation regions and anvils than have been available previously. Collocated radar observations provided additional information on the vertical structure of the cloud layers sampled. Most of the spirals began at cloud top, with temperatures (T) as low as -50 C, and ended at cloud base or below the melting layer (ML). The PSDs broadened from cloud top towards cloud base, with the largest particles increasing in size from several millimeters at cloud top to one centimeter or larger towards cloud base. Some continued growth was noted in the upper part of the ML. Concentrations of particles less than 1 mm in size decreased with decreasing height. The result was a consistent change in the PSDs in the vertical. Similarly, systematic changes in the size dependence of the particle cross-sectional area was noted with decreasing height. Aggregation-as ascertained from both the changes in the PSDs and evolution of particle habits as observed in high detail with the cloud particle imager (CPI) probe-was responsible for these trends. The PSDs were generally well-represented by gamma distributions of the form N = N0 gamma D microns e- lambda gamma D that were fitted to the PSDs over 1-km horizontal intervals throughout the spirals. The intercept (N0 gamma), slope (lambda gamma), and dispersion (microns) values were derived for each PSD. Exponential curves (N = N0e- lambdaD; micron = 0) were also fitted to the distributions. The lambda gamma values for given spirals varied systematically with temperature as did the values of lambda (exponential), and the data generally conformed to values found in previous studies involving exponential fits to size distributions in mid-latitude frontal and cirrus layers. Considerable variability often noted in the PSD properties during the loops of individual spirals was manifested primarily in large changes in N0 gamma and N0, but micron, lambda gamma and lambda remained fairly stable. Temperature is not found to be the sole factor controlling lambda gamma or lambda but is a primary one. Direct relationships were found between lambda gamma and N0 gamma or lambda gamma and micron for the gamma distributions and lambda and N0 for the exponential. The latter relationship was not found as distinctly in earlier studies; observed PSDs in this study had better fidelity with less scatter. The micron values changed monotonically with T over the range of temperatures and were directly related to N0 gamma or lambda gamma, thereby reducing the number of variables in the PSD functional equation to two. In the upper part of the ML, N0 and lambda continued to decrease, and in the lower part these values began to increase as the largest particles melted. We developed general expressions relating various bulk microphysical, radar, and radiative transfer-related variables to N0 gamma and lambda gamma, useful for both tropical and mid-latitude clouds. These relationships facilitate the specification of a number of bulk properties in cloud and climate models. The results presented in this paper apply best to temperatures between 0 and -40 C, for which the measured radar reflectivities fall in the range of 0 to 25 dBZe.

Heymsfield, Andrew J.↗

Spatial and Temporal Extrapolation of Disdrometer Size Distributions Based on a Lagrangian Trajectory Model of Falling Rain

Methodologies to improve disdrometer processing, loosely based on mathematical techniques common to the field of particle flow and fluid mechanics, are examined and tested. The inclusion of advection and vertical wind field estimates appear to produce significantly improved results in a Lagrangian hydrometeor trajectory model, in spite of very strict assumptions of noninteracting hydrometeors, constant vertical air velocity, and time independent advection during the scan time interval. Wind field data can be extracted from each radar elevation scan by plotting and analyzing reflectivity contours over the disdrometer site and by collecting the radar radial velocity data to obtain estimates of advection. Specific regions of disdrometer spectra (drop size versus time) often exhibit strong gravitational sorting signatures, from which estimates of vertical velocity can be extracted. These independent wind field estimates become inputs and initial conditions to the Lagrangian trajectory simulation of falling hydrometeors.

Lane, John E.↗

Light Scalars at FASER

FASER, the ForwArd Search ExpeRiment, is a currently operating experiment at the Large Hadron Collider (LHC) that can detect light long-lived particles produced in the forward region of the LHC interacting point. In this paper, we study the prospect of detecting light CP-even and CP-odd scalars at FASER and FASER 2. Considering a model-independent framework describing the most general interactions between a CP-even or CP-odd scalar and SM particles using the notation of coupling modifiers in the effective Lagrangian, we develop the general formalism for the scalar production and decay. We then analyze the FASER and FASER 2 reaches of light scalars in the large tan β region of the Type-I two Higgs double model as a case study, in which light scalars with relatively long lifetime could be accommodated. In the two benchmark scenarios we considered, the light (pseudo)scalar decay length varies in (10 –8 , 10 5 ) meters. Both FASER and FASER 2 can probe a large part of the parameter space in the large tan β region up to 10 7 , extending beyond the constraints of the other existing experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles is investigated using a novel transient three-dimensional computational fluid dynamics model. The model couples particle–gas hydrodynamics of solid–gas flows laden with polydisperse particles, radiative heat transfer in non-grey absorbing, emitting and anisotropically-scattering multi-component participating media, conduction heat transfer in the gas phase, and interfacial convection heat transfer. The multiphase particle-in-cell method is used to predict high-fidelity solid–gas flow characteristics, such as the local discrete particle size distribution, with increased computational efficiency by combining the advantages of both Eulerian and Lagrangian methods. The multi-component radiative transfer model is implemented using an advanced collision-based Monte Carlo ray-tracing method. The number of the prescribed discrete particle components is found to be the key parameter affecting the computational accuracy and efficiency, which primarily depends on the size distribution of the particles. For the model particle–gas flow featuring free-falling Gamma-distributed ceramic particles exposed to concentrated solar irradiation, the particle volume fraction, radiative, fluid flow and thermal characteristics appear to converge with the increasing number of the discrete particle components. Five particle components are sufficient to obtain physically meaningful results. A further increase in the number of the particle components only slightly increases the accuracy of the numerical predictions at the expense of a rapidly increasing computational time. For five particle components, the particle vertical velocity at the receiver exit for particles with the diameter of 43.4 μm is 57% of that for the particles with the diameter of 202.8 μm. The temperatures of these two particle components increase from the initial ambient values by factors of 2 and 1.2, respectively, during the simulation time. In conclusion, the model developed allows for increased fidelity of particle–gas flow simulations with significant radiative effects.

97 MATHEMATICS AND COMPUTING↗

Improved treatment of multi-material cells in thermal radiation transport codes

High-energy-density physics simulations with non-conformal meshes of materials require a multi-material (MM) closure that affects the thermal radiation transport (TRT). We propose a set of novel closures that work for an arbitrary number of materials, both grey and multigroup energy discretizations, and any angular discretization (such as Sn, IMC, or diffusion). For each spatial cell, our closures let each species (ion and electron) of each material have its own temperature, density, and internal energy, but use a single radiation distribution that interacts with all materials within the cell. Our closures maintain energy conservation, do not incur increased computational cost in the TRT solve itself, are compatible with single-material TRT solvers, and do not make any temperature-equilibrium assumptions. Here we test our closures on a wide range of increasingly realistic problems and find them to be robust.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Towards a theory of hadron resonances

In this review, we present the current state of the art of our understanding of the spectrum of excited strongly interacting particles and discuss methods that allow for a systematic and model-independent calculation of the hadron spectrum. These are lattice QCD and effective field theories. Synergies between both approaches can be exploited allowing for a deeper understanding of the hadron spectrum. Results based on effective field theories and hadron–hadron scattering in lattice QCD or combinations thereof are presented and discussed. We also show that the often used Breit–Wigner parametrization is at odds with chiral symmetry and should not be used in case of strongly coupled channels.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A macroscopic plasma Lagrangian and its application to wave interactions and resonances

The derivation of a macroscopic plasma Lagrangian is considered, along with its application to the description of nonlinear three-wave interaction in a homogeneous plasma and linear resonance oscillations in a inhomogeneous plasma. One approach to obtain the Lagrangian is via the inverse problem of the calculus of variations for arbitrary first and second order quasilinear partial differential systems. Necessary and sufficient conditions for the given equations to be Euler-Lagrange equations of a Lagrangian are obtained. These conditions are then used to determine the transformations that convert some classes of non-Euler-Lagrange equations to Euler-Lagrange equation form. The Lagrangians for a linear resistive transmission line and a linear warm collisional plasma are derived as examples. Using energy considerations, the correct macroscopic plasma Lagrangian is shown to differ from the velocity-integrated low Lagrangian by a macroscopic potential energy that equals twice the particle thermal kinetic energy plus the energy lost by heat conduction.

Peng, Y. K. M.↗

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↗

Evaluation of an efficient statistical transport model for turbulent droplet dispersion in dilute combusting sprays

Evaluation of an efficient statistical transport model for turbulent droplet dispersion is made for a dilute spray of methanol droplets injected into a turbulent, axisymmnetric methane-fuelled diffusion flame burning in stagnant air. In the dispersion model, a computational parcel representing a group of physical particles (droplets) is considered to have a normal probability density function (pdf) in space. The mean is determined by Lagrangian tracking through a sequence of stochastically generated turbulent eddies and the variance is evaluated from a statistical formulation based on the linearized particle equations of motion. The basic validity of this model is established through comparison with available experimental data and with theoretical predictions using a conventional stochastic direct modeling approach. The conclusion of the evaluation is that the proposed dispersion model compares favorably with experimental data and provides a valid technique for simulating turbulent combusting sprays with significant computational savings over conventional methods.

Litchford, Ron J.↗

Fluxes of MeV particles at Earth's orbit and their relationship with the global structure of the solar corona: Observations from SOHO

The SOHO satellite, launched on 2 December 1995, combines a unique set of instruments which allow comparative studies of the interior of the sun, the outer corona and solar to be carried out. In its halo orbit around the L1 Lagrangian point of the sun-earth system, SOHO's comprehensive suprathermal and energetic particle analyzer (COSTEP) measures in situ energetic particles in the energy range of 44 keV/particle to greater than 53 MeV/n. The MeV proton, electron and helium nuclei measurements from the COSTEP electron proton helium instrument (EPHIN) were used to investigate the relationships of intensity increases of these particle species with the large-scale structures of the solar corona and heliosphere, including temporal variations. Coronal observatons are provided by the large angle spectroscopic coronagraph (LASCO) and the extreme ultraviolet imaging telescope (EIT). It was found that during times of minimum solar activity, intensity increases of the particles have two well defined sources: corotating interaction regions (CIRs) in the heliosphere related to coronal holes at the sun and coronal mass ejections.

Posner, A.↗

Lagrangian methods in plasma dynamics. II - Construction of Lagrangians for plasmas

Consideration of the construction of suitable Lagrangian functions for the dynamics of a cold plasma in such a way as to retain the relativistically covariant formalism. In one method, this is achieved by the introduction of a set of three variables which label the world lines of the particles. A second method results in a Clebsch-type representation. Sturrock's relativistic Lagrangian and Low's hot plasma Lagrangian are also briefly discussed in the context of the present work. The behavior of the canonical stress tensor is considered. The applicability of many of the general results in part I (Dougherty, 1970) is ensured by establishing the existence of the Lagrangian function.

Dougherty, J. P.↗