Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Calculation of Unsteady Transonic Flow Over an Airfoil

An implicit finite-difference solver for either the Euler equations or the "thin-layer" Navier-Stokes equations was used to calculate a transonic flow over the NACA 64A010 airfoil pitching about its one-quarter chord. An unsteady automatic grid-generation procedure that will improve significantly the computational efficiency of various unsteady flow problems is described. The calculated results for both inviscid and viscous flows at Mach number 0.8 over the airfoil oscillating with reduced frequency referenced to one-half chord, 0.2, are compared with experimental data measured in the Ames 11 x 11 ft Transonic Wind Tunnel. Nonlinear, unsteady effects of the flow on the surface pressure variations, shock-wave excursions, and overall airloads are examined. Good agreements between the results of computations and experiments were obtained. In the shock-wave region, however, the results of the viscous-flow computations showed closer agreement with the experimental data.

Chyu, W. J.↗
NASA NTRS1982

Finite analytic numerical solution of axisymmetric Navier-Stokes and energy equations

Convective heat transfer for steady state laminar flow in axisymmetric coordinates is considered. Numerical solutions for flow pattern and temperature distribution are obtained by the finite analytic numerical method applied to the Navier-Stokes equations expressed in terms of vorticity and stream function, and the energy equation. The finite analytic numerical method differs from other numerical methods in that it utilizes a local analytic solution in an element of the problem to construct the total numerical solution. Finite analytic solutions of vorticity, stream function, temperature and heat transfer coefficients for flow with Reynolds number of 5.0, 100.0, 1000.0, and 2000.0, and Prandtl number of 0.1, 1.0, and 10.0 with uniform grid sizes are reported for axisymmetric pipe with sudden expansion and contraction. The wall temperature is considered to be isothermal and differs from the inlet temperature. It is shown that the finite analytic solution is stable, converges rapidly, and simulates the convection of fluid flow accurately since the local analytic solution is capable of simulating automatically the influence of skewed convection through the element boundary on the interior nodal values thereby minimizing the false numerical diffusion.

Chen, C.-J.↗
NASA NTRS1997

Development of a Low-Reynolds Number, Nonlinear kappa-epsilon Model for the Reduced Navier-Stokes Equations

Previous work at NASA LeRC has shown that flow distortions in aircraft engine inlet ducts can be significantly reduced by mounting vortex generators, or small wing sections, on the inside surface of the engine inlet. The placement of the vortex generators is an important factor in obtaining the optimal effect over a wide operating envelope. In this regard, the only alternative to a long and expensive test program which would search out this optimal configuration is a good prediction procedure which could narrow the field of search. Such a procedure has been developed in collaboration with NASA LeRC, and results obtained by NASA personnel indicate that it shows considerable promise for predicting the viscous turbulent flow in engine inlet ducts in the presence of vortex generators. The prediction tool is a computer code which numerically solves the reduced Navier-Stokes equations and so is commonly referred to as RNS3D. Obvious deficiencies in RNS3D have been addressed in previous work. Primarily, it is known that the predictions of the mean velocity field of a turbulent boundary layer flow approaching separation are not in good agreement with data. It was suggested that the use of an algebraic mixing-length turbulence model in RNS3D is at least partly to blame for this. Additionally, the current turbulence model includes an assumption of isotropy which will ultimately fail to capture turbulence-driven secondary flow known to exist in noncircular ducts.

Boger, David A.↗
NASA NTRS1991

Numerical investigations in three-dimensional internal flows

In previous efforts, a two-dimensional full Navier-Stokes (FNS) code (SCRAM2D) was used in a design process that involved parametric modifications of the inlet geometry to arrive at what appeared to be an optimum inlet flowfield that produced a uniform flow at the exit in a very short distance. In these previous studies, the technologies for determining the contours with a 'man-in-the-loop' approach for both the ramp and cowl of the inlet were demonstrated, and nearly shock-free exiting flowfields were shown to be obtainable. The resulting two-dimensional compression contours were then used with swept sidewalls to form a three-dimensional inlet. Then the three-dimensional Navier-Stokes code (SCRAM3D) was used to investigate the inlet's three-dimensional flow. One of the major difficulties encountered in the previous studies was that associated with the relatively long time required to obtain a solution using even the 2D FNS code in the design process. Since one of the goals of high-speed inlet design is to produce inputs to the overall aircraft design in a timely manner, it was proposed for this year's research to examine 2D and 3D viscous flow solver techniques alternative to the NFS codes used to date. Areas of the inlet particularly identified for code speed up are those associated with the forebody and external flow ramp systems of the inlet. In these areas, parabolized, or space-marched, Navier-Stokes codes were proposed to be investigated for their applicability in the design process developed previously. This report describes the results of an investigation into the use of two other codes for analyzing the forebody and inlet ramp systems of high-speed inlets.

