Engineering Papers⌕ Search

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.

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

At least 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A critical comparison of two-equation turbulence models

Several two-equation models were proposed and tested against benchmark flows by various researchers. For each study, different numerical methods or codes were used to obtain the results which were reported to be an improvement over other models. However, these comparisons may be overshadowed by the different numerical schemes used to obtain the results. With this in mind, several existing two-equation turbulence models, including k-epsilon, k-tau, k-omega, and q-omega models, are implemented into a common flow solver code for near wall turbulent flows. The quality of each model is based on several criteria, including robustness and accuracy of predicting the turbulent quantities.

Lang, N. J.↗
NASA NTRS1991

Unstructured-grid methods development for unsteady aerodynamic and aeroelastic analyses

The current status of unstructured grid methods development in the Unsteady Aerodynamics Branch at NASA-Langley is described. These methods are being developed for unsteady aerodynamic and aeroelastic analyses. The flow solvers are highlighted which were developed for the solution of the unsteady Euler equations and selected results are given which show various features of the capability. The results demonstrate 2-D and 3-D applications for both steady and unsteady flows. Comparisons are also made with solutions obtained using a structured grid code and with experimental data to determine the accuracy of the unstructured grid methodology. These comparisons show good agreement which thus verifies the accuracy.

Batina, John T.↗
NASA NTRS1991

Parabolized Navier-Stokes methods for hypersonic flows

A representative sampling of the techniques used in the integration of the Parabolized Navier-Stokes (PNS) equations is presented. Special atention is given to recent algorithms developed specifically for application to high speed flows, characterized by the presence of strong embedded shock waves and real gas effects. It is shown that PNS solvers are being used in the analysis of sonic boom signatures. Methods for modeling physical effects are discussed, including an overview of commonly used turbulence models and a more detailed discussion of techniques for including equilibrium and finite rate real gas effects.

Lawrence, Scott L.↗
NASA NTRS1991

An approximate viscous shock layer technique for calculating chemically reacting hypersonic flows about blunt-nosed bodies

An approximate axisymmetric method was developed which can reliably calculate fully viscous hypersonic flows over blunt nosed bodies. By substituting Maslen's second order pressure expression for the normal momentum equation, a simplified form of the viscous shock layer (VSL) equations is obtained. This approach can solve both the subsonic and supersonic regions of the shock layer without a starting solution for the shock shape. The approach is applicable to perfect gas, equilibrium, and nonequilibrium flowfields. Since the method is fully viscous, the problems associated with a boundary layer solution with an inviscid layer solution are avoided. This procedure is significantly faster than the parabolized Navier-Stokes (PNS) or VSL solvers and would be useful in a preliminary design environment. Problems associated with a previously developed approximate VSL technique are addressed before extending the method to nonequilibrium calculations. Perfect gas (laminar and turbulent), equilibrium, and nonequilibrium solutions were generated for airflows over several analytic body shapes. Surface heat transfer, skin friction, and pressure predictions are comparable to VSL results. In addition, computed heating rates are in good agreement with experimental data. The present technique generates its own shock shape as part of its solution, and therefore could be used to provide more accurate initial shock shapes for higher order procedures which require starting solutions.

Cheatwood, F. Mcneil↗
NASA NTRS1990

On the prediction of far field computational aeroacoustics of advanced propellers

A numerical method for determining the acoustic far field generated by a high-speed subsonic aircraft propeller was developed. The approach used in this method was to generate the entire three-dimensional pressure field about the propeller (using an Euler flowfield solver) and then to apply a solution of the wave equation on a cylindrical surface enveloping the propeller. The method is applied to generate the three-dimensional flowfield between two blades of an advanced propeller. The results are compared with experimental data obtained in a wind-tunnel test at a Mach number of 0.6.

Jaeger, Stephen M.↗
NASA NTRS1990

Computation of Space Shuttle high-pressure cryogenic turbopump ball bearing two-phase coolant flow

A homogeneous two-phase fluid flow model, implemented in a three-dimensional Navier-Stokes solver using computational fluid dynamics methodology is described. The application of the model to the analysis of the pump-end bearing coolant flow of the high-pressure oxygen turbopump of the Space Shuttle main engine is studied. Results indicate large boiling zones and hot spots near the ball/race contact points. The extent of the phase change of the liquid oxygen coolant flow due to the frictional and viscous heat fluxes near the contact areas has been investigated for the given inlet conditions of the coolant.

Chen, Yen-Sen↗
NASA NTRS1990

A patched-grid algorithm for complex aircraft configurations

