SEARCH · Search NASA

Results for “Lagrangian methods”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS1974

Experimental determination of turbulence in a GH2-GOX rocket combustion chamber

The intensity of turbulence and the Lagrangian correlation coefficient for a gaseous rocket combustion chamber have been determined from the experimental measurements of the tracer gas diffusion. A combination of Taylor's turbulent diffusion theory and Spalding's 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. An exponential function was used to represent the Lagrangian correlation coefficient. The results indicate that the maximum value of the intensity of turbulence is about 15% and the Lagrangian correlation coefficient drops to about 0.12 in one inch of the chamber length.

Tou, P.↗

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 NTRS1986

A spectral element method for the simulation of unsteady incompressible flows with heat transfer

The spectral element method is a high-order finite element technique for solution of the Navier-Stokes and energy equations. In the isoparametric spectral element discretization, the domain is broken up into general brick elements, and the dependent and independent variables represented as high-order tensor-product Lagrangian interpolants through Chebyshev collocation points. The nonlinear and convective terms in the governing equations are treated with explicit collocation, while the pressure and diffusive contributions are handled implicitly using variational projection operators. The method is applied to flow past a cylinder, flow in grooved channels, and natural convection in an enclosure.

Karniadakis, George E.↗

NASA NTRS1995

Full-Scale Direct Numerical Simulation of Two- and Three-Dimensional Instabilities and Rivulet Formulation in Heated Falling Films

A thin film draining on an inclined plate has been studied numerically using finite element method. Three-dimensional governing equations of continuity, momentum and energy with a moving boundary are integrated in an arbitrary Lagrangian Eulerian frame of reference. Kinematic equation is solved to precisely update interface location. Rivulet formation based on instability mechanism has been simulated using full-scale computation. Comparisons with long-wave theory are made to validate the numerical scheme. Detailed analysis of two- and three-dimensional nonlinear wave formation and spontaneous rupture forming rivulets under the influence of combined thermocapillary and surface-wave instabilities is performed.

Krishnamoorthy, S.↗

NASA NTRS1994

Eulerian-Lagrangian Simulations of Transonic Flutter Instabilities

This paper presents an overview of recent applications of Eulerian-Lagrangian computational schemes in simulating transonic flutter instabilities. This approach, the fluid-structure system is treated as a single continuum dynamics problem, by switching from an Eulerian to a Lagrangian formulation at the fluid-structure boundary. This computational approach effectively eliminates the phase integration errors associated with previous methods, where the fluid and structure are integrated sequentially using different schemes. The formulation is based on Hamilton's Principle in mixed coordinates, and both finite volume and finite element discretization schemes are considered. Results from numerical simulations of transonic flutter instabilities are presented for isolated wings, thin panels, and turbomachinery blades. The results suggest that the method is capable of reproducing the energy exchange between the fluid and the structure with significantly less error than existing methods. Localized flutter modes and panel flutter modes involving traveling waves can also be simulated effectively with no a priori knowledge of the type of instability involved.

Bendiksen, Oddvar O.↗

NASA NTRS2011

A New Limit on Planck Scale Lorentz Violation from Gamma-ray Burst Polarization

Constraints on possible Lorentz invariance violation (UV) to first order in E/M(sub Plank) for photons in the framework of effective field theory (EFT) are discussed, taking cosmological factors into account. Then. using the reported detection of polarized soft gamma-ray emission from the gamma-ray burst GRB041219a that is indicative' of an absence of vacuum birefringence, together with a very recent improved method for estimating the redshift of the burst, we derive constraints on the dimension 5 Lorentz violating modification to the Lagrangian of an effective local QFT for QED. Our new constraints are more than five orders of magnitude better than recent constraints from observations of the Crab Nebula.. We obtain the upper limit on the Lorentz violating dimension 5 EFT parameter absolute value of zeta of 2.4 x 10(exp -15), corresponding to a constraint on the dimension 5 standard model extension parameter. Kappa (sup 5) (sub (v)oo) much less than 4.2 X 10(exp -3)4 / GeV.

