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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Calculations of hot gas ingestion for a STOVL aircraft model

Hot gas ingestion problems for Short Take-Off, Vertical Landing (STOVL) aircraft are typically approached with empirical methods and experience. In this study, the hot gas environment around a STOVL aircraft was modeled as multiple jets in crossflow with inlet suction. The flow field was calculated with a Navier-Stokes, Reynolds-averaged, turbulent, 3D computational fluid dynamics code using a multigrid technique. A simple model of a STOVL aircraft with four choked jets at 1000 K was studied at various heights, headwind speeds, and thrust splay angles in a modest parametric study. Scientific visualization of the computed flow field shows a pair of vortices in front of the inlet. This and other qualitative aspects of the flow field agree well with experimental data.

Fricker, David M.↗
NASA NTRS1993

An adaptive grid/Navier-Stokes methodology for the calculation of nozzle afterbody base flows with a supersonic freestream

This paper describes an adaptive grid method for base flows in a supersonic freestream. The method is based on the direct finite-difference statement of the equidistribution principle. The weighting factor is a combination of the Mach number, density, and velocity first-derivative gradients in the radial direction. Two key ideas of the method are to smooth the weighting factor by using a type of implicit smoothing and to allow boundary points to move in the grid adaptation process. An AGARD nozzle afterbody base flow configuration is used to demonstrate the performance of the adaptive grid methodology. Computed base pressures are compared to experimental data. The adapted grid solutions offer a dramatic improvement in base pressure prediction compared to solutions computed on a nonadapted grid. A total-variation-diminishing (TVD) Navier-Stokes scheme is used to solve the governing flow equations.

Williams, Morgan↗
NASA NTRS2002

Decoupling the Role of Particle Inertia and Gravity on Particle Dispersion

Particle dispersion and the influence that particle momentum exchange has on the properties of a turbulent carrier flow in micro-gravity environments challenge present understanding and predictive schemes. The objective of this effort has been to develop and assess high-fidelity simulation tools for predicting particle transport within micro-gravity environments suspended in turbulent flows. The computational technique is based on Direct Numerical Simulation (DNS) of the incompressible Navier-Stokes equations. The particular focus of the present work is on the class of dilute flows in which particle volume fractions and inter-particle collisions are negligible. Particle motion is assumed to be governed by drag with particle relaxation times ranging from the Kolmogorov scale to the Eulerian timescale of the turbulence and particle mass loadings up to one. The velocity field was made statistically stationary by forcing the low wavenumbers of the flow. The calculations were performed using 96(exp 3) collocation points and the Taylor-scale Reynolds number for the stationary flow was 62. The effect of particles on the turbulence was included in the Navier-Stokes equations using the point-force approximation in which 96(exp 3) particles were used in the calculations. DNS results show that particles increasingly dissipate fluid kinetic energy with increased loading, with the reduction in kinetic energy being relatively independent of the particle relaxation time. Viscous dissipation in the fluid decreases with increased loading and is larger for particles with smaller relaxation times. Fluid energy spectra show that there is a non-uniform distortion of the turbulence with a relative increase in small-scale energy. The non-uniform distortion significantly affects the transport of the dissipation rate, with the production and destruction of dissipation exhibiting completely different behaviors. The spectrum of the fluid-particle energy exchange rate shows that the fluid drags particles at low wavenumbers while the converse is true at high wavenumbers for small particles. A spectral analysis shows that the increase of the high wavenumber portion of the fluid energy spectrum can be attributed to transfer of the fluid-particle covariance by the fluid turbulence. This in turn explains the relative increase of small-scale energy caused by small particles observed in the present simulations as well as those of others.

Squires, Kyle D.↗
NASA NTRS1991

Unsteady vortex-dominated flow around wings with oscillating leading-edge flaps

Unsteady vortex-dominated flow around delta wings with oscillating leading-edge flaps represents an important classs of problems for supermaneuverability and flow control of advanced aircraft. The problem is solved using time accurate integration of the unsteady, compressible, thin-layer Navier-Stokes equations in conjunction with the unsteady, linearized, Navier-displacement equations. Starting with an initial configuration of the wing and its flaps, the Navier-Stokes equations are solved on an initial structured grid for the steady flow. The forced oscillation of the flaps is then applied, and the problem is solved accurately in time. The Navier-displacement equations are solved for the grid deformation and the Navier-Stokes equations are solved for the flowfield. Symmetric and anti-symmetric flaps oscillations are presented to study the effect of the flaps oscillation on the leading-edge vortical flow.