A patched-grid algorithm for the analysis of complex configurations with an implicit, upwind-biased Navier-Stokes solver is presented. Through the use of a generalized coordinate transformation at the zonal interface between two or more blocks, the algorithm can be applied to highly stretched viscous grids and to arbitrarily-shaped patch boundaries. Applications are made to the SR71 reconnaissance aircraft in a high-altitude environment at a supersonic speed and to the F/A-18 forebody-strake configuration at subsonic, high-alpha conditions, in support of the NASA High-Alpha Research Program.

Walters, Robert W.↗
NASA NTRS1990

Prediction of blade-vortex interaction noise using airloads generated by a finite-difference technique

The present numerical finite-difference scheme for helicopter blade-load prediction during realistic, self-generated three-dimensional blade-vortex interactions (BVI) derives the velocity field through a nonlinear superposition of the rotor flow-field yielded by the full potential rotor flow solver RFS2 for BVI, on the one hand, over the rotational vortex flow field computed with the Biot-Savart law. Despite the accurate prediction of the acoustic waveforms, peak amplitudes are found to have been persistently underpredicted. The inclusion of BVI noise source in the acoustic analysis significantly improved the perceived noise level-corrected tone prediction.

Tadghighi, Hormoz↗
NASA NTRS1990

An assessment of full potential and Euler solutions for self-generated rotor blade-vortex interactions

An assessment is made of the accuracy and suitability of the three-dimensional full potential and Euler equations in the modeling of subcritical and supercritical rotor blade-vortex interactions (BVI). The interaction velocity fields are obtained through a nonlinear superposition of the rotor and the vortex wake flow fields. Vortex effects are simulated using the velocity 'transpiration' approach. A modified surface condition is prescribed and enforced at each time step of the computations to satisfy the tangency boundary condition. Potential blade-vortex encounters are identified and tracked in time at equal increments of rotor azimuth using the lifting-line helicopter/rotor trim code CAMRAD. This information is then utilized in interpolation routines within the flow solvers to compute the instantaneous position(s) of the interaction vortex element(s) with respect to the blade for the time-accurate calculations. Results for subcritical self-generated BVI and for supercritical interactions resulting from an upstream-generated vortex are presented.

Hassan, A. A.↗
NASA NTRS1990

On the far-field stream function condition for two-dimensional incompressible flows

The present demonstration of the usefulness of the integral series expansion of the stream function as a far-field computational boundary condition shows the method to require only a 10-percent/time-step increase in computational effort over alternative boundary conditions, in the case of implementation of unsteady problems using a direct elliptic solver. So long as the vorticity was encompassed within the computational domain, the method proved sufficiently accurate to yield virtually identical results for two widely different domains. While the integral-series condition yielded the best results for periodic flow, the Neumann condition gave comparable accuracy with less computation time for the steady-flow case despite its inability to treat periodic flow with vortex shedding.

Sa, Jong-Youb↗
NASA NTRS1991

A multi-block multigrid method for the solution of the Euler and Navier-Stokes equations for three-dimensional flows

A multi-block multigrid method for the solution of the three-dimensional Euler and Navier-Stokes equations is presented. The basic flow solver is a cell-vertex method which employs central-difference spatial approximations and Runge-Kutta time stepping. The use of local time stepping, implicit residual smoothing, multigrid techniques, and variable-coefficient numerical smoothing results in an efficient and robust scheme. The multi-block strategy places the block loop within the Runge-Kutta loop such that accuracy and convergence are not affected by block boundaries. This has been verified by comparing the results of one- and two-block calculations in which the two-block grid is generated by splitting the one-block grid. Results are presented for both Euler and Navier-Stokes computations of wings and wing-fuselage combinations.

Atkins, H. L.↗
NASA NTRS1991

Scramjet mixing establishment times for a pulse facility

A numerical simulation of the temporally developing flow through a generic scramjet combustor duct is presented for stagnation conditions typical of flight at Mach 13 as produced by a shock tunnel pulse facility. The particular focus is to examine the start up transients and to determine the time required for certain flow parameters to become established. The calculations were made with a Navier-Stokes solver SPARK with temporally relaxing inflow conditions derived from operation of the T4 shock tunnel at the University of Queensland in Australia. Calculations at nominal steady inflow conditions were made for comparison. The generic combustor geometry includes the injection of hydrogen fuel from the base of a centrally located strut. In both cases, the flow was assumed laminar and fuel combustion was not included. The establishment process is presented for viscous parameters in the boundary layer and for parameters related to the fuel mixing.

Rogers, R. Clayton↗
NASA NTRS1991

Conical Euler methodology for unsteady vortical flows about rolling delta wings