Stecker, Floyd W.↗

NASA NTRS1980

Solutions of contact problems by the assumed stress hybrid model

A method was developed for contact problems which may be either frictional or frictionless and may involve extensive sliding between deformable bodies. It was based on an assumed stress hybrid approach and on an incremental variational principle for which the Euler's equations of the functional include the equilibrium and compatibility conditions at the contact surface. The tractions at an assumed contact surface were introduced as Lagrangian multipliers in the formulation. It was concluded from the results of several example solutions that the extensive sliding contact between deformable bodies can be solved by the present method.

Kubomura, K.↗

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 NTRS2003

Global Aspects of Charged Particle Motion in Axially Symmetric Multipole Magnetic Fields

The motion of a single charged particle in the space outside of a compact region of steady currents is investigated. The charged particle is assumed to produce negligible electromagnetic radiation, so that its energy is conserved. The source of the magnetic field is represented as a point multipole. After a general description, attention is focused on magnetic fields with axial symmetry. Lagrangian dynamical theory is utilized to identify constants of the motion as well as the equations of motion themselves. The qualitative method of Stonner is used to examine charged particle motion in axisymmetric multipole fields of all orders. Although the equations of motion generally have no analytical solutions and must be integrated numerically to produce a specific orbit, a topological examination of dynamics is possible, and can be used, d la Stonner, to completely describe the global aspects of the motion of a single charged particle in a space with an axisymmetric multipole magnetic field.

Shebalin, John V.↗

NASA NTRS1988

On the simultaneous improvement of a satellite orbit and determination of sea surface topography using altimeter data

A method is presented in satellite altimetry that attempts to simultaneously determine the geoid and sea surface toography with minimum wavelengths of about 500 km and to reduce the radial orbit errors caused by geopotential uncertainties. The modeling of these errors is made using the linearized Lagrangian perturbation theory. Observation equations are developed using sea surface heights and crossover discrepancies as observables. A minimum variance solution with prior information can then provide estimates of parametrs representing the sea surface topography and corrections to the orbit. The potential of the method is demonstrated in a solution where simulated geopotential errors and the Levitus sea surface topography are used to generate the observables for a Seasat 3 day arc. The simulation results suggest that the method can be used to efficiently process real altimeter data.

Engelis, Theodossios↗

NASA NTRS1983

Structural optimization - Challenges and opportunities

A review of developments in structural optimization techniques and their interface with growing computer capabilities is presented. Structural design steps comprise functional definition of an object, an evaluation phase wherein external influences are quantified, selection of the design concept, material, object geometry, and the internal layout, and quantification of the physical characteristics. Optimization of a fully stressed design is facilitated by use of nonlinear mathematical programming which permits automated definition of the physics of a problem. Design iterations terminate when convergence is acquired between mathematical and physical criteria. A constrained minimum algorithm has been formulated using an Augmented Lagrangian approach and a generalized reduced gradient to obtain fast convergence. Various approximation techniques are mentioned. The synergistic application of all the methods surveyed requires multidisciplinary teamwork during a design effort.

Sobieszczanski-Sobieski, J.↗

NASA NTRS1984

Structural optimization: Challenges and opportunities

Sobieszczanski-Sobieski, J.↗

NASA NTRS1987

PIC numerical scheme for a two-dimensional Eulerian code

The present work uses the PIC method with operator splitting to hydroelastoplastic multimaterial flow problems in 2D Eulerian coordinates. The operator splitting procedures are carried out in both Lagrangian and remap phases with alternating direction for each time advancement to maintain the accuracy of one-dimensional scheme. A detail description of the numerical approximations is given with a sample calculation of a tungsten rod penetrating an aluminum target.

Kwak, Dochan↗

NASA NTRS1988

Predicting the velocity and azimuth of fragments generated by the range destruction or random failure of rocket casings and tankage

