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,603 records · Page 89

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 cleaninlet flows, and the present author then studied this fan further in the presence of off-design distorted in-flows. Continuing this research, a three-dimensional, unsteady, Navier-Stokes computational fluid dynamics code is again applied to analyze and corroborate fan performance with clean inlet flow and now with a simplified, sinusoidal distortion of total pressure at the aerodynamic interface plane. This code, already validated in its application to assess aerodynamic damping of vibrating blades at various flow conditions using a one-way coupled energy-exchange approach, is modified to include a two-way coupled time-marching aeroelastic simulation capability. The two coupling methods are compared in their evaluation of flutter stability in the presence of distorted in-flows.

aeroelasticity↗

Assessing Fan Flutter Stability in 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. Continuing this research, a three-dimensional, unsteady, Navier-Stokes computational fluid dynamics code is again applied to analyze and corroborate fan performance with clean inlet flow and now with a simplified, sinusoidal distortion of total pressure at the aerodynamic interface plane. This code, already validated in its application to assess aerodynamic damping of vibrating blades at various flow conditions using a one-way coupled energy-exchange approach, is modified to include a two-way coupled timemarching aeroelastic simulation capability. The two coupling methods are compared in their evaluation of flutter stability in the presence of distorted in-flows.

fluid-structure interaction↗

Shock-wave boundary layer interactions

Presented is a comprehensive, up-to-date review of the shock-wave boundary-layer interaction problem. A detailed physical description of the phenomena for transonic and supersonic speed regimes is given based on experimental observations, correlations, and theoretical concepts. Approaches for solving the problem are then reviewed in depth. Specifically, these include: global methods developed to predict sudden changes in boundary-layer properties; integral or finite-difference methods developed to predict the continuous evolution of a boundary-layer encountering a pressure field induced by a shock wave; coupling methods to predict entire flow fields; analytical methods such as multi-deck techniques; and finite-difference methods for solving the time-dependent Reynolds-averaged Navier-Stokes equations used to predict the development of entire flow fields. Examples are presented to illustrate the status of the various methods and some discussion is devoted to delineating their advantages and shortcomings. Reference citations for the wide variety of subject material are provided for readers interested in further study.

Delery, J.↗

Numerical solution of the Navier-Stokes equations about three-dimensional configurations: A survey

The numerical solution of the Navier-Stokes equations about three-dimensional configurations is reviewed. Formulational and computational requirements for the various Navier-Stokes approaches are examined for typical problems including the viscous flow field solution about a complete aerospace vehicle. Recent computed results, with experimental comparisons when available, are presented to highlight the presentation. The future of Navier-Stokes applications in three-dimensions is seen to be rapidly expanding across a broad front including internal and external flows, and flows across the entire speed regime from incompressible to hypersonic applications. Prospects for the future are described and recommendations for areas of concentrated research are indicated.

Holst, Terry L.↗

Transition to turbulence in an elliptic vortex

We study the three dimensional instability and nonlinear growth of the two dimensional flow described by the stream function Psi = (A sin b1 x sin b2 y)/(b1(exp 2) + b2(exp 2)) where b1 = pi/L1, b2 = pi/L2. This is a swirling flow in a box which is bounded by 0 less than x less than L1, 0 less than y less than L2 and is infinite in the z direction. This flow is a solution of the Navier-Stokes equation with A = exp(-v(b1(exp 2) + b2(exp 2))t) which slowly decays. We seek a viscous solution which starts near this one and slips along but does not penetrate the bounding walls. The vorticity of the basic flow is w(sub z) = A sin b1 x sin b2 y which has maximum value A at the center of the box and drops to zero at the boundaries. We can think of the resulting flow as that of a captive vortex.

Lundgren, Thomas S.↗

Numerical simulation of steady supersonic flow over spinning bodies of revolution

A recently reported parabolized Navier-Stokes code has been employed to compute the supersonic flowfield about a spinning cone and spinning and nonspinning ogive cylinder and boattailed bodies of revolution at moderate incidence. The computations were performed for flow conditions where extensive measurements for wall pressure, boundary-layer velocity profiles, and Magnus force had been obtained. Comparisons between the computational results and experiment indicate excellent agreement for angles of attack up to 6 deg. At angles greater than 6 deg discrepancies are noted which are tentatively attributed to turbulence modeling errors. The comparisons for Magnus effects show that the code accurately predicts the effects of body shape for the selected models.

