SEARCH · Engineering Papers

Results for “solvers”

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 1,243 records · Page 69

All sources NASA NTRS DOE OSTI

NASA NTRS1986

The calculation of supersonic viscous flows using the parabolized Navier-Stokes equations

Solution of the parabolic Navier-Stokes (PNS) equations for supersonic flows is discussed, and compatibility of the PNS method with the triple-deck theory of Stewartson (1974) is demonstrated. Characteristic and stability analyses show that use of an appropriate filter on the pressure term in the x-momentum equation can suppress the departure solutions, giving the usual desired weak interaction solution. An Alternating Direction Explicit procedure, with minimal computer storage requirements compared to the full Navier-Stokes solvers, is proposed to calculate strongly interacting flows using a global iteration procedure for the PNS equations. The PNS equations are used to solve the hypersonic viscous interaction problem, and good agreement is found with experimental results.

Davis, R. T.↗

NASA NTRS1987

Multigrid solution of inviscid transonic flow through rotating blade passages

A fast Euler solver for three dimensional inviscid transonic flow in rotating domains is described. The time dependent Euler equations are discretized spatially with finite volumes, and are advanced temporally with a multiple stage time stepping scheme. A dramatic increase in the rate of convergence for steady solutions is achieved with a multigrid algorithm that employs the multistage scheme as its smoothing procedure. The effectiveness of the multistage scheme as a multigrid driver is enhanced by the utilization of analytically determined combinations of the governing parameters.

Smith, Wayne A.↗

NASA NTRS1987

A generalized procedure for constructing an upwind-based TVD scheme

A generalized formulation for constructing second- and higher-order accurate TVD (total variation diminishing) schemes is presented. A given scheme is made TVD by limiting antidiffusive flux differences with some nonlinear functions, so-called limiters. The general idea of the formulation and its mathematical proof of Harten's TVD conditions is shown by applying the Lax-Wendroff method to a scalar nonlinear equation and constant-coefficient system of conservation laws. For the system of equations, several definitions are derived for the argument used in the limiter function and present their performance to numerical experiments. Then the formulation is formally extended to the nonlinear system of equations. It is demonstrated that use of the present procedure allows easy conversion of existing central or upwind, and second- or higher-order differencing schemes so as to preserve monotonicity and to yield physically admissible solutions. The formulation is simple mathematically as well as numerically; neither matrix-vector multiplication nor Riemann solver is required. Roughly twice as much computational effort is needed as compared to conventional scheme. Although the notion of TVD is based on the initial value problem, application to the steady Euler equations of the formulation is also made. Numerical examples including various ranges of problems show both time- and spatial-accuracy in comparison with exact solutions.

Liou, Meng-Sing↗

NASA NTRS1987

Analysis of viscous transonic flow over airfoil sections

A full Navier-Stokes solver has been used to model transonic flow over three airfoil sections. The method uses a two-dimensional, implicit, conservative finite difference scheme for solving the compressible Navier-Stokes equations. Results are presented as prescribed for the Viscous Transonic Airfoil Workshop to be held at the AIAA 25th Aerospace Sciences Meeting. The NACA 0012, RAE 2822 and Jones airfoils have been investigated for both attached and separated transonic flows. Predictions for pressure distributions, loads, skin friction coefficients, boundary layer displacement thickness and velocity profiles are included and compared with experimental data when possible. Overall, the results are in good agreement with experimental data.

Huff, Dennis L.↗

NASA NTRS1986

Fast algorithm for calculating chemical kinetics in turbulent reacting flow

This paper addresses the need for a fast batch chemistry solver to perform the kinetics part of a split operator formulation of turbulent reacting flows, with special attention focused on the solution of the ordinary differential equations governing a homogeneous gas-phase chemical reaction. For this purpose, a two-part predictor-corrector algorithm which incorporates an exponentially fitted trapezoidal method was developed. The algorithm performs filtering of ill-posed initial conditions, automatic step-size selection, and automatic selection of Jacobi-Newton or Newton-Raphson iteration for convergence to achieve maximum computational efficiency while observing a prescribed error tolerance. The new algorithm, termed CREK1D (combustion reaction kinetics, one-dimensional), compared favorably with the code LSODE when tested on two representative problems drawn from combustion kinetics, and is faster than LSODE.

Radhakrishnan, K.↗

NASA NTRS1986

Multigrid method for nearly singular and slightly indefinite problems

This paper deals with nearly singular, possibly indefinite problems for which the usual multigrid solvers converge very slowly or even diverge. The main difficulty is related to some badly approximated smooth functions which correspond to eigenfunctions with nearly zero eigenvalues. A correction to the usual coarse-grid equations is derived, both in the correction scheme and in the full approximation scheme. The performance of the new algorithm using this correction is essentially as that of usual multigrid for definite problems.

