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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Aerodynamic optimization studies on advanced architecture computers

The approach to carrying out multi-discipline aerospace design studies in the future, especially in massively parallel computing environments, comprises of choosing (1) suitable solvers to compute solutions to equations characterizing a discipline, and (2) efficient optimization methods. In addition, for aerodynamic optimization problems, (3) smart methodologies must be selected to modify the surface shape. In this research effort, a 'direct' optimization method is implemented on the Cray C-90 to improve aerodynamic design. It is coupled with an existing implicit Navier-Stokes solver, OVERFLOW, to compute flow solutions. The optimization method is chosen such that it can accomodate multi-discipline optimization in future computations. In the work , however, only single discipline aerodynamic optimization will be included.

Chawla, Kalpana↗
NASA NTRS1992

Unsteady spot heating of a drop in a microgravity environment

The unsteady localized spot heating of a liquid drop under zero-g conditions is examined theoretically. This pertains to space experiments to measure thermal properties of materials and the purpose here is to predict the thermal behavior of such systems. Spot heating can be achieved by a laser beam focused on a small region of the drop surface. The present theoretical model deals with situations of weak Marangoni flows, whereby the thermal transport is conduction dominated. The heat flow in the drop is treated as unsteady while the surrounding gaseous region is considered to be quasisteady. The ensuing thermally driven flow is analyzed in the Stokes regime.

Sadhal, Satwindar Singh↗
NASA NTRS1997

An Efficient Multiblock Method for Aerodynamic Analysis and Design on Distributed Memory Systems

The work presented in this paper describes the application of a multiblock gridding strategy to the solution of aerodynamic design optimization problems involving complex configurations. The design process is parallelized using the MPI (Message Passing Interface) Standard such that it can be efficiently run on a variety of distributed memory systems ranging from traditional parallel computers to networks of workstations. Substantial improvements to the parallel performance of the baseline method are presented, with particular attention to their impact on the scalability of the program as a function of the mesh size. Drag minimization calculations at a fixed coefficient of lift are presented for a business jet configuration that includes the wing, body, pylon, aft-mounted nacelle, and vertical and horizontal tails. An aerodynamic design optimization is performed with both the Euler and Reynolds Averaged Navier-Stokes (RANS) equations governing the flow solution and the results are compared. These sample calculations establish the feasibility of efficient aerodynamic optimization of complete aircraft configurations using the RANS equations as the flow model. There still exists, however, the need for detailed studies of the importance of a true viscous adjoint method which holds the promise of tackling the minimization of not only the wave and induced components of drag, but also the viscous drag.

Reuther, James↗
NASA NTRS1995

Navier-Stokes Aerodynamic Simulation of the V-22 Osprey on the Intel Paragon MPP

The paper will describe the Development of a general three-dimensional multiple grid zone Navier-Stokes flowfield simulation program (ENS3D-MPP) designed for efficient execution on the Intel Paragon Massively Parallel Processor (MPP) supercomputer, and the subsequent application of this method to the prediction of the viscous flowfield about the V-22 Osprey tiltrotor vehicle. The flowfield simulation code solves the thin Layer or full Navier-Stoke's equation - for viscous flow modeling, or the Euler equations for inviscid flow modeling on a structured multi-zone mesh. In the present paper only viscous simulations will be shown. The governing difference equations are solved using a time marching implicit approximate factorization method with either TVD upwind or central differencing used for the convective terms and central differencing used for the viscous diffusion terms. Steady state or Lime accurate solutions can be calculated. The present paper will focus on steady state applications, although time accurate solution analysis is the ultimate goal of this effort. Laminar viscosity is calculated using Sutherland's law and the Baldwin-Lomax two layer algebraic turbulence model is used to compute the eddy viscosity. The Simulation method uses an arbitrary block, curvilinear grid topology. An automatic grid adaption scheme is incorporated which concentrates grid points in high density gradient regions. A variety of user-specified boundary conditions are available. This paper will present the application of the scalable and superscalable versions to the steady state viscous flow analysis of the V-22 Osprey using a multiple zone global mesh. The mesh consists of a series of sheared cartesian grid blocks with polar grids embedded within to better simulate the wing tip mounted nacelle. MPP solutions will be shown in comparison to equivalent Cray C-90 results and also in comparison to experimental data. Discussions on meshing considerations, wall clock execution time, load balancing, and scalability will be provided.

