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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

NASA NTRS1993

A computer-controlled apparatus for Seebeck inhomogeneity testing of sheathed thermocouples

Mineral-insulated metal-sheathed (MIMS) thermocouple assemblies are used throughout industry and research facilities as a method of temperature measurement where requirements for either harsh environmental conditions exist, or where rigidity of the measurement probe is required. Seebeck inhomogeneity is the abnormal variation of the Seebeck coefficient from point to point in a material. It is not disclosed in conventional calibration. A standardized method of measuring thermoelectric inhomogeneity along the thermocouple probe length is not available. Therefore, calibration for sheathed probes normally does not include testing of probe inhomogeneity. The measurement accuracy would be severely impacted if significant inhomogeneity and a temperature gradient were present in the same region of the probe. A computer-controlled system for determining inhomogeneities was designed, fabricated, and tested. This system provides an accurate method for the identification of the location of inhomogeneity along the length of a sheathed thermocouple and for the quantification of the inhomogeneity. This paper will discuss the apparatus and procedure used to perform these tests and will present data showing tests performed on sheathed thermocouple probes.

Burkett, Cecil G., Jr.↗

NASA NTRS1985

A Residuals Approach to Filtering, Smoothing and Identification for Static Distributed Systems

An approach for state estimation and identification of spatially distributed parameters embedded in static distributed (elliptic) system models is advanced. The method of maximum likelihood is used to find parameter values that maximize a likelihood functional for the system model, or equivalently, that minimize the negative logarithm of this functional. To find the minimum, a Newton-Raphson search is conducted that from an initial estimate generates a convergent sequence of parameter estimates. For simplicity, a Gauss-Markov approach is used to approximate the Hessian in terms of products of first derivatives. The gradient and approximate Hessian are computed by first arranging the negative log likelihood functional into a form based on the square root factorization of the predicted covariance of the measurement process. The resulting data processing approach, referred to here by the new term of predicted data covariance square root filtering, makes the gradient and approximate Hessian calculations very simple. A closely related set of state estimates is also produced by the maximum likelihood method: smoothed estimates that are optimal in a conditional mean sense and filtered estimates that emerge from the predicted data covariance square root filter.

Rodriguez, G.↗

NASA NTRS1991

A nonrecursive 'Order N' preconditioned conjugate gradient/range space formulation of MDOF dynamics

This paper addresses the requirements of present-day mechanical system simulations of algorithms that induce parallelism on a fine scale and of transient simulation methods which must be automatically load balancing for a wide collection of system topologies and hardware configurations. To this end, a combination range space/preconditioned conjugage gradient formulation of multidegree-of-freedon dynamics is developed, which, by employing regular ordering of the system connectivity graph, makes it possible to derive an extremely efficient preconditioner from the range space metric (as opposed to the system coefficient matrix). Because of the effectiveness of the preconditioner, the method can achieve performance rates that depend linearly on the number of substructures. The method, termed 'Order N' does not require the assembly of system mass or stiffness matrices, and is therefore amenable to implementation on work stations. Using this method, a 13-substructure model of the Space Station was constructed.

Kurdila, A. J.↗

NASA NTRS1992

Adaptive mesh strategies for the spectral element method

An adaptive spectral method was developed for the efficient solution of time dependent partial differential equations. Adaptive mesh strategies that include resolution refinement and coarsening by three different methods are illustrated on solutions to the 1-D viscous Burger equation and the 2-D Navier-Stokes equations for driven flow in a cavity. Sharp gradients, singularities, and regions of poor resolution are resolved optimally as they develop in time using error estimators which indicate the choice of refinement to be used. The adaptive formulation presents significant increases in efficiency, flexibility, and general capabilities for high order spectral methods.

Mavriplis, Catherine↗

NASA NTRS1988

Investigation of the influence of wind shear on the aerodynamic characteristics of aircraft using a vortex-lattice method

The objective was to investigate and characterize the aerodynamic effect of shear flow through a series of sensitivity studies of the wind velocity gradients and wing planform geometry parameters. The wind shear effect was computed using a modified vortex-lattice computer program and characterized through the formulation of wind shear aerodynamic coefficients. The magnitude of the aerodynamic effect was demonstrated by computing the resultant change in the aerodynamics of a conventional wing and tail combination on a fixed flight path through a simulated microburst. The results of the study indicate that a significant amount of the control authority of an airplane may be required to counteract the wind shear induced forces and moments in the microburst environment.

Vicroy, Dan D.↗

NASA NTRS1983

An integral analysis of transonic normal shock wave/turbulent boundary layer interactions in internal flow

