Search NASA⌕ Search

SEARCH · Search NASA

Results for “EULER-LAGRANGE EQUATION”

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 37 records · Page 2

Macroscopic Lagrangian description of warm plasmas. I Formulation of the Lagrangian

A macroscopic Lagrangian is derived which includes a pressure tensor, heat conduction, and elastic collisions. Its Euler-Lagrange equations are shown to be the Maxwell equations and the macroscopic force law. The corresponding Hamiltonian is derived, and Hamilton's canonical equations are also demonstrated to lead to the Maxwell equations and the macroscopic force law. The treatment is facilitated by working in a mixture of Eulerian coordinates (for the fields) and Lagrangian coordinates (for the particle motions), and the introduction of a macroscopic potential expressed in terms of the plasma thermal energy and the energy losses by heat conduction.

Peng, Y.-K. M.↗

Minimum dissipation rates in magnetohydrodynamics

Minimum dissipation rate states are explored for a current-carrying channel of magnetofluid, supported by a dc magnetic field and driven by an applied electric field. The minimization is carried out subject to the constraints of constant axial (toroidal) magnetic flux and constant time-averaged rate of supply of magnetic helicity. The solutions of the resulting Euler-Lagrange equations are sensitive to boundary conditions on the current density j. One set of boundary conditions on j leads to the same consequences as Taylor's 'minimum-energy' theory. A different set leads to significantly different consequences, including a departure from the 'force-free' magnetic profile and a toroidal component of current density that does not reverse at the wall when the toroidal magnetic field reverses.

Montgomery, David↗

Energy stability of thermocapillary convection in a model of the float-zone crystal-growth process

Energy stability theory has been applied to a basic state of thermocapillary convection occurring in a cylindrical half-zone of finite length to determine conditions under which the flow will be stable. Because of the finite length of the zone, the basic state must be determined numerically. Instead of obtaining stability criteria by solving the related Euler-Lagrange equations, the variational problem is attacked directly by discretization of the integrals in the energy identity using finite differences. Results of the analysis are values of the Marangoni number below which axisymmetric disturbances to the basic state will decay, for various values of the other parameters governing the problem.

Shen, Y.↗

Structured adaptive grid generation using algebraic methods

The accuracy of the numerical algorithm depends not only on the formal order of approximation but also on the distribution of grid points in the computational domain. Grid adaptation is a procedure which allows optimal grid redistribution as the solution progresses. It offers the prospect of accurate flow field simulations without the use of an excessively timely, computationally expensive, grid. Grid adaptive schemes are divided into two basic categories: differential and algebraic. The differential method is based on a variational approach where a function which contains a measure of grid smoothness, orthogonality and volume variation is minimized by using a variational principle. This approach provided a solid mathematical basis for the adaptive method, but the Euler-Lagrange equations must be solved in addition to the original governing equations. On the other hand, the algebraic method requires much less computational effort, but the grid may not be smooth. The algebraic techniques are based on devising an algorithm where the grid movement is governed by estimates of the local error in the numerical solution. This is achieved by requiring the points in the large error regions to attract other points and points in the low error region to repel other points. The development of a fast, efficient, and robust algebraic adaptive algorithm for structured flow simulation applications is presented. This development is accomplished in a three step process. The first step is to define an adaptive weighting mesh (distribution mesh) on the basis of the equidistribution law applied to the flow field solution. The second, and probably the most crucial step, is to redistribute grid points in the computational domain according to the aforementioned weighting mesh. The third and the last step is to reevaluate the flow property by an appropriate search/interpolate scheme at the new grid locations. The adaptive weighting mesh provides the information on the desired concentration of points to the grid redistribution scheme. The evaluation of the weighting mesh is accomplished by utilizing the weight function representing the solution variation and the equidistribution law. The selection of the weight function plays a key role in grid adaptation. A new weight function utilizing a properly weighted boolean sum of various flowfield characteristics is defined. The redistribution scheme is developed utilizing Non-Uniform Rational B-Splines (NURBS) representation. The application of NURBS representation results in a well distributed smooth grid by maintaining the fidelity of the geometry associated with boundary curves. Several algebraic methods are applied to smooth and/or nearly orthogonalize the grid lines. An elliptic solver is utilized to smooth the grid lines if there are grid crossings. Various computational examples of practical interest are presented to demonstrate the success of these methods.

Yang, Jiann-Cherng↗

A survey of parametrized variational principles and applications to computational mechanics

