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 703 records · Page 39

A diagnostic study of baroclinic disturbances in polar air streams

Quasi-Lagrangian budgets of mass, vorticity and heat are calculated following disturbances that form within polar air streams. Observed cases are extracted from the European Centre for Medium-range Weather Forecasts analyses during the First GARP Global Experiment. Model-generated cases are extracted from the simulations of extratropical cyclogenesis by Sandgathe. These polar lows grow primarily through basic baroclinic instability processes and exhibit many features of larger maritime extratropical cyclones. Polar lows that originate on the poleward (or Cyclonic - Type C) side of the jet and have considerable midtropospheric positive vorticity advection at formation time are contrasted with lows that form on the equatorward (or Anticyclonic - Type A) side of a nearly straight upper-level jet. The midtropospheric positive vorticity advection must be present to enhance the vertical circulation when the large surface fluxes that are associated with strong outbreaks act to damp the thermal wave amplification. Although latent heat release is an important factor in both types, it is an essential energy source for the Type A low developments on the equatorward side. Although the vorticity balance is initially different for the two types of polar lows, the vorticity budgets during later stages are similar. The heat budget and the thickness tendency equation demonstrated that the self-development process that is present in larger maritime cyclones is also important for polar low intensification. The absence of favorable coupling to a jet stream is the missing factor in a model-generated Type A polar low that failed to develop. Consequently, the mid- and upper-tropospheric wind fields determine which polar lows will intensify to significant amplitudes.

Sinclair, Mark R.↗

Equatorial long-wave characteristics determined from satellite sea surface temperature and in situ data

SST maps and imagery derived from the NOAA 6 satellite AVHRR for June and July 1981 in the eastern tropical Pacific portray the wavelike structure of the cool water along the equator from 93 deg W to 125 deg W. Cusped waves of approximately 1000-km zonal wavelength and 25-day period propagated westward with a phase speed of 40 km/day. The observed meridional extent between the crest and trough of the waves is about 300 km. Details in the imagery show cooler water at the cusps advected north and then east with the north equatorial countercurrent (NECC), consistent with the suggestion of a series of anticyclonic eddies occupying the shear zone between the NECC and the westward flowing south equatorial current. Absolute SST estimates from the AVHRR data agree to within 0.6 C with shipboard data taken along 110 deg W between 5 deg N and 5 deg S. The wavelike structures in the SST maps are also in agreement at the surface with a vertical expendable bathythermograph temperature section made along the equator between 93 deg W and 125 deg W, which shows the phase of the waves tilting westward with increasing depth over the upper 75 m. Such a phase shift, if it extended 100-200 km meridionally in either direction from the equator, would be associated with an equatorward flux of heat. Similar phase shifts appear in temperature time series at depths of 20 and 50 m, from a mooring at 0-deg 33-min N, 110-deg 30-min W. Near-surface currents measured at this and a second mooring on the equator at 109-deg 40-min W indicate a regular pattern of northward advection when wave cusps pass them, followed by southwest flow during the passage of wave troughs, again consistent with an equatorward flux of heat, as well as with earlier theoretical and drift buoy findings.

Pullen, Patricia E.↗

Soil moisture estimation using GOES-VISSR infrared data - A case study with a simple statistical method

Five days of clear sky observations of Kansas and Nebraska are used to examine the statistical relationship between soil moisture and infrared surface temperature observations taken from a geosynchronous satellite. Linear regression is used to relate soil moisture to surface temperature and other variables that represent wind speed, vegetation cover, and low-level temperature advection. Results show good agreement between estimated and observed soil moisture features on each of the 5 days. The average coefficient of determination for five pseudoindependent tests in which the test day is held out of the regression is 0.71. It is shown that a depletion coefficient of 0.92, when used to compute antecedent precipitation index (API), produces the best correlation between API and soil moisture as inferred from GOES thermal infrared data. By averaging daily predicted values over the 5-day rain-free case study period, 92 percent of the variance of the morning surface temperature change is explained by a simple multiple linear regression with all independent variables, or, alternatively, 85 percent of the observed variance in API is explained. It is concluded that this approach can distinguish at least four classes of soil wetness, but the necessity for measurement of surface advection may limit its usefulness in remote areas.

Wetzel, Peter J.↗

Measurements of surface energy budgets in the rocky mountains of Colorado

