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 1,387 records · Page 77

Towards a theory of tropical/midlatitude mass exchange from the earth's surface through the stratosphere

The main focus of this work is to understand the dynamics of mass exchange between the tropics and the midlatitudes and to determine any links between tropospheric exchange and that in the stratosphere. We have approached this problem from two different perspectives. The first is aimed towards understanding the troposphere's role in inducing lower stratospheric tropical/midlatitude exchange. For this project we focus on observational analysis of the lower stratosphere to assess the key regions of transport in/out of the tropics and to what extent this transport is driven by tropospheric processes. The second approach is to determine the extent to which stratospheric processes influence the troposphere. In this project we are performing potential vorticity (PV) inversions to assess the winds induced near the tropopause when the stratospheric polar vortex is displaced equatorward. These are each discussed in more detail in the subsections below. Also, we have organized a session on Tropical/Midlatitude Interaction and Transport at the Fall AGU where we will be showing our latest results.

Hartley, Dana↗

Evidence for Arctic Ozone Depletion in Late February and early March 1994

Significant chemical ozone (O3 ) loss in the 1993-94 Arctic winter occurred mainly during an unusually late cold spell of approximately 10 days in late Feb/early Mar. Over the 30 d period studied (including the cold spell), observed vortex-averaged O3 at 465 K (approximately 40 hPa) decreased by approximately 10%. New three-dimensional, diabatic trajectory calculations show that this observed decrease represents only about half of the net chemical loss (approximately 20%) during the 30 day period. The resupply of lower stratospheric O3 by transport in Feb 1994 was considerably greater than in 1993, when transport masked only about a quarter of the chemical loss in Feb/Mar. The net estimated chemical loss over 30 days in 1994 was comparable to that over the same 30 days in 1993, but mainly occurred at a faster rate during the brief cold spell. These results highlight the impact of Arctic interannual variability on the relative roles of chemistry and dynamics in O3 evolution during recent Arctic winters.

Manney, G. L.↗

Interactive Computer Graphics

Aerospace data analysis tools that significantly reduce the time and effort needed to analyze large-scale computational fluid dynamics simulations have emerged this year. The current approach for most postprocessing and visualization work is to explore the 3D flow simulations with one of a dozen or so interactive tools. While effective for analyzing small data sets, this approach becomes extremely time consuming when working with data sets larger than one gigabyte. An active area of research this year has been the development of data mining tools that automatically search through gigabyte data sets and extract the salient features with little or no human intervention. With these so-called feature extraction tools, engineers are spared the tedious task of manually exploring huge amounts of data to find the important flow phenomena. The software tools identify features such as vortex cores, shocks, separation and attachment lines, recirculation bubbles, and boundary layers. Some of these features can be extracted in a few seconds; others take minutes to hours on extremely large data sets. The analysis can be performed off-line in a batch process, either during or following the supercomputer simulations. These computations have to be performed only once, because the feature extraction programs search the entire data set and find every occurrence of the phenomena being sought. Because the important questions about the data are being answered automatically, interactivity is less critical than it is with traditional approaches.

Kenwright, David↗

The Dynamics of Miscible Interfaces: Simulations

The goal of this experimental/computational investigation (joint with Prof Maxworthy at USC) has been to study the dynamics of miscible interfaces, both from a scientific and a practical point of view, and to prepare a related experiment to be flown on the International Space Station. In order to address these effects, we have focused experimental and computational investigations on miscible displacements in cylindrical capillary tubes, as well as in Hele-Shaw cells. Regarding the flow in a capillary tube, the question was addressed as to whether Korteweg stresses and/or divergence effects can potentially account for discrepancies observed between conventional Stokes flow simulations and experiments for miscible flows in capillary tubes. An estimate of the vorticity and streamfunction fields induced by the Kortewegs stresses was derived, which shows these stresses to result in the formation of a vortex ring structure near the tip of the concentration front. Through this mechanism the propagation velocity of the concentration front is reduced, in agreement with the experimental observations. Divergence effects, on the other hand, were seen to be very small, and they have a negligible influence on the tip velocity. As a result, it can be concluded that they are not responsible for the discrepancies between experiments and conventional Stokes simulations. A further part of our investigation focussed on the development of high-accuracy three-dimensional spectral element simulation techniques for miscible flows in capillary tubes, including the effects of variable density and viscosity. Towards this end, the conservation equations are treated in cylindrical coordinates.

