SEARCH · Search NASA

Results for “FLOW EQUATION”

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

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Implicit conservative characteristic modeling schemes for the Euler equations - A new approach

An implicit characteristic-modeling solution scheme for the Euler equations is presented. The scheme does not require the governing equations to be written in characteristic variables or the flux terms to be split into positive and negative contributions. For the two-dimensional problem of a shock wave reflecting from a flat plate, this feature and the simple solution algorithm combine to reduce the computational work per mesh point by 40 percent from that required by a standard, central-difference, implicit, solution algorithm. Application of the method to the quasi-one-dimensional nozzle flow equations for subsonic and supersonic flows without shocks shows the method to be well-conditioned for large time steps.

Wornom, S. F.↗

NASA NTRS1992

Flow prediction over a transport multi-element high-lift system and comparison with flight measurements

Accurate prediction of surface-pressure distributions, merging boundary-layers, and separated-flow regions over multi-element high-lift airfoils is required to design advanced high-lift systems for efficient subsonic transport aircraft. The availability of detailed measurements of pressure distributions and both averaged and time-dependent boundary-layer flow parameters at flight Reynolds numbers is critical to evaluate computational methods and to model the turbulence structure for closure of the flow equations. Several detailed wind-tunnel measurements at subscale Reynolds numbers were conducted to obtain detailed flow information including the Reynolds-stress component. As part of a subsonic-transport high-lift research program, flight experiments are conducted using the NASA-Langley B737-100 research aircraft to obtain detailed flow characteristics for support of computational and wind-tunnel efforts. Planned flight measurements include pressure distributions at several spanwise locations, boundary-layer transition and separation locations, surface skin friction, as well as boundary-layer profiles and Reynolds stresses in adverse pressure-gradient flow.

Vijgen, P. M. H. W.↗

NASA NTRS1993

A controlled variation scheme for convection treatment in pressure-based algorithm

Convection effect and source terms are two primary sources of difficulties in computing turbulent reacting flows typically encountered in propulsion devices. The present work intends to elucidate the individual as well as the collective roles of convection and source terms in the fluid flow equations, and to devise appropriate treatments and implementations to improve our current capability of predicting such flows. A controlled variation scheme (CVS) has been under development in the context of a pressure-based algorithm, which has the characteristics of adaptively regulating the amount of numerical diffusivity, relative to central difference scheme, according to the variation in local flow field. Both the basic concepts and a pragmatic assessment will be presented to highlight the status of this work.

Shyy, Wei↗

NASA NTRS1984

On the use of parametric differentiation to predict weak shock waves

The method of parametric differentiation is shown to be an effective means to transform both the time-linearized and the fully nonlinear transonic small disturbance potential flow equations to a linear system of equations. The globally linearized equations are then solved in phase space by finite difference methods and integral equation methods. The predicted birth, growth, and motion of weak shock waves in two-dimensional transonic flows are compared with existing theory and experiment.

Harris, W. L.↗

NASA NTRS1998

Deterministic Stress Modeling of Hot Gas Segregation in a Turbine

Simulation of unsteady viscous turbomachinery flowfields is presently impractical as a design tool due to the long run times required. Designers rely predominantly on steady-state simulations, but these simulations do not account for some of the important unsteady flow physics. Unsteady flow effects can be modeled as source terms in the steady flow equations. These source terms, referred to as Lumped Deterministic Stresses (LDS), can be used to drive steady flow solution procedures to reproduce the time-average of an unsteady flow solution. The goal of this work is to investigate the feasibility of using inviscid lumped deterministic stresses to model unsteady combustion hot streak migration effects on the turbine blade tip and outer air seal heat loads using a steady computational approach. The LDS model is obtained from an unsteady inviscid calculation. The LDS model is then used with a steady viscous computation to simulate the time-averaged viscous solution. Both two-dimensional and three-dimensional applications are examined. The inviscid LDS model produces good results for the two-dimensional case and requires less than 10% of the CPU time of the unsteady viscous run. For the three-dimensional case, the LDS model does a good job of reproducing the time-averaged viscous temperature migration and separation as well as heat load on the outer air seal at a CPU cost that is 25% of that of an unsteady viscous computation.

Busby, Judy↗

NASA NTRS1992

Computation of unsteady flows over airfoils

Two methods are described for calculating unsteady flows over rapidly pitching airfoils. The first method is based on an interactive scheme in which the inviscid flow is obtained by a panel method. The boundary layer flow is computed by an interactive method that makes use of the Hilbert integral to couple the solutions of the inviscid and viscous flow equations. The second method is based on the solution of the compressible Navier-Stokes equations. The solution of these equations is obtained with an approximately factorized numerical algorithm, and with single block or multiple grids which enable grid embedding to enhance the resolution at isolated flow regions. In addition, the attached flow region can be computed by the numerical solution of compressible boundary layer equations. Unsteady pressure distributions obtained with both methods are compared with available experimental data.

