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 883 records · Page 49

New Capabilities and Improvements to the High-Order Glenn Flux Reconstruction Code

The Glenn Flux Reconstruction (GFR) code is a computational fluid dynamics (CFD) code under development at NASA Glenn Research Center. GFR is based on the high-order flux reconstruction (FR) method and provides a large-eddy simulation (LES) capability that is both accurate and efficient for complex aeropropulsion flows. Three significant new capabilities have been added to the code that improve its performance and functionality. First, a variety of explicit Runge-Kutta methods, including some with adaptive time stepping, were added to GFR with two methods offering a 33% improvement in time-to-solution. Second, GFR can now utilize fully unstructured, mixed-element meshes to more easily facilitate the grid generation process for complex geometries. Finally, a rotating reference frame capability has been added to GFR for solving rotating turbomachinery problems. A selection of results demonstrating these new capabilities are presented in this work. The Taylor-Green vortex problem is used to verify the new unstructured capability by showing similar accuracy and resolution for all element types. LES of the Turbulent Heat Flux Phase III (THX3) experiment with comparison to another high-order LES code and a popular Reynolds-averaged Navier-Stokes (RANS) code demonstrates the accuracy of the code for complex aeropropulsion flows. Finally, LES of a spacecraft cabin ventilation fan shows the ability of GFR to efficiently establish a fan performance map and identify operating points for further analysis at high orders of accuracy.

High-Order Methods↗

New Capabilities and Improvements to the High-Order Glenn Flux Reconstruction Code

The Glenn Flux Reconstruction (GFR) code is a computational fluid dynamics (CFD) code under development at NASA Glenn Research Center. GFR is based on the high-order flux reconstruction (FR) method and provides a large-eddy simulation (LES) capability that is both accurate and efficient for complex aeropropulsion flows. Three significant new capabilities have been added to the code that improve its performance and functionality. First, a variety of explicit Runge-Kutta methods, including some with adaptive time stepping, were added to GFR with two methods offering a 33% improvement in time-to-solution. Second, GFR can now utilize fully unstructured, mixed-element meshes to more easily facilitate the grid generation process for complex geometries. Finally, a rotating reference frame capability has been added to GFR for solving rotating turbomachinery problems. A selection of results demonstrating these new capabilities are presented in this work. The Taylor-Green vortex problem is used to verify the new unstructured capability by showing similar accuracy and resolution for all element types. LES of the Turbulent Heat Flux Phase III (THX3) experiment with comparison to another high-order LES code and a popular Reynolds-averaged Navier-Stokes (RANS) code demonstrates the accuracy of the code for complex aeropropulsion flows. Finally, LES of a spacecraft cabin ventilation fan shows the ability of GFR to efficiently establish a fan performance map and identify operating points for further analysis at high orders of accuracy.

Direct Numerical Simulations↗

Experiments with phase change thermal energy storage canisters for Space Station Freedom

The solar dynamic power module proposed for the Space Station Freedom (SSF) uses the heat of fusion of a phase change material (PCM) to efficiently store thermal energy for use during eclipse periods. The PCM, a LiF-20CaF2 salt, is contained in annular, metal canisters located in a heat receiver at the focus of a solar concentrator. PCM canister ground-based experiments and analytical heat transfer studies are discussed. The hardware, test procedures, and test results from these experiments are discussed. After more than 900 simulated SSF orbital cycles, no canister cracks or leaks were observed and all data were successfully collected. The effect of 1-g test orientation on canister wall temperatures was generally small while void position was strongly dependent on test orientation and canister cooling. In one test orientation, alternating wall temperature data were measured that supports an earlier theory of oscillating vortex flow in the PCM melt. Analytical canister wall temperatures compared very favorably with experimental temperature data. This illustrates that ground-based canister thermal performance can be predicted well by analyses that employ straight-forward, engineering models of void behavior and liquid PCM free convection.

Kerslake, Thomas W.↗

