Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational Calculus”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Feedback control for unsteady flow and its application to the stochastic Burgers equation

The study applies mathematical methods of control theory to the problem of control of fluid flow with the long-range objective of developing effective methods for the control of turbulent flows. Model problems are employed through the formalism and language of control theory to present the procedure of how to cast the problem of controlling turbulence into a problem in optimal control theory. Methods of calculus of variations through the adjoint state and gradient algorithms are used to present a suboptimal control and feedback procedure for stationary and time-dependent problems. Two types of controls are investigated: distributed and boundary controls. Several cases of both controls are numerically simulated to investigate the performances of the control algorithm. Most cases considered show significant reductions of the costs to be minimized. The dependence of the control algorithm on the time-descretization method is discussed.

Choi, Haecheon↗
NASA NTRS1951

Torsion and transverse bending of cantilever plates

The problem of combined bending and torsion of cantilever plates of variable thickness, such as might be considered for solid thin high-speed airplane or missile wings, is considered in this paper. The deflections of the plate are assumed to vary linearly across the chord; minimization of the potential energy by means of the calculus of variations then leads to two ordinary linear differential equations for the bending deflections and the twist of the plate. Because the cantilever is analyzed as a plate rather than as a beam, the effect of constraint against axial warping in torsion is inherently included. The application of this method to specific problems involving static deflection, vibration, and buckling of cantilever plates is presented. In the static-deflection problems, taper and sweep are considered.

WINGS, SOLID↗
NASA NTRS1999

Human Mars Mission: Launch Window from Earth Orbit

The determination of orbital window characteristics is of major importance in the analysis of human interplanetary missions and systems. The orbital launch window characteristics are directly involved in the selection of mission trajectories, the development of orbit operational concepts, and the design of orbital launch systems. The orbital launch window problem arises because of the dynamic nature of the relative geometry between outgoing (departure) asymptote of the hyperbolic escape trajectory and the earth parking orbit. The orientation of the escape hyperbola asymptotic relative to the earth is a function of time. The required hyperbola energy level also varies with time. In addition, the inertial orientation of the parking orbit is a function of time because of the perturbations caused by the Earth's oblateness. Thus, a coplanar injection onto the escape hyperbola can be made only at a point in time when the outgoing escape asymptote is contained by the plane of parking orbit. Even though this condition may be planned as a nominal situation, it will not generally represent the more probable injection geometry. The general case of an escape injection maneuver performed at a time other than the coplanar time will involve both a path angle and plane change and, therefore, a delta V penalty. Usually, because of the delta V penalty the actual departure injection window is smaller in duration than that determined by energy requirement alone. This report contains the formulation, characteristics, and test cases for five different launch window modes for Earth orbit. These modes are: 1) One impulsive maneuver from a Highly Elliptical Orbit (HEO); 2) Two impulsive maneuvers from a Highly Elliptical Orbit (HEO); 3) One impulsive maneuver from a Low Earth Orbit (LEO); 4) Two impulsive maneuvers form LEO; and 5) Three impulsive maneuvers form LEO. The formulation of these five different launch window modes provides a rapid means of generating realistic parametric data for space exploration studies. Also the formulation provides vector and geometrical data sufficient for use as a good starting point in detail trajectory analysis based on calculus of variations, steepest descent, or parameter optimization program techniques.

Young, Archie↗
NASA NTRS1997

Preliminary Design of Low-Thrust Interplanetary Missions

