Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Calculus of Variation”

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

Generation of NEP heliocentric trajectory data

A study, designed to generate representative nuclear electric propulsion data for rendezvous missions to the comet Encke using the variational calculus program HILTOP, is presented. Other purposes of the study include a comparison of the HILTOP data with equivalent data generated with QUICKTOP program and to propose approaches for storing and subsequently accessing the optimum trajectory and performance data in the QUICKLY program.

Horsewood, J. L.

A study of unmanned mission opportunities to comets and asteroids

Several unmanned multiple-target mission opportunities to comets and asteroids were studied. The targets investigated include Grigg-Skjellerup, Giacobini-Zinner, Tuttle-Giacobini-Kresak, Borrelly, Halley, Schaumasse, Geographos, Eros, Icarus, and Toro, and the trajectories consist of purely ballistic flight, except that powered swingbys and deep space burns are employed when necessary. Optimum solar electric rendezvous trajectories to the comets Giacobini-Zinner/85, Borrelly/87, and Temple (2)/83 and /88 employing the 8.67 kw Sert III spacecraft modified for interplanetary flight were also investigated. The problem of optimizing electric propulsion heliocentric trajectories, including the effects of geocentric launch asymptote declination on launch vehicle performance capability, was formulated, and a solution developed using variational calculus techniques. Improvements were made to the HILTOP trajectory optimization computer program. An error analysis of high-thrust maneuvers involving spin-stabilized spacecraft was developed and applied to a synchronous meteorological satellite mission.

Mann, F. I.

Heliocentric interplanetary low thrust trajectory optimization program, supplement 1, part 2

The improvements made to the HILTOP electric propulsion trajectory computer program are described. A more realistic propulsion system model was implemented in which various thrust subsystem efficiencies and specific impulse are modeled as variable functions of power available to the propulsion system. The number of operating thrusters are staged, and the beam voltage is selected from a set of five (or less) constant voltages, based upon the application of variational calculus. The constant beam voltages may be optimized individually or collectively. The propulsion system logic is activated by a single program input key in such a manner as to preserve the HILTOP logic. An analysis describing these features, a complete description of program input quantities, and sample cases of computer output illustrating the program capabilities are presented.

Mann, F. I.

An algorithm to generate input data from meteorological and space shuttle observations to validate a CH4-CO model

Objective procedures to analyze data from meteorological and space shuttle observations to validate a three dimensional model were investigated. The transport and chemistry of carbon monoxide and methane in the troposphere were studied. Four aspects were examined: (1) detailed evaluation of the variational calculus procedure, with the equation of continuity as a strong constraint, for adjustment of global tropospheric wind fields; (2) reduction of the National Meteorological Center (NMC) data tapes for data input to the OSTA-1/MAPS Experiment; (3) interpolation of the NMC Data for input to the CH4-CO model; and (4) temporal and spatial interpolation procedures of the CO measurements from the OSTA-1/MAPS Experiment to generate usable contours of the data.

Peters, L. K.

A variational dynamic programming approach to robot-path planning with a distance-safety criterion

An approach to robot-path planning is developed by considering both the traveling distance and the safety of the robot. A computationally-efficient algorithm is developed to find a near-optimal path with a weighted distance-safety criterion by using a variational calculus and dynamic programming (VCDP) method. The algorithm is readily applicable to any factory environment by representing the free workspace as channels. A method for deriving these channels is also proposed. Although it is developed mainly for two-dimensional problems, this method can be easily extended to a class of three-dimensional problems. Numerical examples are presented to demonstrate the utility and power of this method.

Suh, Suk-Hwan

Frequency optimization of repetitive lattice beam-like structures using a continuum model

A new method for obtaining the maximum frequency design of a beam-like repetitive lattice structure is presented. Using existing techniques, the lattice is first modeled as an equivalent anisotropic Timoshenko beam. The computation of the stiffness and inertial properties of the beam, determined by matching the strain and kinetic energies of the beam with those of the lattice, is facilitated by the repetitive nature of the lattice. The optimum design is obtained by maximizing Rayleigh's quotient using methods of variational calculus. For the problem selected, results show excellent agreement with those obtained by traditional finite-element methods. Moreover, unlike FE methods, cpu time is relatively unaffected by the size of the truss.

Reiss, Robert

Nonprincipal plane scattering of flat plates and pattern control of horn antennas

