Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gradient methods”

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 649 records · Page 36

Heat-Transfer and Pressure Measurements on a Flat-Face Cylinder at a Mach Number Range of 2.49 to 4.44

Heat-transfer coefficients and pressure distributions were obtained on a 4-inch-diameter flat-face cylinder in the Langley Unitary Plan wind tunnel. The measured stagnation heat-transfer coefficient agrees well with 55 percent of the theoretical value predicted by the modified Sibulkin method for a hemisphere. Pressure measurements indicated the dimensionless velocity gradient parameter r du\ a(sub t) dx, where x=0 at the stagnation point was approximately 0.3 and invariant throughout the Mach number range from 2.49 to 4.44 and the Reynolds number range from 0.77 x 10(exp 6) to 1.46 x 10(exp 6). The heat-transfer coefficients on the cylindrical afterbody could be predicted with reasonable accuracy by flat-plate theory at an angle of attack of 0 deg. At angles of attack the cylindrical afterbody stagnation-line heat transfer could be computed from swept-cylinder theory for large distances back of the nose when the Reynolds number is based on the distance from the flow reattachment points.

Burbank, Paige B.↗

Development of Metal Matrix Composites for NASA'S Advanced Propulsion Systems

The state-of-the-art development of several aluminum and copper based Metal Matrix Composites (MMC) for NASA's advanced propulsion systems will be presented. The presentation's goal is to provide an overview of NASA-Marshall Space Flight Center's planned and on-going activities in MMC for advanced liquid rocket engines such as the X-33 vehicle's Aerospike and X-34 Fastrac engine. The focus will be on lightweight and environmental compatibility with oxygen and hydrogen of key MMC materials, within each NASA's new propulsion application, that will provide a high payoff for NASA's reusable launch vehicle systems and space access vehicles. Advanced MMC processing techniques such as plasma spray, centrifugal casting, pressure infiltration casting will be discussed. Development of a novel 3D printing method for low cost production of composite preform, and functional gradient MMC to enhanced rocket engine's dimensional stability will be presented.

Lee, Jonathan A.↗

Computational Relativistic Astrophysics Using the Flowfield-Dependent Variation Theory

Theoretical models, observations and measurements have preoccupied astrophysicists for many centuries. Only in recent years, has the theory of relativity as applied to astrophysical flows met the challenges of how the governing equations can be solved numerically with accuracy and efficiency. Even without the effects of relativity, the physics of magnetohydrodynamic flow instability, turbulence, radiation, and enhanced transport in accretion disks has not been completely resolved. Relativistic effects become pronounced in such cases as jet formation from black hole magnetized accretion disks and also in the study of Gamma-Ray bursts (GRB). Thus, our concern in this paper is to reexamine existing numerical simulation tools as to the accuracy and efficiency of computations and introduce a new approach known as the flowfield-dependent variation (FDV) method. The main feature of the FDV method consists of accommodating discontinuities of shock waves and high gradients of flow variables such as occur in turbulence and unstable motions. In this paper, the physics involved in the solution of relativistic hydrodynamics and solution strategies of the FDV theory are elaborated. The general relativistic astrophysical flow and shock solver (GRAFSS) is introduced, and some simple example problems for Computational Relativistic Astrophysics (CRA) are demonstrated.

Richardson, G. A.↗

Computational Relativistic Astrophysics Using the Flow Field-Dependent Variation Theory

We present our method for solving general relativistic nonideal hydrodynamics. Relativistic effects become pronounced in such cases as jet formation from black hole magnetized accretion disks which may lead to the study of gamma-ray bursts. Nonideal flows are present where radiation, magnetic forces, viscosities, and turbulence play an important role. Our concern in this paper is to reexamine existing numerical simulation tools as to the accuracy and efficiency of computations and introduce a new approach known as the flow field-dependent variation (FDV) method. The main feature of the FDV method consists of accommodating discontinuities of shock waves and high gradients of flow variables such as occur in turbulence and unstable motions. In this paper, the physics involved in the solution of relativistic hydrodynamics and solution strategies of the FDV theory are elaborated. The general relativistic astrophysical flow and shock solver (GRAFSS) is introduced, and some simple example problems for computational relativistic astrophysics (CRA) are demonstrated.

