Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lagrangian”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

An investigation of turbulent transport in the extreme lower atmosphere

A model in which the Lagrangian autocorrelation is expressed by a domain integral over a set of usual Eulerian autocorrelations acquired concurrently at all points within a turbulence box is proposed along with a method for ascertaining the statistical stationarity of turbulent velocity by creating an equivalent ensemble to investigate the flow in the extreme lower atmosphere. Simultaneous measurements of turbulent velocity on a turbulence line along the wake axis were carried out utilizing a longitudinal array of five hot-wire anemometers remotely operated. The stationarity test revealed that the turbulent velocity is approximated as a realization of a weakly self-stationary random process. Based on the Lagrangian autocorrelation it is found that: (1) large diffusion time predominated; (2) ratios of Lagrangian to Eulerian time and spatial scales were smaller than unity; and, (3) short and long diffusion time scales and diffusion spatial scales were constrained within their Eulerian counterparts.

Koper, C. A., Jr.↗
NASA NTRS1975

Turbulence in a gaseous hydrogen-liquid oxygen rocket combustion chamber

The intensity of turbulence and the Lagrangian correlation coefficient for a LOX-GH2 rocket combustion chamber was determined from experimental measurements of tracer gas diffusion. A combination of Taylor's turbulent diffusion theory and a numerical method for solving the conservation equations of fluid mechanics was used to calculate these quantities. Taylor's theory was extended to consider the inhomogeneity of the turbulence field in the axial direction of the combustion chamber, and an exponential function was used to represent the Lagrangian correlation coefficient. The results indicate that the value of the intensity of turbulence reaches a maximum of 14% at a location about 7" downstream from the injector. The Lagrangian correlation coefficient associated with this value is given by the above exponential expression where alpha = 10,000/sec.

Lebas, J.↗
NASA NTRS1978

Measurement of the turbulence in a gaseous rocket combustion chamber

The intensity of turbulence and the Lagrangian correlation coefficient in a gaseous hydrogen (GH2)-gaseous oxygen (GOX) rocket combustion chamber are studied. A helium tracer gas was injected at an upstream point on the combustion chamber centerline. Samples collected downstream were analyzed for the spreading of the tracer gas concentration. Using Taylor's turbulent diffusion theory and Gosman's numerical method, turbulence parameters were found. In order to evaluate turbulence field inhomogeneity in the axial direction, Taylor's theory was extended. An exponential function was used to approximate the Lagrangian correlation coefficient. It was found that turbulence intensity varies from 15% near the injector to 4% at the nozzle entrance, and that the Lagrangian coefficient rapidly decreases when the dispersion time difference increases.

Tou, P. P.↗
NASA NTRS1978

On the mean meridional mass motions of the stratosphere and mesosphere

Using a simplified, approximate 'Lagrangian-mean' dynamical formulation, the mean meridional mass circulation of the stratosphere and mesosphere is discussed. Under solstice conditions, it is shown that this Lagrangian-mean circulation may be inferred, as a first approximation, from the Eulerian-mean diabatic heating. Diabatic heating rates for the solstices, originally derived by Murgatroyd and Goody (1958), result in Lagrangian-mean rising motion at the tropical tropopause, subsidence across the extra-tropical tropopause, and a very strong summer-to-winter pole flow in the mesosphere. This circulation is exactly that obtained by Murgatroyd and Singleton (1961) for the solstices. Those authors, however, attempted to identify this circulation as the Eulerian-mean motion, and were later criticized for their neglect of the meridional eddy heat flux in the calculation, which proved to be extremely important in the winter hemisphere. The present study, nevertheless, indicates that Murgatroyd and Singleton's circulation may in fact be representative of actual air parcel motions in the stratosphere and mesosphere.

Dunkerton, T.↗
NASA NTRS1982

Generalized Roche potential for misaligned binary systems - Properties of the critical lobe

