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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Large Eddy Simulations of a Single-Injector Cooling Flow Using the High-Order Flux Reconstruction Method

A single-injector cooling flow into a heated crossflow was used as a validation case for large eddy simulations (LES) from two high-order CFD codes with comparisons to a state-of-the-art Reynolds averaged Navier-Stokes (RANS) turbulence model. The focus of this paper is on the validation of LES for one of the high-order CFD codes, namely the GFR (Glenn Flux Reconstruction) code that is being developed at NASA Glenn Research Center. Fourth-order LES were performed for the single-injector cooling flow configuration at blowing ratios of both 1.0 and 2.0, and a fifth-order LES was also performed with a blowing ratio of 1.0. The GFR simulations agree very well with the experiment mean quantities and reasonably well with the experiment turbulence quantities. The LES solutions from GFR and the other high-order CFD code, FDL3DI, agreed very closely for nearly every flow quantity examined, despite the fact that these two codes use completely different numerical methods, different grid geometries and domains, and different inflow boundary conditions. Finally, both LES show a significant accuracy improvement over RANS turbulence models for flows with large temperature gradients, where the standard gradient diffusion approximation is insufficient for temperature predictions.

Large Eddy Simulation↗
NASA NTRS1980

Discrete-continuous variable structural synthesis using dual methods

Approximation concepts and dual methods are extended to solve structural synthesis problems involving a mix of discrete and continuous sizing type of design variables. Pure discrete and pure continuous variable problems can be handled as special cases. The basic mathematical programming statement of the structural synthesis problem is converted into a sequence of explicit approximate primal problems of separable form. These problems are solved by constructing continuous explicit dual functions, which are maximized subject to simple nonnegativity constraints on the dual variables. A newly devised gradient projection type of algorithm called DUAL 1, which includes special features for handling dual function gradient discontinuities that arise from the discrete primal variables, is used to find the solution of each dual problem. Computational implementation is accomplished by incorporating the DUAL 1 algorithm into the ACCESS 3 program as a new optimizer option. The power of the method set forth is demonstrated by presenting numerical results for several example problems, including a pure discrete variable treatment of a metallic swept wing and a mixed discrete-continuous variable solution for a thin delta wing with fiber composite skins.

Schmit, L. A.↗
NASA NTRS1984

Special Test Methods for Batteries

Various methods are described for measuring heat generation in primary and secondary batteries as well as the specific heat of batteries and cell thermal conductance. Problems associated with determining heat generation in large batteries are examined. Special attention is given to monitoring temperature gradients in nickel cadmium cells, the use of auxiliary electrodes for conducting tests on battery charge control, evaluating the linear sweep of current from charge to discharge, and determining zero current voltage. The fast transient behavior of batteries in the microsecond range, and the electrical conductance of nickel sinters in the thickness direction are also considered. Mechanical problems experienced in the vibration of Ni-Cd batteries and tests to simulate cyclic fatigue of the steel table connecting the plates to the comb are considered. Methods of defining the distribution of forces when cells are compressed during battery packaging are also explored.

Gross, S.↗
NASA NTRS2017

High-order Accurate Finite-Volume Formulations for the Pressure Gradient Force in Layered Ocean Models

Discretisation of the horizontal pressure gradient force in layered ocean models is a challenging task, with non-trivial interactions between the thermodynamics of the fluid and the geometry of the layers often leading to numerical difficulties. We present two new finite-volume schemes for the pressure gradient operator designed to address these issues. In each case, the horizontal acceleration is computed as an integration of the contact pressure force that acts along the perimeter of an associated momentum control-volume. A pair of new schemes are developed by exploring different control-volume geometries. Non-linearities in the underlying equation-of-state definitions and thermodynamic profiles are treated using a high-order accurate numerical integration framework, designed to preserve hydrostatic balance in a non-linear manner. Numerical experiments show that the new methods achieve high levels of consistency, maintaining hydrostatic and thermobaric equilibrium in the presence of strongly-sloping layer geometries, non-linear equations-of-state and non-uniform vertical stratification profiles. These results suggest that the new pressure gradient formulations may be appropriate for general circulation models that employ hybrid vertical coordinates andor terrain-following representations.

pressure gradient force↗
NASA NTRS2002

Computational AeroAcoustics for Fan Noise Prediction

