Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Vortex dynamics”

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 865 records · Page 48

Influence of disk leakage path on labyrinth seal inlet swirl ratio

The results of numerous investigators have shown the importance of labyrinth seal inlet swirl on the calculated dynamic stiffness of labyrinth seals. These results have not included any calculation of inlet leakage swirl as a function of geometry and sealing conditions of the given seal. This paper outlines a method of calculating the inlet swirl at a given seal by introducing a radial chamber to predict the gas swirl as it goes from the stage tip down to the seal location. For a centrifugal compressor, this amounts to including the flow path from the impeller discharge, down the back of the disk or front of the cover, then into the shaft seal or eye packing, respectively. The solution includes the friction factors of both the disk and stationary wall with account for mass flow rate and calculation of radial pressure gradients by a free vortex solution. The results of various configurations are discussed and comparisons made to other published results of disk swirl.

Kirk, R. Gordon↗

Development of a High-Order Space-Time Matrix-Free Adjoint Solver

The growth in computational power and algorithm development in the past few decades has granted the science and engineering community the ability to simulate flows over complex geometries, thus making Computational Fluid Dynamics (CFD) tools indispensable in analysis and design. Currently, one of the pacing items limiting the utility of CFD for general problems is the prediction of unsteady turbulent ows.1{3 Reynolds-averaged Navier-Stokes (RANS) methods, which predict a time-invariant mean flowfield, struggle to provide consistent predictions when encountering even mild separation, such as the side-of-body separation at a wing-body junction. NASA's Transformative Tools and Technologies project is developing both numerical methods and physical modeling approaches to improve the prediction of separated flows. A major focus of this e ort is efficient methods for resolving the unsteady fluctuations occurring in these flows to provide valuable engineering data of the time-accurate flow field for buffet analysis, vortex shedding, etc. This approach encompasses unsteady RANS (URANS), large-eddy simulations (LES), and hybrid LES-RANS approaches such as Detached Eddy Simulations (DES). These unsteady approaches are inherently more expensive than traditional engineering RANS approaches, hence every e ort to mitigate this cost must be leveraged. Arguably, the most cost-effective approach to improve the efficiency of unsteady methods is the optimal placement of the spatial and temporal degrees of freedom (DOF) using solution-adaptive methods.

Adjoint↗

An Assessment of CFD Effectiveness for Vortex Flow Simulation to Meet Preliminary Design Needs

The low-speed flight and transonic maneuvering characteristics of combat air vehicles designed for efficient supersonic flight are significantly affected by the presence of free vortices. At moderate-to-high angles of attack, the flow invariably separates from the leading edges of the swept slender wings, as well as from the forebodies of the air vehicles, and rolls up to form free vortices. The design of military vehicles is heavily driven by the need to simultaneously improve performance and affordability.1 In order to meet this need, increasing emphasis is being placed on using Modeling & Simulation environments employing the Integrated Product & Process Development (IPPD) concept. The primary focus is on expeditiously providing design teams with high-fidelity data needed to make more informed decisions in the preliminary design stage. Extensive aerodynamic data are needed to support combat air vehicle design. Force and moment data are used to evaluate performance and handling qualities; surface pressures provide inputs for structural design; and flow-field data facilitate system integration. Continuing advances in computational fluid dynamics (CFD) provide an attractive means of generating the desired data in a manner that is responsive to the needs of the preliminary design efforts. The responsiveness is readily characterized as timely delivery of quality data at low cost.

Raj, P.↗

Dynamic Mode Decomposition of Unsteady Pressure-Sensitive Paint Measurements for the NASA Unitary Plan Wind Tunnel Tests

This paper describes the Dynamic Mode Decomposition (DMD) of the pressures on the scale model of the Space Launch System (SLS) Block 1 cargo vehicle with the Unsteady Pressure-Sensitive Paint (uPSP) measurements, which were collected in the Ascent Transient Aerodynamics Tests with the Unitary Plan Wind Tunnel 11-by-11-foot Transonic Wind Tunnel in September 2019 at NASA Ames Research Center. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. The conventional DMD algorithm is based on the Singular Value Decomposition (SVD) of the data matrix. For the matrix of the uPSP measurements of the SLS ATAT, the number of rows is equal to the number of nodes in the grid of the scale model, and the number of columns is equal to the number of frames in the videos taken with 4 Phantom high-speed cameras. In this paper, it is verified that, for the time series with zero mean value, the DMD is equivalent to the decomposition with the Discrete Fourier Transform (DFT). Considering the uPSP is mainly used in the assessment of the unsteady, aerodynamic phenomena, the DMD of the uPSP measurements can be implemented in two steps: (1) subtract the mean value from the uPSP measurement on each of the grid nodes; (2) apply the Fast Fourier Transform (FFT) on the resulting zero-mean time series. The DMD of the uPSP measurements with FFT has two advantages: (1) the computational complexity of FFT is O(N*logN), where N is the length of the time series; (2) compared to the SVD-based DMD algorithm, the DMD with FFT can be easily implemented in parallel processing. A sample matrix of uPSP measurements, at the size of 341 grid nodes and 128 frames, is generated. Figures 1 and 2 show the eigenvalues and the ratios of the eigenvectors, respectively, of the sample matrix, without and with the mean value removed on each of the grid nodes, computed with the SVD-based DMD and the FFT. The figures demonstrate the equivalence of the SVD-based DMD and the decomposition with DFT/FFT for the time series with zero mean value. The results of DMD of the uPSP measurements of the SLS ATAT in September 2019 are presented in the paper. The DMD modes at different frequencies are shown, the aerodynamic phenomena (e.g. shockwave and vortex shedding) are demonstrated and the correlation of the DMD modes with the test configuration parameter (e.g., the Mach Number) is discussed. Figure 3 shows a software tool to visualize the DMD modes. The code to implement the algorithm described in this paper was written in C, with libraries of FFTW for FFT and MPI/OpenMP for parallel processing, and executed on the NASA Pleiades supercomputer. Funding for this research was provided by the NASA Aerosciences Evaluation and Test Capabilities Project.

Pressure-Sensitive Paint↗

Ionospheric traveling convection vortices observed near the polar cleft - A triggered response to sudden changes in the solar wind

Analysis of 20-second resolution magnetometer data from an array of temporary stations operated around Sondre Stromfjord, Greenland, during the summer of 1986 shows the signatures of localized ionospheric traveling convection vortices. An example of an isolated event of this kind observed near 08 local time is presented in detail. This event consists of a twin vortex pattern of convection consistent with the presence of two field-aligned current filaments separated by about 600 km in the east-west direction. This system of currents is observed to move westward (tailward) past the array of stations at about 4 km/sec. The event is associated with relative quiet time ionospheric convection and occurs during an interval of northward IMF. It is, however, associated with a large fluctuation in both the Z and Y components of the IMF and with a large sudden decrease in the solar wind number density. The propagation of the system is inconsistent with existing models of FTE current systems, but nevertheless appears to be related to a readjustment of the magnetopause boundary to a sudden change in the solar wind dynamic pressure and/or to a change in reconnection brought about by a sudden reorientation of the IMF.

Friis-Christensen, E.↗

Diagnostic Comparison of Meteorological Analyses during the 2002 Antarctic Winter

Several meteorological datasets, including U.K. Met Office (MetO), European Centre for Medium-Range Weather Forecasts (ECMWF), National Centers for Environmental Prediction (NCEP), and NASA's Goddard Earth Observation System (GEOS-4) analyses, are being used in studies of the 2002 Southern Hemisphere (SH) stratospheric winter and Antarctic major warming. Diagnostics are compared to assess how these studies may be affected by the meteorological data used. While the overall structure and evolution of temperatures, winds, and wave diagnostics in the different analyses provide a consistent picture of the large-scale dynamics of the SH 2002 winter, several significant differences may affect detailed studies. The NCEP-NCAR reanalysis (REAN) and NCEP-Department of Energy (DOE) reanalysis-2 (REAN-2) datasets are not recommended for detailed studies, especially those related to polar processing, because of lower-stratospheric temperature biases that result in underestimates of polar processing potential, and because their winds and wave diagnostics show increasing differences from other analyses between similar to 30 and 10 hPa (their top level). Southern Hemisphere polar stratospheric temperatures in the ECMWF 40-Yr Re-analysis (ERA-40) show unrealistic vertical structure, so this long-term reanalysis is also unsuited for quantitative studies. The NCEP/Climate Prediction Center (CPC) objective analyses give an inferior representation of the upper-stratospheric vortex. Polar vortex transport barriers are similar in all analyses, but there is large variation in the amount, patterns, and timing of mixing, even among the operational assimilated datasets (ECMWF, MetO, and GEOS-4). The higher-resolution GEOS-4 and ECMWF assimilations provide significantly better representation of filamentation and small-scale structure than the other analyses, even when fields gridded at reduced resolution are studied. The choice of which analysis to use is most critical for detailed transport studies (including polar process modeling) and studies involving synoptic evolution in the upper stratosphere. The operational assimilated datasets are better suited for most applications than the NCEP/CPC objective analyses and the reanalysis datasets.

stratosphere↗

TPSAS-NF1676L-23432-DND

The POLSTRACC (POlar STRAtosphere in a Changing Climate) field campaign was conducted in the Arctic during December 2015 - March 2016 to investigate the chemical, microphysical, and dynamical processes of the Arctic lowermost stratosphere and upper troposphere. The primary measurement platform for POLSTRACC was the German HALO (High Altitude LOng range) research aircraft carrying a large suite of in situ and remote sensing instruments to measure key chemical species, tracers, as well as aerosol and cloud particles and meteorological parameters. Two primary science objectives of POLSTRACC are to improve our understanding of polar stratospheric cloud (PSC) particle characteristics and formation processes and investigate the impact of Arctic cirrus clouds on radiative forcing and chlorine activation. To complement the more focused measurements from the POLSTRACC field campaign, we have used spaceborne lidar measurements from CALIPSO to characterize PSC occurrence and composition, as well as the occurrence of Arctic cirrus during the 2015-2016 season on vortex-wide scales. In this paper, we present a general overview of the 2015-2016 winter, examine in detail the evolution of PSCs and cirrus during the season, and explore the unique aspects of this season in attempt to understand the underlying physical mechanisms.

Michael C Pitts↗

Characterizing Aeroelastic Systems Using Eigenanalysis, Explicitly Retaining The Aerodynamic Degrees of Freedom

Discrete time aeroelastic models with explicitly retained aerodynamic modes have been generated employing a time marching vortex lattice aerodynamic model. This paper presents analytical results from eigenanalysis of these models. The potential of these models to calculate the behavior of modes that represent damped system motion (noncritical modes) in addition to the simple harmonic modes is explored. A typical section with only structural freedom in pitch is examined. The eigenvalues are examined and compared to experimental data. Issues regarding the convergence of the solution with regard to refining the aerodynamic discretization are investigated. Eigenvector behavior is examined; the eigenvector associated with a particular eigenvalue can be viewed as the set of modal participation factors for that particular mode. For the present formulation of the equations of motion, the vorticity for each aerodynamic element appears explicitly as an element of each eigenvector in addition to the structural dynamic generalized coordinates. Thus, modal participation of the aerodynamic degrees of freedom can be assessed in M addition to participation of structural degrees of freedom.

Heeg, Jennifer↗

Computational Fluid Dynamics Simulations of the Transonic Dynamics Tunnel Airstream Oscillator System

This paper presents computational fluid dynamics simulations of the flow in the NASA Langley Research Center Transonic Dynamics Tunnel (TDT) with the airstream oscillator system (AOS) in operation. Test section flow angle can be modulated using sinusoidally oscillating vanes which are located in pairs within the contraction section of the tunnel circuit. A compilation of previously published experimental data for validation has been reviewed that contains hot-wire anemometer and fast-response probe data acquired in an empty test section with the vanes oscillating in-phase and out-of-phase. These experimental data have not previously been compared with computational results. The present results provide an opportunity for validation of the simulations of the sinusoidal AOS flow field in the TDT test section. Computed tunnel centerline, lateral and cross-section contour data of flow angle and vorticity show that resolving the vane tip vortex is critical to accurately predicting the test section center flow angle. The computed results show, as do the experimental data, what appears to be a resonance at higher frequencies. It is currently not known what is the cause of that resonance. Further refinement of the mesh is expected to improve the comparisons with experimental data

Robert E Bartels↗

Investigation of Transonic Truss-Braced Wing Aircraft Transonic Wing-Strut Interference Effects Using FUN3D

This paper presents a computational study of transonic wing-strut interference effects of Transonic Truss-Braced Wing (TTBW) aircraft using the high-fidelity CFD (Computational Fluid Dynamics) code FUN3D (Fully Unstructured Navier-Stokes Three Dimensional). The study is conducted for the wing-strut and the wing-alone configurations at different Mach numbers and Reynolds numbers. The interference effects are calculated by comparing the wing aerodynamics along the spanwise direction between the wing-strut and the wing-alone configurations. The presence of the strut underneath the wing induces a suction peak on the lower surface of the wing, which causes changes in aerodynamic forces and moments, as well as the aerodynamic center location. The interference effects become more pronounced as the Mach number increases. The Reynolds number has less impact on the interference effects. A transonic wing-strut interference aerodynamic correction method is developed for use in a lower-fidelity tool, VSPAERO [a vortex lattice flow solver], coupled to a finite-element model for rapid flutter analysis.

Xiong, Juntao↗

Reduction of Helicopter BVI Noise, Vibration, and Power Consumption Through Individual Blade Control

A wind tunnel test was conducted with a full-scale BO 105 helicopter rotor to evaluate the potential of open-loop individual blade control (IBC) to improve rotor performance, to reduce blade vortex interaction (BVI) noise, and to alleviate helicopter vibrations. The wind tunnel test was an international collaborative effort between NASA/U.S. Army AFDD, ZF Luftfahrttechnik, Eurocopter Deutschland, and the German Aerospace Laboratory (DLR) and was conducted under the auspices of the U.S./German MOU on Rotorcraft Aeromechanics. In this test the normal blade pitch links of the rotor were replaced by servo-actuators so that the pitch of each blade could be controlled independently of the other blades. The specially designed servoactuators and IBC control system were designed and manufactured by ZF Luftfahrttechnik, GmbH. The wind tunnel test was conducted in the 40- by 80-Foot Wind Tunnel at the NASA Ames Research Center. An extensive amount of measurement information was acquired for each IBC data point. These data include rotor performance, static and dynamic hub forces and moments, rotor loads, control loads, inboard and outboard blade pitch motion, and BVI noise data. The data indicated very significant (80 percent) simultaneous reductions in both BVI noise and hub vibrations could be obtained using multi-harmonic input at the critical descent (terminal approach) condition. The data also showed that performance improvements of up to 7 percent could be obtained using 2P input at high-speed forward flight conditions.

Jacklin, Stephen A.↗

Flowfield dynamics in blunt fin-induced shock wave/turbulent boundary layer interactions

Fluctuating wall pressure measurements have been made on centerline upstream of a blunt fin in a Mach 5 turbulent boundary layer. By examining the ensemble averaged wall pressure distributions for different separation shock foot positions, it has been shown that local fluctuating wall pressure measurements are due to a distinct pressure distribution, Rho(sub i), which undergoes a stretching and flattening effect as its upstream boundary translates aperiodically between the upstream influence and separation lines. The locations of the maxima and minima in the wall pressure standard deviation can be accurately predicted using this distribution, providing quantitative confirmation of the model. This model also explains the observed cross-correlations and ensemble average measurements within the interaction. Using the Rho(sub i) model, wall pressure signals from under the separated flow region were used to reproduce the position-time history of the separation shock foot. Further, the negative time delay peak in the cross-correlation between the predicted and actual shock foot histories suggests that the separated region fluctuations precede shock foot motion. The unsteady behavior of the primary horseshoe vortex and its relation to the unsteady separation shock are described.

Dolling, David S.↗

Investigation to advance prediction techniques of the low-speed aerodynamics of V/STOL aircraft

A computer program, VSAERO, has been applied to a number of V/STOL configurations with a view to advancing prediction techniques for the low-speed aerodynamic characteristics. The program couples a low-order panel method with surface streamline calculation and integral boundary layer procedures. The panel method--which uses piecewise constant source and doublet panels-includes an iterative procedure for wake shape and models boundary layer displacement effect using the source transpiration technique. Certain improvements to a basic vortex tube jet model were installed in the code prior to evaluation. Very promising results were obtained for surface pressures near a jet issuing at 90 deg from a flat plate. A solid core model was used in the initial part of the jet with a simple entrainment model. Preliminary representation of the downstream separation zone significantly improve the correlation. The program accurately predicted the pressure distribution inside the inlet on the Grumman 698-411 design at a range of flight conditions. Furthermore, coupled viscous/potential flow calculations gave very close correlation with experimentally determined operational boundaries dictated by the onset of separation inside the inlet. Experimentally observed degradation of these operational boundaries between nacelle-alone tests and tests on the full configuration were also indicated by the calculation. Application of the program to the General Dynamics STOL fighter design were equally encouraging. Very close agreement was observed between experiment and calculation for the effects of power on pressure distribution, lift and lift curve slope.

Maskew, B.↗

Visualization of Secondary Flow Development in High Aspect Ratio Channels with Curvature

The results of an experimental project to visually examine the secondary flow structure that develops in curved, high aspect-ratio rectangular channels are presented. The results provide insight into the fluid dynamics within high aspect ratio channels. A water flow test rig constructed out of plexiglass, with an adjustable aspect ratio, was used for these experiments. Results were obtained for a channel geometry with a hydraulic diameter of 10.6 mm (0.417 in.), an aspect ratio of 5.0, and a hydraulic radius to curvature radius ratio of 0.0417. Flow conditions were varied to achieve Reynolds numbers up to 5,100. A new particle imaging velocimetry technique was developed which could resolve velocity information from particles entering and leaving the field of view. Time averaged secondary flow velocity vectors, obtained using this velocimetry technique, are presented for 30 degrees, 60 degrees, and 90 degrees into a 180 degrees bend and at a Reynolds number of 5,100. The secondary flow results suggest the coexistence of both the classical curvature induced vortex pair flow structure and the eddies seen in straight turbulent channel flow.

Meyer, Michael L.↗

Prediction of Wind Tunnel Effects on the Installed F/A-18A Inlet Flow Field at High Angles-of-attack

NASA Lewis is currently engaged in a research effort as a team member of the High Alpha Technology Program (HATP) within NASA. This program utilizes a specially equipped F/A-18A, the High Alpha Research Vehicle (HARV), in an ambitious effort to improve the maneuverability of high performance military aircraft at low subsonic speed, high angle of attack conditions. The overall objective of the Lewis effort is to develop inlet technology that will ensure efficient airflow delivery to the engine during these maneuvers. One part of the Lewis approach utilizes computational fluid dynamics codes to predict the installed performance of inlets for these highly maneuverable aircraft. Wind tunnel tests were a major component of the Lewis program. Since the available wind tunnel was small (9 x 15 ft) as compared to the scale of the model of the F/A-18A (19.78 percent), there were questions about the capability to obtain useful inlet performance data. The blockage effects were expected to be very large. This report represents the results of an analysis to determine how the wind tunnel walls effect inlet performance at several angles of attack. The predictions for the external particle traces along the fuselage indicate the influence of the wind tunnel side wall under the model is greater at 30 deg angle of attack than at 50 deg angle of attack on the under Leading Edge Extension (LEX) vortex trajectory. The side wall above the model appears to have negligible influence on the under LEX vortex. This may be due to the LEX acting as 'shield' to the upper wall effects. As expected, the wind tunnel has a significant influence on the external forces. The lift and drag coefficients increase significantly for the wind tunnel model as compared to free stream conditions. The wind tunnel had a small effect on the inlet recovery and on inlet total pressure distortion patterns. The predicted recoveries for the wind tunnel model are within one percentage point of the model recoveries in free stream conditions.

Smith, Crawford F.↗

Software Aids Visualization of Computed Unsteady Flow

Unsteady Flow Analysis Toolkit (UFAT) is a computer program that synthesizes motions of time-dependent flows represented by very large sets of data generated in computational fluid dynamics simulations. Prior to the development of UFAT, it was necessary to rely on static, single-snapshot depictions of time-dependent flows generated by flow-visualization software designed for steady flows. Whereas it typically takes weeks to analyze the results of a largescale unsteady-flow simulation by use of steady-flow visualization software, the analysis time is reduced to hours when UFAT is used. UFAT can be used to generate graphical objects of flow visualization results using multi-block curvilinear grids in the format of a previously developed NASA data-visualization program, PLOT3D. These graphical objects can be rendered using FAST, another popular flow visualization software developed at NASA. Flow-visualization techniques that can be exploited by use of UFAT include time-dependent tracking of particles, detection of vortex cores, extractions of stream ribbons and surfaces, and tetrahedral decomposition for optimal particle tracking. Unique computational features of UFAT include capabilities for automatic (batch) processing, restart, memory mapping, and parallel processing. These capabilities significantly reduce analysis time and storage requirements, relative to those of prior flow-visualization software. UFAT can be executed on a variety of supercomputers.

Kao, David↗

Flow Field Measurement of Mixing Driven by Buoyancy

Mixing driven by buoyancy-induced flows inside a cavity consists of stretching and folding of an interface. Measurement of the flow field using particle imaging velocimetry shows that during stretching the flow field has a single elliptic point, thus dominated by a single vortex. However, global bifurcation that results in folding introduces a hyperbolic point whereby the flow field degenerates to multiple vortex interactions. The short-lived coherent structure observed during mixing which results in the Rayleigh- Taylor morphology is attributed to vortex interactions. The mixing characteristics of non-homogeneous fluids driven by buoyancy are important towards understanding transport phenomenon in a microgravity environment. Mixing consists of stretching and folding of an interface due to a flow field whose intensity depends on the body force. For miscible liquids, the characteristic of the flow field determines whether mass transport is governed by diffusion or bulk stirring which induces mixing. For technologically important processes, transport of mass is governed by the coupling of the body force to scalar gradients such as concentration and or temperature' 2 3 . In order to lend insight into these classes of problems we consider a model experimental system to study mixing driven by buoyancy-induced flows. The characteristics of mixing is addressed from detail measurements of the flow field using particle imaging velocimetry (PIV), and its corresponding interface dynamics using image processing techniques.

Batur, C.↗

Explicit and implicit compact high-resolution shock-capturing methods for multidimensional Euler equations 1: Formulation

Two classes of explicit compact high-resolution shock-capturing methods for the multidimensional compressible Euler equations for fluid dynamics are constructed. Some of these schemes can be fourth-order accurate away from discontinuities. For the semi-discrete case their shock-capturing properties are of the total variation diminishing (TVD), total variation bounded (TVB), total variation diminishing in the mean (TVDM), essentially nonoscillatory (ENO), or positive type of scheme for 1-D scalar hyperbolic conservation laws and are positive schemes in more than one dimension. These fourth-order schemes require the same grid stencil as their second-order non-compact cousins. One class does not require the standard matrix inversion or a special numerical boundary condition treatment associated with typical compact schemes. Due to the construction, these schemes can be viewed as approximations to genuinely multidimensional schemes in the sense that they might produce less distortion in spherical type shocks and are more accurate in vortex type flows than schemes based purely on one-dimensional extensions. However, one class has a more desirable high-resolution shock-capturing property and a smaller operation count in 3-D than the other class. The extension of these schemes to coupled nonlinear systems can be accomplished using the Roe approximate Riemann solver, the generalized Steger and Warming flux-vector splitting or the van Leer type flux-vector splitting. Modification to existing high-resolution second- or third-order non-compact shock-capturing computer codes is minimal. High-resolution shock-capturing properties can also be achieved via a variant of the second-order Lax-Friedrichs numerical flux without the use of Riemann solvers for coupled nonlinear systems with comparable operations count to their classical shock-capturing counterparts. The simplest extension to viscous flows can be achieved by using the standard fourth-order compact or non-compact formula for the viscous terms.

Yee, H. C.↗