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 991 records · Page 55

All sources NASA NTRS DOE OSTI

NASA NTRS1987

A low dimensional dynamical system for the wall layer

Low dimensional dynamical systems which model a fully developed turbulent wall layer were derived.The model is based on the optimally fast convergent proper orthogonal decomposition, or Karhunen-Loeve expansion. This decomposition provides a set of eigenfunctions which are derived from the autocorrelation tensor at zero time lag. Via Galerkin projection, low dimensional sets of ordinary differential equations in time, for the coefficients of the expansion, were derived from the Navier-Stokes equations. The energy loss to the unresolved modes was modeled by an eddy viscosity representation, analogous to Heisenberg's spectral model. A set of eigenfunctions and eigenvalues were obtained from direct numerical simulation of a plane channel at a Reynolds number of 6600, based on the mean centerline velocity and the channel width flow and compared with previous work done by Herzog. Using the new eigenvalues and eigenfunctions, a new ten dimensional set of ordinary differential equations were derived using five non-zero cross-stream Fourier modes with a periodic length of 377 wall units. The dynamical system was integrated for a range of the eddy viscosity prameter alpha. This work is encouraging.

Aubry, N.↗

NASA NTRS1993

Triangular spectral elements for incompressible fluid flow

We discuss the use of triangular elements in the spectral element method for direct simulation of incompressible flow. Triangles provide much greater geometric flexibility and are better conditioned and more accurate than quadrilateral elements when small angles arise. We employ a family of tensor product algorithms for triangles, allowing triangular elements to be handled with comparable arithmetic complexity to quadrilateral elements. The triangular discretizations are applied to the Poisson equation and are validated. The triangular discretizations are then applied to the incompressible Navier-Stokes equations, and a laminar channel flow solution is given. The new triangular spectral elements can be combined with standard quadrilateral elements, yielding a general and flexible high order method for complex geometries in two dimensions.

Mavriplis, Catherine↗

NASA NTRS1992

Methodology for sensitivity analysis, approximate analysis, and design optimization in CFD for multidisciplinary applications

Fundamental equations of aerodynamic sensitivity analysis and approximate analysis for the two dimensional thin layer Navier-Stokes equations are reviewed, and special boundary condition considerations necessary to apply these equations to isolated lifting airfoils on 'C' and 'O' meshes are discussed in detail. An efficient strategy which is based on the finite element method and an elastic membrane representation of the computational domain is successfully tested, which circumvents the costly 'brute force' method of obtaining grid sensitivity derivatives, and is also useful in mesh regeneration. The issue of turbulence modeling is addressed in a preliminary study. Aerodynamic shape sensitivity derivatives are efficiently calculated, and their accuracy is validated on two viscous test problems, including: (1) internal flow through a double throat nozzle, and (2) external flow over a NACA 4-digit airfoil. An automated aerodynamic design optimization strategy is outlined which includes the use of a design optimization program, an aerodynamic flow analysis code, an aerodynamic sensitivity and approximate analysis code, and a mesh regeneration and grid sensitivity analysis code. Application of the optimization methodology to the two test problems in each case resulted in a new design having a significantly improved performance in the aerodynamic response of interest.

Taylor, Arthur C., III↗

NASA NTRS1993

The dependence of the dissipation on rotation

A form of the evolution equation for the dissipation showing the effects of rotation is derived directly from the Navier-Stokes equations. Rotation affects the dissipation if the flow is inhomogeneous or compressible. This is due to the fact that pressure scales with inverse Rossby number and certain pressure correlations usually neglected have to be retained when Re Ro is approximately equal to 1. The term reflecting the effects of rotation on the dissipation is a volume integral of the two-point velocity correlation similar in structure to the 'rapid-pressure' correlation appearing in the second-moment equations. Some of its properties are indicated.

Ristorcelli, J. R.↗

NASA NTRS2002

The Delta x B = 0 Constraint Versus Minimization of Numerical Errors in MHD Simulations

