Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Tracer equation”

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 55 records · Page 3

The equilibrium dayside magnetosphere

A method is presented of computing the dayside global earth magnetic field which is in equilibrium with the plasma pressure, based on satellite observations at a local region of the magnetosphere. The method, which utilizes a perturbation around a dipole magnetic field, involves computation of the global plasma pressure profile based on the equatorial (anisotropic) pressure data, derivation of the current profile which satisfies the equilibrium condition, and computation of the magnetic field using the current profile and the boundary current produced by the solar wind. The method is applied for the Active Magnetospheric Particle Tracer Explorers data, and the result of the computation is found to compare reasonably well with the observed magnetic field profile near the geomagnetic equator.

Zavriyev, Anton↗

Evidence for eddy formation in the eastern Arabian Sea during the northeast monsoon

The seasonal formation of a large (500-800 km diameter) anticyclonic eddy in the upper 300-400 m of the eastern Arabian Sea during the northest monsoon period (December-April) is indicated fom hydrographic and satellite altimetry sea level observations, as well as from numerical model experiments. The center of the eddy circulation is approximately 10 deg N, 70 deg E, just to the west of the north-south Laccadive Island chain. In this paper the eddy is called the Laccadive High (LH). In some ways it is like a mirrorlike counterpart to the Great Whirl that develops during the southwest monsoon of the Somali coast (western Arabian Sea). The LH occurs at the same latitude but on the opposite side of the basin during the reversed monsoon. It is different from the Great Whirl, however, in its formation process, its intensity, and its decay. The hydrographic data obtained from surveys all during a single season give sufficiently close station spacing to allow reasonable contouring of the geopotential surfaces and of the properties within and around the LH region with minimum time aliasing. The Geostat altimeter record extends over 4 years, during which the seasonal variability of the LH indicates a dynamic relief of approximately 15-20 cm, which is in good agreement with the hydrographics observations. The altimetry time series also suggests a westward translation of the LH by January with a subsequent dissipation in midbasin. The model used is a wind-forced three-layer primitive equation model which depicts a LH agreement with the timing, position, and amplitude of both the hydrographic and altimetric measurements. The numerical simulation includes a passive tracer located in the Western Bay of Bengal; the western advection of the tracer around the south coasts of Sri Lanka and India in December and January is consistent with the appearance of low-salinity water observed to extend into the Arabian Sea during this period. The modeling studies suggest that both local and remote forcing are important in formation of the LH.

Bruce, John G.↗

The NASA Ames Research Center one- and two-dimensional stratospheric models. Part 2: The two-dimensional model

The two-dimensional model of stratospheric constituents is presented in detail. The derivation of pertinent transport parameters and the numerical solution of the species continuity equations, including a technique for treating the stiff differential equations that represent the chemical kinetic terms, and appropriate methods for simulating the diurnal variations of the solar zenith angle and species concentrations are discussed. Predicted distributions of tracer constituents (ozone, carbon 14, nitric acid) are compared with observed distributions.

Whitten, R. C.↗

A study of the accuracy of neutrally buoyant bubbles used as flow tracers in air

Research has been performed to determine the accuracy of neutrally buoyant and near neutrally buoyant bubbles used as flow tracers in air. Theoretical, computational, and experimental results are presented to evaluate the dynamics of bubble trajectories and factors affecting their ability to trace flow-field streamlines. The equation of motion for a single bubble was obtained and evaluated using a computational scheme to determine the factors which affect a bubble's trajectory. A two-dimensional experiment was also conducted to experimentally determine bubble trajectories in the stagnation region of NACA 0012 airfoil at 0 deg angle of attack using a commercially available helium bubble generation system. Physical properties of the experimental bubble trajectories were estimated using the computational scheme. These properties included the density ratio and diameter of the individual bubbles. the helium bubble system was then used to visualize and document the flow field about a 30 deg swept semispan wing with simulated glaze ice. Results were compared to Navier-Stokes calculations and surface oil flow visualization. The theoretical and computational analysis have shown that neutrally buoyant bubbles will trace even the most complex flow patterns. Experimental analysis revealed that the use of bubbles to trace flow patterns should be limited to qualitative measurements unless care is taken to ensure neutral buoyancy. This is due to the difficulty in the production of neutrally buoyant bubbles.