Richardson, G. A.↗

Thermal Model Correlation for Mars Reconnaissance Orbiter

The Mars Reconnaissance Orbiter (MRO) launched on August 12, 2005 and began aerobraking at Mars in March 2006. In order to save propellant, MRO used aerobraking to modify the initial orbit at Mars. The spacecraft passed through the atmosphere briefly on each orbit; during each pass the spacecraft was slowed by atmospheric drag, thus lowering the orbit apoapsis. The largest area on the spacecraft, most affected by aeroheating, was the solar arrays. A thermal analysis of the solar arrays was conducted at NASA Langley Research Center to simulate their performance throughout the entire roughly 6-month period of aerobraking. A companion paper describes the development of this thermal model. This model has been correlated against many sets of flight data. Several maneuvers were performed during the cruise to Mars, such as thruster calibrations, which involve large abrupt changes in the spacecraft orientation relative to the sun. The data obtained from these maneuvers allowed the model to be well-correlated with regard to thermal mass, conductive connections, and solar response well before arrival at the planet. Correlation against flight data for both in-cruise maneuvers and drag passes was performed. Adjustments made to the model included orientation during the drag pass, solar flux, Martian surface temperature, through-array resistance, aeroheating gradient due to angle of attack, and aeroheating accommodation coefficient. Methods of correlation included comparing the model to flight temperatures, slopes, temperature deltas between sensors, and solar and planet direction vectors. Correlation and model accuracy over 400 aeroheating drag passes were determined, with overall model accuracy better than 5 C.

Amundsen, Ruth M.↗

Flow Separation Control Over a Ramp Using Sweeping Jet Actuators

Flow separation control on an adverse-pressure-gradient ramp model was investigated using various flow-control methods in the NASA Langley 15-Inch Wind Tunnel. The primary flow-control method studied used a sweeping jet actuator system to compare with more classic flow-control techniques such as micro-vortex generators, steady blowing, and steady- and unsteady-vortex generating jets. Surface pressure measurements and a new oilflow visualization technique were used to characterize the effects of these flow-control actuators. The sweeping jet actuators were run in three different modes to produce steady-straight, steady-angled, and unsteady-oscillating jets. It was observed that all of these flow-control methods are effective in controlling the separated flows on the ramp model. The steady-straight jet energizes the boundary layer by momentum addition and was found to be the least effective method for a fixed momentum coefficient. The steady-angled jets achieved better performance than the steady-straight jets because they generate streamwise vortices that energize the boundary layer by mixing high-momentum fluid with near wall low-momentum fluid. The unsteady-oscillating jets achieved the best performance by increasing the pressure recovery and reducing the downstream flow separation. Surface flow visualizations indicated that two out-of-phase counter-rotating vortices are generated per sweeping jet actuator, while one vortex is generated per vortex-generating jets. The extra vortex resulted in increased coverage, more pressure recovery, and reduced flow separation.

Koklu, Mehti↗

Efficient On-Orbit Singularity-Free Geopotential Estimation

For vehicles in orbital cruise, either in orbit around a body or in a transfer orbit, gravity is the primary external force; therefore, the complexity of the geopotential model heavily influences the accuracy of the navigation state. As this type of spacecraft requires a highly accurate geopotential model and performs the computation numerous times a second, this calculation needs to be efficient. This calculation needs to include the first and second derivatives of the geopotential because navigation systems can require both the gravitational acceleration and gradient for inertial integration and state filtering. The most efficient method for calculating the geopotential, the forward column method, is discussed in detail, as well as a method to avoid the singularities that exist when using this method. This is shown to decrease the computation time of the geopotential compared to other popular methods. In addition, methods for first and second order propagation are discussed, which decrease the rate at which the full geopotential model needs to be calculated while maintaining the accuracy required of navigation systems. These estimation methods are shown to decrease the necessary computation of the gravity model by multiple orders of magnitude. This method decreases the computation time of the Exploration Upper Stage navigation system geopotential model by potentially an order of magnitude compared with the previous model without affecting the navigation error. A properly implemented geopotential model can have the accuracy of an 8 by 8 model while approaching the computational requirement as using only the J234 coefficients.