Using the geometrical theory of diffraction, the traditional method of high frequency scattering analysis, the prediction of the radar cross section of a perfectly conducting, flat, rectangular plate is limited to principal planes. Part A of this report predicts the radar cross section in nonprincipal planes using the method of equivalent currents. This technique is based on an asymptotic end-point reduction of the surface radiation integrals for an infinite wedge and enables nonprincipal plane prediction. The predicted radar cross sections for both horizontal and vertical polarizations are compared to moment method results and experimental data from Arizona State University's anechoic chamber. In part B, a variational calculus approach to the pattern control of the horn antenna is outlined. The approach starts with the optimization of the aperture field distribution so that the control of the radiation pattern in a range of directions can be realized. A control functional is thus formulated. Next, a spectral analysis method is introduced to solve for the eigenfunctions from the extremal condition of the formulated functional. Solutions to the optimized aperture field distribution are then obtained.

Balanis, Constantine A.

Minimal Joule dissipation models of magnetospheric convection

This paper gives a topical review of theoretical models of magnetospheric convection based on the concept of minimal Joule dissipation. A two-dimensional slab model of the ionosphere featuring an enhanced conductivity auroral oval is used to compute high-latitude electric fields and currents. Mathematical methods used in the modeling include Fourier analysis, fast Fourier transforms, and variational calculus. Also, conformal transformations are introduced in the analysis, which enable the auroral oval to be represented as a nonconcentric, crescent-shaped figure. Convection patterns appropriate to geomagnetic quiet and disturbed conditions are computed, the differentiating variable being the relative amount of power dissipated in the magnetospheric ring current. When ring current dissipation is small, the convection electric field is restricted to high latitudes (shielding regime), and when it is large, a significant penetration of the field to low latitudes occurs, accompanied by an increase in the ratio of the region I current to the region 2 current.

Barbosa, D. D.

Optimum single modal and bimodal buckling design of symmetric laminates

Variational calculus is used to determine the design that maximizes the resistance of classical symmetric laminates against buckling. The orientations of the constituent orthotropic laminae with respect to the principal axes of the laminate are the design variables. It is shown that the optimal design may not be a point of analyticity of the buckling load. Local analytic extrema are obtained from the design derivatives of the buckling load. Nonanalytic extrema occur whenever the buckling load is a repeated eigenvalue. A novel approach, using a directional design derivative, is employed to determine nonanalytic extrema. Specific examples are presented for biaxial buckling for several different boundary conditions.

Qian, B.

Aeroacoustic and aerodynamic applications of the theory of nonequilibrium thermodynamics

Recent developments in the field of nonequilibrium thermodynamics associated with viscous flows are examined and related to developments to the understanding of specific phenomena in aerodynamics and aeroacoustics. A key element of the nonequilibrium theory is the principle of minimum entropy production rate for steady dissipative processes near equilibrium, and variational calculus is used to apply this principle to several examples of viscous flow. A review of nonequilibrium thermodynamics and its role in fluid motion are presented. Several formulations are presented of the local entropy production rate and the local energy dissipation rate, two quantities that are of central importance to the theory. These expressions and the principle of minimum entropy production rate for steady viscous flows are used to identify parallel-wall channel flow and irrotational flow as having minimally dissipative velocity distributions. Features of irrotational, steady, viscous flow near an airfoil, such as the effect of trailing-edge radius on circulation, are also found to be compatible with the minimum principle. Finally, the minimum principle is used to interpret the stability of infinitesimal and finite amplitude disturbances in an initially laminar, parallel shear flow, with results that are consistent with experiment and linearized hydrodynamic stability theory. These results suggest that a thermodynamic approach may be useful in unifying the understanding of many diverse phenomena in aerodynamics and aeroacoustics.

Horne, W. Clifton

SeGRAm - A practical and versatile tool for spacecraft trajectory optimization

An implementation of the Sequential Gradient/Restoration Algorithm, SeGRAm, is presented along with selected examples. This spacecraft trajectory optimization and simulation program uses variational calculus to solve problems of spacecraft flying under the influence of one or more gravitational bodies. It produces a series of feasible solutions to problems involving a wide range of vehicles, environments and optimization functions, until an optimal solution is found. The examples included highlight the various capabilities of the program and emphasize in particular its versatility over a wide spectrum of applications from ascent to interplanetary trajectories.

Rishikof, Brian H.

Sufficient conditions for a local minimum of the Bolza problem with multiple terminal point constraints

Sufficient conditions for a weak relative minimum are derived for a form of the Bolza problem of variational calculus. The derivation ties together first-order and second-order conditions in a unified consistent manner, addresses controllability considerations in some detail, and clarifies some small inconsistencies in earlier work. The resulting second-order conditions of optimality involve the integration of fewer backward-sweep matrix elements than the standard conditions in the literature for problems with an unspecified final time. As a result, the backward-sweep matrices have the same general structure and dimensionality whether the final time is specified or free, with the terminal values of two of three sweep matrices being more complicated in the latter case.