Sturek, W. B.↗

Highly Accurate Experimental Heave Decay Tests with a Floating Sphere: A Public Benchmark Dataset for Model Validation of Fluid–Structure Interaction

Highly accurate and precise heave decay tests on a sphere with a diameter of 300 mm were completed in a meticulously designed test setup in the wave basin in the Ocean and Coastal Engineering Laboratory at Aalborg University, Denmark. The tests were dedicated to providing a rigorous benchmark dataset for numerical model validation. The sphere was ballasted to half submergence, thereby floating with the waterline at the equator when at rest in calm water. Heave decay tests were conducted, wherein the sphere was held stationary and dropped from three drop heights: a small drop height, which can be considered a linear case, a moderately nonlinear case, and a highly nonlinear case with a drop height from a position where the whole sphere was initially above the water. The precision of the heave decay time series was calculated from random and systematic standard uncertainties. At a 95% confidence level, uncertainties were found to be very low—on average only about 0.3% of the respective drop heights. Physical parameters of the test setup and associated uncertainties were quantified. A test case was formulated that closely represents the physical tests, enabling the reader to do his/her own numerical tests. The paper includes a comparison of the physical test results to the results from several independent numerical models based on linear potential flow, fully nonlinear potential flow, and the Reynolds-averaged Navier–Stokes (RANS) equations. A high correlation between physical and numerical test results is shown. The physical test results are very suitable for numerical model validation and are public as a benchmark dataset.

42 ENGINEERING↗

The low frequency oscillation in the flow over a NACA0012 airfoil with an iced leading edge

The unusually low frequency oscillation in the wake of an airfoil is explored experimentally as well as computationally for a NACA0012 airfoil with a glaze ice accretion at the leading edge. Experimentally, flow oscillations were observed at low frequencies that correspond to a Strouhal number of about 0.02. This occurred in the angle of attack range of 8 to 9 deg, near the onset of static stall for this airfoil. With a Navier-Stokes computation, limit-cycle oscillations in the flow and in the aerodynamic forces were also observed at low Strouhal numbers. However, the occurrence of the oscillation is found to depend on the turbulence model in use as well as the Reynolds number.

Zaman, K. B. M. Q.↗

The low frequency oscillation in the flow over a NACA0012 airfoil with an 'iced' leading edge

The unusually low frequency oscillation in the wake of an airfoil is explored experimentally as well as computationally for a NACA0012 airfoil with a glaze ice accretion at the leading edge. Experimentally, flow oscillations were observed at low frequencies that correspond to a Strouhal number of about 0.02. This occurred in the angle of attack range of 8 to 9 deg, near the onset of static stall for this airfoil. With a Navier-Stokes computation, limit-cycle oscillations in the flow and in the aerodynamic forces were also observed at low Strouhal numbers. However, the occurrence of the oscillation is found to depend on the turbulence model in use as well as the Reynolds number.

Zaman, K. B. M. Q.↗

An Isolated Circular Synthetic Jet in Cross-Flow at Low Momentum-Flux Ratio

A joint experimental and computational investigation was carried out for a round synthetic jet issuing normal to a turbulent boundary layer at a momentum-flux ratio of one. Distributions of velocity and turbulence intensity were measured by hot-wire anemometry. Numerical results were obtained using unsteady Reynolds-averaged Navier-Stokes (URANS) computations. Time and phase-averaged flow properties were compared on the cross sectional plane at x/D = 0.53, 5 and 10 as well as on the axial plane of symmetry. Overall, the numerical results agreed well with the experimental data. CFD predicted a somewhat larger velocity deficit in regions of low-momentum fluid pulled up from the boundary layer. Phase- averaged velocity contours at the plane of symmetry indicated good match between experiments and CFD regarding the size and the position of the periodic flow structure. However, some differences occurred in details such as the shape and inclination of the low-speed flow structure.

Milanovic, Ivana M.↗

Line relaxation methods for the solution of 2D and 3D compressible flows

An implicit finite element based algorithm for the compressible Navier-Stokes equations is outlined, and the solution of the resulting equation by a line relaxation on general meshes of triangles or tetrahedra is described. The problem of generating and adapting unstructured meshes for viscous flows is reexamined, and an approach for both 2D and 3D simulations is proposed. An efficient approach appears to be the use of an implicit/explicit procedure, with the implicit treatment being restricted to those regions of the mesh where viscous effects are known to be dominant. Numerical examples demonstrating the computational performance of the proposed techniques are given.

