SEARCH · Engineering Papers

Results for “Stokes flow”

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 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Numerical solutions of Navier-Stokes equations for a Butler wing

The flow field is simulated on the surface of a given delta wing (Butler wing) at zero incident in a uniform stream. The simulation is done by integrating a set of flow field equations. This set of equations governs the unsteady, viscous, compressible, heat conducting flow of an ideal gas. The equations are written in curvilinear coordinates so that the wing surface is represented accurately. These equations are solved by the finite difference method, and results obtained for high-speed freestream conditions are compared with theoretical and experimental results. In this study, the Navier-Stokes equations are solved numerically. These equations are unsteady, compressible, viscous, and three-dimensional without neglecting any terms. The time dependency of the governing equations allows the solution to progress naturally for an arbitrary initial initial guess to an asymptotic steady state, if one exists. The equations are transformed from physical coordinates to the computational coordinates, allowing the solution of the governing equations in a rectangular parallel-piped domain. The equations are solved by the MacCormack time-split technique which is vectorized and programmed to run on the CDC VPS 32 computer.

Abolhassani, J. S.↗

NASA NTRS1991

Hypersonic waverider analysis - A comparison of numerical and experimental results

In this study two waverider configurations are analyzed numerically and experimentally. The geometries are designed using an inverse design and optimization code with the optimization procedure driven by maximized L/D and limited by geometric constraints. Euler and Navier-Stokes numerical simulations are performed with a partially flux-split finite-difference code developed at NASA Ames, F3D. Experimental data is obtained at NASA Langley in the Unitary Plan Wind Tunnel (UPWT) and the 20-Inch Mach 6 wind tunnel for flow visualization, force/moment, and surface pressure measurements at both on- and off-design conditions. Euler simulations are compared with analytic results obtained from conical flow theory, and Navier-Stokes simulations are compared with experimental results both qualitatively and quantitatively for on- and off-design conditions.

Jones, K. D.↗

DOE OSTI2021

Beyond Taylor’s hypothesis: a novel volumetric reconstruction of velocity and density fields for variable‑density and shear flows

Abstract This work presents a novel numerical procedure for reconstructing volumetric density and velocity fields from planar laser-induced fluorescence (PLIF) and stereoscopic particle image velocimetry (SPIV) data. This new method is theoretically and practically demonstrated to provide more accurate 3D vortical structures and density fields in high shear flows than reconstruction methods based on the mean convective velocity. While Taylor’s hypothesis of frozen turbulence is commonly applied by using the local mean streamwise velocity, the proposed algorithm uses the measured local instantaneous velocity for data convection. It consists of a step-by-step reconstruction based on a mixed Lagrangian–Eulerian solver that includes the 3D interpolation of scattered flow data and that relaxes the Taylor’s hypothesis by iterative enforcement of the incompressibility constraint on the velocity field. This methodology provides 3D fields with temporal resolution, spatial resolution, and accuracy comparable to that of real 3D snapshots, thus providing a practical alternative to tomographic measurements. The procedure is validated using numerical data of the constant-density channel flow available on the Johns Hopkins University Turbulence Database (JHTDB), showing the accurate reconstruction of the 3D velocity field. The algorithm is applied to an experimental dataset of PLIF and SPIV measurements of a variable-density jet flow, demonstrating its capability to provide 3D velocity and density fields that are more consistent with the Navier–Stokes equations compared to the mean flow convective method. Graphic abstract

42 ENGINEERING↗

NASA NTRS1975

Asymptotic analytical methods in fluid mechanics related to drag prediction

Some recent theoretical work of a purely analytical nature is described which promises to provide engineering predictions for the important drag-related phenomena of flow in the stall regime. This analytical work deals with rigorous asymptotic studies of the complete Navier-Stokes equations that govern the viscous flow around any aerodynamic body under conditions where boundary layer separation takes place from the body surface.

Inger, G. R.↗

NASA NTRS1981

Numerical simulation of a compressible homogeneous, turbulent shear flow