An overview of the current state-of-the-art in computational aeroacoustics as applied to fan noise prediction at NASA Glenn is presented. Results from recent modeling efforts using three dimensional inviscid formulations in both frequency and time domains are summarized. In particular, the application of a frequency domain method, called LINFLUX, to the computation of rotor-stator interaction tone noise is reviewed and the influence of the background inviscid flow on the acoustic results is analyzed. It has been shown that the noise levels are very sensitive to the gradients of the mean flow near the surface and that the correct computation of these gradients for highly loaded airfoils is especially problematic using an inviscid formulation. The ongoing development of a finite difference time marching code that is based on a sixth order compact scheme is also reviewed. Preliminary results from the nonlinear computation of a gust-airfoil interaction model problem demonstrate the fidelity and accuracy of this approach. Spatial and temporal features of the code as well as its multi-block nature are discussed. Finally, latest results from an ongoing effort in the area of arbitrarily high order methods are reviewed and technical challenges associated with implementing correct high order boundary conditions are discussed and possible strategies for addressing these challenges ore outlined.

Envia, Ed↗
NASA NTRS1990

Development of an integrated BEM approach for hot fluid structure interaction

A comprehensive boundary element method is presented for transient thermoelastic analysis of hot section Earth-to-Orbit engine components. This time-domain formulation requires discretization of only the surface of the component, and thus provides an attractive alternative to finite element analysis for this class of problems. In addition, steep thermal gradients, which often occur near the surface, can be captured more readily since with a boundary element approach there are no shape functions to constrain the solution in the direction normal to the surface. For example, the circular disc analysis indicates the high level of accuracy that can be obtained. In fact, on the basis of reduced modeling effort and improved accuracy, it appears that the present boundary element method should be the preferred approach for general problems of transient thermoelasticity.

Dargush, G. F.↗
NASA NTRS1986

A comparison of experimental and theoretical results for leakage, pressure gradients, and rotordynamic coefficients for tapered annular gas seal

A brief review of current annular seal theory and a discussion of the predicted effect on stiffness of tapering the seal stator are presented. An outline of Nelson's analytical-computational method for determining rotordynamic coefficients for annular compressible-flow seals is included. Modifications to increase the maximum rotor speed of an existing air-seal test apparatus at Texas A&M University are described. Experimental results, including leakage, entrance-loss coefficients, pressure distributions, and normalized rotordynamic coefficients, are presented for four convergent-tapered, smooth-rotor, smooth-stator seals. A comparison of the test results shows that an inlet-to-exit clearance ratio of 1.5 to 2.0 provides the maximum direct stiffness, a clearance ratio of 2.5 provides the greatest stability, and a clearance ratio of 1.0 provides the least stability. The experimental results are compared to theoretical results from Nelson's analysis with good agreement. Test results for cross-coupled stiffness show less sensitivity of fluid prerotation than predicted.

Elrod, D. A.↗
NASA NTRS1994

Multigrid one shot methods for optimal control problems: Infinite dimensional control

The multigrid one shot method for optimal control problems, governed by elliptic systems, is introduced for the infinite dimensional control space. ln this case, the control variable is a function whose discrete representation involves_an increasing number of variables with grid refinement. The minimization algorithm uses Lagrange multipliers to calculate sensitivity gradients. A preconditioned gradient descent algorithm is accelerated by a set of coarse grids. It optimizes for different scales in the representation of the control variable on different discretization levels. An analysis which reduces the problem to the boundary is introduced. It is used to approximate the two level asymptotic convergence rate, to determine the amplitude of the minimization steps, and the choice of a high pass filter to be used when necessary. The effectiveness of the method is demonstrated on a series of test problems. The new method enables the solutions of optimal control problems at the same cost of solving the corresponding analysis problems just a few times.

Arian, Eyal↗
NASA NTRS2003

Risk of Skin Cancer from Space Radiation

We review the methods for estimating the probability of increased incidence of skin cancers from space radiation exposure, and describe some of the individual factors that may contribute to risk projection models, including skin pigment, and synergistic effects of combined ionizing and UV exposure. The steep dose gradients from trapped electrons, protons, and heavy ions radiation during EVA and limitations in EVA dosimetry are important factors for projecting skin cancer risk of astronauts. We estimate that the probability of increased skin cancer risk varies more than 10-fold for individual astronauts and that the risk of skin cancer could exceed 1 % for future lunar base operations for astronauts with light skin color and hair. Limitations in physical dosimetry in estimating the distribution of dose at the skin suggest that new biodosimetry methods be developed for responding to accidental overexposure of the skin during future space missions.

Cucinotta, Francis A.↗
NASA NTRS1992

Stress Intensity Factors for Part-Through Surface Cracks in Hollow Cylinders

