Search NASA⌕ Search

SEARCH · Search NASA

Results for “time-dependent Navier–Stokes”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Flux-vector splitting algorithm for chain-rule conservation-law form

A flux-vector splitting algorithm with Newton-Raphson iteration was developed for the 'full compressible' Navier-Stokes equations cast in chain-rule conservation-law form. The algorithm is intended for problems with deforming spatial domains and for problems whose governing equations cannot be cast in strong conservation-law form. The usefulness of the algorithm for such problems was demonstrated by applying it to analyze the unsteady, two- and three-dimensional flows inside one combustion chamber of a Wankel engine under nonfiring conditions. Solutions were obtained to examine the algorithm in terms of conservation error, robustness, and ability to handle complex flows on time-dependent grid systems.

Shih, T. I.-P.↗
NASA NTRS1993

Computational fluid dynamic modelling of cavitation

Models in sheet cavitation in cryogenic fluids are developed for use in Euler and Navier-Stokes codes. The models are based upon earlier potential-flow models but enable the cavity inception point, length, and shape to be determined as part of the computation. In the present paper, numerical solutions are compared with experimental measurements for both pressure distribution and cavity length. Comparisons between models are also presented. The CFD model provides a relatively simple modification to an existing code to enable cavitation performance predictions to be included. The analysis also has the added ability of incorporating thermodynamic effects of cryogenic fluids into the analysis. Extensions of the current two-dimensional steady state analysis to three-dimensions and/or time-dependent flows are, in principle, straightforward although geometrical issues become more complicated. Linearized models, however offer promise of providing effective cavitation modeling in three-dimensions. This analysis presents good potential for improved understanding of many phenomena associated with cavity flows.

Deshpande, Manish↗
NASA NTRS1994

An Experimental Study of the Flowfield on a Semispan Rectangular Wing with a Simulated Glaze Ice Accretion

Wind tunnel experiments were conducted in order to study the effect of a simulated glaze ice accretion on the flowfield of a semispan, reflection-plane, rectangular wing at Re = 1.5 million and M = 0.12. A laser Doppler velocimeter was used to map the flowfield on the upper surface of the model in both the clean and iced configurations at alpha = 0, 4, and 8 degrees angle of attack. At low angles of attack, the massive separation bubble aft of the leading edge ice horn was found to behave in a manner similar to laminar separation bubbles. At alpha = 0 and 4 degrees, the locations of transition and reattachment, as deduced from momentum thickness distributions, were found to be in good agreement with transition and reattachment locations in laminar separation bubbles. These values at y/b = 0.470, the centerline measurement location, matched well with data obtained on a similar but two dimensional model. The measured velocity profiles on the iced wing compared reasonably with the predicted profiles from Navier-Stokes computations. The iced-induced separation bubble was also found to have features similar to the recirculating region aft of rearward-facing steps. At alpha = 0 degrees and 4 degrees, reverse flow magnitudes and turbulence intensity levels were typical of those found in the recirculating region aft of rearward-facing steps. The calculated separation streamline aft of the ice horn at alpha = 4 degrees, y/b = 0.470 coincided with the locus of the maximum Reynolds normal stress. The maximum Reynolds normal stress peaked at two locations along the separation streamline. The location of the first peak-value coincided with the transition location, as deduced from the momentum thickness distributions. The location of the second peak was just upstream of reattachment, in good agreement with measurements of flows over similar obstacles. The intermittency factor in the vicinity of reattachment at alpha = 4 degrees, y/b = 0.470, revealed the time-dependent nature of the reattachment process. The size and extent of the separation bubble were found to be a function of angle of attack and the spanwise location. Three dimensional effects were found to be strongest at alpha = 8 degrees. The calculated separation and stagnation streamlines were found to vary little with spanwise location at alpha = 0 degrees. The calculated separation streamlines at alpha = 4 degrees revealed that the bubble was largest near the centerline measurement plane, whereas the tip-induced vortex flow and the model root-tunnel wall boundary-layer interaction reduced the size of the bubble. These effects were found to be most dramatic at alpha = 8 degrees.

Khodadoust, Abdollah↗
DOE OSTI2025