Meiburg, Eckart↗

Effect of azimuthal flow fluctuations on flow and flame dynamics of axisymmetric swirling flames

We report recent studies have clearly shown the important role of swirl fluctuations (or, more precisely, fluctuations in axial vorticity) in the response of premixed flames to flow oscillations. An important implication of this mechanism is that the axial location of the swirler plays a key role in the phase between the acoustic flow excitation source and the resulting axial vorticity fluctuation at the flame. Similar to the previously well recognized role of azimuthal vorticity fluctuations, these swirl fluctuations are vortical and convect at the mean flow velocity, unlike the acoustic flow fluctuations. However, there is a fundamental difference between axial and azimuthal vorticity disturbances in terms of the flow oscillations they induce on the flame. Specifically, azimuthal vorticity disturbances excite radial and axial flow disturbances, while axial vorticity oscillations, in general induce both radial and azimuthal flow fluctuations, but in the axisymmetric case, they only directly excite azimuthal flow fluctuations. The axial vorticity fluctuations do, however, indirectly excite axial and radial velocity fluctuations when the axial vortex tube is tilted off-axis, such as at locations of area expansion. This difference is significant because axisymmetric flames are disturbed only by the velocity component normal to it, which stem from axial and radial velocity components only. This implies that axisymmetric mean flames are not directly affected by azimuthal flow fluctuations, since they are tangential to it. Thus, it is the extent to which the axial vorticity is tilted and rotated that controls the strength of the flow oscillations normal to the flame and, in turn, lead to heat release oscillations. This coupling process is not easily amenable to analytical calculations and, as such, we report here a computational study of the role of these different flow fluctuations on the flame response in an axisymmetric framework. In conclusion, the results indicate that the swirl fluctuations can act as significant source of flame heat release disturbances, due to azimuthal and radial vortex tubes tilting in mixing passage boundary layers and at the expansion into the combustor.

42 ENGINEERING↗

Investigation of the structure and dynamics of the ozone layer in the Eastern Arctic region during EASOE Campaign

Balloon measurements of the vertical distribution of ozone and aerosol were made at Dickson Island (73 deg N, 81 deg E), Kiruna (68 deg N, 20 deg E) and Heiss Island (81 deg N, 58 deg E) from December 1991 to March 1992. To acquire data on the seasonal variability of the vertical ozone distribution, electrochemical ozonesondes ECC-4A were flown three times a week. With ozonesondes on the same balloons, backscattersondes were flown on the average of two or three times per month. Using these instruments, altitude profiles of backscatter ratio were measured at two wavelengths centered at 490 nm and 940 nm. Additionally, at Heiss Island, Dickson, and Yakutsk (63 deg N, 130 deg E) regular total ozone measurements were obtained using Brewer spectrophotometers. Based on measurements of backscatter ratio it was found that after the Pinatubo eruption in June 1991 significant amount of stratospheric aerosols were formed and transported to the Arctic before the polar vortex was well developed. Analysis of ozone data has shown a deep decrease of ozone concentration in the lower stratosphere in times of intensive transportation of air masses from low latitudes to the polar region in the second half of January and also for some periods in February and March of 1992. When the values of backscatter ratio beta were more than 8-10 at a wavelength of 940 nm strong anticorrelation occurred between aerosol loading and ozone concentration in the lower stratosphere. At 50-70 deg N, the mean monthly values of total ozone in winter-spring 1992 proved to be much lower than the climatic mean values.

Khattatov, V.↗

Hydrodynamic and Structural Response Modeling of a Prototype Floating Membrane Reservoir System for Pumped Storage Hydropower