Kerho, Michael F.↗

The Origin of Antarctic Precipitation: A Modeling Approach

Isotope concentrations in polar ice cores have long been used to estimate paleotemperatures. Underlying the use of this "isotope paleothermometer" is the assumption that the relationship between surface temperature and isotope concentration over time at a single geographical point is the same as that observed over space during the present-day climate. The validity of this assumption may in fact be compromised by several factors related to climate change. The specific factor studied in this paper involves the evaporative sources for polar precipitation. Climatic changes in the relative strengths of these sources would imply a need for a recalibration of the paleothermometer. To quantify such changes, we performed two GCM simulations, one of present-day climate and the other of the climate during the Last Glacial Maximum (LGM), roughly 18000 years ago. Evaporative sources of Antarctic precipitation were established using special tracer diagnostics. Results suggest that polar precipitation during the LGM does indeed consist of (relatively) more water from tropical oceans, a direct reflection of the LGM's increased equator-to-pole temperature gradient and its increased sea ice extent, which reduces high latitude evaporation. This result implies that an uncalibrated ice core paleothermometer would produce LGM temperatures that are biased slightly low. Because LGM boundary conditions are still under debate, we performed a third GCM simulation using a modified set of LGM boundary conditions. Using this simulation gives some qualitatively similar results, though the tropical contribution is not quite as high. Uncertainties in the LGM boundary conditions does hamper success in calibrating the paleothermometer.

Delaygue, Gilles↗

Finite element procedures for time-dependent convection-diffusion-reaction systems

New finite element procedures based on the streamline-upwind/Petrov-Galerkin formulations are developed for time-dependent convection-diffusion-reaction equations. These procedures minimize spurious oscillations for convection-dominated and reaction-dominated problems. The results obtained for representative numerical examples are accurate with minimal oscillations. As a special application problem, the single-well chemical tracer test (a procedure for measuring oil remaining in a depleted field) is simulated numerically. The results show the importance of temperature effects on the interpreted value of residual oil saturation from such tests.

Tezduyar, T. E.↗

Missing matter in the vicinity of the sun

The Poisson and Vlasov equations are solved numerically for realistic Galaxy models which include multiple disk components, a Population II spheroid, and an unseen massive halo. The total amount of matter in the vicinity of the sun is determined by comparing the observed distributions of tracer stars, samples of F dwarfs, and K giants with the predictions of the Galaxy models. Results are obtained for a number of different assumed distributions of the unseen disk mass. For all the observed samples, typical models imply that about half of the mass in the solar vicinity must be in the form of unobserved matter. The volume density of unobserved material near the sun is about 0.1 solar mass/cu pc; the corresponding column density is about 30 solar mass/sq pc. This so far unseen material must be in a disk with an exponential scale height of less than 0.7 kpc.

Bahcall, John N.↗

The NASA Ames Research Center one- and two-dimensional stratospheric models. Part 1: The one-dimensional model

A one-dimensional model of stratospheric trace constituents is described in detail. Specifically, the numerical solution of the species continuity equations, including a technique for treating the stiff differential equations representing the chemical kinetic terms, and an appropriate method for simulating the diurnal variations of the species concentrations, are discussed. A specialized treatment of atmospheric photodissociation rates is outlined in the text. The choice of a vertical eddy diffusivity profile and its success in predicting the vertical tracer distributions (carbon 14, methane, and nitrous oxide) are also discussed.

Turco, R. P.↗

EOS Interdisciplinary Investigation: Observational and Modeling Studies of Radiative, Chemical, and Dynamical Interactions in the Earth's Atmosphere