Rose, William C.↗
DOE OSTI2024

ExaWind: Open‐source CFD for hybrid‐RANS/LES geometry‐resolved wind turbine simulations in atmospheric flows

Abstract Predictive high‐fidelity modeling of wind turbines with computational fluid dynamics, wherein turbine geometry is resolved in an atmospheric boundary layer, is important to understanding complex flow accounting for design strategies and operational phenomena such as blade erosion, pitch‐control, stall/vortex‐induced vibrations, and aftermarket add‐ons. The biggest challenge with high‐fidelity modeling is the realization of numerical algorithms that can capture the relevant physics in detail through effective use of high‐performance computing. For modern supercomputers, that means relying on GPUs for acceleration. In this paper, we present ExaWind, a GPU‐enabled open‐source incompressible‐flow hybrid‐computational fluid dynamics framework, comprising the near‐body unstructured grid solver Nalu‐Wind, and the off‐body block‐structured‐grid solver AMR‐Wind, which are coupled using the Topology Independent Overset Grid Assembler. Turbine simulations employ either a pure Reynolds‐averaged Navier–Stokes turbulence model or hybrid turbulence modeling wherein Reynolds‐averaged Navier–Stokes is used for near‐body flow and large eddy simulation is used for off‐body flow. Being two‐way coupled through overset grids, the two solvers enable simulation of flows across a huge range of length scales, for example, 10 orders of magnitude going from O(μm) boundary layers along the blades to O(10 km) across a wind farm. In this paper, we describe the numerical algorithms for geometry‐resolved turbine simulations in atmospheric boundary layers using ExaWind. We present verification studies using canonical flow problems. Validation studies are presented using megawatt‐scale turbines established in literature. Additionally presented are demonstration simulations of a small wind farm under atmospheric inflow with different stability states.

17 WIND ENERGY↗
NASA NTRS1989

Computation of turbulent incompressible wing-body junction flow

A three-dimensional incompressible Reynolds-averaged Navier-Stokes solver is presently used in conjunction with a mixing-length turbulence model to characterize the flow around a wing that is mounted on a flat plate, in a wind tunnel, as well as the flow around a support strut within a turnaround duct. Good agreement is found between predicted and observed values of flat-plate static pressure, horseshoe vortex system size, and mean flow velocities in the case of the wing; the case of the strut in a duct is noted to exhibit many of the same overall flow features as the wing/plate.

Burke, R. W.↗
NASA NTRS1989

Experimental and computational studies of dynamic stall

A review of dynamic stall research in progress under the Navy-NASA Joint Institute of Aeronautics is presented. This effort, which includes both experimental and computational studies of the dynamic stall process, is directed toward better understanding and modeling of the fluid flow that occurs on helicopters and aircraft flying in conditions that induce dynamic stall. The results of research now in progress are presented, with discussion of the experimental program on compressibility effects on dynamic stall, related CFD studies of the stall process based on Navier-Stokes modeling, and viscous-inviscid flow modeling of the incipient stall process.

Carr, L. W.↗
NASA NTRS1995

Viscous aerodynamic design using the adjoint variable approach

The use of classical optimal control methods, in particular variational methods, to solve the airfoil optimization problem, by deriving a set of adjoint (costate) equations and boundary conditions has already been done for inviscid (potential and Euler flows) and two dimensional, steady state, incompressible flow governed by the Navier-Stokes equations. The interior and boundary terms of the volume integral have been derived (in this work) for the steady Navier-Stokes equations in three dimensions for a viscous, compressible heat conducting fluid. This can be used to derive the adjoint equations and numerical boundary conditions for general classes of problems and hence paves the way for a solution to the aerodynamic optimization problem for compressible viscous flows. The next steps to the realization of that goal are projected as below. The usual square integral pressure functional as an objective function is being replaced by a more realistic drag functional subject to a lift constraint. The feasibility of attempting the more difficult time dependent problem is being investigated. It remains to get the full system of adjoint equations and boundary conditions with the new functional. The state and adjoint equations must be discretized and coded. An appropriate optimization program must be used (steepest descents seems inadequate) and various known airfoil shapes should be recovered in test cases of the computer program.

