Search NASA⌕ Search

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.

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

At least 19 records

Nonunique solutions to the transonic potential flow equation

Steinhoff and Jameson (1981) have shown that within a certain range of angle of attack and freestream Mach number, numerical solutions of the full-potential equation for flow past an airfoil are not unique. This study was mainly concerned with showing that the anomaly is inherent to the partial-differential equation governing the flow and not a result of its discrete representation. Steinhoff and Jameson conjectured that the anomaly may have a physical basis. The present investigation has two objectives. Results are to be presented which indicate that the anomaly is due to a breakdown in the potential approximation, rather than a phenomenon associated with the inviscid flowfield. The second objective is to show that the lift coefficient, predicted by the potential equation, is a smooth but multivalued function of the angle of attack.

Salas, M. D.↗

Thin airfoil theory based on approximate solution of the transonic flow equation

A method is presented for the approximate solution of the nonlinear equations transonic flow theory. Solutions are found for two-dimensional flows at a Mach number of 1 and for purely subsonic and purely supersonic flows. Results are obtained in closed analytic form for a large and significant class of nonlifting airfoils. At a Mach number of 1 general expressions are given for the pressure distribution on an airfoil of specified geometry and for the shape of an airfoil having a prescribed pressure distribution. Extensive comparisons are made with available data, particularly for a Mach number of 1, and with existing solutions.

WING-SECTION THEORY↗

A rule for selecting analytical boundary conditions for the conservative quasi-one-dimensional nozzle flow equations

For the one-dimensional nozzle flow equations, the number of analytical boundary conditions that may be applied corresponds to the number of exterior characteristics reaching the boundary points. Since the slopes of the characteristics are the characteristic speeds u, u+c, and u-c, one concludes that two analytical boundary conditions are required at a subsonic inflow boundary and one at a subsonic outflow boundary. However, using this guideline alone can lead to nonunique solutions in subsonic cases. This paper presents a rule for selecting analytical boundary conditions for the dependent variables which will yield unique solutions if they exist.

Wornom, S. F.↗

Accurate solutions, parameter studies and comparisons for the Euler and potential flow equations

Parameter studies are conducted using the Euler and potential flow equation models for unsteady and steady flows in both two and three dimensions. The Euler code is an implicit, upwind, finite volume code which uses the Van Leer method of flux-vector-splitting which has been recently extended for use on dynamic meshes and maintain all the properties of the original splitting. The potential flow code is an implicit, finite difference method for solving the transonic small disturbance equations and incorporates both entropy and vorticity corrections into the solution procedures thereby extending its applicability into regimes where shock strength normally precludes its use. Parameter studies resulting in benchmark type calculations include the effects of spatial and temporal refinement, spatial order of accuracy, far field boundary conditions for steady flow, frequency of oscillation, and the use of subiterations at each time step to reduce linearization and factorization errors. Comparisons between Euler and potential flows results are made as well as with experimental data where available.

Anderson, W. Kyle↗

Accurate solutions, parameter studies and comparisons for the Euler and potential flow equations

Parameter studies are conducted using the Euler and potential flow equation models for steady and unsteady flows in both two and three dimensions. The Euler code is an implicit, upwind, finite volume code which uses the Van Leer method of flux vector splitting which has been recently extended for use on dynamic meshes and maintain all the properties of the original splitting. The potential flow code is an implicit, finite difference method for solving the transonic small disturbance equations and incorporates both entropy and vorticity corrections into the solution procedures thereby extending its applicability into regimes where shock strength normally precludes its use. Parameter studies resulting in benchmark type calculations include the effects of spatial and temporal refinement, spatial order of accuracy, far field boundary conditions for steady flow, frequency of oscillation, and the use of subiterations at each time step to reduce linearization and factorization errors. Comparisons between Euler and potential flow results are made, as well as with experimental data where available.

Anderson, W. Kyle↗

A Highly Accurate Technique for the Treatment of Flow Equations at the Polar Axis in Cylindrical Coordinates using Series Expansions

Numerical methods for solving the flow equations in cylindrical or spherical coordinates should be able to capture the behavior of the exact solution near the regions where the particular form of the governing equations is singular. In this work we focus on the treatment of these numerical singularities for finite-differences methods by reinterpreting the regularity conditions developed in the context of pseudo-spectral methods. A generally applicable numerical method for treating the singularities present at the polar axis, when nonaxisymmetric flows are solved in cylindrical, coordinates using highly accurate finite differences schemes (e.g., Pade schemes) on non-staggered grids, is presented. Governing equations for the flow at the polar axis are derived using series expansions near r=0. The only information needed to calculate the coefficients in these equations are the values of the flow variables and their radial derivatives at the previous iteration (or time) level. These derivatives, which are multi-valued at the polar axis, are calculated without dropping the accuracy of the numerical method using a mapping of the flow domain from (0,R)*(0,2pi) to (-R,R)*(0,pi), where R is the radius of the computational domain. This allows the radial derivatives to be evaluated using high-order differencing schemes (e.g., compact schemes) at points located on the polar axis. The proposed technique is illustrated by results from simulations of laminar-forced jets and turbulent compressible jets using large eddy simulation (LES) methods. In term of the general robustness of the numerical method and smoothness of the solution close to the polar axis, the present results compare very favorably to similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r=0, following the method proposed recently by Mohseni and Colonius (1). Extension of the method described here for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical or spherical coordinates with finite-differences schemes of various level of accuracy is immediate.

Constantinescu, George S.↗

A New Method for Accurate Treatment of Flow Equations in Cylindrical Coordinates Using Series Expansions

