SEARCH · Engineering Papers

Results for “Stokes equations”

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 811 records · Page 45

All sources NASA NTRS DOE OSTI

NASA NTRS2010

Development and Application of Agglomerated Multigrid Methods for Complex Geometries

We report progress in the development of agglomerated multigrid techniques for fully un- structured grids in three dimensions, building upon two previous studies focused on efficiently solving a model diffusion equation. We demonstrate a robust fully-coarsened agglomerated multigrid technique for 3D complex geometries, incorporating the following key developments: consistent and stable coarse-grid discretizations, a hierarchical agglomeration scheme, and line-agglomeration/relaxation using prismatic-cell discretizations in the highly-stretched grid regions. A signi cant speed-up in computer time is demonstrated for a model diffusion problem, the Euler equations, and the Reynolds-averaged Navier-Stokes equations for 3D realistic complex geometries.

Nishikawa, Hiroaki↗

NASA NTRS2022

Boundary Layer Instabilities Over a Cone-Cylinder-Flare Model at Mach 6

Computations are performed to investigate the boundary-layer instabilities over a sharp cone-cylinder-flare model at zero degrees angle of attack. The model geometry and the flow conditions are selected to match the experiments conducted in the Boeing/AFOSR Mach 6 Quiet Tunnel (BAM6QT) at Purdue University. The geometry consists of a nominally sharp 5-degree half-angle cone, followed by a cylindrical segment and then a 10-degree flare. An axisymmetric separation bubble is generated as a result of the laminar shock/boundary-layer interaction in the cylinder-flare region. The comparison of the laminar flow solution and the schlieren images shows a remarkable agreement between the respective locations of both the boundarylayer edge and the reattachment shock. The predicted heat flux distribution is also in agreement with the measured values downstream of the reattachment location. The analysis of convective and global instabilities is performed for flare half angles equal to 8, 10, and 12 degrees and nosetip radii equal to 0.1, 1, and 5 mm. The linear amplification of first and second Mack mode instabilities that begin to amplify in the cone region are computed with a combination of the parabolized stability equations (PSE) and the harmonic linearized Navier-Stokes equations (HLNSE). The predicted frequency spectra of the surface pressure fluctuations associated with both planar and oblique instability waves are compared with the measured spectra at the various locations of the PCB and Kulite sensors. The comparison shows that the computational analysis captures the distinct lobes within the disturbance amplification spectra measured in the experiments, but some differences in amplification characteristics are noted at low frequencies. Overall, the oblique disturbances are found to be more amplified than the planar disturbances. To our knowledge, this represents the first successful comparison between convective instability analysis and measured surface pressure fluctuations for a hypersonic configuration with a separation bubble. Finally, the global instability analysis shows that the laminar flow becomes supercritical for flare half angles larger than 8 degrees. The unstable global mode for the experimental configuration of a 10 degrees flare and a sharp nosetip cone corresponds to a stationary three-dimensional disturbance that is concentrated in the recirculation region and achieves its maximum growth rate for an azimuthal wavenumber of 5.

Hydrodynamic instability↗

NASA NTRS1982

Finite analytic numerical solution of axisymmetric Navier-Stokes and energy equations

Convective heat transfer for steady state laminar flow in axisymmetric coordinates is considered. Numerical solutions for flow pattern and temperature distribution are obtained by the finite analytic numerical method applied to the Navier-Stokes equations expressed in terms of vorticity and stream function, and the energy equation. The finite analytic numerical method differs from other numerical methods in that it utilizes a local analytic solution in an element of the problem to construct the total numerical solution. Finite analytic solutions of vorticity, stream function, temperature and heat transfer coefficients for flow with Reynolds number of 5.0, 100.0, 1000.0, and 2000.0, and Prandtl number of 0.1, 1.0, and 10.0 with uniform grid sizes are reported for axisymmetric pipe with sudden expansion and contraction. The wall temperature is considered to be isothermal and differs from the inlet temperature. It is shown that the finite analytic solution is stable, converges rapidly, and simulates the convection of fluid flow accurately since the local analytic solution is capable of simulating automatically the influence of skewed convection through the element boundary on the interior nodal values thereby minimizing the false numerical diffusion.

Chen, C.-J.↗

NASA NTRS1985

Turbulent solutions of equations of fluid motion

