SEARCH · Search NASA

Results for “Scale resolving simulations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS2001

Warm-Core Intensification Through Horizontal Eddy Heat Transports into the Eye

The mechanism for the formation and intensification of the hurricane warm core is not well understood. The generally accepted explanation is that the warm core forms as a result of gentle subsidence of air within the eye that warms as a result of adiabatic compression. Malkus suggested that this subsidence is part of a deep circulation in which air begins descent at high levels in the eye, acquires cyclonic angular momentum as it descends to lower levels, and then diverges at low levels, where it is entrained back into the eyewall. Inward mixing from the eyewall is hypothesized to force the subsidence and maintain the moisture and momentum budgets of the subsiding air. Willoughby suggested that air within the eye has remained so since it was first enclosed during the formation of the eyewall and that it subsides at most only a few kilometers rather than through the depth of the troposphere. He relates the subsidence to the low-level divergence and entrainment into the eyewall noted by Malkus, but suggests that shrinkage of the eye's volume is more than adequate to account for the air lost to the eyewall or converted to cloudy air by turbulent mixing across the eye boundary. Smith offered an alternative view of the subsidence forcing, suggesting that vertical motion in a mature hurricane eye is generated largely by imbalances between the downward vertical pressure gradient force and the upward buoyancy force. The vertical pressure gradient force is associated with the decay and/or radial spread of the tangential wind field with height at those levels were the winds are in approximate gradient wind balance. The rate of subsidence is just that required to warm the air sufficiently such that the buoyancy remains in close hydrostatic balance with an increasing vertical pressure gradient force. In this study, a very high-resolution simulation of Hurricane Bob using a cloud-resolving grid scale of 1.3 km is used to examine the heat budget within the storm with particular emphasis on the mechanisms for warming of the eye.

Braun, Scott A.↗

NASA NTRS2001

Cloud Microphysics Budget in the Tropical Deep Convective Regime

Cloud microphysics budgets in the tropical deep convective regime are analyzed based on a 2-D cloud resolving simulation. The model is forced by the large-scale vertical velocity and zonal wind and large-scale horizontal advections derived from TOGA COARE for a 20-day period. The role of cloud microphysics is first examined by analyzing mass-weighted mean heat budget and column-integrated moisture budget. Hourly budgets show that local changes of mass-weighted mean temperature and column-integrated moisture are mainly determined by the residuals between vertical thermal advection and latent heat of condensation and between vertical moisture advection and condensation respectively. Thus, atmospheric thermodynamics depends on how cloud microphysical processes are parameterized. Cloud microphysics budgets are then analyzed for raining conditions. For cloud-vapor exchange between cloud system and its embedded environment, rainfall and evaporation of raindrop are compensated by the condensation and deposition of supersaturated vapor. Inside the cloud system, the condensation of supersaturated vapor balances conversion from cloud water to raindrop, snow, and graupel through collection and accretion processes. The deposition of supersaturated vapor balances conversion from cloud ice to snow through conversion and riming processes. The conversion and riming of cloud ice and the accretion of cloud water balance conversion from snow to graupel through accretion process. Finally, the collection of cloud water and the melting of graupel increase raindrop to compensate the loss of raindrop due to rainfall and the evaporation of raindrop.

Li, Xiao-Fan↗

NASA NTRS2007

OMI Tropospheric NO2 from Lightning in Observed Convective Events

