Search NASA⌕ Search

SEARCH · Search NASA

Results for “implicit schemes”

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 181 records · Page 10

Compressible studies on dynamic stall

The purpose of this work is to investigate the effects of compressibility on the dynamic stall phenomenon by numerical simulation of the unsteady flow. The full two-dimensional unsteady compressible Navier-Stokes equations are solved for flows over oscillating airfoils and airfoils pitching rapidly to high angles of attack. The free-stream speeds vary from low subsonic with mild compressibility effects, to moderate subsonic where strong compressibility effects appear close to the leading edge at high angles of attack. An Alternating Direction Implicit scheme is implemented for the numerical solution with the viscous terms retained in both directions. The numerical results are compared with available experimental data for a Sikorsky airfoil for compressible high Reynolds number flows. There is good agreement between the computed and measured unsteady lift and pitching moment coefficient time histories. The computed high-speed subsonic unsteady results give a good picture of the entire flow field, and the dynamic stall progression in the compressible flow regime. It was observed that compressibility effects are more severe close to the leading edge at moderate angles of attack, and that the dynamic stall vortex appears at lower angles of attack as the free-stream speed increases.

Ekaterinaris, J. A.↗

Application of a full-potential solver to bending-torsion flutter in cascades

A two-dimensional, unsteady, full-potential cascade aerodynamic code has been used for flutter calculations. The code, which uses a time-marching implicit scheme with internal Newton-iterations at each time step, allows oscillatory blade motions with phase lag between adjacent blades to be modeled using a single inter-blade passage. This aerodynamic code has been used with a two-degrees-of-freedom typical section structural model to determine the aeroelastic stability of a bladed disk. The unsteady pressure results calculated from the aerodynamic code have been verified by comparison with results from the linear and linearized potential theories. Flutter results have been presented for two cases. A comparison with linear theory shows identical trends; however, the present calculations show a lower flutter frequency.

Bakhle, Milind A.↗

LU-SGS implicit algorithm for three-dimensional incompressible Navier-Stokes equations with source term

A numerical method is developed for solving the incompressible Navier-Stokes equations using the concept of pseudocompressibility. A lower-upper symmetric-Gauss-Seidel implicit scheme is developed for three-dimensional incompressible viscous flow computations. The present algorithm offers additional advantages when solving the flow equations with source terms. Complete vectorizability of the algorithm on oblique planes of sweep in three-dimensions is accomplished in a new flow solver, INS3D-LU code. Spatial differencing is a second-order accurate semi-discrete finite-volume method augmented by a third-order accurate numerical dissipation model which is based on spectral-radii. Comparison of numerical solutions for a curved duct with experimental data shows good agreement. The method is applied to calculate the inducer flow of the Space Shuttle Main Engine turbopump.

Yoon, Seokkwan↗

Numerical simulations in combustion

This paper reviews numerical simulations in reacting flows in general and combustion phenomena in particular. It is shown that use of implicit schemes and/or adaptive mesh strategies can improve convergence, stability, and accuracy of the solution. Difficulties increase as turbulence and multidimensions are considered, particularly when finite-rate chemistry governs the given combustion problem. Particular attention is given to the areas of solid-propellant combustion dynamics, turbulent diffusion flames, and spray droplet vaporization.

Chung, T. J.↗

Implicit and explicit computations of flows past cavities with and without yaw

Implicit and explicit computations are presented to study the supersonic, turbulent flows past three-dimensional, rectangular cavities at 0 and 45 deg yaw angle. The Reynolds-averaged, unsteady, compressible, complete Navier-Stokes equations are solved time-accurately, including the Reynolds stresses. The number of major vortices inside and the separation characteristics on the cavity floor are different for cavities with length-to-depth ratios of 3.0 and 6.7. Due to the oscillations of the shear layer, which randomly exposes the cavity flow to the external flow, vortices are shed into the mainstream. When the upstream flow approaches a cavity at yaw, the width-to-depth ratio also becomes an important parameter. The dominant rotational planes of the major vortices incline at an angle to the vertical plane due to yaw. The execution of the implicit scheme is about three times faster, but it requires about four times more computer memory.

Baysal, Oktay↗

Preconditioned conjugate gradient methods for the compressible Navier-Stokes equations