The surface energy balance of a mountain valley and a mountain peak were continuously monitored for periods of 126 and 34 days, respectively. The data yielded temporally resolved, high-quality measurements of the radiative, soil, and total atmospheric turbulent heat fluxes. Comparative data for the two sites reveal the effects of terrain setting and surface cover in modulating the heat balance. A procedure is suggested for accommodating local advective affects due to variable terrain and surface cover to allow partitioning of the total turbulent flux into values for sensible and latent heat. This procedure relies on a stability-dependent coefficient of turbulent transfer which is locally calibrated to systematically minimize advection-related errors. Scaling functions for the transfer coeffient are based on the effective blackbody temperature of the surface. Consequently, the procedure offers an opportunity for developing reliable bulk representations of heat transfer at the surface for use with remote sensing data and with models of the soil and surface cover in complex terrain.

Sheaffer, John D.↗

Cellular automaton formulation of passive scalar dynamics

Cellular automata modeling of the advection of a passive scalar in a two-dimensional flow is examined in the context of discrete lattice kinetic theory. It is shown that if the passive scalar is represented by tagging or 'coloring' automation particles a passive advection-diffusion equation emerges without use of perturbation expansions. For the specific case of the hydrodynamic lattice gas model of Frisch et al. (1986), the diffusion coefficient is calculated by perturbation.

Chen, Hudong↗

Simulated sudden stratospheric warming - Synoptic evolution

An analysis is presented of a sudden stratospheric warming event which occurred spontaneously during a general circulation model simulation of the global atmospheric circulation. Two separate warming pulses exhibit the same dynamical evolution with a 'cycle' of about two weeks. Two distinct phases of the warming cycle are apparent: (1) the generation of an intense localized warm cell in conjunction with significant adiabatic heating associated with cross-isobar flow which has been induced by vertically propagating long wave disturbances; and (2) the northward transport of that warm cell via advection by the essentially geostrophic windfield corresponding to an intense, offset polar cyclone, in conjunction with a strong Aleutian anticyclone. During the first warming pulse in January, a moderate Aleutian anticyclone was in place prior to the warming cycle and was intensified by interaction with an eastward traveling anticyclone induced by the differential advection of the warm cell. The second warming pulse occurred in early February with a strong Aleutian anticyclone already established. In contrast to the January event, the warming in February culminated with reversal of the zonal westerlies to easterlies over a significant depth of the stratosphere.

Blackshear, W. T.↗

Recent progress in the development and understanding of SUPG methods with special reference to the compressible Euler and Navier-Stokes equations

The current status of streamline-upwind/Petrov-Galerkin (SUPG) methods for the analysis of flow problems is surveyed in an analytical review. Problem areas addressed include classical Galerkin, upwind, artificial-diffusion, SUPG, discontinuous Galerkin, space-time FEM, and discontinuity-capturing approaches to the scalar advection-diffusion equation; incompressible flows; advective-diffusive systems; and the compressible Euler and Navier-Stokes equations. Graphs and diagrams are provided, and the good stability properties of state-of-the-art SUPG methods are pointed out.

Hughes, Thomas J. R.↗

Diagnostic study of explosive cyclogenesis during FGGE

Two rapidly developing extratropical maritime cyclones (one, developed during January 13-15, 1979, along an intense frontal zone south of Japan, the other, of January 26-27, 1979, in a polar mass over the North Atlantic) were investigated using the FGGE data and ECMWF level IIIb analyses to describe the structure and dynamics of these events. Although the cyclones evolved from a strong low-level baroclinic zone without initial large midtropospheric vorticity advection, thus resembling the Petterssen type A development, rapid deepening occurred in both cases when an approaching upper tropospheric jet with appreciable shear vorticity advection became favorably superposed over the surface low. During the development period, stability decreased in the low troposphere, aiding in the rapid development of an intense mass-circulation and low tropospheric vorticity production by the divergence term. The results suggest that upper-level forcing plays a greater role in the initiation of explosive oceanic development than is suggested by the Petterssen and Smebye (1971) description.

Wash, Carlyle H.↗

Third-order multi-dimensional Euler/Navier-Stokes solver

A nonoscillatory multidimensional advection code is presented which is applicable to the study of very high Reynolds number steady flow, and results are obtained for the oblique advection of a step profile. The scheme is based on an exponential upwinding or linear extrapolation refinement of a previous QUICK scheme, allowing sharp monotonic resolution of the step while retaining third-order accuracy. The resulting simple flux-limiting criteria have been applied to a higher order method and a formally third-order supercompressive algorithm. Results are presented for first-order, second-order, and third-order upwinding.

Leonard, B. P.↗

Rossby wave activity in a two-dimensional model - Closure for wave driving and meridional eddy diffusivity