A direct, low Reynolds number, numerical simulation was performed on a homogeneous turbulent shear flow. The full compressible Navier-Stokes equations were used in a simulation on the ILLIAC IV computer with a 64,000 mesh. The flow fields generated by the code are used as an experimental data base, to examine the behavior of the Reynols stresses in this simple, compressible flow. The variation of the structure of the stresses and their dynamic equations as the character of the flow changed is emphasized. The structure of the tress tensor is more heavily dependent on the shear number and less on the fluctuating Mach number. The pressure-strain correlation tensor in the dynamic uations is directly calculated in this simulation. These correlations are decomposed into several parts, as contrasted with the traditional incompressible decomposition into two parts. The performance of existing models for the conventional terms is examined, and a model is proposed for the 'mean fluctuating' part.

Feiereisen, W. J.↗

NASA NTRS1987

Stability of time dependent and spatially varying flows; Proceedings of the Symposium, Hampton, VA, Aug. 19-23, 1985

Papers are presented on the application of stability theory to laminar flow control, secondary instabilities in boundary layers, a Floquet analysis of secondary instability in shear flows, and the generation of Tollmien-Schlichting waves by long wavelength free stream disturbances. Also considered are numerical experiments on boundary-layer receptivity, short-scale inviscid instabilities in the flow past surface-mounted obstacles, wave phenomena in a high Reynolds number compressible boundary layer, and instability of time-periodic flows. Other topics include high frequency Rayleigh instability of Stokes layers, stability and resonance in grooved-channel flows, finite length Taylor Couette flow, and vortical structures in the breakdown stage of transition.

Dwoyer, Douglas L.↗

NASA NTRS1991

Pressure flux-split viscous solutions for subsonic diffusers

A method for predicting flow in diffusers with inlet swirl has been developed. Solutions of the pressure flux-split Reduced Navier Stokes (RNS) equations are obtained for flows in axisymmetric diffusers. Viscous-Inviscid interactions and flow fields with toroidal recirculation regions are efficiently captured. The computational model is verified by comparision with experimental data and other computations. The extreme sensitivity to grid, turbulent closure model and inlet profiles is discussed.

Pordal, H. S.↗

NASA NTRS2016

Enhancement of the Open National Combustion Code (OpenNCC) and Initial Simulation of Energy Efficient Engine Combustor

In this paper, we present the recent enhancement of the Open National Combustion Code (OpenNCC) and apply the OpenNCC to model a realistic combustor configuration (Energy Efficient Engine (E3)). First, we perform a series of validation tests for the newly-implemented advection upstream splitting method (AUSM) and the extended version of the AUSM-family schemes (AUSM+-up). Compared with the analytical/experimental data of the validation tests, we achieved good agreement. In the steady-state E3 cold flow results using the Reynolds-averaged Navier-Stokes(RANS), we find a noticeable difference in the flow fields calculated by the two different numerical schemes, the standard Jameson- Schmidt-Turkel (JST) scheme and the AUSM scheme. The main differences are that the AUSM scheme is less numerical dissipative and it predicts much stronger reverse flow in the recirculation zone. This study indicates that two schemes could show different flame-holding predictions and overall flame structures.

Combustor-turbine interaction↗

NASA NTRS2016

Enhancement of the Open National Combustion Code (OpenNCC) and Initial Simulation of Energy Efficient Engine Combustor

Code validation↗

NASA NTRS1993

Superfluid turbulence

At low temperatures (below 5 Kelvin), helium is a liquid with a very low kinematic viscosity. It was proposed that wind tunnels could be built using liquid helium as the test fluid. The primary advantages of such wind tunnels would be a combination of large Reynolds numbers and a relatively small apparatus. It is hoped that this combination will allow the study of high Reynolds number flows in an academic setting. There are two basic types of liquid helium wind tunnels that can be built, corresponding to the two phases of liquid helium. The high temperature phase (between approximately 2 to 5 Kelvin) is called helium 1 and is a Navier-Stokes fluid. There are no unanswered scientific questions about the design or operation of a wind tunnel in the helium 1 phase. The low temperature phase (below approximately 2 Kelvin) of liquid helium is called helium 2. This is a quantum fluid, meaning that there are some properties of helium 2 which are directly due to quantum mechanical effects and which are not observed in Navier-Stokes fluids. The quantum effects that are relevant to this paper are: (1) helium 2 is well described as a superposition of two separate fluids called the superfluid and the normal fluid. The normal-fluid component is a Navier-Stokes fluid and the superfluid is an irrotational Euler fluid; and (2) circulation in the superfluid exists only in quantized vortex filaments. All quantized vortex filaments have identical circulations kappa and core size a. The objective of the research at CTR was to develop an understanding of the microscopic processes responsible for the observed Navier-Stokes behavior of helium 2 flows.