Analysis of the weighted shifted boundary method for the Poisson and Stokes problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. Accuracy is maintained by properly shifting the location and values of the boundary conditions. This avoids integration over cut cells and the associated implementation issues. Recently, the Weighted SBM (WSBM) was proposed for the Navier-Stokes equations with free surfaces and the Stokes flow with moving boundaries. The attribute “weighted” in the name WSBM stems from the fact that its variational form is weighted with the elemental volume fraction of active fluid. The motivation for the development of the WSBM was the preservation of the volume of active fluid to a higher degree of accuracy, which in turn resulted in improved stability and robustness characteristics in moving-boundary, time-dependent simulations. In this article, we present the numerical analysis of the WSBM formulations for the Poisson and Stokes problems. We give mathematical conditions under which the bilinear forms defining the discrete variational formulations are uniformly coercive (Poisson problem) or inf-sup stable (Stokes problem). By these results, stability and optimal convergence is proven in the natural norm; L2-error estimates can also be derived.

Approximate domain boundaries↗
DOE OSTI2025

Optical neural engine for solving scientific partial differential equations

Abstract Solving partial differential equations (PDEs) is the cornerstone of scientific research and development. Data-driven machine learning (ML) approaches are emerging to accelerate time-consuming and computation-intensive numerical simulations of PDEs. Although optical systems offer high-throughput and energy-efficient ML hardware, their demonstration for solving PDEs is limited. Here, we present an optical neural engine (ONE) architecture combining diffractive optical neural networks for Fourier space processing and optical crossbar structures for real space processing to solve time-dependent and time-independent PDEs in diverse disciplines, including Darcy flow equation, the magnetostatic Poisson’s equation in demagnetization, the Navier-Stokes equation in incompressible fluid, Maxwell’s equations in nanophotonic metasurfaces, and coupled PDEs in a multiphysics system. We numerically and experimentally demonstrate the capability of the ONE architecture, which not only leverages the advantages of high-performance dual-space processing for outperforming traditional PDE solvers and being comparable with state-of-the-art ML models but also can be implemented using optical computing hardware with unique features of low-energy and highly parallel constant-time processing irrespective of model scales and real-time reconfigurability for tackling multiple tasks with the same architecture. The demonstrated architecture offers a versatile and powerful platform for large-scale scientific and engineering computations.

Tang, Yingheng (ORCID:0009000153622546)↗
NASA NTRS1974

Numerical calculations of velocity and pressure distribution around oscillating airfoils

An analytical procedure based on the Navier-Stokes equations was developed for analyzing and representing properties of unsteady viscous flow around oscillating obstacles. A variational formulation of the vorticity transport equation was discretized in finite element form and integrated numerically. At each time step of the numerical integration, the velocity field around the obstacle was determined for the instantaneous vorticity distribution from the finite element solution of Poisson's equation. The time-dependent boundary conditions around the oscillating obstacle were introduced as external constraints, using the Lagrangian Multiplier Technique, at each time step of the numerical integration. The procedure was then applied for determining pressures around obstacles oscillating in unsteady flow. The obtained results for a cylinder and an airfoil were illustrated in the form of streamlines and vorticity and pressure distributions.

Bratanow, T.↗
NASA NTRS1974

Analysis of three-dimensional unsteady flow around oscillating wings

A method based on the Navier-Stokes equations was developed for determining analytically the three-dimensional unsteady flow patterns around oscillating wings. The Helmholz vorticity transport equations were discretized in three-dimensional finite element form from a variational formulation and integrated numerically. At each time step of the numerical integration the velocity field was calculated from the representation of the three-dimensional wing by a system of optimized distribution of vortices in space. During the numerical integration of the vorticity transport equations the time-dependent boundary conditions on the wing were specified as external constraint conditions. Examples of obtained results describing the three-dimensional unsteady flow around a wing were presented.

Bratanow, T.↗
NASA NTRS1977

An experimental and computational investigation of the steady and unsteady transonic flow fields about an airfoil in a solid-wall test channel

An experimental and computational investigation of the steady and unsteady transonic flow field about a thick airfoil is described. An operational computer code for solving the two-dimensional, compressible Navier-Stokes equations for flow over airfoils was modified to include solid-wall, slipflow boundary conditions to properly assess the code and help guide the development of improved turbulence models. Steady and unsteady flow fields about an 18% thick circular arc airfoil at Mach numbers of 0.720, 0.754, and 0.783 and a chord Reynolds number of 11 x 10 to the 6th are predicted and compared with experiment. For the first time, computed results for unsteady turbulent flows with separation caused by a shock wave were obtained which qualitatively reproduce the time-dependent aspects of experiments. Features such as the intensity and reduced frequency of airfoil surface-pressure fluctuations, oscillatory regions of trailing-edge and shock-induced separation, and the Mach number range for unsteady flows were all qualitatively reproduced.