Hassan, O.↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper presents the results of an ongoing study into the linear stability characteristics of the boundary layer flow over the common research model with natural laminar flow (CRMNLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work involving parallel stability computations of a boundary layer flow based on the conical flow approximation has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien- Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. Here, we examine the effects of the simplifying approximations in both basic state computation and the stability analysis, with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, the basic states are computed by using full Navier-Stokes equations and the stability analysis is performed by using a nonorthogonal coordinate system that allows a clear distinction between planar TS and CF instabilities. Furthermore, the effects of curvature and nonparallel mean flow have been included in the stability computations based on the parabolized stability equations (PSE). The fully turbulent Reynolds-Averaged-Navier-Stokes (RANS) mean flow solutions show good agreement with the measured wall pressure distribution. Viscous-inviscid interactive effects are observed to be important because the shock fronts along the suction surface are influenced by the imposed transition front. The stability results confirm the previous findings related to TS amplification within the inboard region of the wing and the dominance of stationary CF modes in the outboard region. However, given the close proximity of the measured transition front and the dual shock system within the outer part of the wing, the onset of transition may well be shock limited within the outboard region. In general, the transition criterion based on the dual N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at Re MAC = 15 million and AoA = 1.5, 2 degrees; however, the low values of the correlating N-factors at Re MAC = 17.5 million support the hypothesis that the measured transition at the higher Reynolds number may have been strongly influenced by the merging of turbulent wedges that originate from surface imperfections near the leading edge.

Boundary layer transition↗

Improved Models for the Porous Surface with Passive Control

A computational investigation of afterbody flow using a passive control method is conducted. The passive control method consists of a porous surface placed over a plenum. The purpose of the passive control method is to exploit the adverse pressure gradient present in afterbody flow in an attempt to reduce boundary layer separation and afterbody drag. Four different porous wall models are used to model the transpiration velocity in the region of passive control. A three-dimensional, time-dependent, Reynolds-averaged, simplified Navier-Stokes solver, PAB3D, is used to simulate afterbody flow with and without passive control. Three afterbody configurations with boat-tail angles of 10, 20, and 30 deg. are used to obtain two-dimensional solutions with a freestream Mach number of 0.6 and nozzle pressure ratio of 6. The region of passive control was initially placed from 20-60% of the nozzle length. The effect of the porous placement and porous extent is also studied. Baseline (no porosity) two-dimensional solutions are qualitatively similar to experimental data but under-predict the magnitude of the pressure recovery. Results for the subsonic solutions show losses in the pressure recovery for some cases with passive control. Three-dimensional effects are also investigated and seen to be very significant. Three-dimensional baseline solutions, for both sub- and super-critical freestream Mach numbers, compare very favorably with the experimental data in comparison to the two-dimensional solution. Future work is required to examine three-dimensional afterbody flows with passive porosity.

Kraushaar, Sandra L.↗

Large Photospheric Doppler Shift in Solar Active Region 12673: I. Field-Aligned Flows

Delta (δ) sunspots sometimes host fast photospheric flows along the central magnetic polarity inversion line (PIL). Here we study the strong Doppler shift signature in the central penumbral light bridge of solar active region NOAA 12673. Observations from the Helioseismic and Magnetic Imager (HMI) indicate highly sheared, strong magnetic fields. Large Doppler shifts up to 3.2 km s −1 appeared during the formation of the light bridge and persisted for about 16 hours. A new velocity estimator, called DAVE4VMwDV, reveals fast converging and shearing motion along the PIL from HMI vector magnetograms, and recovers the observed Doppler signal much better than an old version of the algorithm. The inferred velocity vectors are largely (anti-)parallel to the inclined magnetic fields, suggesting that the observed Doppler shift contains significant contribution from the projected, field-aligned flows. High-resolution observations from the Hinode/Spectro-Polarimeter (SP) further exhibit a clear correlation between the Doppler velocity and the cosine of the magnetic inclination, which is in agreement with HMI results and consistent with a field-aligned flow of about 9.6 km s −1 . The complex Stokes profiles suggest significant gradients of physical variables along the line of sight. We discuss the implications on the δ-spot magnetic structure and the flow-driving mechanism.

Solar active region magnetic fields↗

On computations of the integrated space shuttle flowfield using overset grids

Numerical simulations using the thin-layer Navier-Stokes equations and chimera (overset) grid approach were carried out for flows around the integrated space shuttle vehicle over a range of Mach numbers. Body-conforming grids were used for all the component grids. Testcases include a three-component overset grid - the external tank (ET), the solid rocket booster (SRB) and the orbiter (ORB), and a five-component overset grid - the ET, SRB, ORB, forward and aft attach hardware, configurations. The results were compared with the wind tunnel and flight data. In addition, a Poisson solution procedure (a special case of the vorticity-velocity formulation) using primitive variables was developed to solve three-dimensional, irrotational, inviscid flows for single as well as overset grids. The solutions were validated by comparisons with other analytical or numerical solution, and/or experimental results for various geometries. The Poisson solution was also used as an initial guess for the thin-layer Navier-Stokes solution procedure to improve the efficiency of the numerical flow simulations. It was found that this approach resulted in roughly a 30 percent CPU time savings as compared with the procedure solving the thin-layer Navier-Stokes equations from a uniform free stream flowfield.

Chiu, I-T.↗

Two-dimensional computations of multi-stage compressor flows using a zonal approach

A clear understanding of the fluid dynamics associated with rotor/stator configurations can be very helpful when optimizing the performance of turbomachinery. In this study, a two-dimensional, implicit, thin-layer, Navier-Stokes zonal approach has been used to investigate the flow within a 2 1/2-stage compressor. Relative motion between the rotor and stator airfoils is made possible with the use of systems of patched and overlaid grids that move with respect to each other. The treatment of multistage turbomachines with arbitrary numbers of airfoils per row is made possible by the use of a flexible database system. Results in the form of instantaneous pressure and entropy contours and time-averaged pressures are presented for the 2 1/2-stage compressor. Time-averaged pressures and pressure amplitudes for a single-stage turbine configuration are also presented. The numerical results compare well with experimental data.

Gundy-Burlet, Karen L.↗

Critical issues in viscous flow computations

In developing computer programs to numerically solve the Navier-Stokes equations, the purpose of the computation must be clearly kept in mind. In the Air Force, the purpose is to provide design information on non-linear aerodynamic phenomenon for aircraft that perform throughout the flight corridor. This translates into the requirement for a computer program which can solve the time averaged compressible Navier-Stokes equations (with a turbulence model) in three dimensions for generalized geometries. The intended application of the results then controls the priorities in addressing critical issues. Recurrent problem areas encountered in the study of viscous flow include: (1) grid generation for arbitrary geometry; (2) numerical difficulties; (3) turbulence models; (4) accuracy and efficiency; and (5) smearing of discontinuities.

Hankey, W. L., Jr.↗

NCC: A Physics-Based Design and Analysis Tool for Combustion Systems

The National Combustion Code (NCC) is an integrated system of computer codes for physics-based design and analysis of combustion systems. It uses unstructured meshes and runs on parallel computing platforms. The NCC is composed of a set of distinct yet closely related modules. They are: (1) a gaseous flow module solving 3-D Navier-Stokes equations; (2) a turbulence module containing the non-linear k-epsilon models; (3) a chemistry module using either the conventional reduced kinetics approach of solving species equations or the Intrinsic Low Dimensional Manifold (ILDM) kinetics approach of table looking up in conjunction with solving the equations of the progressive variables; (4) a turbulence-chemistry interaction module including the option of solving the joint probability density function (PDF) for species and enthalpy; and (5) a spray module for solving the liquid phase equations. In early 1995, an industry-government team was formed to develop the NCC. In July 1998, the baseline beta version was completed and presented in two NCC sessions at the 34th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, July 1998. An overview of this baseline beta version was presented at the NASA HPCCP/CAS Workshop 98, August 1998. Since then, the effort has been focused on the streamlining, validation, and enhancement of the th baseline beta version. The progress is presented in two NCC sessions at the AIAA 38 Aerospace Sciences Meeting & Exhibit, January 2000. At this NASA HPCCP/CAS Workshop 2000, an overview of the NCC papers presented at the AIAA 38 th Aerospace Sciences Meeting & Exhibit is presented, with emphasis on the reduction of analysis time of simulating the (gaseous) reacting flows in full combustors. In addition, results of NCC simulation of a modern turbofan combustor will also be reported.

Liu, Nan-Suey↗