Search NASA⌕ Search

SEARCH · Search NASA

Results for “trajectory simulation”

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 199 records · Page 11

Inviscid and Viscous CFD Analysis of Booster Separation for the Space Launch System Vehicle

This paper presents details of Computational Fluid Dynamic (CFD) simulations of the Space Launch System during solid-rocket booster separation using the Cart3D inviscid and Overflow viscous CFD codes. The discussion addresses the use of multiple data sources of computational aerodynamics, experimental aerodynamics, and trajectory simulations for this critical phase of flight. Comparisons are shown between Cart3D simulations and a wind tunnel test performed at NASA Langley Research Center's Unitary Plan Wind Tunnel, and further comparisons are shown between Cart3D and viscous Overflow solutions for the flight vehicle. The Space Launch System (SLS) is a new exploration-class launch vehicle currently in development that includes two Solid Rocket Boosters (SRBs) modified from Space Shuttle hardware. These SRBs must separate from the SLS core during a phase of flight where aerodynamic loads are nontrivial. The main challenges for creating a separation aerodynamic database are the large number of independent variables (including orientation of the core, relative position and orientation of the boosters, and rocket thrust levels) and the complex flow caused by exhaust plumes of the booster separation motors (BSMs), which are small rockets designed to push the boosters away from the core by firing partially in the direction opposite to the motion of the vehicle.

Viscous↗

The Development and Evaluation of an Operational Aerobraking Strategy for the Mars 2001 Odyssey Orbiter

The Mars 2001 Odyssey Orbiter successfully completed the aerobraking phase of its mission on January 11, 2002. This paper discusses the support provided by NASA's Langley Research Center to the navigation team at the Jet Propulsion Laboratory in the planning and operational support of Mars Odyssey Aerobraking. Specifically, the development of a three-degree-of-freedom aerobraking trajectory simulation and its application to pre-flight planning activities as well as operations is described. The importance of running the simulation in a Monte Carlo fashion to capture the effects of mission and atmospheric uncertainties is demonstrated, and the utility of including predictive logic within the simulation that could mimic operational maneuver decision-making is shown. A description is also provided of how the simulation was adapted to support flight operations as both a validation and risk reduction tool and as a means of obtaining a statistical basis for maneuver strategy decisions. This latter application was the first use of Monte Carlo trajectory analysis in an aerobraking mission.

Tartabini, Paul V.↗

The effect of interplanetary trajectory options on a manned Mars aerobrake configuration

Manned Mars missions originating in low Earth orbit (LEO) in the time frame 2010 to 2025 were analyzed to identify preferred mission opportunities and their associated vehicle and trajectory characteristics. Interplanetary and Mars atmospheric trajectory options were examined under the constraints of an initial manned exploration scenario. Two chemically propelled vehicle options were considered: (1) an all propulsive configuration, and (2) a configuration which employs aerobraking at Earth and Mars with low lift/drag (L/D) shapes. Both the interplanetary trajectory options as well as the Mars atmospheric passage are addressed to provide a coupled trajectory simulation. Direct and Venus swingby interplanetary transfers with a 60 day Mars stopover are considered. The range and variation in both Earth and Mars entry velocity are also defined. Two promising mission strategies emerged from the study: (1) a 1.0 to 2.0 year Venus swingby mission, and (2) a 2.0 to 2.5 year direct mission. Through careful trajectory selection, 11 mission opportunities are identified in which the Mars entry velocity is between 6 and 10 km/sec and Earth entry velocity ranges from 11.5 to 12.5 km/sec. Simulation of the Earth return aerobraking maneuver is not performed. It is shown that a low L/D configuration is not feasible for Mars aerobraking without substantial improvements in the interplanetary navigation system. However, even with an advanced navigation system, entry corridor and aerothermal requirements restrict the number of potential mission opportunities. It is also shown that for a large blunt Mars aerobrake configuration, the effects of radiative heating can be significant at entry velocities as low as 6.2 km/sec and will grow to dominate the aerothermal environment at entry velocities above 8.5 km/sec. Despite the additional system complexity associated with an aerobraking vehicle, the use of aerobraking was shown to significantly lower the required initial LEO weight. In comparison with an all propulsive mission, savings between 19 and 59 percent were obtained depending upon launch date.