Levy, L. L., Jr.↗
NASA NTRS1981

VNAP2: A Computer Program for Computation of Two-dimensional, Time-dependent, Compressible, Turbulent Flow

A computer program, VNAP2, for calculating turbulent (as well as laminar and inviscid), steady, and unsteady flow is presented. It solves the two dimensional, time dependent, compressible Navier-Stokes equations. The turbulence is modeled with either an algebraic mixing length model, a one equation model, or the Jones-Launder two equation model. The geometry may be a single or a dual flowing stream. The interior grid points are computed using the unsplit MacCormack scheme. Two options to speed up the calculations for high Reynolds number flows are included. The boundary grid points are computed using a reference plane characteristic scheme with the viscous terms treated as source functions. An explicit artificial viscosity is included for shock computations. The fluid is assumed to be a perfect gas. The flow boundaries may be arbitrary curved solid walls, inflow/outflow boundaries, or free jet envelopes. Typical problems that can be solved concern nozzles, inlets, jet powered afterbodies, airfoils, and free jet expansions. The accuracy and efficiency of the program are shown by calculations of several inviscid and turbulent flows. The program and its use are described completely, and six sample cases and a code listing are included.

Cline, M. C.↗
NASA NTRS1980

Solution of 3-dimensional time-dependent viscous flows. Part 2: Development of the computer code

There is considerable interest in developing a numerical scheme for solving the time dependent viscous compressible three dimensional flow equations to aid in the design of helicopter rotors. The development of a computer code to solve a three dimensional unsteady approximate form of the Navier-Stokes equations employing a linearized block emplicit technique in conjunction with a QR operator scheme is described. Results of calculations of several Cartesian test cases are presented. The computer code can be applied to more complex flow fields such as these encountered on rotating airfoils.

Weinberg, B. C.↗
NASA NTRS1982

Solution of 3-dimensional time-dependent viscous flows. Part 3: Application to turbulent and unsteady flows

A numerical scheme is developed for solving the time dependent, three dimensional compressible viscous flow equations to be used as an aid in the design of helicopter rotors. In order to further investigate the numerical procedure, the computer code developed to solve an approximate form of the three dimensional unsteady Navier-Stokes equations employing a linearized block implicit technique in conjunction with a QR operator scheme is tested. Results of calculations are presented for several two dimensional boundary layer flows including steady turbulent and unsteady laminar cases. A comparison of fourth order and second order solutions indicate that increased accuracy can be obtained without any significant increases in cost (run time). The results of the computations also indicate that the computer code can be applied to more complex flows such as those encountered on rotating airfoils. The geometry of a symmetric NACA four digit airfoil is considered and the appropriate geometrical properties are computed.

Weinberg, B. C.↗
NASA NTRS1989

On the Nature of Navier-stokes Turbulence

Several turbulent and nonturbulent solutions of the Navier-Stokes equations are obtained. The unaveraged equations are used numerically in conjunction with tools and concepts from nonlinear dynamics, including time series, phase portraits, Poincare sections, largest Liapunov exponents, power spectra, and strange attractors. Initially neighboring solutions for a low-Reynolds-number fully developed turbulence are compared. The solutions, separate exponentially with time, having a positive Liapunov exponent. Thus the turbulence is characterized as chaotic. In a search for solutions which contrast with the turbulent ones, the Reynolds number is reduced. Several qualitatively different flows are noted. These are, fully chaotic, complex period, weakly chaotic, simple periodic, and fixed-point. Of these, only the fully chaotic flows are classified as turbulent. Those flows have both a positive Liapunov exponent and Poincare sections without pattern. By contrast, the weakly chaotic flows have some pattern in their Poincare sections. The fixed-point and periodic flows are nonturbulent, since turbulence, is both time-dependent and aperiodic. Turbulent solutions are obtained in which energy cascades from large to small-scale motions. In general, the spectral energy transfer takes place between wavenumber bands that are considerably separated. The special transfer can occur either as a result of nonlinear turbulence self-interaction or by interaction of turbulence with mean gradients. Turbulent systems are compared with those studied in kinetic theory. The two types of systems are fundamentally different (continuous and dissipative as opposed to discrete and conservative), but there are similarities. For instance, both are nonlinear and show sensitive dependence on initial conditions. Also, the turbulent and molecular stress tensors are identical if the macroscopic velocities for the turbulent stress are replaced by molecular velocities.