The MHD equations are a system of non-strictly hyperbolic conservation laws. The non-convexity of the inviscid flux vector resulted in corresponding Jacobian matrices with undesirable properties. It has previously been shown by Powell et al. (1995) that an 'almost' equivalent MHD system in non-conservative form can be derived. This non-conservative system has a better conditioned eigensystem. Aside from Powell et al., the MHD equations can be derived from basic principles in either conservative or non-conservative form. The Delta x B = 0 constraint of the MHD equations is only an initial condition constraint, it is very different from the incompressible Navier-Stokes equations in which the divergence condition is needed to close the system (i.e., to have the same number of equations and the same number of unknown). In the MHD formulations, if Delta x B = 0 initially, all one needs is to construct appropriate numerical schemes that preserve this constraint at later time evolutions. In other words, one does not need the Delta x B condition to close the MHD system. We formulate our new scheme together with the Cargo & Gallice (1997) form of the MHD approximate Riemann solver in curvilinear grids for both versions of the MHD equations. A novel feature of our new method is that the well-conditioned eigen-decomposition of the non-conservative MHD equations is used to solve the conservative equations. This new feature of the method provides well-conditioned eigenvectors for the conservative formulation, so that correct wave speeds for discontinuities are assured. The justification for using the non-conservative eigen-decomposition to solve the conservative equations is that our scheme has a better control of the numerical error associated with the divergence of the magnetic condition. Consequently, computing both forms of the equations with the same eigen-decomposition is almost equivalent. It will be shown that this approach, using the non-conservative eigensystem when solving the conservative equations, also works well in the context of standard shock-capturing schemes.

Yee, H. C.↗

NASA NTRS2009

Wind-US User's Guide, Version 2.0

Wind-US is a computational platform which may be used to numerically solve various sets of equations governing physical phenomena. Currently, the code supports the solution of the Euler and Navier-Stokes equations of fluid mechanics, along with supporting equation sets governing turbulent and chemically reacting flows. Wind-US is a product of the NPARC Alliance, a partnership between the NASA Glenn Research Center (GRC) and the Arnold Engineering Development Center (AEDC) dedicated to the establishment of a national, applications-oriented flow simulation capability. The Boeing Company has also been closely associated with the Alliance since its inception, and represents the interests of the NPARC User's Association. The "Wind-US User's Guide" describes the operation and use of Wind-US, including: a basic tutorial; the physical and numerical models that are used; the boundary conditions; monitoring convergence; the files that are read and/or written; parallel execution; and a complete list of input keywords and test options.

Towne, Charles E.↗

DOE OSTI2022

PsDNS

PsDNS is a Python package which makes it easy to solve partial differential equation using a pseudo-spectral method. The code uses Python classes and Numpy-like arrays to provide a simple interface with which users can construct equations, integrators, and diagnostics. It also includes implementation for common problems, in particular, the incompressible Navier-Stokes equation in a periodic domain, which is commonly used as a benchmark problem for turbulence research. PsDNS uses MPI to allow massively parallel computation for large problems.

Israel, Daniel↗

NASA NTRS1980

Prediction and measurement of turbulent aerodynamic trailing edge flows

A viscous-inviscid interaction algorithm is developed for prediction of two-dimensional mean and fluctuating velocity distributions in the wake immediately downstream of an airfoil trailing edge. A composite pressure field is defined, and a Poisson equation solved for transverse pressure variations. A parabolized form of the time-averaged steady Navier-Stokes equations are solved in conjunction with a viscous-augmented two-dimensional inviscid potential flow analysis. A tensor constitutive equation is employed to predict Reynolds stress distributions from solutions of a turbulence kinetic energy two equation closure model. Numerical predictions compared favorably with detailed experimental data for mean and fluctuating velocities, and Reynolds shear stress distributions, in the trailing edge region of a NACA 63-012 airfoil.

Baker, A. J.↗

NASA NTRS1980

A simple numerical orthogonal coordinate generator for fluid dynamic applications

