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 37 records · Page 2

An experimental and computational investigation of the flow field about a transonic airfoil in supercritical flow with turbulent boundary-layer separation

A combined experimental and computational research program is described for testing and guiding turbulence modeling within regions of separation induced by shock waves incident in turbulent boundary layers. Specifically, studies are made of the separated flow the rear portion of an 18%-thick circular-arc airfoil at zero angle of attack in high Reynolds number supercritical flow. The measurements include distributions of surface static pressure and local skin friction. The instruments employed include highfrequency response pressure cells and a large array of surface hot-wire skin-friction gages. Computations at the experimental flow conditions are made using time-dependent solutions of ensemble-averaged Navier-Stokes equations, plus additional equations for the turbulence modeling.

Rubesin, M. W.↗

An experimental and computational investigation of the flow field about a transonic airfoil in supercritical flow with turbulent boundary-layer separation

A combined experimental and computational research program for testing and guiding turbulence modeling within regions of separation induced by shock waves incident on turbulent boundary layers is described. Specifically, studies are made of the separated flow over the rear portion of an 18%-thick circular-arc airfoil at zero angle of attack in high Reynolds number supercritical flow. The measurements include distributions of surface static pressure and local skin friction. The instruments employed include high-frequency response pressure cells and a large array of surface hot-wire skin-friction gages. Computations at the experimental flow conditions are made using time-dependent solutions of ensemble-averaged Navier-Stokes equations, plus additional equations for the turbulence modeling.

Rubesin, M. W.↗

Viscous reacting flows with wall slip and catalysis applied to spheres in arc jets and flight

The influence of wall slip and catalytic atom-recombination on the flow field and wall heat flux are calculated for high altitude flight and arc jet flow conditions. Boundary equations, which include velocity slip, temperature jump, and wall catalytic atom recombination, are coupled to the viscous reacting multicomponent Navier-Stokes equations. These equations are solved using a time-dependent finite difference technique applied to spheres in an arc jet flow (Reynolds number of 550) and a high altitude flight case representative of the Space Shuttle Orbiter (Reynolds number of 450). The results indicate that catalysis strongly influences the temperature jump, but not the velocity slip. Slip increases the atom fraction and temperature at both the wall and the flow field. Likewise, the shock stand-off distance, the wall heat flux, and friction coefficient are increased over the nonslip cases. The reacting gas calculations confirm the chemically frozen nature of the shock layer in arc jet flows.

Scott, C. D.↗

Interaction of strong shock wave with turbulent boundary layer

The reported investigation represents an extension of the time-dependent solution of separated laminar flows based on the complete Navier-Stokes equations reported by MacCormack (1971) and Carter (1973). The current study includes turbulence models in conjunction with the compressible flow equations. The calculations start with a uniform flow except for values imposed along the upstream and outer boundaries. The basic numerical method is discussed along with questions concerning the exponential accuracy and the resolution of the viscous sublayer in a compressed region.

Baldwin, B. S.↗

Numerical computation of two dimensional viscous blunt body flows with an impinging shock, part 2

Two-dimensional viscous blunt body flows with an impinging shock have been computed using a time-dependent finite-difference method which solves the complete set of Navier-Stokes equations for a compressible flow. For low Reynolds number flows, the entire flow field, including the bow shock and impinging shock, has been captured in the computation. For higher Reynolds number flows, the bow shock is treated as a discontinuity across which the Rankine-Hugoniot equations are applied, while the boundary layer and interaction regions are captured as before. Using this latter shock-fitting approach, a Type III shock interaction flow field has been computed with flow conditions corresponding to the space shuttle orbiter freestream conditions at 61 km (200,000 ft).

Holst, T. L.↗

Numerical computation of viscous blunt body flows with a planar impinging shock

Two- and three-dimensional, viscous blunt body flows with planar impinging shocks are computed using an explicit, time-dependent, finite-difference method to solve the complete set of Navier-Stokes equations. The bow shock is treated as a discontinuity, while all interior shock layer detail such as shear layers, shock waves, jets and the wall boundary layer are automatically captured in the solution. Numerical results are presented for cases in which planar shock waves of different strengths and orientations are allowed to impinge on the flow field surroundings an infinite cylinder resulting in two- and three-dimensional shock interference patterns. The numerical results are compared with experiment.

