Engineering PapersSearch

SEARCH · Engineering Papers

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

Determination of the temperature profile in an atmosphere from its outgoing radiation

A highly convergent relaxation method for the inversion of the full radiative transfer equation was developed. The results of the iterative solution indicate that convergence can be achieved over a wide range of initial guesses, enabling the temperature profile of a relatively unknown atmosphere to be unambiguously determined. The method is illustrated by examples for the outgoing radiance in the earth's atmosphere for the region of the 4.3 microns CO2 band, but can be similarly applied in other frequency ranges.

Chahine, M. T.

Energy distribution among reaction products. V.

Discussion of three reactions, one point of theoretical interest being the predicted correlation between barrier height and barrier location. The H + Br 2 reaction having a lower activation barrier than H + Cl 2, should have an earlier barrier, and hence a greater percentage attractive energy release and higher efficiency of vibrational excitation. Information is developed concerning the effect of isotopic substitution in the pair of reactions H + Cl 2 and D + Cl 2. The 'arrested relaxation' method was used. Essentially, the method involves reacting two diffuse reagent beams in a reaction vessel with background pressure less than 0.001 torr, and with walls cooled by liquid nitrogen or liquid helium.

Anlauf, K. G.

High-Q toroidal cavities for high frequency klystrons.

A toroidal cavity developed for a 4-KW 12 GHz satellite-borne television transmitter klystron is described. The cavity has an internal Q 40% higher than that of a conventional cylindrical doubly reentrant cavity, thus yielding higher circuit efficiency and conserving energy which cannot be recovered in multistage depressed potential beam collectors. As a result of optimization studies with a digital computer program for obtaining cavity field distributions by a relaxation method and for computing the intrinsic cavity parameters, a particular cavity configuration with conical reentrant tunnel tips and toroidal walls is shown to provide good thermal characteristics and mechanical rigidity as well as low internal losses.

Branch, G. M.

Relaxation solutions for inviscid axisymmetric transonic flow over blunt or pointed bodies.

A finite-difference relaxation method is presented for numerical solution of the full potential equation and exact boundary conditions for general axisymmetric bodies is inviscid, steady transonic flow. Body-normal coordinates are used in the nose region and sheared cylindrical coordinates are used on the afterbody to accommodate corners such as boattails and flares. An improved difference scheme is used which does not require that the flow be nearly alined with a coordinate direction in supersonic regions, and which treats either subsonic or supersonic free streams. Numerical results are illustrated for some simple classical shapes such as spheres and ellipsoids, and for more practical shapes like tangent-ogives with boattails. Special attention is given to bodies which have been studied for area-rule applications. Agreement with available experimental results is good in cases where viscous effects and wind-tunnel wall interference are not important.

South, J. C., Jr.

Calculation of separated flows at subsonic and transonic speeds

A boundary-layer integral approach is combined with a finite-difference relaxation method to calculate viscous interactions between separated flows at subsonic and transonic velocities. Results are obtained for separated laminar flows on circular-arc airfoils at zero angle of attack and are compared with data of Collins (1972). Inviscid and viscous flows are covered.

Klineberg, J. M.

Remote sounding of cloudy atmospheres. I - The single cloud layer

The relaxation method for the inverse solution of the radiative transfer equation is applied in a dual-frequency scheme for the determination of complete vertical temperature profiles in cloudy atmospheres from radiance observations alone, without any additional information related to the expected solutions. The dual-frequency principle employs to advantage a property in the Planck function of the dependence of intensity on frequency. This property leads to the formulation of a new convergence criterion for the selection of cloud-sounding frequencies to be used for reconstructing the clear column radiance from observations made in the presence of a broken cloud layer in all fields of view. The principle is applied to the case of observations in two adjacent or partially overlapping fields of view and to the case of observations in a single field of view. The solutions are illustrated by numerical examples in the dual-frequency ranges of the 4.3 and 15-micron CO2 bands of the terrestrial atmosphere.

Chahine, M. T.

An investigation on the effect of second-order additional thickness distributions to the upper surface of an NACA 64 sub 1-212 airfoil

An investigation was conducted on a CDC 7600 digital computer to determine the effects of additional thickness distributions to the upper surface of an NACA 64 sub 1 - 212 airfoil. Additional thickness distributions employed were in the form of two second-order polynomial arcs which have a specified thickness at a given chordwise location. The forward arc disappears at the airfoil leading edge, the aft arc disappears at the airfoil trailing edge. At the juncture of the two arcs, x = x, continuity of slope is maintained. The effect of varying the maximum additional thickness and its chordwise location on airfoil lift coefficient, pitching moment, and pressure distribution was investigated. Results were obtained at a Mach number of 0.2 with an angle-of-attack of 6 degrees on the basic NACA 64 sub 1 - 212 airfoil, and all calculations employ the full potential flow equations for two dimensional flow. The relaxation method of Jameson was employed for solution of the potential flow equations.

Hague, D. S.

An investigation on the effect of second-order additional thickness distributions to the upper surface of an NACA 64-206 airfoil

An investigation was conducted on a CDC 7600 digital computer to determine the effects of additional thickness distributions to the upper surface of an NACA 64-206 airfoil. Additional thickness distributions employed were in the form of two second-order polynomial arcs which have a specified thickness at a given chordwise location. The forward arc disappears at the airfoil leading edge, the aft arc disappears at the airfoil trailing edge. At the juncture of the two arcs, x = x, continuity of slope is maintained. The effect of varying the maximum additional thickness and its chordwise location on airfoil lift coefficient, pitching moment, and pressure distribution was investigated. Results were obtained at a Mach number of 0.2 with an angle-of-attack of 6 degrees on the basic NACA 64-206 airfoil, and all calculations employ the full potential flow equations for two dimensional flow. The relaxation method of Jameson was employed for solution of the potential flow equations.

Merz, A. W.

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 quadratically 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 the NACA 64-206 and 64 sub 1 - 212 airfoils. The additional thickness distribution had the form of a continuous mathematical function which disappears at both the leading edge and the trailing edge. The function behaves as a polynomial of order epsilon sub 1 at the leading edge, and a polynomial of order epsilon sub 2 at the trailing edge. Epsilon sub 2 is a constant and epsilon sub 1 is varied over a range of practical interest. The magnitude of the additional thickness, y, is a second input parameter, and the effect of varying epsilon sub 1 and y on the aerodynamic performance of the airfoil was investigated. Results were obtained at a Mach number of 0.2 with an angle-of-attack of 6 degrees on the basic airfoils, and all calculations employ the full potential flow equations for two dimensional flow. The relaxation method of Jameson was employed for solution of the potential flow equations.

Merz, A. W.

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.