Ekaterinaris, J. A.↗

NASA NTRS1989

Numerical modeling of a vortex stabilized arcjet thruster

A numerical method to solve the equations governing a vortex stabilized arcjet thruster is being developed. The model will allow the effects of swirling flow and geometry on arcjet thruster performance to be determined. The propellant flow equations, i.e., the axisymmetric, thin layer, Navier-Stokes equations, are solved using a Gauss-Siedel line-relaxation procedure. An implicit FTCS method is used to solve the electromagnetic field equations. A grid generation scheme was developed for an arbitrary arcjet geometry. The model will allow the radial and axial components of velocity and current distributions to be determined from a region upstream of the cathode, through the constrictor, to the exit plane of the nozzle.

Keith, Theo G., Jr.↗

NASA NTRS1991

An inverse inviscid method for the design of quasi-three dimensional rotating turbomachinery cascades

A new inverse inviscid method suitable for the design of rotating blade sections lying on an arbitrary axisymmetric stream-surface with varying streamtube width is presented. The geometry of the axisymmetric stream-surface and the streamtube width variation with meridional distance, the number of blades, the inlet flow conditions, the rotational speed and the suction and pressure side velocity distributions as functions of the normalized arc-length are given. The flow is considered irrotational in the absolute frame of reference and compressible. The output of the computation is the blade section that satisfies the above data. The method solves the flow equations on a (phi 1, psi) potential function-streamfunction plane for the velocity modulus, W and the flow angle beta; the blade section shape can then be obtained as part of the physical plane geometry by integrating the flow angle distribution along streamlines. The (phi 1, psi) plane is defined so that the monotonic behavior of the potential function is guaranteed, even in cases with high peripheral velocities. The method is validated on a rotating turbine case and used to design new blades. To obtain a closed blade, a set of closure conditions were developed and referred.

Bonataki, E.↗

NASA NTRS1985

A theoretical basis for extending surface-paneling methods to transonic flow

The surface integral terms in Green's third identity are often used to solve the Prandtl-Glauert (linear potential-flow) equation with panel methods. This can be done, as in the PAN AIR code, for either subsonic or supersonic flow about complete aircraft. The extension to transonic flow is suggested by the volume integral terms of Green's third identity. The mathematical basis for this extension, without the use of body-fitted grids, is presented. Supercritical transonic results computed from a two-dimensional transonic PAN AIR research code demonstrate the method.

Erickson, L. L.↗

NASA NTRS1989

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↗

NASA NTRS1999

Practical Aerodynamic Design Optimization Based on the Navier-Stokes Equations and a Discrete Adjoint Method

The technical details are summarized below: Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. . An efficient surface parameterization based on a free-form deformation technique has been utilized and the resulting codes have been integrated with an optimization package. Lastly, sample optimizations have been shown for inviscid and turbulent flow over an ONERA M6 wing. Drag reductions have been demonstrated by reducing shock strengths across the span of the wing.

Grossman, Bernard↗

NASA NTRS1973

Numerical calculation of the three dimensional transonic flow over a yawed wing.

Results are presented of calculations of the three dimensional steady transonic flow over a finite yawed wing. The full potential flow equation is solved in a transformed coordinate system which permits the boundary conditions to be satisfied exactly. The correct differential properties are enforced by rotating the difference scheme to conform with the flow direction, and fast convergence is assured by simulating a time dependent equation designed to settle quickly to a steady state. Computed lift drag ratios are consistent with the results of wind tunnel tests of a yawed wing conducted by R. T. Jones (1972).-

Jameson, A.↗

NASA NTRS1990

Investigation of vortex breakdown on a delta wing using Euler and Navier-Stokes equations

A numerical investigation of leading-edge vortex breakdown on a delta wing at high angles of attack is presented. The analysis has been restricted to low speed flows on a flat plate wing with sharp leading edges. Both Euler and Navier-Stokes (assuming fully laminar and fully turbulent flows) equations have been used in this study and the results are compared against experimental data. Predictions of vortex breakdown progression with angle of attack with both Euler and Navier-Stokes equations are shown to be consistent with the experimental data. However, the Navier-Stokes predictions show significant improvements in breakdown location at angles of attack where the vortex breakdown approaches the wing apex. The location of the primary vortex and the level of vorticity in the pre-breakdown regions in these flowfield solutions are affected very little by the viscous effects, even though the Navier-Stokes solutions exhibit viscous phenomena such as secondary and tertiary vortices which are not present in the Euler solutions. In the post-breakdown regions, however, the levels of vorticity in the primary vortex have increased differences between the Euler and Navier-Stokes solutions at comparable locations. Navier-Stokes indicates the presence of a secondary vortex even after the primary vortex is burst. The predicted trajectories of the primary vortex are in very good agreement with the test data, the laminar solutions providing the overall best comparison. The Euler shows a small displacement of the primary vortex, relative to experiment, due to the lack of secondary vortices. The turbulent Navier-Stokes, in general, fall between the Euler and laminar solutions. These findings are based on solutions from meshes that are not usually considered fine enough for resolving vortical flows. To further understand the vortex breakdown phenomenon and the fine details of the vortical flow structure, further calculations with the finer meshes and improved turbulence models are necessary.