An approximate integral viscous-inviscid interaction method is presented for calculating the development of a turbulent boundary layer subjected to a normal shock wave induced adverse pressure gradient in an internal axisymmetric flow. The inflow conditions and the downstream pressure are provided for the computation. In the supersonic region of shock pressure rise, the Prandtl-Meyer function is used to couple the viscous and inviscid flows. An analytical model for the coupling process is postulated and appropriate equations are defined. Downstream of the sonic point, one-dimensional inviscid flow is assumed for coupling with the viscous flow. The turbulent boundary layer is calculated using Green's integral lag-entrainment method. Comparisons of the solutions with the experimental data are made for interactions which are unseparated, near separation and separated. For comparison purposes, solutions to the time-dependent, mass-averaged, Navier-Stokes equations incorporating a two-equation, Wilcox-Rubesin turbulence model are also shown. The computed results from the integral method show good agreement with experimental data for unseparated interactions and reasonable agreement with the trend of the viscous effects when the interaction becomes increasingly separated.

Om, D.↗

NASA NTRS1950

Approximate Turbulent Boundary-Layer Development in Plane Compressible Flow Along Thermally Insulated Surfaces with Application to Supersonic-Tunnel Contour Correction

Numerical solutions of the differential equation obtained from the momentum theorem for the development of a turbulent boundary layer along a thermally insulated surface in two-dimensional and in radial shock-free flow are presented in tabular form for a range of Mach numbers from 0.100 to 10. The solution can be used in a step-wise procedure with any given distribution of favorable pressure gradients and for zero pressure gradients. Solutions are also given for use with moderate adverse pressure gradients. The mean velocity in the boundary layer is approximated by a power-law profile. In view of the stepwise integration methods to be used, the exponent designated the profile shape can be varied along the surface between the integral fraction limits 1/5 and 1/11 through interpolation. Agreement obtained between theoretical and experimental boundary-layer development in a supersonic nozzle at a nominal Mach number of 2 indicates the general validity of the approximations used in the analysis - in particular, the method of extrapolating low-speed skin-friction relations to high Mach number flows. The extrapolation method used assumes that the skin-friction coefficient depend primarily on Reynolds number, provided that the density and the kinematic viscosity are evaluated at surface conditions.

BOUNDARY LAYER - WIND TUNNEL WALLS↗

NASA NTRS1977

A novel concept for subsonic inlet boundary-layer control

A self-bleeding method for boundary layer control is described and tested for a subsonic inlet designed to operate in the flowfield generated by high angles of attack. Naturally occurring surface static pressure gradients are used to remove the boundary layer from a separation-prone region of the inlet and to reinject it at a less critical location with a net performance gain. The results suggest that this self-bleeding method for boundary-layer control might be successfully applied to other inlets operating at extreme aerodynamic conditions.

Miller, B. A.↗

NASA NTRS1971

An analysis and comparison of several trajectory optimization methods

The sensitivities of the convergence characteristics of the methods to initially assumed parameters and trial solution, convergence times, computer logic, and storage requirements are discussed. Numerical comparison of the convergence characteristics is made by considering a minimum time, low thrust, Earth-Mars transfer trajectory. A modified quasi-linearization method reduces convergence time by approximately 70% when compared with the generalized Newton-Raphson method and allows the terminal boundary to be specified by a general function of the problem variables. A uniquely specified and easily determined, time dependent weighting matrix for the gradient techniques accelerates the shaping of the optimal control program and improves the convergence characteristics during the terminal iterations. Convergence envelopes, indicating how sensitive the convergence characteristics are to initially assumed parameters, are plotted for the perturbation and quasi-linearization methods. Several iteration schemes are proposed which increase the size of the convergence envelopes and decrease the sensitivity of the method to initially assumed parameters.

Lewallen, J. M.↗

NASA NTRS1974

Preparative electrophoresis of living lymphocytes

Vertical liquid columns containing low molecular weight dextran density gradients can be used for preparative lymphocyte electrophoresis on earth, in simulation of 0 gravity conditions. Another method that has been tested at 1 G, is the electrophoresis of lymphocytes in a upward direction in vertical columns. By both methods up to 10 to the 7th power lymphocytes can be separated at one time in a 30 cm glass column of 8 mm inside diameter, at 12 v/cm, in 2 hours. Due to convection and sedimentation problems, the separation at 1 G is less than ideal, but it is expected that at 0 gravity electrophoresis will prove to be a uniquely powerful cell separation tool. The technical feasibility of electrophoresing inert particles at 0 G has been proven earlier, during the flight of Apollo 16.

Vanoss, C. J.↗

NASA NTRS1974

Preparative liquid column electrophoresis of T and B lymphocytes at gravity = 1

Vertical liquid columns containing low-molecular-weight dextran density gradients can be used for preparative lymphocyte electrophoresis on earth, in simulation of zero gravity conditions. Another method that has been tested at 1 g, is the electrophoresis of lymphocytes in an upward direction in vertical columns. By both methods up to 100 million lymphocytes can be separated at one time in a 30-cm glass column of 8-mm inside diameter, at 12 V/cm, in two hours. Due to convection and sedimentation problems, the separation at 1 g is less than ideal, but it is expected that at zero gravity electrophoresis will probe to be a uniquely powerful cell separation tool.

Van Oss, C. J.↗

NASA NTRS1976

