Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solve time”

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 1,117 records · Page 62

Evolution of Orion Mission Design for Exploration Mission 1 and 2

The evolving mission design and concepts of NASA's next steps have shaped Orion into the spacecraft that it is today. Since the initial inception of Orion, through the Constellation Program, and now in the Exploration Mission frame-work with the Space Launch System (SLS), each mission design concept and program goal have left Orion with a set of capabilities that can be utilized in many different mission types. Exploration Missions 1 and 2 (EM-1 and EM-2) have now been at the forefront of the mission design focus for the last several years. During that time, different Design Reference Missions (DRMs) were built, analyzed, and modified to solve or mitigate enterprise level design trades to ensure a viable mission from launch to landing. The resulting DRMs for EM-1 and EM-2 were then expanded into multi-year trajectory scans to characterize vehicle performance and Earth-Moon geometry trends. This provides Orion's subsystems with stressing reference trajectories to help design their system. Now that Orion has progressed through the Preliminary and Critical Design Re-views (PDR and CDR) there is a general shift in the focus of mission design from aiding the vehicle design to providing mission specific products needed for pre-flight and real time operations. Some of the mission specific products need-ed include analysis of steering law performance, inputs into navigational accura-cy assessments, abort options at any point in the mission for each valid trajecto-ry in the launch window, recontact avoidance between the upper stage and Orion post nominal separation, etc.

Gutkowski, Jeffrey P.↗

Evolution of Orion Mission Design for Exploration Mission 1 and 2

The evolving mission design and concepts of NASA’s next steps have shaped Orion into the spacecraft that it is today. Since the initial inception of Orion, through the Constellation Program, and now in the Exploration Mission frame-work with the Space Launch System (SLS), each mission design concept and pro-gram goal have left Orion with a set of capabilities that can be utilized in many different mission types. Exploration Missions 1 and 2 (EM-1 and EM-2) have now been at the forefront of the mission design focus for the last several years. During that time, different Design Reference Missions (DRMs) were built, analyzed, and modified to solve or mitigate enterprise level design trades to ensure a viable mission from launch to landing. The resulting DRMs for EM-1 and EM-2 were then expanded into multi-year trajectory scans to characterize vehicle performance as affected by variations in Earth-Moon geometry. This provides Orion’s subsystems with stressing reference trajectories to help design their system. Now that Orion has progressed through the Preliminary and Critical Design Reviews (PDR and CDR), there is a general shift in the focus of mission design from aiding the vehicle design to providing mission specific products needed for pre-flight and real time operations. Some of the mission specific products needed include, large quantities of nominal trajectories for multiple monthly launch periods and abort options at any point in the mission for each valid trajectory in the launch window.

Gutkowski, Jeffrey P.↗

Time dependent computation of the Euler equations for designing 2-D cascades, including the case of transonic shock free design

A numerical method to solve the inverse problem for airfoil cascades is described. The case of compressible, inviscid flow is considered. The proposed method is based on the finite difference approximation of the time-dependent Euler equations. Numerical examples show the capability of the method to design subsonic or transonic shockless blades.

Zannetti, L.↗

Scattering of impact wave by a crack in composite plate

The surface responses due to impact load on an infinite uniaxial graphite/epoxy plate with and without delamination cracks are investigated both in time and frequency domain by using a hybrid method combining the finite element discretization of the near-field with boundary integral representation of the field outside a contour enclosing completely the crack. This combined method leads to a set of linear unsymmetric complex matrix equations which are solved to obtain the response in the frequency domain by biconjugate gradient method. The time domain response is then obtained by using an FFT. In order to capture the time-domain characteristics accurately, high order finite elements have been used. Also, both the six node singular elements and eight node transition elements are used around the crack tips to model the crack-tip singularity. It is shown that from the numerical results for surface responses both depth and length of this crack can be identified.

Ju, T. H.↗

Optimal K/S impulsive transfer

The necessary equations in K/S variables (E-formulation) are developed for an optimal N-impulse solution of the time-open and rendezvous problems. The equations are particularly simple when only two impulses are assumed. The developed algorithm is extremely fast and reliable because the equations that must be solved are smooth, well-behaved (almost sinusoidal) functions. Numerical results for a two-impulse time-open and rendezvous problem are presented.

Jezewski, D. J.↗

Measuring Greenland Ice Mass Variation With Gravity Recovery and the Climate Experiment Gravity and GPS