Seasonal Evolution of Saturn's Polar Temperatures and Composition

The seasonal evolution of Saturn's polar atmospheric temperatures and hydrocarbon composition is derived from a decade of Cassini Composite Infrared Spectrometer (CIRS) 7-16 micrometers thermal infrared spectroscopy. We construct a near-continuous record of atmospheric variability poleward of 60 deg from northern winter/southern summer (2004, Ls = 293 deg) through the equinox (2009, Ls= 0 deg) to northern spring/southern autumn (2014, Ls = 56 deg). The hot tropospheric polar cyclones that are entrained by pro-grade jets within 2-3 deg of each pole, and the hexagonal shape of the north polar belt, are both persistent features throughout the decade of observations. The hexagon vertices rotated westward by approx. equal to 30 deg longitude between March 2007 and April 2013, confirming that they are not stationary in the Voyager-defined System III longitude system as previously thought. Tropospheric temperature contrasts between the cool polar zones (near 80-85 deg) and warm polar belts (near 75-80 deg) have varied in both hemispheres, resulting in changes to the vertical wind shear on the zonal jets in the upper troposphere and lower stratosphere. The extended region of south polar stratospheric emission has cooled dramatically poleward of the sharp temperature gradient near 75 deg S (by approximately -5 K/yr), coinciding with a depletion in the abundances of acetylene (0030 +/- 0.005 ppm/yr) and ethane (0.35 +/- 0.1 ppm/yr), and suggestive of stratospheric upwelling with vertical wind speeds of w approx. equal to +0.1 mm/s. The upwelling appears most intense within 5 deg latitude of the south pole. This is mirrored by a general warming of the northern polar stratosphere (+5 K/yr) and an enhancement in acetylene (0.030 +/- 0.003 ppm/yr) and ethane (0.45 +/- 0.1 ppm/yr) abundances that appears to be most intense poleward of 75 deg N, suggesting subsidence at w approx. equal to -0.15 mm/ s. However, the sharp gradient in stratospheric emission expected to form near 75 deg N by northern summer solstice (2017, Ls = 90 deg) has not yet been observed, so we continue to await the development of a northern summer stratospheric vortex. The peak stratospheric warming in the north occurs at lower pressure levels (p less than 1 mbar) than the peak stratospheric cooling in the south (p greater than 1 mbar). Vertical motions are derived from both the temperature field (using the measured rates of temperature change and the deviations from the expectations of radiative equilibrium models) and hydrocarbon distributions (solving the continuity equation). Vertical velocities tend towards zero in the upper troposphere where seasonal temperature contrasts are smaller, except within the tropospheric polar cyclones where w approx. equal to +0.02 mm/s. North polar minima in tropospheric and stratospheric temperatures were detected in 2008-2010 (lagging one season, or 6-8 years, behind winter solstice); south polar maxima appear to have occurred before the start of the Cassini observations (1-2 years after summer solstice), consistent with the expectations of radiative climate models. The influence of dynamics implies that the coldest winter temperatures occur in the 75-80 deg region in the stratosphere, and in the cool polar zones in the troposphere, rather than at the poles themselves. In addition to vertical motions, we propose that the UV-absorbent polar stratospheric aerosols entrained within Saturn's vortices contribute significantly to the radiative budget at the poles, adding to the localized enhancement in the south polar cooling and north polar warming poleward of +/-75 deg.

polar temperatures↗

Control of Meridional Flow by a Non-Uniform Rotational Magnetic Field