Lightning is responsible for an estimated 15 percent of total NO emissions, and is one of the most prominent sources in the upper troposphere. In this study, we present evidence of lightning-generated NO2 (LNO2) using data from the Ozone Monitoring Instrument (OMI), which has observed tropospheric NO2 since its launch in 2004. Although LNO2 has been also reported in previous satellite studies from the Global Ozone Monitoring Experiment (GOME) and SCIAMACHY, OMI is better suited for such measurements by virtue of its higher spatial resolution and daily global coverage. We will present data clearly showing the LNO2 signal in the OMI tropospheric NO2 product on two days over and downwind of specific convective systems in the US Midwest. Gridded monthly mean tropospheric NO 2 data are subtracted from the daily gridded data to obtain the presumed LNO2 signal. Observed cloud-to-ground (CG) lightning flashes from the National Lightning Detection Network (NLDN) were counted along middle and upper tropospheric back trajectories that were run from the regions containing the LNO2 signal. A vertically-weighted average number of upwind CG flashes was obtained using a profile of LNO(x) mass obtained from a series of midlatitude cloud-resolved storm chemistry simulations. The number of CG flashes was scaled up to total flashes (intracloud (IC) flashes plus CG) using a climatological IC/CG ratio. The number of moles of LNO(x) in the region considered was estimated by assuming that LNO2 is 30 percent of LNO(x). This value was divided by the number of upwind flashes to obtain an average estimate of the number of moles produced per flash. Results yield values in the range obtained through other estimation techniques (e.g., aircraft measurements, models). We will also present a similar analysis over northern Australia during the SCOUT-O3/ACTIVE field campaigns in November and December 2005, in which we will compare the OMI LNOx signals with aircraft observations from the storm anvils.

Pickering, Kenneth↗

NASA NTRS1993

On the large eddy simulation of turbulent flows in complex geometry

Application of the method of Large Eddy Simulation (LES) to a turbulent flow consists of three separate steps. First, a filtering operation is performed on the Navier-Stokes equations to remove the small spatial scales. The resulting equations that describe the space time evolution of the 'large eddies' contain the subgrid-scale (sgs) stress tensor that describes the effect of the unresolved small scales on the resolved scales. The second step is the replacement of the sgs stress tensor by some expression involving the large scales - this is the problem of 'subgrid-scale modeling'. The final step is the numerical simulation of the resulting 'closed' equations for the large scale fields on a grid small enough to resolve the smallest of the large eddies, but still much larger than the fine scale structures at the Kolmogorov length. In dividing a turbulent flow field into 'large' and 'small' eddies, one presumes that a cut-off length delta can be sensibly chosen such that all fluctuations on a scale larger than delta are 'large eddies' and the remainder constitute the 'small scale' fluctuations. Typically, delta would be a length scale characterizing the smallest structures of interest in the flow. In an inhomogeneous flow, the 'sensible choice' for delta may vary significantly over the flow domain. For example, in a wall bounded turbulent flow, most statistical averages of interest vary much more rapidly with position near the wall than far away from it. Further, there are dynamically important organized structures near the wall on a scale much smaller than the boundary layer thickness. Therefore, the minimum size of eddies that need to be resolved is smaller near the wall. In general, for the LES of inhomogeneous flows, the width of the filtering kernel delta must be considered to be a function of position. If a filtering operation with a nonuniform filter width is performed on the Navier-Stokes equations, one does not in general get the standard large eddy equations. The complication is caused by the fact that a filtering operation with a nonuniform filter width in general does not commute with the operation of differentiation. This is one of the issues that we have looked at in detail as it is basic to any attempt at applying LES to complex geometry flows. Our principal findings are summarized.

Ghosal, Sandip↗

NASA NTRS1993

A normal stress subgrid-scale eddy viscosity model in large eddy simulation

The Smagorinsky subgrid-scale eddy viscosity model (SGS-EVM) is commonly used in large eddy simulations (LES) to represent the effects of the unresolved scales on the resolved scales. This model is known to be limited because its constant must be optimized in different flows, and it must be modified with a damping function to account for near-wall effects. The recent dynamic model is designed to overcome these limitations but is compositionally intensive as compared to the traditional SGS-EVM. In a recent study using direct numerical simulation data, Horiuti has shown that these drawbacks are due mainly to the use of an improper velocity scale in the SGS-EVM. He also proposed the use of the subgrid-scale normal stress as a new velocity scale that was inspired by a high-order anisotropic representation model. The testing of Horiuti, however, was conducted using DNS data from a low Reynolds number channel flow simulation. It was felt that further testing at higher Reynolds numbers and also using different flows (other than wall-bounded shear flows) were necessary steps needed to establish the validity of the new model. This is the primary motivation of the present study. The objective is to test the new model using DNS databases of high Reynolds number channel and fully developed turbulent mixing layer flows. The use of both channel (wall-bounded) and mixing layer flows is important for the development of accurate LES models because these two flows encompass many characteristic features of complex turbulent flows.