For interplanetary missions, highly efficient electric propulsion systems can be used to increase the mass delivered to the destination and/or reduce the trip time over typical chemical propulsion systems. This technology is being demonstrated on the Deep Space 1 mission - part of NASA's New Millennium Program validating technologies which can lower the cost and risk and enhance the performance of future missions. With the successful demonstration on Deep Space 1, future missions can consider electric propulsion as a viable propulsion option. Electric propulsion systems, while highly efficient, produce only a small amount of thrust. As a result, the engines operate during a significant fraction of the trajectory. This characteristic makes it much more difficult to find optimal trajectories. The methods for optimizing low-thrust trajectories are typically categorized as either indirect, or direct. Indirect methods are based on calculus of variations, resulting in a two-point boundary value problem that is solved by satisfying terminal constraints and targeting conditions. These methods are subject to extreme sensitivity to the initial guess of the variables - some of which are not physically intuitive. Adding a gravity assist to the trajectory compounds the sensitivity. Direct methods parameterize the problem and use nonlinear programming techniques to optimize an objective function by adjusting a set of variables. A variety of methods of this type have been examined with varying results. These methods are subject to the limitations of the nonlinear programming techniques. In this paper we present a direct method intended to be used primarily for preliminary design of low-thrust interplanetary trajectories, including those with multiple gravity assists. Preliminary design implies a willingness to accept limited accuracy to achieve an efficient algorithm that executes quickly.

Sims, Jon A.↗
NASA NTRS2000

Guidance Concept for a Mars Ascent Vehicle First Stage

This paper presents a guidance concept for use on the first stage of a Mars Ascent Vehicle (MAV). The guidance is based on a calculus of variations approach similar to that used for the final phase of the Apollo Earth return guidance. A three degree-of-freedom (3DOF) Monte Carlo simulation is used to evaluate performance and robustness of the algorithm.

Queen, Eric M.↗
NASA NTRS2001

Human Exploration Missions Study Launch Window from Earth Orbit

The determination of orbital launch window characteristics is of major importance in the analysis of human interplanetary missions and systems. The orbital launch window characteristics are directly involved in the selection of mission trajectories, the development of orbit operational concepts, and the design of orbital launch systems. The orbital launch window problem arises because of the dynamic nature of the relative geometry between outgoing (departure) asymptote of the hyperbolic escape trajectory and the earth parking orbit. The orientation of the escape hyperbola asymptotic relative to earth is a function of time. The required hyperbola energy level also varies with time. In addition, the inertial orientation of the parking orbit is a function of time because of the perturbations caused by the Earth's oblateness. Thus, a coplanar injection onto the escape hyperbola can be made only at a point in time when the outgoing escape asymptote is contained by the plane of parking orbit. Even though this condition may be planned as a nominal situation, it will not generally represent the more probable injection geometry. The general case of an escape injection maneuver performed at a time other than the coplanar time will involve both a path angle and plane change and, therefore, a Delta(V) penalty. Usually, because of the Delta(V) penalty the actual departure injection window is smaller in duration than that determined by energy requirement alone. This report contains the formulation, characteristics, and test cases for five different launch window modes for Earth orbit. These modes are: (1) One impulsive maneuver from a Low Earth Orbit (LEO), (2) Two impulsive maneuvers from LEO, (3) Three impulsive maneuvers from LEO, (4) One impulsive maneuvers from a Highly Elliptical Orbit (HEO), (5) Two impulsive maneuvers from a Highly Elliptical Orbit (HEO) The formulation of these five different launch window modes provides a rapid means of generating realistic parametric data for space exploration studies. Also the formulation provides vector and geometrical data sufficient for use as a good starting point in detail trajectory analysis based on calculus of variations, steepest descent, or parameter optimization program techniques.

Young, Archie↗
NASA NTRS2002

Trajectories for High Specific Impulse High Specific Power Deep Space Exploration

Flight times and deliverable masses for electric and fusion propulsion systems are difficult to approximate. Numerical integration is required for these continuous thrust systems. Many scientists are not equipped with the tools and expertise to conduct interplanetary and interstellar trajectory analysis for their concepts. Several charts plotting the results of well-known trajectory simulation codes were developed and are contained in this paper. These charts illustrate the dependence of time of flight and payload ratio on jet power, initial mass, specific impulse and specific power. These charts are intended to be a tool by which people in the propulsion community can explore the possibilities of their propulsion system concepts. Trajectories were simulated using the tools VARITOP and IPOST. VARITOP is a well known trajectory optimization code that involves numerical integration based on calculus of variations. IPOST has several methods of trajectory simulation; the one used in this paper is Cowell's method for full integration of the equations of motion. An analytical method derived in the companion paper was also evaluated. The accuracy of this method is discussed in the paper.