Samuels, David C.↗

NASA NTRS1983

Computation of supersonic viscous flows around pointed bodies at large incidence

A recently reported parabolized Navier-Stokes method has been extended to compute turbulent supersonic flows around cones and an ogive-cylinder body at large incidence. The algebraic eddy-viscosity turbulence model contained in the code was modified to properly account for the large regions of cross-flow separation that occur in these flows. Extensive comparisons between computed results and experimentally measured flow fields are presented. The results show good agreement for viscous-layer profiles and details of the external leeward-side vortex structure at angles of attack up to three times the cone half angles. Details of the modified turbulence model are presented and discussed.

Degani, D.↗

NASA NTRS1994

Turbulence models and Reynolds analogy for two-dimensional supersonic compression ramp flow

Results of the application of turbulence models and the Reynolds analogy to the Navier-Stokes computations of Mach 2.9 two-dimensional compression ramp flows are presented. The Baldwin-Lomax eddy viscosity model and the kappa-epsilon turbulence transport equations for the turbulent momentum flux modeling in the Navier-Stokes equations are studied. The Reynolds analogy for the turbulent heat flux modeling in the energy equation was also studied. The Navier-Stokes equations and the energy equation were numerically solved for the flow properties. The Reynolds shear stress, the skin friction factor, and the surface heat transfer rate were calculated and compared with their measurements. It was concluded that with a hybrid kappa-epsilon turbulence model for turbulence modeling, the present computations predicted the skin friction factors of the 8 deg and 16 deg compression ramp flows and with the turbulent Prandtl number Pr(sub t) = 0.93 and the ratio of the turbulent thermal and momentum transport coefficients mu(sub q)/mu(sub t) = 2/Prt, the present computations also predicted the surface heat transfer rates beneath the boundary layer flow of the 16 compression ramp.

Wang, Chi R.↗

NASA NTRS1992

Calculation of two-dimensional shear-driven cavity flows at high Reynolds numbers

The time-dependent Navier-Stokes equations are numerically integrated for two-dimensional incompressible viscous flow in a shear-driven square cavity. Using a time-splitting method and finite differences on a staggered mesh, the momentum and pressure equations are directly solved by a tensor product method where one finite difference direction is diagonalized by eigenvalue decomposition. The effects of increasing Reynolds number are studied and the developing boundary layer is captured by using a finely clustered mesh. At Re = 30,000 the flow is in a continuously developing unsteady regime. Power spectrum plots indicate that the unsteady flow oscillates with one fundamental frequency and exhibits some characteristics of transition between laminar and turbulent states.

Huser, A.↗

NASA NTRS2022

Combined Bluntness and Roughness Effects on Cones at Hypersonic Speeds

This computational study investigates the effects of discrete roughness elements on a blunt cone at zero degrees angle of attack in a Mach 6 flow. Motivation was provided by experiments conducted in the Air Force Research Laboratory Mach 6 High Reynolds Number facility on a 7-degree half-angle cone with a roughness array located at 45 degrees from the apex on two nose tips of different blutness but equivalent roughness Reynolds number. Transition was only affected on the blunter cone, indicating that the transition onset is associated with the combined effects of bluntness and roughness. The present study investigates the 15.24 mm nose radius, 420 azimuthal wavenumber case via Navier-Stokes computations of the laminar base flow and instability analysis. Plane-marching parabolized stability equations (PSE) and inflow-resolvent analysis based on the three-dimensional, harmonic linearized Navier-Stokes equations (HLNSE) are used to calculate the amplification of disturbances along the roughness wake as well as over the roughness nearfield. Results show that the roughness shape can have a great impact on the characteristics of the most amplified wake instabilities. For the experimental configuration with cubic roughness elements of 15 μ m height, the flow is marginally unstable. For prismatic elements of 20 μ m height, the PSE predicts a logarithmic disturbance amplification ratio of N = 5.6 along its wake, but this ratio increases to N = 9.6 when the amplification over the roughness and separation regions is included in the inflow-resolvent analysis.