Vadyak, Joseph↗
NASA NTRS1997

Static Aeroelastic Analysis of Transonic Wind Tunnel Models Using Finite Element Methods

A computational method for accurately predicting the static aeroelastic deformations of typical transonic transport wind tunnel models is described. The method utilizes a finite element method (FEM) for predicting the deformations. Extensive calibration/validation of this method was carried out using a novel wind-off wind tunnel model static loading experiment and wind-on optical wing twist measurements obtained during a recent wind tunnel test in the National Transonic Facility (NTF) at NASA LaRC. Further validations were carried out using a Navier-Stokes computational fluid dynamics (CFD) flow solver to calculate wing pressure distributions about several aeroelastically deformed wings and comparing these predictions with NTF experimental data. Results from this aeroelastic deformation method are in good overall agreement with experimentally measured values. Including the predicted deformations significantly improves the correlation between CFD predicted and experimentally measured wing & pressures.

Hooker, John R.↗
NASA NTRS2005

Receptivity of Supersonic Boundary Layers to Acoustic Disturbances

Boundary layer receptivity to two-dimensional slow and fast acoustic waves is investigated by solving Navier-Stokes equations for Mach 4.5 flow over a flat plate with a finite-thickness leading edge. Higher order spatial and temporal schemes are employed to obtain the solution whereby the flat-plate leading edge region is resolved by providing a sufficiently refined grid. The results show that the instability waves are generated in the leading edge region and that the boundary-layer is much more receptive to slow acoustic waves (by almost a factor of 20) as compared to the fast waves. Hence, this leading-edge receptivity mechanism is expected to be more relevant in the transition process for high Mach number flows. The effect of acoustic wave incidence angle is also studied and it is found that the receptivity of the boundary layer on the windward side (with respect to the acoustic forcing) decreases by more than a factor of 4 when the incidence angle is increased from 0 to 45 deg. However, the receptivity coefficient for the leeward side is found to vary relatively weakly with the incidence angle. The effect of leading-edge thickness is also studied and bluntness is found to stabilize the boundary layer. The relative significance of fast acoustic waves is enhanced in the presence of bluntness.

Malik, Mujeeb R.↗
NASA NTRS2006

Numerical Issues for Circulation Control Calculations

Steady-state and time-accurate two-dimensional solutions of the compressible Reynolds-averaged Navier- Stokes equations are obtained for flow over the Lockheed circulation control (CC) airfoil and the General Aviation CC (GACC) airfoil. Numerical issues in computing circulation control flows such as the effects of grid resolution, boundary and initial conditions, and unsteadiness are addressed. For the Lockheed CC airfoil computed solutions are compared with detailed experimental data, which include velocity and Reynolds stress profiles. Three turbulence models, having either one or two transport equations, are considered. Solutions are obtained on a sequence of meshes, with mesh refinement primarily concentrated on the airfoil circular trailing edge. Several effects related to mesh refinement are identified. For example, sometimes sufficient mesh resolution can exclude nonphysical solutions, which can occur in CC airfoil calculations. Also, sensitivities of the turbulence models with mesh refinement are discussed. In the case of the GACC airfoil the focus is on the difference between steady-state and time-accurate solutions. A specific objective is to determine if there is self-excited vortex shedding from the jet slot lip.

Swanson, Roy C., Jr.↗
NASA NTRS2007

Acoustic Receptivity of Mach 4.5 Boundary Layer with Leading- Edge Bluntness