Hydrodynamic, computational fluid dynamics, and finite-element modeling were performed for a novel floating membrane reservoir system design for closed-loop pumped storage hydropower application. The conceptual design, which is now protected under an invention disclosure with a patent pending, offers a potential low-cost, low-impact solution to address the high costs, long investment return periods, and environmental disruptions encountered with traditional pumped storage development while offering modularity to enable replication at many locations. Prior to physical tests on a prototype, this paper documents the response analysis and modeling completed to simulate hydraulic and structural performance under different reservoir deployment and alignment arrangements, evaluate system stability, and refine the conceptual design. The results indicate that the excitation frequency from vortex shedding is at least an order of magnitude lower than the water sloshing frequencies in the reservoir, the structural natural frequency of the entire reservoir, and the vibrational frequencies of the side membranes, although excitation frequencies from other sources could cause mechanical resonance under certain conditions. To control destabilizing effects and prevent rocking motion and vibration, the authors conclude that the design could be improved by including a support structure around the floating reservoir, which could also provide floating-membrane containment, improve safety, and facilitate vertical reservoir movement. These conclusions, based on hydrodynamic and structural response modeling, help define specifications for upcoming full-scale prototype construction, deployment, and testing. This study demonstrates standard modeling technique application to an innovative design for which similar applications have not been previously evaluated. As a result, the refined design is capable of improving the scalability and feasibility of pumped storage hydropower in the United States and will be considered for commercialization following prototype testing.

13 HYDRO ENERGY↗

Flight, Wind-Tunnel, and Computational Fluid Dynamics Comparison for Cranked Arrow Wing (F-16XL-1) at Subsonic and Transonic Speeds

Geometrical, flight, computational fluid dynamics (CFD), and wind-tunnel studies for the F-16XL-1 airplane are summarized over a wide range of test conditions. Details are as follows: (1) For geometry, the upper surface of the airplane and the numerical surface description compare reasonably well. (2) For flight, CFD, and wind-tunnel surface pressures, the comparisons are generally good at low angles of attack at both subsonic and transonic speeds, however, local differences are present. In addition, the shock location at transonic speeds from wind-tunnel pressure contours is near the aileron hinge line and generally is in correlative agreement with flight results. (3) For boundary layers, flight profiles were predicted reasonably well for attached flow and underneath the primary vortex but not for the secondary vortex. Flight data indicate the presence of an interaction of the secondary vortex system and the boundary layer and the boundary-layer measurements show the secondary vortex located more outboard than predicted. (4) Predicted and measured skin friction distributions showed qualitative agreement for a two vortex system. (5) Web-based data-extraction and computational-graphical tools have proven useful in expediting the preceding comparisons. (6) Data fusion has produced insightful results for a variety of visualization-based data sets.

John E Lamar↗

Full-Span Tiltrotor Aeroacoustic Model (TRAM) Overview and 40- by 80-Foot Wind Tunnel Test

Most helicopter data trends cannot be extrapolated to tiltrotors because blade geometry and aerodynamic behavior, as well as rotor and fuselage interactions, are significantly different for tiltrotors. A tiltrotor model has been developed to investigate the aeromechanics of tiltrotors, to develop a comprehensive database for validating tiltrotor analyses, and to provide a research platform for supporting future tiltrotor designs. The Full-Span Tiltrotor Aeroacoustic Model (FS TRAM) is a dual-rotor, powered aircraft model with extensive instrumentation for measurement of structural and aerodynamic loads. This paper will present the Full-Span TRAM test capabilities and the first set of data obtained during a 40- by 80-Foot Wind Tunnel test conducted in late 2000 at NASA Ames Research Center. The Full-Span TRAM is a quarter-scale representation of the V-22 Osprey aircraft, and a heavily instrumented NASA and U.S. Army wind tunnel test stand. Rotor structural loads are monitored and recorded for safety-of-flight and for information on blade loads and dynamics. Left and right rotor balance and fuselage balance loads are monitored for safety-of-flight and for measurement of vehicle and rotor aerodynamic performance. Static pressure taps on the left wing are used to determine rotor/wing interactional effects and rotor blade dynamic pressures measure blade airloads. All of these measurement capabilities make the FS TRAM test stand a unique and valuable asset for validation of computational codes and to aid in future tiltrotor designs. The Full-Span TRAM was tested in the NASA Ames Research Center 40- by 80-Foot Wind Tunnel from October through December 2000. Rotor and vehicle performance measurements were acquired in addition to wing pressures, rotor acoustics, and Laser Light Sheet (LLS) flow visualization data. Hover, forward flight, and airframe (rotors off) aerodynamic runs were performed. Helicopter-mode data were acquired during angle of attack and thrust sweeps for a variety of tunnel speeds. Wake geometry images were acquired using LLS photographs and suggest dual tip vortex formation at low thrust conditions. The full paper will include comparisons to isolated-rotor TRAM data acquired at the Duits-Nederlandse Windtunnel (DNW) in 1998. The FS TRAM has been established as a valuable national asset for tiltrotor research. Data reduction and analysis of the 40- by 80-Foot Wind Tunnel test results are underway. Follow-on testing of the FS TRAM is currently being planned for the NASA Ames 80- by 120-Foot Wind Tunnel in late 2001.