The response of the Greenland ice sheet to climate change could significantly alter sea level. The ice sheet was much thicker at the last glacial maximum. To gain insight into the global change process and the future trend, it is important to evaluate the ice mass variation as a function of time and space. The Gravity Recovery and Climate Experiment (GRACE) mission to fly in 2001 for 5 years will measure gravity changes associated with the current ice variation and the solid earth's response to past variations. Our objective is to assess the separability of different change sources, accuracy and resolution in the mass variation determination by the new gravity data and possible Global Positioning System (GPS) bedrock uplift measurements. We use a reference parameter state that follows a dynamic ice model for current mass variation and a variant of the Tushingham and Peltier ICE-3G deglaciation model for historical deglaciation. The current linear trend is also assumed to have started 5 kyr ago. The Earth model is fixed as preliminary reference Earth model (PREM) with four viscoelastic layers. A discrete Bayesian inverse algorithm is developed employing an isotropic Gaussian a priori covariance function over the ice sheet and time. We use data noise predicted by the University of Texas and JPL for major GRACE error sources. A 2 mm/yr uplift uncertainty is assumed for GPS occupation time of 5 years. We then carry out covariance analysis and inverse simulation using GRACE geoid coefficients up to degree 180 in conjunction with a number of GPS uplift rates. Present-day ice mass variation and historical deglaciation are solved simultaneously over 146 grids of roughly 110 km x 110 km and with 6 time increments of 3 kyr each, along with a common starting epoch of the current trend. For present-day ice thickness change, the covariance analysis using GRACE geoid data alone results in a root mean square (RMS) posterior root variance of 2.6 cm/yr, with fairly large a priori uncertainties in the parameters and a Gaussian correlation length of 350 km. Simulated inverse can successfully recover most features in the reference present-day change. The RMS difference between them over the grids is 2.8 cm/yr. The RMS difference becomes 1.1 cm/yr when both are averaged with a half Gaussian wavelength of 150 km. With a fixed Earth model, GRACE alone can separate the geoid signals due to past and current load fairly well. Shown are the reference geoid signatures of direct and elastic effects of the current trend, the viscoelastic effect of the same trend starting from 5 kyr ago, the Post Glacial Rebound (PGR), and the predicted GRACE geoid error. The difference between the reference and inverse modeled total viscoelastic signatures is also shown. Although past and current ice mass variations are allowed the same spatial scale, their geoid signals have different spatial patterns. GPS data can contribute to the ice mass determination as well. Additional information is contained in the original.

Wu, Xiao-Ping↗

Characteristics of the boundary-layer equations of the minimum time-to-climb problem

In many singular perturbation solutions of optimal control problems, the most difficult numerical task is to solve the boundary-layer equations. However, these equations have a special structure that may often be used to expedite their solution. This paper begins by noting the general nature of the boundary-layer equations for optimal control problems. These results are then applied to the aircraft minimum time-to-climb problem. A specific numerical example is considered to illustrate the characteristics of the solution of the boundary-layer equations for this problem.

Ardema, M. D.↗

Aerodynamic Characteristics of an Aerospace Vehicle During a Subsonic Pitch-Over Maneuver

Time-dependent CFD has been used to predict aerospace vehicle aerodynamics during a subsonic rotation maneuver. The inviscid 3D3U code is employed to solve the 3-D unsteady flow field using an unstructured grid of tetrahedra. As this application represents a challenge to time-dependent CFD, observations concerning spatial and temporal resolution are included. It is shown that even for a benign rotation rate, unsteady aerodynamic effects are significant during the maneuver. Possibly more significant, however, the rotation maneuver creates ow asymmetries leading to yawing moment, rolling moment, and side force which are not present in the quasi-steady case. A series of steady solutions at discrete points in the maneuver are also computed for comparison with wind tunnel measurements and as a means of quantifying unsteady effects.

Kleb, William L.↗

An efficient finite element technique for sound propagation in axisymmetric hard wall ducts carrying high subsonic Mach number flows

The general non-linear three-dimensional equation for acoustic potential is derived by using a perturbation technique. The linearized axisymmetric equation is then solved by using a finite element algorithm based on the Galerkin formulation for a harmonic time dependence. The solution is carried out in complex number notation for the acoustic velocity potential. Linear, isoparametric, quadrilateral elements with non-uniform distribution across the duct section are implemented. The resultant global matrix is stored in banded form and solved by using a modified Gauss elimination technique. Sound pressure levels and acoustic velocities are calculated from post element solutions. Different duct geometries are analyzed and compared with experimental results.

Tag, I. A.↗

Numerical solutions of the complete Navier-Stokes equations

The objective of this study is to compare the use of assumed pdf (probability density function) approaches for modeling supersonic turbulent reacting flowfields with the more elaborate approach where the pdf evolution equation is solved. Assumed pdf approaches for averaging the chemical source terms require modest increases in CPU time typically of the order of 20 percent above treating the source terms as 'laminar.' However, it is difficult to assume a form for these pdf's a priori that correctly mimics the behavior of the actual pdf governing the flow. Solving the evolution equation for the pdf is a theoretically sound approach, but because of the large dimensionality of this function, its solution requires a Monte Carlo method which is computationally expensive and slow to coverage. Preliminary results show both pdf approaches to yield similar solutions for the mean flow variables.

Hassan, H. A.↗

Computer simulation for time optimal or energy optimal attitude control of spin-stabilized spacecraft.

