Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Lagrangian function”

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 73 records · Page 4

Flows in the space of interacting chiral boson theories

We study interacting theories of N left-moving and N ¯ right-moving Floreanini-Jackiw bosons in two dimensions. A parametrized family of such theories is shown to enjoy (nonmanifest) Lorentz invariance if and only if its Lagrangian obeys a flow equation driven by a function of the energy-momentum tensor. We discuss the canonical quantization of such theories along classical stress tensor flows, focusing on the case of the root- T T ¯ deformation, where we obtain perturbative results for the deformed spectrum in a certain large-momentum limit. In the special case N = N ¯ , we consider the quantum effective action for the root- T T ¯ -deformed theory by expanding around a general classical background, and we find that the one-loop contribution vanishes for backgrounds with constant scalar gradients. Our analysis can also be interpreted via dual U ( 1 ) Chern-Simons theories in three dimensions, which might be used to describe deformations of charged AdS 3 black holes or quantum Hall systems. Published by the American Physical Society 2024

Ebert, Stephen (ORCID:0000000154670155)

The restricted three-body problem including radiation pressure

The force of radiation pressure is included in the restricted three-body problem by considering the major radial components of the pressure force for the case of a particle in the vicinity of two luminous massive bodies, as well as by introducing the Poynting-Robertson effect for the case of one luminous body. The positions of the Lagrangian points L4 and L5 are found as functions of the ratio of radiation to gravitational forces. The Poynting-Robertson effect renders the L4 and L5 points unstable on a time scale that is long compared to the period of rotation of the two massive bodies. Implications for space colonization and a mechanism for producing azimuthal asymmetries in the interplanetary dust complex are discussed.

Schuerman, D. W.

Numerical solution of periodic vortical flows about a thin airfoil

A numerical method is developed for computing periodic, three-dimensional, vortical flows around isolated airfoils. The unsteady velocity is split into a vortical component which is a known function of the upstream flow conditions and the Lagrangian coordinates of the mean flow, and an irrotational field whose potential satisfies a nonconstant-coefficient, inhomogeneous, convective wave equation. Solutions for thin airfoils at zero degrees incidence to the mean flow are presented in this paper. Using an elliptic coordinate transformation, the computational domain is transformed into a rectangle. The Sommerfeld radiation condition is applied to the unsteady pressure on the grid line corresponding to the far field boundary. The results are compared with a Possio solver, and it is shown that for maximum accuracy the grid should depend on both the Mach number and reduced frequency. Finally, in order to assess the range of validity of the classical thin airfoil approximation, results for airfoils with zero thickness are compared with results for airfoils with small thickness.

Scott, James R.

Numerical solution of periodic vortical flows about a thin airfoil

A numerical method is developed for computing periodic, three-dimensional, vortical flows around isolated airfoils. The unsteady velocity is split into a vortical component which is a known function of the upstream flow conditions and the Lagrangian coordinates of the mean flow, and an irrotational field whose potential satisfies a nonconstant-coefficient, inhomogeneous, convective wave equation. Solutions for thin airfoils at zero degrees incidence to the mean flow are presented in this paper. Using an elliptic coordinate transformation, the computational domain is transformed into a rectangle. The Sommerfield radiation condition is applied to the unsteady pressure on the grid line corresponding to the far field boundary. The results are compared with a Possion solver, and it is shown that for maximum accuracy the grid should depend on both the Mach number and reduced frequency. Finally, in order to assess the range of validity of the classical thin airfoil approximation, results for airfoils with zero thickness are compared with results for airfoils with small thickness.

Scott, James R.

A note on a Lagrangian formulation for motion about the collinear points

A Lagrangian formulation for the three-dimensional motion of a satellite in the vicinity of the collinear points of the circular-restricted problem is reconsidered. It is shown that the influence of the primaries can be expressed in the form of two third-body disturbing functions. By use of this approach, the equations for the Lagrangian and for the motion itself are readily developed into highly compact expressions. All orders of the non-linear developments are shown to be easily obtainable using well-known recursive relationships. The resulting forms for these equations are well suited for use in the initial phase of canonical or non-canonical investigations.

Richardson, D. L.

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.

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.

Truss topology optimization with simultaneous analysis and design

Strategies for topology optimization of trusses for minimum weight subject to stress and displacement constraints by Simultaneous Analysis and Design (SAND) are considered. The ground structure approach is used. A penalty function formulation of SAND is compared with an augmented Lagrangian formulation. The efficiency of SAND in handling combinations of general constraints is tested. A strategy for obtaining an optimal topology by minimizing the compliance of the truss is compared with a direct weight minimization solution to satisfy stress and displacement constraints. It is shown that for some problems, starting from the ground structure and using SAND is better than starting from a minimum compliance topology design and optimizing only the cross sections for minimum weight under stress and displacement constraints. A member elimination strategy to save CPU time is discussed.

