Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “advection”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Simulation of the modern arctic climate by the NCAR CCM1

The National Center of Atmospheric Research (NCAR) Community Climate Model Version 1 (CCM1's) simulation of the modern arctic climate is evaluated by comparing a five-year seasonal cycle simulation with the European Center for Medium-Range Weather Forecasts (ECMWF) global analyses. The sea level pressure (SLP), storm tracks, vertical cross section of height, 500-hPa height, total energy budget, and moisture budget are analyzed to investigate the biases in the simulated arctic climate. The results show that the model simulates anomalously low SLP, too much storm activity, and anomalously strong baroclinicity to the west of Greenland and vice versa to the east of Greenland. This bias is mainly attributed to the model's topographic representation of Greenland. First, the broadened Greenland topography in the model distorts the path of cyclone waves over the North Atlantic Ocean. Second, the model oversimulates the ridge over Greenland, which intensifies its blocking effect and steers the cyclone waves clockwise around it and hence produces an artificial circum-Greenland trough. These biases are significantly alleviated when the horizontal resolution increases to T42. Over the Arctic basin, the model simulates large amounts of low-level (stratus) clouds in winter and almost no stratus in summer, which is opposite to the observations. This bias is mainly due to the location of the simulated SLP features and the negative anomaly of storm activity, which prevent the transport of moisture into this region during summer but favor this transport in winter. The moisture budget analysis shows that the model's net annual precipitation (P-E) between 70 deg N and the North Pole is 6.6 times larger than the observations and the model transports six times more moisture into this region. The bias in the advection term is attributed to the positive moisture fixer scheme and the distorted flow pattern. However, the excessive moisture transport into the Arctic basin does not solely result from the advection term. The contribution by the moisture fixer is as large as from advection. By contrast, the semi-Lagrangian transport scheme used in the CCM2 significantly improves the moisture simulation for this region; however, globally the error is as serious as for the positive moisture fixer scheme. Finally, because the model has such serious problems in simulating the present arctic climate, its simulations of past and future climate change for this region are questionable.

Bromwich, David H.↗
NASA NTRS1994

Surface gravity waves at equilibrium with a steady wind

Observations of wave fields' spatial evolution and of gravity wave spectra S(sub omega) are analyzed on the basis of the data reported by several research groups as well as on a 2-year data set of wind and wave measurements by stationary National Oceanic and Atmospheric Administration buoys near the Hawaiian Islands. We seek to clarify the role of the wave energy advection (with the wave group velocity) in the overall energy balance. This advective transfer appears to be no less important than the local (breaking wave induced) dissipation as a factor of wind-wave equilibrium. The advection is found to manifest itself in the shape of wave spectra by reducing the rate at which the spectra density of the wave energy, S(sub omega approximately omega (exp -p), falls off as the frequency increases away from the spectra peak. This and other conclusions are derived by comparing the field observations with theoretical predictions of the weak turbulence theory for a spatially inhomogeneous, statistically stationary, wave field. The observations also indicate that the typical wave age xi = C(sub 0)/U in the open ocean is much greater than the limiting value 1.2 attributed to the 'fully developed sea.' Although the observed spectra can be approximated by a power law with a single 'effective' exponent, this apparent exponent, p, is found to depend on the wave age. At high xi and at frequencies below the generation range, -p tends to -3 rather than the value of -11/3 predicted by the Zakharov-Zaslavskii theory. This deviation is interpreted as pointing to a nonconservative nature of the inverse cascade, the latter including a leakage of energy to low-wavenumber modes. Dependence of the overall effective exponent on xi is shown to be responsible for variation in the coefficients b, B, c, C appearing in empirical fetch laws, such as xi = C bar-x (exp c) and e = B bar-x (exp B), where bar-x and e are the dimensionless fetch and wave energy, respectively.

Glazman, Roman E.↗
NASA NTRS1994

Synoptic interpretation of measurements from HALOE

The degree to which the Southern Hemisphere polar vortex is isolated against horizontal (isentropic) mixing is investigated using data from the Halogen Occultation Experiment (HALOE), U.K. Meteorological Office (UKMO) potential vorticity (PV), and contour advection diagnostics. Measurements of methane and water vapor taken by HALOE during a disturbed period in the Southern Hemisphere springtime (21 September - 15 October 1992) are interpreted in light of the prevailing synoptic meteorology. Daily fields of winds and PV are shown to be essential in the interpretation of the data. A climatological high pressure region is responsible for a distorted vortex, and a substantial 'vortex stripping' event is present, associated with the early stages of vortex breakdown. This leads to significant temporal, zonal, and altitudinal variations in the distribution of tracers. The authors point out the difficulties this presents for the interpretation of solar occultation data, especially with regard to the use of zonal average time series. Longitude-height methane distributions from two days during the period are examined. Both days show substantial variations in abundance around a latitude circle. In particular, the authors investigate HALOE measurements at 77 deg S on 15 October 1992, which indicate an abundance of methane in the height region 600-2000 K (approximately 30-1 mb) that is more typical of midlatitude air. Similar distributions, observed in the 1991 HALOE data, have previously been interpreted as evidence for the penetration of midlatitude air into the vortex. Gradients of potential vorticity and contour advection diagnositcs are employed to examine whether the UKMO winds are consistent with this hypothesis in 1992. Although midlatitude air is able to penetrate poleward of the main jet core by advection processes alone, an essentially intact inner core of vortex air remains, which does not mix to any great extent with air from lower latitudes. The authors show that the high-latitude HALOE abundances that are typical of midlatitude air were observed in a region of extensive filamentation and mixing, rather than within the inner, more isolated, core.

Bithell, M.↗
NASA NTRS1994

New production in the warm waters of the tropical Pacific Ocean

The average depth-integrated rate of new production in the tropical Pacific Ocean was estimated from a calculation of horizontal and vertical nitrate balance over the region enclosed by the climatological 26 C isotherm. The net turbulent flux of nitrate into the region was computed in terms of the climatological net surface heat flux and the nitrate-temperature relationship at the base of the 26 C isotherm. The net advective transport of nitrate into the region was estimated using the mean nitrate distribution obtained from the analysis of historical data and previous results of a general circulation model of the tropical Pacific. The rate of new production resulting from vertical turbulent fluxes of nitrate was found to be similar in magnitude to that due to advective transport. Most (about 75%) of the advective input of nitrate was due to the horizontal transport of nutrient-rich water from the eastern equatorial region rather than from equatorial upwelling. An average rate of new production of 14.5 - 16 g C/sq m/yr was found for the warm waters of the tropical Pacific region. These values are in good agreement with previous estimates for this region and are almost five times less than is estimated for the eastern equatorial Pacific, where most of the nutrient upwelling occurs.

Pena, M. Angelica↗
NASA NTRS1993

Water Surface Currents, Short Gravity-Capillary Waves and Radar Backscatter

Despite their importance for air-sea interaction and microwave remote sensing of the ocean surface, intrinsic properties of short gravity-capillary waves are not well established. This is largely due to water surface currents and their effects on the direct measurements of wave parameters conducted at a fixed point. Frequencies of small scale waves propagating on a surface which itself is in motion, are subject to Doppler shifts. Hence, the high frequency tail of the wave spectra obtained from such temporal observations is smeared. Conversion of this smeared measured-frequency spectra to intrinsic-frequency (or wavenumber) spectra requires corrections for the Doppler shifts. Such attempts in the past have not been very successful in particular when field data were used. This becomes evident if the amplitude modulation of short waves by underlying long waves is considered. Microwave radar studies show that the amplitude of a short wave component attains its maximum value near the crests and its minimum in the troughs of the long waves. Doppler-shifted wave data yield similar results but much larger in modulation magnitude, as expected. In general, Doppler shift corrections reduce the modulation magnitude. Overcorrection may result in a negligible modulation or even in a strong modulation with the maximum amplitude in the wave troughs. The latter situation is clearly contradictory to our visual observations as well as the radar results and imply that the advection by currents is overestimated. In this study, a differential-advection approach is used in which small scale waves are advected by the currents evaluated not at the free surface, but at a depth proportional to their wavelengths. Applicability of this approach is verified by the excellent agreement in phase and magnitude of short-wave modulation between results based on radar and on wave-gauge measurements conducted on a lake.

Atakturk, Serhad S.↗
NASA NTRS1995

A diagnosis of the development of a winter anticyclone over North America

This paper examines the 48-h life cycle of a winter anticyclone occurring over North America from 18 to 20 January 1979 using Goddard Laboratory for Atmospheres FGGE level 3b (SOP 1) global analyses on a 4 deg latitude by 5 deg longitude grid. Applying the relatively new methodology of the Zwack-Okossi equation, results show that anticyclonic vorticity advection and cold-air advection acted to develop the anticyclone, while adiabatic warming in the descending air opposed development. Other forcing processes made only small contributions to anticyclone changes. Vertical profiles of the development quantities reveal that vorticity and temperature advections, as well as the adiabatic warming, maximized in the 200-300-mb layer.

King, Melinda L.↗
NASA NTRS1994

Convection in the Physical Vapor Transport Process: Thermosolutal - Part 2

We consider the effect of an inert gas on the diffusive-convective physical vapor transport process. We investigate the case when the temperature gradient is stabilizing and the concentration gradient is destabilizing for a wide parametric range. When an inert gas is present, the thermal and solutal convection oppose each other. The solutal field is destabilizing while the thermal field and the advective-diffusive flux stabilize the flow field. When the pressure of the inert component is increased, the stabilizing effect of the advective-diffusive flux is decreased. The intensity of convection as well as the oscillatory transient time increases. Below, the critical Rayleigh number, the nonlinear dynamics of the flow field show an oscillatory approach to steady state. For parametric values in the neighborhood of the critical Rayleigh number, the flow field undergoes a chaotic transient which settles to a periodic state. The asymptotic state of the flow field shows that growth and amalgamation of cells yields an overturning motion which results in an asymmetric cellular structure. The low gravity environment yields the stabilizing advective-diffusive flow which results in uniform temperature and concentration gradients near the crystal interface.

Duval, Walter M. B.↗
NASA NTRS1996

An Atmospheric Tape Recorder: The Imprint of Tropical Tropopause Temperatures on Stratospheric Water Vapor

We describe observations of tropical stratospheric water vapor q that show clear evidence of large-scale upward advection of the signal from annual fluctuations in the effective 'entry mixing ratio' q(sub E) of air entering the tropical stratosphere. In other words, air is 'marked,' on emergence above the highest cloud tops, like a signal recorded on an upward moving magnetic tape. We define q(sub E) as the mean water vapor mixing ratio, at the tropical tropopause, of air that will subsequently rise and enter the stratospheric 'overworld' at about 400 K. The observations show a systematic phase lag, increasing with altitude, between the annual cycle in q(sub E) and the annual cycle in q at higher altitudes. The observed phase lag agrees with the phase lag calculated assuming advection by the transformed Eulerian-mean vertical velocity of a q(sub E) crudely estimated from 100-hPa temperatures, which we use as a convenient proxy for tropopause temperatures. The phase agreement confirms the overall robustness of the calculation and strongly supports the tape recorder hypothesis. Establishing a quantitative link between q(sub E) and observed tropopause temperatures, however, proves difficult because the process of marking the tape depends subtly on both small- and large-scale processes. The tape speed, or large-scale upward advection speed, has a substantial annual variation and a smaller variation due to the quasi-biennial oscillation, which delays or accelerates the arrival of the signal by a month or two in the middle stratosphere. As the tape moves upward, the signal is attenuated with an e-folding time of about 7 to 9 months between 100 and 50 hPa and about 15 to 18 months between 50 and 20 hPa, constraining possible orders of magnitude both of vertical diffusion K(sub z) and of rates of mixing in from the extratropics. For instance, if there were no mixing in, then K(sub z) would be in the range 0.03-0.09 m(exp 2)/s; this is an upper bound on K(sub z).

Mote, Philip W.↗
NASA NTRS1997

On the Behavior of Velocity Fluctuations in Rapidly Rotating Flows

The behavior of velocity fluctuations subjected to rapid rotation is examined. The rapid rotation considered is any arbitrary combination of two basic forms of rotation, reference frame rotation and mean flow rotation. It is recognized that the two types of rotating flows differ in the manner in which the fluctuating fields are advected. The first category is comprised of flows in rotating systems of which synoptic scale geophysical flows are a good example. In this class of flows the fluctuating velocity field advects and rotates with the mean flow. In the rapid rotation limit, the Taylor-Proudman theorem describes the behavior of this class of fluctuations. Velocity fluctuations that are advected without rotation by the mean flow constitute the second category which includes vortical flows of aerodynamic interest. The Taylor-Proudman theorem is not pertinent to I his class flows and a new result appropriate to this second category of fluctuations is derived. The present development demonstrates that the fluctuating velocity fields are rendered two-dimensional and horizontally non-divergent in the limit of any large combination of reference frame rotation and mean-flow rotation. The concommitant 'geostrophic' balance of the momentum equation is, however, dependent upon the form of rapid rotation. It is also demonstrated that the evolution equations of a two-dimensional fluctuating velocity fields are frame-indifferent with any imposed mean-flow rotation. The analyses and results of this paper highlight many fundamental aspects of rotating flows and have important consequences for their turbulence closures in inertial and non-inertial frames.

Girimaji, S. S.↗
NASA NTRS1997

Spectral and Temporal Properties of Galactic Black Hole Systems

Kusunose, Mineshige & Yamada (1996; hereafter KMY) extended the model of Kusunose & Mineshige (1995) to the Galactic black hole candidates by considering nonthermal electron injection with gamma(EQ\0(,\s\up2(less than),\s\do-l(_))) 10. The effects of pair escape and advection on the disk structure and general relativistic effects on the emission spectrum were also examined. They found that the energy spectral index (alpha)(sub x) of the power law X-rays is about-0.8 and-2.0 when 1(sub soft)/1 = 0.2 and 2, respectively, where 1(sub soft)/1 is the ratio of the compactness of the injected soft photons to that of the gravitational energy. The power law index was found to be nearly independent of the mass accretion which is consistent with the observed luminosity independence. The model with small 1(sub soft)/1 (less than 1) shows promise for explaining the low state observed in Galactic black hole candidates. Model fits were provided for GX339-4 and Cyg X- 1 data from COMPTEL and OSSE on the Compton Gamma Ray Observatory. The difference in emission spectra between thermal disks and the model of KMY appears only in the energy range greater than 100 keV. Li, Kusunose and Liang (1996) studied stochastic particle acceleration to produce nonthermal particle distributions which then were used in the model of Kusunose & Mineshige (1995) to model the spectrum above 1 Mev from GBHC's. Under certain conditions, stochastic electron acceleration overcame Coulomb and Compton losses resulting in a suprathermal electron population. Good fits were obtained by COMPTEL and OSSE observations of Cyg X-1 and GRO J0422+23. Kusunose & Mineshige (1996a) examined the role of electron-positron pairs in advection-dominated disks. They found that the results for advection-dominated disks without pairs are not qualitatively changed by including pairs. Summaries of work sponsored by this grant are given in Wheeler, Kim, Moscoso, Kusunose & Mineshige (1996) and Kusunose (1996) Work was also done on developing a model for an e(+-) pair wind from the inner disk region of a black hole. The model consists of three zones: a pair production/annihilation zone at the base of the wind, a pair annihilation zone slightly further out from the disk and a wind zone extending to infinity where no annihilation occurs. The model assumes an input X-ray / gamma-ray power-law spectrum as a function of photon energy and angular distribution. Pairs are created in the pair production/annihilation zone via photon-photon, photonparticle and particle-particle collisions. The bulk velocity of the pairs is obtained from the momentum component perpendicular to the disk taking into account the radiation pressure acceleration of the pairs. Energy balance in the pair production/annihilation zone is calculated by assuming that the momentum component parallel to the disk is thermalized and by taking into account heating/cooling via anisotropic thermal Compton scattering along with the cooling processes of bremsstrahlung and synchrotron radiation. The pair density is also calculated. The outer boundary of the pair production/annihilation zone is where the optical depth for photon-photon pair produciton is unity. Beyond this point only pair annihilation will occur in the pair annihilation zone. When the pair density becomes small, pair annihilation will be negligible and the pairs will flow freely to infinity. In this model we have found that the X-ray power-law spectral index, alpha(sub x) is the primary parameter which determines the density, temperature and velocity of the pair production/annihilation zone (Moscoso, Kusunose & Wheeler 1996).

Wheeler, J. Craig↗
NASA NTRS1998

Bulk-Flow Analysis of Hybrid Thrust Bearings for Advanced Cryogenic Turbopumps

A bulk-flow analysis and computer program for prediction of the static load performance and dynamic force coefficients of angled injection, orifice-compensated hydrostatic/hydrodynamic thrust bearings have been completed. The product of the research is an efficient computational tool for the design of high-speed thrust bearings for cryogenic fluid turbopumps. The study addresses the needs of a growing technology that requires of reliable fluid film bearings to provide the maximum operating life with optimum controllable rotordynamic characteristics at the lowest cost. The motion of a cryogenic fluid on the thin film lands of a thrust bearing is governed by a set of bulk-flow mass and momentum conservation and energy transport equations. Mass flow conservation and a simple model for momentum transport within the hydrostatic bearing recesses are also accounted for. The bulk-flow model includes flow turbulence with fluid inertia advection, Coriolis and centrifugal acceleration effects on the bearing recesses and film lands. The cryogenic fluid properties are obtained from realistic thermophysical equations of state. Turbulent bulk-flow shear parameters are based on Hirs' model with Moody's friction factor equations allowing a simple simulation for machined bearing surface roughness. A perturbation analysis leads to zeroth-order nonlinear equations governing the fluid flow for the thrust bearing operating at a static equilibrium position, and first-order linear equations describing the perturbed fluid flow for small amplitude shaft motions in the axial direction. Numerical solution to the zeroth-order flow field equations renders the bearing flow rate, thrust load, drag torque and power dissipation. Solution to the first-order equations determines the axial stiffness, damping and inertia force coefficients. The computational method uses well established algorithms and generic subprograms available from prior developments. The Fortran9O computer program hydrothrust runs on a Windows 95/NT personal computer. The program, help files and examples are licensed by Texas A&M University Technology License Office. The study of the static and dynamic performance of two hydrostatic/hydrodynamic bearings demonstrates the importance of centrifugal and advection fluid inertia effects for operation at high rotational speeds. The first example considers a conceptual hydrostatic thrust bearing for an advanced liquid hydrogen turbopump operating at 170,000 rpm. The large axial stiffness and damping coefficients of the bearing should provide accurate control and axial positioning of the turbopump and also allow for unshrouded impellers, therefore increasing the overall pump efficiency. The second bearing uses a refrigerant R134a, and its application in oil-free air conditioning compressors is of great technological importance and commercial value. The computed predictions reveal that the LH2 bearing load capacity and flow rate increase with the recess pressure (i.e. increasing orifice diameters). The bearing axial stiffness has a maximum for a recess pressure rati of approx. 0.55. while the axial damping coefficient decreases as the recess pressure ratio increases. The computer results from three flow models are compared. These models are a) inertialess, b) fluid inertia at recess edges only, and c) full fluid inertia at both recess edges and film lands. The full inertia model shows the lowest flow rates, axial load capacity and stiffness coefficient but on the other hand renders the largest damping coefficients and inertia coefficients. The most important findings are related to the reduction of the outflow through the inner radius and the appearance of subambient pressures. The performance of the refrigerant hybrid thrust bearing is evaluated at two operating speeds and pressure drops. The computed results are presented in dimensionless form to evidence consistent trends in the bearing performance characteristics. As the applied axial load increases, the bearing film thickness and flow rate decrease while the recess pressure increases. The axial stiffness coefficient shows a maximum for a certain intermediate load while the damping coefficient steadily increases. The computed results evidence the paramount of centrifugal fluid inertia at low recess pressures (i.e. low loads), and where there is actually an inflow through the bearing inner diameter, accompanied by subambient pressures just downstream of the bearing recess edge. These results are solely due to centrifugal fluid inertia and advection transport effects. Recommendations include the extension of the computer program to handle flexure pivot tilting pad hybrid bearings and the ability to calculate moment coefficients for shaft angular misalignments.

