Search NASA⌕ Search

SEARCH · Search NASA

Results for “Relaxation method”

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 127 records · Page 7

Theoretical effect of modifications to the upper surface of two NACA airfoils using smooth polynomial additional thickness distributions which emphasize leading edge profile and which vary linearly at the trailing edge

An investigation was conducted on a CDC 7600 digital computer to determine the effects of additional thickness distributions to the upper surface of airfoils. The additional thickness distribution had the form of a continuous mathematical function which disappears at both the leading edge and the trailing edge. Results were obtained at a Mach number of 0.2 with an angle of attack of 6 deg. All calculations employed the full potential flow equations for two dimensional flow. The relaxation method of Jameson was used for solution of the potential flow equations. It is shown that increasing the thickness and variations in shape increases the lift and the adverse pitching moment coefficients.

Hague, D. S.↗

Computer simulation of transonic flow past airfoils with boundary layer correction

A computer program has been developed to solve the compressible flow equation for the velocity potential. The exterior of the airfoil is mapped onto the unit circle and the flow is computed on a grid in the circle plane. A relaxation method using backward differencing in the flow direction at supersonic points permits solutions for large supersonic areas. The pressure distribution resulting from the flow becomes the input to the von Karman momentum equation which when integrated gives the displacement thickness. This displacement thickness is smoothed and added to the airfoil to account for the turbulent boundary layer. The boundary layer correction is computed iteratively with the flow. Results from this program and test data agree well.

Bauer, F.↗

Time-dependent transonic flow solutions for axial turbomachinery

Three-dimensional unsteady transonic flow through an axial turbomachine stage is described in terms of a pair of two-dimensional formulations pertaining to orthogonal surfaces, namely, a blade-to-blade surface and a hub-to-casing surface. The resulting systems of nonlinear, inviscid, compressible equations of motion are solved by an explicit finite-difference technique. The blade-to-blade program includes the periodic interaction between rotor and stator blade rows. Treatment of the boundary conditions and of the blade slipstream motion by a characteristic type procedure is discussed in detail. Harmonic analysis of the acoustic far field produced by the blade row interaction, including an arbitrary initial transient, is outlined. Results from the blade-to-blade program are compared with experimental measurements of the rotating pressure field at the tip of a high-speed fan. The hub-to-casing program determines circumferentially averaged flow properties on a meridional plane. Blade row interactions are neglected in this formulation, but the force distributions over the entire blade surface for both the rotor and stator are obtained. Results from the hub-to-casing program are compared with a relaxation method solution for a subsonic rotor. Results are also presented for a quiet fan stage which includes transonic flow in both the rotor and stator and a normal shock in the stator.

Erdos, J.↗

Approximate method for calculating transonic flow about lifting wing-body configurations

The three-dimensional problem of transonic flow about lifting wing-body configurations is reduced to a two-variable computational problem with the method of matched asymptotic expansions. The computational problem is solved with the method of relaxation. The method accounts for leading-edge separation, the presence of shock waves, and the presence of solid, slotted, or porous tunnel walls. The Mach number range of the method extends from zero to the supersonic value at which the wing leading edge becomes sonic. A modified form of the transonic area rule which accounts for the effect of lift is developed. This effect is explained from simple physical considerations.

Barnwell, R. W.↗

Numerical studies of three-dimensional breakdown in trailing vortex wakes

The development of dangerous vortex wakes trailing between aircraft for appreciable distance downstream, and posing hazards to following aircraft, is calculated using finite-element three-dimensional relaxation methods. Near-wake, middle-wake, and far-wake conditions are delineated. Vortex system finite lifetimes can be calculated with greater facility owing to the inclusion of a self-induction term in the solution dependent upon local curvature and vortex core radius. Looping and convergence phenomena are studied for single vortex pairs, twin pairs, and multiple pairs. Proportions of flap span and wing span are shown to affect the time to convergence of vortex pairs. The self-induction term renders the convergence time of vortices a sensitive function of wavelength. Wakes trailing between wide-body ('Jumbo') aircraft are also examined.