The paper considers the Roche potential for binary systems where the stellar rotation axis is not aligned with the orbital revolution axis. It is shown that, as the degree of misalignment varies, internal Lagrangian points and external Lagrangian points may switch their roles. A systematic method to identify the internal Lagrangian point and to calculate the volume of the critical lobe is developed, and numerical results for a wide range of parameters of binary systems with circular orbits are presented. For binary systems with large enough misalignment, discrete changes occur in the topological structure of the equipotential surfaces as the orbital phase varies. The volume of the critical lobe has minima, as a function of orbital phase, at the two instances when the secondary crosses the equatorial plane of the primary. In semidetached systems, mass transfer may be confined to the vicinity of these two instances.

Avni, Y.↗
NASA NTRS1993

An Eulerian/Lagrangian coupling procedure for three-dimensional vortical flows

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of 3D vortical flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method, added to the Eulerian time-marching procedure, provides a correction of the Eulerian solution. In turn, the Eulerian solution is used to integrate the Lagrangian state-vector along the particles trajectories. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers describe accurately the convection properties and enhance the vorticity and entropy capturing capabilities of the Eulerian solver. The Eulerian/Lagrangian coupling strategies are discussed and the combined scheme is tested on a constant stagnation pressure flow in a 90 deg bend and on a swirling pipe flow. As the numerical diffusion is reduced when using the Lagrangian correction, a vorticity gradient augmentation is identified as a basic problem of this inviscid calculation.

Felici, Helene M.↗
NASA NTRS1993

Reduction of numerical diffusion in three-dimensional vortical flows using a coupled Eulerian/Lagrangian solution procedure

A new approach based on the coupling of an Eulerian and a Lagrangian solver, aimed at reducing the numerical diffusion errors of standard Eulerian time-marching finite-volume solvers, is presented. The approach is applied to the computation of the secondary flow in two bent pipes and the flow around a 3D wing. Using convective point markers the Lagrangian approach provides a correction of the basic Eulerian solution. The Eulerian flow in turn integrates in time the Lagrangian state-vector. A comparison of coarse and fine grid Eulerian solutions makes it possible to identify numerical diffusion. It is shown that the Eulerian/Lagrangian approach is an effective method for reducing numerical diffusion errors.

Felici, Helene M.↗
NASA NTRS1953

Diffusion of Heat from a Line Source in Isotropic Turbulence

An experimental and analytical study has been made of some features of the turbulent heat diffusion behind a line heated wire stretched perpendicular to a flowing isotropic turbulence. The mean temperature distributions have been measured with systematic variations in wind speed, size of turbulence-producing grid, and downstream location of heat source. The nature of the temperature fluctuation field has been studied. A comparison of Lagrangian and Eulerian analyses for diffusion in a nondecaying turbulence yields an expression for turbulent-heat-transfer coefficient in terms of turbulence velocity and a Lagrangian "scale." the ratio of Eulerian to Lagrangian microscale has been determined theoretically by generalization of a result of Heisenberg and with arbitrary constants taken from independent sources, shows rough agreement with experimental results. A convenient form has been deduced for the criterion of interchangeability of instantaneous space and time derivatives in a flowing turbulence.

Uberoi, Mahinder S↗
NASA NTRS1994

Numerical methods for multibody systems

This article gives a brief summary of some results obtained by Nasser on modeling and simulation of inequality problems in multibody dynamics. In particular, the augmented Lagrangian method discussed here is applied to a constrained motion problem with impulsive inequality constraints. A fundamental characteristic of the multibody dynamics problem is the lack of global convexity of its Lagrangian. The problem is transformed into a convex analysis problem by localization (piecewise linearization), where the augmented Lagrangian has been successfully used. A model test problem is considered and a set of numerical experiments is presented.

Glowinski, Roland↗
NASA NTRS1994

Mixing processes within the polar night jet