SanAndres, Luis↗
NASA NTRS1999

Physically-Derived Dynamical Cores in Atmospheric General Circulation Models

The algorithm chosen to represent the advection in atmospheric models is often used as the primary attribute to classify the model. Meteorological models are generally classified as spectral or grid point, with the term grid point implying discretization using finite differences. These traditional approaches have a number of shortcomings that render them non-physical. That is, they provide approximate solutions to the conservation equations that do not obey the fundamental laws of physics. The most commonly discussed shortcomings are overshoots and undershoots which manifest themselves most overtly in the constituent continuity equation. For this reason many climate models have special algorithms to model water vapor advection. This talk focuses on the development of an atmospheric general circulation model which uses a consistent physically-based advection algorithm in all aspects of the model formulation. The shallow-water model is generalized to three dimensions and combined with the physics parameterizations of NCAR's Community Climate Model. The scientific motivation for the development is to increase the integrity of the underlying fluid dynamics so that the physics terms can be more effectively isolated, examined, and improved. The expected benefits of the new model are discussed and results from the initial integrations will be presented.

Rood, Richard B.↗
NASA NTRS1999

Physically-Derived Dynamical Cores in Atmospheric General Circulation Models

The algorithm chosen to represent the advection in atmospheric models is often used as the primary attribute to classify the model. Meteorological models are generally classified as spectral or grid point, with the term grid point implying discretization using finite differences. These traditional approaches have a number of shortcomings that render them non-physical. That is, they provide approximate solutions to the conservation equations that do not obey the fundamental laws of physics. The most commonly discussed shortcomings are overshoots and undershoots which manifest themselves most overtly in the constituent continuity equation. For this reason many climate models have special algorithms to model water vapor advection. This talk focuses on the development of an atmospheric general circulation model which uses a consistent physically-based advection algorithm in all aspects of the model formulation. The shallow-water model of Lin and Rood (QJRMS, 1997) is generalized to three dimensions and combined with the physics parameterizations of NCAR's Community Climate Model. The scientific motivation for the development is to increase the integrity of the underlying fluid dynamics so that the physics terms can be more effectively isolated, examined, and improved. The expected benefits of the new model are discussed and results from the initial integrations will be presented.