Boundary-layer transition↗

NASA NTRS1981

Computation of hypersonic viscous flow around three-dimensional bodies at high angles of attack

A parabolized Navier-Stokes code capable of predicting steady viscous supersonic flows with cross-flow separation is applied to three-dimensional arbitrary geometries at high angles of attack. The numerical procedure, which is implicit, noniterative, and of second-order accuracy in the marching direction, has been used to compute complicated flow fields containing a relatively thick sonic layer and regions of strong viscous-inviscid interaction. A consistent and accurate procedure has also been developed to provide the necessary starting data through timewise integration of the equations of motion near the nose-tip region of the body. Numerical results obtained from the present method compare well with experiment for both the surface pressures and heat transfer.

Rizk, Y. M.↗

NASA NTRS1983

Transonic, Axisymmetric Flow Over Nozzle Afterbodies With Supersonic Jet Exhausts

Predictions require less computation than Navier-Stokes solutions. RAXJET computer program predicts transonic, axisymmetric flow over nozzle afterbodies with supersonic jet exhausts and includes effects of boundarylayer displacement, separation, jet entrainment, and inviscid jet plume blockage. RAXJET written in FORTRAN IV.

Wilmoth, R. G.↗

NASA NTRS1991

Wake interference behind two flat plates normal to the flow - A finite-element study

A finite-element model of the Navier-Stokes equations is used for numerical simulation of flow past two normal flat plates arranged side by side at Reynolds number 80 and 160. The results from this simulation indicate that when the gap between the plates is twice the width of a single plate, the individual wakes of the plates behave independently, with the antiphase vortex shedding being dominant. At smaller gap sizes, the in-phase vortex shedding, with strong wake interaction, is favored. The gap flow in those cases becomes biased, with one of the wakes engulfing the other. The direction of the biased flow was found to be switching at irregular intervals, with the time histories of the indicative flow parameters and their power spectra resembling those of a chaotic system.

Behr, M.↗

DOE OSTI2020

Evaluating the Effect of Vane Trailing Edge Flow on Turbine Rim Sealing

Abstract Modern gas turbine development continues to move toward increased overall efficiency, driven in part by higher firing temperatures that point to a need for more cooling air to prevent catastrophic component failure. However, using additional cooling flow bled from the upstream compressor causes a corresponding detriment to overall efficiency. A primary candidate for cooling flow optimization is purge flow, which contributes to sealing the stator–rotor cavity and prevents ingestion of hot main gas path (MGP) flow into the wheelspace. Previous research has identified that the external main gas path flow physics play a significant role in driving rim seal ingestion. However, the potential impact of other cooling flow features on ingestion behavior, such as vane trailing edge (VTE) flow, is absent in the open literature. This paper presents experimental measurements of rim cavity cooling effectiveness collected from a one-stage turbine operating at engine-representative Reynolds and Mach numbers. Carbon dioxide (CO2) was used as a tracer gas in both the purge flow and vane trailing edge flow to investigate flow migration into and out of the wheelspace. Results show that the vane trailing edge flow does in fact migrate into the rim seal and that there is a superposition relationship between individual cooling flow contributions. Computational fluid dynamics (CFD) simulations using unsteady Reynolds-averaged Navier–Stokes (URANS) were used to confirm VTE flow ingestion into the rim seal cavity. Radial and circumferential traverse surveys were performed to quantify cooling flow radial migration through the main gas path with and without vane trailing edge flow. The surveys confirmed that vane trailing edge flow is entrained into the wheelspace as purge flow is reduced. Local CO2 measurements also confirmed the presence of VTE flow deep in the wheelspace cavity.

Engineering↗