Brandt, A.↗

NASA NTRS1987

Studies of vortex dominated flows; Proceedings of the Symposium, Hampton, VA, July 9-11, 1985

Papers are presented on waves and bifurcations in vortex filaments, a ring-vortex representation of an axisymmetric vortex sheet, and comparison of experiment with the dynamics of the von Karman vortex trail. Also considered are force-free and loss-free transitions between vortex flow states, a vortex breakdown simulation based on a nonlinear inviscid method, and the prediction of highly vortical flows using an Euler equation model. Other topics include the theory of high-Reynolds-number flow past a blunt body, progress on the calculation of large-scale separation at high Reynolds numbers, and viscous-inviscid interaction solvers and computation of highly separated flows. Papers are also presented on simulation studies of vortex dynamics of a leading edge vortex flap, methods for numerical simulation of leading edge vortex flow, and comparison of measured and computed pitot pressures in a leading edge vortex from a delta wing.

Hussaini, M. Y.↗

NASA NTRS1987

Comparison of measured and computed pitot pressures in a leading edge vortex from a delta wing

Calculations are presented for a 75-deg swept flat plate wing tested at a freestream Mach number of 1.95 and 10 degrees angle of attack. Good agreement is found between computational data and previous experimental pitot pressure measurements in the core of the vortex, suggesting that the total pressure losses predicted by the Euler equation solvers are not errors, but realistic predictions. Data suggest that the magnitude of the total pressure loss is related to the circumferential velocity field through the vortex, and that it increases with angle of attack and varies with Mach number and sweep angle.

Murman, Earll M.↗

NASA NTRS1987

Three dimensional mesh generation by triangulation of arbitrary point sets

A method for generating an unstructured mesh is described. The approach is quite general and joins an arbitrary set of points to produce a covering of three dimensional space by tetrahedra. After removing the tetrahedra that connect surface points, a mesh suitable for a finite element based flow solver is obtained. Details of the triangulation algorithm are provided together with an analysis of the algorithm efficiency and validity.

Baker, Timothy J.↗

NASA NTRS1987

Grid generation and inviscid flow computation about cranked-winged airplane geometries

An algebraic grid generation procedure that defines a patched multiple-block grid system suitable for fighter-type aircraft geometries with fuselage and engine inlet, canard or horizontal tail, cranked delta wing and vertical fin has been developed. The grid generation is based on transfinite interpolation and requires little computational power. A finite-volume Euler solver using explicit Runge-Kutta time-stepping has been adapted to this grid system and implemented on the VPS-32 vector processor with a high degree of vectorization. Grids are presented for an experimental aircraft with fuselage, canard, 70-20-cranked wing, and vertical fin. Computed inviscid compressible flow solutions are presented for Mach 2 at 3.79, 7 and 10 deg angles of attack. Conmparisons of the 3.79 deg computed solutions are made with available full-potential flow and Euler flow solutions on the same configuration but with another grid system. The occurrence of an unsteady solution in the 10 deg angle of attack case is discussed.

Eriksson, L.-E.↗

NASA NTRS1987

An adaptive finite element scheme for the Euler and Navier-Stokes equations

The implementation of two explicit finite element schemes for the solution of the compressible Euler and Navier-Stokes equations is presented. The schemes can be employed with general unstructured triangular meshes in two dimensions. Either scheme can therefore be used as the basic solver in a solution adaptive mesh procedure in a direct manner. The particular adaptive approach which is advocated here is intended for the solution of steady state problems only and involves an adaptive regeneration of the grid at prescribed stages during the false transient. The grid regeneration is accomplished by a mesh generator which has the capability of generating triangular grids over computational domains of arbitrary shape. The procedure is illustrated by solving transonic flows over multi-airfoil configurations and high speed flows, involving shock interactions, past circular cylinders.

Morgan, K.↗

NASA NTRS1987

Aerodynamic design modification of a hypersonic wind tunnel nozzle by CSCM with high order accuracy

An improved version of the CSCM implicit Navier-Stokes solver with a flexible data structure has been applied, in conjunction with a sophisticated patched grid system, to the aerodynamic analysis of hypersonic axisymmetric contoured nozzles of the NASA-Ames 3.5-foot hypersonic wind tunnel. The results obtained are found to be in good agreement with experiment. Based on the analysis, a newly designed throat contour for the Mach 14 nozzle has been recommended.

Yang, J. Y.↗

NASA NTRS1987

A theoretical evaluation of laser-sustained plasma thruster performance