Horiuti, K.↗

NASA NTRS1992

Large eddy simulation of compressible turbulent channel flows

The development of turbulence within rocket propulsion chamber flows remains a difficult problem to predict. Within solid propellant rockets, the flow can exhibit multiple regions of transition to turbulence, and is susceptible to various modes of aero-acoustic interaction. The objective is to formulate, develop, and validate, a large eddy simulation (LES) method for compressible channel flows. The LES method separates the resolvable scale motions from the unresolvable scales by applying a spatial filter to the compressible Navier-Stokes equations. The equations are solved numerically using a modified four-step Runge-Kutta procedure in time and second or fourth order finite differences in space. The current simulations have validated the algorithm in the incompressible limit and have demonstrated the ability of the method to simulate high subsonic Mach number flows.

Beddini, Robert A.↗

NASA NTRS1992

Search for subgrid scale parameterization by projection pursuit regression

The dependence of subgrid-scale stresses on variables of the resolved field is studied using direct numerical simulations of isotropic turbulence, homogeneous shear flow, and channel flow. The projection pursuit algorithm, a promising new regression tool for high-dimensional data, is used to systematically search through a large collection of resolved variables, such as components of the strain rate, vorticity, velocity gradients at neighboring grid points, etc. For the case of isotropic turbulence, the search algorithm recovers the linear dependence on the rate of strain (which is necessary to transfer energy to subgrid scales) but is unable to determine any other more complex relationship. For shear flows, however, new systematic relations beyond eddy viscosity are found. For the homogeneous shear flow, the results suggest that products of the mean rotation rate tensor with both the fluctuating strain rate and fluctuating rotation rate tensors are important quantities in parameterizing the subgrid-scale stresses. A model incorporating these terms is proposed. When evaluated with direct numerical simulation data, this model significantly increases the correlation between the modeled and exact stresses, as compared with the Smagorinsky model. In the case of channel flow, the stresses are found to correlate with products of the fluctuating strain and rotation rate tensors. The mean rates of rotation or strain do not appear to be important in this case, and the model determined for homogeneous shear flow does not perform well when tested with channel flow data. Many questions remain about the physical mechanisms underlying these findings, about possible Reynolds number dependence, and, given the low level of correlations, about their impact on modeling. Nevertheless, demonstration of the existence of causal relations between sgs stresses and large-scale characteristics of turbulent shear flows, in addition to those necessary for energy transfer, provides important insight into the relation between scales in turbulent flows.

Meneveau, C.↗

NASA NTRS2009

Combining Satellite Microwave Radiometer and Radar Observations to Estimate Atmospheric Latent Heating Profiles

In this study, satellite passive microwave sensor observations from the TRMM Microwave Imager (TMI) are utilized to make estimates of latent + eddy sensible heating rates (Q1-QR) in regions of precipitation. The TMI heating algorithm (TRAIN) is calibrated, or "trained" using relatively accurate estimates of heating based upon spaceborne Precipitation Radar (PR) observations collocated with the TMI observations over a one-month period. The heating estimation technique is based upon a previously described Bayesian methodology, but with improvements in supporting cloud-resolving model simulations, an adjustment of precipitation echo tops to compensate for model biases, and a separate scaling of convective and stratiform heating components that leads to an approximate balance between estimated vertically-integrated condensation and surface precipitation. Estimates of Q1-QR from TMI compare favorably with the PR training estimates and show only modest sensitivity to the cloud-resolving model simulations of heating used to construct the training data. Moreover, the net condensation in the corresponding annual mean satellite latent heating profile is within a few percent of the annual mean surface precipitation rate over the tropical and subtropical oceans where the algorithm is applied. Comparisons of Q1 produced by combining TMI Q1-QR with independently derived estimates of QR show reasonable agreement with rawinsonde-based analyses of Q1 from two field campaigns, although the satellite estimates exhibit heating profile structure with sharper and more intense heating peaks than the rawinsonde estimates. 2