The diffusive mass transfer of species during crystal growth in vertical ampoules is significantly affected by fluid flow in the liquid mother phase (melt). For electrically conductive melts, an elegant way of remotely inducing and controlling this flow is by utilizing a uniform rotational magnetic field (RMF) in the transverse direction. It induces an azimuthal flow which tends to homogenize the thermal and solutal fields. The rotating field also reduces the diffusion boundary layer, stabilizes temperature fluctuations, and promotes better overall crystal growth. For moderate strengths of the applied magnetic field (2-20 m Tesla) with frequencies of up to 400 Hz, the induced secondary meridional flow becomes significant. It typically consists of one roll at the bottom of the liquid column and a second roll (vortex) at the top. The flow along the centerline (ampoule axis) is directed from the growing solid (interface) towards the liquid (melt). In case of convex interfaces (e.g. in floating zone crystal growth) such flow behavior is beneficial since it suppresses diffusion at the center. However, for concave interfaces (e.g. vertical Bridgman crystal growth) such a flow tends to exacerbate the situation in making the interface shape more concave. It would be beneficial to have some control of this meridional flow- for example, a single recirculating cell with controllable direction and flow magnitude will make this technique even more attractive for crystal growth. Such flow control is a possibility if a non-uniform PNE field is utilized for this purpose. Although this idea has been proposed earlier, it has not been conclusively demonstrated so far. In this work, we derive the governing equations for the fluid dynamics for such a system and obtain solutions for a few important cases. Results from parallel experimental measurements of fluid flow in a mercury column subjected to non-uniform RMF will also be presented.

Mazuruk, Konstantin↗

HART-II Acoustic Predictions using a Coupled CFD/CSD Method

This paper documents results to date from the Rotorcraft Acoustic Characterization and Mitigation activity under the NASA Subsonic Rotary Wing Project. The primary goal of this activity is to develop a NASA rotorcraft impulsive noise prediction capability which uses first principles fluid dynamics and structural dynamics. During this effort, elastic blade motion and co-processing capabilities have been included in a recent version of the computational fluid dynamics code (CFD). The CFD code is loosely coupled to computational structural dynamics (CSD) code using new interface codes. The CFD/CSD coupled solution is then used to compute impulsive noise on a plane under the rotor using the Ffowcs Williams-Hawkings solver. This code system is then applied to a range of cases from the Higher Harmonic Aeroacoustic Rotor Test II (HART-II) experiment. For all cases presented, the full experimental configuration (i.e., rotor and wind tunnel sting mount) are used in the coupled CFD/CSD solutions. Results show good correlation between measured and predicted loading and loading time derivative at the only measured radial station. A contributing factor for a typically seen loading mean-value offset between measured data and predictions data is examined. Impulsive noise predictions on the measured microphone plane under the rotor compare favorably with measured mid-frequency noise for all cases. Flow visualization of the BL and MN cases shows that vortex structures generated in the prediction method are consist with measurements. Future application of the prediction method is discussed.

Boyd, D. Douglas, Jr.↗

Arcjet thruster research and technology. Phase 1: Executive summary

The principal objective of this two phase program is to conduct the development research required to make the low power arcjet a flight ready technology. Many important results were obtained during Phase 1 to move closer to this objective. Fundamental analyses were performed of the arcjet nozzle, the gas kinetic reaction effects, the thermal environment, and the arc stabilizing vortex. These aided the conceptual understanding of the arcjet and guided design work. A hydrazine (N2H4) arcjet was designed that combined a flight qualified catalyst bed with a modular arcjet. Extensive testing was performed which demonstrated the feasibility of using this propellant in an arcjet for the first time. Startup techniques were developed, stability maintained, material compatibility tests conducted, and performance mapping tests performed. Specific impulse values from 400 to 730 seconds were produced with a non-optimized design. These levels are higher than were originally thought possible and proved that extremely high enthalpy values can be obtained with constricted arc technology. Erosion rate data are promising for lifetime extensions to meet flight application requirements. Power control unit (PCU) development was started with the design and fabrication of a laboratory high switching frequency supply. Valuable data were obtained on PCU operation and on the interaction with the dynamic arc. Phase 2 efforts presently underway are resolving key issues for multi-hundred hour lifetimes, are continuing to investigate arcjet/PCU interactions, and will demonstrate duty cycle N2H4 arcjet/PCU operation in a simulated flight mode for lifetimes consistent with initial applications.