Rood, Richard B.↗
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 NTRS2001

A Comparison of Cumulus Parameterizations in Idealized Sea-Breeze Simulations

Four cumulus parameterizations in the Penn State-NCAR model MM5 are compared in idealized sea-breeze simulations, with the aim of discovering why they work as they do. The most realistic results appear to be those using the Kain-Fritsch scheme. Rainfall is significantly delayed with the Betts-Miller-Janjic scheme, due to the method of computing the reference sounding. This method can be corrected, but downdrafts should be added in a physically realistic manner. Even without downdrafts, a corrected version of the BMJ scheme produces nearly the same timing and location of deep convection as the KF scheme, despite the very different physics. In order to simulate the correct timing of the rainfall, a minimum amount of mass is required in the layer that is the source of a parameterized updraft. The Grell parameterization, in the present simulation, always derives the updraft from the top of the mixed layer, where vertical advection predominates over horizontal advection in increasing the moist static energy. This makes the application of the quasi-equilibrium closure more correct than it would be if the updrafts were always derived from the most unstable layer, but it evades the question of whether or not horizontal advection generates instability. Using different physics, the parameterizations produce significantly different cloud-top heights.

Cohen, Charles↗
NASA NTRS2003

Use of GTE, Mozaic, Sonex, and UARS-MLS Data in Understanding Tropospheric Processes Critical to Model Development