Sankaranarayanan, S.

Real-time trajectory optimization on parallel processors

A parallel algorithm has been developed for rapidly solving trajectory optimization problems. The goal of the work has been to develop an algorithm that is suitable to do real-time, on-line optimal guidance through repeated solution of a trajectory optimization problem. The algorithm has been developed on an INTEL iPSC/860 message passing parallel processor. It uses a zero-order-hold discretization of a continuous-time problem and solves the resulting nonlinear programming problem using a custom-designed augmented Lagrangian nonlinear programming algorithm. The algorithm achieves parallelism of function, derivative, and search direction calculations through the principle of domain decomposition applied along the time axis. It has been encoded and tested on 3 example problems, the Goddard problem, the acceleration-limited, planar minimum-time to the origin problem, and a National Aerospace Plane minimum-fuel ascent guidance problem. Execution times as fast as 118 sec of wall clock time have been achieved for a 128-stage Goddard problem solved on 32 processors. A 32-stage minimum-time problem has been solved in 151 sec on 32 processors. A 32-stage National Aerospace Plane problem required 2 hours when solved on 32 processors. A speed-up factor of 7.2 has been achieved by using 32-nodes instead of 1-node to solve a 64-stage Goddard problem.

Psiaki, Mark L.

A vortex wake capturing method for potential flow calculations

A method is presented for modifying finite difference solutions of the potential equation to include the calculation of non-planar vortex wake features. The approach is an adaptation of Baker's 'cloud in cell' algorithm developed for the stream function-vorticity equations. The vortex wake is tracked in a Lagrangian frame of reference as a group of discrete vortex filaments. These are distributed to the Eulerian mesh system on which the velocity is calculated by a finite difference solution of the potential equation. An artificial viscosity introduced by the finite difference equations removes the singular nature of the vortex filaments. Computed examples are given for the two-dimensional time dependent roll-up of vortex wakes generated by wings with different spanwise loading distributions.

Murman, E. M.

Tropical Cyclogenesis in a Tropical Wave Critical Layer: Easterly Waves

The development of tropical depressions within tropical waves over the Atlantic and eastern Pacific is usually preceded by a "surface low along the wave" as if to suggest a hybrid wave-vortex structure in which flow streamlines not only undulate with the waves, but form a closed circulation in the lower troposphere surrounding the low. This structure, equatorward of the easterly jet axis, is identified herein as the familiar critical layer of waves in shear flow, a flow configuration which arguably provides the simplest conceptual framework for tropical cyclogenesis resulting from tropical waves, their interaction with the mean flow, and with diabatic processes associated with deep moist convection. The recirculating Kelvin cat's eye within the critical layer represents a sweet spot for tropical cyclogenesis in which a proto-vortex may form and grow within its parent wave. A common location for storm development is given by the intersection of the wave's critical latitude and trough axis at the center of the cat's eye, with analyzed vorticity centroid nearby. The wave and vortex live together for a time, and initially propagate at approximately the same speed. In most cases this coupled propagation continues for a few days after a tropical depression is identified. For easterly waves, as the name suggests, the propagation is westward. It is shown that in order to visualize optimally the associated Lagrangian motions, one should view the flow streamlines, or stream function, in a frame of reference translating horizontally with the phase propagation of the parent wave. In this co-moving frame, streamlines are approximately equivalent to particle trajectories. The closed circulation is quasi-stationary, and a dividing streamline separates air within the cat's eye from air outside.

Dunkerton, T. J.

Analysis of liquid rocket combustion chamber turbulence levels from diffusion data

The intensity of turbulence and the Lagrangian correlation coefficient in a liquid-rocket combustion chamber have been analytically determined from an analysis of experimental diffusion data obtained in a small rocket engine which operated at 300-psia chamber pressure and produced approximately 250 pounds thrust. Results of gas-sample measurements obtained by Orsat and gas-chromatograph techniques to determine helium-concentration profiles were analyzed on the basis of Taylor's (1921) turbulent diffusion theory to obtain turbulence flow-field parameters. The results of the analysis indicate that turbulent diffusion in a combustion chamber can be adequately modeled by the one-dimensional Taylor theory, which assumes that the intensity of turbulence is a function only of axial distance in the chamber and that the Lagrangian correlation coefficient is expressed by a power law. The results indicate a higher intensity of turbulence and lower correlation than previously expected.