An application of a simple numerical technique which allows for the rapid construction of orthogonal coordinate systems about two dimensional and axisymmetric bodies is presented. The technique is based on a predictor corrector numerical method. It can be used to generate orthogonal meshes which have unequally spaced points in two directions. These orthogonal meshes in their transformed computational plane are, however, equally spaced so that the differencing for the metric coefficients and the fluid dynamic equation terms can be easily determined using equally spaced central finite differences. Solutions to the Navier-Stokes equations for flow over blunt bodies with reverse curvature are presented. The coupling of the time dependent fluid dynamic equations and the coordinate generator worked well with no undesirable effects noted.

Graves, R. A., Jr.↗

NASA NTRS1987

A theoretical investigation of laser-sustained plasma thruster

A numerical code has been successfully developed for the investigation of thruster performance using a laser-sustained hydrogen plasma as the propellant. The plasma was sustained using a 10.6-micron CO2 laser beam focused at different positions within the thruster. The physical model assumed that plasma is in thermodynamical equilibrium (LTE), and geometric ray tracing was adopted to describe the laser beam. The steady-state, axisymmetric, Navier-Stokes equations coupled with the laser power absorption process have been solved numerically. A pressure based Navier-Stokes numerical solver using body-fitted coordinates was used to calculate the laser-supported rocket flow which includes both subsonic and supersonic flow regions. From the limited parametric study, which did not try to optimize the rocket performance, it was found that better performance was obtained when the laser beam was focused closer to the rocket throat.

Jeng, San-Mou↗

NASA NTRS1987

A theoretical evaluation of laser-sustained plasma thruster performance

An extensive numerical experiment has been conducted to evaluate rocket thruster performance using a laser-sustained hydrogen plasma as the propellant. The plasma was sustained using a 30 kW CO2 laser beam operated at 10.6 microns focused inside the thruster. The steady-state Navier-Stokes equations coupled with the laser power absorption process have been solved numerically. A pressure based Navier-Stokes solver using body-fitted coordinate was used to calculate the laser-supported rocket flow which included both recirculating and transonic flow regions. The local thermodynamic equilibrium (LTE) assumption was used for the plasma thermophysical and optical properties. Geometric ray tracing was adopted to describe the laser beam. Several different throat size thrusters operated at 150 and 300 kPa chamber stagnation pressure were studied. It was found that the thruster performance (vacuum specific impulse) was highly dependent on the operating conditions, and a properly designed laser supported thruster can attain a specific impulse around 1500 secs. The heat loading on the thruster wall was also estimated and was in the range of that for a conventional chemical rocket.

Jeng, San-Mou↗

NASA NTRS1987

Navier-Stokes cascade analysis with a stiff Kappa-Epsilon turbulence solver

The two dimensional, compressible, thin layer Navier-Stokes equations with the Baldwin-Lomax turbulence model and the kinetic energy-energy dissipation (k-epsilon) model are solved numerically to simulate the flow through a cascade. The governing equations are solved for the entire flow domain, without the boundary layer assumptions. The stiffness of the k-epsilon equations is discussed. A semi-implicit, Runge-Kutta, time-marching scheme is developed to solve the k-epsilon equations. The impact of the k-epsilon solver on the explicit Runge-Kutta Navier-Stokes solver is discussed. Numerical solutions are presented for two dimensional turbulent flow over a flat plate and a double circular arc cascade and compared with experimental data.

Liu, Jong-Shang↗

NASA NTRS1988

Navier-Stokes cascade analysis with a stiff k-epsilon turbulence solver

Liu, Jong-Shang↗

NASA NTRS1994

Low-dimensional description of the dynamics in separated flow past thick airfoils