A parameterization of the effects of Rossby waves in the middle atmosphere is proposed for use in two-dimensional models. By adding an equation for conservation of Rossby wave activity, closure is obtained for the meridional eddy fluxes and body force due to Rossby waves. Rossby wave activity is produced in a climatological fashion at the tropopause, is advected by a group velocity which is determined solely by model zonal winds, and is absorbed where it converges. Absorption of Rossby wave activity causes both an easterly torque and an irreversible mixing of potential vorticity, represented by the meridional eddy diffusivity, K(yy). The distribution of Rossby wave driving determines the distribution of K(yy), which is applied to all of the chemical constituents. This provides a self-consistent coupling of the wave activity with the winds, tracer distributions and the radiative field. Typical winter stratospheric values for K(yy) of 2 million sq m/sec are obtained. Poleward tracer advection is enhanced and meridional tracer gradients are reduced where Rossby wave activity is absorbed in the model.

Hitchman, Matthew H.↗

Cut-and-connect of two antiparallel vortex tubes

Motivated by an early conjecture that vortex cut-and-connect plays a key role in mixing and production of turbulence, helicity and aerodynamic noise, the cross-linking of two antiparallel viscous vortex tubes via direct numerical simulation is studied. The Navier-Stokes equations are solved by a dealiased pseudo-spectral method with 64 cubed grid points in a periodic domain for initial Reynolds numbers Re up to 1000. The vortex tubes are given an initial sinusoidal perturbation to induce a collision and keep the two tubes pressed against each other as annihilation continues. Cross-sectional and wire plots of various properties depict three stages of evolution: (1) Inviscid induction causing vortex cores to first approach and form a contact zone with a dipole cross-section, and then to flatten and stretch; (2) Vorticity annihilation in the contact zone accompanied by bridging between the two vortices at both ends of the contact zone due to a collection of cross-linked vortex lines, now orthogonal to the initial vortex tubes. The direction of dipole advection in the contact zone reverses; and (3) Threading of the remnants of the original vortices in between the bridges as they pull apart. The crucial stage 2 is shown to be a simple consequence of vorticity annihilation in the contact zone, link-up of the un-annihilated parts of vortex lines, and stretching and advection by the vortex tube swirl of the cross-linked lines, which accumulate at stagnation points in front of the annihilating vortex dipole. It is claimed that bridging is the essence of any vorticity cross-linking and that annihilation is sustained by stretching of the dipole by the bridges. Vortex reconnection details are found to be insensitive to asymmetry. Modeling of the reconnection process is briefly examined. The 3D spatial details of scalar transport (at unity Schmidt number), enstrophy production, dissipation and helicity are also examined.

Melander, Mogens V.↗

Analysis of sea ice dynamics

The ongoing work has established the basis for using multiyear sea ice concentrations from SMMR passive microwave for studies of largescale advection and convergence/divergence of the Arctic sea ice pack. Comparisons were made with numerical model simulations and buoy data showing qualitative agreement on daily to interannual time scales. Analysis of the 7-year SMMR data set shows significant interannual variations in the total area of multiyear ice. The scientific objective is to investigate the dynamics, mass balance, and interannual variability of the Arctic sea ice pack. The research emphasizes the direct application of sea ice parameters derived from passive microwave data (SMMR and SSMI) and collaborative studies using a sea ice dynamics model. The possible causes of observed interannual variations in the multiyear ice area are being examined. The relative effects of variations in the large scale advection and convergence/divergence within the ice pack on a regional and seasonal basis are investigated. The effects of anomolous atmospheric forcings are being examined, including the long-lived effects of synoptic events and monthly variations in the mean geostrophic winds. Estimates to be made will include the amount of new ice production within the ice pack during winter and the amount of ice exported from the pack.

Zwally, J.↗

Numerical simulation of tropical cumulus clouds and their interaction with the subcloud layer

A two-dimensional numerical model, with an emphasis on turbulent processes in the boundary layer and in clouds, was developed for simulating an ensemble of cumulus clouds. The model was used to investigate the response of tropical cumulus clouds to imposed large-scale vertical advection by performing two simulations which were identical except that one was with upward large-scale vertical velocity (the 'disturbed' case) and the other was without any advection (the 'undisturbed' case). The results showed that deep cumulus clouds formed in the disturbed case, while only shallow clouds formed in the undisturbed case. The subcloud layer (SCL) in the disturbed case was strongly affected by cumulus circulations and rain evaporation, with both cumulus updrafts and downdrafts being important contributors to the SCL sensible and latent heat fluxes. Cumulus-scale circulation in the disturbed SCL determined where new clouds formed by creating convergence zones and moisture anomalies.

Krueger, Steven K.↗

Nonlinear Hadley circulation driven by asymmetric differential heating