Wood, Lincoln J.

A simple analytical aerodynamic model of Langley Winged-Cone Aerospace Plane concept

A simple three DOF analytical aerodynamic model of the Langley Winged-Coned Aerospace Plane concept is presented in a form suitable for simulation, trajectory optimization, and guidance and control studies. The analytical model is especially suitable for methods based on variational calculus. Analytical expressions are presented for lift, drag, and pitching moment coefficients from subsonic to hypersonic Mach numbers and angles of attack up to +/- 20 deg. This analytical model has break points at Mach numbers of 1.0, 1.4, 4.0, and 6.0. Across these Mach number break points, the lift, drag, and pitching moment coefficients are made continuous but their derivatives are not. There are no break points in angle of attack. The effect of control surface deflection is not considered. The present analytical model compares well with the APAS calculations and wind tunnel test data for most angles of attack and Mach numbers.

Pamadi, Bandu N.

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.

Inner Core Rotation from Geomagnetic Westward Drift and a Stationary Spherical Vortex in Earth's Core

The idea that geomagnetic westward drift indicates convective leveling of the planetary momentum gradient within Earth's core is pursued in search of a differentially rotating mean state, upon which various oscillations and secular effects might be superimposed. The desired state conforms to roughly spherical boundary conditions, minimizes dissipative interference with convective cooling in the bulk of the core, yet may aid core cooling by depositing heat in the uppermost core and lower mantle. The variational calculus of stationary dissipation applied to a spherical vortex within the core yields an interesting differential rotation profile, akin to spherical Couette flow bounded by thin Hartmann layers. Four boundary conditions are required. To concentrate shear induced dissipation near the core-mantle boundary, these are taken to be: (i) no-slip at the core-mantle interface; (ii) geomagnetically estimated bulk westward flow at the base of the core-mantle boundary layer; (iii) no-slip at the inner-outer core interface; and, to describe magnetic locking of the inner core to the deep outer core; (iv) hydrodynamically stress-free at the inner-outer core boundary. By boldly assuming the axial core angular momentum anomaly to be zero, the super-rotation of the inner core relative to the mantle is calculated to be at most 1.5 deg./yr.

Voorhies, Coerte V.

Inner Core Rotation from Geomagnetic Westward Drift and a Stationary Spherical Vortex in Earth's Core

The idea that geomagnetic westward drift indicates convective leveling of the planetary momentum gradient within Earth's core is pursued in search of a differentially rotating mean state, upon which various oscillations and secular effects might be superimposed. The desired state conforms to roughly spherical boundary conditions, minimizes dissipative interference with convective cooling in the bulk of the core, yet may aide core cooling by depositing heat in the uppermost core and lower mantle. The variational calculus of stationary dissipation applied to a spherical vortex within the core yields an interesting differential rotation profile akin to spherical Couette flow bounded by thin Hartmann layers. Four boundary conditions are required. To concentrate shear induced dissipation near the core-mantle boundary, these are taken to be: (i) no-slip at the core-mantle interface; (ii) geomagnetically estimated bulk westward flow at the base of the core-mantle boundary layer; (iii) no-slip at the inner-outer core interface; and, to describe magnetic locking of the inner core to the deep outer core, (iv) hydrodynamically stress-free at the inner-outer core boundary. By boldly assuming the axial core angular momentum anomaly to be zero, the super-rotation of the inner core is calculated to be at most 1.5 degrees per year.

Voorhies, C. V.

An historical survey of computational methods in optimal control.

Review of some of the salient theoretical developments in the specific area of optimal control algorithms. The first algorithms for optimal control were aimed at unconstrained problems and were derived by using first- and second-variation methods of the calculus of variations. These methods have subsequently been recognized as gradient, Newton-Raphson, or Gauss-Newton methods in function space. A much more recent addition to the arsenal of unconstrained optimal control algorithms are several variations of conjugate-gradient methods. At first, constrained optimal control problems could only be solved by exterior penalty function methods. Later algorithms specifically designed for constrained problems have appeared. Among these are methods for solving the unconstrained linear quadratic regulator problem, as well as certain constrained minimum-time and minimum-energy problems. Differential-dynamic programming was developed from dynamic programming considerations. The conditional-gradient method, the gradient-projection method, and a couple of feasible directions methods were obtained as extensions or adaptations of related algorithms for finite-dimensional problems. Finally, the so-called epsilon-methods combine the Ritz method with penalty function techniques.

Polak, E.