McCluer, Megan S.↗

Integrated Physics-based Modeling and Experiments for Improved Prediction of Combustion Dynamics in Low-Emission Systems

Concurrent simulation and experiment was undertaken to assess the ability of a hybrid RANS-LES model to predict combustion dynamics in a single-element lean direct-inject (LDI) combustor showing self-excited instabilities. High frequency pressure modes produced by Fourier and modal decomposition analysis were compared quantitatively, and trends with equivalence ratio and inlet temperature were compared qualitatively. High frequency OH PLIF and PIV measurements were also taken. Submodels for chemical kinetics and primary and secondary atomization were also tested against the measured behavior. For a point-wise comparison, the amplitudes matched within a factor of two. The dependence on equivalence ratio was matched. Preliminary results from simulation using an 18-reaction kinetics model indicated instability amplitudes closer to measurement. Analysis of the simulations suggested a band of modes around 1400 Hz were due to a vortex bubble breakdown and a band of modes around 6 kHz were due to a precessing vortex core hydrodynamic instability. The primary needs are directly coupled and validated ab initio models of the atomizer free surface flow and the primary atomization processes, and more detailed study of the coupling between the 3D swirling flow and the local thermoacoustics in the diverging venturi section.

combustion↗

A Fluid Solver for Studying Torsional Galloping in Solar-Tracking PV Panel Arrays

Solar-tracking photovoltaic arrays are susceptible to aeroelastic fluttering during high-wind events. This dynamic fluttering behavior can grow in amplitude until the panels enter an unstable mode known as torsional galloping which can lead to panel failure or total array destruction. To better understand the physics of the torsional galloping phenomenon, and to inform the discussion around panel design and recommended panel stow positions during high wind events, a fluid-structure interaction solver composed of a simulated atmospheric boundary layer with simplified panel structural responses was designed. The simulation choices and features of this solver were informed by the geometry and physical properties of an experimental panel array known to exhibit torsional galloping behavior during hind-wind events. These simulations revealed that the torsional galloping instability is driven by a combination of cyclic vortex shedding from the sun-facing side of the panel and the elastic properties of the torque tube linking the panel assemblies. Testing different stow angles across a range of wind speeds indicates that panels are generally more stable when stowed at negative angles where the leading edge is closer to the ground, hypothesized to be due to ground-blocking effects.

41 EE - Solar Energy Technologies Office (EE-4S)↗

A time‐varying formulation of the curled wake model within the FAST.Farm framework

Abstract In this article, we present a time‐varying formulation of the curled wake model that we implemented in FAST.Farm. The curled wake model, originally developed for steady‐state conditions, is used to produce realistic wake profiles behind a wind turbine in yawed (or skewed) conditions. We begin by introducing the key elements of the FAST.Farm framework. Then, after briefly summarizing the original wake dynamics formulation of FAST.Farm based on a polar wake profile, we present the new time‐varying formulation of the curled wake model, compare the two, and highlight the differences with the original curled wake model. After discussing some implementation details, we present different applications with increasing levels of complexity: single turbine with uniform and turbulent inflow, fixed and transient yaw, and multiple turbines. We verify our results using the original FAST.Farm implementation and large‐eddy simulations. The results with the new curled wake model are improved compared to the original implementation, as they include cross‐flow velocities and wake asymmetry. Yet, large‐eddy simulation results show a more pronounced lateral convection of the wake and a stronger concentration of vorticity at the top vortex. The new curled wake implementation in FAST.Farm should enable the calculation of not only generator power but also wind turbine structural loads for applications involving intentional or unintentional skewed flow and wind‐farm control involving wake steering.

