SEARCH · Engineering Papers
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.
The operation of the BGRID computer program is described for generating block-structured grids. Examples are provided to illustrate the code input and output. The application of a fully implicit AF (approximation factorization)-based computer code, called TWINGB (Transonic WING), for solving the 3D transonic full potential equation in conservation form on block-structured grids is also discussed.
An evaluation of the transonic-wing-analysis computer code TWING is presented. TWING utilizes a fully implicit, approximate-factorization iteration scheme to solve the full-potential equation in conservative form. A numerical elliptic-solver grid-generation scheme is used to generate the required finite-difference mesh. Several wing configurations have been analyzed, and comparisons of computed results have been made with available experimental data. Results indicate that the code is robust, accurate (when significant viscous effects are not present), and efficient. TWING generally produces solutions an order of magnitude faster than other conservative, full-potential codes using successive-line overrelaxation. The present method is applicable to a wide range of isolated wing configurations, including high-aspect-ratio transport wings and low-aspect-ratio, high-sweep, fighter configurations.
Engquist and Osher (1980) have introduced a finite difference scheme for solving the transonic small disturbance equation, taking into account cases in which only compression shocks are admitted. Osher et al. (1983) studied a class of schemes for the full potential equation. It is proved that these schemes satisfy a new discrete 'entropy inequality' which rules out expansion shocks. However, the conducted analysis is restricted to steady two-dimensional flows. The present investigation is concerned with the adoption of a heuristic approach. The full potential equation in conservation form is solved with the aid of a modified artificial density method, based on flux biasing. It is shown that, with the current scheme, expansion shocks are not possible.
Self-consistent transonic multifluid solutions for the dynamical and thermal structure of the two-phase dusty gas atmosphere of a comet containing a dirty H2O-clathrate ice nucleus have been obtained. This was done by solving a simultaneous set of differential equations representing conservation of number density, momentum, and energy, as well as the transfer of solar radiation responsible for the major photolytic processes and nuclear heating. The results presented include velocity and Mach number profiles of heavy species and the velocity of dust within the subsonic and supersonic regions, temperature profiles of the dust and gas in the supersonic region, and radial profiles of electron and ion densities in the collision-dominated region.
NASCRIN analyzes two-dimensional flow fields in supersonic combustion ramjet (scramjet) inlets. Solves two-dimensional Euler or Navier-Strokes equations in conservative form by unsplit, explicit, two-step finitedifference method. More recent explicit/implicit, two-step scheme also incorporated for viscous flow analysis. Algebraic, two-layer eddy-viscosity model used for turbulent flow calculations.
A computer program NASCRIN has been developed for analyzing two-dimensional flow fields in high-speed inlets. It solves the two-dimensional Euler or Navier-Stokes equations in conservation form by an explicit, two-step finite-difference method. An explicit-implicit method can also be used at the user's discretion for viscous flow calculations. For turbulent flow, an algebraic, two-layer eddy-viscosity model is used. The code is operational on the CDC CYBER 203 computer system and is highly vectorized to take full advantage of the vector-processing capability of the system. It is highly user oriented and is structured in such a way that for most supersonic flow problems, the user has to make only a few changes. Although the code is primarily written for supersonic internal flow, it can be used with suitable changes in the boundary conditions for a variety of other problems.
During the last eight years, the subject of computational unsteady transonic aerodynamics has undergone a period of rapid growth. A brief survey of some of the fundamental advances during this period is provided, and some of the research done at NASA Ames Research Center is described in detail. The small-distance potential equation is considered, taking into account features of the unsteady transonic flow, the governing equation, the pressure coefficient, the wake condition, the airfoil tangency condition, the downstream boundary condition, the alternating direction implicit (ADI) algorithm, a grid system, boundary conditions, an airfoil with oscillating trailing edge flap, three types of shock-wave motion, and an example of a simple aeroelastic problem. A description is given of an implicit algorithm for the unsteady full potential equation in conservation form.
Explore the source record for details and available documents.
NASCRIN program analyzes two-dimensional flow fields in supersoniccombustion ramjet (scramjet) inlets. Solves two-dimensional Euler or Navier-Stokes equations in conservative form by unsplit, explicit, two-step finite-difference method. More recent explicit/implicit, two-step scheme incorporated by analysis of viscous flow. Algebraic, two-layer eddy-viscosity model used for turbulent-flow calculations. Vectorized version, written for CDC CYBER 205, whereas scalar version, can be run on CRAY or other scalar computers.
Supersonic inlet flows with mixed external-internal compressions were computed using a combined implicit-explicit (Beam-Warming-Steger/MacCormack) method for solving the three-dimensional unsteady, compressible Navier-Stokes equations in conservation form. Numerical calculations were made of various flows related to such inlet operations as the shock-wave intersections, subsonic spillage around the cowl lip, and inlet started versus unstarted conditions. Some of the computed results were compared with wind tunnel data.
Supersonic inlet flows with mixed external-internal compressions of an axisymmetric inlet model were computed using a combined implicit-explicit (Beam-Warming-Steger/MacCormack) method for solving the three-dimensional unsteady, compressible Navier-Stokes equations in conservation form. Numerical calculations were made of various flows typically found in supersonic inlets such as shock-wave intersections, flow spillage around the cowl lip, shock-wave/boundary-layer interactions, control of shock-induced flow separation by means of boundary layer bleed, internal normal (terminal) shocks, and the effects of flow incidence. Computed results were compared with available wind tunnel data.
A nongeostrophic theory of zonally averaged circulation is formulated using the nonlinear primitive equations (mass conservation, thermodynamics, and zonal momentum) on a sphere. The relationship between the mean meridional circulation and diabatic heating rate is studied. Differences between results of nongeostropic theory and the geostrophic formulation concerning the role of eddy forcing of the diabatic circulation and the nonlinear nearly inviscid limit versus the geostrophic limit are discussed. Consideration is given to the Eliassen-Palm flux divergence, the Eliassen-Palm pseudodivergence, the nonacceleration theorem, and the nonlinear nongeostrophic Taylor relationship.
High power laser beams can drive plasma particles to relativistic energies. An accurate description of strong waves requires the inclusion of ion dynamics in the analysis. The equations governing the propagation of relativistic electromagnetic waves in a cold electron-ion plasma can be reduced to two equations expressing conservation of energy-momentum of the system. The two conservation constants are functions of the plasma stream velocity, the wave velocity, the wave amplitude, and the electron-ion mass ratio. The dynamic parameter, expressing electron-ion momentum conversation in the laboratory frame, can be regarded as an adjustable quantity, a suitable choice of which will yield self-consistent solutions when other plasma parameters were specified. Circularly polarized electromagnetic waves and electrostatic plasma waves are used as illustrations.
The steady form of the full potential equation, in conservative form, is employed to analyze and design a wide variety of complex aerodynamic shapes. The nonlinear method is based on the theory of characteristic signal propagation coupled with novel flux biasing concepts and body-fitted mapping procedures. The resulting codes are vectorized for the CRAY XMP and the VPS-32 supercomputers. Use of the full potential nonlinear theory is demonstrated for a single-point supersonic wing design and a multipoint design for transonic maneuver/supersonic cruise/maneuver conditions. Achievement of high aerodynamic efficiency through numerical design is verified by wind tunnel tests. Other studies reported include analyses of a canard/wing/nacelle fighter geometry.
The introduction of artificial time derivatives which permit faster convergence to the steady state is used to precondition steady state incompressible equation-solving acceleration methods. An optimal value is obtained for beta, given constant alpha, where alpha and beta are arbitrary functions in the general preconditioning stage. A symmetric hyperbolic system is formed by the resulting incompressible equations. Compressible equations in conservation form with slow flow are also considered.
A new time-accurate Euler method is presented for calculating the transonic flow over a lifting or nonlifting rotor blade in both hover and forward flight. The approach is to solve the conservative Euler equations in a rotor-fixed frame of reference using a finite-volume method. The discretized equations are solved by using a Runge-Kutta multistage scheme with a new higher-order implicit residual smoothing procedure. Preliminary results were compared with wind-tunnel data. In all the cases considered, good agreement was found with experimental data.
Three-dimensional steady and unsteady vortex-dominated flows around sharp-edged delta wings are considered in this paper. The problem is formulated by using the unsteady conservative Euler equations for the flow relative motion with respect to a moving frame of reference. An implicit approximately-factored finite volume scheme is used to solve the resulting equations on a three-dimensional computational grid which is generated by using a modified Joukowski transformation in cross-flow planes at the grid chord stations. The scheme is applied to a delta wing undergoing pitching oscillation around a large angle of attack. The initial conditions correspond to a steady flow around a delta wing of aspect ratio of one, freestream Mach number of 0.3 and mean angle of attack of 20.5. The steady flow results are compared with those of an explicit computational scheme and the experimental data, and they are in good agreement.
Supersonic inlet flows with mixed external-internal compressions were computed using a combined implicit-explicit (Beam-Warming-Steger/MacCormack) method for solving the three-dimensional unsteady, compressible Navier-Stokes equations in conservation form. Numerical calculations were made of various flows related to such inlet operations as the shock-wave intersections, subsonic spillage around the cowl lip, and inlet started versus unstarted conditions. Some of the computed results were compared with wind tunnel data.