During the first year we focused on the analysis of data collected on over 7600 commercial aircraft flights (the MOZAIC program). The aim was to further our understanding of the fundamental dynamical processes that drive mesoscale phenomena in the upper troposphere and lower stratosphere, and their effects on the advection of passive scalars. Through these studies we made the following findings. 2001]: We derived the Kolmogorov equation for the third-order velocity structure function on an f-plane. We showed how the sign of the function yields the direction of the energy cascade. The remarkable linearity of the measured off diagonal third-order structure function was studied. We suggested that the Coriolis term, which appears explicitly in this equation, may be crucial in understanding the observed kinetic energy spectra at scales larger than 100 km, instead of the nonlinear advection term as previously assumed. Also, we showed that decreases with latitude for mid- to high-latitudes, and tabulated the values. Ozone concentration structure functions were calculated by restricting the data points to approximately isentropic pairings. In this way we were able to make comparisons with scalar advection theories. We found that, at even the largest scales, there was no evidence for the simple scaling predicted for smooth advection/diffusion.

Newell, Reginald E.↗
NASA NTRS2004

Influence of Sub-grid-Scale Isentropic Transports on McRAS Evaluations using ARM-CART SCM Datasets

In GCM-physics evaluations with the currently available ARM-CART SCM datasets, McRAS produced very similar character of near surface errors of simulated temperature and humidity containing typically warm and moist biases near the surface and cold and dry biases aloft. We argued it must have a common cause presumably rooted in the model physics. Lack of vertical adjustment of horizontal transport was thought to be a plausible source. Clearly, debarring such a freedom would force the incoming air to diffuse into the grid-cell which would naturally bias the surface air to become warm and moist while the upper air becomes cold and dry, a characteristic feature of McRAS biases. Since, the errors were significantly larger in the two winter cases that contain potentially more intense episodes of cold and warm advective transports, it further reaffirmed our argument and provided additional motivation to introduce the corrections. When the horizontal advective transports were suitably modified to allow rising and/or sinking following isentropic pathways of subgrid scale motions, the outcome was to cool and dry (or warm and moisten) the lower (or upper) levels. Ever, crude approximations invoking such a correction reduced the temperature and humidity biases considerably. The tests were performed on all the available ARM-CART SCM cases with consistent outcome. With the isentropic corrections implemented through two different numerical approximations, virtually similar benefits were derived further confirming the robustness of our inferences. These results suggest the need for insentropic advective transport adjustment in a GCM due to subgrid scale motions.

Sud, Y. C.↗
NASA NTRS1999

On Multi-Dimensional Unstructured Mesh Adaption

Anisotropic unstructured mesh adaption is developed for a truly multi-dimensional upwind fluctuation splitting scheme, as applied to scalar advection-diffusion. The adaption is performed locally using edge swapping, point insertion/deletion, and nodal displacements. Comparisons are made versus the current state of the art for aggressive anisotropic unstructured adaption, which is based on a posteriori error estimates. Demonstration of both schemes to model problems, with features representative of compressible gas dynamics, show the present method to be superior to the a posteriori adaption for linear advection. The performance of the two methods is more similar when applied to nonlinear advection, with a difference in the treatment of shocks. The a posteriori adaption can excessively cluster points to a shock, while the present multi-dimensional scheme tends to merely align with a shock, using fewer nodes. As a consequence of this alignment tendency, an implementation of eigenvalue limiting for the suppression of expansion shocks is developed for the multi-dimensional distribution scheme. The differences in the treatment of shocks by the adaption schemes, along with the inherently low levels of artificial dissipation in the fluctuation splitting solver, suggest the present method is a strong candidate for applications to compressible gas dynamics.

Wood, William A.↗