Search NASA⌕ Search

SEARCH · Search NASA

Results for “Navier Stokes algorithm”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Vortex shedding from solid rocket propellant inhibitors

Vortex shedding frequency caused by the protrusion of inhibitors into the flow field of a solid rocket motor is investigated by experimental and mathematical models. The time dependent Navier-Stokes equations are solved using a finite difference technique assuming incompressible, two-dimensional flow under both laminar and turbulent flow conditions. For laminar flow, explicit solutions are obtained using a vorticity-transport equation in place of the Navier-Stokes equations. For turbulent flow, a two-equation (k-epsilon) model is used for turbulent modeling and the SIMPLE algorithm is employed as the computational scheme. Cold flow tests were conducted to confirm the basic flow structure and to determine the vortex shedding frequency under both laminar and turbulent flow conditions. The vortex shedding frequencies were determined using a stroboscope to measure the oscillating frequency of yarn tufts which were fastened to one inhibitor in the models. A hot-film anemometer established the velocity history behind the inhibitor. Good agreement between the theoretical results and measurements of the vortex shedding frequencies is demonstrated.

Shu, P. H.↗
NASA NTRS1986

Solution of separated flows using a simultaneous iteration technique

Generalization and improvement of an earlier work developed for studying separated flows using boundary layer type equations are presented. The improvements include extensions to a general coordinate system and use of a more general zonal technique for solving the coupled equations. In order to be able to consider arbitrary geometries, second order accurate (in space) conservative differences are generated by considering the integral formulation of the governing equations in a general coordinate system. The general coordinate system is handled in as general a manner as possible to allow for the use of either analytically or numerically generated coordinate systems. A marching procedure was used for solving the Partially Parabolized Navier-Stokes (PPNS) equations in the viscous region coupled in a fully implicit manner with the elliptic inviscid equation. To test the algorithm and compare to other solutions, solutions for flow over a flat plate and flow past the symmetrical 12 percent thick Joukowski airfoil (J012) at zero angle of attack were presented.

Source record↗
NASA NTRS1984

Feasibility of using the Massively Parallel Processor for large eddy simulations and other Computational Fluid Dynamics applications

The results of an investigation into the feasibility of using the MPP for direct and large eddy simulations of the Navier-Stokes equations is presented. A major part of this study was devoted to the implementation of two of the standard numerical algorithms for CFD. These implementations were not run on the Massively Parallel Processor (MPP) since the machine delivered to NASA Goddard does not have sufficient capacity. Instead, a detailed implementation plan was designed and from these were derived estimates of the time and space requirements of the algorithms on a suitably configured MPP. In addition, other issues related to the practical implementation of these algorithms on an MPP-like architecture were considered; namely, adaptive grid generation, zonal boundary conditions, the table lookup problem, and the software interface. Performance estimates show that the architectural components of the MPP, the Staging Memory and the Array Unit, appear to be well suited to the numerical algorithms of CFD. This combined with the prospect of building a faster and larger MMP-like machine holds the promise of achieving sustained gigaflop rates that are required for the numerical simulations in CFD.

Bruno, John↗
NASA NTRS1990

Preconditioned conjugate gradient methods for the compressible Navier-Stokes equations

The compressible Navier-Stokes equations are solved for a variety of two-dimensional inviscid and viscous problems by preconditioned conjugate gradient-like algorithms. Roe's flux difference splitting technique is used to discretize the inviscid fluxes. The viscous terms are discretized by using central differences. An algebraic turbulence model is also incorporated. The system of linear equations which arises out of the linearization of a fully implicit scheme is solved iteratively by the well known methods of GMRES (Generalized Minimum Residual technique) and Chebyschev iteration. Incomplete LU factorization and block diagonal factorization are used as preconditioners. The resulting algorithm is competitive with the best current schemes, but has wide applications in parallel computing and unstructured mesh computations.

Venkatakrishnan, V.↗
NASA NTRS1990

Shock interference studies on a circular cylinder at Mach 16