Lagrangian material line simulations are performed using U.K. Meteorological Office simulated winds and temperatures to examine mixing processes in the middle- and lower-stratospheric polar night jet during the 1992 Southern Hemisphere spring and Northern Hemisphere winter. The Lagrangian simulations are undertaken to provide insight into the effects of mixing within the polar night jet on observations of the polar vortex made by instruments onboard the Upper Atmosphere Research Satellite (UARS) during these periods. A moderate to strong kinematic barrier to large-scale isentropic exchange, similar to the barrier identified in General Circulation Model (GCM) simulations, is identified during both of these periods. Characteristic timescales for mixing by large-scale isentropic motions within the polar night jet range from 20 days in the Southern Hemisphere lower stratosphere to years in the Northern Hemisphere middle stratosphere. The long mixing timescales found in the Northern Hemisphere polar night jet do not persist. Instead, the Northern Hemisphere kinematic barriers are broken down as part of the large-scale stratospheric response to a strong tropospheric blocking event. A series of Lagrangian experiments are conducted to investigate the sensitivity of the kinematic barrier to diabatic effects and to small-scale inertial gravity wave motions. Differential diabatic descent is found to have a significant impact on mixing processes within the Southern Hemisphere middle-stratospheric jet core. The interaction between small-scale displacements by idealized, inertial gravity waves and the large-scale flow is found to have a significant impact on mixing within the polar night jet in both hemispheres. These sensitivity experiments suggest that scales of motion that are unresolved in global assimilated datasets may contribute to mass exchange across the kinematic barrier to large-scale isentropic motion.

Pierce, R. Bradley↗
NASA NTRS1997

Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design

Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.

Fahrenthold, Eric P.↗
NASA NTRS1999

Dynamical and Chemical Behavior of the Lower Stratosphere and Interactions with the Troposphere

Equivalent-barotropic calculations, in tandem with Lagrangian analyses, reveal how changes of total ozone follow from vertical and horizontal transport by planetary waves. Those calculations also throw light on how diabatic motions comprising the Brewer-Dobson circulation develop from quasi-horizontal advection by planetary waves. Potential temperature along a material surface indicates organized subsidence inside the polar-night vortex, resembling tracer observations from UARS. Lagrangian histories illustrate that this sinking motion follows in large part from parcels being driven out of thermodynamic equilibrium by planetary waves, especially at high latitudes. Irreversible heat transfer then produces a net drift of air across isentropic surfaces as parcels orbit about the displaced vortex. By driving mean-meridional overturning in the stratosphere, this downward drift is ultimately responsible for transferring ozone from the tropics to the extratropical lower stratosphere. It also introduces horizontal structure into the distribution of total ozone, which surfaces clearly in ozone trends. High-resolution global cloud imagery constructed from 6 satellites simultaneously observing the Earth was used to investigate the spectrum of equatorial waves generated by tropical convection and propagating vertically into the stratosphere. The results indicate that temperature variability is dominated by planetary-scale equatorial waves like the Kelvin mode, which agrees with satellite observations of the tropical stratosphere. However, the Kelvin mode accounts for only about 30 - 50% of the eastward momentum flux radiating into the stratosphere, the remainder coming from gravity waves. An algorithm was developed to determine 3-dimensional atmospheric motion from satellite tracer measurements. Based on Lagrangian constraints, the algorithm circumvents limitations of the traditional scheme for inferring motion from temperature measurements and determines the circulation in the tropics as reliably as elsewhere. A study of deep convection revealed that the highest towers (those penetrating into stratospheric air and controlling tropopause height and composition through convective mixing) occur in close association with the diurnal cycle of convection. Clouds colder than 220 K develop almost entirely in association with the diurnal cycle of convection over tropical landmasses and substantially in association with it even over maritime regions.

Salby, Murry L.↗
NASA NTRS2003

Water Impact Test and Simulation of a Composite Energy Absorbing Fuselage Section