A calculational method is described which provides a powerful tool for predicting solid rocket motor (SRM) casing and liquid rocket tankage fragmentation response. The approach properly partitions the available impulse to each major system-mass component. It uses the Pisces code developed by Physics International to couple the forces generated by an Eulerian-modeled gas flow field to a Lagrangian-modeled fuel and casing system. The details of the predictive analytical modeling process and the development of normalized relations for momentum partition as a function of SRM burn time and initial geometry are discussed. Methods for applying similar modeling techniques to liquid-tankage-overpressure failures are also discussed. Good agreement between predictions and observations are obtained for five specific events.

Eck, Marshall↗

NASA NTRS1993

Non-recursive augmented Lagrangian algorithms for the forward and inverse dynamics of constrained flexible multibodies

A technique is presented for solving the inverse dynamics of flexible planar multibody systems. This technique yields the non-causal joint efforts (inverse dynamics) as well as the internal states (inverse kinematics) that produce a prescribed nominal trajectory of the end effector. A non-recursive global Lagrangian approach is used in formulating the equations for motion as well as in solving the inverse dynamics equations. Contrary to the recursive method previously presented, the proposed method solves the inverse problem in a systematic and direct manner for both open-chain as well as closed-chain configurations. Numerical simulation shows that the proposed procedure provides an excellent tracking of the desired end effector trajectory.

Bayo, Eduardo↗

DOE OSTI2025

Parallel-in-Time Solution of Scalar Nonlinear Conservation Laws

Here, we consider the parallel-in-time solution of scalar nonlinear conservation laws in one spatial dimension. The equations are discretized in space with a conservative finite-volume method using weighted essentially nonoscillatory (WENO) reconstructions, and in time with high-order explicit Runge–Kutta methods. The solution of the global, discretized space-time problem is sought via a nonlinear iteration that uses a novel linearization strategy in cases of nondifferentiable equations. Under certain choices of discretization and algorithmic parameters, the nonlinear iteration coincides with Newton’s method, although, more generally, it is a preconditioned residual correction scheme. At each nonlinear iteration, the linearized problem takes the form of a certain discretization of a linear conservation law over the space-time domain in question. An approximate parallel-in-time solution of the linearized problem is computed with a single multigrid reduction-in-time (MGRIT) iteration; however, any other effective parallel-in-time method could be used in its place. The MGRIT iteration employs a novel coarse-grid operator that is a modified conservative semi-Lagrangian discretization and generalizes those we have developed previously for nonconservative scalar linear hyperbolic problems. Numerical tests are performed for the inviscid Burgers and Buckley–Leverett equations. For many test problems, the solver converges in just a handful of iterations with a convergence rate independent of mesh resolution, including problems with (interacting) shocks and rarefactions.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1987

Hamiltonian theory of nonlinear waves in planetary rings

The derivation of a Hamiltonian field theory for nonlinear density waves in Saturn's rings is discussed. Starting with a Hamiltonian for a discrete system of gravitating streamlines, an averaged Hamiltonian is obtained by successive applications of Lie transforms. The transformation may be carried out to any desired order in q, where q is the nonlinearity parameter defined in the work of Shu, et al (1985) and Borderies et al (1985). Subsequent application of the Wentzel-Kramer-Brillouin Method approximation yields an asymptotic field Hamiltonian. Both the nonlinear dispersion relation and the wave action transport equation are easily derived from the corresponding Lagrangian by the standard variational principle.

Stewart, G. R.↗

NASA NTRS1990

The pdf approach to turbulent flow

This paper provides a detailed discussion of the theory and application of probability density function (pdf) methods, which provide a complete statistical description of turbulent flow fields at a single point or a finite number of points. The basic laws governing the flow of Newtonian fluids are set up in the Eulerian and the Lagrangian frame, and the exact and linear equations for the characteristic functionals in those frames are discussed. Pdf equations in both frames are derived as Fourier transforms of the equations of the characteristic functions. Possible formulations for the nonclosed terms in the pdf equation are discussed, their properties are assessed, and closure modes for the molecular-transport and the fluctuating pressure-gradient terms are reviewed. The application of pdf methods to turbulent combustion flows, supersonic flows, and the interaction of turbulence with shock waves is discussed.

Kollmann, W.↗