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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Characterization of Martian near-subsurface materials by determination of cohesion and angle of internal friction

Back-analysis (reconstruction) of the stability of thirty avalanche chutes was performed in the very limited areas where high resolution imaging overlapped with available 1:500 K topographic map coverage. A new technique was developed to incorporate the third dimension (width) of an avalanche chute in stability back-analysis in order to yield unambiguous values of cohesion and angle of internal friction. The procedure is based upon extending the ordinary method of slices to three dimensions, in order to construct avalanche chute cross-sections whose widths and depths vary as a function of gradient, gravity, density of material, and phi and c. Applying the technique to the well documented slide at Lodalen, Norway as a test produces excellent correspondence with reality. Generally, the technique reveals that the width:depth ratio of any avalanche chute decreases with increasing contrast between the average slope angle and the angle of internal friction. Applying this technique to the martian avalanche chute yields results consistent with indications from earlier work, but with greater certainty. Values of cohesion and angle of internal friction identify the materials at the time of failure as moderately cohesive debris. If Sharp's identification of these features as avalanche chutes is correct, then the results here imply that weathering processes have had a significant effect to depths of tens of meters (where failure has occured) below the martian surface. It is also implied that on relatively steep slopes within Valles Marineris, sizable, unaltered, unmantled bedrock exposures for high resolution spectral and spatial scanning by Mars Observer may be scarce.

Sullivan, R. J.↗
NASA NTRS1992

Microwave temperature profiler for clear air turbulence prediction

A method is disclosed for determining Richardson Number, Ri, or its reciprocal, RRi, for clear air prediction using measured potential temperature and determining the vertical gradient of potential temperature, d(theta)/dz. Wind vector from the aircraft instrumentation versus potential temperature, dW/D(theta), is determined and multiplies by d(theta)/dz to obtain dW/dz. Richardson number or its reciprocal is then determined from the relationship Ri = K(d theta)/dz divided by (dW/dz squared) for use in detecting a trend toward a threshold value for the purpose of predicting clear air turbulence. Other equations for this basic relationship are disclosed together with the combination of other atmospheric observables using multiple regression techniques.

Gary, Bruce L.↗
NASA NTRS1988

Design of a lunar transportation system

The establishment of lunar bases is the next logical step in the exploration of space. Permanent lunar bases will support scientific investigation, the industrialization of space, and the development of self-sufficiency on the Moon. Scientific investigation and research and development would lead to applications utilizing lunar material resources. By utilizing these resources, the industrialization of space can become a reality. The above two factors coupled with the development of key and enabling technologies would lead to achievement of self-sufficiency of the lunar base. Attention was focused on specific design(s) to be pursued during subsequent stages in advanced courses. Some of the objectives in the project included: (1) minimizing the transportation of construction material and fuel from earth, or maximizing the use of the lunar material; (2) use of novel materials and light weight structures; (3) use of new manufacturing methods and technology such as magnetically levitated, or superconducting materials; and (4) innovative concepts of effectively utilizing the exotic lunar conditions, i.e. high thermal gradients, lack of atmosphere, zero wind forces, and lower gravity, etc.

Source record↗
NASA NTRS1989

Design of a lunar transportation system

The development of a good transportation infrastructure is a major requirement for the establishment of a permanent lunar base. Transportation is characterized by the technology available in a specific time frame and the need to transport personnel and cargo between Earth and Moon, and between lunar bases. In our study, attention was first focused on developing a transportation system for the first generation lunar base. As a first step, a tracked-type multipurpose lunar transportation vehicle was considered as a possible mode of transportation and a detailed study was conducted on the various aspects of the vehicle. Since the vehicle is composed of many moving parts, exposing it to the environment of the Moon, where fine dust particles are prevalent, can cause problems associated with lubrication and friction. The vehicle also posed problems concerning weight and power. Hence, several modifications were made to the above design ideas conceptually, and a Lunar Articulated Remote Transportation System (Lunar ARTS) is proposed as a more effective alternative with the following objectives: (1) minimizing the transportation of construction material and fuel from Earth or maximizing the use of the lunar material; (2) use of novel materials and light-weight structures; (3) use of new manufacturing methods and technology such as magnetic levitation using superconducting materials; and (4) innovative concepts of effectively utilizing the exotic lunar conditions, i.e., high thermal gradients, lack of atmosphere, lower gravity, etc. To achieve the above objectives of designing transportation systems from concept to operation, the project was planned in three phases: (1) conceptual design; (2) detailed analysis and synthesis; and (3) construction, testing, evaluation, and operation. In this project, both phases 1 and 2 have been carried out and work on phase 3 is in progress. In this paper, the details of the Lunar ARTS are discussed and the future work on the vehicle are also outlined.