The dynamical state of the stratosphere influenced by radiative heating, with no internal sources or sinks of angular momentum, is examined. It is shown that there exists a nonlinear Hadley regime driven by antisymmetric (or more generally, asymmetric) thermal equilibria typical of the middle atmosphere at the solstices. This regime consists of a single mean meridional cell, equatorial easterlies, and strong winter westerlies. Outside of the circulation region the flow is in thermal equilibrium. The effect of one-sided friction, acting as a drag on midlatitude westerlies only, is to expand the Hadley cell into the winter hemisphere and increase the magnitude of cross-equatorial flow. This result is possible even in the steady state when the advection of angular momentum in the tropics is made small by reducing the gradient of angular momentum in this region instead of the advecting velocity.

Dunkerton, Timothy J.↗

The sensitivity of total ozone and ozone perturbation scenarios in a two-dimensional model due to dynamical inputs

The sensitivity of the total ozone distribution in the Douglass et al. (1989) photochemical model to the dynamical inputs of ozone and the zone perturbation scenarios was investigated by computing the residual circulation from three sets of heating rates (used in computing the advection field). The model ozone results were compared to solar backscattered UV ozone data (SBUV). It was found that, depending on the dynamical inputs, the modeled global average total ozone varied up to 10 percent (with even more substantial seasonal and latitudinal variabilities), indicating that total ozone in the middle to high latitudes is very sensitive to the advection field below 100 mbar, the region where the heating rates are most uncertain. The modeled total ozone showed better agreement with solar SBUV climatology when computed horizontal eddy diffusion was used.

Jackman, Charles H.↗

On the mechanisms responsible for high-latitude thermospheric composition variations during the recovery phase of a geomagnetic storm

The causal mechanisms for the recovery of the perturbed high-latitude thermospheric composition to the unperturbed state in the period following a geomagnetic storm are investigated. Model runs of the NCAR thermosphere/ionosphere GCM (TIGCM) and thermosphere GCM are used to calculate the averaged mass mixing ratio variations and the forcing terms responsible for these variations during the recovery phase of a geomagnetic storm. High latitude compositional recovery is found to occur in the NCAR TIGCM on a time scale of about 12 hr to 1 day. This time scale is in agreement with previously observed time scales for typical poststorm F region electron density recoveries. Neither molecular diffusion nor large-scale horizontal advection is the dominant process in determining the compositional state during the recovery period. Thermospheric compositional recovery at high geomagnetic latitudes is driven primarily by vertical advection.

Burns, A. G.↗

A study of numerical methods for hyperbolic conservation laws with stiff source terms

In the present study of the behavior of typical numerical methods in the case of a model advection equation having a parameter-dependent source term, two approaches to the incorporation of the source terms are used: MacCormack-type predictor-corrector methods with flux limiters, and splitting methods in which the fluid dynamics and chemistry are handled in separate steps. The latter are found to perform slightly better. The model scalar equation is used to show that the incorrectness of the propagation speeds of discontinuities observed in the stiff case is due to the introduction of nonequilibrium values through numerical dissipation in the advection step.

Leveque, R. J.↗

Synoptic/planetary-scale interactions and blocking over the North Atlantic Ocean

Work was completed on the height tendency diagnoses of two extratropical cyclones that occurred upstream from the blocking event studied previously. One developed explosively over water 60 to 36 hours before the block first appeared, while the second developed explosively over the southeastern United States during the time of block formation. In both cases, both vorticity and temperature advection were consistently important forcing mechanisms. This is in contrast to the block itself, in which vorticity advection was easily the dominant forcing mechanism. Latent heat release was also significant, accounting for about 50 percent of the total height falls in the cyclone below 850 mb. Estimates of latent heat release were greatly enhanced by coupling parameterized estimates with values derived from GOES IR data using an algorithm developed by Marshall's F. R. Robertson. Among the difficulties encountered in this work was the identification of an appropriate lower boundary condition for the solution of the height tendency equation. The zero value currently used tends to yield underestimates of the lower troposphere height tendencies. To address this problem a new diagnostic technique was developed in cooperation with Dr. Peter Zwack of the University of Quebec at Montreal. Based on an equation Dr. Zwack had previously developed (the Zwack-Okossi development equation), researchers now have a relationship that is completely consistent with the height tendency equation and provides estimates of lower boundary geostrophic vorticity or height tendencies. Finally, comparison of the SAT (satellite data) and NOSAT (no satellite data) analyses is progressing well. The present focus is on both the new diagnostic technique and the SAT/NOSAT comparisons. The former is being tested on the southeastern United States cyclone case previously mentioned and compared with the height tendency diagnoses already completed. The latter are being examined for the blocking case described in the publications cited in this summary. In addition to obtaining statistics that will allow general comparison of the two analyses, it will be possible to determine whether conclusions about the dynamics of the block development are influenced by the analysis set used.

Smith, Phillip J.↗