Flaws resulting from improper welding and forging are usually modeled as cracks in flat plates, hollow cylinders or spheres. The stress intensity factor solutions for these crack cases are of great practical interest. This report describes some recent efforts at improving the stress intensity factor solutions for cracks in such geometries with emphasis on hollow cylinders. Specifically, two crack configurations for cylinders are documented. One is that of a surface crack in an axial plane and the other is a part-through thumb-nail crack in a circumferential plane. The case of a part-through surface crack in flat plates is used as a limiting case for very thin cylinders. A combination of the two cases for cylinders is used to derive a relation for the case of a surface crack in a sphere. Solutions were sought which cover the entire range of the geometrical parameters such as cylinder thickness, crack aspect ratio and crack depth. Both the internal and external position of the cracks are considered for cylinders and spheres. The finite element method was employed to obtain the basic solutions. Power-law form of loading was applied in the case of flat plates and axial cracks in cylinders and uniform tension and bending loads were applied in the case of circumferential (thumb-nail) cracks in cylinders. In the case of axial cracks, the results for tensile and bending loads were used as reference solutions in a weight function scheme so that the stress intensity factors could be computed for arbitrary stress gradients in the thickness direction. For circumferential cracks, since the crack front is not straight, the above technique could not be used. Hence for this case, only the tension and bending solutions are available at this time. The stress intensity factors from the finite element method were tabulated so that results for various geometric parameters such as crack depth-to-thickness ratio (a/t), crack aspect ratio (a/c) and internal radius-to-thickness ratio (R/t) or the crack length-to-width ratio (2c/W) could be obtained by interpolation and extrapolation. Such complete tables were then incorporated into the NASA/FLAGRO computer program which is widely used by the aerospace community for fracture mechanics analysis.

Mettu, Sambi R.↗
NASA NTRS2004

Optical Diagnostics for Thermal Barrier Coatings

The translucent nature of ceramic oxide thermal barrier coatings (TBCs) provides an opportunity to employ optical probes to monitor temperature gradients and buried damage propagation within the coating. An important advantage of noncontact optical diagnostics is that they are amendable to health monitoring of TBCs in service. In this paper, two optical diagnostic approaches, operating in different wavelength regimes, are discussed. The first approach is the use of mid-infrared reflectance (MIR) to monitor the progression of TBC delamination produced by thermal cycling. This approach takes advantage of the maximum transparency of the TBCs at mid-infrared wavelengths, in particular, between 3 and 5 microns. Recent progress in extending the MIR method to a more practical visual inspection tool will be presented. A second approach, using visible wavelengths, is the embedding of thermographic phosphors within the TBC to add sensing functions to the coating that can provide depth-selective information about temperature gradients and TBC integrity. Emphasis will be given to the use of fluorescence decay time measurements to provide temperature readings from a thermographic phosphor layer residing beneath the TBC.

Eldridge, J. I.↗
NASA NTRS1972

A theoretical and experimental investigation of the hypersonic turbulent boundary layer subject to normal and longitudinal pressure gradients and cross flow along a windward plane of symmetry.

A theoretical and experimental investigation for the hypersonic turbulent boundary layer undergoing both normal and longitudinal pressure gradients and cross flow along the plane of symmetry on the windward side of a body of revolution is presented. The purpose of this investigation is to provide a method for calculating boundary layers that occur on the centerline of symmetry of an inlet. The three-dimensional compressible integral equations are written along the symmetry plane and integrated numerically, and the effect of the cross flow on the behavior of the boundary layer is also studied. The experiments were performed at Mach number 6 and Reynolds number of 44,000,000/ft for the axisymmetric case and for the case with cross flow. The velocity profiles and the boundary layer parameters agree relatively well with theoretical values with and without cross flow effects. The presence of the cross flow alters significantly the value of the profiles. Total temperature-velocity correlations deviate from the linear Crocco distribution. Heat transfer rates agree with theoretical values and values obtained by the flat plate reference enthalpy method. Cross flow influence on the heat transfer rates is of minor relevance.

Zakkay, V.↗
NASA NTRS1976

Studies of noble-metal thermocouple stability at high temperatures

