Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time-accurate”

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

Massively parallel computing for the simulation of unsteady flows in turbomachinery

This paper deals with evaluating the capabilities of the massively parallel Connection Machine CM2 in predicting unsteady flows in turbomachines. The implementation on the CM2 of an implicit, time-accurate, zonal algorithm for the Navier-Stokes equations in two dimensions is described. Programming issues and modifications made to the original sequential algorithm to improve performance on the CM2 are briefly discussed. Performance is compared to a functionally equivalent code for the Cray YMP.

Madavan, Nateri K.↗

Suppression of pressure oscillations in an open cavity by passive pneumatic control

A computational investigation has been conducted to determine the effectiveness of a passive control technique in suppressing the oscillations in an open cavity exposed to a supersonic flow. Time-accurate solutions of the unsteady, Reynolds-averaged, Navier-Stokes equations were obtained with an explicit predictor-corrector algorithm; the passive control was implemented through the use of a simple linear pressure-velocity law along the porous cavity floor. The computational code was validated by comparisons with experimental data for the cavity flow without control. The computational results with control demonstrate that the oscillations are suppressed; the fluid dynamic mechanism of the control is seen to be a stabilization of the free-shear layer. Spectral analysis of the unsteady data shows that the resonant frequencies are essentially unchanged by the control. Beneficial reductions in the cavity resonant drag are also observed.

Chokani, N.↗

Prediction of unsteady aerodynamic loads in cascades using the linearized Euler equations on deforming grids

A linearized Euler solver for calculating unsteady flows in turbomachinery blade rows due to both incident gusts and blade motion is presented. Using the linearized Euler technique, one decomposes the flow into a mean (or steady) flow plus an unsteady, harmonically varying, small disturbance flow. Linear variable coefficient equations describe the small disturbance behavior of the flow, and are solved using a pseudo-time marching Lax-Wendroff scheme. For the blade motion problem, a harmonically deforming computational rid that conforms to the motion of vibrating blades eliminates large error producing mean flow gradient terms that would otherwise appear in the unsteady flow tangency boundary condition. The paper also presents a new, numerically exact, nonreflecting far-field boundary condition based on an eigenanalysis of the discretized equations. Computed flow solutions demonstrate the computational accuracy and efficiency of the present method. The solution of the linearized Euler equations requires one to two orders of magnitude less computer time than solution of the nonlinear Euler equations using traditional time-accurate time-marching techniques. In addition, the deformable grid significantly improves the accuracy of the solution.

Hall, Kenneth C.↗

Numerical analysis of a nuclear fuel element for nuclear thermal propulsion

A computational fluid dynamics model with porosity and permeability formulations in the transport equations has been developed to study the concept of nuclear thermal propulsion through the analysis of a pulsed irradiation of a particle bed element (PIPE). The numerical model is a time-accurate pressure-based formulation. An adaptive upwind scheme is employed for spatial discretization. The upwind scheme is based on second- and fourth-order central differencing with adaptive artificial dissipation. Multiblocked porosity regions have been formulated to model the cold frit, particle bed, and hot frit. Multiblocked permeability regions have been formulated to describe the flow shaping effect from the thickness-varying cold frit. Computational results for several zero-power density PIPEs and an elevated-particle-temperature PIPE are presented. The implications of the computational results are discussed.

Wang, Ten-See↗

Investigation of vortex development on a pitching slender body of revolution

A computational study of the unsteady flow about a pitching 3.5 caliber tangent ogive forebody is presented. The flow is simulated using the full 3D unsteady Navier-Stokes equations and a time-accurate implicit algorithm. Comparison to available experimental data for a steady 20 deg case is presented as validation. Effects of grid resolution and a comparison of solutions using full Navier-Stokes and the thin-layer approximation are included. The forebody is simulated in a 'pitch-up to 20 deg and hold' maneuver, and two different pitch axis locations are used in the study. Examination of the unsteady vorticity field for the pitch-up cases reveals the formation of strong shear layers as the body decelerates, and their roll-up into vortical structures in a process similar to that observed in airfoil dynamic stall. Pronounced vortex/surface interactions are seen which produce multiple secondary separation regions, ejection of vorticity from the surface, and embedded regions of high suction.

Stanek, M. J.↗

Asymmetric turbulent vortical flows over slender bodies

Time-accurate numerical solutions have been obtained of equations modeling turbulent subsonic flows over a slender ogive-cylinder body of revolution in the high-angle-of-attack regime where a large asymmetry in the mean flow has been observed experimentally. A modified algebraic eddy-viscosity turbulence model was utilized to correctly compute the effects of the asymmetric vortices on the underlying viscous layers. In order to reproduce any one of the experimentally observed asymmetric flow-fields, it was found necessary to add a small geometrical disturbance near the body apex. By determining an appropriate size of the disturbance, it was possible to obtain excellent agreement between numerical results and experimental data for angles of attack of 30 and 40 deg, Reynolds numbers of 3.0 x 10 to the 6th and = 4.0 x 10 to the 6th, and several roll angles. When the disturbance was removed, the flow field returned to its original symmetric shape. These results are similar in behavior to solutions obtained previously for laminar flows. Just as in the laminar case, results suggest that the origin of the asymmetry is a convective-type instability of an originally symmetric flow.

