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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Random element method for numerical modeling of diffusional processes

The random element method is a generalization of the random vortex method that was developed for the numerical modeling of momentum transport processes as expressed in terms of the Navier-Stokes equations. The method is based on the concept that random walk, as exemplified by Brownian motion, is the stochastic manifestation of diffusional processes. The algorithm based on this method is grid-free and does not require the diffusion equation to be discritized over a mesh, it is thus devoid of numerical diffusion associated with finite difference methods. Moreover, the algorithm is self-adaptive in space and explicit in time, resulting in an improved numerical resolution of gradients as well as a simple and efficient computational procedure. The method is applied here to an assortment of problems of diffusion of momentum and energy in one-dimension as well as heat conduction in two-dimensions in order to assess its validity and accuracy. The numerical solutions obtained are found to be in good agreement with exact solution except for a statistical error introduced by using a finite number of elements, the error can be reduced by increasing the number of elements or by using ensemble averaging over a number of solutions.

Ghoniem, A. F.↗
NASA NTRS1987

Numerical experiments on boundary-layer receptivity

The incompressible laminar flow over an infinitely thin flat plate is obtained using a Navier-Stokes code in vorticity-velocity variables. The flow at and near the leading edge of the plate is an integral part of the solution algorithm which requires no special treatment; thus allowing for the flow field in this region to be studied in detail. An incident plane sound wave is imposed in the free-stream flow and the receptivity of the boundary layer is studied with particular emphasis to the flow near and at the leading edge.

Gatski, Thomas B.↗
NASA NTRS1989

Spectral solution of the incompressible Navier-Stokes equations on the Connection Machine 2

The issue of solving the time-dependent incompressible Navier-Stokes equations on the Connection Machine 2 is addressed, for the problem of transition to turbulence of the steady flow in a channel. The spectral algorithm used serially requires O(N(4)) operations when solving the equations on an N x N x N grid; using the massive parallelism of the CM, it becomes an O(N(2)) problem. Preliminary timings of the code, written in LISP, are included and compared with a corresponding code optimized for the Cray-2 for a 128 x 128 x 101 grid.

Tomboulian, Sherryl↗
NASA NTRS1990

Characteristic-based algorithms for flows in thermo-chemical nonequilibrium

A generalized finite-rate chemistry algorithm with Steger-Warming, Van Leer, and Roe characteristic-based flux splittings is presented in three-dimensional generalized coordinates for the Navier-Stokes equations. Attention is placed on convergence to steady-state solutions with fully coupled chemistry. Time integration schemes including explicit m-stage Runge-Kutta, implicit approximate-factorization, relaxation and LU decomposition are investigated and compared in terms of residual reduction per unit of CPU time. Practical issues such as code vectorization and memory usage on modern supercomputers are discussed.

Walters, Robert W.↗
NASA NTRS1989

Spectral solution of the incompressible Navier-Stokes equations on the Connection Machine2

The issue of solving the time-dependent incompressible Navier-Stokes equations on the Connection Machine 2 is addressed, for the problem of transition to turbulence of the steady flow in a channel. The spectral algorithm used serially requires O(N4) operations when solving the equations on an N x N x N grid; using the massive parallelism of the CM, it becomes an O(N2) problem. Preliminary timings of the code, written in LISP, are included and compared with a corresponding code optimized for the Cray-2 for a 128 x 128 x 101 grid.

Tomboulian, Sherryl↗
NASA NTRS1990

Numerical study of supersonic combustors by multi-block grids with mismatched interfaces

A three dimensional, finite rate chemistry, Navier-Stokes code was extended to a multi-block code with mismatched interface for practical calculations of supersonic combustors. To ensure global conservation, a conservative algorithm was used for the treatment of mismatched interfaces. The extended code was checked against one test case, i.e., a generic supersonic combustor with transverse fuel injection, examining solution accuracy, convergence, and local mass flux error. After testing, the code was used to simulate the chemically reacting flow fields in a scramjet combustor with parallel fuel injectors (unswept and swept ramps). Computational results were compared with experimental shadowgraph and pressure measurements. Fuel-air mixing characteristics of the unswept and swept ramps were compared and investigated.

Moon, Young J.↗
NASA NTRS1991

Unsteady shock-vortex interaction on a flexible delta wing

