SEARCH · Engineering Papers

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 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Turbulent flow calculations using unstructured and adaptive meshes

A method of efficiently computing turbulent compressible flow over complex two dimensional configurations is presented. The method makes use of fully unstructured meshes throughout the entire flow field, thus enabling the treatment of arbitrarily complex geometries and the use of adaptive meshing techniques throughout both viscous and inviscid regions of the flow field. Mesh generation is based on a locally mapped Delaunay technique in order to generate unstructured meshes with highly-stretched elements in the viscous regions. The flow equations are discretized using a finite element Navier-Stokes solver, and rapid convergence to steady-state is achieved using an unstructured multigrid algorithm. Turbulence modeling is performed using an inexpensive algebraic model, implemented for use on unstructured and adaptaive meshes. Compressible turbulent flow solutions about multiple-element airfoil geometries are computed and compared with expeirmental data.

Mavriplis, Dimitri J.↗

NASA NTRS1996

Modeling the effect of unsteady flows on the time average flow field of a blade row embedded in an axial flow multistage turbomachine

The role of unsteady flow processes in establishing the performance of axial flow turbomachinery was investigated. The development of the flow model, as defined by the time average flow equations associated with the passage of the blade row embedded in a multistage configuration, is presented. The mechanisms for unsteady momentum and energy transport is outlined. The modeling of the unsteady momentum and energy transport are discussed. The procedure for simulating unsteady multistage turbomachinery flows is described.

Adamczyk, John J.↗

NASA NTRS1977

Analytical study of viscous effects on transonic flow over boattail nozzles

An analytical study has been made of the effects of viscous/inviscid interactions on the transonic flow over boattail nozzles. The theoretical method couples a relaxation solution of the full potential transonic flow equations with a conventional boundary-layer solution to account for displacement-thickness effects. Surface pressures calculated on circular-arc boattails with solid jet-plume simulators show good agreement with experiment for free-stream Mach numbers less than 0.90. For separated flows, an empirical discriminating streamline model of the separation bubble gives good results until the onset of shock-induced separation, which occurs typically at a Mach number of about 0.90.

Wilmoth, R. G.↗

NASA NTRS1993

Preliminary experimental results for a cryogenic brush seal configuration

Preliminary fluid nitrogen flow data are reported for a five-brush, ceramic-coated-rub-runner brush seal system, where the brushes and the rub runner were placed at each end of a centrally pressurized multifunction tester ('back-to-back' set of brushes) and tested at rotor speeds of 0, 10, 18, and 22.5 krpm. After testing, both the brushes and the ceramic-coated rub runner appeared pristine. The coating withstood both the thermomechanical and dynamic loadings with minor wear track scarring. The bristle tips showed some indication of material shearing (smearing) wear. The Ergun porous flow equation was applied to the brush seal data. The Ergun relation, which required heuristic information to characterize the coefficients, fit the gaseous data but was in poor agreement with the fluid results. The brush seal exit conditions were two phase. Two-phase, choked-flow design charts were applied but required one data point at each rotor speed to define the (C(sub f)A x Constant) flow and area coefficients. Reasonable agreement between prediction and data was found, as expected, but such methods are not to be construed as two-phase-flow brush seal analyses.

Carlile, J. A.↗

NASA NTRS1994

Benchmark problems in computational aeroacoustics

A recent directive at NASA Langley is aimed at numerically predicting principal noise sources. During my summer stay, I worked with high-order ENO code, developed by Dr. Harold Atkins, for solving the unsteady compressible Navier-Stokes equations, as it applies to computational aeroacoustics (CAA). A CAA workshop, composed of six categories of benchmark problems, has been organized to test various numerical properties of code. My task was to determine the robustness of Atkins' code for these test problems. In one category, we tested the nonlinear wave propagation of the code for the one-dimensional Euler equations, with initial pressure, density, and velocity conditions. Using freestream boundary conditions, our results were plausible. In another category, we solved the linearized two-dimensional Euler equations to test the effectiveness of radiation boundary conditions. Here we utilized MAPLE to compute eigenvalues and eigenvectors of the Jacobian given variable and flux vectors. We experienced a minor problem with inflow and outflow boundary conditions. Next, we solved the quasi one dimensional unsteady flow equations with an incoming acoustic wave of amplitude 10(exp -6). The small amplitude sound wave was incident on a convergent-divergent nozzle. After finding a steady-state solution and then marching forward, our solution indicated that after 30 periods the acoustic wave had dissipated (a period is time required for sound wave to traverse one end of nozzle to other end).