Boundary layer receptivity to two-dimensional slow and fast acoustic waves is investigated by solving Navier-Stokes equations for Mach 4.5 flow over a flat plate with a finite-thickness leading edge. Higher order spatial and temporal schemes are employed to obtain the solution whereby the flat-plate leading edge region is resolved by providing a sufficiently refined grid. The results show that the instability waves are generated in the leading edge region and that the boundary-layer is much more receptive to slow acoustic waves (by almost a factor of 20) as compared to the fast waves. Hence, this leading-edge receptivity mechanism is expected to be more relevant in the transition process for high Mach number flows where second mode instability is dominant. Computations are performed to investigate the effect of leading-edge thickness and it is found that bluntness tends to stabilize the boundary layer. Furthermore, the relative significance of fast acoustic waves is enhanced in the presence of bluntness. The effect of acoustic wave incidence angle is also studied and it is found that the receptivity of the boundary layer on the windward side (with respect to the acoustic forcing) decreases by more than a factor of 4 when the incidence angle is increased from 0 to 45 deg. However, the receptivity coefficient for the leeward side is found to vary relatively weakly with the incidence angle.

Malik, Mujeeb R.↗
NASA NTRS2011

Receptivity of Hypersonic Boundary Layers to Acoustic and Vortical Disturbances

Boundary layer receptivity to two-dimensional acoustic disturbances at different incidence angles and to vortical disturbances is investigated by solving the Navier-Stokes equations for Mach 6 flow over a 7deg half-angle sharp-tipped wedge and a cone. Higher order spatial and temporal schemes are employed to obtain the solution. The results show that the instability waves are generated in the leading edge region and that the boundary layer is much more receptive to slow acoustic waves as compared to the fast waves. It is found that the receptivity of the boundary layer on the windward side (with respect to the acoustic forcing) decreases when the incidence angle is increased from 0 to 30 degrees. However, the receptivity coefficient for the leeward side is found to vary relatively weakly with the incidence angle. The maximum receptivity is obtained when the wave incident angle is about 20 degrees. Vortical disturbances also generate unstable second modes, however the receptivity coefficients are smaller than that for the acoustic waves. Vortical disturbances first generate the fast acoustic modes and they switch to the slow mode near the continuous spectrum.

Balakamar, P.↗
NASA NTRS2013

Assessments of a Turbulence Model Based on Menter's Modification to Rotta's Two-Equation Model

The main objective of this paper is to construct a turbulence model with a more reliable second equation simulating length scale. In the present paper, we assess the length scale equation based on Menter s modification to Rotta s two-equation model. Rotta shows that a reliable second equation can be formed in an exact transport equation from the turbulent length scale L and kinetic energy. Rotta s equation is well suited for a term-by-term modeling and shows some interesting features compared to other approaches. The most important difference is that the formulation leads to a natural inclusion of higher order velocity derivatives into the source terms of the scale equation, which has the potential to enhance the capability of Reynolds-averaged Navier-Stokes (RANS) to simulate unsteady flows. The model is implemented in the PAB3D solver with complete formulation, usage methodology, and validation examples to demonstrate its capabilities. The detailed studies include grid convergence. Near-wall and shear flows cases are documented and compared with experimental and Large Eddy Simulation (LES) data. The results from this formulation are as good or better than the well-known SST turbulence model and much better than k-epsilon results. Overall, the study provides useful insights into the model capability in predicting attached and separated flows.

Abdol-Hamid, Khaled S.↗
NASA NTRS2015

The Prediction of Scattered Broadband Shock-Associated Noise

A mathematical model is developed for the prediction of scattered broadband shock-associated noise. Model arguments are dependent on the vector Green's function of the linearized Euler equations, steady Reynolds-averaged Navier-Stokes solutions, and the two-point cross-correlation of the equivalent source. The equivalent source is dependent on steady Reynolds-averaged Navier-Stokes solutions of the jet flow, that capture the nozzle geometry and airframe surface. Contours of the time-averaged streamwise velocity component and turbulent kinetic energy are examined with varying airframe position relative to the nozzle exit. Propagation effects are incorporated by approximating the vector Green's function of the linearized Euler equations. This approximation involves the use of ray theory and an assumption that broadband shock-associated noise is relatively unaffected by the refraction of the jet shear layer. A non-dimensional parameter is proposed that quantifies the changes of the broadband shock-associated noise source with varying jet operating condition and airframe position. Scattered broadband shock-associated noise possesses a second set of broadband lobes that are due to the effect of scattering. Presented predictions demonstrate relatively good agreement compared to a wide variety of measurements.

