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,459 records · Page 81

Subseasonal to Seasonal Prediction of Wintertime Northern Hemisphere Extratropical Cyclone Activity by S2S and NMME Models

In this study, the prediction of wintertime extratropical cyclone activity from subseasonal to seasonal timescale in current dynamical models' reforecasts is investigated. On seasonal time scales, the North American Multi-Model Ensemble (NMME) models show skillful predictions over the eastern North Pacific, North America and the western North Atlantic with at least 5 months lead. The prediction skill is highly related to El Nino-Southern Oscillation (ENSO), as using the ENSO-related SST pattern gives rise to prediction skill with very similar spatial pattern and amplitude. On subseasonal time scales, models in the Seasonal to Sub-seasonal Prediction (S2S) dataset have skillful predictions up to 4 weeks lead over regions from the eastern North Pacific to the western Atlantic, as well as northern Europe, the eastern Atlantic and East Asia. Generally, forecast skill improves with a larger ensemble size. The subseasonal prediction skill from the Pacific to the western Atlantic is related to ENSO, and that over eastern Atlantic, Europe and East Asia are associated with stratospheric polar vortex anomalies. Current models do not show much skill from the Madden-Julian Oscillation (MJO), as the MJO impact on extratropical cyclone activity is not well captured by the models. European Centre for Medium-Range Weather Forecasts (ECMWF) model has the highest single model subseasonal prediction skill. The prediction skill in the ECMWF model is higher than its estimated potential predictability, likely because the signal-to-noise ratio is too low in the model hindcasts.

54 ENVIRONMENTAL SCIENCES↗

Effect of Controlled Artificial Disorder on the Magnetic Properties of EuFe 2 (As 1– x P x ) 2 Ferromagnetic Superconductor

Static (DC) and dynamic (AC, at 14 MHz and 8 GHz) magnetic susceptibilities of single crystals of a ferromagnetic superconductor, EuFe 2 ( As 1 – x P x ) 2 ( x = 0.23), were measured in pristine state and after different doses of 2.5 MeV electron or 3.5 MeV proton irradiation. The superconducting transition temperature, T c ( H ) , shows an extraordinarily large decrease. It starts at T c ( H = 0 ) ≈ 24 K in the pristine sample for both AC and DC measurements, but moves to almost half of that value after moderate irradiation dose. Remarkably, after the irradiation not only T c moves significantly below the FM transition, its values differ drastically for measurements at different frequencies, ≈16 K in AC measurements and ≈12 K in a DC regime. We attribute such a large difference in T c to the appearance of the spontaneous internal magnetic field below the FM transition, so that the superconductivity develops directly into the mixed spontaneous vortex-antivortex state where the onset of diamagnetism is known to be frequency-dependent. We also examined the response to the applied DC magnetic fields and studied the annealing of irradiated samples, which almost completely restores the superconducting transition. Overall, our results suggest that in EuFe 2 ( As 1 – x P x ) 2 superconductivity is affected by local-moment ferromagnetism mostly via the spontaneous internal magnetic fields induced by the FM subsystem. Another mechanism is revealed upon irradiation where magnetic defects created in ordered Eu 2 + lattice act as efficient pairbreakers leading to a significant T c reduction upon irradiation compared to other 122 compounds. On the other hand, the exchange interactions seem to be weakly screened by the superconducting phase leading to a modest increase of T m (less than 1 K) after the irradiation drives T c to below T m . Our results suggest that FM and SC phases coexist microscopically in the same volume.

36 MATERIALS SCIENCE↗

The effect of non-linear mutual friction on pulsar glitch sizes and rise times

ABSTRACT Observations of pulsar glitches have the potential to provide constraints on the dynamics of the high density interior of neutron stars. However, to do so, realistic glitch models must be constructed and compared to the data. We take a step towards this goal by testing non-linear models for the mutual friction force, which is responsible for the exchange of angular momentum between the neutron superfluid and the observable normal component in a glitch. In particular, we consider a non-linear dependence of the drag force on the relative velocity between superfluid vortices and the normal component, in which the contributions of both kelvin and phonon excitations are included. This non-linear model produces qualitatively new features, and is able to reproduce the observed bimodal distribution of glitch sizes in the pulsar population. The model also suggests that the differences in size distributions in individual pulsars may be due to the glitches being triggered in regions with different pinning strengths, as stronger pinning leads to higher vortex velocities and a qualitatively different mutual friction coupling with respect to the weak pinning case. Glitches in pulsars that appear to glitch quasi-periodically with similar sizes may thus be due to the same mechanisms as smaller events in pulsars that have no preferred glitch size, but simply originate in stronger pinning regions, possibly in the core of the star.