Porter-Locklear, Freda↗

NASA NTRS2004

Procedure for Determining One-Dimensional Flow Distributions in Arbitrarily Connected Passages Without the Influence of Pumping

A calculation procedure is presented which allows the one-dimensional determination of flow distributions in arbitrarily connected (branching) flow passages having multiple inlets and exits. The procedure uses an adaptation of the finite element technique, iteratively coupled with an accurate one-dimensional flow solver. The procedure eliminates the usual restrictions inherent with finite element flow calculations. Unlike existing one-dimensional methods, which require simplifications to the flow equations (uncoupling the momentum and energy equations), to allow for arbitrary branching and multiple inlets and exits, the only limitation of the described methodology is that, at present, it can only accommodate non-rotating configurations (no pumping effects). The calculation procedure is robust, and will always converge for physically possible flow. The procedure is described, and its use is illustrated by an example.

Meitner, Peter L.↗

NASA NTRS1974

Development of inverse inviscid transonic solution methods

A numerical method suitable for the analysis and/or design of supercritical transonic airfoils is reported. In order to achieve accuracy, the method utilizes the full inviscid potential flow equations; and in order to remain simple it solves the problem in a stretched Cartesian grid system. The resulting computer program has several advantages over others of its type -- its use in either the direct analysis mode in which the airfoil shape is prescribed and the flow field and surface pressures are determined, or in the inverse mode in which the surface pressures are given and the airfoil shape and flow field are computed. Other advantages of the program include its use in a design program, the rotated finite difference scheme and its determination of the airfoil shape simultaneously with the flow field relaxation solution.

Carlson, L. A.↗

NASA NTRS1996

Transport Coefficients in weakly compressible turbulence

A theory of transport coefficients in weakly compressible turbulence is derived by applying Yoshizawa's two-scale direct interaction approximation to the compressible equations of motion linearized about a state of incompressible turbulence. The result is a generalization of the eddy viscosity representation of incompressible turbulence. In addition to the usual incompressible eddy viscosity, the calculation generates eddy diffusivities for entropy and pressure, and an effective bulk viscosity acting on the mean flow. The compressible fluctuations also generate an effective turbulent mean pressure and corrections to the speed of sound. Finally, a prediction unique to Yoshizawa's two-scale approximation is that terms containing gradients of incompressible turbulence quantities also appear in the mean flow equations. The form these terms take is described.

Rubinstein, Robert↗

NASA NTRS1978

Computational transonic flow about realistic aircraft configurations

A numerical method has been developed to compute transonic flows about realistic wing-fuselage configurations. The finite difference scheme employs an improved small disturbance flow equation. A unique grid embedding technique, which was heretofore applied to airfoils and wings, has been extended to include the treatment of both body and wing-body shapes. The resulting high-density mesh is shown to be a valuable asset in resolving details of the three-dimensional flow. A mathematical modeling system is used to process arbitrary fuselage geometries for body boundary conditions. Correlations with experimental data for simple isolated bodies, an isolated fuselage, and wing-fuselage combinations are included.

Boppe, C. W.↗

NASA NTRS1978

Forebody and afterbody solutions of the Navier-Stokes equations for supersonic flow over blunt bodies in a generalized orthogonal coordinate system