A 3-dimensional model of atmospheric dynamics and photochemistry has been developed from the primitive equations in isentropic coordinates. The model extends from the upper troposphere through the middle atmosphere and is driven by observed dynamical and chemical structure in the troposphere. The model's formulation is entirely spectral: Horizontal structure is represented in terms of vector Hough functions, which explicitly partition the motion into its rotational and divergent components. This formulation enables the model's computational performance to be increased dramatically by filtering high-frequency gravity waves, which do not affect PV conservation. Vertical structure is represented in terms of eigenfunctions that follow from the primitive equations in isentropic coordinates. The model's fully spectral formulation enables scale-selective dissipation, necessary for numerical stability, to be applied at 6th order-in all 3 coordinate directions. This feature leaves all but the shortest vertical scales undamped and, consequently, allows potential vorticity and chemical tracers to be conserved quite accurately, distinctly better than in the model's layered counterpart. These dynamical features are complemented by a basic but fairly complete treatment of gas phase photochemistry, which accounts for some 48 chemical species, diurnally-varying SW absorption by O2 and O3, and LW cooling calculated via a detailed band calculation of radiative transfer. Driven by observed tropospheric structure, the 3D model reproduces observed structure in the middle atmosphere, inclusive of transport by the Brewer-Dobson circulation and accompanying chemical variations. Calculated residual motion is consistent with diabatic cooling rates, poleward transport, and descent in the polar night derived from satellite measurements of chemical tracers made from LIMS, SAGE, and UARS. The model also reproduces the observed structure of chemical species like O3, HNO3, and ClO. At the same time, it provides a respectable 3D simulation of the wintertime increase of midlatitude total ozone, as well as its geographical structure observed by Nimbus-7 TOMS. The 3D model has been used to explore the origin of diabatic upwelling in the tropical lower stratosphere, which has been related to planetary wave absorption in the winter- time stratosphere, as well as to influences from tropospheric convection at the tropical tropopause.

Salby, Murry L.↗

Motions and Initial Conditions in Star-Forming Dense Cores

The main focus was the study of star-forming regions through high spectral- and spatial resolution observations of mm-wavelength lines, and through models of the observations. The main accomplishments were a) demonstration that more than 15 starless cores show substantial evidence of extended inward motion at about half the sound speed; b) observations of infall asymmetry in several cores, in lines of N2H(+) and DCO(+), low- depletion tracers of the "inner core"; c) observation of "infall asymmetry" of spectral lines over approx. 0.5 pc in the NGC1333 cluster-forming region; d) observations indicating that cores are nearly at rest with respect to their envelopes; and e) development of analytic, power-series solutions to the equations of motions for condensing 1-D systems (layers, cylinders and spheres).

Myers, Philip C.↗

Variability in the Speed of the Brewer-Dobson Circulation as Observed by Aura/MLS

We use Aura/MLS stratospheric water vapour (H2O) measurements as tracer for dynamics and infer interannual variations in the speed of the Brewer-Dobson circulation (BDC) from 2004 to 2011. We correlate one-year time series of H2O in the lower stratosphere at two subsequent pressure levels (68 hPa, approx.18.8 km and 56 hPa, approx 19.9 km at the Equator) and determine the time lag for best correlation. The same calculation is made on the horizontal on the 100 hPa (approx 16.6 km) level by correlating the H2O time series at the Equator with the ones at 40 N and 40 S. From these lag coefficients we derive the vertical and horizontal speeds of the BDC in the tropics and extra-tropics, respectively. We observe a clear interannual variability of the vertical and horizontal branch. The variability reflects signatures of the Quasi Biennial Oscillation (QBO). Our measurements confirm the QBO meridional circulation anomalies and show that the speed variations in the two branches of the BDC are out of phase and fairly well anti-correlated. Maximum ascent rates are found during the QBO easterly phase. We also find that transport of H2O towards the Northern Hemisphere (NH) is on the average two times faster than to the Southern Hemisphere (SH) with a mean speed of 1.15m/s at 100 hPa. Furthermore, the speed towards the NH shows much more interannual variability with an amplitude of about 21% whilst the speed towards the SH varies by only 10 %. An amplitude of 21% is also observed in the variability of the ascent rate at the Equator which is on the average 0.2mm/s.