Polsgrove, Tara↗
NASA NTRS2002

Simulation of Trajectories for High Specific Impulse Deep Space Exploration

Difficulties in approximating flight times and deliverable masses for continuous thrust propulsion systems have complicated comparison and evaluation of proposed propulsion concepts. These continuous thrust propulsion systems are of interest to many groups, not the least of which are the electric propulsion and fusion communities. Several charts plotting the results of well-known trajectory simulation codes were developed and are contained in this paper. These charts illustrate the dependence of time of flight and payload ratio on jet power, initial mass, specific impulse and specific power. These charts are intended to be a tool by which people in the propulsion community can explore the possibilities of their propulsion system concepts. Trajectories were simulated using the tools VARITOP and IPOST. VARITOP is a well known trajectory optimization code that involves numerical integration based on calculus of variations. IPOST has several methods of trajectory simulation; the one used in this paper is Cowell's method for full integration of the equations of motion. The analytical method derived in the companion paper was also used to simulate the trajectory. The accuracy of this method is discussed in the paper.

Polsgrove, Tara↗
NASA NTRS2002

Atmospheric Ascent Guidance for Rocket-Powered Launch Vehicles

An advanced ascent guidance algorithm for rocket- powered launch vehicles is developed. This algorithm cyclically solves the calculus-of-variations two-point boundary-value problem starting at vertical rise completion through main engine cutoff. This is different from traditional ascent guidance algorithms which operate in a simple open-loop mode until high dynamic pressure (including the critical max-Q) portion of the trajectory is over, at which time guidance operates under the assumption of negligible aerodynamic acceleration (i.e., vacuum dynamics). The initial costate guess is corrected based on errors in the terminal state constraints and the transversality conditions. Judicious approximations are made to reduce the order and complexity of the state/costate system. Results comparing guided launch vehicle trajectories with POST open-loop trajectories are given verifying the basic formulation of the algorithm. Multiple shooting is shown to be a very effective numerical technique for this application. In particular, just one intermediate shooting point, in addition to the initial shooting point, is sufficient to significantly reduce sensitivity to the guessed initial costates. Simulation results from a high-fidelity trajectory simulation are given for the case of launch to sub-orbital cutoff conditions as well as launch to orbit conditions. An abort to downrange landing site formulation of the algorithm is presented.

Dukeman, Greg A.↗
NASA NTRS2002

Primer Vector Optimization: Survey of Theory, New Analysis and Applications

In this paper, a summary of primer vector theory is presented. The applicability of primer vector theory is examined in an effort to understand when and why the theory can fail. For example, since the Calculus of Variations is based on "small" variations, singularities in the linearized (variational) equations of motion along the arcs must be taken into account. These singularities are a recurring problem in analyse that employ small variations. Two examples, the initialization of an orbit and a line of apsides rotation, are presented. Recommendations, future work, and the possible addition of other optimization techniques are also discussed.

Guzman, J. J.↗
NASA NTRS2003

Ultrasail