Hackett, J. E.↗

Numerical Studies of Three-dimensional Breakdown in Trailing Vortex Wakes

Finite element, three dimensional relaxation methods are used to calculate the development of vortex wakes behind aircraft for a considerable downstream distance. The inclusion of a self-induction term in the solution, dependent upon local curvature and vortex core radius, permits calculation of finite lifetimes for systems for which infinite life would be predicted two dimensionally. The associated computer program is described together with single-pair, twin-pair, and multiple-pair studies carried out using it. It is found, in single-pair studies, that there is a lower limit to the wavelengths at which the Crow-type of instability can occur. Below this limit, self-induction effects cause the plane of the disturbance waves to rotate counter to the vortex direction. Self induction in two dimensionally generated twin spiral waves causes an increase in axial length which becomes more marked with decreasing initial wavelength. The time taken for vortex convergence toward the center plane is correspondingly increased. The limited parametric twin-pair study performed suggests that time-to-converge increases with increasing flap span. Limited studies of Boeing 747 configurations show correct qualitative response to removal of the outer flap and to gear deployment, as compared with wind tunnel and flight test experience.

Evans, P. F.↗

Transonic flow past swept wings

FLO-22 computer program aids in numerical analysis of transonic potential flow past lifting, swept wing. FLO-22 uses relaxation method to solve finite-difference approximation of full-potential equation for transonic flow past configuration consisting of wing or arbitrary planform and dihedral extending from symmetry plane or wall. Comparisons of FLO-22 calculated results with experimental data for both conventional and super-critical transport wings show good agreement.

Caughey, D. A.↗

Transonic flow calculations

The development of relaxation methods to calculate transonic flows is discussed. Rather accurate predictions can be made for a number of flows of interest using the transonic potential flow equation, which may be derived from the Euler equations for inviscid compressible flow by introducing the assumption that the flow is irrotational. Once the choice of a mathematical model has been settled, the numerical procedure for actually computing a solution contains two main elements: the construction of a discrete approximation which converges to the solution of the continuous problem in the limit as the mesh width is reduced to zero, and the solution of the resulting set of nonlinear difference equations by a convergent iterative scheme. The choice of an appropriate coordinate system and its influence on the accuracy of the discrete approximation is discussed. Several applications of the general method are described.

Jameson, A.↗

Inviscid transonic flow over axisymmetric bodies

Axisymmetric transonic flow is of interest not only because of its practical application to missile and launch vehicle aerodynamics but also because of its relation, in terms of area rule, to fully three dimensional flow. RAXBOD computer program analyzes steady, inviscid, irrotational, transonic flow over axisymmetric bodies in free air. RAXBOD uses finite-difference relaxation method to solve numerically exact formulation of disturbance velocity potential with exact surface boundary conditions. Agreement with available experimental results has been good in cases where viscous effects and wind-tunnel wall interference are not important.

South, J. C., Jr.↗

Stellar convection 2: A multi-mode numerical solution for convection in spheres

The convective flow of a self gravitating sphere of Boussinesq fluid for small Reynolds and Peclet numbers is numerically determined. The decomposition of the equations of motion into modes is reviewed and a relaxation method is developed and presented to compute the solutions to these equations. The stable equilibrium flow for a Rayleigh number of 10 to the 4th power and a Prandtl number of 10 is determined. The 2 and 3 dimensional spectra of the kinetic and thermal energies and the convective flux as a function of wavelengths are calculated in terms of modes. The anisotropy of the flow as a function of wavelength is defined.

Marcus, P. S.↗

Optimal large angle maneuvers with simultaneous shape control/vibration arrest

A relaxation method is demonstrated which reliably solves the nonlinear two point boundary value problem which arises when optimal control theory is applied to determination of large angle maneuvers of flexible spacecraft. The basic ideas are summarized and several idealized maneuvers are determined. The emphasis is upon demonstrating the basic ideas and practical aspects of the methodology.

Turner, J. D.↗

Nonequilibrium material effects on the behavior of polymeric composite matrices and their related composites