Deissler, Robert G.↗
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 NTRS1992

Fractional-Step, Finite-Volume Computation Of Flow

Method of solving Navier-Stokes equations of incompressible flow in general nonorthogonal curvilinear coordinates incorporates fractional-step and finite-volume approaches. Developed in continuing effort to achieve accuracy without need for excessive computation time in numerical simulation of time-dependent, three-dimensional flows bounded by surfaces of complicated shape.

Kwak, Dochan↗
NASA NTRS1991

Time-dependent solutions of viscous incompressible flows in moving co-ordinates

A time-accurate solution method for the incompressible Navier-Stokes equations in generalized moving coordinates is presented. A finite-volume discretization method that satisfies the geometric conservation laws for time-varying computational cells is used. The discrete equations are solved by a fractional-step solution procedure. The solution is second-order-accurate in space and first-order-accurate in time. The pressure and the volume fluxes are chosen as the unknowns to facilitate the formulation of a consistent Poisson equation and thus to obtain a robust Poisson solver with favorable convergence properties. The method is validated by comparing the solutions with other numerical and experimental results. Good agreement is obtained in all cases.

Rosenfeld, Moshe↗
NASA NTRS1992

A finite-volume numerical method to calculate fluid forces and rotordynamic coefficients in seals

A numerical method to calculate rotordynamic coefficients of seals is presented. The flow in a seal is solved by using a finite-volume formulation of the full Navier-Stokes equations with appropriate turbulence models. The seal rotor is perturbed along a diameter such that the position of the rotor is a sinusoidal function of time. The resulting flow domain changes with time, and the time-dependent flow in the seal is solved using a space conserving moving grid formulation. The time-varying fluid pressure reaction forces are then linked with the rotor center displacement, velocity and acceleration to yield the rotordynamic coefficients. Results for an annular seal are presented, and compared with experimental data and other more simplified numerical methods.

Athavale, M. M.↗
NASA NTRS1994

Boundary Conditions for Jet Flow Computations

Ongoing activities are focused on capturing the sound source in a supersonic jet through careful large eddy simulation (LES). One issue that is addressed is the effect of the boundary conditions, both inflow and outflow, on the predicted flow fluctuations, which represent the sound source. In this study, we examine the accuracy of several boundary conditions to determine their suitability for computations of time-dependent flows. Various boundary conditions are used to compute the flow field of a laminar axisymmetric jet excited at the inflow by a disturbance given by the corresponding eigenfunction of the linearized stability equations. We solve the full time dependent Navier-Stokes equations by a high order numerical scheme. For very small excitations, the computed growth of the modes closely corresponds to that predicted by the linear theory. We then vary the excitation level to see the effect of the boundary conditions in the nonlinear flow regime.

Hayder, M. E.↗
NASA NTRS1994

Three-dimensional fully spectral numerical method for mantle convection with depth-dependent properties

A semi-implicit fully spectral collocation method for the simulation of three-dimensional mantle convection with depth-dependent thermo-dynamic and transport properties is presented. The variable property Navier-Stokes equation expressed in terms of the primitive variable velocity and pressure is solved with the mass continuity and temperature equations. The periodic horizontal boundary conditions allow a Fourier expansion for the two horizontal directions. The stress-free, impermeable isothermal boundary conditions along with the depth dependent coefficients are handled with a Chebyshev expansion in the vertical direction. In the limit of an infinite Prandtl number appropriate to mantle convection, the inertial terms in the momentum equation are unimportant. In this case an explicit solution of a Poisson equation for pressure can be avoided; instead a fourth-order equation for vertical velocity can be solved. Simultaneous imposition of both impermeable and continuity boundary conditions during the vertical velocity evaluation is discussed. The pressure distributions on the top and bottom bounding planes were determined by means of an influence matrix technique. The numerical method employed here avoids time-splitting errors and enforces velocity boundary conditions and continuity over the entire domain, including the boundaries, to machine accuracy. Strongly time-dependent three-dimensional solutions up to a surface Rayleigh number of 1 x 10(exp 7) have been obtained. Strong upwellings, pulsating chaotically, are formed by the collective merging of cylindrical plumes.

Balachandar, S.↗