Amert, Joel↗

Approach for Input Uncertainty Propagation and Robust Design in CFD Using Sensitivity Derivatives

An implementation of the approximate statistical moment method for uncertainty propagation and robust optimization for quasi 3-D Euler CFD code is presented. Given uncertainties in statistically independent, random, normally distributed input variables, first- and second-order statistical moment procedures are performed to approximate the uncertainty in the CFD output. Efficient calculation of both first- and second-order sensitivity derivatives is required. In order to assess the validity of the approximations, these moments are compared with statistical moments generated through Monte Carlo simulations. The uncertainties in the CFD input variables are also incorporated into a robust optimization procedure. For this optimization, statistical moments involving first-order sensitivity derivatives appear in the objective function and system constraints. Second-order sensitivity derivatives are used in a gradient-based search to successfully execute a robust optimization. The approximate methods used throughout the analyses are found to be valid when considering robustness about input parameter mean values.

Putko, Michele M.↗

Analyses of large quasistatic deformations of inelastic bodies by a new hybrid-stress finite element algorithm - Applications

A new hybrid-stress finite element algorithm suitable for analyzing large quasistatic deformations of inelastic solids is presented and its feasibility and performance are demonstrated with examples. The algorithm provides extremely accurate bifurcation analysis which is stable with respect to variation in the finite element mesh, so long as the same type of element is used in every mesh. When the mesh element is varied, the result changes in a predictable manner. The method does not necessarily lead to an upper or lower bound for the critical load. An explicit forward gradient scheme is used to improve stability and is shown to be useful also for elongation-dominated deformations. The application of the method to the onset of necking in plane extension and to deformation and stress in plane extension of an elasticoviscous fluid with an array of cylindrical voids is given in detail.

Reed, K. W.↗

Charts and Tables for Estimating the Stability of the Compressible Laminar Boundary Layer with Heat Transfer and Arbitrary Pressure Gradient

The minimum critical Reynolds numbers for the similar solutions of the compressible laminar boundary layer computed by Cohen and Reshotko and also for the Falkner and Skan solutions as recomputed by Smith have been calculated by Lin's rapid approximate method for two-dimensional disturbances. These results enable the stability of the compressible laminar boundary layer with heat transfer and pressure gradient to be easily estimated after the behavior of the boundary layer has been computed by the approximate method of Cohen and Reshotko. The previously reported unusual result (NACA Technical Note 4037) that a highly cooled stagnation point flow is more unstable than a highly cooled flat-plate flow is again encountered. Moreover, this result is found to be part of the more general result that a favorable pressure gradient is destabilizing for very cool walls when the Mach number is less than that for complete stability. The minimum critical Reynolds numbers for these wall temperature ratios are, however, all larger than any value of the laminar-boundary-layer Reynolds number likely to be encountered. For Mach numbers greater than those for which complete stability occurs a favorable pressure gradient is stabilizing, even for very cool walls.

Tetervin, Neal↗

Approach for Uncertainty Propagation and Robust Design in CFD Using Sensitivity Derivatives

This paper presents an implementation of the approximate statistical moment method for uncertainty propagation and robust optimization for a quasi 1-D Euler CFD (computational fluid dynamics) code. Given uncertainties in statistically independent, random, normally distributed input variables, a first- and second-order statistical moment matching procedure is performed to approximate the uncertainty in the CFD output. Efficient calculation of both first- and second-order sensitivity derivatives is required. In order to assess the validity of the approximations, the moments are compared with statistical moments generated through Monte Carlo simulations. The uncertainties in the CFD input variables are also incorporated into a robust optimization procedure. For this optimization, statistical moments involving first-order sensitivity derivatives appear in the objective function and system constraints. Second-order sensitivity derivatives are used in a gradient-based search to successfully execute a robust optimization. The approximate methods used throughout the analyses are found to be valid when considering robustness about input parameter mean values.

Putko, Michele M.↗

Effects of specularly reflected radiation on spacecraft temperatures and thermal gradients

This paper describes the effect that specularly reflected solar energy has upon the heating load imposed upon orbiting spacecraft. Because this reflection may increase the total heating by factors of two or three, it is important that it can be computed accurately. An efficient method for treating multiple reflections is given and demonstrated by computing the temperatures and thermal gradients in a reflecting spherical radiator.