Some turbulent solutions of the unaveraged Navier-Stokes equations (equations of fluid motion) are reviewed. Those equations are solved numerically in order to study the nonlinear physics of incompressible turbulent flow. The three components of the mean-square velocity fluctuations are initially equal for the conditions chosen. The resulting solutions show characteristics of turbulence, such as the linear and nonlinear excitation of small-scale fluctuations. For the stronger fluctuations the initially nonrandom flow develops into an apparently random turbulence. The cases considered include turbulence that is statistically homogeneous or inhomogeneous and isotropic or anisotropic. A statistically steady-state turbulence is obtained by using a spatially periodic body force. Various turbulence processes, including the transfer of energy between eddy sizes and between directional components and the production, dissipation, and spatial diffusion of turbulence, are considered. It is concluded that the physical processes occurring in turbulence can be profitably studied numerically.

Deissler, R. G.↗

NASA NTRS1992

Runge-Kutta methods combined with compact difference schemes for the unsteady Euler equations

Recent development using compact difference schemes to solve the Navier-Stokes equations show spectral-like accuracy. A study was made of the numerical characteristics of various combinations of the Runge-Kutta (RK) methods and compact difference schemes to calculate the unsteady Euler equations. The accuracy of finite difference schemes is assessed based on the evaluations of dissipative error. The objectives are reducing the numerical damping and, at the same time, preserving numerical stability. While this approach has tremendous success solving steady flows, numerical characteristics of unsteady calculations remain largely unclear. For unsteady flows, in addition to the dissipative errors, phase velocity and harmonic content of the numerical results are of concern. As a result of the discretization procedure, the simulated unsteady flow motions actually propagate in a dispersive numerical medium. Consequently, the dispersion characteristics of the numerical schemes which relate the phase velocity and wave number may greatly impact the numerical accuracy. The aim is to assess the numerical accuracy of the simulated results. To this end, the Fourier analysis is to provide the dispersive correlations of various numerical schemes. First, a detailed investigation of the existing RK methods is carried out. A generalized form of an N-step RK method is derived. With this generalized form, the criteria are derived for the three and four-step RK methods to be third and fourth-order time accurate for the non-linear equations, e.g., flow equations. These criteria are then applied to commonly used RK methods such as Jameson's 3-step and 4-step schemes and Wray's algorithm to identify the accuracy of the methods. For the spatial discretization, compact difference schemes are presented. The schemes are formulated in the operator-type to render themselves suitable for the Fourier analyses. The performance of the numerical methods is shown by numerical examples. These examples are detailed. described. The third case is a two-dimensional simulation of a Lamb vortex in an uniform flow. This calculation provides a realistic assessment of various finite difference schemes in terms of the conservation of the vortex strength and the harmonic content after travelling a substantial distance. The numerical implementation of Giles' non-refelctive equations coupled with the characteristic equations as the boundary condition is discussed in detail. Finally, the single vortex calculation is extended to simulate vortex pairing. For the distance between two vortices less than a threshold value, numerical results show crisp resolution of the vortex merging.

Yu, Sheng-Tao↗

NASA NTRS1988

Application of unsteady aerodynamic methods for transonic aeroelastic analysis

Aerodynamic methods for aeroelastic analysis are applied to various flow problems. These methods include those that solve the three dimensional transonic small disturbance (TSD) potential equation, the two dimensional (2-D) full potential (FP) equation, and the 2-D thin layer Navier-Stokes equations. Flutter analysis performed using TSD aerodynamics show that such methods can be used to analyze some aeroelastic phenomena. For thicker bodies and larger amplitude motions, a nonisentropic FP method is presented. The unsteady FP equation is modified to model the entropy jumps across shock waves. The conservative form of the modified equation is solved in generalized coordinates using an implicit, approximate factorization method. Pressures calculated on the NLR 7301 and NACA 64A010A airfoils using the nonisentropic FP method are presented. It is shown that modeling shock generated entropy extends the range of validity of the FP method. A Navier-Stokes code is correlated with pressures measured on a supercritical airfoil at transonic speeds. When corrections are made for wind tunnel wall effects, the calculations correlate well with the measured data.

Whitlow, Woodrow, Jr.↗

NASA NTRS1982

A viscous-inviscid interaction algorithm for three-dimensional turbulent subsonic aerodynamic juncture region flow

An order-of magnitude analysis of the subsonic, three-dimensional, steady time-averaged Navier-Stokes equations for semibounded aerodynamic juncture geometries yields the parabolic Navier-Stokes simplification. The numerical solution of the resultant pressure Poisson equation is cast into complementary and particular parts, yielding an iterative interaction algorithm with an exterior three-dimensional potential flow solution. A parabolic transverse momentum equation set is constructed, wherein robust enforcement first-order continuity effects is accomplished using a penalty differential constraint concept within a finite element solution algorithm. A Reynolds stress constitutive equation, with low turbulence Reynolds number wall functions, is employed for closure, using parabolic forms of the two-equation turbulent kinetic energy-dissipation equation system. The algorithm is employed to predict the three-dimensional turbulent flowfield in the juncture region formed by two intersecting parabolic arcs at a free-stream Mach-number of 0.08, and Re/C = 600,000/m.