Braun, Robert D.↗

Post2 End-to-End Descent and Landing Simulation for ALHAT Design Analysis Cycle 2

The ALHAT project is an agency-level program involving NASA centers, academia, and industry, with a primary goal to develop a safe, autonomous, precision-landing system for robotic and crew-piloted lunar and planetary descent vehicles. POST2 is used as the 6DOF descent and landing trajectory simulation for determining integrated system performance of ALHAT landing-system models and lunar environment models. This paper presents updates in the development of the ALHAT POST2 simulation, as well as preliminary system performance analysis for ALDAC-2 used for the testing and assessment of ALHAT system models. The ALDAC-2 POST2 Monte Carlo simulation results have been generated and focus on HRN model performance with the fully integrated system, as well performance improvements of AGNC and TSAR model since the previous design analysis cycle

Davis, Jody L.↗

Mission Analysis Program for Solar Electric Propulsion (MAPSEP). Volume 2: User's manual for earth orbital MAPSEP

The trajectory simulation mode (SIMSEP) requires the namelist SIMSEP to follow TRAJ. The SIMSEP contains parameters which describe the scope of the simulation, expected dynamic errors, and cumulative statistics from previous SIMSEP runs. Following SIMSEP are a set of GUID namelists, one for each guidance correction maneuver. The GUID describes the strategy, knowledge or estimation uncertainties and cumulative statistics for that particular maneuver. The trajectory display mode (REFSEP) requires only the namelist TRAJ followed by scheduling cards, similar to those used in GODSEP. The fixed field schedule cards define: types of data displayed, span of interest, and frequency of printout. For those users who can vary the amount of blank common storage in their runs, a guideline to estimate the total MAPSEP core requirements is given. Blank common length is related directly to the dimension of the dynamic state (NDIM) used in transition matrix (STM) computation, and, the total augmented (knowledge) state (NAUG). The values of program and blank common must be added to compute the total decimal core for a CDC 6500. Other operating systems must scale these requirements appropriately.

Source record↗

Improvement of Automated POST Case Success Rate Using Support Vector Machines

During early conceptual design of complex systems, concept down selection can have a large impact upon program life-cycle cost. Therefore, any concepts selected during early design will inherently commit program costs and affect the overall probability of program success. For this reason it is important to consider as large a design space as possible in order to better inform the down selection process. For conceptual design of launch vehicles, trajectory analysis and optimization often presents the largest obstacle to evaluating large trade spaces. This is due to the sensitivity of the trajectory discipline to changes in all other aspects of the vehicle design. Small deltas in the performance of other subsystems can result in relatively large fluctuations in the ascent trajectory because the solution space is non-linear and multi-modal. In order to help capture large design spaces for new launch vehicles, the authors have performed previous work seeking to automate the execution of the industry standard tool, Program to Optimize Simulated Trajectories (POST). This work initially focused on implementation of analyst heuristics to enable closure of cases in an automated fashion, with the goal of applying the concepts of design of experiments (DOE) and surrogate modeling to enable near instantaneous throughput of vehicle cases.3 As noted in [4] work was then completed to improve the DOE process by utilizing a graph theory based approach to connect similar design points.

Zwack, Mathew R.↗

Improvement of Automated POST Case Success Rate Using Support Vector Machines