Ultrasail is a complete sail system for the launch, deployment, stabilization and control of very large solar sails enabling reduced mission times for interplanetary and deep space spacecraft. Ultrasail is an innovative, non-traditional approach to propulsion technology achieved by combining propulsion and control systems developed for formation-flying microsatellites with an innovative solar sail architecture to achieve sq km-class controllable sail areas, sail subsystem area densities of 1 gm per sq m, and thrust levels equivalent to 400 kW ion thruster systems used for comparable deep space missions. Ultrasail can conceivably even achieve outer planetary rendezvous, a deep space capability now reserved for high-mass nuclear and chemical systems. Ultrasail is a Delta IV-launched multi-blade spin-stabilized system with blade lengths as long as 50 km, reminiscent of the MacNeal Heliogyro. The primary innovation is the near-elimination of sail supporting structures by attaching the sail tip to a rigid formation-flying microsatellite truss which deploys the sail blade, and which then articulates the blade to provide attitude control, including spin stabilization and precession of the spin axis. These tip microsatellites are controlled by a solar-powered 3-axis microthruster system (electric or cold gas) to maintain proper sail film tension during deployment and spin-up. The satellite mass also provides a stabilizing centrifugal force on the blade while in rotation. Understanding the dynamics of individual blades is key to the overall dynamics of Ultrasail. Forces and torques that must be modeled include those due to solar pressure, those generated by the microsatellite at the blade tip and by torques applied at the blade root. Centrifugal forces also play a significant role in the deployment and maintenance of the sail configuration. To capture the dynamics of the overall system, the equations of motion for the blades have been derived. Using these differential equations, a control law will be derived to maneuver Ultrasail. This law involves the pitching of the individual blades thereby moving the distribution of the radiation pressure on each individual blade and inducing a resultant torque on the system. The direction of the angular momentum vector and its rate of precession can be controlled through the pitch angle of the blades. The Ultrasail trajectory is also being studied. Optimal or near-optimal trajectories are being generated to showcase Ultrasail performance. Various missions, e.g. outer planet and solar polar missions for observation of the Sun, are currently being investigated to demonstrate the performance enhancements generated by Ultrasail technology. Calculus-of-variations-based optimization software is used to produce optimal Ultrasail trajectories. The performance of these trajectories is being compared to optimal results generated with other propulsion models, including chemical propulsion, ion propulsion, and competing solar sail concepts. Results of these studies will quantify the performance of Ultrasail compared to existing solar sail concepts for high energy missions.

Burton, R.↗
NASA NTRS2004

Angle-of-Attack-Modulated Terminal Point Control for Neptune Aerocapture

An aerocapture guidance algorithm based on a calculus of variations approach is developed, using angle of attack as the primary control variable. Bank angle is used as a secondary control to alleviate angle of attack extremes and to control inclination. The guidance equations are derived in detail. The controller has very small onboard computational requirements and is robust to atmospheric and aerodynamic dispersions. The algorithm is applied to aerocapture at Neptune. Three versions of the controller are considered with varying angle of attack authority. The three versions of the controller are evaluated using Monte Carlo simulations with expected dispersions.

Queen, Eric M.↗
NASA NTRS2014

DUKSUP: A Computer Program for High Thrust Launch Vehicle Trajectory Design and Optimization

From the late 1960's through 1997, the leadership of NASA's Intermediate and Large class unmanned expendable launch vehicle projects resided at the NASA Lewis (now Glenn) Research Center (LeRC). One of LeRC's primary responsibilities --- trajectory design and performance analysis --- was accomplished by an internally-developed analytic three dimensional computer program called DUKSUP. Because of its Calculus of Variations-based optimization routine, this code was generally more capable of finding optimal solutions than its contemporaries. A derivation of optimal control using the Calculus of Variations is summarized including transversality, intermediate, and final conditions. The two point boundary value problem is explained. A brief summary of the code's operation is provided, including iteration via the Newton-Raphson scheme and integration of variational and motion equations via a 4th order Runge-Kutta scheme. Main subroutines are discussed. The history of the LeRC trajectory design efforts in the early 1960's is explained within the context of supporting the Centaur upper stage program. How the code was constructed based on the operation of the Atlas/Centaur launch vehicle, the limits of the computers of that era, the limits of the computer programming languages, and the missions it supported are discussed. The vehicles DUKSUP supported (Atlas/Centaur, Titan/Centaur, and Shuttle/Centaur) are briefly described. The types of missions, including Earth orbital and interplanetary, are described. The roles of flight constraints and their impact on launch operations are detailed (such as jettisoning hardware on heating, Range Safety, ground station tracking, and elliptical parking orbits). The computer main frames on which the code was hosted are described. The applications of the code are detailed, including independent check of contractor analysis, benchmarking, leading edge analysis, and vehicle performance improvement assessments. Several of DUKSUP's many major impacts on launches are discussed including Intelsat, Voyager, Pioneer Venus, HEAO, Galileo, and Cassini.