Baker, A. J.↗

NASA NTRS1982

The CMC:3DPNS computer program for prediction of three-dimensional, subsonic, turbulent aerodynamic juncture region flow. Volume 1: Theoretical

An order-of-magnitude analysis of the subsonic three dimensional steady time averaged Navier-Stokes equations, for semibounded aerodynamic juncture geometries, yields the parabolic Navier-Stokes simplification. The numerical solution of the resultant pressure Poisson equation is cast into complementary and particular parts, yielding an iterative interaction algorithm with an exterior three dimensional potential flow solution. A parabolic transverse momentum equation set is constructed, wherein robust enforcement of first order continuity effects is accomplished using a penalty differential constraint concept within a finite element solution algorithm. A Reynolds stress constitutive equation, with low turbulence Reynolds number wall functions, is employed for closure, using parabolic forms of the two-equation turbulent kinetic energy-dissipation equation system. Numerical results document accuracy, convergence, and utility of the developed finite element algorithm, and the CMC:3DPNS computer code applied to an idealized wing-body juncture region. Additional results document accuracy aspects of the algorithm turbulence closure model.

Baker, A. J.↗

NASA NTRS1993

Recent developments in high order K-exact reconstruction on unstructured meshes

This paper presents recent improvements in high-order K-exact reconstruction on unstructured meshes. The new reconstruction procedures are incorporated into a basic upwind finite-volume scheme suitable for solving scalar advection-diffusion equations as well as the Euler and Navier-Stokes equations. Numerical calculations are performed comparing the present method with lower order accurate reconstruction procedures (piecewise constant and piecewise linear) and various competing technologies such as the fluctuation splitting method of Roe (1987) and Deconinck et al. (1992) and a system-variant of the streamline diffusion Petrov-Galerkin method developed by Hansbo (1991) and Hansbo and Johnson (1991). Five test problems are used in the numerical comparisons: scalar circular advection, transonic and supersonic Euler flow, laminar boundary-layer flow, and general compressible Navier-Stokes flow.

Barth, Timothy J.↗

NASA NTRS1980

The dynamics of oceanic fronts. I - The Gulf Stream

The establishment and maintenance of the mean hydrographic properties of large-scale density fronts in the upper ocean is considered. The dynamics is studied by posing an initial value problem starting with a near-surface discharge of buoyant water with a prescribed density deficit into an ambient stationary fluid of uniform density; full time dependent diffusion and Navier-Stokes equations are then used with constant eddy diffusion and viscosity coefficients, together with a constant Coriolis parameter. Scaling analysis reveals three independent scales of the problem including the radius of deformation of the inertial length, buoyancy length, and diffusive length scales. The governing equations are then suitably scaled and the resulting normalized equations are shown to depend on the Ekman number alone for problems of oceanic interest. It is concluded that the mean Gulf Stream dynamics can be interpreted in terms of a solution of the Navier-Stokes and diffusion equations, with the cross-stream circulation responsible for the maintenance of the front; this mechanism is suggested for the maintenance of the Gulf Stream dynamics.

Kao, T. W.↗

NASA NTRS1982

A numerical solution of the supersonic flow over a rearward facing step with transverse non-reacting hydrogen injection

The present investigation represents an application of computational fluid dynamics to a problem associated with the flow in the combustor region of a supersonic combustion ramjet engine (scramjet). The governing equations are considered, taking into account the Navier-Stokes equations, a molecular viscosity calculation, the molecular thermal conductivity, molecular diffusion, and a turbulence model. The employed numerical solution is patterned after the explicit, time-dependent, unsplit, predictor-corrector, finite-difference method given by MacCormack (1969). The calculation is concerned with the supersonic flow over a rearward-facing step with transverse H2 injection at conditions germane to the combustor region of a scramjet engine. The H2 jet acts as an effective body which essentially shields the primary flow from the rearward-facing step, thus substantially changing the wave pattern in the primary flow.

Berman, H. A.↗

NASA NTRS1983

Assessment of a 3-D boundary layer code to predict heat transfer and flow losses in a turbine