During early conceptual design of complex systems, concept down selection can have a large impact upon program life-cycle cost. Therefore, any concepts selected during early design will inherently commit program costs and affect the overall probability of program success. For this reason it is important to consider as large a design space as possible in order to better inform the down selection process. For conceptual design of launch vehicles, trajectory analysis and optimization often presents the largest obstacle to evaluating large trade spaces. This is due to the sensitivity of the trajectory discipline to changes in all other aspects of the vehicle design. Small deltas in the performance of other subsystems can result in relatively large fluctuations in the ascent trajectory because the solution space is non-linear and multi-modal [1]. In order to help capture large design spaces for new launch vehicles, the authors have performed previous work seeking to automate the execution of the industry standard tool, Program to Optimize Simulated Trajectories (POST). This work initially focused on implementation of analyst heuristics to enable closure of cases in an automated fashion, with the goal of applying the concepts of design of experiments (DOE) and surrogate modeling to enable near instantaneous throughput of vehicle cases [2]. Additional work was then completed to improve the DOE process by utilizing a graph theory based approach to connect similar design points [3]. The conclusion of the previous work illustrated the utility of the graph theory approach for completing a DOE through POST. However, this approach was still dependent upon the use of random repetitions to generate seed points for the graph. As noted in [3], only 8% of these random repetitions resulted in converged trajectories. This ultimately affects the ability of the random reps method to confidently approach the global optima for a given vehicle case in a reasonable amount of time. With only an 8% pass rate, tens or hundreds of thousands of reps may be needed to be confident that the best repetition is at least close to the global optima. However, typical design study time constraints require that fewer repetitions be attempted, sometimes resulting in seed points that have only a handful of successful completions. If a small number of successful repetitions are used to generate a seed point, the graph method may inherit some inaccuracies as it chains DOE cases from the non-global-optimal seed points. This creates inherent noise in the graph data, which can limit the accuracy of the resulting surrogate models. For this reason, the goal of this work is to improve the seed point generation method and ultimately the accuracy of the resulting POST surrogate model. The work focuses on increasing the case pass rate for seed point generation.

Zwack, Matthew R.↗

Autolanding Trajectory Design for the X-34

An Autolanding I-load Program (ALIP) is developed to design unpowered autolanding trajectories for the X-34 Mach 8 vehicle. The trajectory is comprised of geometric flight segments that are based on the shuttle approach and landing design (steep glideslope, circular flare, exponential flare to shallow glideslope). Enforcing physical constraints such as loads, vertical descent rate, continuity and smoothness reduces the design problem to a two point boundary value problem with conditions on the initial and final dynamic pressure. Finding a solution required the development of trajectory simulation techniques that constrained the flight profile to a prescribed geometry. The design methodology can be extended beyond the autolanding flight regime by repeating the series of geometric segments and solving multiple two-point boundary values problems (one for each series). The techniques described in this paper facilitate rapid design of reference trajectories.

Barton, Gregg H.↗

Hybrid 3-D rocket trajectory program. Part 1: Formulation and analysis. Part 2: Computer programming and user's instruction

Models utilizing various sub-sets of the six degrees of freedom are used in trajectory simulation. A 3-D model with only linear degrees of freedom is especially attractive, since the coefficients for the angular degrees of freedom are the most difficult to determine and the angular equations are the most time consuming for the computer to evaluate. A computer program is developed that uses three separate subsections to predict trajectories. A launch rail subsection is used until the rocket has left its launcher. The program then switches to a special 3-D section which computes motions in two linear and one angular degrees of freedom. When the rocket trims out, the program switches to the standard, three linear degrees of freedom model.

Huang, L. C. P.↗

DAVINCI: Venus Atmospheric Model Comparisons

The Deep Atmosphere Venus Investigation of Noble gases, Chemistry, and Imaging(DAVINCI) mission aims to answer long-standing questions regarding Venus’ origin using Zephyr, an atmospheric descent probe. Zephyr will be the first probe to take high-resolution aerial photographs of a mountainous tesserae surface as it descends over the Alpha Regio highlands region, which has the oldest surfaces of Venus. The Zephyr’s descent trajectory that determines the touchdown in the Alpha Regio, which is crucial for the DAVINCI mission, depends on Venus’ atmospheric properties and winds. Unfortunately, the atmospheric data for Venus from previous missions is sparse. Therefore, it is essential to consider various atmospheric models and scenarios from past flight data to predict Zephyr’s flight performance, specifically landing ellipse. To this end, this work compares three atmospheric models: the Venus Global Reference Atmospheric Model (Venus-GRAM), the Venus Climate Database(VCD), and an empirical wind model developed by Ralph Lorenz for the DAVINCI trajectory simulation and modeling. This paper compares the mean and variations of different atmospheric properties and winds from these atmospheric models. In addition, this work combines the atmospheric properties and the wind variability from the Venus-GRAM with the winds from the Lorenz-based model to have more stressing Venus wind dispersions that allow for more conservative trajectory analysis. Furthermore, this work relies on the DAVINCI landing ellipse size as a metric to measure how robust the trajectory analysis will be to the change in the atmospheric properties and winds of the Venus atmosphere.