Miller, Steven A. E.↗
NASA NTRS2016

Receptivity and Forced Response to Acoustic Disturbances in High-Speed Boundary Layers

Supersonic boundary-layer receptivity to freestream acoustic disturbances is investigated by solving the Navier-Stokes equations for Mach 3.5 flow over a sharp flat plate and a 7-deg half-angle cone. The freestream disturbances are generated from a wavy wall placed at the nozzle wall. The freestream acoustic disturbances radiated by the wavy wall are obtained by solving the linearized Euler equations. The results for the flat plate show that instability modes are generated at all the incident angles ranging from zero to highly oblique. However, the receptivity coefficient decreases by about 20 times when the incident angle increases from zero to a highly oblique angle of 68 degrees. The results for the cone show that no instability modes are generated when the acoustic disturbances impinge the cone obliquely. The results show that the perturbations generated inside the boundary layer by the acoustic disturbances are the response of the boundary layer to the external forcing. The amplitude of the forced disturbances inside the boundary layer are about 2.5 times larger than the incoming field for zero azimuthal wavenumber and they are about 1.5 times for large azimuthal wavenumbers.

Balakumar, P.↗
NASA NTRS2017

Receptivity and Forced Response to Acoustic Disturbances in High-Speed Boundary Layers

Supersonic boundary-layer receptivity to freestream acoustic disturbances is investigated by solving the Navier–Stokes equations for Mach 3.5 flow over a 7 deg half-angle cone. The freestream disturbances are generated from a wavy wall placed at the nozzle wall. The freestream acoustic disturbances radiated by the wavy wall are obtained by solving the linearized Euler equations. The results show that no noticeable instability modes are generated when the acoustic disturbances impinge the cone obliquely. The results show that the perturbations generated inside the boundary layer by the acoustic disturbances are the response of the boundary layer to the external forcing. The amplitude of the forced disturbances inside the boundary layer are about 2.5 times larger than the incoming field for zero azimuthal wave number, and they are about 1.5 times for large azimuthal wave numbers.

Balakumar, P.↗
NASA NTRS2000

Towards LES Models of Jets and Plumes

As pointed out by Rodi standard integral solutions for jets and plumes developed for discharge into infinite, quiescent ambient are difficult to extend to complex situations, particularly in the presence of boundaries such as the sea floor or ocean surface. In such cases the assumption of similarity breaks down and it is impossible to find a suitable entrainment coefficient. The models are also incapable of describing any but the most slowly varying unsteady motions. There is therefore a need for full time-dependent modeling of the flow field for which there are three main approaches: (1) Reynolds averaged numerical simulation (RANS), (2) large eddy simulation (LES), and (3) direct numerical simulation (DNS). Rodi applied RANS modeling to both jets and plumes with considerable success, the test being a match with experimental data for time-averaged velocity and temperature profiles as well as turbulent kinetic energy and rms axial turbulent velocity fluctuations. This model still relies on empirical constants, some eleven in the case of the buoyant jet, and so would not be applicable to a partly laminar plume, may have limited use in the presence of boundaries, and would also be unsuitable if one is after details of the unsteady component of the flow (the turbulent eddies). At the other end of the scale DNS modeling includes all motions down to the viscous scales. Boersma et al. have built such a model for the non-buoyant case which also compares well with measured data for mean and turbulent velocity components. The model demonstrates its versatility by application to a laminar flow case. As its name implies, DNS directly models the Navier-Stokes equations without recourse to subgrid modeling so for flows with a broad spectrum of motions (high Re) the cost can be prohibitive - the number of required grid points scaling with Re(exp 9/4) and the number of time steps with Re(exp 3/4). The middle road is provided by LES whereby the Navier-Stokes equations are formally filtered with the filter chosen to only exclude the smallest turbulent motions. If successful, LES should provide much of the detail available to DNS but at more bearable cost. Fatica et al. in comparing LES with DNS for a low Reynolds number jet showed that the LES could simulate the temporally evolving behavior including growth of the jet thickness. It is the intention of this report to explore the application of an LES model to jets and plumes. As always, before tackling complex situations, the model must be tested for the simplest of cases and so we address only two, a non-buoyant axisymmetric jet issuing steadily from an orifice into a semi-infinite stationary environment and a buoyant jet in the same environment. The work is a continuation of Basu and Mansour.