A coordinate transformation, which can approximate many different two-dimensional and axisymmetric body shapes with an analytic function, is used as a basis for solving the Navier-Stokes equations for the purpose of predicting 0 deg angle of attack supersonic flow fields. The transformation defines a curvilinear, orthogonal coordinate system in which coordinate lines are perpendicular to the body and the body is defined by one coordinate line. This system is mapped in to a rectangular computational domain in which the governing flow field equations are solved numerically. Advantages of this technique are that the specification of boundary conditions are simplified and, most importantly, the entire flow field can be obtained, including flow in the wake. Good agreement has been obtained with experimental data for pressure distributions, density distributions, and heat transfer over spheres and cylinders in supersonic flow. Approximations to the Viking aeroshell and to a candidate Jupiter probe are presented and flow fields over these shapes are calculated.

Gnoffo, P. A.↗

NASA NTRS1992

Heat transfer in rotating serpentine passages with trips skewed to the flow

Johnson, B. V.↗

NASA NTRS1986

A Constitutive Model for the Inelastic Multiaxial Cyclic Response of a Nickel Base Superalloy Rene 80

The objective was to develop unified constitutive equations which can model a variety of nonlinear material phenomena observed in Rene 80 at elevated temperatures. A constitutive model was developed based on back stress and drag stress. The tensorial back stress was used to model directional effects; whereas, the scalar drag stress was used to model isotropic effects and cyclic hardening or softening. A flow equation and evolution equations for the state variables were developed in multiaxial form. Procedures were developed to generate the material parameters. The model predicted very well the monotonic tensile, cyclic, creep, and stress relaxation behavior of Rene 80 at 982 C. The model was then extended to 871, 760, and 538 C. It was shown that strain rate dependent behavior at high temperatures and strain rate independent behavior at the lower temperatures could be predicted very well. A large number of monotonic tensile, creep, stress relation, and cyclic experiments were predicted. The multiaxial capabilities of the model were verified extensively for combined tension/torsion experiments. The prediction of the model agreed very well for proportional, nonproportional, and pure shear cyclic loading conditions at 982 and 871 C.

Ramaswamy, V. G.↗

NASA NTRS1992

A finite-difference approximate-factorization algorithm for solution of the unsteady transonic small-disturbance equation

A time-accurate approximate-factorization (AF) algorithm is described for solution of the three-dimensional unsteady transonic small-disturbance equation. The AF algorithm consists of a time-linearization procedure coupled with a subiteration technique. The algorithm is the basis for the Computational Aeroelasticity Program-Transonic Small Disturbance (CAP-TSD) computer code, which was developed for the analysis of unsteady aerodynamics and aeroelasticity of realistic aircraft configurations. The paper describes details on the governing flow equations and boundary conditions, with an emphasis on documenting the finite-difference formulas of the AF algorithm.

Batina, John T.↗

NASA NTRS1972

An improved design method and experimental performance of two dimensional curved wall diffusers

A computer design program was developed to incorporate the suction slots in solving the potential flow equations with prescribed boundary conditions. Using the contour generated from this program two Griffith diffusers were tested having area ratios AR = 3 and 4. The inlet Reynolds number ranged from 600,000 to 7 million. It was found that the slot suction required for metastable operation depends on the sidewall suction applied. Values of slot suction of 8% of the inlet flow rate was required for AR = 4 with metastable condition, provided that enough sidewall suction was applied. For AR = 3, the values of slot suction was about 25% lower than those required for AR = 4. For nearly all unseparated test runs, the effectiveness was 100% and the exit flow was uniform. In addition to the Griffith diffusers, dump and cusp diffusers of comparable area ratios were built and tested. The results obtained from these diffusers were compared with those of the Griffith diffusers. Flow separation occurred in all test runs with the dump and cusp diffusers.

Yang, T.↗

NASA NTRS1988

A constitutive material model for nonlinear finite element structural analysis using an iterative matrix approach