An investigation of flap-lag stability of wind turbine rotors in the presence of velocity gradients and helicopter rotors in forward flight

The flap-lag equations of motion of a torsionally rigid and centrally hinged spring-restrained rigid blade are developed using quasi-steady blade-element aerodynamic theory. These equations have periodic coefficients and are applicable to both wind-turbine rotors with velocity gradients and helicopter rotors in forward flight. By solving these equations both by the Floquet-Liapunov method and by an approximate method, the range of applicability of the latter method is established. Flap-lag and pure flapping stability boundaries illustrating the effects of flow and rotor parameters are presented. Finally, some discussion on the techniques for generating the Floquet transition matrix and on the relative merits of the coordinate axes system is presented.

Kaza, K. R. V.↗

NASA NTRS1976

Calculation methods for compressible turbulent boundary layers

Calculation procedures for non-reacting compressible two- and three-dimensional turbulent boundary layers were reviewed. Integral, transformation and correlation methods, as well as finite difference solutions of the complete boundary layer equations summarized. Alternative numerical solution procedures were examined, and both mean field and mean turbulence field closure models were considered. Physics and related calculation problems peculiar to compressible turbulent boundary layers are described. A catalog of available solution procedures of the finite difference, finite element, and method of weighted residuals genre is included. Influence of compressibility, low Reynolds number, wall blowing, and pressure gradient upon mean field closure constants are reported.

Bushnell, D. M.↗

NASA NTRS1981

Transfer of line radiation in differentially expanding atmospheres. VI The plane parallel atmosphere with expanding and contracting regions

The non-LTE radiative transfer problem for a two level atom with complete redistribution over a Doppler profile is solved for a plane parallel slab (overlying a radiating photosphere) that has a velocity field which rises symmetrically from zero at either face to a central maximum. Since the velocity gradient reverses, distant layers of the slab become coupled by radiation that jumps intervening layers. The Feautrier method is used, but an iterative variant is also employed as a check in cases where poorly conditioned matrices are encountered. Approximations are developed to explain some of the principal features. It is found that the source function S tends to have two plateaus with values near 2/3 I sub 0 and 1/3 I sub 0, where I sub 0 is the photospheric continuum incident from below; the larger value lies nearer the photosphere. The upper layers sometimes exhibit a rise in S owing to interconnection by radiation to the base. It is noted that the radiation force is largest at the two faces and the midplane. Some line profiles are found to have unusually steep absorptions at rest frequency because of the low excitation in the uppermost, stationary layers.

Noerdlinger, P. D.↗

NASA NTRS1981

Optical observations of thermospheric dynamics at Arecibo

A description is given of wind and temperature data obtained with a 15-cm Fabry-Perot interferometer. The observations were made between February and August 1980 and confirm the details of the conceptual model of Herrero and Meriwether (1980). Attention is given to the method used to analyze the neutral wind field, the zonal winds, the meridional wind, the horizontal wind gradient, the vertical wind velocity, the full wind field, and the neutral temperature. The most regular feature of the neutral wind at Arecibo is the zonal component, which has a peak velocity in the eastward direction near sunset, and decreases in velocity as the night progresses. Both the wind and temperature measurements reported provide ample confirmation of the importance of the equatorial midnight pressure bulge in modifying the dynamic behavior of the nighttime thermosphere at Arecibo.

Burnside, R. G.↗

NASA NTRS1986

An explicit finite-volume time-marching procedure for turbulent flow calculations

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite-volume, time-marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier-Stokes equations. The entire calculation is performed in the physical domain. Control volumes are chosen so that smoothing of flow properties, typically required for stability, is not required. Different time steps are used in the different governing equations. A new pressure interpolation scheme is introduced which improves the shock capturing ability of the method. A multi-volume method for pressure changes in the boundary layer allows calculations which use very long and thin control volumes (length/height - 1000). The method is compared with two test cases. Essentially incompressible turbulent boundary layer flow in an adverse pressure gradient is calculated and the computed distributions of mean velocity and shear are in good agreement with the measurements. Transonic viscous flow in a converging diverging nozzle is calculated; the Mach number upstream of the shock is approximately 1.25. The agreement between the calculated and measured shock strength and total pressure losses is good.

Nicholson, Stephen↗

NASA NTRS1987

An explicit finite-volume time-marching procedure for turbulent flow calculations

Nicholson, Stephen↗

NASA NTRS1992

A study of equation solvers for linear and non-linear finite element analysis on parallel processing computers

Concurrent computing environments provide the means to achieve very high performance for finite element analysis of systems, provided the algorithms take advantage of multiple processors. The authors have examined several algorithms for both linear and nonlinear finite element analysis. The performance of these algorithms on an Alliant FX/80 parallel supercomputer has been studied. For single load case linear analysis, the optimal solution algorithm is strongly problem dependent. For multiple load cases or nonlinear analysis through a modified Newton-Raphson method, decomposition algorithms are shown to have a decided advantage over element-by-element preconditioned conjugate gradient algorithms.

Watson, Brian C.↗