DeRise, George↗
NASA NTRS1981

The effect of modeled turbulence on a hypersonic shock layer with massive ablation injection

For the first time, candidate turbulence models are analyzed for shock-layer conditions as severe as those of a Jovian entry. The various models investigated are two standard models developed primarily by Cebeci (1970) and Kendall (1972), and two modifications of the standard models developed primarily by Baldwin and Lomax (1978) and by Nicolet (1979). The analysis is conducted within the context of a new, viscous shock-layer flow-field code, to properly isolate differences in performance among the various models. This code retains all terms in the Navier-Stokes equations necessary to describe the unique flow events, and incorporates computational features that greatly facilitate the analysis of such events. The relative superiority of the turbulence models cannot be determined from this study alone; however, insight is provided into the relative severity of the various models. In addition, attention is focused on an important research area for experimental and theoretical investigation to develop realistic turbulence models.

Menees, G. P.↗
NASA NTRS1997

Three-Dimensional Viscous Alternating Direction Implicit Algorithm and Strategies for Shape Optimization

A gradient-based shape optimization based on quasi-analytical sensitivities has been extended for practical three-dimensional aerodynamic applications. The flow analysis has been rendered by a fully implicit, finite-volume formulation of the Euler and Thin-Layer Navier-Stokes (TLNS) equations. Initially, the viscous laminar flow analysis for a wing has been compared with an independent computational fluid dynamics (CFD) code which has been extensively validated. The new procedure has been demonstrated in the design of a cranked arrow wing at Mach 2.4 with coarse- and fine-grid based computations performed with Euler and TLNS equations. The influence of the initial constraints on the geometry and aerodynamics of the optimized shape has been explored. Various final shapes generated for an identical initial problem formulation but with different optimization path options (coarse or fine grid, Euler or TLNS), have been aerodynamically evaluated via a common fine-grid TLNS-based analysis. The initial constraint conditions show significant bearing on the optimization results. Also, the results demonstrate that to produce an aerodynamically efficient design, it is imperative to include the viscous physics in the optimization procedure with the proper resolution. Based upon the present results, to better utilize the scarce computational resources, it is recommended that, a number of viscous coarse grid cases using either a preconditioned bi-conjugate gradient (PbCG) or an alternating-direction-implicit (ADI) method, should initially be employed to improve the optimization problem definition, the design space and initial shape. Optimized shapes should subsequently be analyzed using a high fidelity (viscous with fine-grid resolution) flow analysis to evaluate their true performance potential. Finally, a viscous fine-grid-based shape optimization should be conducted, using an ADI method, to accurately obtain the final optimized shape.

Pandya, Mohagna J.↗
NASA NTRS2001

Development of Flow and Heat Transfer Models for the Carbon Fiber Rope in Nozzle Joints of the Space Shuttle Reusable Solid Rocket Motor

Methodologies have been developed for modeling both gas dynamics and heat transfer inside the carbon fiber rope (CFR) for applications in the space shuttle reusable solid rocket motor joints. Specifically, the CFR is modeled using an equivalent rectangular duct with a cross-section area, friction factor and heat transfer coefficient such that this duct has the same amount of mass flow rate, pressure drop, and heat transfer rate as the CFR. An equation for the friction factor is derived based on the Darcy-Forschheimer law and the heat transfer coefficient is obtained from pipe flow correlations. The pressure, temperature and velocity of the gas inside the CFR are calculated using the one-dimensional Navier-Stokes equations. Various subscale tests, both cold flow and hot flow, have been carried out to validate and refine this CFR model. In particular, the following three types of testing were used: (1) cold flow in a RSRM nozzle-to-case joint geometry, (2) cold flow in a RSRM nozzle joint No. 2 geometry, and (3) hot flow in a RSRM nozzle joint environment simulator. The predicted pressure and temperature history are compared with experimental measurements. The effects of various input parameters for the model are discussed in detail.

Wang, Q.↗
DOE OSTI2020

Blood flow imaging by optimal matching of computational fluid dynamics to 4D-flow data