The prediction of the complete flow field in a turbine passage is an extremely difficult task due to the complex three dimensional pattern which contains separation and attachment lines, a saddle point and horseshoe vortex. Whereas, in principle such a problem can be solved using full Navier-Stokes equations, in reality methods based on a Navier-Stokes solution procedure encounter difficulty in accurately predicting surface quantities (e.g., heat transfer) due to grid limitations imposed by the speed and size of the existing computers. On the other hand the overall problem is strongly three dimensional and too complex to be analyzed by the current design methods based on inviscid and/or viscous strip theories. Thus there is a strong need for enhancing the current prediction techniques through inclusion of 3-D viscous effects. A potentially simple and cost effective way to achieve this is to use a prediction method based on three dimensional boundary layer (3-DBL) theory. The major objective of this program is to assess the applicability of such a 3-DBL approach for the prediction of heat loads, boundary layer growth, pressure losses and streamline skewing in critical areas of a turbine passage. A brief discussion of the physical problem addressed here along with the overall approach is presented.

Vatsa, V. N.↗

NASA NTRS1987

A technique for the prediction of airfoil flutter characteristics in separated flow

A solution procedure is described for determining the two-dimensional, one- or two-degree-of-freedom flutter characteristics of arbitrary airfoils at large angles of attack. The same procedure is used to predict stall flutter. This procedure requires a simultaneous integration in time of the solid and fluid equations of motion. The fluid equations of motion are the unsteady compressible Navier-Stokes equations, solved in a body-fitted moving coordinate system using an approximate factorization scheme. The solid equations of motion are integrated in time using an Euler implicit scheme. Flutter is said to occur if small disturbances imposed on the airfoil attitude lead to divergent oscillatory motions at subsequent times. Results for a number of special cases are presented to demonstrate the suitability of this scheme to predict flutter at large mean angles of attack. Some stall flutter applications are also presented.

Wu, Jiunn-Chi↗

NASA NTRS1988

A decade of computer simulations for space shuttle aerodynamics

Ten years ago computer simulations of the flow field around the Space Shuttle Orbiter were limited to inviscid calculations for the windward side of the forebody and viscous calculations for selected 2-D problems. Advances in computer hardware and numerical methods during the past ten years made it possible to calculate viscous flow over the complete orbiter configuration at angle of attack. The equations solved are the Reynolds-averaged, Navier-Stokes equation, simplified by either the thin-layer or parabolized approximation. An algebraic eddy viscosity model is used for turbulent flow. The free stream is assumed to be a perfect gas for wind tunnel conditions and a real gas in thermodynamic equilibrium for flight conditions. Four examples of recent computer simulations are presented. Flow field results include oil flow patterns on the surface and Mach number contours, isobars, and cross-flow velocity vectors in the shock layer.

Inouye, Mamoru↗

NASA NTRS1988

A decade of computer simulations for Space Shuttle aerodynamics

Inouye, Mamoru↗

NASA NTRS1989

Galerkin/least-squares procedures in computational fluid dynamics

The application of Galerkin/least-squares methodology to computational fluid dynamics is reviewed. Applications to advective-diffusive systems of equations and incompressible and compressible Euler and Navier-Stokes equations are described. Space-time finite-element methodology for time-dependent problems is emphasized and the status and role of so-called 'discontinuity-capturing operators' is reviewed.

Hughes, Thomas J. R.↗

NASA NTRS1991

Numerical investigation of an internal layer in turbulent flow over a curved hill

A numerical investigation of incompressible and compressible turbulent flows over strongly curved surfaces is presented. The turbulent flow equations are solved by a pressure based Navier-Stokes equations solver. In the method, the conservation of mass equation is replaced by a pressure correction equation applicable for both compressible and incompressible flows. The turbulence is described by a multiple time scale turbulence model supplemented with a near-wall turbulence model. The numerical results show that the internal layer is a strong turbulence field which is developed beneath the external boundary layer and is located very close to the wall. The development of the internal layer is attributed to the enormous mean flow strain rate caused by the streamline curvature. The external boundary layer flow responds rather slowly to the streamline curvature. Thus, the turbulence field of the forward corner of the curved hill is characterized by two turbulence fields interacting with each other. The turbulence intensity of the internal layer is much stronger than that of the external boundary layer, so that the development of a new boundary layer in the downstream region of the curved hill depends mostly on the internal layer. These numerical results are in good agreement with the measured data, and show that the turbulence model can resolve the turbulence field subjected to the strong streamline curvature.

Kim, S.-W.↗

NASA NTRS1990

Computation of instability and transition

The state of transition calculations in 1989 and desirable future trends are analyzed within the framework of aeronautical applications with emphasis on wall-bounded flows. A shift to studies of the early stages of transition, particularly to receptivity, from studies of the late stages in excessively confined periodic domains is proposed. It is pointed out that early-stage studies require a larger computational domain and more difficult boundary conditions but not a very fine grid. Intermediate strategies, more involved than the classical Orr-Sommerfeld equation but much more economical than full Navier-Stokes equations, are emphasized, along with the need for compressible calculations.

Spalart, P. R.↗