Sankaravelu, A.↗
NASA NTRS1994

Thermocapillary motion of deformable drops

The thermocapillary motion of initially spherical drops/bubbles driven by a constant temperature gradient in an unbounded liquid medium is simulated numerically. Effects of convection of momentum and energy, as well as shape deformations, are addressed. The method used is based on interface tracking on a base cartesian grid, and uses a smeared color or indicator function for the determination of the surface topology. Quad-tree adaptive refinement of the cartesian grid is implemented to enhance the fidelity of the surface tracking. It is shown that convection of energy results in a slowing of the drop, as the isotherms get wrapped around the front of the drop. Shape deformation resulting from inertial effects affect the migration velocity. The physical results obtained are in agreement with the existing literature. Furthermore, remarks are made on the sensitivity of the calculated solutions to the smearing of the fluid properties. Analysis and simulations show that the migration velocity depends very strongly on the smearing of the interfacial force whereas it is rather insensitive to the smearing of other properties, hence the adaptive grid.

Haj-Hariri, Hossein↗
NASA NTRS2002

Application of Accelerometer Data to Mars Odyssey Aerobraking and Atmospheric Modeling

Aerobraking was an enabling technology for the Mars Odyssey mission even though it involved risk due primarily to the variability of the Mars upper atmosphere. Consequently, numerous analyses based on various data types were performed during operations to reduce these risk and among these data were measurements from spacecraft accelerometers. This paper reports on the use of accelerometer data for determining atmospheric density during Odyssey aerobraking operations. Acceleration was measured along three orthogonal axes, although only data from the component along the axis nominally into the flow was used during operations. For a one second count time, the RMS noise level varied from 0.07 to 0.5 mm/s2 permitting density recovery to between 0.15 and 1.1 kg per cu km or about 2% of the mean density at periapsis during aerobraking. Accelerometer data were analyzed in near real time to provide estimates of density at periapsis, maximum density, density scale height, latitudinal gradient, longitudinal wave variations and location of the polar vortex. Summaries are given of the aerobraking phase of the mission, the accelerometer data analysis methods and operational procedures, some applications to determining thermospheric properties, and some remaining issues on interpretation of the data. Pre-flight estimates of natural variability based on Mars Global Surveyor accelerometer measurements proved reliable in the mid-latitudes, but overestimated the variability inside the polar vortex.

Tolson, R. H.↗
NASA NTRS2004

Employing Sensitivity Derivatives for Robust Optimization under Uncertainty in CFD

A robust optimization is demonstrated on a two-dimensional inviscid airfoil problem in subsonic flow. Given uncertainties in statistically independent, random, normally distributed flow parameters (input variables), an approximate first-order statistical moment method is employed to represent the Computational Fluid Dynamics (CFD) code outputs as expected values with variances. These output quantities are used to form the objective function and constraints. The constraints are cast in probabilistic terms; that is, the probability that a constraint is satisfied is greater than or equal to some desired target probability. Gradient-based robust optimization of this stochastic problem is accomplished through use of both first and second-order sensitivity derivatives. For each robust optimization, the effect of increasing both input standard deviations and target probability of constraint satisfaction are demonstrated. This method provides a means for incorporating uncertainty when considering small deviations from input mean values.

Newman, Perry A.↗
NASA NTRS2016

Temperature Dependence of Electrical Resistance of Woven Melt-Infiltrated SiCf/SiC Ceramic Matrix Composites