Degani, David↗

Navier-Stokes calculations of transonic flows past cavities

Presented in this paper is a computational investigation of subsonic and transonic flows past three-dimensional deep and transitional cavities. Simulations of these self-induced oscillatory flows have been generated through time-accurate solutions of the Reynolds averaged, full Navier-Stokes equations, using the explicit MacCormack scheme. The Reynolds stresses have been included through the Baldwin-Lomax algebraic turbulence model with certain modifications. The computational results include instantaneous and time averaged flow properties. The results of an experimental investigation have been used not only to validate the time-averaged results, but also to investigate the effects of varying the Mach number and the incoming boundary-layer thickness. Time series analyses have been performed for the instantaneous pressure values on the cavity floor and compared with the results obtained by a predictive formula. While most of the comparisons have been favorable, some discrepancies have been observed, particularly on the rear face. The present results help understanding the three-dimensional and unsteady features of the separations, vortices, the shear layer, as well as some of the aeroacoustic phenomena of compressible cavity flows.

Baysal, O.↗

Leading-edge receptivity to a vortical freestream disturbance: A numerical analysis

The receptivity to freestream vorticity of the boundary layer over a flat plate with an elliptic leading edge is investigated numerically. The flow is simulated by solving the incompressible Navier-Stokes system in general curvilinear coordinates with the vorticity and stream function as dependent variables. A finite-difference scheme which is second-order accurate in both space and time is used. As a first step, the steady basic-state solution is computed. Then a small amplitude vortical disturbance is introduced at the upstream boundary and the governing equations are solved time-accurately to evaluate the spatial and temporal growth of the perturbations leading to instability waves (Tollmien-Schlichting waves) inside the boundary layer. Preliminary results for a symmetric, 2-D disturbance reveal the presence of Tollmien-Schlichting waves aft of the flat-plate/ellipse juncture.

Buter, Thomas A.↗

Application of program LAURA to perfect gas shock tube flows: A parametric study

The Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) was originally developed to solve steady-flow problems. The desire to validate the algorithm with shock tube experimental data motivated the development of a time-accurate version of the LAURA code. The current work presents a test of the Algorithm. Computational results are compared with the exact solution for a simple shock tube case. The parameters examined are Courant number, relaxation sweeps, grid spacing, and the inviscid relaxation factor. The results of the study indicate that LAURA is capable of producing accurate solutions when appropriate values are used for each parameter.

Mitterer, K. F.↗

Computation of incompressible viscous flows through turbopump components

A finite-difference, three-dimensional, incompressible Navier-Stokes formulation for calculating the flow through turbopump components is presented. The solution method is based on the pseudocompressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line-relaxation method. Both steady and unsteady flow calculations can be performed using the presented algorithm. In this paper, the equations are solved in steadily rotating reference frames by using the steady-state formulation in order to simulate the flow through a turbopump inducer. Eddy viscosity is computed by using the Baldwin-Lomax model. Numerical results are compared with experimental measurements and good agreement is found between the two. Time-accurate calculations will be reported in future publications.

Kiris, Cetin↗

Transient cooling of a square region of radiating medium

A time-accurate numerical solution was carried out for transient radiative cooling of a gray emitting and absorbing medium in a square two-dimensional region. The integro-differential energy equation for transient temperature distributions was solved in two stages. At each time increment, the local radiative source term was obtained by numerical integration of the temperature field using two-dimensional Gaussian integration over rectangular subregions. Then the differential portion of the equation was integrated forward in time by use of the local first and second time derivatives. The results were compared with available limiting case, and excellent agreement was obtained. Transient results are given for a wide range of optical thicknesses of the region. Optimum transient cooling is obtained when the optical side length is about 4.

Siegel, Robert↗

Prediction of unsteady rotor-surface pressure and heat transfer from wake passings

The research described in this paper is a numerical investigation of the effects of unsteady flow on gas turbine heat transfer, particularly on a rotor blade surface. The unsteady flow in a rotor blade passage and the unsteady heat transfer on the blade surface as a result of wake/blade interaction are modeled by the inviscid flow/boundary layer approach. The Euler equations which govern the inviscid flow are solved using a time-accurate marching scheme. The unsteady flow in the blade passage is induced by periodically moving a wake model across the passage inlet. Unsteady flow solutions in the passage provide pressure gradients and boundary conditions for the boundary-layer equations which govern the viscous flow adjacent to the blade surface. Numerical solutions of the unsteady turbulent boundary layer yield surface heat flow values which can then be compared to experimental data. Comparisons with experimental data show that unsteady heat flux on the blade suction surface is well predicted, but the predictions of unsteady heat flux on the blade pressure surface do not agree.

Tran, Le T.↗

Time-dependent solutions of viscous incompressible flows in moving co-ordinates