Grecu, Mircea↗

NASA NTRS2000

Ion Dynamics Model for Collisionless Radio Frequency Sheaths

Full scale reactor model based on fluid equations is widely used to analyze high density plasma reactors. It is well known that the submillimeter scale sheath in front of a biased electrode supporting the wafer is difficult to resolve in numerical simulations, and the common practice is to use results for electric field from some form of analytical sheath model as boundary conditions for full scale reactor simulation. There are several sheath models in the literature ranging from Child's law to a recent unified sheath model [P. A. Miller and M. E. Riley, J. Appl. Phys. 82, 3689 (1997)l. In the present work, the cold ion fluid equations in the radio frequency sheath are solved numerically to show that the spatiotemporal variation of ion flux inside the sheath, commonly ignored in analytical models, is important in determining the electric field and ion energy at the electrode. Consequently, a semianalytical model that includes the spatiotemporal variation of ion flux is developed for use as boundary condition in reactor simulations. This semianalytical model is shown to yield results for sheath properties in close agreement with numerical solutions.

Bose, Deepak↗

NASA NTRSDate not supplied

Aeroacoustic Study of a Subscale Large Civil Transport (STAR) Model – Part 2: Validation of Simulated Results

Aeroacoustic measurements of the 26%-scale, semispan Boeing 777-200 Subsonic Transport Aeroacoustic Research (STAR) model tested in the NASA Ames Research Center 40- by 80-foot wind tunnel were used to ascertain the efficacy of high-fidelity simulations to accurately predict noise from the landing gear of large commercial transports. The simulations, conducted with the lattice Boltzmann solver PowerFLOW®, used a digital replica of the STAR model with or without main landing gear deployed and slats and flaps set to their highest deflection angles to represent aircraft during landing. The computations were performed at a Mach number of 0.21, Reynolds number of 8.2 million based on the model mean aerodynamic chord, and other conditions prevalent during the STAR model test. Measured and computed surface pressures were in very good agreement at most port locations on the model, as were global force coefficients, indicating that the simulations captured the impact of main gear deployment on inboard flap loading. Noise sources produced by the main landing gear and high-lift devices were determined via source localization maps generated with CLEAN from synthetic and experimental data. In general, very good agreement between predicted and measured acoustic source location and relative strength was observed in the maps. Comparisons of far-field noise spectra obtained from the CLEAN deconvolution maps showed remarkable agreement between synthetic and experimental broadband noise at low and medium frequencies. Main landing gear sources for model-scale frequencies above 7,000 Hz could not be resolved with the spatial resolution used during the simulations.

airframe noise↗

NASA NTRS2001

Precipitation Efficiency in the Tropical Deep Convective Regime

Precipitation efficiency in the tropical deep convective regime is analyzed based on a 2-D cloud resolving simulation. The cloud resolving model is forced by the large-scale vertical velocity and zonal wind and large-scale horizontal advections derived from TOGA COARE for a 20-day period. Precipitation efficiency may be defined as a ratio of surface rain rate to sum of surface evaporation and moisture convergence (LSPE) or a ratio of surface rain rate to sum of condensation and deposition rates of supersaturated vapor (CMPE). Moisture budget shows that the atmosphere is moistened (dryed) when the LSPE is less (more) than 100 %. The LSPE could be larger than 100 % for strong convection. This indicates that the drying processes should be included in cumulus parameterization to avoid moisture bias. Statistical analysis shows that the sum of the condensation and deposition rates is bout 80 % of the sum of the surface evaporation rate and moisture convergence, which ads to proportional relation between the two efficiencies when both efficiencies are less han 100 %. The CMPE increases with increasing mass-weighted mean temperature and creasing surface rain rate. This suggests that precipitation is more efficient for warm environment and strong convection. Approximate balance of rates among the condensation, deposition, rain, and the raindrop evaporation is used to derive an analytical solution of the CMPE.

Li, Xiaofan↗

NASA NTRS2012

TRMM Latent Heating Retrieval and Comparisons with Field Campaigns and Large-Scale Analyses