Brewer-Dobson circulation↗

Remarks on the overall distribution of hydrogen in the galactic disk

Several current problems concerning the overall distribution of hydrogen in the galaxy are discussed in general terms. These problems include the degree of saturation characterizing low-latitude emission observations of HI, and the optical-depth corrections to the derived column and volume densities; the amount of fine-scale velocity and spatial structure diluted by the instrumental limitations of the presently available surveys; and the general problem of detailed mapping of the HI in the galaxy. Comparison is made between the distribution of HI and that of CO and several other galactic tracers. The galactic disk as defined by atomic hydrogen has a diameter fully twice as large as that defined by the ionized and molecular states of hydrogen, as well as by other molecules, supernova remnants, pulsars, gamma radiation, synchroton radiation, and the youngest stars. It is also less confined to the galactic equator than most of the other constituents. The degree of small-scale structure apparent in the molecular observations is much greater than in the HI observations. Parameters describing the small-scale structure were determined using Monte Carlo techniques to simulate the observations.

Burton, W. B.↗

Performance limitations in parallel processor simulations

A jet-engine model is partitioned and simulated on a parallel processor system consisting of five 8086/8087 floating-point computers. The simulation uses Heun's integration method. A near-optimal parallel simulation (in the sense of minimum execution time) achieves speedup of only 2.13 and efficiency of 42.6 percent, in effect wasting 57.4 percent of the available processing power. A detailed analysis identifies and graphically demonstrates why the system fails to achieve ideal performance (viz., speedup of 5 and efficiency of 100 percent). Inherent characteristics of the problem equations and solution algorithm account for the loss of nearly half of the available processing power. Overheads associated with interprocessor communication and processor synchronization account for only a small fraction of the lost processing power. The effects of these and other factors which limit parallel processor performance are illustrated through real-time timing-analyzer tracers describing the run/idle status of the parallel processors during the simulation.

O'Grady, E. Pearse↗

Low Ozone in the Marine Boundary Layer of the Tropical Pacific Ocean

Aircraft measurements of ozone, its key precursors, and a variety of chemical tracers were made in the troposphere of the western and central Pacific in October 1991. These data are presented and analyzed to examine the occurrence of low ozone concentrations in the remote marine boundary layer of the tropical and equatorial Pacific Ocean. The data from these flights out of Guam, covering an area extending from the equator to 20 N and from south of the Philippines to Hawaii, show average O3 concentrations as low as 8-9 ppb (ppb=10(exp-9)v/v) at altitudes of 0.3-0.5 km in the boundary layer. Individual measurements as low as 2-5 ppb were recorded. Low O3 concentrations do not always persist in space and time. High O3, generally associated with the transport of upper tropospheric air, was also encountered in the boundary layer. In practically all cases, O3 increased to values as large as 25-30 ppb within 2 km above the boundary layer top. Steady state model computations are used to suggest that these low O3 concentrations are a result of net photochemical O3 destruction in a low NO environment, sea-surface deposition, and extremely low net entrainment rates (1-2 mm per second) from the free troposphere. Day/night measurements of ethane, propane, gaseous and aerosol Cl suggest that daytime (morning) Cl atom concentrations in the vicinity of 10(exp 5) molecules per cubic centimeter may be present in the marine boundary layer. This Cl atom abundance can be rationalized only if sea salt aerosols can release free chlorine (Cl2) to the gas phase in the presence of sun light (and possibly O3). These Cl atom concentrations, however, are still insufficient and Cl (or Br) chemistry is not likely to be an important cause of the observed low O3.

Singh, Hanwant B.↗

High Resolution Numerical Simulations and Diagnostic Studies Atmospheric Transport During TRACE-A