Emery, A. F.↗

Aerodynamic Shape Optimization of Supersonic Aircraft Configurations via an Adjoint Formulation on Parallel Computers

This work describes the application of a control theory-based aerodynamic shape optimization method to the problem of supersonic aircraft design. The design process is greatly accelerated through the use of both control theory and a parallel implementation on distributed memory computers. Control theory is employed to derive the adjoint differential equations whose solution allows for the evaluation of design gradient information at a fraction of the computational cost required by previous design methods. The resulting problem is then implemented on parallel distributed memory architectures using a domain decomposition approach, an optimized communication schedule, and the MPI (Message Passing Interface) Standard for portability and efficiency. The final result achieves very rapid aerodynamic design based on higher order computational fluid dynamics methods (CFD). In our earlier studies, the serial implementation of this design method was shown to be effective for the optimization of airfoils, wings, wing-bodies, and complex aircraft configurations using both the potential equation and the Euler equations. In our most recent paper, the Euler method was extended to treat complete aircraft configurations via a new multiblock implementation. Furthermore, during the same conference, we also presented preliminary results demonstrating that this basic methodology could be ported to distributed memory parallel computing architectures. In this paper, our concern will be to demonstrate that the combined power of these new technologies can be used routinely in an industrial design environment by applying it to the case study of the design of typical supersonic transport configurations. A particular difficulty of this test case is posed by the propulsion/airframe integration.

Reuther, James↗

Aerodynamic Shape Optimization of Supersonic Aircraft Configurations via an Adjoint Formulation on Parallel Computers

This work describes the application of a control theory-based aerodynamic shape optimization method to the problem of supersonic aircraft design. The design process is greatly accelerated through the use of both control theory and a parallel implementation on distributed memory computers. Control theory is employed to derive the adjoint differential equations whose solution allows for the evaluation of design gradient information at a fraction of the computational cost required by previous design methods (13, 12, 44, 38). The resulting problem is then implemented on parallel distributed memory architectures using a domain decomposition approach, an optimized communication schedule, and the MPI (Message Passing Interface) Standard for portability and efficiency. The final result achieves very rapid aerodynamic design based on higher order computational fluid dynamics methods (CFD). In our earlier studies, the serial implementation of this design method (19, 20, 21, 23, 39, 25, 40, 41, 42, 43, 9) was shown to be effective for the optimization of airfoils, wings, wing-bodies, and complex aircraft configurations using both the potential equation and the Euler equations (39, 25). In our most recent paper, the Euler method was extended to treat complete aircraft configurations via a new multiblock implementation. Furthermore, during the same conference, we also presented preliminary results demonstrating that the basic methodology could be ported to distributed memory parallel computing architectures [241. In this paper, our concem will be to demonstrate that the combined power of these new technologies can be used routinely in an industrial design environment by applying it to the case study of the design of typical supersonic transport configurations. A particular difficulty of this test case is posed by the propulsion/airframe integration.

Reuther, James↗

Surface, Water, and Air Biocharacterization (SWAB) Flight Experiment

The determination of risk from infectious disease during spaceflight missions is composed of several factors including both the concentration and characteristics of the microorganisms to which the crew are exposed. Thus, having a good understanding of the microbial ecology aboard spacecraft provides the necessary information to mitigate health risks to the crew. While preventive measures are taken to minimize the presence of pathogens on spacecraft, medically significant organisms have been isolated from both the Mir and International Space Station (ISS). Historically, the method for isolation and identification of microorganisms from spacecraft environmental samples depended upon their growth on culture media. Unfortunately, only a fraction of the organisms may grow on a specific culture medium, potentially omitting those microorganisms whose nutritional and physical requirements for growth are not met. To address this bias in our understanding of the ISS environment, the Surface, Water, and Air Biocharacterization (SWAB) Flight Experiment was designed to investigate and develop monitoring technology to provide better microbial characterization. For the SWAB flight experiment, we hypothesized that environmental analysis using non-culture-based technologies would reveal microorganisms, allergens, and microbial toxins not previously reported in spacecraft, allowing for a more complete health assessment. Key findings during this experiment included: a) Generally, advanced molecular techniques were able to reveal a few organisms not recovered using culture-based methods; however, there is no indication that current monitoring is "missing" any medically significant bacteria or fungi. b) Molecular techniques have tremendous potential for microbial monitoring, however, sample preparation and data analysis present challenges for spaceflight hardware. c) Analytical results indicate that some molecular techniques, such as denaturing gradient gel electrophoresis (DGGE), can be much less sensitive than culture-based methods. d) More sensitive molecular techniques, such as quantitative polymerase chain reaction (QPCR), were able to identify viral DNA from ISS environments, suggesting potential transfer of the organism between crewmembers. In addition, the hardware selected for this experiment represented advances for next-generation sample collection. The advanced nature of this collection hardware was noted, when the Sartorius MD8 Air Port air sampler from the SWAB experiment remained on board ISS at the request of JAXA investigators, who intend to use it in completion of their microbial ecology experiment.