A conical Euler methodology was developed to study unsteady vortex-dominated flows about rolling highly swept delta wings undergoing either forced or free-to-roll motions. The flow solver of the code involves a multistage Runge-Kutta time-stepping scheme which uses a finite-volume spatial discretization of the Euler equations on an unstructured grid of triangles. Results are presented for a 75-deg swept sharp-leading-edge delta wing at a freestream Mach number of 1.2 and at alpha = 10, 20, and 30 deg. At the 10 and 20 deg, forced harmonic analyses indicate that the rolling moment coefficients provide a positive damping which is verified by free-to-roll calculations. In contrast, at 30 deg, a forced harmonic analysis indicates that the rolling moment coefficient provides a negative damping at the small roll amplitudes. A free-to-roll calculation for this case produces an initially divergent response, but as the amplitude of motion grows with time, the response transitions to a wing-rock type of limit cycle oscillation.

Lee, Elizabeth M.↗
NASA NTRS1991

Cascade flutter analysis with transient response aerodynamics

Two methods for calculating linear frequency domain aerodynamic coefficients from a time-marching Full-Potential cascade solver are developed and verified. In the first method, the Influence Coefficient method, solutions to elemental problems are superposed to obtain the solutions for a cascade in which all blades are vibrating with a constant interblade phase angle. The elemental problem consists of a single blade in the cascade oscillating while the other blades remain stationary. In the second method, the Pulse Response method, the response to the transient motion of a blade is used to calculate influence coefficients. This is done by calculating the Fourier transforms of the blade motion and the response. Both methods are validated by comparison with the Harmonic Oscillation method and give accurate results. The aerodynamic coefficients obtained from these methods are used for frequency domain flutter calculations involving a typical section blade structural model. An eigenvalue problem is solved for each interblade phase angle mode and the eigenvalues are used to determine aeroelastic stability. Flutter calculations are performed for two examples over a range of subsonic Mach numbers using both flat plates and actual airfoils.

Bakhle, Milind A.↗
NASA NTRS1991

An efficient advancing front algorithm for Delaunay triangulation

There has been some recent interest in fluid dynamics calculations on unstructured meshes. One method of unstructured mesh generation involves Delaunay triangulation. This method has certain advantages but it can be expensive to implement. Furthermore, there can be problems with crossing grid lines near boundaries. A method shown here avoids many of the robustness and efficiency problems previously associated with Delaunay triangulation. As an added feature, a simple algorithm is shown which allows removal of diagonal edges from cells that are nearly rectangular. This can result in significant savings in the cost per iteration of a flow solver using this grid.

Merriam, Marshal L.↗
NASA NTRS1990

Flow studies in close-coupled ventral nozzles for STOVL aircraft

Flow in a generic ventral nozzle system was studied experimentally and analytically with the PARC3D computational fluid dynamics program (a full Navier-Stokes equations solver) in order to evaluate the program's ability to predict system performance and internal flow patterns. A generic model of a tailpipe with a rectangular ventral nozzle, about one-third of full size, was tested with unheated air at steady-state pressure ratios up to 4.0. Measurements showed about 5.5 percent flow-turning loss and reasonable nozzle performance coefficients. The flow turned more than the designed 90 deg, causing an aftward axial component in the total thrust. Flow behavior into and through the ventral duct is discussed and illustrated with paint streak flow visualization photographs. PARC3D graphic images are shown for comparison with the experiment photographs. The program successfully predicted internal flow patterns; it also computed thrust and discharge coefficients within 1 percent of measured values.

Mcardle, Jack G.↗
NASA NTRS1991

Development of a solution adaptive unstructured scheme for quasi-3D inviscid flows through advanced turbomachinery cascades

The objective of the present research is to develop a general solution adaptive scheme for the accurate prediction of inviscid quasi-three-dimensional flow in advanced compressor and turbine designs. The adaptive solution scheme combines an explicit finite-volume time-marching scheme for unstructured triangular meshes and an advancing front triangular mesh scheme with a remeshing procedure for adapting the mesh as the solution evolves. The unstructured flow solver has been tested on a series of two-dimensional airfoil configurations including a three-element analytic test case presented here. Mesh adapted quasi-three-dimensional Euler solutions are presented for three spanwise stations of the NASA rotor 67 transonic fan. Computed solutions are compared with available experimental data.

Usab, William J., Jr.↗
NASA NTRS1991

Unstructured Euler flow solutions using hexahedral cell refinement

An attempt is made to extend grid refinement into three dimensions by using unstructured hexahedral grids. The flow solver is developed using the TIGER (topologically Independent Grid, Euler Refinement) as the starting point. The program uses an unstructured hexahedral mesh and a modified version of the Jameson four-stage, finite-volume Runge-Kutta algorithm for integration of the Euler equations. The unstructured mesh allows for local refinement appropriate for each freestream condition, thereby concentrating mesh cells in the regions of greatest interest. This increases the computational efficiency because the refinement is not required to extend throughout the entire flow field.

Melton, John E.↗