79 ASTRONOMY AND ASTROPHYSICS↗

Temporal coherent control of resonant two-photon double ionization of the hydrogen molecule via doubly excited states

Here, we use time-delayed, counter-rotating, circularly polarized few-cycle attosecond nonoverlapping pulses to study the temporal coherent control of the resonant process of two-photon double ionization (TPDI) of hydrogen molecule via doubly excited states for pulse propagation direction along $\hat k$ either parallel or perpendicular to the molecular axis $\hat R$. For $\hat k$∥ $\hat R$ and a pulse carrier frequency of 36 eV resonantly populating the Q 2 1 Π$^+_u$ (1) doubly excited state as well as other 1 Π$^+_u$ doubly excited states, we find that the indirect ionization pathway through these doubly excited states changes the character of the kinematical vortex-shaped momentum distribution produced by the two direct ionization pathways from fourfold to twofold rotational symmetry. This result is similar to what found in TPDI of the He atom involving 1 P$^ o_ {±1}$ doubly excited states; however, angular distributions exhibiting a quantum beat effect between the ground state and a doubly excited state seen for the He atom are observed here for its molecular counterpart with an anomaly in shape and magnitude, not in frequency. The sixfold differential probability integrated over the azimuthal angle of the photoelectron pair shows that this anomaly is due to autoionization decays and quantum beats between doubly excited states. For $\hat k$⟂$\hat R$ and a broadband pulse carrier frequency of 30 eV populating the Q 1 1 Π$^+_u$(1), Q 1 1 Σ$^+_u$(1), Q 2 1 Π$^+_u$(1) , and Q 1 1 Σ$^+_u$ (2) doubly excited states, the momentum distribution is shown to exhibit dynamical electron vortices with four spiral arms, which originates from the interplay between the 1 Δ$^+_ g$ , 1 Π$^+_ g$, and 1 Σ$^+_ g$ dynamical ionization amplitudes. Our treatment within either the adiabatic-nuclei approximation or fixed-nuclei approximation shows that the latter provides a very good account for this correlated process.

74 ATOMIC AND MOLECULAR PHYSICS↗

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.↗

The Large-Scale Vorticity Balance of the Antarctic Continental Margin in a Fine-Resolution Global Simulation

The depth-integrated vorticity budget of a global, eddy-permitting ocean/sea ice simulation over the Antarctic continental margin (ACM) is diagnosed to understand the physical mechanisms implicated in meridional transport. The leading-order balance is between the torques due to lateral friction, nonlinear effects, and bottom vortex stretching, although details vary regionally. Maps of the time-averaged depth-integrated vorticity budget terms and time series of the spatially averaged, depth-integrated vorticity budget terms reveal that the flow in the Amundsen, Bellingshausen, and Weddell Seas and, to a lesser extent, in the western portion of East Antarctica, is closer to an approximate topographic Sverdrup balance (TSB) compared to other segments of the ACM. Correlation and coherence analyses further support these findings, and also show that inclusion of the vorticity tendency term in the response (the planetary vorticity advection and the bottom vortex stretching term) increases the correlation with the forcing (the vertical net stress curl), and also increases the coherence between forcing and response at high frequencies across the ACM, except for the West Antarctic Peninsula. Finally, these findings suggest that the surface stress curl, imparted by the wind and the sea ice, has the potential to contribute to the meridional, approximately cross-slope, transport to a greater extent in the Amundsen, Bellingshausen, Weddell, and part of the East Antarctic continental margin than elsewhere in the ACM.

58 GEOSCIENCES↗

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↗

Review: knots and other new topological effects in liquid crystals and colloids