Two investigatory studies on performance characteristics of noble-metal thermocouples are described. (1) thermoelectric stability as affected by preferential oxidation of iridium in the system iridium-40% rhodium versus iridium, and (2) the effects of temperature gradients on the emf stability of the systems platinum-13% rhodium versus platinum and iridium-40% rhodium versus iridium, operating in air. The stability investigation was carried out at three temperatures - 1700, 1850, and 2000 C - by comparing the output of the test thermocouple in air with the output of an identically constructed reference thermocouple in nitrogen. The results show that no calibration shift was observed producing a change in output greater than that corresponding to a 2.0% change in the indicated temperature for all samples tested. The investigation of gradient effects was carried out by subjecting test thermocouples to both severe and mild gradients for periods up to 200 hours. For the platinum system, the operating temperature was 1500 C with gradients of 1475 and 700 C/cm; for the iridium system, 2000 C with gradients of 700, 1500, and 1975 C/cm. Exposure to temperature gradients was found to introduce significant changes in calibration for both systems. In both investigations, the thermoelements were examined by means of electron-probe analysis and by metallographic methods to detect chemical and structural changes. Data and micrographs are presented.

Freeze, P.↗
NASA NTRS1977

A numerical investigation of the finite element method in compressible primitive variable Navier-Stokes flow

The results of a comprehensive numerical investigation of the basic capabilities of the finite element method (FEM) for numerical solution of compressible flow problems governed by the two-dimensional and axis-symmetric Navier-Stokes equations in primitive variables are presented. The strong and weak points of the method as a tool for computational fluid dynamics are considered. The relation of the linear element finite element method to finite difference methods (FDM) is explored. The calculation of free shear layer and separated flows over aircraft boattail afterbodies with plume simulators indicate the strongest assets of the method are its capabilities for reliable and accurate calculation employing variable grids which readily approximate complex geometry and capably adapt to the presence of diverse regions of large solution gradients without the necessity of domain transformation.

Cook, C. H.↗
NASA NTRS1978

Aeroelastic stability and response of horizontal axis wind turbine blades

The coupled flap-lag-torsion equations of motion of an isolated horizontal axis wind turbine blade are formulated. Quasi-steady blade-element strip theory was applied to derive the aerodynamic operator which includes boundary layer type gradient winds. The final equations which have periodic coefficients were solved in order to obtain the aeroelastic response and stability of large horizontal axis wind turbine blade. A new method of generating an appropriate time-dependent equilibrium position (required for the stability analysis) has been implemented. Representative steady-state responses and stability boundaries, applicable mainly to an existing blade design (NASA/-ERDA MOD-0), are presented. The results indicate that the MOD-0 configuration is a basically stable design and that blade stability is not sensitive to offsets between blade elastic axis and aerodynamic center. Blade stability appears to be sensitive to precone. The tower shadow (or wake) has a considerable effect on the flap response but leaves blade stability unchanged. Finally, it was found that non linear terms in the equations of motion can significantly affect the linearized stability boundaries, however, these terms have a negligible effect on blade response at operating conditions.

Kottapalli, S. B. R.↗
NASA NTRS1982

Interactive initialization of heat flux parameters for numerical models using satellite temperature measurements

A method for obtaining patterns of moisture availability (and net evaporation) from satellite infrared measurements employs Carlson's boundary layer model and a variety of image processing routines executed by a minicomputer. To test the method with regard to regional scale moisture analyses, two case studies were chosen because of the availability of HCMM data and because of the presence of a large horizontal gradient in antecedent precipitation and crp moisture index. Results show some correlation in both cases between antecedent precipitation and derived moisture availability. Apparently, regional-scale moisture availability patterns can be determined with some degree of fidelity but the values themselves may be useful only in the relative sense and significant to within plus or minus one category of dryness over a range of 4 or 5 categories between absolutely dry and field saturation. Preliminary results suggest that the derived moisture values correlate best with longer-term precipitation totals, suggesting that the infrared temperatures respond more sensitively to a relatively deep substrate layer.

Carlson, T. N.↗
NASA NTRS1985

Computational aspects of unsteady flows

The calculation of unsteady flows and the development of numerical methods for solving unsteady boundary layer equations and their application to the flows around important configurations such as oscillating airfoils are presented. A brief review of recent work is provided with emphasis on the need for numerical methods which can overcome possible problems associated with flow reversal and separation. The zig-zag and characteristic box schemes are described in this context, and when embodied in a method which permits interaction between solutions of inviscid and viscous equations, the characteristic box scheme is shown to avoid the singularity associated with boundary layer equations and prescribed pressure gradient. Calculations were performed for a cylinder started impulsively from rest and oscillating airfoils. The results are presented and discussed. It is conlcuded that turbulence models based on an algebraic specification of eddy viscosity can be adequate, that location of translation is important to the calculation of the location of flow separation and, therefore, to the overall lift of an oscillating airfoil.

Cebeci, T.↗