The enclosed publications constitute our final report. These are publications completed in referred journals. The completed work includes: Meteorological conditions associated with vertical distributions of aerosols off the west coast of Africa. TRACE-A trajectory intercomparison 2. Isentropic and kinematic methods. Chemical characteristics of tropospheric air over the tropical South Atlantic Ocean: Relationship to trajectory history. Passive Tracer Transport Relevant to the TRACE-A Experiment. The meteorological environment of the tropospheric ozone maximum over the tropical south Atlantic Ocean. : Influence of a middle-latitude cyclone on tropospheric ozone distributions during a period of TRACE-A. All of this work provided the meteorological background for the Transport of Atmospheric Chemistry near the Equator (TRACE-A) project sponsored by Dr,, Joe McNeal of NASA. Our principal findings are that bio-mass burning sources from Africa and South America did contribute to the accumulation of tropospheric ozone over the tropical Atlantic Ocean. Other findings relate to circulations, advections of Ozone and other plausible sources for the TOMS based ozone maximum over the Atlantic Ocean.

Krishnamurti, T. N.↗

Mixed Layer Sub-Mesoscale Parameterization - Part 1: Derivation and Assessment

Several studies have shown that sub-mesoscales (SM 1km horizontal scale) play an important role in mixed layer dynamics. In particular, high resolution simulations have shown that in the case of strong down-front wind, the re-stratification induced by the SM is of the same order of the de-stratification induced by small scale turbulence, as well as of that induced by the Ekman velocity. These studies have further concluded that it has become necessary to include SM in ocean global circulation models (OGCMs), especially those used in climate studies. The goal of our work is to derive and assess an analytic parameterization of the vertical tracer flux under baroclinic instabilities and wind of arbitrary directions and strength. To achieve this goal, we have divided the problem into two parts: first, in this work we derive and assess a parameterization of the SM vertical flux of an arbitrary tracer for ocean codes that resolve mesoscales, M, but not sub-mesoscales, SM. In Part 2, presented elsewhere, we have used the results of this work to derive a parameterization of SM fluxes for ocean codes that do not resolve either M or SM. To carry out the first part of our work, we solve the SM dynamic equations including the non-linear terms for which we employ a closure developed and assessed in previous work. We present a detailed analysis for down-front and up-front winds with the following results: (a) down-front wind (blowing in the direction of the surface geostrophic velocity) is the most favorable condition for generating vigorous SM eddies; the de-stratifying effect of the mean flow and re-stratifying effect of SM almost cancel each other out,

Canuto, V. M.↗

Predictions and measurements of isothermal flowfields in axisymmetric combustor geometries

Numerical predictions, flow visualization experiments and time-mean velocity measurements were obtained for six basic nonreacting flowfields (with inlet swirl vane angles of 0 (swirler removed), 45 and 70 degrees and sidewall expansion angles of 90 and 45 degrees) in an idealized axisymmetric combustor geometry. A flowfield prediction computer program was developed which solves appropriate finite difference equations including a conventional two equation k-epsilon eddy viscosity turbulence model. The wall functions employed were derived from previous swirling flow measurements, and the stairstep approximation was employed to represent the sloping wall at the inlet to the test chamber. Recirculation region boundaries have been sketched from the entire flow visualization photograph collection. Tufts, smoke, and neutrally buoyant helium filled soap bubbles were employed as flow tracers. A five hole pitot probe was utilized to measure the axial, radial, and swirl time mean velocity components.

Rhodes, D. L.↗

Remarks on the overall distribution of hydrogen in the galactic disk

The degree of saturation characterizing low altitude emission observations of H(I) and the optical depth corrections to the derived column and volume densities are discussed. The amount of fine-scale velocity and spatial structure diluted by the instrumental limitations of the presently available surveys are described. The general problem of detailed mapping of H(I) in the galaxy is explored. Comparison is made between the distribution of H(I) and that of CO and several other galactic tracers. Atomic hydrogen is unique in its distribution, instead of being typical of many Population I constituents. As defined by atomic hydrogen, the galactic disk has a diameter fully twice as large as that defined by the ionized and molecular states of hydrogen, as well as by other molecules, supernova remnants, pulsars, gamma-radiation, synchroton radiation, and the youngest stars. It is also less confined to the galactic equator than most of the other constituents. The degree of small scale structure apparent in the molecular observations is much greater than that in the H(I) observations.

Burton, W. B.↗