Abstract Humankind has been obsessed with knots in religion, culture and daily life for millennia, while physicists like Gauss, Kelvin and Maxwell already involved them in models centuries ago. Nowadays, colloidal particles can be fabricated to have shapes of knots and links with arbitrary complexity. In liquid crystals, closed loops of singular vortex lines can be knotted by using colloidal particles and laser tweezers, as well as by confining nematic fluids into micrometer-sized droplets with complex topology. Knotted and linked colloidal particles induce knots and links of singular defects, which can be interlinked (or not) with colloidal particle knots, revealing the diversity of interactions between topologies of knotted fields and topologically nontrivial surfaces of colloidal objects. Even more diverse knotted structures emerge in nonsingular molecular alignment and magnetization fields in liquid crystals and colloidal ferromagnets. The topological solitons include hopfions, skyrmions, heliknotons, torons and other spatially localized continuous structures, which are classified based on homotopy theory, characterized by integer-valued topological invariants and often contain knotted or linked preimages, nonsingular regions of space corresponding to single points of the order parameter space. A zoo of topological solitons in liquid crystals, colloids and ferromagnets promises new breeds of information displays and a plethora of data storage, electro-optic and photonic applications. Their particle-like collective dynamics echoes coherent motions in active matter, ranging from crowds of people to schools of fish. This review discusses the state of the art in the field, as well as highlights recent developments and open questions in physics of knotted soft matter. We systematically overview knotted field configurations, the allowed transformations between them, their physical stability and how one can use one form of knotted fields to model, create and imprint other forms. The large variety of symmetries accessible to liquid crystals and colloids offer insights into stability, transformation and emergent dynamics of fully nonsingular and singular knotted fields of fundamental and applied importance. The common thread of this review is the ability to experimentally visualize these knots in real space. The review concludes with a discussion of how the studies of knots in liquid crystals and colloids can offer insights into topologically related structures in other branches of physics, with answers to many open questions, as well as how these experimentally observable knots hold a strong potential for providing new inspirations to the mathematical knot theory.

Physics↗

Gaseous jet through an outward opening injector: Details of mixing characteristic and turbulence scales

Direct injection (DI) strategy of natural gas (NG) into internal combustion engines (ICE) has led to higher thermal efficiency and lower exhaust emissions. In order to thoroughly understand the most relevant phenomena affecting the performances of such engines, computational fluid dynamics (CFD) plays a key role as an accurate description of the jet evolution and interaction within the combustion chamber is required to that aim. Accurate description of high-pressure gaseous jets is rather challenging at high Mach numbers, as the injected gas is strongly under-expanded once in the ambient, giving room to shocks due to compressibility effects. Also the interaction between shock waves and mixing layers needs to be carefully represented with a multi-dimensional model, calling for substantial computational resources requirements. In this paper a numerical investigation of the behavior of a gaseous jet (Argon) through an outward opening injector has been carried out. A Large Eddy Simulation (LES) approach has been used in order to track the structures derived by the interaction of the injected fuel with the surrounding ambient. Although already good results were obtained using a Reynolds Averaged Navier-Stokes (RANS) approach, the adoption of LES is required to characterize more accurately the jet properties in terms of vortex structures and mixing effectiveness. Finally, the effect of the Nozzle Pressure Ratio (NPR) on the jet evolution has been highlighted in the paper, showing how a higher NPR would give a faster injection process, compromising however the homogeneity of the mixture.

33 ADVANCED PROPULSION SYSTEMS↗

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↗

Simulation of the Multi-Wake Evolution of Two Sandia National Labs/National Rotor Testbed Turbines Operating in a Tandem Layout

The future of wind power systems deployment is in the form of wind farms comprised of scores of such large turbines, most likely at offshore locations. Individual turbines have grown in span from a few tens of meters to today’s large turbines with rotor diameters that dwarf even the largest commercial aircraft. These massive dynamical systems present unique challenges at scales unparalleled in prior applications of wind science research. Fundamental to this effort is the understanding of the wind turbine wake and its evolution. Furthermore, the optimization of the entire wind farm depends on the evolution of the wakes of different turbines and their interactions within the wind farm. In this article, we use the capabilities of the Common ODE Framework (CODEF) model for the analysis of the effects of wake–rotor and wake-to-wake interactions between two turbines situated in a tandem layout fully and partially aligned with the incoming wind. These experiments were conducted in the context of a research project supported by the National Rotor Testbed (NRT) program of Sandia National Labs (SNL). Results are presented for a layout which emulates the turbine interspace and relative turbine emplacement found at SNL’s Scaled Wind Technologies Facility (SWiFT), located in Lubbock, Texas. The evolution of the twin-wake interaction generates a very rich series of secondary transitions in the vortex structure of the combined wake. These ultimately affect the wake’s axial velocity patterns, altering the position, number, intensity, and shape of localized velocity-deficit zones in the wake’s cross-section. This complex distribution of axial velocity patterns has the capacity to substantially affect the power output, peak loads, fatigue damage, and aeroelastic stability of turbines located in subsequent rows downstream on the farm.