Kandil, Osama A.↗
NASA NTRS1992

Numerical modeling of fluid and electromagnetic phenomena in an arcjet

An explicit numerical technique is used to solve the axisymmetric reduced electromagnetic field equation. The effect of an electrical arc on a viscous, axisymmetric flow is approximated using an implicit thin layer Navier-Stokes solver with additional electromagnetic source terms in conjunction with the explicit finite difference code.

Flowe, Anita C.↗
NASA NTRS2011

Multidisciplinary Modeling Software for Analysis, Design, and Optimization of HRRLS Vehicles

The concept for Highly Reliable Reusable Launch Systems (HRRLS) under the NASA Hypersonics project is a two-stage-to-orbit, horizontal-take-off / horizontal-landing, (HTHL) architecture with an air-breathing first stage. The first stage vehicle is a slender body with an air-breathing propulsion system that is highly integrated with the airframe. The light weight slender body will deflect significantly during flight. This global deflection affects the flow over the vehicle and into the engine and thus the loads and moments on the vehicle. High-fidelity multi-disciplinary analyses that accounts for these fluid-structures-thermal interactions are required to accurately predict the vehicle loads and resultant response. These predictions of vehicle response to multi physics loads, calculated with fluid-structural-thermal interaction, are required in order to optimize the vehicle design over its full operating range. This contract with ResearchSouth addresses one of the primary objectives of the Vehicle Technology Integration (VTI) discipline: the development of high-fidelity multi-disciplinary analysis and optimization methods and tools for HRRLS vehicles. The primary goal of this effort is the development of an integrated software system that can be used for full-vehicle optimization. This goal was accomplished by: 1) integrating the master code, FEMAP, into the multidiscipline software network to direct the coupling to assure accurate fluid-structure-thermal interaction solutions; 2) loosely-coupling the Euler flow solver FEFLO to the available and proven aeroelasticity and large deformation (FEAP) code; 3) providing a coupled Euler-boundary layer capability for rapid viscous flow simulation; 4) developing and implementing improved Euler/RANS algorithms into the FEFLO CFD code to provide accurate shock capturing, skin friction, and heat-transfer predictions for HRRLS vehicles in hypersonic flow, 5) performing a Reynolds-averaged Navier-Stokes computation on an HRRLS configuration; 6) integrating the RANS solver with the FEAP code for coupled fluid-structure-thermal capability; and 7) integrating the existing NASA SRGULL propulsion flow path prediction software with the FEFLO software for quasi-3D propulsion flow path predictions, 8) improving and integrating into the network, an existing adjoint-based design optimization code.

Spradley, Lawrence W.↗
NASA NTRS1985

Numerical study of Reynolds stress in compressible flows

A second order closure has been implemented in an implicit Navier-Stokes solver to study the behavior of the Reynolds stresses under the influence of severe pressure gradients. In the boundary layer zone, the strongly sheared character of the mean flow dominates the turbulence generation mechanisms. However, the pressure gradients also play a very important role for these processes, but at different locations within the boundary layer.

Vandromme, D.↗
NASA NTRS1996

Challenges and Accomplishments of Viscous Incompressible Flow Computations in Aerospace and Beyond

Recent progress in incompressible Navier-Stokes solution methods will be presented. Discussions are focused on the methods designed for complex geometry applications in three dimensions, and thus are limited to primitive variable formulation. Both steady- and unsteady-solution algorithms and their salient features are discussed. A summary of our recent progress in flow solver development is given followed by numerical studies of a few example problems of our current interest. Solvers discussed here are based on structured-grid approach using finite-difference or finite-volume frame work. This short course will be delivered in three one-our lectures. The material in the course are collected from the work performed by the Incompressible Navier-Stokes group at NASA Ames Research Center over the past several years, and can be found in our publications widely disseminated in the US and abroad. This short course is sponsored by AGARD Consultant and Exchange Program under Support Project P-110.

Kwak, Dochan↗
NASA NTRS1985

Unsteady turbulent boundary layers - Some effects of abrupt free-stream velocity changes