17 WIND ENERGY↗

PyroCbs from Australia Fires and its Impact Using Satellite Observations from CrIS and TROPOMI and Reanalysis Data

Pyrocumulonimbus (pyroCb) clouds are thunder clouds created by intense heat from the Earth’s surface. They are formed similarly to cumulonimbus clouds, but the intense heat that results in the vigorous updraft comes from fire, either large wildfires or volcanic eruptions. Australia’s unprecedented fire disasters at the end of 2019 to early 2020 emitted huge amounts of carbon monoxide (CO) and fire aerosol particles to the atmosphere, particularly during the pyroCb outbreak that occurred in southeast Australia between 29 December 2019 and 4 January 2020. It was estimated that at least 18 pyroCbs were generated during this episode, and some of them injected ice, smoke, and biomass burning gases above the local tropopause. An unprecedented abundance of H2O and CO in the stratosphere, and the displacement of background ozone (O3) and N2O from rapid ascent of air from the troposphere and lower stratosphere were found from satellite observations. Some other studies also found that the fire emissions and their long-range transport resulted in stratospheric aerosol, temperature, and O3 anomalies after the 2020 Australian bushfires and altered the Antarctic ozone and vortex, posing great impact to local air quality and climate change. Further study on the atmospheric thermodynamic status of atmosphere associated with these pyroCbs, and the change of the cloud properties and trace gases during this unprecedented Australia fires will be made using a new single Field of View (SFOV) Sounder Atmospheric Products (SiFSAP). SiFSAP was developed by NASA using the Cross-track Infrared Sounder (CrIS) and Advanced Technology Microwave Sounder (ATMS) onboard SNPP and JPSS-1, and will soon be available to the public at NASA DAAC. Since this product has a spatial resolution of 15 km at nadir, which is better than most global weather and climate models and other current operational sounding products, a process-oriented analysis of the dynamic transport of CO and fire plumes during this unprecedented fire disasters will be made in this study. Based on a Principal Component Radiative Transfer Model (PCRTM) and an optimized estimation retrieval algorithm, a simultaneously retrieval is made using the whole spectral information measured by CrIS, and the derived SiFSAP include temperature, water vapor, trace gases (such as O3, CO2, CO, CH4 and N2O), cloud properties and surface properties. Use of ATMS together with CrIS allows SiFSAP to get accurate retrieval products under thick pyroCb conditions, and an algorithm to detect pyroCb based on the hyperspectral infrared sounder spectrum from CrIS will be developed and verified. In addition to SiFSAP sounding products, other products like CO, O3, NO2 from TROPOMI, O3 from OMPS will be used for retrospective analysis. The wind fields from the NASA’s Modern-Era Retrospective Analysis for Research and Applications Version-2 (MERRA-2) and ERA5 will be used to characterize the transport, and the SiFSAP temperature and water vapor profiles within and around pyroCbs will be compared with MERRA-2 and ERA5 products.

Xiaozhen (Shawn) Xiong↗

Some numerical and physical aspects of unsteady Navier-Stokes computations over airfoils using dynamic meshes

An upwind-biased implicit approximate factorization algorithm is applied to several unsteady flows on dynamic meshes. The thin-layer form of the compressible Navier-Stokes equations is used to solve both laminar and turbulent flows over airfoils pitching about the quarter chord. Numerical aspects of the solutions are investigated, including grid and time step effects. Two methods for determining fluxes - flux-vector splitting and flux-difference splitting - are compared. Flux-difference splitting predicts results more accurately than flux-vector splitting on a coarse mesh, but both methods agree on a fine mesh. Physical aspects of the computations are also examined. An equilibrium turbulent boundary layer model computes generally better unsteady results in comparison with experiment than a nonequilibrium model for the transonic case analyzed. Also, the size and location of the primary shed vortex for an airfoil pitching up at a constant rate is calculated in good agreement with experiment for two pitch rates.

Rumsey, Christopher L.↗

Overview of the NASA Glenn Flux Reconstruction Based High-Order Unstructured Grid Code