Holst, T. L.↗

Numerical computation of two-dimensional viscous blunt body flows with an impinging shock

Two-dimensional, viscous, blunt body flows with an impinging shock wave are computed using a time-dependent, finite-difference method to solve the complete set of Navier-Stokes equations. The bow shock wave is treated as a discontinuity, while all interior shock layer detail such as shear layers, shock waves, jets, and the wall boundary layer are automatically captured in the solution. Numerical results are presented for cases in which shock waves of different strengths are allowed to impinge on the flow field surrounding a circular cylinder resulting in different shock interference patterns. The two-dimensional results are compared qualitatively with existing three-dimensional experiments.

Tannehil, J. C.↗

Comparison of a two-dimensional shock impingement computation with experiment

Results of computations of two-dimensional viscous blunt-body flowfields with an impinging shock wave, with a time-dependent finite-difference method employed to solve the complete set of Navier-Stokes equations, are compared with experimental results. The experimental results were obtained in a 20-inch hypersonic tunnel with a planar shock impinging on the cylindrical leading edge of a fin, hence with the shock parallel to the centerline of the leading edge, so that type III and type IV interference patterns were generated. Close agreement is found. The overall effects of smoothing and grid size on the calculations are determined. A 31 x 51 mesh is adequate for wall pressure values (except in peaked regions).

Tannehill, J. C.↗

Calculation of the flow on a blunted cone at high angle of attack

A new technique for calculating the entire flow-field on spherically blunted circular cones at high angles of attack and high Reynolds numbers is described. The calculations are based on a single-layer system of three-dimensional parabolic equations which are approximations to the full steady Navier-Stokes equations. Initial conditions at the sphere-cone tangency plane are provided by using an inviscid time-dependent solution added to a viscous nonsimilar boundary layer solution. Calculated results are compared with experimental heat transfer and pressure data for a 15 deg half-angle cone with a 1-in. spherical nose at 15 deg angle of attack. The free-stream Mach number is 10.6, and the free-stream Reynolds number is 1,200,000 per foot. Excellent agreement between the calculated and experimental data for both pressure and heat transfer is obtained.

Lubard, S. C.↗

An experimental and analytical study of slip and catalytic boundary conditions applied to spheres in low Reynolds number arc jet flows

Measurements of low-density dissociated arc-heated flows around low and high surface catalycity spheres are compared with numerical solutions to the full Navier-Stokes equations for reacting gas flow fields, taking into account slip and finite surface catalytic boundary conditions. A time-dependent finite difference solution technique developed by Li (1973, 1974) is used in the investigation and a simple transport model is employed for the flows considered. The correctness of the solution approach is demonstrated by the very good agreement obtained between experimental and computational data.

Scott, C. D.↗

Convection in the tanks of a rotating spacecraft

A general study of convection and mixing of a stratified fluid in a rotating container is described, with application to the special problem of fluid heating and convection in a spacecraft tank. The analysis was based on a set of approximate equations for the Navier-Stokes description of fluid convection with small density variations in a rotating system, including effects of body forces due to temperature stratification (caused by a heater) and arbitrary time-dependent rotation of the tank about a noncentral axis. An efficient numerical finite difference scheme and computational method are described for the convection of vorticity and energy in a two-dimensional tank. Special procedures were developed for analysis of the thermodynamic states resulting from the approximate flow equations derived for small density variations. Results of the numerical simulation are presented for studying the effectiveness of rotation maneuvers in mixing stratified oxygen in the tanks of an Apollo spacecraft. Significant effects of the rotation maneuvers are discussed, including the reduction of the potential for pressure collapse.

Source record↗

Finite element analysis of unsteady incompressible flow around an oscillating obstacle of arbitrary shape.

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

Bratanow, T.↗

Calculation of a separated turbulent boundary layer

The properties of a Navier-Stokes solution of a shock-separated turbulent flow over a flat wall are investigated. Refinements of an algebraic relaxation turbulence model previously shown to be of value for the simulation of separated flows are presented. A simplified analysis applicable near an adiabatic wall is developed and used to help verify the accuracy of the numerical solution. Features of the time-dependent response of a turbulent boundary layer to shock impingement are presented.

Baldwin, B.↗

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↗

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)↗

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.↗

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.↗