S. Agrawal↗

NASA NTRS1974

Thermodynamic analysis of an axial-flow laser-powered thruster

A one-dimensional model of a laser-powered thruster is analyzed. The model consists of two flow regions: a laser-heated region and a coaxial region outside the laser region. Different bulk velocities and temperatures are assigned for each region, and a common pressure is assigned for both regions. The one-dimensional inviscid flow equations are solved for each region with the same initial conditions. The results indicate that in order to prevent choked flow while increasing the anthalpy of the inlet gas by a factor of 10, the initial Mach number must be less than 0.14. For hydrogen gas seeded with 10 percent tungsten and entering the thruster tube at 10 atm and 300 K and a laser beam one-half the area of the tube, a tube length between 20 and 50 cm is required to vaporize the particles.

Kascak, A. F.↗

NASA NTRS1977

Computation of transonic boattail flow with separation

The relaxation procedure of South and Jameson for the full potential transonic flow equation was coupled to a modified Reshotko-Tucker integral boundary-layer technique with an empirical model for separated flow. The viscous and inviscid flows were solved iteratively until convergence was obtained. This iterative method was then applied to the subsonic and transonic flow over a series of axisymmetric circular-arc boattails with solid jet plume simulators. Comparisons of theoretical and experimental surface pressures and boattail drag are presented over a free-stream Mach numbers below 0.90. The qualitative variation of boattail drag with free-stream Mach number and boattail angle well into the region of transonic drag rise was correctly predicted; however, the absolute drag levels were significantly underpredicted. For separated flows, the empirical discriminating streamline model gives good results up to a free-stream Mach number of about 0.90 and allows reasonable predictions for shock-induced separation if the proper separation location and separation turning angle are known.

Wilmoth, R. G.↗

NASA NTRS1989

Calculation of steady and unsteady pressures at supersonic speeds with CAP-TSD

A finite difference technique is used to solve the transonic small disturbance flow equation making use of shock capturing to treat wave discontinuities. Thus the nonlinear effects of thickness and angle of attack are considered. Such an approach is made feasible by the development of a new code called CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance), and is based on a fully implicit approximate factorization (AF) finite difference method to solve the time dependent transonic small disturbance equation. The application of the CAP-TSD code to the calculation of low to moderate supersonic steady and unsteady flows is presented. In particular, comparisons with exact linear theory solutions are made for steady and unsteady cases to evaluate shock capturing and other features of the current method. In addition, steady solutions obtained from an Euler code are used to evaluate the small disturbance aspects of the code. Steady and unsteady pressure comparisons are made with measurements for an F-15 wing model and for the RAE tailplane model.

Bennett, Robert M.↗

NASA NTRS1989

Three-dimensional multigrid Navier-Stokes computations for turbomachinery applications

The fully three-dimensional, time-dependent compressible Navier-Stokes equations in cylindrical coordinates are presently used, in conjunction with the multistage Runge-Kutta numerical integration scheme for solution of the governing flow equations, to simulate complex flowfields within turbomechanical components whose pertinent effects encompass those of viscosity, compressibility, blade rotation, and tip clearance. Computed results are presented for selected cascades, emphasizing the code's capabilities in the accurate prediction of such features as airfoil loadings, exit flow angles, shocks, and secondary flows. Computations for several test cases have been performed on a Cray-YMP, using nearly 90,000 grid points.

Subramanian, S. V.↗

NASA NTRS1994

A modified restricted Euler equation for turbulent flows with mean velocity gradients

The restricted Euler equation captures many important features of the behavior of the velocity gradient tensor observed in direct numerical simulations (DNS) of isotropic turbulence. However, in slightly more complex flows the agreement is not good, especially in regions of low dissipation. In this paper, it is demonstrated that the Reynolds-averaged restricted Euler equation violates the balance of mean momentum for virtually all homogeneous turbulent flows with only two major exceptions: isotropic and homogeneously-sheared turbulence. A new model equation which overcomes this shortcoming and is more widely applicable is suggested. This modele is derived from the Navier-Stokes equation with a restricted Euler type approximation made on the fluctuating velocity gradient field. Analytical solutions of the proposed modified restricted Euler equation appear to be difficult to obtain. Hence, a strategy for numerically calculating the velocity gradient tensor is developed. Preliminary calculations tend to indicate that the modified restricted Euler equation captures many important aspects of the behavior of the fluctuating velocity gradients in anisotropic homogeneous turbulence.

Girimaji, Sharath S.↗