An original technique for determining the optimal magnetic torque strategy for control of the attitude of spin stabilized spacecraft is presented. By employing Lagrange multipliers and the Calculus of Variations, optimal control equations are derived which define minimum time and minimum energy attitude maneuvers. Computer program algorithms to numerically solve these optimal control equations are also described. The performance of this technique is compared with a commonly employed planning method.

Woolley, R. D.↗

An efficiency study on obtaining the minimum weight of a thermal protection system

Three minimizing techniques are evaluated to determine the most efficient method for minimizing the weight of a thermal protection system and for reducing computer usage time. The methods used (numerical optimization and nonlinear least squares) for solving the minimum-weight problem involving more than one material and more than one constraint are discussed. In addition, the one material and one constraint problem is discussed.

Williams, S. D.↗

Numerical solutions of supersonic and hypersonic laminar flows over a two-dimensional compression corner

An efficient time-splitting, second-order accurate, numerical scheme is used to solve the complete Navier-Stokes equations for supersonic and hypersonic laminar flow over a two-dimensional compression corner. A fine, exponentially stretched mesh spacing is used in the region near the wall for resolving the viscous layer. Good agreement is obtained between the present computed results and experimental measurement for a Mach number of 14.1, a Reynolds number of 104,000, and wedge angles of 15, 18, and 24 deg. The details of the pressure variation across the boundary layer are given, and a correlation between the leading edge shock and the peaks in surface pressure and heat transfer is observed.

Hung, C. M.↗

A coupled fluid-structure dynamic analysis for water impact loads

In this paper a nonlinear transient, hydroelastic method is developed for response analysis of the Space Shuttle recoverable booster on water impact. At each time step, coupled equations of fluid flow and structure dynamics are solved with an iterative, self-correcting process. Rigid body motions are represented by large angle equations. Local deformations are represented by a NASTRAN-generated model of the three-dimensional structure. The fluid is represented by series solutions of the potential flow equation which include the effects of the local structure motions. The resulting water loads and structural stresses obtained will aid in qualification of the final structural design. Comparisons are made to experimental impact data to validate the method.

Herting, D. N.↗

Predicting surface heat flux

Report presents technique involving single embedded thermocouple used to predict flux and temperature for high-or-low conductivity materials that have temperature-and-pressure dependent properties. Technique solves for heat rate and temperature at given surface at each time step, rather than for entire history. Newton-Raphson technique is used for temperature solutions; then quadratic fit is employed. Report contains representative graphs and tables.

Curry, D. M.↗

Some calculations of transonic potential flow for the NACA 64A006 airfoil with oscillating flap

A method for calculating the transonic flow over steady and oscillating airfoils was developed by Isogai. It solves the full potential equation with a semi-implicit, time-marching, finite difference technique. Steady flow solutions are obtained from time asymptotic solutions for a steady airfoil. Corresponding oscillatory solutions are obtained by initiating an oscillation and marching in time for several cycles until a converged periodic solution is achieved. In this paper the method is described in general terms, and results are compared with experimental data for both steady flow and for oscillations at several values of reduced frequency. Good agreement for static pressures is shown for subcritical speeds, with increasing deviation as Mach number is increased into the supercritical speed range. Fair agreement with experiment was obtained at high reduced frequencies with larger deviations at low reduced frequencies.

Bennett, R. M.↗

Numerical calculation of transonic axial turbomachinery flows

A numerical method and the results of a computer program are presented for solving an exact, three dimensional, full potential equation that models rotating and nonrotating inviscid, absolutely irrotational, homentropic flows. Besides calculating the flows through an arbitrarily shaped rotor or stator blade row mounted on an axisymmetric hub and confined in an axisymmetric duct, the computer program is also capable of analyzing flow fields about arbitrarily shaped wing body combinations, propellers, helicopter rotors in hover, and wind turbine rotors. The governing equation is solved numerically in a fully conservative form by using an artificial time concept, a finite volume technique, rotated type dependent differencing, successive line overrelaxation, and sequential boundary conforming grid refinement. An artificial viscosity is added in fully conservative form, and an initial guess for the potential field is applied, as determined by a two dimensional cascade analysis.

Dulikravich, D. S.↗

Comparison between optical measurements and a numerical solution of the flow field within a transonic axial-flow compressor rotor

A comparison between numerical and experimental results is presented for the flowfield within a transonic axial-flow compressor rotor. The rotor was tested at design speed and a wide open throttle discharge condition. The relative tip Mach number was 1.4. A laser anemometer system was used to measure velocity and flow angle upstream, within, and downstream of the rotor. A holographic interferometer was used to visualize the rotor shock system near the tip. The computational procedure solves the full three-dimensional Euler equations using a time-marching technique. Shock location and shape determined from the two optical systems are compared. Calculated relative Mach number and flow angle contours, shock locations, and shock strength are compared to values measured with the laser anemometer.

Strazisar, A. J.↗