Results are presented for the numerical simulation of unsteady viscous incompressible flow past thick airfoils. Specifically, flow past a NACA 4424 at an angle of attack of 2.5 deg and Reynolds numbers in the range of 1700-4000 has been simulated using the spectral element method. At these conditions the flow is separatedd and an unsteady wake is formed. Application of the method of empirical eigenfunction reveals the structure of the most energetic components of the flow. These are found to occur in pairs that, through phase exchange, are responsible for the vortex shedding. A set of ordinary differential equations is obtained for the amplitudes of these eigenfunctions by a Galerkin projection of the Navier-Stokes equations. The solutions of the model system are compared with the full simulation. The work is of relevance to the transition process and observed routes to chaos in airfoil wakes.

Deane, Anil E.↗

NASA NTRS1996

Model equations for simulating flows in multistage turbomachinery

A steady, three dimensional average-passage equation system was derived. The purpose was to simulate multistage turbomachinery flows. These equations describe a steady, viscous flow that is periodic from blade passage to blade passage. Moreover, these equations have a closure problem that is similar to that of the Reynolds-average Navier-Stokes equations. A scaled form of the average-passage equation system could provide an improved mathematical model for simulating the flow in the design and in the off-design conditions of a multistage machine.

Adamczyk, John J.↗

NASA NTRS2002

An Evaluation of Linear Instability Waves as Sources of Sound in a Supersonic Turbulent Jet

Mach wave radiation from supersonic jets is revisited to better quantify the extent to which linearized equations represent the details of the actual mechanism. To this end, we solve the linearized Navier-Stokes equations (LNS) with precisely the same mean flow and inflow disturbances as a previous direct numerical simulation (DNS) of a perfectly expanded turbulent M = 1.92 jet. We restrict our attention to the first two azimuthal modes, n = 0 and n = 1, which constitute most of the acoustic field. The direction of peak radiation and the peak Strouhal number matches the DNS reasonably well, which is in accord with previous experimental justification of the linear theory. However, it is found that the sound pressure level predicted by LNS is significantly lower than that from DNS. In order to investigate the discrepancy, individual frequency components of the solution are examined. These confirm that near the peak Strouhal number, particularly for the first helical mode n = 1, the amplification of disturbances in the LNS closely matches the DNS. However, away from the peak frequency (and generally for the azimuthal mode n = 0), modes in the LNS are damped while those in the DNS grow at rates comparable to those at the peak Strouhal number.

Mohseni, Kamran↗

NASA NTRS2004

An Investigation of Wave Propagations in Discontinuous Galerkin Method

Analysis of the discontinuous Galerkin method has been carried out for one- and two-dimensional system of hyperbolic equations. Analytical, as well as numerical, properties of wave propagation in a DGM scheme are derived and verified with direct numerical simulations. In addition to a systematic examination of the dissipation and dispersion errors, behaviours of a DG scheme at an interface of two different grid topologies are also studied. Under the same framework, a quantitative discrete analysis of various artificial boundary conditions is also conducted. Progress has been made in numerical boundary condition treatment that is closely related to the application of DGM in aeroacoustics problems. Finally, Fourier analysis of DGM for the Convective diffusion equation has also be studied in connection with the application of DG schemes for the Navier-Stokes equations. This research has resulted in five(5) publications, plus one additional manuscript in preparation, four(4) conference presentations, and three(3) departmental seminars, as summarized in part II. Abstracts of papers are given in part 111 of this report.

Hu, Fang Q.↗

NASA NTRS2019

Streamlined Convergence Acceleration for CFD Codes

Enigma, a simplified interface to the PETSc library, is shown to enable the rapid solution of discrete partial differential equations. Two CFD codes, LAURA and HyperSolve, use Enigma to compute steady solutions of the Navier-Stokes equations. Using PETSc, Enigma is shown to provide a Jacobian-Free Newton-Krylov method (JFNK), globalized with pseudotransient continuation, that improves efficiency over the point-implicit relaxation method traditionally used by LAURA. It is shown that iterative error has a large impact on surface heat transfer predicted by LAURA on an axisymmetric sphere-cone geometry. Also, the convergence rate of HyperSolve simulating subsonic flow over a delta wing geometry with the JFNK method is shown to be more efficient than employing a defect correction method as the nonlinear solver.

Thompson, Kyle B.↗