A unified numerical method for the integration of stiff time dependent constitutive equations is presented. The solution process is directly applied to a constitutive model proposed by Bodner. The theory confronts time dependent inelastic behavior coupled with both isotropic hardening and directional hardening behaviors. Predicted stress-strain responses from this model are compared to experimental data from cyclic tests on uniaxial specimens. An algorithm is developed for the efficient integration of the Bodner flow equation. A comparison is made with the Euler integration method. An analysis of computational time is presented for the three algorithms.

Koenig, Herbert A.↗

NASA NTRS1994

A flux splitting scheme with high-resolution and robustness for discontinuities

A flux splitting scheme is proposed for the general nonequilibrium flow equations with an aim at removing numerical dissipation of Van-Leer-type flux-vector splittings on a contact discontinuity. The scheme obtained is also recognized as an improved Advection Upwind Splitting Method (AUSM) where a slight numerical overshoot immediately behind the shock is eliminated. The proposed scheme has favorable properties: high-resolution for contact discontinuities; conservation of enthalpy for steady flows; numerical efficiency; applicability to chemically reacting flows. In fact, for a single contact discontinuity, even if it is moving, this scheme gives the numerical flux of the exact solution of the Riemann problem. Various numerical experiments including that of a thermo-chemical nonequilibrium flow were performed, which indicate no oscillation and robustness of the scheme for shock/expansion waves. A cure for carbuncle phenomenon is discussed as well.

Wada, Yasuhiro↗

NASA NTRS2006

Reentry-Vehicle Shape Optimization Using a Cartesian Adjoint Method and CAD Geometry

A DJOINT solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (e.g., geometric parameters that control the shape). Classic aerodynamic applications of gradient-based optimization include the design of cruise configurations for transonic and supersonic flow, as well as the design of high-lift systems. are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric computer-aided design (CAD). In previous work on Cartesian adjoint solvers, Melvin et al. developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the two-dimensional Euler equations using a ghost-cell method to enforce the wall boundary conditions. In Refs. 18 and 19, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm were the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The accuracy of the gradient computation was verified using several three-dimensional test cases, which included design variables such as the free stream parameters and the planform shape of an isolated wing. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Factors under consideration include the computation of mesh sensitivities that provide a reliable approximation of the objective function gradient, as well as the computation of surface shape sensitivities based on a direct-CAD interface. We present detailed gradient verification studies and then focus on a shape optimization problem for an Apollo-like reentry vehicle. The goal of the optimization is to enhance the lift-to-drag ratio of the capsule by modifying the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design.

Nemec, Marian↗

NASA NTRS1992

Heat transfer in rotating serpentine passages with trips skewed to the flow

Experiments were conducted to determine the effects of buoyancy and Coriolis forces on heat transfer in turbine blade internal coolant passages. The experiments were conducted with a large scale, multi-pass, heat transfer model with both radially inward and outward flow. Trip strips, skewed at 45 deg to the flow direction, were machined on the leading and trailing surfaces of the radial coolant passages. An analysis of the governing flow equations showed that four parameters influence the heat transfer in rotating passages: coolant-to-wall temperature, rotation number, Reynolds number, and radius-to-passage hydraulic diameter ratio. The first three of these four parameters were varied over ranges which are typical of advanced gas turbine engine operating conditions. Results were correlated and compared to previous results from similar stationary and rotating models with smooth walls and with trip strips normal to the flow direction. The heat transfer coefficients on surfaces, where the heat transfer decreased with rotation and buoyancy, decreased to as low as 40 percent of the value without rotation. However, the maximum values of the heat transfer coefficients with high rotation were only slightly above the highest levels previously obtained with the smooth wall models. It was concluded that (1) both Coriolis and buoyancy effects must be considered in turbine blade cooling designs with trip strips, (2) the effects of rotation are markedly different depending upon the flow direction, and (3) the heat transfer with skewed trip strips is less sensitive to buoyancy than the heat transfer in models with either smooth or normal trips. Therefore, skewed trip strips rather than normal trip strips are recommended and geometry-specific tests are required for accurate design information.

Johnson, B. V.↗