high thrust trajectory design↗
NASA NTRS2015

DUKSUP: A Computer Program for High Thrust Launch Vehicle Trajectory Design and Optimization

From the late 1960s through 1997, the leadership of NASAs Intermediate and Large class unmanned expendable launch vehicle projects resided at the NASA Lewis (now Glenn) Research Center (LeRC). One of LeRCs primary responsibilities --- trajectory design and performance analysis --- was accomplished by an internally-developed analytic three dimensional computer program called DUKSUP. Because of its Calculus of Variations-based optimization routine, this code was generally more capable of finding optimal solutions than its contemporaries. A derivation of optimal control using the Calculus of Variations is summarized including transversality, intermediate, and final conditions. The two point boundary value problem is explained. A brief summary of the codes operation is provided, including iteration via the Newton-Raphson scheme and integration of variational and motion equations via a 4th order Runge-Kutta scheme. Main subroutines are discussed. The history of the LeRC trajectory design efforts in the early 1960s is explained within the context of supporting the Centaur upper stage program. How the code was constructed based on the operation of the AtlasCentaur launch vehicle, the limits of the computers of that era, the limits of the computer programming languages, and the missions it supported are discussed. The vehicles DUKSUP supported (AtlasCentaur, TitanCentaur, and ShuttleCentaur) are briefly described. The types of missions, including Earth orbital and interplanetary, are described. The roles of flight constraints and their impact on launch operations are detailed (such as jettisoning hardware on heating, Range Safety, ground station tracking, and elliptical parking orbits). The computer main frames on which the code was hosted are described. The applications of the code are detailed, including independent check of contractor analysis, benchmarking, leading edge analysis, and vehicle performance improvement assessments. Several of DUKSUPs many major impacts on launches are discussed including Intelsat, Voyager, Pioneer Venus, HEAO, Galileo, and Cassini.

Launch vehicle performance optimization↗
NASA NTRS2020

Shape Validation and RF Performance of Inflatable Antennas

Inflatable aperture antennas are an emerging technology that is being investigated for potential use in science and exploration missions. In particular, for missions to Mars and beyond, large deployable aperture antennas can provide the antenna gain required for high data rate communications, where the necessary antenna diameter exceeds the available volume of typical launch vehicle platforms. As inflatable aperture antennas have not been proven fully qualified for space missions, the author's Master's Thesis assessed the Ruze equation in characterizing this antenna technology. Inflatable aperture antennas do not follow a parabolic shape, and so the Ruze equation is not applicable due to the macroscopic shape errors of this technology. Therefore, geometric evaluations of the surface profile cannot simply correlate antenna gain degradation with the root-mean-square shape error with a parabolic surface. Consequently, the focus of this work was to derive an accurate mathematical model of an inflatable aperture antenna in order to characterize its Radio Frequency (RF) performance. Calculus of Variations methodologies were used to derive the surface profile shape of the inflatable aperture antenna. Physical Optics techniques were used to generate the antenna pattern profile. Validation testing of the predicted inflatable antenna shape model was performed through use of Laser Radar metrology measurements on an inflatable test article. Assessments of the RF performance of the inflatable aperture antenna, compared with nominally shaped paraboloidal antennas, were obtained through simulations of both technologies using a common diameter, depth, and arc length. Assessments of the RF performance of the inflatable aperture antenna was also performed against itself for changes in distance of the antenna feed location in the axial direction. Whereas the Ruze equation is limited to assessing gain reduction, this effort will also assess beam spreading and first side lobe angle and magnitude. The ability to characterize the RF response of this antenna will provide for an improved understanding of this technology. The accurate representation of the shape of this type of antenna technology will help to identify the most appropriate ways in which this technology could be utilized in planning future communication architectures for NASA missions to Mars and beyond.