Source record↗

Features of steady magnetospheric convection

The large-scale patterns of ionospheric convection and particle precipitation are described during two intervals of steady magnetospheric convection (SMC) on November 24, 1981. The unique data set used in the analysis includes recordings from the worldwide network of magnetometers and all-sky cameras, global auroral images from the Dynamics Explorer (DE) 1 spacecraft, and particle precipitation data from low-altitude National Oceanic and Atmospheric Administration (NOAA) 6 and NOAA 7 spacecraft. The data show that intense magnetospheric convection continued during more than 10 hours under the steady southward interplanetary magnetic field without any distinct substorm signatures. All data sets available confirmed the stable character of the large-scale magnetospheric configuration during this period. In particular, the magnetic flux threading the polar cap was stable (within 10%) during 3.5 hours of continued DE 1 observations. The dayside cusp was located at an unusually low latitude (70 deg CGL). The nightside auroral pattern consisted of two distinct regions. The diffuse aurora in the equatorward half of the expanded (10 deg wide) auroral oval was well-separated from the bright, active auroral forms found in the vicinity of the poleward boundary of the oval. The twin-vortex convection pattern had no signature of the Harang discontinuity; its nightside 'convection throat' was spatially coincident with the poleward active auroras. This region of the auroral oval was identified as the primary site of the short-lived transient activations during the SMC intervals. The energetic particle observations show that the auroral precipitation up to its high-latitude limit is on closed field lines and that particle acceleration up to greater than 30-keV energy starts close to this limit. The isotropic boundaries of the greater than 30-keV protons and electrons were found close to each other, separating regions of discrete and diffuse precipitation. This suggests that these precipitation types originate on the very taillike and very dipolelike field lines, respectively.

Yahnin, A.↗

The Evolution of Friction Stir Welding Theory at Marshall Space Flight Center

From 1995 to the present the friction stir welding (FSW) process has been under study at Marshall Space Flight Center (MSFC). This is an account of the progressive emergence of a set of conceptual tools beginning with the discovery of the shear surface, wiping metal transfer, and the invention of a kinematic model and making possible a treatment of both metallurgical structure formation and process dynamics in friction stir welding from a unified point of view. It is generally observed that the bulk of the deformation of weld metal around the FSW pin takes place in a very narrow, almost discontinuous zone with high deformation rates characteristic of metal cutting. By 1999 it was realized that this zone could be treated as a shear surface like that in simple metal cutting models. At the shear surface the seam is drawn out and compressed and pressure and flow conditions determine whether or not a sound weld is produced. The discovery of the shear surface was followed by the synthesis of a simple 3- flow kinematic model of the FSW process. Relative to the tool the flow components are: (1) an approaching translational flow at weld speed V, (2) a rotating cylindrical plug flow with the angular velocity of the tool , and (3) a relatively slow ring vortex flow (like a smoke ring) encircling the tool and driven by shoulder scrolls and pin threads. The rotating plug flow picks up an element of weld metal, rotates it around with the tool, and deposits it behind the tool ( wiping metal transfer ); it forms plan section loops in tracers cut through by the tool. Radially inward flow from the ring vortex component retains metal longer in the rotating plug and outward flow expels metal earlier; this interaction forms the looping weld seam trace and the tongue and groove bimetallic weld contour. The radial components of the translational and ring vortex flows introduce parent metal intrusions into the small grained nugget material close to the tool shoulder; if this feature is pronounced, nugget collapse may result. Certain weld features, in particular internal banding seen in transverse section as onion rings and associated surface ridges called tool marks , have long implied an oscillation flow component, but have only recently been attributed in the literature to tool eccentricity. Rotating plug shape, typically a hollow cylinder flared at the end where it sticks to the shoulder, varies as pressure distribution on the tool determines where sticking occurs. Simplified power input estimates balanced against heat loss estimates give reasonable temperature estimates, explain why the power requirement changes hardly at all over a wide range of RPM s, and yield isotherms that seem to fall along boundaries of parameter windows of operation.