The response of a turbulent boundary layer to abrupt changes in free-stream velocity is studied. Streamwise velocity measurements are examined during the abrupt events and after them, as the flow relaxes to its equilibrium condition. The measured velocity field is found to be in good agreement with the analytic solution of the Stokes equation for boundary conditions describing abrupt changes. While the flow recovers to its equilibrium state, changes in integral, boundary-layer parameters scale on the local free-stream time of flight.

Brereton, G. J.↗
NASA NTRS1998

Numerical Simulations of Steady and Unsteady Transition in Low-Pressure Turbine Blade Rows

Transition plays an important role in the prediction of losses and performance in low-pressure turbines. The transition location on a turbine blade may vary significantly because of the wakes from upstream blade rows, and intermittent flow separation can also affect the transition process in an unsteady flow environment. In the present investigation, an unsteady Navier-Stokes analysis is used to predict transition in a low-pressure turbine cascade and a low-pressure turbine stage. The numerical flow analysis is third-order spatially accurate and second-order temporally accurate, and the equations of motion are integrated using an implicit time-marching procedure. The Baldwin-Lomax and k-epsilon turbulence models, in conjunction with several algebraic transition models, have been used to predict the location of transition. Predicted results include unsteady blade loadings, time-histories of the pressure, transition locations and boundary layer quantities, as well as performance quantities and comparisons with the available experimental/design data.

Dorney, Daniel↗
NASA NTRS2011

Steady and Unsteady Flow Effects of Circumferential Grooves Casing Treatment in a Transonic Compressor Rotor

The current paper reports on an investigation of steady and unsteady flow effects of circumferential grooves casing treatment in a transonic compressor rotor. Circumferential grooves casing treatment is used mainly to increase stall margin in axial compressors with a relatively small decrease in aerodynamic efficiency. It is widely believed that flow mechanisms of circumferential grooves casing treatment near stall conditions are not yet well understood even though this treatment has been used widely in real engines. Numerical analysis based on steady Reynolds-averaged Navier-Stokes (RANS) has been the primary tool used to understand flow mechanism for circumferential grooves casing treatment. Although steady RANS explains some flow effects of circumferential grooves casing treatment, it does not calculate all the measured changes in the compressor characteristics. Therefore, design optimization of circumferential grooves with steady RANS has not been very successful. As a compressor operates toward the stall condition, the flow field becomes transient. Major sources of self-generated flow unsteadiness are shock oscillation and interaction between the passage shock and the tip leakage vortex. In the present paper, an unsteady Reynolds-averaged Navier-Stokes (URANS) approach is applied to study the effects of circumferential grooves in a transonic compressor. The results from URANS are compared with the results from RANS and measured data. The current investigation shows that there are significant unsteady flow effects on the performance of the circumferential grooves casing treatment. For the currently investigated rotor, the unsteady effects are of the same magnitude as the steady effects in terms of extending the compressor stall margin.

Hah, Chunill↗
NASA NTRS1993

Prediction of incidence and surface roughness effects on turbine performance

The results of a Navier-Stokes analysis for predicting the change in turbine efficiency due to a change in either incidence or surface roughness is discussed. It was experimentally determined by Boynton, Tabibzadeh, and Hudson that polishing the SSME high pressure fuel turbine blades improved turbine efficiency by about 2 points over a wide range of operating conditions. These conditions encompassed the range of incidence seen by the turbine blading during flight. It is also necessary to be able to predict turbine performance at various operating points for future rocket turbopump applications. The code RVCQ3D, developed by Rod Chima, was used to determine the effects of changes in incidence angle on turbine blade row efficiency. The midspan Navier-Stokes results were used in conjunction with an inviscid flow analysis code to predict the efficiency of the two stage SSME over a wide range of operating conditions for smooth and rough turbine blades. The use of the Navier-Stokes analysis to predict changes in turbine efficiency due to variation in incidence angles was found to be superior to other incidence loss correlations available in the literature. The sensitivity of the Navier-Stokes results to grid parameters is discussed. The effects of the surface roughness were accounted for using the Cebeci-Chang rough wall turbulence model. This model was implemented in the code RVCQ3D. The implementation of this model for predicting the change in efficiency is also discussed.

Boyle, R. J.↗
DOE OSTI2023

Moment Representation of Regularized Lattice Boltzmann Methods on NVIDIA and AMD GPUs