A computational fluid dynamics code based on the flux reconstruction (FR) method is currently being developed at NASA Glenn Research Center to ultimately provide a large- eddy simulation capability that is both accurate and efficient for complex aeropropulsion flows. The FR approach offers a simple and efficient method that is easy to implement and accurate to an arbitrary order on common grid cell geometries. The governing compressible Navier-Stokes equations are discretized in time using various explicit Runge-Kutta schemes, with the default being the 3-stage/3rd-order strong stability preserving scheme. The code is written in modern Fortran (i.e., Fortran 2008) and parallelization is attained through MPI for execution on distributed-memory high-performance computing systems. An h- refinement study of the isentropic Euler vortex problem is able to empirically demonstrate the capability of the FR method to achieve super-accuracy for inviscid flows. Additionally, the code is applied to the Taylor-Green vortex problem, performing numerous implicit large-eddy simulations across a range of grid resolutions and solution orders. The solution found by a pseudo-spectral code is commonly used as a reference solution to this problem, and the FR code is able to reproduce this solution using approximately the same grid resolution. Finally, an examination of the code's performance demonstrates good parallel scaling, as well as an implementation of the FR method with a computational cost/degree- of-freedom/time-step that is essentially independent of the solution order of accuracy for structured geometries.

High-Order Methods↗

Particle concentration at planet-induced gap edges and vortices. I. Inviscid three-dimensional hydro disks

We perform a systematic study of the dynamics of dust particles in protoplanetary disks with embedded planets using global two-dimensional and three-dimensional inviscid hydrodynamic simulations. Lagrangian particles have been implemented into the magnetohydrodynamic code Athena with cylindrical coordinates. We find two distinct outcomes depending on the mass of the embedded planet. In the presence of a low-mass planet (8 M {sub ⊕}), two narrow gaps start to open in the gas on each side of the planet where the density waves are shocked. These shallow gaps can dramatically affect particle drift speed and cause significant, roughly axisymmetric dust depletion. On the other hand, a more massive planet (>0.1 M{sub J} ) carves out a deeper gap with sharp edges, which are subject to Rossby wave instability leading to vortex formation. Particles with a wide range of sizes (0.02 < Ωt{sub s} < 20) are trapped and settle to the midplane in the vortex, with the strongest concentration for particles with Ωt{sub s} ∼ 1. The dust concentration is highly elongated in the φ direction, and can be as wide as four disk scale heights in the radial direction. Dust surface density inside the vortex can be increased by more than a factor of 10{sup 2} in a very non-axisymmetric fashion. For very big particles (Ωt{sub s} >> 1) we find strong eccentricity excitation, in particular around the planet and in the vicinity of the mean motion resonances, facilitating gap openings there. Our results imply that in weakly turbulent protoplanetary disk regions (e.g., the {sup d}ead zone{sup )} dust particles with a very wide range of sizes can be trapped at gap edges and inside vortices induced by planets with M{sub p} < M{sub J} , potentially accelerating planetesimal and planet formation there, and giving rise to distinctive features that can be probed by ALMA and the Extended Very Large Array.

79 ASTRONOMY AND ASTROPHYSICS↗

High resolution numerical simulations of methane pool fires using adaptive mesh refinement

The ability to accurately predict the structure and dynamics of pool fires using computational simulations is of great interest in a wide variety of applications, including accidental and wildland fires. However, the presence of physical processes spanning a broad range of spatial and temporal scales poses a significant challenge for simulations of such fires, particularly at conditions near the transition between laminar and turbulent flow. Here, in this study, we examine the transition to turbulence in methane pool fires using high-resolution simulations with multi-step finite rate chemistry, where adaptive mesh refinement (AMR) is used to directly resolve small-scale flow phenomena. We perform three simulations of methane pool fires, each with increasing diameter, corresponding to increasing inlet Reynolds and Richardson numbers. As the diameter increases, the flow transitions from organized vortex roll-up via the puffing instability to much more chaotic mixing associated with finger formation along the shear layer and core collapse near the inlet. These effects combine to create additional mixing close to the inlet, thereby enhancing fuel consumption and causing more rapid acceleration of the fluid above the pool. We also make comparisons between the transition to turbulence and core collapse in the present pool fires and in inert helium plumes, which are often used as surrogates for the study of buoyant reacting flows.

42 ENGINEERING↗