Nunes, Arthur C.↗

Numerical design of advanced multi-element airfoils

The current study extends the application of computational fluid dynamics to three-dimensional high-lift systems. Structured, overset grids are used in conjunction with an incompressible Navier-Stokes flow solver to investigate flow over a two-element high-lift configuration. The computations were run in a fully turbulent mode using the one-equation Baldwin-Barth turbulence model. The geometry consisted of an unswept wing which spanned a wind tunnel test section. Flows over full and half-span Fowler flap configurations were computed. Grid resolution issues were investigated in two dimensional studies of the flapped airfoil. Results of the full-span flap wing agreed well with experimental data and verified the method. Flow over the wing with the half-span was computed to investigate the details of the flow at the free edge of the flap. The results illustrated changes in flow streamlines, separation locations, and surface pressures due to the vortex shed from the flap edge.

Mathias, Donovan L.↗

Computational Aeroacoustic Analysis System Development

Many industrial and commercial products operate in a dynamic flow environment and the aerodynamically generated noise has become a very important factor in the design of these products. In light of the importance in characterizing this dynamic environment, Rocketdyne has initiated a multiyear effort to develop an advanced general-purpose Computational Aeroacoustic Analysis System (CAAS) to address these issues. This system will provide a high fidelity predictive capability for aeroacoustic design and analysis. The numerical platform is able to provide high temporal and spatial accuracy that is required for aeroacoustic calculations through the development of a high order spectral element numerical algorithm. The analysis system is integrated with well-established CAE tools, such as a graphical user interface (GUI) through PATRAN, to provide cost-effective access to all of the necessary tools. These include preprocessing (geometry import, grid generation and boundary condition specification), code set up (problem specification, user parameter definition, etc.), and postprocessing. The purpose of the present paper is to assess the feasibility of such a system and to demonstrate the efficiency and accuracy of the numerical algorithm through numerical examples. Computations of vortex shedding noise were carried out in the context of a two-dimensional low Mach number turbulent flow past a square cylinder. The computational aeroacoustic approach that is used in CAAS relies on coupling a base flow solver to the acoustic solver throughout a computational cycle. The unsteady fluid motion, which is responsible for both the generation and propagation of acoustic waves, is calculated using a high order flow solver. The results of the flow field are then passed to the acoustic solver through an interpolator to map the field values into the acoustic grid. The acoustic field, which is governed by the linearized Euler equations, is then calculated using the flow results computed from the flow solver.

Hadid, A.↗

2nd-Order CESE Results For C1.4: Vortex Transport by Uniform Flow

The Conservation Element and Solution Element (CESE) method was used as implemented in the NASA research code ez4d. The CESE method is a time accurate formulation with flux-conservation in both space and time. The method treats the discretized derivatives of space and time identically and while the 2nd-order accurate version was used, high-order versions exist, the 2nd-order accurate version was used. In regards to the ez4d code, it is an unstructured Navier-Stokes solver coded in C++ with serial and parallel versions available. As part of its architecture, ez4d has the capability to utilize multi-thread and Messaging Passage Interface (MPI) for parallel runs.

Space-time CE/SE method↗

The Atmospheric Dynamics of Jupiter, Saturn, and Titan