This paper presents the results of a computational study on the shock interference problems on a cylindrical body typical of an engine inlet cowl leading edge at a nominal Mach number of 16. The two-dimensional Navier-Stokes equations are solved assuming the flow to be in chemical and thermal equilibrium and using a finite element method. The algorithm employs a cell-centered fully implicit upwind scheme. Adaptively generated unstructured meshes of triangles and quadrilaterals are employed. Under certain conditions the finite element code resulted in oscillatory solutions for shock interference at Mach 16. Some of the causes of the unsteady behavior are identified, and to the extent possible, such situations are avoided in the present application. Two shock interference conditions involving a Type IV (supersonic jet) and a Type III (attaching shear layer) are solved. The results are compared with available experimental data and reasonable agreement is observed. A semi-empirical method is also used to estimate the maximum surface heat flux and static pressure.

Prabhu, Ramadas K.↗
NASA NTRS1992

Automatic differentiation of advanced CFD codes for multidisciplinary design

Automated multidisciplinary design of aircraft and other flight vehicles requires the optimization of complex performance objectives with respect to a number of design parameters and constraints. The effect of these independent design variables on the system performance criteria can be quantified in terms of sensitivity derivatives which must be calculated and propagated by the individual discipline simulation codes. Typical advanced CFD analysis codes do not provide such derivatives as part of a flow solution; these derivatives are very expensive to obtain by divided (finite) differences from perturbed solutions. It is shown that sensitivity derivatives can be obtained accurately and efficiently using the ADIFOR source translator for automatic differentiation. In particular, it is demonstrated that the 3-D, thin-layer Navier-Stokes, multigrid flow solver called TLNS3D is amenable to automatic differentiation in the forward mode even with its implicit iterative solution algorithm and complex turbulence modeling. It is significant that by using computational differentiation, consistent discrete nongeometric sensitivity derivatives have been obtained from an aerodynamic 3-D CFD code in a relatively short time, e.g., O(man-week) not O(man-year).

Bischof, C.↗
NASA NTRS1994

Multi-dimensional combustor flowfield analyses in gas-gas rocket engine

The objectives of the present research are to improve design capabilities for low thrust rocket engines through understanding of the detailed mixing and combustions processes. Of particular interest is a small gaseous hydrogen-oxygen thruster which is considered as a coordinated part of an on-going experimental program at NASA LeRC. Detailed computational modeling requires the application of the full three-dimensional Navier Stokes equations, coupled with species diffusion equations. The numerical procedure is performed on both time-marching and time-accurate algorithms and using an LU approximate factorization in time, flux split upwinding differencing in space. The emphasis in this paper is focused on using numerical analysis to understand detailed combustor flowfields, including the shear layer dynamics created between fuel film cooling and the core gas in the vicinity on the nearby combustor wall; the integrity and effectiveness of the coolant film; three-dimensional fuel jets injection/mixing/combustion characteristics; and their impacts on global engine performance.

Tsuei, Hsin-Hua↗
NASA NTRS1994

Optimal control of thermally coupled Navier Stokes equations

The optimal boundary temperature control of the stationary thermally coupled incompressible Navier-Stokes equation is considered. Well-posedness and existence of the optimal control and a necessary optimality condition are obtained. Optimization algorithms based on the augmented Lagrangian method with second order update are discussed. A test example motivated by control of transport process in the high pressure vapor transport (HVPT) reactor is presented to demonstrate the applicability of our theoretical results and proposed algorithm.

Ito, Kazufumi↗
NASA NTRS1995

A fast Poisson solver for unsteady incompressible Navier-Stokes equations on the half-staggered grid

In this paper, a fast Poisson solver for unsteady, incompressible Navier-Stokes equations with finite difference methods on the non-uniform, half-staggered grid is presented. To achieve this, new algorithms for diagonalizing a semi-definite pair are developed. Our fast solver can also be extended to the three dimensional case. The motivation and related issues in using this second kind of staggered grid are also discussed. Numerical testing has indicated the effectiveness of this algorithm.

Golub, G. H.↗
DOE OSTI2026

A fourth order sharp immersed method for the incompressible Navier-Stokes equations with stationary and moving boundaries and interfaces