A time-accurate solution method for the incompressible Navier-Stokes equations in generalized moving coordinates is presented. A finite-volume discretization method that satisfies the geometric conservation laws for time-varying computational cells is used. The discrete equations are solved by a fractional-step solution procedure. The solution is second-order-accurate in space and first-order-accurate in time. The pressure and the volume fluxes are chosen as the unknowns to facilitate the formulation of a consistent Poisson equation and thus to obtain a robust Poisson solver with favorable convergence properties. The method is validated by comparing the solutions with other numerical and experimental results. Good agreement is obtained in all cases.

Rosenfeld, Moshe↗

Time accurate Navier-Stokes computations of solid rocket motor internal field joint cavity flows

The low-speed flow of combustion gases in the Space Shuttle Solid Rocket Motor (SRM) following ignition are investigated with a time-accurate CFD technique. The flow of the combustion gases is examined in the forward joint section of the SRM at 80 s after ignition with attention given to the amplification of existing 'organ-pipe' modes of 15, 30, and 45 Hz resulting from frequency lock-in with vortex-shedding modes. The CFD code employed is a multizone compressible Navier-Stokes flow solver using a finite-volume, upwind-biased, total-variation-diminishing scheme. The vortex-generated oscillations generated by the inhibitor and the cavity at the forward field joint are described in terms of amplitude and frequency and related to the flow dynamics in the SRM. Oscillations are shown to occur spontaneously when the trailing edge of the cavity at the joint is sharp and can be amplified if the unburned portion of the inhibitor protrudes into the flow.

Farr, Rebecca A.↗

Reduction of vortex-induced vibration in vane geometries

Computations using a time-accurate, compressible Navier-Stokes flow model were conducted to analyze both unsteady laminar and turbulent flows over two curved configurations of vanes which were shaped to treat a 4.0 kHz vibration problem which has occurred in several SSMEs. This problem involves vortex shedding from vanes which causes excessive vane vibration and cracking. The original vane configuration exhibited strong flow-induced vibrations at a Strouhal number near 0.19 for a the first bending mode excitation. Scalloping of the leading edge raised the frequency of the vane's first torsional mode, which in turn increased the onset flow velocity at lock-in. Beveling the vane's trailing edge eliminated vortex shedding at the trailing edge, which decoupled the flow oscillation from the vibrating wave. A modified vane configuration with a beveled trailing edge at a 30 deg angle was also studied.

Liu, B. L.↗

Three dimensional flows in multistage compressors

Current multi-stage turbomachinery design/analysis methods are based on a time-averaged, axisymmetric representation of the flow field. The actual flow field is asymmetric and unsteady due to blade row interactions. The Reynolds averaged Navier-Stokes solvers are limited to single-stage machines for existing computers. Therefore, advanced multi-stage compressors will operate far off-design for portions of the flight regime. The objectives are to provide an experimentally validated average-passage calculation of multistage compressor blade row interactions and an experimentally validated time-accurate calculation of multi-stage compressor blade row interactions. Various aspects of this investigation are presented in viewgraph form.

Hathaway, Michael D.↗

Development of generalized pressure velocity coupling scheme for the analysis of compressible and incompressible combusting flows

The objective of this investigation has been to develop an algorithm (or algorithms) for the improvement of the accuracy and efficiency of the computer fluid dynamics (CFD) models to study the fundamental physics of combustion chamber flows, which are necessary ultimately for the design of propulsion systems such as SSME and STME. During this three year study (May 19, 1978 - May 18, 1992), a unique algorithm was developed for all speed flows. This newly developed algorithm basically consists of two pressure-based algorithms (i.e. PISOC and MFICE). This PISOC is a non-iterative scheme and the FICE is an iterative scheme where PISOC has the characteristic advantages on low and high speed flows and the modified FICE has shown its efficiency and accuracy to compute the flows in the transonic region. A new algorithm is born from a combination of these two algorithms. This newly developed algorithm has general application in both time-accurate and steady state flows, and also was tested extensively for various flow conditions, such as turbulent flows, chemically reacting flows, and multiphase flows.

Chen, C. P.↗

Investigation of advanced counterrotation blade configuration concepts for high speed turboprop systems. Task 2: Unsteady ducted propfan analysis computer program users manual

The primary objective of this study was the development of a time-dependent three-dimensional Euler/Navier-Stokes aerodynamic analysis to predict unsteady compressible transonic flows about ducted and unducted propfan propulsion systems at angle of attack. The computer codes resulting from this study are referred to as Advanced Ducted Propfan Analysis Codes (ADPAC). This report is intended to serve as a computer program user's manual for the ADPAC developed under Task 2 of NASA Contract NAS3-25270, Unsteady Ducted Propfan Analysis. Aerodynamic calculations were based on a four-stage Runge-Kutta time-marching finite volume solution technique with added numerical dissipation. A time-accurate implicit residual smoothing operator was utilized for unsteady flow predictions. For unducted propfans, a single H-type grid was used to discretize each blade passage of the complete propeller. For ducted propfans, a coupled system of five grid blocks utilizing an embedded C-grid about the cowl leading edge was used to discretize each blade passage. Grid systems were generated by a combined algebraic/elliptic algorithm developed specifically for ducted propfans. Numerical calculations were compared with experimental data for both ducted and unducted propfan flows. The solution scheme demonstrated efficiency and accuracy comparable with other schemes of this class.

Hall, Edward J.↗