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 127 records · Page 7

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.↗

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.↗

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.↗

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.↗

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.↗

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↗

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↗

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.↗

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.↗

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.↗

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.↗

Evolution of light domain walls interacting with dark matter, part 1

The evolution of domain walls generated in the early Universe is discussed considering an interaction between the walls and a major gaseous component of the dark matter. The walls are supposed able to reflect the particles elastically and with a reflection coefficient of unity. A toy Lagrangian that could give rise to such a phenomenon is discussed. In the simple model studied, highly non-relativistic and slowly varying speeds are obtained for the domain walls (approximately 10 (exp -2)(1+z)(exp -1)) and negligible distortions of the microwave background. In addition, these topological defects may provide a mechanism of forming the large scale structure of the Universe, by creating fluctuations in the dark matter (delta rho/rho approximately O(1)) on a scale comparable with the distance the walls move from the formation (in the model d less than 20 h(exp -1) Mpc). The characteristic scale of the wall separation can be easily chosen to be of the order of 100 Mpc instead of being restricted to the horizon scale, as usually obtained.

Massarotti, Alessandro↗