Pardha Sai Chadalavada↗

Simulation reduction using the Taguchi method

A large amount of engineering effort is consumed in conducting experiments to obtain information needed for making design decisions. Efficiency in generating such information is the key to meeting market windows, keeping development and manufacturing costs low, and having high-quality products. The principal focus of this project is to develop and implement applications of Taguchi's quality engineering techniques. In particular, we show how these techniques are applied to reduce the number of experiments for trajectory simulation of the LifeSat space vehicle. Orthogonal arrays are used to study many parameters simultaneously with a minimum of time and resources. Taguchi's signal to noise ratio is being employed to measure quality. A compromise Decision Support Problem and Robust Design are applied to demonstrate how quality is designed into a product in the early stages of designing.

Mistree, Farrokh↗

Development of a Mars Airplane Entry, Descent, and Flight Trajectory

An entry, descent, and flight (EDF) trajectory profile for a Mars airplane mission is defined as consisting of the following elements: ballistic entry of an aeroshell; supersonic deployment of a decelerator parachute; subsonic release of a heat shield; release, unfolding, and orientation of an airplane to flight attitude; and execution of a pull up maneuver to achieve trimmed, horizontal flight. Using the Program to Optimize Simulated Trajectories (POST) a trajectory optimization problem was formulated. Model data representative of a specific Mars airplane configuration, current models of the Mars surface topography and atmosphere, and current estimates of the interplanetary trajectory, were incorporated into the analysis. The goal is to develop an EDF trajectory to maximize the surface-relative altitude of the airplane at the end of a pull up maneuver, while subject to the mission design constraints. The trajectory performance was evaluated for three potential mission sites and was found to be site-sensitive. The trajectory performance, examined for sensitivity to a number of design and constraint variables, was found to be most sensitive to airplane mass, aerodynamic performance characteristics, and the pull up Mach constraint. Based on the results of this sensitivity study, an airplane-drag optimized trajectory was developed that showed a significant performance improvement.

Murray, James E.↗

Comparison of Pre- and Post-Flight Estimations of Pitch Damping Coefficients for Unguided Entry Vehicles