We propose a fourth order Navier-Stokes solver based on the immersed interface method (IIM), for flow problems with stationary and one-way coupled moving boundaries and interfaces. Our algorithm employs a Runge-Kutta-based projection method that maintains high-order temporal accuracy in both velocity and pressure for steady and unsteady velocity boundary conditions. Fourth order spatial accuracy is achieved through a novel fifth order IIM discretization scheme for the advection term, as well as existing high-order interface-corrected finite difference schemes for the other differential operators. Using a set of manufactured flow problems with stationary and moving boundaries, we demonstrate fourth order convergence of velocity and pressure in the infinity norm, both inside the domain and on the immersed boundaries. The solver’s performance is further validated through a range of practical flow simulations, highlighting its efficiency over a second order scheme. Finally, we showcase the ability of our immersed discretization scheme to handle interface-coupled multiphysics problems by solving a conjugate heat transfer problem with multiple immersed solids. Overall, the proposed approach robustly combines the efficiency of high order discretization schemes with the flexibility of immersed discretizations for flow problems with complex, moving boundaries and interfaces.

42 ENGINEERING↗
NASA NTRS1993

Advanced adaptive computational methods for Navier-Stokes simulations in rotorcraft aerodynamics

A phase 2 research and development effort was conducted in area transonic, compressible, inviscid flows with an ultimate goal of numerically modeling complex flows inherent in advanced helicopter blade designs. The algorithms and methodologies therefore are classified as adaptive methods, which are error estimation techniques for approximating the local numerical error, and automatically refine or unrefine the mesh so as to deliver a given level of accuracy. The result is a scheme which attempts to produce the best possible results with the least number of grid points, degrees of freedom, and operations. These types of schemes automatically locate and resolve shocks, shear layers, and other flow details to an accuracy level specified by the user of the code. The phase 1 work involved a feasibility study of h-adaptive methods for steady viscous flows, with emphasis on accurate simulation of vortex initiation, migration, and interaction. Phase 2 effort focused on extending these algorithms and methodologies to a three-dimensional topology.

Stowers, S. T.↗
NASA NTRS2004

Automated Euler and Navier-Stokes Database Generation for a Glide-Back Booster

The past two decades have seen a sustained increase in the use of high fidelity Computational Fluid Dynamics (CFD) in basic research, aircraft design, and the analysis of post-design issues. As the fidelity of a CFD method increases, the number of cases that can be readily and affordably computed greatly diminishes. However, computer speeds now exceed 2 GHz, hundreds of processors are currently available and more affordable, and advances in parallel CFD algorithms scale more readily with large numbers of processors. All of these factors make it feasible to compute thousands of high fidelity cases. However, there still remains the overwhelming task of monitoring the solution process. This paper presents an approach to automate the CFD solution process. A new software tool, AeroDB, is used to compute thousands of Euler and Navier-Stokes solutions for a 2nd generation glide-back booster in one week. The solution process exploits a common job-submission grid environment, the NASA Information Power Grid (IPG), using 13 computers located at 4 different geographical sites. Process automation and web-based access to a MySql database greatly reduces the user workload, removing much of the tedium and tendency for user input errors. The AeroDB framework is shown. The user submits/deletes jobs, monitors AeroDB's progress, and retrieves data and plots via a web portal. Once a job is in the database, a job launcher uses an IPG resource broker to decide which computers are best suited to run the job. Job/code requirements, the number of CPUs free on a remote system, and queue lengths are some of the parameters the broker takes into account. The Globus software provides secure services for user authentication, remote shell execution, and secure file transfers over an open network. AeroDB automatically decides when a job is completed. Currently, the Cart3D unstructured flow solver is used for the Euler equations, and the Overflow structured overset flow solver is used for the Navier-Stokes equations. Other codes can be readily included into the AeroDB framework.

Chaderjian, Neal M.↗
DOE OSTI2024

Shape effects on the local dynamics of suspensions of spheroidal particles

