SEARCH · Engineering Papers
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.
During the session on modeling of the planetary boundary layer (PBL) processes the following topics were addressed: current PBL parameterizations used in mesoscale models; critical issues in improving PBL parameterizations for mesoscale models; suggestions for future modeling efforts; and suggestions for observing system strategies.
This project involves collaborative research with Stephen W. Pacala and Simon A. Levin of Princeton University to calibrate, test, and analyze models of heterogeneous forested landscapes containing a diverse array of habitats. The project is an extension of previous, NASA-supported research to develop a spatially-explicit model of forest dynamics at the scale of an individual forest stand (hectares to square kilometer spatial scales). That model (SORTIE) has been thoroughly parameterized from field studies in the modal upland environment of western Connecticut. Under our current funding, we are scaling-up the model and parameterizing it for the broad range of upland environments in the region. Our most basic goal is to understand the linkages between stand-level dynamics (as revealed in our previous research) and landscape-level dynamics of forest composition and structure.
A parameterized reaction mechanism has been created to examine the interactions between isoprene and other tropospheric gas-phase chemicals. Tests of the parameterization have shown that its results match those of a more complex reaction set to a high degree of accuracy. Comparisons between test runs have shown that the presence of isoprene at the start of a six day interval can enhance later ozone concentrations by as much as twenty-nine percent. The test cases used no input fluxes beyond the initial time, implying that a single input of a biogenic hydrocarbon to an airmass can alter its ozone chemistry over a time scale on the order of a week.
A parameterization of NAT (nitric acid trihydrate) clouds is developed for use in 2D models of the stratosphere. The parameterization uses model distributions of HNO3 and H2O to determine critical temperatures for NAT formation as a function of latitude and pressure. National Meteorological Center temperature fields are then used to determine monthly temperature frequency distributions, also as a function of latitude and pressure. The fractions of these distributions which fall below the critical temperatures for NAT formation are then used to determine the NAT cloud surface area density for each location in the model grid. By specifying heterogeneous reaction rates as functions of the surface area density, it is then possible to assess the effects of the NAT clouds on model constituent distributions. We also consider the increase in the NAT cloud formation in the presence of a fleet of stratospheric aircraft. The stratospheric aircraft NO(x) and H2O perturbations result in increased HNO3 as well as H2O. This increases the probability of NAT formation substantially, especially if it is assumed that the aircraft perturbations are confined to a corridor region.
A fundamental problem in the study of the terrestrial middle atmosphere is to calculate accurately the local heating due to the absorption of solar radiation. Knowledge of the heat budget is essential to understanding the atmospheric thermal structure, atmospheric motions, atmospheric chemistry, and their coupling. The evaluation of heating rates is complicated (especially above the stratopause) by the fact that the heating is not a simple one-step process. That is, the absorbed solar energy does not all immediately appear as heat. Rather, substantial portions of the incident energy may appear as internal energy of excited photolysis products (e.g., O(1D) or O2(1 delta)) or as chemical potential energy of product species such as atomic oxygen. The ultimate disposition of the internal and chemical energy possessed by the photolysis products determines the efficiency and thus the rate at which the middle atmosphere is heated. In studies of the heat budget, it is also vitally important to consider transport of long lived chemical species such as atomic oxygen above approximately 80 km. In such cases, the chemical potential energy may be transported great distances (horizontally or vertically) before undergoing a reaction to release the heat. Atomic oxygen influences the heating not only by reactions with itself and with O2 but also by reactions with odd-hydrogen species, especially those involving OH (Mlynczak and Solomon, 1991a). Consequently, absorbed solar energy may finally by converted to heat a long time after and at a location far from the original deposition. The purpose of this paper is to examine the solar and chemical heating processes and to present parameterizations for the heating efficiencies readily applicable for use in numerical models and heat budget studies. In the next two sections the processes relevant to the heating efficiencies for ozone and molecular oxygen will be reviewed. In section 4 the processes for the exothermic reactions will be reviewed and parameterizations for the heating efficiencies for both the solar and chemical processes will be presented in Section 5.
During the year of 1994, the engineering problems in aircraft modeling were studied. The initial concern was to obtain a surface model with desirable geometric characteristics. Much of the effort during the first half of the year was to find an efficient way of solving a computationally difficult optimization model. Since the smoothing technique in the proposal 'Surface Modeling and Optimization Studies of Aerodynamic Configurations' requires solutions of a sequence of large-scale quadratic programming problems, it is important to design algorithms that can solve each quadratic program in a few interactions. This research led to three papers by Dr. W. Li, which were submitted to SIAM Journal on Optimization and Mathematical Programming. Two of these papers have been accepted for publication. Even though significant progress has been made during this phase of research and computation times was reduced from 30 min. to 2 min. for a sample problem, it was not good enough for on-line processing of digitized data points. After discussion with Dr. Robert E. Smith Jr., it was decided not to enforce shape constraints in order in order to simplify the model. As a consequence, P. Dierckx's nonparametric spline fitting approach was adopted, where one has only one control parameter for the fitting process - the error tolerance. At the same time the surface modeling software developed by Imageware was tested. Research indicated a substantially improved fitting of digitalized data points can be achieved if a proper parameterization of the spline surface is chosen. A winning strategy is to incorporate Dierckx's surface fitting with a natural parameterization for aircraft parts. The report consists of 4 chapters. Chapter 1 provides an overview of reverse engineering related to aircraft modeling and some preliminary findings of the effort in the second half of the year. Chapters 2-4 are the research results by Dr. W. Li on penalty functions and conjugate gradient methods for quadratic programming problems.
The evolution of synoptic-scale dynamics associated with a middle and upper tropospheric cloud event that occurred on 26 November 1991 is examined. The case under consideration occurred during the FIRE CIRRUS-II Intensive Field Observing Period held in Coffeyville, KS during Nov. and Dec., 1991. Using data from the wind profiler demonstration network and a temporally and spatially augmented radiosonde array, emphasis is given to explaining the evolution of the kinematically-derived ageostrophic vertical circulations and correlating the circulation with the forcing of an extensively sampled cloud field. This is facilitated by decomposing the horizontal divergence into its component parts through a natural coordinate representation of the flow. Ageostrophic vertical circulations are inferred and compared to the circulation forcing arising from geostrophic confluence and shearing deformation derived from the Sawyer-Eliassen Equation. It is found that a thermodynamically indirect vertical circulation existed in association with a jet streak exit region. The circulation was displaced to the cyclonic side of the jet axis due to the orientation of the jet exit between a deepening diffluent trough and building ridge. The cloud line formed in the ascending branch of the vertical circulation with the most concentrated cloud development occurring in conjunction with the maximum large-scale vertical motion. The relationship between the large scale dynamics and the parameterization of middle and upper tropospheric clouds in large-scale models is discussed and an example of ice water contents derived from a parameterization forced by the diagnosed vertical motions and observed water vapor contents is presented.
We investigate the thermal history of an Earth-like planet with the help of a parameterized mantle convection model including the volatile exchange between mantle and surface reservoirs. The weakening of mantle silicates by dissolved volatiles is described by a functional relationship between creep rate and water fugacity. We use flow law parameters of diffusion creep in olivine under dry and wet conditions. The mantle degassing rate is considered as directly proportional to the seafloor spreading rate, which is also dependent on the mantle heat flow. To calculate the spreading rate, we assume that the heat flow under the mid-ocean ridges is double the average mantle heat flow. The rate of regassing also depends on the seafloor spreading rate as well as on other factors like the efficiency of volatile recycling through island arc volcanism. Both mechanisms (de- and regassing) are coupled self-consistently with the help of the parameterized convection model under implementation of a temperature and volatile-content-dependent mantle viscosity.
The evolution of synoptic-scale dynamics associated with a middle and upper tropospheric cloud event that occurred on 26 November 1991 is examined. The case under consideration occurred during the FIRE CIRRUS-II Intensive Field Observing Period held in Coffeyville, KS during Nov. and Dec., 1991. Using data from the wind profiler demonstration network and a temporally and spatially augmented radiosonde array, emphasis is given to explaining the evolution of the kinematically-derived ageostrophic vertical circulations and correlating the circulation with the forcing of an extensively sampled cloud field. This is facilitated by decomposing the horizontal divergence into its component parts through a natural coordinate representation of the flow. Ageostrophic vertical circulations are inferred and compared to the circulation forcing arising from geostrophic confluence and shearing deformation derived from the Sawyer-Eliassen Equation. It is found that a thermodynamically indirect vertical circulation existed in association with a jet streak exit region. The circulation was displaced to the cyclonic side of the jet axis due to the orientation of the jet exit between a deepening diffluent trough and building ridge. The cloud line formed in the ascending branch of the vertical circulation with the most concentrated cloud development occurring in conjunction with the maximum large-scale vertical motion. The relationship between the large scale dynamics and the parameterization of middle and upper tropospheric clouds in large-scale models is discussed and an example of ice water contents derived from a parameterization forced by the diagnosed vertical motions and observed water vapor contents is presented.
In the surface hydrologic parameterization of general circulation models (GCMs), it is commonly assumed that the precipitation processes are homogeneous over a GCM grid square and that the precipitation intensity is uniformly distributed. Based on evidence that the spatial distribution of precipitation within a GCM grid square is crucial for the land surface hydrology parameterization, a few researchers have explored the impacts of assuming that the precipitation is exponentially distributed. This paper explores the suitability of the aforementioned assumptions. First, a statistical analysis is conducted of historical precipitation data for three GCM grids in different regions of the United States. The analysis suggests that neither the uniform nor the exponential distribution assumption may be suitable at the GCM grid scale and, that instead, the spatial variability in precipitation is characterized by statistical patterns that are inhomogeneous. These patterns vary from grid to grid and are induced by the interaction between atmospheric conditions and various land surface characteristics, such as topographical features, surface properties, etc. Within the same grid square, however, the statistical patterns are generally constant from year to year. Based on this analysis, a computationally viable (i.e., usable with GCMs) stochastic precipitation disaggregation scheme that utilizes these stable statistical patterns is proposed. The method was used to generate spatially distributed hourly rainfall for a summer season in the southwestern region of the continental United States. Analysis of the results shows that the methodology preserves the seasonal characteristics of spatial variability in precipitation that is observed in the long-term historical data.
The formation, viscous evolution, and gravitational stability of protoplanetary disks are investigated. The formation process is parameterized by the angular velocity of the molecular cloud core omega, while the viscous evolution is parameterized by the viscosity parameter alpha in the disk; in this study we consider a range of (0.4-6) x 10(exp -14)/s for omega and from 10(exp -5) to 10(exp -1) for alpha. The axisymmetric gravitational stabilities of the disks are checked using Toomre's criterion. The resulting disk surface temperature distribution, (d log T(sub s)/d log R) approximately = -0.6 (R is the cylindrical radius), can be attributed to two heating sources: the viscous heating dominant in the inner disk region, and the accretion shock heating dominant in the outer disk region. This surface temperature distribution matches that observed in many disks around young stellar objects. During the infall stage, disks with alpha less than 10(exp -1.5) become gravitationally unstable independent of omega. The gravitational instabilities occur at radii ranging from 5 to 40 AU. The ratio of the disk mass to the central star mass ranges from 0.2 to 0.5 at the times of instability, about 4 x 10(exp -5) x (omega/10(exp -14)/s)(exp -0.67) yr. Most disks with low alpha and high omega become gravitationally unstable during their formation phase.
Two independent datasets for the solar radiation at the surface derived from satellites are compared. The data derived from the Earth Radiation Budget Experiment (ERBE) is for the net solar radiation at the surface whereas the International Satellite Cloud Climatology Project (ISCCP) data is for the downward flux only and was corrected with a space- and time-varying albedo. The ISCCP net flux is at all times higher than the ERBE flux. The difference can be divided into an offset that decreases with latitude and another component that correlates with high tropical cloud cover. With this latter exception the two datasets provide spatial patterns of solar flux that are very similar. A tropical Pacific Ocean model is forced with these two datasets and observed climatological winds. The upward heat flux is parameterized taking into account separately the longwave radiative, latent, and sensible heat fluxes. Best fit values for the uncertain parameters are found using an optimization procedure that seeks to minimize the difference between model and observed SST by varying the parameters within a reasonable range of uncertainty. The SST field the model produces with the best fit parameters is the best the model can do. If the differences between the model and data are larger than can be accounted for by remaining uncertainties in the heat flux parameterization and forcing data then the ocean model must be held to be at fault. Using this method of analysis, a fundamental model fault is identified. Inadequate treatment of mixed layer/entrainment processes in upwelling regions of the eastern tropical Pacific leads to a large and seasonally varying error in the model SST. Elsewhere the model SST is insufficiently different from observed to be able to identify model errors.
The one-dimensional radiative-convective equilibrium model with an explicit hydrologic cycle introduced in part 1 is used to study the sensitivity of the model's atmosphere to large changes in the solar forcing, when various cumulus convection parameterizations are used. When the concentration of the absorbing as in the atmosphere is temperature dependent, equilibrium is impossible for values of the solar forcing larger than a critical value. This result is referred to as a runaway greenhouse. The cumulus convection parameterization schemes currently in use in global climate models (GCMs) employ different assumptions about moistening. This causes the critical solar forcing above which a runaway greenhouse occurs to be very sensitive to the cumulus convection scheme employed. Furthermore, we show that the sensitivity of the equilibrium temperature to changes in the solar forcing depends crucially on the microphysics of cumulus convection. For fixed cloud conditions, the critical forcing for a runaway greenhouse to occur is between approximately 1.22 and 1.49 times the global mean value for the Earth, and for clear sky conditions, it is a few percent lower. The runaway greenhouse in the experiments with the mass flux schemes generally occurs more rapidly than in the experiments with the adjustment schemes. In addition, the inability of the hard convective adjustment scheme to produce an efficient vertical transport of moisture, together with the saturation requirement for convection to occur, leads to the breakdown of the radiative-convective equilibria when other processes are not available to provide the necessary vertical transport of water vapor.
Generated by landscape discontinuities (e.g., sea breezes) mesoscale circulation processes are not represented in large-scale atmospheric models (e.g., general circulation models), which have an inappropiate grid-scale resolution. With the assumption that atmospheric variables can be separated into large scale, mesoscale, and turbulent scale, a set of prognostic equations applicable in large-scale atmospheric models for momentum, temperature, moisture, and any other gaseous or aerosol material, which includes both mesoscale and turbulent fluxes is developed. Prognostic equations are also developed for these mesoscale fluxes, which indicate a closure problem and, therefore, require a parameterization. For this purpose, the mean mesoscale kinetic energy (MKE) per unit of mass is used, defined as E-tilde = 0.5 (the mean value of u'(sub i exp 2), where u'(sub i) represents the three Cartesian components of a mesoscale circulation (the angle bracket symbol is the grid-scale, horizontal averaging operator in the large-scale model, and a tilde indicates a corresponding large-scale mean value). A prognostic equation is developed for E-tilde, and an analysis of the different terms of this equation indicates that the mesoscale vertical heat flux, the mesoscale pressure correlation, and the interaction between turbulence and mesoscale perturbations are the major terms that affect the time tendency of E-tilde. A-state-of-the-art mesoscale atmospheric model is used to investigate the relationship between MKE, landscape discontinuities (as characterized by the spatial distribution of heat fluxes at the earth's surface), and mesoscale sensible and latent heat fluxes in the atmosphere. MKE is compared with turbulence kinetic energy to illustrate the importance of mesoscale processes as compared to turbulent processes. This analysis emphasizes the potential use of MKE to bridge between landscape discontinuities and mesoscale fluxes and, therefore, to parameterize mesoscale fluxes generated by such subgrid-scale landscape discontinuities in large-scale atmospheric models.
We examine the electromagnetic (EM) bias by using retracked TOPEX altimeter data. In contrast to previous studies, we use a parameterization of the EM bias which does not make stringent assumptions about the form of the correction or its global behavior. We find that the most effective single parameter correction uses the altimeter-estimated wind speed but that other parameterizations, using a wave age related parameter of significant wave height, may also significantly reduce the repeat pass variance. The different corrections are compared, and their improvement of the TOPEX height variance is quantified.
The wind speed and significant wave height (H(sub 1/3)) dependencies of the sea state bias in altimeter estimates of sea level, expressed in the form (Delta)h(sub SSB) = bH(sub 1/3), are examined from least squares analysis of 21 cycles of collinear TOPEX data. The bias coefficient b is found to increase in magnitude with increasing wind speed up to about 12 m/s and decrease monotonically in magnitude with increasing H(sub 1/3). A parameterization of b as a quadratic function of wind speed only, as in the formation used to produce the TOPEX geophysical data records (GDRs), is significantly better than a parameterization purely in terms of H(sub 1/3). However, a four-parameter combined wind speed and wave height formulation for b (quadratic in wind speed plus linear in H(sub 1/3)) significantly improves the accuracy of the sea state bias correction. The GDR formulation in terms of wind speed only should therefore be expanded to account for a wave height dependence of b. An attempt to quantify the accuracy of the sea state bias correction (Delta)h(sub SSB) concludes that the uncertainty is a disconcertingly large 1% of H(sub 1/3).
Formulas are presented that parameterize the heating rate and coefficient of turbulent heat conduction produced by saturated internal gravity waves (IGW) in the upper atmosphere. Estimates of these values are made using observational data. The parameterization of IGW influences are introduced into a one-dimensional model of global mean thermal and composition balances of the upper atmosphere. Computations are performed for different values of IGW energy fluxes entering into the upper atmosphere from below. It is shown that realistic vertical profiles of the global mean temperature can be obtained using different values of IGW energy flux into the upper atmosphere. Increasing the IGW intensity leads not only to an increase of the heating rate due to wave enery dissipation, but also to an increase of the heating rate due to wave energy dissipation, but also to an increase in the coefficient of turbulent heat conduction and cooling rate produced by turbulence generated by the wave. So, near an altitude of 100 km the main part of solar heating is compensated by infrared cooling on one hand, and the main part of wave dissipation heating is compensated by turbulent cooling on the other hand. These quasi-balances generally hold for different values of IGW intensity.
At the last International Reference Ionosphere (IRI) Workshop it was decided that future editions of the model should include a representationof the aurroral oval boundaries. In this paper we review the different existing parameterizations of the auroral oval discussing their data base, boundary criteria, matematical formation, and overall usefulness for IRI. As a first candidate for incorporation into IRI we recomment the parameterization of the Feldstein ovals by Holzworth and Meng. Ways of implementing this model into IRI are discussed. We will also address adjustability with user-provided boundaries or boundary-related parameters, to better support strom-related studies.