Castro, V. A.↗

Incompressible spectral-element method: Derivation of equations

A fractional-step splitting scheme breaks the full Navier-Stokes equations into explicit and implicit portions amenable to the calculus of variations. Beginning with the functional forms of the Poisson and Helmholtz equations, we substitute finite expansion series for the dependent variables and derive the matrix equations for the unknown expansion coefficients. This method employs a new splitting scheme which differs from conventional three-step (nonlinear, pressure, viscous) schemes. The nonlinear step appears in the conventional, explicit manner, the difference occurs in the pressure step. Instead of solving for the pressure gradient using the nonlinear velocity, we add the viscous portion of the Navier-Stokes equation from the previous time step to the velocity before solving for the pressure gradient. By combining this 'predicted' pressure gradient with the nonlinear velocity in an explicit term, and the Crank-Nicholson method for the viscous terms, we develop a Helmholtz equation for the final velocity.

Deanna, Russell G.↗

Adjoint Sensitivity Computations for an Embedded-Boundary Cartesian Mesh Method and CAD Geometry

Cartesian-mesh methods are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric Computer-Aided Design (CAD) tools. Our goal is to combine the automation capabilities of Cartesian methods with an eficient computation of design sensitivities. We address this issue using the adjoint method, where the computational cost of the design sensitivities, or objective function gradients, is esseutially indepeudent of the number of design variables. In previous work, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm included the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Central to this development is the computation of volume-mesh sensitivities to obtain a reliable approximation of the objective finction gradient. Motivated by the success of mesh-perturbation schemes commonly used in body-fitted unstructured formulations, we propose an approach based on a local linearization of a mesh-perturbation scheme similar to the spring analogy. This approach circumvents most of the difficulties that arise due to non-smooth changes in the cut-cell layer as the boundary shape evolves and provides a consistent approximation tot he exact gradient of the discretized abjective function. A detailed gradient accurace study is presented to verify our approach. Thereafter, we focus on a shape optimization problem for an Apollo-like reentry capsule. The optimization seeks to enhance the lift-to-drag ratio of the capsule by modifyjing the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design. This abstract presents only a brief outline of the numerical method and results; full details will be given in the final paper.

Nemec, Marian↗

Impedance Eduction in a Duct Using the Linearized Euler Equations

An impedance eduction technique that is based upon a numerical solution to the linearized Euler equations coupled with the minimization of an objective function is presented. The linearized Euler equations are solved numerically using a cubic finite element method and the minimum of the objective function is obtained via a gradient based optimizer. The mean flow is allowed to have a gradient in two coordinate directions and the impedance of the liner may vary arbitrarily over the liner surface. A new regularized wall impedance boundary condition presented by Rienstra which contains the effects of the mean boundary layer thickness is implemented. Impedances are educed for a perforate over honeycomb liner using test data acquired in the Langley Grazing Flow Impedance Tube. The impedance of the liner is educed both with an upstream source and a downstream source. The method is validated by comparing the educed impedance to a benchmark method. The primary conclusion of the study is that the impedances educed for upstream and downstream sources are not well matched and that the primary effect of the boundary thickness is to reduce the resistance of the liner.

Watson, Willie R.↗