Recent studies have successfully shown the use of electrical resistance (ER)measurements to monitor room temperature damage accumulation in SiC fiber reinforced SiC matrix composites (SiCf/SiC) Ceramic Matrix Composites (CMCs). In order to determine the feasibility of resistance monitoring at elevated temperatures, the present work investigates the temperature dependent electrical response of various MI (Melt Infiltrated)-CVI (Chemical Vapor Infiltrated) SiC/SiC composites containing Hi-Nicalon Type S, Tyranno ZMI and SA reinforcing fibers. Test were conducted using a commercially available isothermal testing apparatus as well as a novel, laser-based heating approach developed to more accurately simulate thermomechanical testing of CMCs. Secondly, a post-test inspection technique is demonstrated to show the effect of high-temperature exposure on electrical properties. Analysis was performed to determine the respective contribution of the fiber and matrix to the overall composite conductivity at elevated temperatures. It was concluded that because the silicon-rich matrix material dominates the electrical response at high temperature, ER monitoring would continue to be a feasible method for monitoring stress dependent matrix cracking of melt-infiltrated SiC/SiC composites under high temperature mechanical testing conditions. Finally, the effect of thermal gradients generated during localized heating of tensile coupons on overall electrical response of the composite is determined.

Ceramic Matrix Composites↗
NASA NTRS2017

Constrained Burn Optimization for the International Space Station

In long-term trajectory planning for the International Space Station (ISS), translational burns are currently targeted sequentially to meet the immediate trajectory constraints, rather than simultaneously to meet all constraints, do not employ gradient-based search techniques, and are not optimized for a minimum total deltav (v) solution. An analytic formulation of the constraint gradients is developed and used in an optimization solver to overcome these obstacles. Two trajectory examples are explored, highlighting the advantage of the proposed method over the current approach, as well as the potential v and propellant savings in the event of propellant shortages.

Brown, Aaron J.↗
NASA NTRSDate not supplied

Nonlinear Homogenization of Finitely Deformed Viscoelastic-Viscoplastic Composites Using Mechanics of Structure Genome

The objective of this paper is to develop a micromechanics approach to homogenizing finitely deformed viscoelastic-viscoplastic composites using the mechanics of structure genome. The incremental constitutive relation for glassy polymers, formulated in the spatial configuration, is implemented in the present approach.This involves (1) pulling-back the constitutive model to the material configuration and (2)choosing the deformation gradient tensor and the first Piola–Kirchhoff stress tensor as the strain and the stress measures during homogenization, respectively. An Euler–Newton predictor–corrector method is developed for homogenization. Each step involves formulating a variational statement using the mechanics of structure genome, discretizing the statement in a finite-dimensional space, and solving the problem using an Euler/multilevel Newton method. The present approach is demonstrated by homogenizing fiber- and particle-reinforced composites undergoing uniaxial, biaxial, or shear deformation, at different stain rates.

Multi-scale modeling, High Strain Composites, Visc↗
NASA NTRS1999

Demonstration of Automatically-Generated Adjoint Code for Use in Aerodynamic Shape Optimization

Gradient-based optimization requires accurate derivatives of the objective function and constraints. These gradients may have previously been obtained by manual differentiation of analysis codes, symbolic manipulators, finite-difference approximations, or existing automatic differentiation (AD) tools such as ADIFOR (Automatic Differentiation in FORTRAN). Each of these methods has certain deficiencies, particularly when applied to complex, coupled analyses with many design variables. Recently, a new AD tool called ADJIFOR (Automatic Adjoint Generation in FORTRAN), based upon ADIFOR, was developed and demonstrated. Whereas ADIFOR implements forward-mode (direct) differentiation throughout an analysis program to obtain exact derivatives via the chain rule of calculus, ADJIFOR implements the reverse-mode counterpart of the chain rule to obtain exact adjoint form derivatives from FORTRAN code. Automatically-generated adjoint versions of the widely-used CFL3D computational fluid dynamics (CFD) code and an algebraic wing grid generation code were obtained with just a few hours processing time using the ADJIFOR tool. The codes were verified for accuracy and were shown to compute the exact gradient of the wing lift-to-drag ratio, with respect to any number of shape parameters, in about the time required for 7 to 20 function evaluations. The codes have now been executed on various computers with typical memory and disk space for problems with up to 129 x 65 x 33 grid points, and for hundreds to thousands of independent variables. These adjoint codes are now used in a gradient-based aerodynamic shape optimization problem for a swept, tapered wing. For each design iteration, the optimization package constructs an approximate, linear optimization problem, based upon the current objective function, constraints, and gradient values. The optimizer subroutines are called within a design loop employing the approximate linear problem until an optimum shape is found, the design loop limit is reached, or no further design improvement is possible due to active design variable bounds and/or constraints. The resulting shape parameters are then used by the grid generation code to define a new wing surface and computational grid. The lift-to-drag ratio and its gradient are computed for the new design by the automatically-generated adjoint codes. Several optimization iterations may be required to find an optimum wing shape. Results from two sample cases will be discussed. The reader should note that this work primarily represents a demonstration of use of automatically- generated adjoint code within an aerodynamic shape optimization. As such, little significance is placed upon the actual optimization results, relative to the method for obtaining the results.