The effects of physical aging on the material properties of some linear and network macromolecular glasses are discussed. The free volume concept is used to describe this behavior. The effect of physical aging on properties of some uniaxial graphite/fiber epoxy resin composites is investigated using stress relaxation in both tensile and flexural modes. The matrix polymers used were resins both of which are based on a 4,4-methylenedianiline derivative of epichlorohydrin with diamino diphenyl sulfone (DDS) as the curing agent. The matrix resin, as used in the practical application in composites, not fully cured and the glass transition of the network was dependent on the curing schedule. The physical aging of the bulk crosslinked epoxy was found to depend on the annealing temperature, and the T sub g of the resin. The physical aging of the composite, monitored by the stress relaxation method, was found to be dependent on the testing direction.

Wilkes, G. L.↗

The global distribution of water vapor in the middle atmosphere of Venus

Near-IR measurements are presented of the mean vertical and horizontal distribution of water vapor in the Venus clouds as measured by the Pioneer Venus Orbiter IR radiometer, and comparisons are made with previous data. Six thermal channels were used to generate several hundred thousand readings for determination of the mean mixing ratio. Averaging was performed as a function of the solar zenith angle, with profiles retrieved with a relaxation method applied to radiance data at 45 microns. Consideration was given to mean cloud models and temperature profiles obtained from the five temperature sounding channels scanning from 11.5-15 microns. Laboratory tests were effected to validate the transmission functions. The results included a maximum column abundances above the cloud optical depth in the early afternoon in the equatorial regions. Mixing ratio enhancement was highest on the dayside and at high altitudes, with a mean ratio of 0.0001 at a 40% uncertainty level. Day-to-day fluctuations in the pressure level at 11.5 microns was larger than 10%, far below the factors of 2-3 determined by other investigators.

Schofield, J. T.↗

Error analysis and correction of discrete solutions from finite element codes

Many structures are an assembly of individual shell components. Therefore, results for stresses and deflections from finite element solutions for each shell component should agree with the equations of shell theory. This paper examines the problem of applying shell theory to the error analysis and the correction of finite element results. The general approach to error analysis and correction is discussed first. Relaxation methods are suggested as one approach to correcting finite element results for all or parts of shell structures. Next, the problem of error analysis of plate structures is examined in more detail. The method of successive approximations is adapted to take discrete finite element solutions and to generate continuous approximate solutions for postbuckled plates. Preliminary numerical results are included.

Thurston, G. A.↗

Computing with high-resolution upwind schemes for hyperbolic equations

Computational aspects of modern high-resolution upwind finite-difference schemes for hyperbolic systems of conservation laws are examined. An operational unification is demonstrated for constructing a wide class of flux-difference-split and flux-split schemes based on the design principles underlying total variation diminishing (TVD) schemes. Consideration is also given to TVD scheme design by preprocessing, the extension of preprocessing and postprocessing approaches to general control volumes, the removal of expansion shocks and 'glitches', relaxation methods for implicit TVD schemes, and a new family of high-accuracy TVD schemes.

Chakravarthy, S. R.↗

Crack Growth in Single-Crystal Silicon

Report describes experiments on crack growth in single-crystal silicon at room temperature in air. Crack growth in (111) cleavage plane of wafers, 50 by 100 by 0.76 mm in dimension, cut from Czochralski singlecrystal silicon studied by double-torsion load-relaxation method and by acoustic-emission measurements. Scanning electron microscopy and X-ray topography also employed. Results aid in design and fabrication of silicon photovoltaic and microelectronic devices.

Chen, C. P.↗

Improvement of the process for electrochemical impregnation of nickel hydroxide electrodes

Nickel hydroxide electrodes containing 11g/dsqm hydroxide, with capacities of 3.6 to 3.8 Ah/dsqm were prepared at 353 K by electrochemical impregnation. The reproducibility of the results is obtained by readjusting the pH before each preparation. The control of each electrode is done during two cycles of charge and discharge following the manufacture by a potential relaxation method.

Comtat, M.↗