The effect of shape on the dynamics of suspensions of non-spherical heavy particles is examined by fully resolved numerical simulations of oblate and prolate spheroids, as well as spheres, for a density ratio of ten, volume fractions ranging from 0.5% to 5%, and Reynolds numbers between 20 and 30. The dynamics is determined both by the interactions of the particles with the fluid as well as by collisions, with the number and importance of collisions increasing with volume fractions. A single isolated oblate or prolate spheroid falling under gravity generally falls broadside on, for the governing parameters examined here, and at low-volume fractions, the majority of particles in a suspension fall that way. At higher-volume fractions, the orientation is more random. The slip velocity decreases as the volume fraction increases for all shapes, as expected, but the effect of the shape is much less than seen for a single particle. Furthermore, this seems to be due to two effects. For all volume fractions, the most deformed particles cluster more than spheres and less deformed particles, which increases their slip velocity. As the concentration increases, the increased particle interactions also causes more particles to fall short side-on, which reduces the frontal area and the resulting drag, increasing the slip velocity. This second effect is, of course, absent for spherical particles.

42 ENGINEERING↗
NASA NTRS1978

Finite element analysis of turbulent flows

The purpose of the reported investigation is to provide a rather comprehensive statement of the problem of turbulent flow prediction, and to present and place in context representative finite element numerical solutions obtained to date. A general statement of the solution algorithm is presented, and definitive comparisons between finite element-derived and conventional and/or popular finite difference algorithms is made in the few places where available. In a consideration of the theoretical formulation attention is given to the Navier-Stokes equations, Reynolds stress equations, and parabolic Navier-Stokes equations.

Baker, A. J.↗
NASA NTRS1981

Direct numerical simulations of turbulent shear flows

Numerical simulations of wakes of axisymmetric bodies and of turbulent mixing layers are reported. The flows were assumed to be statistically homogeneous in the mean flow direction, in concert with experimental data and the self-similarity theorem. The nonlinear Navier-Stokes equations were solved by a pseudo-spectral numerical method using a 32 x 32 x 33 point grid and an algorithm for fast Fourier transforms and inverse transforms. Leapfrog time differencing was employed on nonlinear terms and time differencing on viscous terms. Towed wakes and wakes behind a self-propelled body were simulated, showing that the towed wakes exhibited a proper temporal behavior after an initial period of adjustment, including the development of a kurtosis near the wake edge, which is experimentally verifiable. The mixing-layer simulation displayed the laboratory demonstrated presence of large scale features such as vortex cores, while the lateral coherence was weak.

Metcalfe, R. W.↗
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 NTRS1989

Numerical Simulation of Flow Through an Artificial Heart

A solution procedure was developed that solves the unsteady, incompressible Navier-Stokes equations, and was used to numerically simulate viscous incompressible flow through a model of the Pennsylvania State artificial heart. The solution algorithm is based on the artificial compressibility method, and uses flux-difference splitting to upwind the convective terms; a line-relaxation scheme is used to solve the equations. The time-accuracy of the method is obtained by iteratively solving the equations at each physical time step. The artificial heart geometry involves a piston-type action with a moving solid wall. A single H-grid is fit inside the heart chamber. The grid is continuously compressed and expanded with a constant number of grid points to accommodate the moving piston. The computational domain ends at the valve openings where nonreflective boundary conditions based on the method of characteristics are applied. Although a number of simplifing assumptions were made regarding the geometry, the computational results agreed reasonably well with an experimental picture. The computer time requirements for this flow simulation, however, are quite extensive. Computational study of this type of geometry would benefit greatly from improvements in computer hardware speed and algorithm efficiency enhancements.

Rogers, Stuart E.↗
NASA NTRS1989

An investigation of cell centered and cell vertex multigrid schemes for the Navier-Stokes equations

Two efficient and robust finite-volume multigrid schemes for solving the Navier-Stokes equations are investigated. These schemes employ either a cell centered or a cell vertex discretization technique. An explicit Runge-Kutta algorithm is used to advance the solution in time. Acceleration techniques are applied to obtain faster steady-state convergence. Accuracy and convergence of the schemes are examined. Computational results for transonic airfoil flows are essentially the same, even for a coarse mesh. Both schemes exhibit good convergence rates for a broad range of artificial dissipation coefficients.

Radespiel, R.↗