Comparative studies of Jupiter and Saturn often emphasize their similarities, but recent observations have highlighted important differences. The stratospheres of both planets exhibit an equatorial oscillation reminiscent of that in Earth's middle atmosphere. Jupiter's oscillation has a 4-5 year period, not linked to its season, and it has been modeled as an analog to the terrestrial quasi-biennial oscillation, driven by the stresses associated with vertically propagating waves. Saturn's equatorial oscillation is nearly semiannual, but wave activity may still be a driver. Jupiter's internal rotation rate is inferred from its steady modulated radio emission. Saturn's internal rotation is more enigmatic. It has been inferred from the modulation of the body's kilometric radio emission, but this period has varied by 1% over the last 25 years. Saturn's equatorial winds are also puzzling, as those inferred from cloud tracking by Cassini and more recent HST observations are weaker than those from Voyager. Whether this is attributable to a difference in altitudes of the tracked clouds in winds with vertical shear or a real temporal change in the winds is not known. Both winter and summer poles of Saturn exhibit very compact circumpolar vortices with warm cores, indicating subsidence. Titan's middle atmosphere is characterized by global cyclostrophic winds, particularly the strong circumpolar vortex in the winter hemisphere. In many ways, the spatial distribution of temperature, gaseous constituents, and condensates is reminiscent of conditions in terrestrial winter vortices, albeit with different chemistry. The meridional contrast in Titan's tropospheric temperatures is small, only a few kelvins.

Flasar, F. M.↗

Observed aspects of atmospheric blocking

An important problem in large-scale atmospheric dynamics is understanding the evolution of these blocking systems. Inspection of a large number of cases of 500 mb cyclonic and anticyclonic blocking vortices during recent winter seasons reveals that the occurrence of these systems is preceded by high amplitude (but not necessarily stationary) planetary waves at 500 mb and by intense synoptic-scale surface cyclone activity. A possible relationship between the cyclones and planetary waves during block evolution is thus hypothesized. Detailed diagnostic study of blocking cases during January 1977, February 1978 and November 1980 discloses that the development of 500 mb anticyclonic blocking vortices is preceded in time by the intense surface cyclogenesis between 500 mb planetary scale troughs and ridges, whereas blocking cyclonic vortices form following cyclogenesis near the axes of 500 mb troughs. A positive feedback between the cyclone and planetary waves is thus implied. The 500 mb height changes accompanying the above processes were analyzed in terms of quasi-geostrophic potential vorticity. During anticyclonic vortex development, there are spatially and temporally persistent calculated quasi-geostrophic height rises due to anticyclonic potential vorticity advections occurring near the cyclone track. Similarly, during cyclonic vortex development, there are spatially and temporally persistent calculated quasi-geostrophic height falls due to cyclonic potential vorticity advections occurring near the cyclone track.

Colucci, S. J.↗

Computations of three-dimensional steady and unsteady viscous incompressible flows

The INS3D family of computational fluid dynamics computer codes is presented. These codes are used to as tools in developing and assessing algorithms for solving the incompressible Navier-Stokes equations for steady-state and unsteady flow problems. This work involves applying the codes to real-world problems involving complex three-dimensional geometries. The algorithms utilized include the method of pseudocompressibility including both central and upwind differencing, several types of artificial dissipation schemes, approximate factorization, and an implicit line-relaxation scheme. These codes have been validated using a wide range of problems including flow over a backward-facing step, driven cavity flow, flow through various types of ducts, and steady and unsteady flow over a circular cylinder. Many diverse flow applications have been solved using these codes including parts of the Space Shuttle Main Engine, problems in naval hydrodynamics, low-speed aerodynamics, and biomedical fluid flows. The presentation details several of these, including the flow through a Space Shuttle Main Engine inducer, vortex shedding behind a circular cylinder, and flow through an artificial heart.

Kwak, Dochan↗

Assessment of Pseudoshock Models Against Experiment in a Low-Aspect-Ratio Isolator