An extensive numerical experiment has been conducted to evaluate rocket thruster performance using a laser-sustained hydrogen plasma as the propellant. The plasma was sustained using a 30 kW CO2 laser beam operated at 10.6 microns focused inside the thruster. The steady-state Navier-Stokes equations coupled with the laser power absorption process have been solved numerically. A pressure based Navier-Stokes solver using body-fitted coordinate was used to calculate the laser-supported rocket flow which included both recirculating and transonic flow regions. The local thermodynamic equilibrium (LTE) assumption was used for the plasma thermophysical and optical properties. Geometric ray tracing was adopted to describe the laser beam. Several different throat size thrusters operated at 150 and 300 kPa chamber stagnation pressure were studied. It was found that the thruster performance (vacuum specific impulse) was highly dependent on the operating conditions, and a properly designed laser supported thruster can attain a specific impulse around 1500 secs. The heat loading on the thruster wall was also estimated and was in the range of that for a conventional chemical rocket.

Jeng, San-Mou↗

NASA NTRS1987

Euler and Navier-Stokes solutions for the leeside flow over delta wings at supersonic speeds

Distinctly different types of leeside flowfields over highly swept sharp leading edge delta wings in supersonic flow were numerically simulated using Euler and Navier-Stokes solvers. The Euler code was seen to be adequate only in predicting primary flow structures (leading edge vortex and cross flow shock) whereas the Navier-Stokes code was capable of predicting secondary flow structures (i.e., secondary vortex). A comparison of laminar and turbulent Navier-Stokes solutions indicated that the turbulent boundary layer model is more accurate in predicting the effect of the boundary layer model on the flowfield. Also, the Navier-Stokes code indicated detailed flow structures not observed in the qualitative experimental data available (i.e., vapor screen photographs) indicating a need for quantitative flow field data.

Mcmillin, S. N.↗

NASA NTRS1987

A strake design method for supersonic speeds and low lift

A method is presented for the design of a strake for a forebody typical of supersonic aircraft at typical cruise (low lift) conditions, in order to ensure attached flow on the strake upper surface. The method systematically varies strake camber and thickness until an attached flow surface pressure distribution is achieved. Surface pressures are computed by the full potential flow solver, NCOREL. The high-alpha flowfield of the strake is also investigated experimentally and computationally; a large vortex is noted on the strake leeside, together with a pair of vortices on the forebody upper centerline. The effect of sideslip on the high-alpha flowfield is also discussed.

Pittman, James L.↗

NASA NTRS1987

Inviscid transonic wing design using inverse methods in curvilinear coordinates

An inverse wing design method has been developed around an existing transonic wing analysis code. The original analysis code, TAWFIVE, has as its core the numerical potential flow solver, FLO30, developed by Jameson and Caughey. Features of the analysis code include a finite-volume formulation; wing and fuselage fitted, curvilinear grid mesh; and a viscous boundary layer correction that also accounts for viscous wake thickness and curvature. The development of the inverse methods as an extension of previous methods existing for design in Cartesian coordinates is presented. Results are shown for inviscid wing design cases in super-critical flow regimes. The test cases selected also demonstrate the versatility of the design method in designing an entire wing or discontinuous sections of a wing.

Gally, Thomas A.↗

NASA NTRS1987

Computation of turbulent flows using an extended k-epsilon turbulence closure model

An extended kappa-epsilon turbulence model is proposed and tested with successful results. An improved transport equation for the rate of dissipation of the turbulent kinetic energy, epsilon, is proposed. The proposed model gives more effective response to the energy production rate than does the standard kappa-epsilon turbulence model. An extra time scale of the production range is included in the dissipation rate equation. This enables the present model to perform equally well for several turbulent flows with different characteristics, e.g., plane and axisymmetric jets, turbulent boundary layer flow, turbulent flow over a backward-facing step, and a confined turbulent swirling flow. A second-order accurate finite difference boundary layer code and a nearly second-order accurate finite difference elliptic flow solver are used for the present numerical computations.

Chen, Y.-S.↗

NASA NTRS1987

A full potential flow analysis with realistic wake influence for helicopter rotor airload prediction

A 3-D, quasi-steady, full potential flow solver was adapted to include realistic wake influence for the aerodynamic analysis of helicopter rotors. The method is based on a finite difference solution of the full potential equation, using an inner and outer domain procedure for the blade flowfield to accommodate wake effects. The nonlinear flow is computed in the inner domain region using a finite difference solution method. The wake is modeled by a vortex lattice using prescribed geometry techniques to allow for the inclusion of realistic rotor wakes. The key feature of the analysis is that vortices contained within the finite difference mesh (inner domain) were treated with a vortex embedding technique while the influence of the remaining portion of the wake (in the outer domain) is impressed as a boundary condition on the outer surface of the finite difference mesh. The solution procedure couples the wake influence with the inner domain solution in a consistent and efficient solution process. The method has been applied to both hover and forward flight conditions. Correlation with subsonic and transonic hover airload data is shown which demonstrates the merits of the approach.

Egolf, T. Alan↗