Rainfall production is a fundamental process within the Earth's hydrological cycle because it represents both a principal forcing term in surface water budgets, and its energetics corollary, latent heating (LH), is one of the principal sources of atmospheric diabatic heating. Latent heat release itself is a consequence of phase changes between the vapor, liquid, and frozen states of water. The vertical distribution of LH has a strong influence on the atmosphere, controlling large-scale tropical circulations, exciting and modulating tropical waves, maintaining the intensities of tropical cyclones, and even providing the energetics of midlatitude cyclones and other mobile midlatitude weather systems. Moreover, the processes associated with LH result in significant non-linear changes in atmospheric radiation through the creation, dissipation and modulation of clouds and precipitation. Yanai et al. (1973) utilized the meteorological data collected from a sounding network to present a pioneering work on thermodynamic budgets, which are referred to as the apparent heat source (Q1) and apparent moisture sink (Q2). Yanai's paper motivated the development of satellite-based LH algorithms and provided a theoretical background for imposing large-scale advective forcing into cloud-resolving models (CRMs). These CRM-simulated LH and Q1 data have been used to generate the look-up tables used in LH algorithms. This paper examines the retrieval, validation, and application of LH estimates based on rain rate quantities acquired from the Tropical Rainfall Measuring Mission satellite (TRMM). TRMM was launched in November 1997 as a joint enterprise between the American and Japanese space agencies -- with overriding goals of providing accurate four-dimensional estimates of rainfall and LH over the global Tropics and subtropics equatorward of 35o. Other literature has acknowledged the achievement of the first goal of obtaining an accurate rainfall climatology. This paper describes the second major goal of obtaining credible LH estimates as well as their applications within TRMM's zone of coverage, the standard TRMM LH products, and areas for further improvement.

latent heating↗

NASA NTRS2006

Unsteady Flow Simulation of High-speed Turbopumps

Computation of high-speed hydrodynamics requires high-fidelity simulation to resolve flow features involving transient flow, cavitation, tip vortex and multiple scales of unsteady fluctuations. One example of this type in aerospace is related to liquid-fueled rocket turbopump. Rocket turbopumps operate under severe conditions at very high rotational speeds typically at thousands of rpm. For example, the Shuttle orbiter low-pressure-fuel-turbopump creates transient flow features associated with reverse flows, tip clearance effects, secondary flows, vortex shedding, junction flows, and cavitation effects. Flow unsteadiness originating from the orbiter Low-Pressure-Fuel-Turbopump (LPFTP) inducer is one of the major contributors to the high frequency cyclic loading that results in high cycle fatigue damage to the flow liners just upstream of the LPFTP. The reverse flow generated at the tip of the inducer blades travels upstream and interacts with the bellows cavity. Simulation procedure for this type high-speed hydrodynamic problems requires a method for quantifying multi-scale and multi-phase flow as well as an efficient high-end computing strategy. The current paper presents a high-fidelity computational procedure for unsteady hydrodynamic problems using a high-speed liquid-fueled rocket turbopump.

Kiris, Cetin C.↗

NASA NTRS1991

Large-eddy simulation of compressible turbulence

The increase in the range of length scales with increasing Reynolds number limits the direct simulation of turbulent flows to relatively simple geometries and low Reynolds numbers. However, since most flows of engineering interest occur at much higher Reynolds number than is currently within the capabilities of full simulation, prediction of these flow fields can only be obtained by solving some suitably-averaged set of governing equations. In the traditional Reynolds-averaged approach, the Navier-Stokes equations are averaged over time. This in turn yields correlations between various turbulence fluctuations. It is these terms, e.g. the Reynolds stresses, for which a turbulence model must be derived. Turbulence modeling of incompressible flows has received a great amount of attention in the literature. An area of research that has received comparatively less attention is the modeling of compressible turbulent flows. An approach to simulating compressible turbulence at high Reynolds numbers is through the use of Large-Eddy Simulation (LES). In LES the dependent variables are decomposed into a large-scale (resolved) component and a sub-grid scale component. It is the small-scale components of the velocity field which are presumably more homogeneous than the large scales and, therefore, more easily modeled. Thus, it seems plausible that simpler models, which should be more universal in character than those employed in second-order closure schemes, may be developed for LES of compressible turbulence. The objective of the present research, therefore, is to explore models for the Large-Eddy Simulation of compressible turbulent flows. Given the recent successes of Zeman in second order closure modeling of compressible turbulence, model development was guided by principals employed in second-order closures.