Three-dimensional, time-resolved blood flow measurement (4D-flow) is a powerful research and clinical tool, but improved resolution and scan times are needed. Therefore, this study aims to (1) present a postprocessing framework for optimization-driven simulation-based flow imaging, called 4D-flow High-resolution Imaging with a priori Knowledge Incorporating the Navier-Stokes equations and the discontinuous Galerkin method (4D-flow HIKING), (2) investigate the framework in synthetic tests, (3) perform phantom validation using laser particle imaging velocimetry, and (4) demonstrate the use of the framework in vivo. An optimizing computational fluid dynamics solver including adjoint-based optimization was developed to fit computational fluid dynamics solutions to 4D-flow data. Synthetic tests were performed in 2D, and phantom validation was performed with pulsatile flow. Reference velocity data were acquired using particle imaging velocimetry, and 4D-flow data were acquired at 1.5 T. In vivo testing was performed on intracranial arteries in a healthy volunteer at 7 T, with 2D flow as the reference. Results Synthetic tests showed low error (0.4%-0.7%). Phantom validation showed improved agreement with laser particle imaging velocimetry compared with input 4D-flow in the horizontal (mean -0.05 vs -1.11 cm/s, P < .001; SD 1.86 vs 4.26 cm/s, P < .001) and vertical directions (mean 0.05 vs -0.04 cm/s, P = .29; SD 1.36 vs 3.95 cm/s, P < .001). In vivo data show a reduction in flow rate error from 14% to 3.5%. Phantom and in vivo results from 4D-flow HIKING show promise for future applications with higher resolution, shorter scan times, and accurate quantification of physiological parameters.

4D-flow MRI↗
NASA NTRS1991

Effects of inlet distortion on the development of secondary flows in a subsonic axial inlet compressor rotor

Detailed flow measurements were taken inside an isolated axial compressor rotor operating subsonically near peak efficiency. Laser anemometer measurements were made with two inlet velocity profiles. One profile consisted of an unmodified baseline flow, and the second profile was distorted by placing axisymmetric screens on the hub and shroud well upstream of the rotor. A primary flow is defined in the rotor and deviations from this primary flow for each inlet flow condition identified. A comparison between the two flow deviations is made to assess the development of a passage vortex due to the distortion of the inlet flow. A comparison of experimental results with computational predictions from a Navier-Stokes solver showed good agreement between predicted and measured flow. Measured results indicate that a distorted inlet profile has minimal effect on the development of the flow in the rotor passage and the resulting passage vortex.

Owen, Albert K.↗
NASA NTRS1993

Direct simulation of compressible wall-bounded turbulence

When analyzing many turbulent flows, the effects of compressibility can be neglected. Even some relatively high-speed flows, such as boundary layers generated by a supersonic aircraft, produce turbulent statistics that are similar to those found for the incompressible case. There are other situations, however, in which the non-zero divergence of the turbulence leads to behavior that is fundamentally different from that found at constant density. Examples include flows created by internal combustion engines, hypersonic flight, and supersonic combustion. It is with instances such as these that this project is concerned. In particular, we are interested in the effects of compressibility on turbulence near a smooth solid constant temperature surface; our primary objective is an increased physical understanding that can be used to improve turbulence models of wall-bounded compressible flows. With this in mind, we have begun a direct numerical simulation (DNS) study of turbulence in a plane channel. Because all of the relevant spatial and temporal scales are to be resolved, the simulations require no subgrid scale parameterization. The DNS code developed by Buell to study compressible plane Couette flow was modified to solve the compressible Navier-Stokes equations in the plane channel. The channel was chosen over the Couette flow for two reasons: (1) to avoid using the very large streamwise domains needed to adequately capture the large Couette vortical structures; and (2) to make use of previous experience by considering the compressible version of a well established case and isolate finite Mach number effects by comparing to the incompressible channel.

Coleman, Gary N.↗
NASA NTRS2003

Asymmetric Stokes-V Profiles at the Penumbral Boundary of a Sunspot