The compressible Navier-Stokes equations are solved for a variety of two-dimensional inviscid and viscous problems by preconditioned conjugate gradient-like algorithms. Roe's flux difference splitting technique is used to discretize the inviscid fluxes. The viscous terms are discretized by using central differences. An algebraic turbulence model is also incorporated. The system of linear equations which arises out of the linearization of a fully implicit scheme is solved iteratively by the well known methods of GMRES (Generalized Minimum Residual technique) and Chebyschev iteration. Incomplete LU factorization and block diagonal factorization are used as preconditioners. The resulting algorithm is competitive with the best current schemes, but has wide applications in parallel computing and unstructured mesh computations.

Venkatakrishnan, V.↗

Interactive grid generation for turbomachinery flow field simulations

The control point form of algebraic grid generation presented provides the means that are needed to generate well structured grids of turbomachinery flow simulations. It uses a sparse collection of control points distributed over the flow domain. The shape and position of coordinate curves can be adjusted from these control points while the grid conforms precisely to all boundaries. An interactive program called TURBO, which uses the control point form, is being developed. Basic features of the code are discussed and sample grids are presented. A finite volume LU implicit scheme is used to simulate flow in a turbine cascade on the grid generated by the program.

Choo, Yung K.↗

An upwind approach to unsteady flowfield simulation

A numerical method to determine unsteady solutions of the laminar, perfect gas Navier-Stokes equations has been developed. The structured finite-volume, approximately factored implicit scheme uses Newton subiterations to obtain the spatially and temporally second-order accurate time history of the interaction of blast-waves with stationary targets. The inviscid flux is evaluated using either of two upwind techniques, while the viscous terms are computed by central differencing. Comparisons of numerical, analytical, and experimental results are made in two and three dimensions. The results show accurate wave speed resolution and nonoscillatory discontinuity capturing.

Atwood, Christopher A.↗

Three-dimensional Euler time accurate simulations of fan rotor-stator interactions

A numerical method useful to describe unsteady 3-D flow fields within turbomachinery stages is presented. The method solves the compressible, time dependent, Euler conservation equations with a finite volume, flux splitting, total variation diminishing, approximately factored, implicit scheme. Multiblock composite gridding is used to partition the flow field into a specified arrangement of blocks with static and dynamic interfaces. The code is optimized to take full advantage of the processing power and speed of the Cray Y/MP supercomputer. The method is applied to the computation of the flow field within a single stage, axial flow fan, thus reproducing the unsteady 3-D rotor-stator interaction.

Boretti, A. A.↗

Navier-Stokes analysis of turbomachinery blade external heat transfer

The two-dimensional, compressible, thin-layer Navier-Stokes and energy equations were solved numerically to obtain heat transfer rates on turbomachinery blades. The Baldwin-Lomax algebraic model and the q-omega low Reynolds number two-equation model were used for modeling of turbulence. For the numerical solution of the governing equations a four-stage Runge-Kutta solver was employed. The turbulence model equations were solved using an implicit scheme. Numerical solutions are presented for two-dimensional flow within two vane cascades. The heat transfer results and the pressure distributions were compared with published experimental data. The agreement between the numerical calculations and the experimental values were found to be generally favorable. The position of transition from laminar to turbulent flow was also predicted accurately.

Gorla, Rama S. R.↗

Combustion of liquid fuel droplets in supercritical conditions

A comprehensive analysis of liquid-fuel droplet combustion in both sub- and super-critical environments has been conducted. The formulation is based on the complete conservation equations for both gas and liquid phases, and accommodates finite-rate chemical kinetics and a full treatment of liquid-vapor phase equilibrium at the droplet surface. The governing equations and the associated interface boundary conditions are solved numerically using a fully coupled, implicit scheme with the dual time-stepping integration technique. The model is capable of treating the entire droplet history, including the transition from the subcritical to the supercritical state. As a specific example, the combustion of n-pentane fuel droplets in air is studied for pressures of 5-140 atm. Results indicate that the ambient gas pressure exerts significant control of droplet gasification and burning processes through its influences on the fluid transport, gas/liquid interface thermodynamics, and chemical reactions. The droplet gasification rate increases progressively with pressure. However, the data for the overall burnout time exhibits a significant variation near the critical burning pressure, mainly as a result of reduced mass-diffusion rate and latent heat of vaporization with increased pressure. The influence of droplet size on the burning characteristics is also noted.

Shuen, J. S.↗

An implicit three-dimensional Navier-Stokes solver for compressible flow