Squires, Kyle D.↗

NASA NTRS1995

Numerical simulation of filling a magnetic flux tube with a cold plasma: Anomalous plasma effects

Large-scale models of plasmaspheric refilling have revealed that during the early stage of the refilling counterstreaming ion beams are a common feature. However, the instability of such ion beams and its effect on refilling remain unexplored. In order to learn the basic effects of ion beam instabilities on refilling, we have performed numerical simulations of the refilling of an artificial magnetic flux tube. (The shape and size of the tube are assumed so that the essential features of the refilling problem are kept in the simulation and at the same time the small scale processes driven by the ion beams are sufficiently resolved.) We have also studied the effect of commonly found equatorially trapped warm and/or hot plasma on the filling of a flux tube with a cold plasma. Three types of simulation runs have been performed.

Singh, Nagendra↗

NASA NTRS1994

Multifractal characterizations of nonstationarity and intermittency in geophysical fields: Observed, retrieved, or simulated

Geophysical data rarely show any smoothness at any scale, and this often makes comparison with theoretical model output difficult. However, highly fluctuating signals and fractal structures are typical of open dissipative systems with nonlinear dynamics, the focus of most geophysical research. High levels of variability are excited over a large range of scales by the combined actions of external forcing and internal instability. At very small scales we expect geophysical fields to be smooth, but these are rarely resolved with available instrumentation or simulation tools; nondifferentiable and even discontinuous models are therefore in order. We need methods of statistically analyzing geophysical data, whether measured in situ, remotely sensed or even generated by a computer model, that are adapted to these characteristics. An important preliminary task is to define statistically stationary features in generally nonstationary signals. We first discuss a simple criterion for stationarity in finite data streams that exhibit power law energy spectra and then, guided by developments in turbulence studies, we advocate the use of two ways of analyzing the scale dependence of statistical information: singular measures and qth order structure functions. In nonstationary situations, the approach based on singular measures seeks power law behavior in integrals over all possible scales of a nonnegative stationary field derived from the data, leading to a characterization of the intermittency in this (gradient-related) field. In contrast, the approach based on structure functions uses the signal itself, seeking power laws for the statistical moments of absolute increments over arbitrarily large scales, leading to a characterization of the prevailing nonstationarity in both quantitative and qualitative terms. We explain graphically, step by step, both multifractal statistics which are largely complementary to each other. The geometrical manifestations of nonstationarity and intermittency, 'roughness' and 'sparseness', respectively, are illustrated and the associated analytical (differentiability and continuity) properties are discussed. As an example, the two techniques are applied to a series of recent measurements of liquid water distributions inside marine stratocumulus decks; these are found to be multifractal over scales ranging from approximately 60 m to approximately 60 km. Finally, we define the 'mean multifractal plane' and show it to be a simple yet comprehensive tool with many applications including data intercomparison, (dynamical or stochastic) model and retrieval validations.

Davis, Anthony↗

NASA NTRS1979

Computational aerodynamics development and outlook /Dryden Lecture in Research for 1979/

Some past developments and current examples of computational aerodynamics are briefly reviewed. An assessment is made of the requirements on future computer memory and speed imposed by advanced numerical simulations, giving emphasis to the Reynolds averaged Navier-Stokes equations and to turbulent eddy simulations. Experimental scales of turbulence structure are used to determine the mesh spacings required to adequately resolve turbulent energy and shear. Assessment also is made of the changing market environment for developing future large computers, and of the projections of micro-electronics memory and logic technology that affect future computer capability. From the two assessments, estimates are formed of the future time scale in which various advanced types of aerodynamic flow simulations could become feasible. Areas of research judged especially relevant to future developments are noted.

Chapman, D. R.↗

NASA NTRS2013

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