This survey paper describes recent developments in the area of parametrized variational principles (PVP's) and selected applications to finite-element computational mechanics. A PVP is a variational principle containing free parameters that have no effect on the Euler-Lagrange equations. The theory of single-field PVP's based on gauge functions (also known as null Lagrangians) is a subset of the inverse problem of variational calculus that has limited value. On the other hand, multifield PVP's are more interesting from theoretical and practical standpoints. Following a tutorial introduction, the paper describes the recent construction of multifield PVP's in several areas of elasticity and electromagnetics. It then discusses three applications to finite-element computational mechanics: the derivation of high-performance finite elements, the development of element-level error indicators, and the constructions of finite element templates. The paper concludes with an overview of open research areas.

Felippa, Carlos A.↗

International Space Station Centrifuge Rotor Models A Comparison of the Euler-Lagrange and the Bond Graph Modeling Approach

The assembly and operation of the International Space Station (ISS) require extensive testing and engineering analysis to verify that the Space Station system of systems would work together without any adverse interactions. Since the dynamic behavior of an entire Space Station cannot be tested on earth, math models of the Space Station structures and mechanical systems have to be built and integrated in computer simulations and analysis tools to analyze and predict what will happen in space. The ISS Centrifuge Rotor (CR) is one of many mechanical systems that need to be modeled and analyzed to verify the ISS integrated system performance on-orbit. This study investigates using Bond Graph modeling techniques as quick and simplified ways to generate models of the ISS Centrifuge Rotor. This paper outlines the steps used to generate simple and more complex models of the CR using Bond Graph Computer Aided Modeling Program with Graphical Input (CAMP-G). Comparisons of the Bond Graph CR models with those derived from Euler-Lagrange equations in MATLAB and those developed using multibody dynamic simulation at the National Aeronautics and Space Administration (NASA) Johnson Space Center (JSC) are presented to demonstrate the usefulness of the Bond Graph modeling approach for aeronautics and space applications.

Nguyen, Louis H.↗

Variational calculations for resonance oscillations of inhomogeneous plasmas

The electrostatic resonance properties of an inhomogeneous plasma column are reported by application of the Rayleigh-Ritz method. A description of the rf equation of motion and pressure term that expresses the system of equations in Euler-Lagrange form is presented. The Rayleigh-Ritz procedure is applied to the corresponding Lagrangian to obtain approximate resonance frequencies and eigenfunctions. An appropriate set of trial coordinate functions is defined, which leads to frequency and eigenfunction estimates.

Peng, Y. K. M.↗

Variational calculations for resonance oscillations of inhomogeneous plasmas

The electrostatic resonance properties of an inhomogeneous plasma column are treated by the Rayleigh-Ritz method. In contrast to Parker, Nickel & Gould (1964), who carried out an exact computation, the present treatment uses a description of the RF equation of motion and pressure term that allows one to express the system of equations in Euler-Lagrange form. The Rayleigh-Ritz procedure is then applied to the corresponding Lagrangian, to obtain approximate resonance frequences and eigenfunctions. An appropriate set of trial coordinate functions is defined, which leads to frequency and eigenfunction estimates in excellent agreement with the work of Parker et al. (1964).

Peng, Y.-K. M.↗

Multi-Impulse to Time Optimal Finite Burn Trajectory Conversion

A novel algorithm is presented that provides an improvement over a traditional parameter optimization method when solving time-optimal, finite-burn pseudo-rendezvous spacecraft trajectory problems. A hybrid optimization procedure is described that converts a set of multiple-impulses, representing high- or low-thrust maneuvers, to an exact time-optimal finite-burn trajectory for a thrust limited, constant exhaust velocity spacecraft. The Hybrid Method applies a control law derived from the Euler-Lagrange system of equations within the classical Indirect Method to a modern Direct Method. An iterative adjoint-control transformation and an evolving constraint vector are introduced to solve the optimal control two-point boundary value problem. This method requires no prior knowledge of the solution, which adds simplicity to the trajectory design process and aids automation. Examples are shown for low-thrust apogee raise maneuvers, non-coplanar Earth orbit transfers, and a modified Deep Space 1 low-thrust trajectory. A numerically significant improvement to objective cost is shown across all application problems compared to a traditional solution method, as well as a significant improvement to convergence speed for a select class of problems.

Joshua A Fogel↗

Variational objective analyses for cyclone studies

The basic analysis equations, i.e., the two horizontal momentum equations, the hydrostatic equation, and the integrated continuity equation were derived for the nonlinear vertical coordinate, nondimensionalized, and expressed in finite differences on a staggered grid. Special care was taken to transform the hydrostatic equation and the pressure gradient terms of the horizontal momentum equations to nearly eliminate truncation error over steeply sloping terrain. This formulation also eliminated explicit reference to orographically induced variations in the thermodynamic variables so that the variational adjustments are on the scale of the meteorological perturbations. The analysis equations were subjected to the Euler-Lagrange operations as expressed for finite differences and an additional set of five partial differential equations was derived, bringing to nine the number of equations in Model I. Higher order terms, terms containing observed quantities, and terms containing none of the variables to be adjusted were grouped into forcing functions and the equations were solved for the zero order terms. Zero order variables were eliminated between these equations and there resulted two diagnostic equations which take the form of general linear second order partial differential equations with nonconstant coefficients.

Achtemeier, G. L.↗

Green's formula and variational principles for cosmic-ray transport with application to rotating and shearing flows

Green's theorem and Green's formula for the diffusive cosmic-ray transport equation in relativistic flows are derived. Green's formula gives the solution of the transport equation in terms of the Green's function of the adjoint transport equation, and in terms of distributed sources throughout the region R of interest, plus terms involving the particle intensity and streaming on the boundary. The adjoint transport equation describes the time-reversed particle transport. An Euler-Lagrange variational principle is then obtained for both the mean scattering frame distribution function f, and its adjoint f(dagger). Variations of the variational functional with respect to f(dagger) yield the transport equation, whereas variations of f yield the adjoint transport equation. The variational principle, when combined with Noether's theorem, yields the conservation law associated with Green's theorem. An investigation of the transport equation for steady, azimuthal, rotating flows suggests the introduction of a new independent variable H to replace the comoving frame momentum variable p'. For the case of rigid rotating flows, H is conserved and is shown to be analogous to the Hamiltonian for a bead on a rigidly rotating wire. The variable H corresponds to a balance between the centrifugal force and the particle inertia in the rotating frame. The physical interpretation of H includes a discussion of nonrelativistic and special relativistic rotating flows as well as the cases of aziuthal, differentially rotating flows about Schwarzs-child and Kerr black holes. Green's formula is then applied to the problem of the acceleration of ultra-high-energy cosmic rays by galactic rotation. The model for galactic rotation assumes an angular velocity law Omega = Omega(sub 0)(omega(sub 0)/omega), where omega denotes radial distance from the axis of rotation. Green's functions for the galactic rotation problem are used to investigate the spectrum of accelerated particles arising from monoenergetic and truncated power-law sources. We conclude that it is possible to accelerate particles beyond the knee by galactic rotation, but not in sufficient number to adequately explain the observed spectrum.

Webb, G. M.↗

Radiation from advanced solid rocket motor plumes

The overall objective of this study was to develop an understanding of solid rocket motor (SRM) plumes in sufficient detail to accurately explain the majority of plume radiation test data. Improved flowfield and radiation analysis codes were developed to accurately and efficiently account for all the factors which effect radiation heating from rocket plumes. These codes were verified by comparing predicted plume behavior with measured NASA/MSFC ASRM test data. Upon conducting a thorough review of the current state-of-the-art of SRM plume flowfield and radiation prediction methodology and the pertinent data base, the following analyses were developed for future design use. The NOZZRAD code was developed for preliminary base heating design and Al2O3 particle optical property data evaluation using a generalized two-flux solution to the radiative transfer equation. The IDARAD code was developed for rapid evaluation of plume radiation effects using the spherical harmonics method of differential approximation to the radiative transfer equation. The FDNS CFD code with fully coupled Euler-Lagrange particle tracking was validated by comparison to predictions made with the industry standard RAMP code for SRM nozzle flowfield analysis. The FDNS code provides the ability to analyze not only rocket nozzle flow, but also axisymmetric and three-dimensional plume flowfields with state-of-the-art CFD methodology. Procedures for conducting meaningful thermo-vision camera studies were developed.

Farmer, Richard C.↗

Spatially homogeneous world models

Einstein field equations for spatially homogeneous spaces, considering rotating matter model in regularized Euler-Lagrange form

Ozsvath, I.↗

Model-Reference Adaptive Control of Distributed Lagrangian Infinite-Dimensional Systems Using Hamilton’s Principle

This paper presents a Hamilton's principle for distributed control of infinite-dimensional systems modeled by a distributed form of the Euler-Lagrange method. The distributed systems are governed by a system of linear partial differential equations in space and time. A generalized potential energy expression is developed that can capture most physical systems including those systems that have no spatial distribution. The Hamilton's principle is applied to derive distributed feedback control methods without resorting to the standard weak-form discretization approach to convert an infinite-dimensional systems to a finite-dimensional systems. It can be shown by the principle of least action that the distributed control synthesized by the Hamilton's principle is a minimum-norm control. A model-reference adaptive control framework is developed for distributed Lagrangian systems in the presence of uncertainty. The theory is demonstrated by an application of adaptive flutter suppression control of a flexible aircraft wing.

Nguyen, Nhan T.↗

Energy Stability of Thermocapillary Convection in Models of the Float Zone Process

The energy-stability of thermocapillary convection in models of the float-zone, crystal-growing process was studied. Stability limits, as functions of pertinent parameters, that will identify conditions which will not allow the existence of an undesirable oscillatory flow instability were determined. Such instabilities may occur in the space processing of semiconductor materials. The determination of the stability limits will involve two sets of numerical computations: (1) solution of the nonlinear governing equations together with the appropriate boundary conditions to determine the basic state (in general, velocity, pressure and temperature fields and the displacement of free surfaces and interfaces); and (2) solution of a nonlinear Euler-Lagrange systems for the energy-stability limit. Both computations, while difficult, should be within the scope of available computer capability and available concepts in numerical analysis. Finite-element methods are attractive candidates for the numerical work.

Neitzel, G. P.↗