Unsteady Navier-Stokes computations have been carried out for simulating transonic flows over a clipped delta wing undergoing oscillatory and ramp motions, including flexibility. The implicit upwind algorithm has been validated by comparing the solutions with experimental data for the oscillatory pitching motion cases. The numerical and experimental results agree well at moderate angles of attack, where a leading-edge vortex develops. The ramp motion cases have demonstrated the effects of unsteadiness of the flow field and structural flexibility on the wing responses. For the 10 deg ramp motion, a vortex breakdown is observed. The interaction with the shock wave plays an essential role in the process of the breakdown observed in the present calculation.

Obayashi, Shigeru↗
NASA NTRS1980

On the numerical solution of time-dependent viscous incompressible fluid flows involving solid boundaries

An inherent numerical problem associated with the fully explicit pseudospectral numerical simulation of the incompressible Navier-Stokes equation for viscous flows with no-slip walls is described. A semi-implicit scheme which circumvents this numerical difficulty is presented. In this algorithm the equation of continuity rather than the Poisson equation for pressure is solved directly. Pseudospectral formulation of the channel flow problem using Fourier series and Chebyshev polynomials expansions is given for this scheme. An example demonstrating the applicability of the method is given.

Moin, P.↗
NASA NTRS1983

Heating analysis of bent-nose biconics at high angles of attack using the parabolized Navier-Stokes equations

A method based on the Parabolized Navier-Stokes equations is used to calculate the flow field and heat transfer of lifting entry vehicles. The method is based on the Bean and Warming implicit algorithm and uses a new procedure for preventing departure solutions. Calculations are carried out for blunt on-axis and bent biconics, assuming a perfect gas and laminar flow, and compared with available heat transfer, surface pressure and shock shape measurements for a range of Mach numbers and angles of attack. In all calculations presented here, the starting solution is obtained from available inviscid and boundary layer codes. Good agreement with experiment is indicated. Thus, the method provides an accurate and rather inexpensive procedure for calculating three-dimensional flows at supersonic Mach numbers.

Stephenson, B. L.↗
NASA NTRS1991

Extension of multigrid methodology to supersonic/hypersonic 3-D viscous flows

A multigrid acceleration technique developed for solving 3-D Navier-Stokes equations for subsonic/transonic flows was extended to supersonic/hypersonic flows. An explicit multistage Runge-Kutta type of time stepping scheme is used as the basic algorithm in conjunction with the multigrid scheme. Solutions were obtained for a blunt conical frustum at Mach 6 to demonstrate the applicability of the multigrid scheme to high speed flows. Computations were performed for a generic High Speed Civil Transport configuration designed to cruise at Mach 3. These solutions show both the efficiency and accuracy of the present scheme for computing high speed viscous flows over configurations of practical interest.

Vatsa, Veer N.↗
NASA NTRS1990

Numerical study of supersonic combustors by multi-block grids with mismatched interfaces

A three-dimensional, finite-rate chemistry, Navier-Stokes code has been extended to a multi-block with mismatched interface code for practical calculations of supersonic combustors. To ensure global conservation, a conservative algorithm was used for the treatment of mismatched interfaces. The extended code was checked against one test case, i.e., a generic supersonic combustor with transverse fuel injection, examining solution accuracy, convergence, and local mass flux error. After testing, the code was used to simulate the chemically reacting flow fields in a scramjet combustor with parallel fuel injectors (unswept and swept ramps). Computational results were compared with experimental shadowgraph and pressure measurements. Fuel-air mixing characteristics of the unswept and swept ramps were compared and investigated.

Moon, Young J.↗
NASA NTRS1993

Time-accurate Navier-Stokes calculations with multigrid acceleration

An efficient method for calculating unsteady flows is presented, with emphasis on a modified version of the thin-layer Navier-Stokes equations. Fourier stability analysis is used to illustrate the effect of treating the source term implicitly instead of explicity, as well as to illustrate other algorithmic choices. A 2D circular cylinder (with a Reynolds number of 1200 and a Mach number of 0.3) is calculated. The present scheme requires only about 10 percent of the computer time required by global minimum time stepping.

Melson, N. D.↗
NASA NTRS1995

Development of numerical methods for overset grids with applications for the integrated Space Shuttle vehicle

Algorithms and computer code developments were performed for the overset grid approach to solving computational fluid dynamics problems. The techniques developed are applicable to compressible Navier-Stokes flow for any general complex configurations. The computer codes developed were tested on different complex configurations with the Space Shuttle launch vehicle configuration as the primary test bed. General, efficient and user-friendly codes were produced for grid generation, flow solution and force and moment computation.

Chan, William M.↗
NASA NTRS1997

On Efficient Parallel Implementation of Moving Body Overset Grid Methods