Smith, L. O., Jr.

Design optimization of space structures

The topology-shape-size optimization of space structures is investigated through Kikuchi's homogenization method. The method starts from a 'design domain block,' which is a region of space into which the structure is to materialize. This domain is initially filled with a finite element mesh, typically regular. Force and displacement boundary conditions corresponding to applied loads and supports are applied at specific points in the domain. An optimal structure is to be 'carved out' of the design under two conditions: (1) a cost function is to be minimized, and (2) equality or inequality constraints are to be satisfied. The 'carving' process is accomplished by letting microstructure holes develop and grow in elements during the optimization process. These holes have a rectangular shape in two dimensions and a cubical shape in three dimensions, and may also rotate with respect to the reference axes. The properties of the perforated element are obtained through an homogenization procedure. Once a hole reaches the volume of the element, that element effectively disappears. The project has two phases. In the first phase the method was implemented as the combination of two computer programs: a finite element module, and an optimization driver. In the second part, focus is on the application of this technique to planetary structures. The finite element part of the method was programmed for the two-dimensional case using four-node quadrilateral elements to cover the design domain. An element homogenization technique different from that of Kikuchi and coworkers was implemented. The optimization driver is based on an augmented Lagrangian optimizer, with the volume constraint treated as a Courant penalty function. The optimizer has to be especially tuned to this type of optimization because the number of design variables can reach into the thousands. The driver is presently under development.

Felippa, Carlos

Distributed Control by Lagrangian Steepest Descent

Often adaptive, distributed control can be viewed as an iterated game between independent players. The coupling between the players mixed strategies, arising as the system evolves from one instant to the next, is determined by the system designer. Information theory tells us that the most likely joint strategy of the players, given a value of the expectation of the overall control objective function, is the minimizer of a function o the joint strategy. So the goal of the system designer is to speed evolution of the joint strategy to that Lagrangian mhimbhgpoint,lowerthe expectated value of the control objective function, and repeat Here we elaborate the theory of algorithms that do this using local descent procedures, and that thereby achieve efficient, adaptive, distributed control.

Wolpert, David H.

An Empirical Quantile Estimation Approach for Chance-Constrained Nonlinear Optimization Problems

We investigate an empirical quantile estimation approach to solve chance-constrained nonlinear optimization problems. Our approach is based on the reformulation of the chance constraint as an equivalent quantile constraint to provide stronger signals on the gradient. In this approach, the value of the quantile function is estimated empirically from samples drawn from the random parameters, and the gradient of the quantile function is estimated via a finite-difference approximation on top of the quantile-function-value estimation. We establish a convergence theory of this approach within the framework of an augmented Lagrangian method for solving general nonlinear constrained optimization problems. The foundation of the convergence analysis is a concentration property of the empirical quantile process, and the analysis is divided based on whether or not the quantile function is differentiable. In contrast to the sampling-and-smoothing approach used in the literature, the method developed in this paper does not involve any smoothing function and hence the quantile-function gradient approximation is easier to implement and there are less accuracy-control parameters to tune. Furthermore, we demonstrate the effectiveness of this approach and compare it with a smoothing method for the quantile-gradient estimation. Numerical investigation shows that the two approaches are competitive for certain problem instances.

Applied Probability

Large-Deformation Displacement Transfer Functions for Shape Predictions of Highly Flexible Slender Aerospace Structures

Large deformation displacement transfer functions were formulated for deformed shape predictions of highly flexible slender structures like aircraft wings. In the formulation, the embedded beam (depth wise cross section of structure along the surface strain sensing line) was first evenly discretized into multiple small domains, with surface strain sensing stations located at the domain junctures. Thus, the surface strain (bending strains) variation within each domain could be expressed with linear of nonlinear function. Such piecewise approach enabled piecewise integrations of the embedded beam curvature equations [classical (Eulerian), physical (Lagrangian), and shifted curvature equations] to yield closed form slope and deflection equations in recursive forms.

Deflection equations

Structural optimization of framed structures using generalized optimality criteria

The application of a generalized optimality criteria to framed structures is presented. The optimality conditions, Lagrangian multipliers, resizing algorithm, and scaling procedures are all represented as a function of the objective and constraint functions along with their respective gradients. The optimization of two plane frames under multiple loading conditions subject to stress, displacement, generalized stiffness, and side constraints is presented. These results are compared to those found by optimizing the frames using a nonlinear mathematical programming technique.

Kolonay, R. M.