A three-dimensional numerical method based on the lower-upper symmetric-Gauss-Seidel implicit scheme in conjunction with the flux-limited dissipation model is developed for solving the compressible Navier-Stokes equations. A new computer code which is based on this method requires only 9 microsec per grid-point per iteration on a single processor of a Cray YMP computer and executes at the sustained rate of 170 MFLOPS. A reduction of 4 orders of magnitude in the residual for a high Reynolds number flow using 230 K grid points is obtained in 24 minutes. The computational results compare well with available experimental data.

Yoon, Seokkwan↗

Diagonal implicit multigrid solution of compressible turbulent flows

The Caughey (1988) multigrid diagonal implicit algorithm for Euler equations is extended to solve the two-dimensional thin layer Navier-Stokes equations for turbulent transonic flows. To further accelerate convergence to a steady state, the implicit scheme is used within the multigrid method. Results are presented for transonic flows past airfoils, and the flow-field results are compared with other computational data and experiments demonstrating the accuracy of the method.

Varma, R. R.↗

Multitasking the INS3D-LU code on the Cray Y-MP

This paper presents the results of multitasking the INS3D-LU code on eight processors. The code is a full Navier-Stokes solver for incompressible fluid in three dimensional generalized coordinates using a lower-upper symmetric-Gauss-Seidel implicit scheme. This code has been fully vectorized on oblique planes of sweep and parallelized using autotasking with some directives and minor modifications. The timing results for five grid sizes are presented and analyzed. The code has achieved a processing rate of over one Gflops.

Fatoohi, Rod↗

The effect of bleed configuration on shock/boundary layer interactions

The flow field characteristics are simulated numerically in an oblique shock wave/laminar boundary layer interaction for three different bleed configurations. The numerical solution for the flow field is obtained for the strong conservation-law form of the two-dimensional compressible Navier-Stokes equations using an implicit scheme. The computations model the flow in the interaction region and inside the bleed slot for an impinging oblique shock of sufficient strength to cause boundary layer separation on a flat plate boundary layer. The computed results are presented for a normal bleed slot at three different locations, upstream, across and downstream of the shock impingement point. The detailed flow characteristics in the interaction zone and inside the bleed slot are compared for the three cases. The resulting surface pressure and shear stress distributions are also compared for the three bleed cases and for the shock wave boundary layer interaction without bleed.

Hamed, A.↗

Three-dimensional incompressible Navier-Stokes solver using lower-upper symmetric-Gauss-Seidel algorithm

A numerical method based on the pseudocompressibility concept is developed for solving the three-dimensional incompressible Navier-Stokes equations using the lower-upper symmetric-Gauss-Seidel implicit scheme. Very high efficiency is achieved in a new flow solver, INS3D-LU code, by accomplishing the complete vectorizability of the algorithm on oblique planes of sweep in three dimensions.

Yoon, Seokkwan↗

Vaporization and combustion of fuel droplets at supercritical conditions

Vaporization and combustion liquid-fuel droplets in both sub- and super-critical environments have been examined. The formulation is based on the complete conservation equations for both gas and liquid phases, and accommodates finite-rate chemical kinetics and a full treatment of liquid-vapor phase equilibrium at the droplet surface. The governing equations and the associated interface boundary conditions are solved numerically using a fully coupled, implicit scheme with the dual time-stepping integration technique. The model is capable of treating the entire droplet history, including the transition from the subcritical to the supercritical state. As a specific example, the combustion of n-pentane fuel droplets in air is studied for pressures of 5-140 atm. In addition, the dynamic responses of droplet vaporization and combustion to ambient-pressure oscillations are investigated. Results indicate that the droplet gasification and burning mechanisms depend greatly on the ambient pressure. In particular, a rapid enlargement of the vaporization and combustion responses occurs when the droplet surface reaches its critical point, mainly due to the strong variations of latent heat of vaporization and thermophysical properties at the critical state.

Yang, Vigor↗

Navier-Stokes analysis of turbomachinery blade external heat transfer

The two-dimensional, compressible, thin-layer Navier-Stokes and energy equations were solved numerically to obtain heat transfer rates on turbomachinery blades. The Baldwin-Lomax algebraic model and the q - omega low Reynolds number, two-equation model were used for modeling of turbulence. For the numerical solution of the governing equations a four-stage Runge-Kutta solver was employed. The turbulence model equations were solved using an implicit scheme. Numerical solutions are presented for two-dimensional flow within two vane cascades. The heat transfer results and the pressure distributions were compared with published experimental data. The agreement between the numerical calculations and the experimental values were found to be generally favorable. The position of transition from laminar to turbulent flow was also predicted accurately.

Ameri, A. A.↗