Webb, A. T.↗
NASA NTRS1986

Computation of separated flow on a ramp using the space marching conservative supra-characteristics method

Steady, hypersonic viscous flows over compression corners with streamwise separation have been computed using the space marching Conservative Supra-Characteristics Method (CSCM-S) of Lombard. The CSCM-S method permits stable space marching of the parabolized Navier-Stokes (PNS) equations through large separated flow regions. The present method has been used to compute surface pressure, heat transfer, and skin friction coefficients for two compression corner cases studied experimentally by Holden and Moselle. The computed results compare favorably with previous Navier-Stokes results and the experimental data. The present method has also been compared with the conventional Beam-Warming scheme for solving the PNS equations. Comparison are made for accuracy, computer time, and computer storage.

Stookesberry, D. C.↗
NASA NTRS1987

A numerical study of three-dimensional separated flow past a hemisphere cylinder

Separated and vortical flow about a hemisphere-cylinder body has been investigated. An algorithm featuring two implicit factors, and partial flux splitting has been used to solve the thin-layer Navier-Stokes equations. In analyzing the complex flow patterns, experimental data and topological concepts are used to complement the numerical results in interpreting the surface-flow patterns as well as the flowfield structures. Basic issues concerning the three-dimensional separation characteristics and the leeward vortical structures are examined.

Ying, Susan X.↗
NASA NTRS1976

Further studies of unsteady boundary layers with flow reversal

One set of calculations was performed using the first order, time dependent turbulent boundary layer equations, and extended earlier work by Nash and Patel to a wider range of flows. Another set of calculations was performed for laminar flow using the time dependent Navier-Stokes equations. The results of the calculations confirm previous conclusions concerning the existence of a regime of unseparated flow, containing an embedded region of reversal, which is accessible to first order boundary layer theory. However, certain doubts are cast on the precise nature of the events which accompany the eventual breakdown of the theory due to singularity onset. The earlier view that the singularity appears as the final event in a sequence involving rapid thickening of the boundary layer and the formation of a localized region of steep gradients is called into question by the present results. It appears that singularity onset is not necessarily preceded by rapid boundary layer thickening, or even necessarily produces immediate thickening.

Nash, J. F.↗
NASA NTRS1982

Burgers approximation for two-dimensional flow past an ellipse

A linearization of the Navier-Stokes equation due to Burgers in which vorticity is transported by the velocity field corresponding to continuous potential flow is examined. The governing equations are solved exactly for the two dimensional steady flow past an ellipse of arbitrary aspect ratio. The requirement of no slip along the surface of the ellipse results in an infinite algebraic system of linear equations for coefficients appearing in the solution. The system is truncated at a point which gives reliable results for Reynolds numbers R in the range 0 R 5. Predictions of the Burgers approximation regarding separation, drag and boundary layer behavior are investigated. In particular, Burgers linearization gives drag coefficients which are closer to observed experimental values than those obtained from Oseen's approximation. In the special case of flow past a circular cylinder, Burgers approximation predicts a boundary layer whose thickness is roughly proportional to R-1/2.

Dorrepaal, J. M.↗