SEARCH · Search NASA
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.
The BFS method for alloys is applied to the study of Co growth on Cu(111). The parameterization of the Co-Cu system is obtained from first-principles calculations, and tested against known experimental features for low coverage Co deposition on Cu(100) and Cu(111). Atomistic simulations are performed to investigate the behavior of Co on Cu(111) as a function of coverage.
Mars accretion is known to be fast compared to Earth. Basaltic samples provide a probe into the interior and allow reconstruction of siderophile element contents of the mantle. These estimates can be used to estimate conditions of core formation, as for Earth. Although many assume that Mars went through a magma ocean stage, and possibly even complete melting, the siderophile element content of Mars mantle is consistent with relatively low pressure and temperature (PT) conditions, implying only shallow melting, near 7 GPa and 2073 K. This is a pressure range where some have proposed a change in siderophile element partitioning behavior. We will examine the databases used for parameterization and split them into a low and higher pressure regime to see if the methods used to reach this conclusion agree for the two sets of data.
Among the many possible norm-based optimization methods, the concept of H-infinity optimal control has gained enormous attention in the past few years. Here the H-infinity framework, based on the Small Gain Theorem and the Youla Parameterization, effectively treats system uncertainties in the control law synthesis. A design approach involving a mixed H(sub 2)/H-infinity norm strives to combine the advantages of both methods. This advantage motivates researchers toward finding solutions to the mixed H(sub 2)/H-infinity control problem. The approach developed in this research is based on a finite time cost functional that depicts an H-infinity bound control problem in a H(sub 2)-optimization setting. The goal is to define a time-domain cost function that optimizes the H(sub 2)-norm of a system with an H-infinity-constraint function.
Among the many possible norm-based optimization methods, the concept of H-infinity optimal control has gained enormous attention in the past few years. Here the H-infinity framework, based on the Small Gain Theorem and the Youla Parameterization, effectively treats system uncertainties in the control law synthesis. A design approach involving a mixed H(sub 2)/H-infinity norm strives to combine the advantages of both methods. This advantage motivates researchers toward finding solutions to the mixed H(sub 2)/H-infinity control problem. The approach developed in this research is based on a finite time cost functional that depicts an H-infinity bound control problem in a H(sub 2)-optimization setting. The goal is to define a time-domain cost function that optimizes the H(sub 2)-norm of a system with an H-infinity-constraint function.
The gap in explicit resolution of phenomena between global climate models and cloud resolving models is shrinking at a steady pace with global integrations of several km in resolution now practical for at least some time scales. In moving toward finer resolution the long standing problem of convective parameterization is being examined along with the assumption of convective quasi-equilibrium and how this can be reconciled with the stochastic and intermittent nature of convection. In this context we examine the nature of parameterized convection in the NASA Goddard Earth Observing System (GEOS-5) Atmospheric General Circulation Model. Our analysis uses both coarse (2.5 degree) and fine scale 0.25 degree spatial resolution integrations. Two basic formulations are compared: the default option is the Relaxed Arakawa-Schubert (RAS) scheme which invokes a sequence of linearly entraining plumes and quasi-equilibrium closure. An optional modification of this method (the "Stochastic Tokioka" constraint) places a random lower limit on plume entrainment. An alternative representation is the Kain-Fritsch parameterization which was originally developed for mesoscale numerical modeling strategies. Here entrainment is determined by a crude buoyancy sorting approach that allows the plume spectrum to be more responsive to ambient vertical stratification of moisture. Diagnostics of the model behavior are referenced to recent observational evidence of continuous phase transition behavior. In particular we examine the relationship between column water vapor and probablility of convective presence and intensity. Sensitivity of the statistics of convective behavior to parcel mixing/entrainment formulations and parcel initial thermodynamics are considered. Observational statistics from A-Train and TRMM sensors provide validation of the model integrations.
The renormalization group (RG) procedure for nonlinear, dissipative systems is now quite standard, and its applications to the problem of hydrodynamic turbulence are becoming well known. In summary, the RG method isolates self similar behavior and provides a systematic procedure to describe scale invariant dynamics in terms of large scale variables only. The parameterization of the small scales in a self consistent manner has important implications for sub-grid modeling. This paper develops the homogeneous, isotropic turbulence and addresses the meaning and consequence of epsilon-expansion. The theory is then extended to include a weak mean flow and application of the RG method to a sequence of models is shown to converge to the Navier-Stokes equations.
The World Climate Research Programme (WCRP) plan is concerned with the need to develop a uniform global cloud climatology as part of a broad research program on climate processes. The International Satellite Cloud Climatology Project (ISCCP) has been approved as the first project of the WCRP. The ISCCP has the basic objective to collect and analyze satellite radiance data to infer the global distribution of cloud radiative properties in order to improve the modeling of cloud effects on climate. Research is conducted to explore an algorithm for retrieving cloud properties by utilizing the available infrared sounder data from polar-orbiting satellites. A numerical method is developed for computing cloud top heights, amount, and emissivity on the basis of a parameterized infrared radiative transfer equation for cloudy atmospheres. Theoretical studies were carried out by considering a synthetic atmosphere.
A summary is presented of the main results obtained during the course of research on output feedback control. The term output feedback is used to denote a controller design approach which does not rely on an observer to estimate the states of the system. Thus, the order of the controller is fixed, and can even be zero order, which amounts to constant gain ouput feedback. The emphasis has been on optimal output feedback. That is, a fixed order controller is designed based on minimizing a suitably chosen quadratic performance index. A number of problem areas that arise in this context have been addressed. These include developing suitable methods for selecting an index of performance, both time domain and frequency domain methods for achieving robustness of the closed loop system, developing canonical forms to achieve a minimal parameterization for the controller, two time scale design formulations for ill-conditioned systems, and the development of convergent numerical algorithms for solving the output feedback problem.
Regional climate models are expected to exhibit improved skill at finer spatial resolutions due to improved representation of land surface heterogeneity. However, at spatial scales between 1 to 10 km (grey scales), these improvements are often illusive due to the competing benefits from spatial resolution and cumulus parameterization. This study provides insights into the impact of model resolution and cumulus parameterization on precipitation prediction in the Central Great Plains by using an object-based evaluation method. Our results show limited improvement solely from finer resolution but larger improvement without using the cumulus scheme at a 4-kmresolution. Compared to traditional evaluation methods, the object-based analysis shows that without the cumulus scheme the spatial properties of precipitation are better represented. In contrast, all model configurations show a dry bias in precipitation days and a tendency to produce widespread precipitation but with fewer hours with precipitationwhich indicates other shortcomings in the model.
A large set of deterministic and ensemble forecasts was produced to identify the optimal spacing for forecasting U.S. East Coast snowstorms. WRF forecasts were produced on cloud-allowing (~1-km grid spacing) and convection-allowing (3–4 km) grids, and compared against forecasts with parameterized convection (>~10 km). Performance diagrams were used to evaluate 19 deterministic forecasts from the winter of 2013–14. Ensemble forecasts of five disruptive snowstorms spanning the years 2015–18 were evaluated using various methods to evaluate probabilistic forecasts. While deterministic forecasts using cloud-allowing grids were not better than convection-allowing forecasts, both had lower bias and higher success ratios than forecasts with parameterized convection. All forecasts were underdispersive. Nevertheless, forecasts on the higher-resolution grids were more reliable than those with parameterized convection. Forecasts on the cloud-allowing grid were best able to discriminate areas that received heavy snow and those that did not, while the forecasts with parameterized convection were least able to do so. It is recommended to use convection-resolving and (if computationally possible) to use cloud-allowing forecast grids when predicting East Coast winter storms.
The effects of nonspherical small ice crystals on the IR radiative properties of cirrus clouds are investigated utilizing light scattering properties of spheroidal particles. Employing the anomalous diffraction theory for spheres and results from the exact spheroid scattering program, efficient parameterization equations are developed for calculations of the absorption and scattering properties for small ice crystals. Parameterization formulas are developed for large ice crystals employing results computed from the geometric ray-tracing method and the Fraunhofer diffraction theory for hexagonal crystals and spheroids.
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.
Evapotranspiration (ET) is critical in linking global water, carbon and energy cycles. However, direct measurement of global terrestrial ET is not feasible. Here, we first reviewed the basic theory and state-of-the-art approaches for estimating global terrestrial ET, including remote-sensing-based physical models, machine-learning algorithms and land surface models (LSMs). We then utilized 4 remote-sensing-based physical models, 2 machine-learning algorithms and 14 LSMs to analyze the spatial and temporal variations in global terrestrial ET. The results showed that the ensemble means of annual global terrestrial ET estimated by these three categories of approaches agreed well, with values ranging from 589.6 mm/yr (6.56×10^4 cu.km/yr) to 617.1 mm/yr (6.87×10^4 cu.km/yr). For the period from 1982 to 2011, both the ensembles of remote-sensing-based physical models and machine-learning algorithms suggested increasing trends in global terrestrial ET (0.62 mm/sq.yr with a significance level of p<0.05 and 0.38 mm yr−2 with a significance level of p<0.05, respectively). In contrast, the ensemble mean of the LSMs showed no statistically significant change (0.23 mm/sq.yr, p>0.05), although many of the individual LSMs reproduced an increasing trend. Nevertheless, all 20 models used in this study showed that anthropogenic Earth greening had a positive role in increasing terrestrial ET. The concurrent small interannual variability, i.e., relative stability, found in all estimates of global terrestrial ET, suggests that a potential planetary boundary exists in regulating global terrestrial ET, with the value of this boundary being around 600 mm/yr. Uncertainties among approaches were identified in specific regions, particularly in the Amazon Basin and arid/semiarid regions. Improvements in parameterizing water stress and canopy dynamics, the utilization of new available satellite retrievals and deep-learning methods, and model–data fusion will advance our predictive understanding of global terrestrial ET.
A method is presented for the unique derivation of both starspot temperatures and effective areas from the range of standardized V and R light curves obtained through the rotation periods, and the method is applied to spots observed on II Pegasi, BY Draconis, and HD 209813. The method involves the determination of the temperature difference between the spot and the surrounding photosphere from V and V-R photometry and the Barnes-Evans visual surface brightness relation, and area determination by the derivation of a parameterized effective area representing the fractional area of disk covered by spots under the assumption that the spot or spot group is adequately described by a single circular spot. By this method, an effective temperature of 3400 K and an area representing 37% of the facing hemisphere are derived for the large, cool spot observed on II Pegasi during 1977, an effective temperature of 3500 K and area of 34% are derived for BY Dra and an effective temperature of less than 3840 K is obtained for HD 209813. The derived starspot temperatures, which are believed to correspond to the starspot penumbrae, imply umbral temperatures of 2600-2900 K, in excellent agreement with the 3 kilogauss umbral spot models of Mullan (1974). The geometric modelling of the formation and migration of the spot on BY Dra also indicates a differential rotation amounting to -9.7 x 10 to the -9th rad/sec per deg, in agreement with the solar value.
Continuing efforts toward validating the buildup factor method and the BRYNTRN code, which use the deterministic approach in solving radiation transport problems and are the candidate engineering tools in space radiation shielding analyses, are presented. A simplified theory of proton buildup factors assuming no neutron coupling is derived to verify a previously chosen form for parameterizing the dose conversion factor that includes the secondary particle buildup effect. Estimates of dose in tissue made by the two deterministic approaches and the Monte Carlo method are intercompared for cases with various thicknesses of shields and various types of proton spectra. The results are found to be in reasonable agreement but with some overestimation by the buildup factor method when the effect of neutron production in the shield is significant. Future improvement to include neutron coupling in the buildup factor theory is suggested to alleviate this shortcoming. Impressive agreement for individual components of doses, such as those from the secondaries and heavy particle recoils, are obtained between BRYNTRN and Monte Carlo results.
This paper explains why high resolution design spaces encourage traditional airfoil optimization algorithms to generate noisy shape modifications, which lead to inaccurate linear predictions of aerodynamic coefficients and potential failure of descent methods. By using auxiliary drag constraints for a simultaneous drag reduction at all design points and the least shape distortion to achieve the targeted drag reduction, an improved algorithm generates relatively smooth optimal airfoils with no severe off-design performance degradation over a range of flight conditions, in high resolution design spaces parameterized by cubic B-spline functions. Simulation results using FUN2D in Euler flows are included to show the capability of the robust aerodynamic shape optimization method over a range of flight conditions.
We solve the problem of determining airfoils that approximate, in a least square sense, given surface pressure distributions in transonic flight regimes. The flow is modeled by means of the Euler equations and the solution procedure is an adjoint- based minimization algorithm that makes use of the inverse Theodorsen transform in order to parameterize the airfoil. Fast convergence to the optimal solution is obtained by means of the pseudo-time method. Results are obtained using three different pressure distributions for several free stream conditions. The airfoils obtained have given a trailing edge angle.
A systematic development of the delta-four-stream approximation for calculations of radiative fluxes in planetary atmosphere is presented. It is shown that an analytic solution for this approximation can be derived explicitly, with minimum computational effort for flux calculations. Relative accuracy checks for reflection, transmission, and absorption for numerous asymmetry factors, single-scattering albedos, optical depths, and solar zenith angles have been performed with respect to the 'exact' results computed from the adding method for radiative transfer. Overall, results from the delta-four-stream approximation yield relative accuracies within about 5 percent. This approximation is well suited to radiative transfer parameterizations involving flux and heating calculations in aerosol and cloudy atmospheres.