Introduction: Pre-flight predictions of pitch damping coefficient (𝐶!!+𝐶!"̇)curves are used in entry simulations to determine flight readiness and if any dynamic remediation strategies are required. Accurate characterization of these pitch damping coefficient curves is crucial for mission success. Historically, pre-flight estimations of these curves are derived experimentally either via ballistic range testing or with 1-degree-of-freedom (DOF) free-or forced-oscillation wind tunnel testing. Methodology:With post-flight best-estimated trajectories, a fitting procedure as described in Karlgaard, et. al. [1] is used to compare pre-and post-flight pitch damping coefficients. This technique relies on an assumed analytical form of the trajectory, a limited 3-DOF representation described in Schoenenbergerand Queen [2], which results in anEuler-Cauchy solution form.This solution assumes constant density,thus limited amplitude multi-peak windows of the trajectory are considered to populate the pitch damping-amplitudes pace. Fig. 1 shows preliminary estimations of the pitch damping coefficients estimated using the post-flight best estimated trajectoryfor MER B, using a3-amplitude-peak window. Important to note is that the Euler-Cauchy solution also assumes planar motion, thus coning or 𝛽contribution is not included. Further,the estimated post-flight pitch damping approximations and the pre-flight predictions shown are effective representations of pitch damping, or the integrated effect of the pitch damping coefficient over one pitch cycle at the corresponding cycle’s peak amplitude. Effective pitch damping coefficient curves provide more intuitive insights about where the vehicle is stable and permits readily processed comparisons between data-sets. It is important to note that most pitch damping curves, such as those reported in aerodatabases and utilized in flight mechanics trajectory simulations, are represented in instantaneous angle-of-attack space distinct from what is presented in Fig 1. The extracted post-flight pitch damping coefficient points show a wide range of values across a narrow band of amplitudes (Fig.1a), albeit centered on the published pre-flight approximation. This also applies to the prediction as a function of Mach in Fig1b, where the post-flight scatter points are centered around the pre-flight predictions, but the pre-flight predictions notably lack data aboveMach 5.Further tuning of the technique, such as changing the number of peaks considered in each stencil, resolving density in higher resolution, and employing more accurate fitting procedures (beyond scipy.curve_fit()), will be performed. PosterFocus:This study will apply this reconstruction method to the as-flown trajectories for the following unguided entry vehicles: MER A, MER B, IRVE3, Phoenix, Insight, Pathfinder, and SIAD. The pitch damping coefficient will be derived from these trajectories and compared with the pre-flight predictions. Comparing across different missions, entry environments, and vehicle geometries can inform where our pre-flight testing methodologies accurately capture the flight dynamics. Results will provide valuable insight into the accuracy of pre-flight predicted pitch damping curves. Beyond that, this method can shed light on how accurate these curves need to be, which can inform uncertainty ranges applied to pre-flight curves for other missions.

pitch damping coefficient↗

Aerocapture Performance Analysis for a Neptune-Triton Exploration Mission

A systems analysis has been conducted for a Neptune-Triton Exploration Mission in which aerocapture is used to capture a spacecraft at Neptune. Aerocapture uses aerodynamic drag instead of propulsion to decelerate from the interplanetary approach trajectory to a captured orbit during a single pass through the atmosphere. After capture, propulsion is used to move the spacecraft from the initial captured orbit to the desired science orbit. A preliminary assessment identified that a spacecraft with a lift to drag ratio of 0.8 was required for aerocapture. Performance analyses of the 0.8 L/D vehicle were performed using a high fidelity flight simulation within a Monte Carlo executive to determine mission success statistics. The simulation was the Program to Optimize Simulated Trajectories (POST) modified to include Neptune specific atmospheric and planet models, spacecraft aerodynamic characteristics, and interplanetary trajectory models. To these were added autonomous guidance and pseudo flight controller models. The Monte Carlo analyses incorporated approach trajectory delivery errors, aerodynamic characteristics uncertainties, and atmospheric density variations. Monte Carlo analyses were performed for a reference set of uncertainties and sets of uncertainties modified to produce increased and reduced atmospheric variability. For the reference uncertainties, the 0.8 L/D flatbottom ellipsled vehicle achieves 100% successful capture and has a 99.87 probability of attaining the science orbit with a 360 m/s V budget for apoapsis and periapsis adjustment. Monte Carlo analyses were also performed for a guidance system that modulates both bank angle and angle of attack with the reference set of uncertainties. An alpha and bank modulation guidance system reduces the 99.87 percentile DELTA V 173 m/s (48%) to 187 m/s for the reference set of uncertainties.

Starr, Brett R.↗

Combined Experimental and Numerical Simulations of Thermal Barrier Coated Turbine Blades Erosion

A combined experimental and computational study was conducted to investigate the erosion of thermal barrier coated (TBC) blade surfaces by alumina particles ingestion in a single stage turbine. In the experimental investigation, tests of particle surface interactions were performed in specially designed tunnels to determine the erosion rates and particle restitution characteristics under different impact conditions. The experimental results show that the erosion rates increase with increased impingement angle, impact velocity and temperature. In the computational simulations, an Euler-Lagrangian two stage approach is used in obtaining numerical solutions to the three-dimensional compressible Reynolds Averaged Navier-Stokes equations and the particles equations of motion in each blade passage reference frame. User defined functions (UDF) were developed to represent experimentally-based correlations for particle surface interaction models which were employed in the three-dimensional particle trajectory simulations to determine the particle rebound characteristics after each surface impact. The experimentally based erosion UDF model was used to predict the TBC erosion rates on the turbine blade surfaces based on the computed statistical data of the particles impact locations, velocities and angles relative to the blade surface. Computational results are presented for the predicted TBC blade erosion in a single stage commercial APU turbine, for a NASA designed automotive turbine, and for the NASA turbine scaled for modern rotorcraft operating conditions. The erosion patterns in the turbines are discussed for uniform particle ingestion and for particle ingestion concentrated in the inner and outer 5 percent of the stator blade span representing the flow cooling the combustor liner.

Hamed, Awate↗

Comparison of Pre- and Post-Flight Estimations of Pitch Damping Coefficients for Unguided Entry Vehicles

Pre-flight predictions of pitch damping coefficient (C M q + C m ȧ ) curves are used in entry simulations to determine flight readiness and if any dynamic remediation strategies are required. Accurate characterization of these pitch damping coefficient curves is crucial for mission success. Historically, pre-flight estimations of these curves are derived experimentally either via ballistic range testing or with 1-degree-of-freedom (DOF) free-or forced-oscillation wind tunnel testing. With post-flight best-estimated trajectories, a fitting procedure as described in Karlgaard, et. al. is used to compare pre-and post-flight pitch damping coefficients. This technique relies on an assumed analytical form of the trajectory, a limited 3-DOF representation described in Schoenenberger and Queen, which results in an Euler-Cauchy solution form. This solution assumes constant density, thus limited amplitude multi-peak windows of the trajectory are considered to populate the pitch damping-amplitudes pace. The estimated post-flight pitch damping approximations and the pre-flight predictions shown are effective representations of pitch damping, or the integrated effect of the pitch damping coefficient over one pitch cycle at the corresponding cycle’s peak amplitude. Effective pitch damping coefficient curves provide more intuitive insights about where the vehicle is stable and permits readily processed comparisons between data-sets. It is important to note that most pitch damping curves, such as those reported in aero databases and utilized in flight mechanics trajectory simulations, are represented in instantaneous angle-of-attack space distinct from what is presented in the paper.

pitch damping coefficient↗

Guidance and Control System for a Satellite Constellation

A distributed guidance and control algorithm was developed for a constellation of satellites. The system repositions satellites as required, regulates satellites to desired orbits, and prevents collisions. 1. Optimal methods are used to compute nominal transfers from orbit to orbit. 2. Satellites are regulated to maintain the desired orbits once the transfers are complete. 3. A simulator is used to predict potential collisions or near-misses. 4. Each satellite computes perturbations to its controls so as to increase any unacceptable distances of nearest approach to other objects. a. The avoidance problem is recast in a distributed and locally-linear form to arrive at a tractable solution. b. Plant matrix values are approximated via simulation at each time step. c. The Linear Quadratic Gaussian (LQG) method is used to compute perturbations to the controls that will result in increased miss distances. 5. Once all danger is passed, the satellites return to their original orbits, all the while avoiding each other as above. 6. The delta-Vs are reasonable. The controller begins maneuvers as soon as practical to minimize delta-V. 7. Despite the inclusion of trajectory simulations within the control loop, the algorithm is sufficiently fast for available satellite computer hardware. 8. The required measurement accuracies are within the capabilities of modern inertial measurement devices and modern positioning devices.

Bryson, Jonathan Lamar↗

Helicopter flight control system design and evaluation for NOE operations using controller inversion techniques

A flight control system design technique is described that allows rapid evaluation of candidate control laws using computer simulated trajectories. The two-step technique first computes the optimum trajectory for the basic aircraft, and then reconstructs the pilot's control displacements necessary to fly this trajectory, given a particular control system architecture. Since the flight control system itself does not add to the acceleration potential of the helicopter, one may evaluate several proposed control laws in terms of the resulting stick inputs for the same maneuver. The method is illustrated through application to implicit and explicit model-following control laws designed for a battlefield helicopter with and without an auxiliary propulsor. It is shown that the technique is a valuable aid for selection of promising control laws for further simulator studies and flight investigations.

Mckillip, Robert M., Jr.↗