An investigation into the parallel performance of moving-body overset grid methods will be presented. Parallel versions of the OVERFLOW flow solver, DCF3D domain connectivity software, and SIXDO six-degree-of-freedom routine are coupled with an automatic load balance routine and tested for 3D Navier-Stokes calculations on the IBM SP2. The primary source of parallel inefficiency in moving and problems are the domain connectivity costs with DCF 3D. Although this algorithm constitutes a relatively low fraction of the total solution cost (e.g. 10-20%) in calculations on serial machines, the consequently cause a significant degradation in the overall parallel performance. The paper will highlight some approaches for improving the scalability of DCF3D. The paper will present results of a proposed new load balancing scheme that seeks more equal distribution of the inter-grid boundary points in order to more evenly load balance the donor search costs associated with DCF3D. Some preliminary results will also be given from a new solution-adaption algorithm coupled with OVERFLOW which incorporates overset cartesian grids with various levels of refinement. The measured parallel performance from a descending delta-wing configuration and a generic store-separation from a wing/pylon case will be presented.

Wissink, Andrew M.↗
DOE OSTI2024

Asymptotic consistency of the WSINDy algorithm in the limit of continuum data

In this work we study the asymptotic consistency of the weak-form sparse identification of nonlinear dynamics algorithm (WSINDy) in the identification of differential equations from noisy samples of solutions. We prove that the WSINDy estimator is unconditionally asymptotically consistent for a wide class of models that includes the Navier–Stokes, Kuramoto–Sivashinsky and Sine–Gordon equations. We thus provide a mathematically rigorous explanation for the observed robustness to noise of weak-form equation learning. Conversely, we also show that, in general, the WSINDy estimator is only conditionally asymptotically consistent, yielding discovery of spurious terms with probability one if the noise level exceeds a critical threshold σ c . We provide explicit bounds on σ c in the case of Gaussian white noise and we explicitly characterize the spurious terms that arise in the case of trigonometric and/or polynomial libraries. Furthermore, we show that, if the data is suitably denoised (a simple moving average filter is sufficient), then asymptotic consistency is recovered for models with locally-Lipschitz, polynomial-growth nonlinearities. Our results reveal important aspects of weak-form equation learning, which may be used to improve future algorithms. We demonstrate our findings numerically using the Lorenz system, the cubic oscillator, a viscous Burgers-growth model and a Kuramoto–Sivashinsky-type high-order PDE.

asymptotic consistency↗
NASA NTRS1984

A compressible Navier-Stokes code for turbulent flow modeling

An implicit, finite volume code for solving two dimensional, compressible turbulent flows is described. Second order upwind differencing of the inviscid terms of the equations is used to enhance stability and accuracy. A diagonal form of the implicit algorithm is used to improve efficiency. Several zero and two equation turbulence models are incorporated to study their impact on overall flow modeling accuracy. Applications to external and internal flows are discussed.

Coakley, T. J.↗
NASA NTRS1993

Parallel execution of a three-dimensional, chemically reacting, Navier-Stokes code on distributed-memory machines

This paper describes the parallel version of the three-dimensional, chemically reacting, computational fluid dynamics (CFD) code, SPARK. This work was performed on the Intel iPSC/860-based parallel computers. The SPARK code utilizes relatively simple explicit numerical algorithms, but models complex chemical reactions. The code solves the equations over a regular structured mesh so a simple dam decomposition is used to assign work to the individual processors. The explicit nature of the algorithm, combined with the computational intensity of the chemistry calculations, results in a very low communication-to-computation ratio when compared to typical CFD codes. The efficiency of the parallel code is examined and shown to be about 65 percent when the problem size is scaled with the number of processors. Two low-angle wall-jet injection cases are solved to demonstrate the capability of the parallel code for solving large problems efficiently.

Otto, John C.↗
NASA NTRS1999

A Pseubo-Temporal Multi-Grid Relaxation Scheme for Solving the Parabolized Navier-Stokes Equations

A multi-grid, flux-difference-split, finite-volume code, VULCAN, is presented for solving the elliptic and parabolized form of the equations governing three-dimensional, turbulent, calorically perfect and non-equilibrium chemically reacting flows. The space marching algorithms developed to improve convergence rate and or reduce computational cost are emphasized. The algorithms presented are extensions to the class of implicit pseudo-time iterative, upwind space-marching schemes. A full approximate storage, full multi-grid scheme is also described which is used to accelerate the convergence of a Gauss-Seidel relaxation method. The multi-grid algorithm is shown to significantly improve convergence on high aspect ratio grids.

Morrison, J. H.↗