In March 2002, a 25-ft/s vertical drop test of a composite fuselage section was conducted onto water. The purpose of the test was to obtain experimental data characterizing the structural response of the fuselage section during water impact for comparison with two previous drop tests that were performed onto a rigid surface and soft soil. For the drop test, the fuselage section was configured with ten 100-lb. lead masses, five per side, that were attached to seat rails mounted to the floor. The fuselage section was raised to a height of 10-ft. and dropped vertically into a 15-ft. diameter pool filled to a depth of 3.5-ft. with water. Approximately 70 channels of data were collected during the drop test at a 10-kHz sampling rate. The test data were used to validate crash simulations of the water impact that were developed using the nonlinear, explicit transient dynamic codes, MSC.Dytran and LS-DYNA. The fuselage structure was modeled using shell and solid elements with a Lagrangian mesh, and the water was modeled with both Eulerian and Lagrangian techniques. The fluid-structure interactions were executed using the fast general coupling in MSC.Dytran and the Arbitrary Lagrange-Euler (ALE) coupling in LS-DYNA. Additionally, the smooth particle hydrodynamics (SPH) meshless Lagrangian technique was used in LS-DYNA to represent the fluid. The simulation results were correlated with the test data to validate the modeling approach. Additional simulation studies were performed to determine how changes in mesh density, mesh uniformity, fluid viscosity, and failure strain influence the test-analysis correlation.

Fasanella, Edwin L.↗
NASA NTRS2004

Adaptive, Distributed Control of Constrained Multi-Agent Systems