A highly confined shock train is investigated in a direct-connect isolator facility with a Mach 2 inflow and a constant-area low-aspect-ratio rectangular test section. High-speed schlieren imaging, wall static pressure measurements, surface oil-flow visualization, and particle image velocimetry from this isolator are synthesized into a three-dimensional schematic of the shock train structure. Against this, the prevailing pseudoshock models in the literature are assessed to evaluate the validity of their underlying assumptions. None of the prevailing pseudoshock models are found to simultaneously model the pressure and Mach number profiles, indicating a gap in the model formation and underlying assumptions when applied to the experimental isolator of interest. The presence of distortion in the isolator flowfield, such as a wall-bounded vortex, is found to skew the structure of the shock train, altering the strength and distribution of the compressive pressure gradient. It is further observed that the separated flow morphology surrounding the shock train is not monolithic, as is typically assumed, adjusting the balance of compressive forces within the shock cells. These findings lead to the conclusion that existing flux-conserved modeling approaches require modification to be effective in distorted and highly confined cases, including closure models that capture the three-dimensional distorted structure of the approach flow and its evolution along the shock train.

Shock Train↗

Aerodynamic properties of turbulent combustion fields

Flow fields involving turbulent flames in premixed gases under a variety of conditions are modeled by the use of a numerical technique based on the random vortex method to solve the Navier-Stokes equations and a flame propagation algorithm to trace the motion of the front and implement the Huygens principle, both due to Chorin. A successive over-relaxation hybrid method is applied to solve the Euler equation for flows in an arbitrarily shaped domain. The method of images, conformal transformation, and the integral-equation technique are also used to treat flows in special cases, according to their particular requirements. Salient features of turbulent flame propagation in premixed gases are interpreted by relating them to the aerodynamic properties of the flow field. Included among them is the well-known cellular structure of flames stabilized by bluff bodies, as well as the formation of the characteristic tulip shape of flames propagating in ducts. In its rudimentary form, the mechanism of propagation of a turbulent flame is shown to consist of: (1) rotary motion of eddies at the flame front, (2) self-advancement of the front at an appropriate normal burning speed, and (3) dynamic effects of expansion due to exothermicity of the combustion reaction. An idealized model is used to illustrate these fundamental mechanisms and to investigate basic aerodynamic features of flames in premixed gases. The case of a confined flame stabilized behind a rearward-facing step is given particular care and attention. Solutions are shown to be in satisfactory agreement with experimental results, especially with respect to global properties such as the average velocity profiles and reattachment length.

Hsiao, C. C.↗

Simulating Virtual Terminal Area Weather Data Bases for Use in the Wake Vortex Avoidance System (Wake VAS) Prediction Algorithm

During the research project, sounding datasets were generated for the region surrounding 9 major airports, including Dallas, TX, Boston, MA, New York, NY, Chicago, IL, St. Louis, MO, Atlanta, GA, Miami, FL, San Francico, CA, and Los Angeles, CA. The numerical simulation of winter and summer environments during which no instrument flight rule impact was occurring at these 9 terminals was performed using the most contemporary version of the Terminal Area PBL Prediction System (TAPPS) model nested from 36 km to 6 km to 1 km horizontal resolution and very detailed vertical resolution in the planetary boundary layer. The soundings from the 1 km model were archived at 30 minute time intervals for a 24 hour period and the vertical dependent variables as well as derived quantities, i.e., 3-dimensional wind components, temperatures, pressures, mixing ratios, turbulence kinetic energy and eddy dissipation rates were then interpolated to 5 m vertical resolution up to 1000 m elevation above ground level. After partial validation against field experiment datasets for Dallas as well as larger scale and much coarser resolution observations at the other 8 airports, these sounding datasets were sent to NASA for use in the Virtual Air Space and Modeling program. The application of these datasets being to determine representative airport weather environments to diagnose the response of simulated wake vortices to realistic atmospheric environments. These virtual datasets are based on large scale observed atmospheric initial conditions that are dynamically interpolated in space and time. The 1 km nested-grid simulated datasets providing a very coarse and highly smoothed representation of airport environment meteorological conditions. Details concerning the airport surface forcing are virtually absent from these simulated datasets although the observed background atmospheric processes have been compared to the simulated fields and the fields were found to accurately replicate the flows surrounding the airport where coarse verification data were available as well as where airport scale datasets were available.

Kaplan, Michael L.↗