The lattice Boltzmann method is a highly scalable Navier-Stokes solver that has been applied to flow problems in a wide array of domains. However, the method is bandwidth-bound on modern GPU accelerators and has a large memory footprint. In this paper, we present new 2D and 3D GPU implementations of two different regularized lattice Boltzmann methods, which are not only able to achieve an acceleration of ∼ 1.4 × w.r.t. reference lattice Boltzmann implementations but also reduce the memory requirements by up to 35% and 47% in 2D and 3D simulations respectively. These new approaches are evaluated on NVIDIA and AMD GPU architectures.

Valero Lara, Pedro↗
DOE OSTI2020

Two dimensional topology optimization of heat exchangers with the volume fraction method

We perform a comparison study of two topology optimizations methods applied to the design of a two fluids heat exchanger modeled with a coupled thermal-flow problem. The flow follows an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation without internal generation and high Peclet number. To keep the two fluid phases separated, we solve two Stokes-Brinkman equations, where the Brinkman term models the other phase as a solid. These two velocity fields are then fed to the heat transfer equation. Our goal is to maximize the enthalpy at the cold outlet while constraining the pressure drop. We first solve the design modeling the solid and fluid phases with a volume fraction variable. A SIMP-like penalization in the Brinkman term drives the optimization to a discrete design. The cost and constraint function derivatives are calculated with the library pyadjoint and the optimization is performed by IPOPT. We present optimized designs in two dimensions and discuss the influence of the parameters.

Beck, VictorA.↗
DOE OSTI2021

Shear thinning rheology in a concentrated suspension of fibers: the role of attractive interactions

Fiber-reinforced composites are ubiquitously encountered in the engineering, automobile, and aerospace industries. Fabrication of these composites requires the mixing of fibers dispersed in a liquid, and the final structure is affected by fiber properties, interactions, and flow fields. A good understanding of the rheology of fiber suspensions can aid in the design and optimization of the fabrication processes. To this end, we simulate the simple shear flow of fiber suspensions by accounting for short range lubrication, van der Waals attractive, electrostatic repulsion interactions, as well as friction between fibers. Direct numerical simulations are performed using the Immersed Boundary Method where the fibers are modeled as flexible slender bodies governed by the Euler-Bernoulli beam theory and the fluid flow is resolved using the Navier-Stokes equations. The simulation results for the suspension viscosity and yield stress are consistent with the experimental data from the literature for polyamide fiber suspensions. The shear-thinning behavior becomes stronger as we increase the magnitude of attractive interactions. Furthermore, we perform a parametric study varying fiber flexibility, volume fraction, aspect ratio, inertia and examine suspension viscosity and yield stress.

Khan, Monsurul↗
NASA NTRS1975

Numerical solutions of the unsteady Navier-Stokes equations for arbitrary bodies using boundary-fitted curvilinear coordinates

A method of automatic body-fitted curvilinear coordinate generation is described and used to construct a finite-difference solution of the full incompressible time-dependent Navier-Stokes equations for the unsteady laminar viscous flow arbitrary two-dimensional airfoils or any other two-dimensional body. A method of controlling the spacing of the coordinate lines encircling the body is developed in order to treat high Reynolds number flows, since the coordinate lines must concentrate near the surface to a greater degree as the Reynolds number increases. Multiple airfoils and submerged hydrofoils are treated as illustrative examples. The solution shows good agreement with the Blasius boundary layer solution for the flow past a semi-infinite flat plate.

Thompson, J. F.↗
NASA NTRS1978

Prediction of unsteady separated flows on oscillating airfoils

Techniques for calculating high Reynolds number flow around an airfoil undergoing dynamic stall are reviewed. Emphasis is placed on predicting the values of lift, drag, and pitching moments. Methods discussed include: the discrete potential vortex method; thin boundary layer method; strong interaction between inviscid and viscous flows; and solutions to the Navier-Stokes equations. Empirical methods for estimating unsteady airloads on oscillating airfoils are also described. These methods correlate force and moment data from wind tunnel tests to indicate the effects of various parameters, such as airfoil shape, Mach number, amplitude and frequency of sinosoidal oscillations, mean angle, and type of motion.

Mccroskey, W. J.↗
NASA NTRS1977

The finite element method in computational fluid mechanics

A finite element solution algorithm is established for a general statement of the time-averaged Navier-Stokes equations governing multi-dimensional turbulent flows. Numerical results are presented which evaluate factors affecting solution accuracy for a broad spectrum of linear and non-linear problem classes.

Baker, A. J.↗