Product Distribution (PO) theory was recently developed as a broad framework for analyzing and optimizing distributed systems. Here we demonstrate its use for adaptive distributed control of Multi-Agent Systems (MASS), i.e., for distributed stochastic optimization using MAS s. First we review one motivation of PD theory, as the information-theoretic extension of conventional full-rationality game theory to the case of bounded rational agents. In this extension the equilibrium of the game is the optimizer of a Lagrangian of the (Probability dist&&on on the joint state of the agents. When the game in question is a team game with constraints, that equilibrium optimizes the expected value of the team game utility, subject to those constraints. One common way to find that equilibrium is to have each agent run a Reinforcement Learning (E) algorithm. PD theory reveals this to be a particular type of search algorithm for minimizing the Lagrangian. Typically that algorithm i s quite inefficient. A more principled alternative is to use a variant of Newton's method to minimize the Lagrangian. Here we compare this alternative to RL-based search in three sets of computer experiments. These are the N Queen s problem and bin-packing problem from the optimization literature, and the Bar problem from the distributed RL literature. Our results confirm that the PD-theory-based approach outperforms the RL-based scheme in all three domains.

Bieniawski, Stefan↗
NASA NTRS2004

Product Distribution Theory for Control of Multi-Agent Systems

Product Distribution (PD) theory is a new framework for controlling Multi-Agent Systems (MAS's). First we review one motivation of PD theory, as the information-theoretic extension of conventional full-rationality game theory to the case of bounded rational agents. In this extension the equilibrium of the game is the optimizer of a Lagrangian of the (probability distribution of) the joint stare of the agents. Accordingly we can consider a team game in which the shared utility is a performance measure of the behavior of the MAS. For such a scenario the game is at equilibrium - the Lagrangian is optimized - when the joint distribution of the agents optimizes the system's expected performance. One common way to find that equilibrium is to have each agent run a reinforcement learning algorithm. Here we investigate the alternative of exploiting PD theory to run gradient descent on the Lagrangian. We present computer experiments validating some of the predictions of PD theory for how best to do that gradient descent. We also demonstrate how PD theory can improve performance even when we are not allowed to rerun the MAS from different initial conditions, a requirement implicit in some previous work.

Lee, Chia Fan↗
NASA NTRS1999

Acoustic Energy Estimates in Inhomogeneous Moving Media

In ducted fan engine noise research, there is a need for defining a simple and easy to use acoustic energy conservation law to help in quantification of noise control techniques. There is a well known conservation law relating acoustic energy and acoustic energy flux in the case of an isentropic irrotational flow. Several different approaches have been taken to generalize this conservation law. For example, Morfey finds an identity by separating out the irrotational part of the perturbed flow. Myers is able to find a series of indentities by observing an algebraic relationship between the basic conservation of energy equation for a background flow and the underlying equations of motion. In an approximate sense, this algebraic relationship is preserved under perturbation. A third approach which seems to have not been pursued in the literature is a result known as Noether's theorem. There is a Lagrangian formulation for the Euler equation of fluid mechanics. Noether's theorem says that any group action that leaves the Lagrangian action invariant leads to a conserved quantity. This presentation will include a survey of current results regarding acoustic energy and preliminary results on the symmetries of the Lagrangian.

Farassat, F.↗
NASA NTRS2004

A Survey Of Earth-Moon Libration Orbits: Stationkeeping Strategies And Intra-Orbit Transfers

Cislunar space is a readily accessible region that may well develop into a prime staging area in the effort to colonize space near Earth or to colonize the Moon. While there have been statements made by various NASA programs regarding placement of resources in orbit about the Earth-Moon Lagrangian locations, there is no survey of the total cost associated with attaining and maintaining these unique orbits in an operational fashion. Transfer trajectories between these orbits required for assembly, servicing, and positioning of these resources have not been extensively investigated. These orbits are dynamically similar to those used for the Sun-Earth missions, but differences in governing gravitational ratios and perturbation sources result in unique characteristics. We implement numerical computations using high fidelity models and linear and nonlinear targeting techniques to compute the various maneuver (Delta)V and temporal costs associated with orbits about each of the Earth-Moon Lagrangian locations (L1, L2, L3, L4, and L5). From a dynamical system standpoint, we speak to the nature of these orbits and their stability. We address the cost of transfers between each pair of Lagrangian locations.

Folta, David↗
NASA NTRS2006

Simulations of Sea-Ice Dynamics Using the Material-Point Method

In recent years, the availability of large volumes of recorded ice motion derived from high-resolution SAR data has provided an amazingly detailed look at the deformation of the ice cover. The deformation is dominated by the appearance of linear kinematic features that have been associated with the presence of leads. These remarkable data put us in a position to begin detailed evaluation of current coupled mechanical and thermodynamic models of sea ice. This presentation will describe the material point method (MPM) for solving these model equations. MPM is a numerical method for continuum mechanics that combines the best aspects of Lagrangian and Eulerian discretizations. The material points provide a Lagrangian description of the ice that models convection naturally. Thus, properties such as ice thickness and compactness are computed in a Lagrangian frame and do not suffer from errors associated with Eulerian advection schemes, such as artificial diffusion, dispersion, or oscillations near discontinuities. This desirable property is illustrated by solving transport of ice in uniform, rotational and convergent velocity fields. Moreover, the ice geometry is represented by unconnected material points rather than a grid. This representation facilitates modeling the large deformations observed in the Arctic, as well as localized deformation along leads, and admits a sharp representation of the ice edge. MPM also easily allows the use of any ice constitutive model. The versatility of MPM is demonstrated by using two constitutive models for simulations of wind-driven ice. The first model is a standard viscous-plastic model with two thickness categories. The MPM solution to the viscous-plastic model agrees with previously published results using finite elements. The second model is a new elastic-decohesive model that explicitly represents leads. The model includes a mechanism to initiate leads, and to predict their orientation and width. The elastic-decohesion model can provide similar overall deformation as the viscous-plastic model; however, explicit regions of opening and shear are predicted. Furthermore, the efficiency of MPM with the elastic-decohesive model is competitive with the current best methods for sea ice dynamics. Simulations will also be presented for an area of the Beaufort Sea, where predictions can be validated against satellite observations of the Arctic.

Sulsky, D.↗