The motivation of this work is the ongoing effort at the Center for Turbulence Research (CTR) to use large eddy simulation (LES) techniques to calculate the noise radiated by jet engines. The focus on engine exhaust noise reduction is motivated by the fact that a significant reduction has been achieved over the last decade on the other main sources of acoustic emissions of jet engines, such as the fan and turbomachinery noise, which gives increased priority to jet noise. To be able to propose methods to reduce the jet noise based on results of numerical simulations, one first has to be able to accurately predict the spatio-temporal distribution of the noise sources in the jet. Though a great deal of understanding of the fundamental turbulence mechanisms in high-speed jets was obtained from direct numerical simulations (DNS) at low Reynolds numbers, LES seems to be the only realistic available tool to obtain the necessary near-field information that is required to estimate the acoustic radiation of the turbulent compressible engine exhaust jets. The quality of jet-noise predictions is determined by the accuracy of the numerical method that has to capture the wide range of pressure fluctuations associated with the turbulence in the jet and with the resulting radiated noise, and by the boundary condition treatment and the quality of the mesh. Higher Reynolds numbers and coarser grids put in turn a higher burden on the robustness and accuracy of the numerical method used in this kind of jet LES simulations. As these calculations are often done in cylindrical coordinates, one of the most important requirements for the numerical method is to provide a flow solution that is not contaminated by numerical artifacts. The coordinate singularity is known to be a source of such artifacts. In the present work we use 6th order Pade schemes in the non-periodic directions to discretize the full compressible flow equations. It turns out that the quality of jet-noise predictions using these schemes is especially sensitive to the type of equation treatment at the singularity axis. The objective of this work is to develop a generally applicable numerical method for treating the singularities present at the polar axis, which is particularly suitable for highly accurate finite-differences schemes (e.g., Pade schemes) on non-staggered grids. The main idea is to reinterpret the regularity conditions developed in the context of pseudo-spectral methods. A set of exact equations at the singularity axis is derived using the appropriate series expansions for the variables in the original set of equations. The present treatment of the equations preserves the same level of accuracy as for the interior scheme. We also want to point out the wider utility of the method, proposed here in the context of compressible flow equations, as its extension for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical coordinates with finite-differences schemes of various level of accuracy is straightforward. The robustness and accuracy of the proposed technique is assessed by comparing results from simulations of laminar forced-jets and turbulent compressible jets using LES with similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r = 0.

Constantinescu, G.S.↗

Thin airfoil theory based on approximate solution of the transonic flow equation

A method is presented for the approximate solution of the nonlinear equations of transonic flow theory. Solutions are found for two-dimensional flows at a Mach number of 1 and for purely subsonic and purely supersonic flows. Results are obtained in closed analytic form for a large and significant class of nonlifting airfoils. At a Mach number of 1 general expressions are given for the pressure distribution on an airfoil of specified geometry and for the shape of an airfoil having a prescribed pressure distribution. Extensive comparisons are made with available data, particularly for a Mach number of 1, and with existing solutions.

Spreiter, John R↗

Multigrid acceleration of the isenthalpic form of the compressible flow equations

A numerical method for solving the isenthalpic form of the governing equations for compressible inviscid flows was developed. The method is based on the concept of flux vector splitting in its implicit form and was tested on several demanding configurations. Time marching to steady state was accelerated by the implementation of the multigrid procedure which very effectively increased the rate of convergence. High quality steady-state results were obtained for various test cases and required only short computational times due to the relative efficiency of the basic method.

Melson, N. Duane↗

Multigrid acceleration of the isenthalpic form of the compressible flow equations

A numerical method for solving the isenthalpic form of the governing equations for compressible inviscid flows was developed. The method is based on the concept of flux vector splitting in its implicit form and was tested on several demanding configurations. Time marching to steady state was accelerated by the implementation of the multigrid procedure which very effectively increased the rate of convergence. High quality steady-state results were obtained for various test cases and required only short computational times due to the relative efficiency of the basic method.

Melson, N. Duane↗

Fluid-structural interactions using Navier-Stokes flow equations coupled with shell finite element structures

A computational procedure is presented to study fluid-structural interaction problems for three-dimensional aerospace structures. The flow is modeled using the three-dimensional unsteady Euler/Navier-Stokes equations and solved using the finite-difference approach. The three dimensional structure is modeled using shell/plate finite-element formulation. The two disciplines are coupled using a domain decomposition approach. Accurate procedures both in time and space are developed to combine the solutions from the flow equations with those of the structural equations. Time accuracy is maintained using aeroelastic configuration-adaptive moving grids that are computed every time step. The work done by aerodynamic forces due to structural deformations is preserved using consistent loads. The present procedure is validated by computing the aeroelastic response of a wing and comparing with experiment. Results are illustrated for a typical wing-body configuration.

Guruswamy, Guru P.↗

A sonic flow equation for electric arc jets

The relationship between total enthalpy and the flow parameters of two types of electric arc jets is discussed. A simple equation for the supersonic arc jet, based on ARCFLO code calculations for mass-average total enthalpy, is presented in terms of a sonic flow parameter. At enthalpies greater than about 25 MJ/kg, this equation shows better agreement with experimental arc jet data than a previous equation.

Shepard, Charles E.↗

Slender-Body Theory Based On Approximate Solution of the Transonic Flow Equation

Approximate solution of the nonlinear equations of the small disturbance theory of transonic flow are found for the pressure distribution on pointed slender bodies of revolution for flows with free-stream, Mach number 1, and for flows that are either purely subsonic or purely supersonic. These results are obtained by application of a method based on local linearization that was introduced recently in the analysis of similar problems in two-dimensional flows. The theory is developed for bodies of arbitrary shape, and specific results are given for cone-cylinders and for parabolic-arc bodies at zero angle of attack. All results are compared either with existing theoretical results or with experimental data.

Spreiter, John R.↗