Green, Lawrence↗
NASA NTRS1995

Acoustic receptivity due to weak surface inhomogeneities in adverse pressure gradient boundary layers

The boundary layer receptivity to free-stream acoustic waves in the presence of localized surface disturbances is studied for the case of incompressible Falkner-Skan flows with adverse pressure gradients. These boundary layers are unstable to both viscous and inviscid (i.e., inflectional) modes, and the finite Reynolds number extension of the Goldstein-Ruban theory provides a convenient method to compare the efficiency of the localized receptivity processes in these two cases. The value of the efficiency function related to the receptivity caused by localized distortions in surface geometry is relatively insensitive to the type of instability mechanism, provided that the same reference length scale is used to normalize the efficiency function for each type of instability. In contrast, when the receptivity is induced by variations in wall suction velocity or in wall admittance distribution, the magnitudes of the related efficiency functions, as well as the resulting coupling coefficients, are smaller for inflectional (i.e., Rayleigh) modes than for the viscous Tollmien-Schlichting waves. The reduced levels of receptivity can be attributed mainly to the shorter wavelengths and higher frequencies of the inflectional modes. Because the most critical band of frequencies shifts toward higher values, the overall efficiency of the wall suction- and the wall admittance-induced receptivity decreases with an increase in the adverse pressure gradient.

Choudhari, Meelan↗
NASA NTRS1993

Application of multigrid methods to the solution of liquid crystal equations on a SIMD computer

We will describe a finite difference code for computing the equilibrium configurations of the order-parameter tensor field for nematic liquid crystals in rectangular regions by minimization of the Landau-de Gennes Free Energy functional. The implementation of the free energy functional described here includes magnetic fields, quadratic gradient terms, and scalar bulk terms through the fourth order. Boundary conditions include the effects of strong surface anchoring. The target architectures for our implementation are SIMD machines, with interconnection networks which can be configured as 2 or 3 dimensional grids, such as the Wavetracer DTC. We also discuss the relative efficiency of a number of iterative methods for the solution of the linear systems arising from this discretization on such architectures.

Farrell, Paul A.↗
NASA NTRSDate not supplied

Edge Based Viscous Method for Node-Centered Formulations

This paper presents a novel, efficient, conservative, edge-based method for evaluation of mean flow viscous fluxes and turbulence-model diffusion terms of the Reynolds-averaged Navier-Stokes equations on tetrahedral grids. The new method is implemented in a practical, node-centered, finite-volume computational fluid dynamics solver. The baseline finite-volume scheme that is equivalent to a second-order accurate finite-element Galerkin approximation of viscous stresses is reformulated. The order of operations to compute the cell-based Green-Gauss gradients is changed to combine the operations by edge, which leads to an equivalent formulation on tetrahedral grids, improves efficiency, and preserves the compact discretization stencil based on the nearest neighbors. The computational results presented in this paper verify the implementation of this edge-based method by comparing its accuracy and iterative convergence with those of the well verified and validated baseline formulation. Efficiency gains for residual and Jacobian evaluations result in significant reduction of time to solution. This novel edge-based formulation on tetrahedra can be seamlessly combined with the baseline formulation on cells of other types for computing solutions on mixed-element grids.

Edge Based↗
NASA NTRS1991

A look-ahead variant of the Lanczos algorithm and its application to the quasi-minimal residual method for non-Hermitian linear systems