Welch, Bryan W.↗
NASA NTRS2020

Shape Validation and RF Performance of Inflatable Antennas

Inflatable aperture antennas are an emerging technology that is being investigated for potential use in science and exploration missions. In particular, for missions to Mars and beyond, large deployable aperture antennas can provide the antenna gain required for high data rate communications, where the necessary antenna diameter exceeds the available volume of typical launch vehicle platforms. As inflatable aperture antennas have not been proven fully qualified for space missions, the author's Master's Thesis assessed the Ruze equation in characterizing this antenna technology. Inflatable aperture antennas do not follow a parabolic shape, and so the Ruze equation is not applicable due to the macroscopic shape errors of this technology. Therefore, geometric evaluations of the surface profile cannot simply correlate antenna gain degradation with the root-mean square shape error with a parabolic surface. Consequently, the focus of this work was to derive an accurate mathematical model of an inflatable aperture antenna to in order to characterize its Radio Frequency (RF) performance. Calculus of Variations methodologies were used to derive the surface profile shape of the inflatable aperture antenna. Physical Optics techniques were used to generate the antenna pattern profile. Validation testing of the predicted inflatable antenna shape model was performed through use of Laser Radar metrology measurements on an inflatable test article. Assessments of the RF performance of the inflatable aperture antenna, compared with nominally shaped paraboloidal antennas, were obtained through simulations of both technologies using a common diameter, depth, and arc length. Assessments of the RF performance of the inflatable aperture antenna was also performed against itself for changes in distance of the antenna feed location in the axial direction. Whereas the Ruze equation is limited to assessing gain reduction, this effort will also assess beam spreading and first side lobe angle and magnitude. The ability to characterize the RF response of this antenna will provide for an improved understanding of this technology. The accurate representation of the shape of this type of antenna technology will help to identify the most appropriate ways in which this technology could be utilized in planning future communication architectures for NASA missions to Mars and beyond.

Welch, Bryan↗
NASA NTRS2019

Investigation of Direct Force Control for Planetary Aerocapture at Neptune

In this work, a direct force control numerical predictor-corrector guidance architecture is developed to enable Neptune aerocapture using flight-heritage blunt body aeroshells. A linear aerodynamics model is formulated for a Mars Science Laboratory-derived aeroshell. The application of calculus of variations shows that the optimal angle of attack and side-slip angle control laws are bang-bang. A closed-loop numerical predictor-corrector direct force control guidance algorithm is developed and numerically simulated using the Program to Optimize Simulated Trajectories II. The Monte Carlo simulated trajectories are demonstrated to be robust to the modeled dispersions in aerodynamics, atmospheric density, and entry state. An aerocapture technology trade study demonstrates that blunt body direct force control aerocapture enables similar performance as slender body bank angle control but halves the peak g-loading.

Deshmukh, Rohan G.↗
NASA NTRS1994

(abstract) Optimal Low Thrust Trajectories Using Differential Inclusion Concepts

Low thrust propulsion systems typically have their greatest benefit for high energy missions or missions with large post-launch maneuver requirements. Missions which have been examined include main belt asteroid rendezvous, comet rendezvous, outer planet and Mercury orbiters, Pluto flyby, and solar probe missions. Low thrust mission design software used to determine these trajectories is based on two distinct formulations of the optimal control problem: the indirect and direct methods. The traditional approach (indirect) is to use the calculus of variations to obtain first order necessary conditions on the states and costates. In contrast, direct methods are conceptually different in that no explicit integration takes place. A direct method based on differential inclusion concepts has been developed and successfully used to compute low thrust trajectories. This new approach removes explicit control dependence from the problem thereby reducing the dimension of the parameter space for the nonlinear programming problem. Also when compared to other direct methods, fewer nonlinear constraints are required to represent the dynamics of the problem.

low↗