Baruah, Apurva (ORCID:0000000252354068)↗

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.↗

Three-dimensional topological magnetic monopoles and their interactions in a ferromagnetic meta-lattice

Topological magnetic monopoles (TMMs), also known as hedgehogs or Bloch points, are three-dimensional (3D) nonlocal spin textures that are robust to thermal and quantum fluctuations due to the topology protection. Although TMMs have been observed in skyrmion lattices, spinor Bose–Einstein condensates, chiral magnets, vortex rings and vortex cores, it has been difficult to directly measure the 3D magnetization vector field of TMMs and probe their interactions at the nanoscale. Here, we report the creation of 138 stable TMMs at the specific sites of a ferromagnetic meta-lattice at room temperature. We further develop soft x-ray vector ptycho-tomography to determine the magnetization vector and emergent magnetic field of the TMMs with a 3D spatial resolution of 10 nm. This spatial resolution is comparable to the magnetic exchange length of transition metals, enabling us to probe monopole-monopole interactions. We find that the TMM and anti-TMM pairs are separated by 18.3±1.6 nm, while the TMM and TMM, anti-TMM and anti-TMM pairs are stabilized at comparatively longer distances of 36.1±2.4 nm and 43.1±2.0 nm, respectively. We also observe virtual TMMs created by magnetic voids in the meta-lattice. This work demonstrates that ferromagnetic meta lattices could be used as a platform to create and investigate the interactions and dynamics of TMMs. Furthermore, we expect that soft x-ray vector ptycho-tomography can be broadly applied to quantitatively image 3D vector fields in magnetic and anisotropic materials at the nanoscale.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Studies of Inviscid Flux Schemes for Acoustics and Turbulence Problems

The last two decades have witnessed tremendous growth in computational power, the development of computational fluid dynamics (CFD) codes which scale well over thousands of processors, and the refinement of unstructured grid-generation tools which facilitate rapid surface and volume gridding of complex geometries. Thus, engineering calculations of 10(exp 7) - 10(exp 8) finite-volume cells have become routine for some types of problems. Although the Reynolds Averaged Navier Stokes (RANS) approach to modeling turbulence is still in extensive and wide use, increasingly large-eddy simulation (LES) and hybrid RANS-LES approaches are being applied to resolve the largest scales of turbulence in many engineering problems. However, it has also become evident that LES places different requirements on the numerical approaches for both the spatial and temporal discretization of the Navier Stokes equations than does RANS. In particular, LES requires high time accuracy and minimal intrinsic numerical dispersion and dissipation over a wide spectral range. In this paper, the performance of both central-difference and upwind-biased spatial discretizations is examined for a one-dimensional acoustic standing wave problem, the Taylor-Green vortex problem, and the turbulent channel fl ow problem.

Morris, C. I.↗

Studies of Inviscid Flux Schemes for Acoustics and Turbulence Problems

The last two decades have witnessed tremendous growth in computational power, the development of computational fluid dynamics (CFD) codes which scale well over thousands of processors, and the refinement of unstructured grid-generation tools which facilitate rapid surface and volume gridding of complex geometries. Thus, engineering calculations of 10(exp 7) - 10(exp 8) finite-volume cells have become routine for some types of problems. Although the Reynolds Averaged Navier Stokes (RANS) approach to modeling turbulence is still in extensive and wide use, increasingly large-eddy simulation (LES) and hybrid RANS-LES approaches are being applied to resolve the largest scales of turbulence in many engineering problems. However, it has also become evident that LES places different requirements on the numerical approaches for both the spatial and temporal discretization of the Navier Stokes equations than does RANS. In particular, LES requires high time accuracy and minimal intrinsic numerical dispersion and dissipation over a wide spectral range. In this paper, the performance of both central-difference and upwind-biased spatial discretizations is examined for a one-dimensional acoustic standing wave problem, the Taylor-Green vortex problem, and the turbulent channel ow problem.

Morris, Christopher I.↗