The Lanczos algorithm can be used both for eigenvalue problems and to solve linear systems. However, when applied to non-Hermitian matrices, the classical Lanczos algorithm is susceptible to breakdowns and potential instabilities. In addition, the biconjugate gradient (BCG) algorithm, which is the natural generalization of the conjugate gradient algorithm to non-Hermitian linear systems, has a second source of breakdowns, independent of the Lanczos breakdowns. Here, we present two new results. We propose an implementation of a look-ahead variant of the Lanczos algorithm which overcomes the breakdowns by skipping over those steps where a breakdown or a near-breakdown would occur. The new algorithm can handle look-ahead steps of any length and requires the same number of matrix-vector products and inner products per step as the classical Lanczos algorithm without look-ahead. Based on the proposed look-ahead Lanczos algorithm, we then present a novel BCG-like approach, the quasi-minimal residual (QMR) method, which avoids the second source of breakdowns in the BCG algorithm. We present details of the new method and discuss some of its properties. In particular, we discuss the relationship between QMR and BCG, showing how one can recover the BCG iterates, when they exist, from the QMR iterates. We also present convergence results for QMR, showing the connection between QMR and the generalized minimal residual (GMRES) algorithm, the optimal method in this class of methods. Finally, we give some numerical examples, both for eigenvalue computations and for non-Hermitian linear systems.

Nachtigal, Noel M.↗
NASA NTRS2010

Receptivity of Hypersonic Boundary Layers over Straight and Flared Cones

The effects of adverse pressure gradients on the receptivity and stability of hypersonic boundary layers were numerically investigated. Simulations were performed for boundary layer flows over a straight cone and two flared cones. The steady and the unsteady flow fields were obtained by solving the two-dimensional Navier-Stokes equations in axi-symmetric coordinates using the 5th order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The mean boundary layer profiles were analyzed using local stability and non-local parabolized stability equations (PSE) methods. After the most amplified disturbances were identified, two-dimensional plane acoustic waves were introduced at the outer boundary of the computational domain and time accurate simulations were performed. The adverse pressure gradient was found to affect the boundary layer stability in two important ways. Firstly, the frequency of the most amplified second-mode disturbance was increased relative to the zero pressure gradient case. Secondly, the amplification of first- and second-mode disturbances was increased. Although an adverse pressure gradient enhances instability wave growth rates, small nose-tip bluntness was found to delay transition due to the low receptivity coefficient and the resulting weak initial amplitude of the instability waves. The computed and measured amplitude-frequency spectrums in all three cases agree very well in terms of frequency and the shape except for the amplitude.

Balakumar, Ponnampalam↗
NASA NTRS1982

Numerical simulations of unsteady transonic flow in diffusers

Forced and naturally occurring, self-sustaining oscillations of transonic flows in two-dimensional diffusers were computed using MacCormack's hybrid method. Depending upon the shock strengths and the area ratios, the flow was fully attached or separated by either the shock or the adverse pressure gradient associated with the enlarging diffuser area. In the case of forced oscillations, a sinusoidal plane pressure wave at frequency 300 Hz was prescribed at the exit. A sufficiently large amount of data were acquired and Fourier analyzed. The distrbutions of time-mean pressures, the power spectral density, and the amplitude with phase angle along the top wall and in the core region were determined. Comparison with experimental results for the forced oscillation generally gave very good agreement; some success was achieved for the case of self-sustaining oscillation despite substantial three-dimensionality in the test. An observation of the sequence of self-sustaining oscillations was given.

Liou, M.-S.↗
NASA NTRS1989

Graphical postprocessing for 3-D mesh quality evaluation

An important objective of three-dimensional graphical finite-element postprocessing is to indicate to the engineer the accuracy of analysis results. The inclusion of mesh quality sensors permits a subjective evaluation of the adequacy of a single analysis being interpreted. For graphical approaches, both strain-energy-density gradients and discontinuities of unsmoothed responses and their gradients have proved to be effective sensors. Interactive graphical tools which can display discontinuity information effectively are described; these are essentially different from the ordinary methods used for the viewing of smoothed results.

Panthaki, M. J.↗