We present the spectropolarimetric measurements of a sunspot in the active region NOAA 6958 (15S03W), situated near the central meridian disk passage. The follower polarity sunspot was somewhat symmetrically round shaped with an elongated penumbra. There were several opposite polarity magnetic elements at, and beyond the penumbral boundary. The H-alpha images of the sunspot show the bright emission regions near the penumbral boundary towards the sun-center, which was of opposite polarity with respect to the main spot. The net-circular polarization (NCP) map shows that NCP is negative in the inner part of the spot and positive at the penumbral boundary and near the H-alpha plage. The Doppler velocities were determined by measuring the center-of-gravity (COG) of the Stokes-I profile and zero-crossing (ZC) wavelength of the Stokes-V profiles. The COG velocity map in general agrees with the Evershed flow. In addition, it shows the up flow in the penumbral region. The ZC velocities show the strong down flow at the penumbral boundary. Double-lobed Stokes-V profiles are observed at the locations, where the penumbral fibrils terminate coinciding the H-alpha plage. The Double lobed profiles had an unshifted component similar to the Stokes-V profiles of the sunspot penumbra and a shifted component with a velocity of about 5 km/s. The amplitude of the second component increases along the penumbral fibril as a function of the distance from the center of the sunspot. In this paper we discuss the role of emerging flux in generating the observed double lobed profiles. Based on our present observations, we propose to observe with the Solar-B Spectropolarimeter for understanding the nature of emerging flux near the sunspots.

Choudhary, Debi Prasad↗
NASA NTRS1993

A computational study of wingtip vortex flowfield

The near-field behavior of a wingtip vortex flow is studied computationally using an incompressible flow solver for the Navier-Stokes equations based on the artificial compressibility method. Inaccuracies in current computational studies are addressed, especially, the role of numerical errors and transition/turbulence modeling. A subset problem is devised in order to make the study of vortex preservation more tractable. As part of the numerical checks, the flow is first run laminar while performing a systematic grid refinement study for the subset problem. Further studies on the numerical errors are conducted with the measured Reynolds stresses introduced into the momentum equations as source terms. As a preliminary study of turbulent flows, the one-equation Baldwin-Barth turbulence model is implemented as well as the approximation of the production term. The full-geometry case is computed using 1.1 million grid points. The results are compared with experiment.

Dacles-Mariani, Jennifer↗
DOE OSTI2022

Modeling of Cube Array Roughness: RANS, Large Eddy Simulation, and Direct Numerical Simulation

Abstract Flow over arrays of cubes is an extensively studied model problem for rough wall turbulent boundary layers. While considerable research has been performed in computationally investigating these topologies using direct numerical simulation (DNS) and large eddy simulation (LES), the ability of sublayer-resolved Reynolds-averaged Navier–Stokes (RANS) to predict the bulk flow phenomena of these systems is relatively unexplored, especially at low and high packing densities. Here, RANS simulations are conducted on six different packing densities of cubes in aligned and staggered configurations. The packing densities investigated span from what would classically be defined as isolated, up to those in the d-type roughness regime, filling in the gap in the present literature. Three different sublayer-resolved turbulence closure models were tested for each case: a low Reynolds number k–ϵ model, the Menter k–ω SST model, and a full Reynolds stress model. Comparisons of the velocity fields, secondary flow features, and drag coefficients are made between the RANS results and existing LES and DNS results. There is a significant degree of variability in the performance of the various RANS models across all comparison metrics. However, the Reynolds stress model demonstrated the best accuracy in terms of the mean velocity profile as well as drag partition across the range of packing densities.

Engineering↗
NASA NTRS2014

Assessing Fan Flutter Stability in the Presence of Inlet Distortion Using One-way and Two-way Coupled Methods

Concerns regarding noise, propulsive efficiency, and fuel burn are inspiring aircraft designs wherein the propulsive turbomachines are partially (or fully)embedded within the airframe; such designs present serious concerns with regard to aerodynamic and aeromechanic performance of the compression system in response to inlet distortion. Previously, a preliminary design of a forward-swept high-speed fan exhibited flutter concerns in clean-inlet flows, and the present author then studied this fan further in the presence of off-design distorted in-flows. A three-dimensional, unsteady, Navier-Stokes computational fluid dynamics code is applied to analyze and corroborate fan performance with clean inlet flow. This code, already validated in its application to assess aerodynamic damping of vibrating blades at various flow conditions using a loosely-coupled approach, is modified to include a tightly-coupled aeroelastic simulation capability, and then loosely-coupled and tightly-coupled methods arecompared in their evaluation of flutter stability in distorted in-flows.

fluid-structure interaction↗