Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parameterization method”

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

A stochastic precipitation disaggregation scheme for GCM applications

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.

Gao, Xiaogang↗
NASA NTRS2004

An Empirical Cumulus Parameterization Scheme for a Global Spectral Model

Realistic vertical heating and drying profiles in a cumulus scheme is important for obtaining accurate weather forecasts. A new empirical cumulus parameterization scheme based on a procedure to improve the vertical distribution of heating and moistening over the tropics is developed. The empirical cumulus parameterization scheme (ECPS) utilizes profiles of Tropical Rainfall Measuring Mission (TRMM) based heating and moistening derived from the European Centre for Medium- Range Weather Forecasts (ECMWF) analysis. A dimension reduction technique through rotated principal component analysis (RPCA) is performed on the vertical profiles of heating (Q1) and drying (Q2) over the convective regions of the tropics, to obtain the dominant modes of variability. Analysis suggests that most of the variance associated with the observed profiles can be explained by retaining the first three modes. The ECPS then applies a statistical approach in which Q1 and Q2 are expressed as a linear combination of the first three dominant principal components which distinctly explain variance in the troposphere as a function of the prevalent large-scale dynamics. The principal component (PC) score which quantifies the contribution of each PC to the corresponding loading profile is estimated through a multiple screening regression method which yields the PC score as a function of the large-scale variables. The profiles of Q1 and Q2 thus obtained are found to match well with the observed profiles. The impact of the ECPS is investigated in a series of short range (1-3 day) prediction experiments using the Florida State University global spectral model (FSUGSM, T126L14). Comparisons between short range ECPS forecasts and those with the modified Kuo scheme show a very marked improvement in the skill in ECPS forecasts. This improvement in the forecast skill with ECPS emphasizes the importance of incorporating realistic vertical distributions of heating and drying in the model cumulus scheme. This also suggests that in the absence of explicit models for convection, the proposed statistical scheme improves the modeling of the vertical distribution of heating and moistening in areas of deep convection.

Rajendran, K.↗
NASA NTRS2012

Estimation of the Relationship Between Remotely Sensed Anthropogenic Heat Discharge and Building Energy Use

This paper examined the relationship between remotely sensed anthropogenic heat discharge and energy use from residential and commercial buildings across multiple scales in the city of Indianapolis, Indiana, USA. The anthropogenic heat discharge was estimated with a remote sensing-based surface energy balance model, which was parameterized using land cover, land surface temperature, albedo, and meteorological data. The building energy use was estimated using a GIS-based building energy simulation model in conjunction with Department of Energy/Energy Information Administration survey data, the Assessor's parcel data, GIS floor areas data, and remote sensing-derived building height data. The spatial patterns of anthropogenic heat discharge and energy use from residential and commercial buildings were analyzed and compared. Quantitative relationships were evaluated across multiple scales from pixel aggregation to census block. The results indicate that anthropogenic heat discharge is consistent with building energy use in terms of the spatial pattern, and that building energy use accounts for a significant fraction of anthropogenic heat discharge. The research also implies that the relationship between anthropogenic heat discharge and building energy use is scale-dependent. The simultaneous estimation of anthropogenic heat discharge and building energy use via two independent methods improves the understanding of the surface energy balance in an urban landscape. The anthropogenic heat discharge derived from remote sensing and meteorological data may be able to serve as a spatial distribution proxy for spatially-resolved building energy use, and even for fossil-fuel CO2 emissions if additional factors are considered.

Multi-scale↗
NASA NTRS1996

A Revised Thermosphere for the Mars Global Reference Atmospheric Model (Mars-GRAM Version 3.4)

This report describes the newly-revised model thermosphere for the Mars Global Reference Atmospheric Model (Mars-GRAM, Version 3.4). It also provides descriptions of other changes made to the program since publication of the programmer's guide for Mars-GRAM Version 3.34. The original Mars-GRAM model thermosphere was based on the global-mean model of Stewart. The revised thermosphere is based largely on parameterizations derived from output data from the three-dimensional Mars Thermospheric Global Circulation Model (MTGCM). The new thermospheric model includes revised dependence on the 10.7 cm solar flux for the global means of exospheric temperature, temperature of the base of the thermosphere, and scale height for the thermospheric temperature variations, as well as revised dependence on orbital position for global mean height of the base of the thermosphere. Other features of the new thermospheric model are: (1) realistic variations of temperature and density with latitude and time of day, (2) more realistic wind magnitudes, based on improved estimates of horizontal pressure gradients, and (3) allowance for user-input adjustments to the model values for mean exospheric temperature and for height and temperature at the base of the thermosphere. Other new features of Mars-GRAM 3.4 include: (1) allowance for user-input values of climatic adjustment factors for temperature profiles from the surface to 75 km, and (2) a revised method for computing the sub-solar longitude position in the 'ORBIT' subroutine.

Justus, C. G.↗
NASA NTRS1983

Data processing techniques used with MST radars: A review

The data processing methods used in high power radar probing of the middle atmosphere are examined. The radar acts as a spatial filter on the small scale refractivity fluctuations in the medium. The characteristics of the received signals are related to the statistical properties of these fluctuations. A functional outline of the components of a radar system is given. Most computation intensive tasks are carried out by the processor. The processor computes a statistical function of the received signals, simultaneously for a large number of ranges. The slow fading of atmospheric signals is used to reduce the data input rate to the processor by coherent integration. The inherent range resolution of the radar experiments can be improved significant with the use of pseudonoise phase codes to modulate the transmitted pulses and a corresponding decoding operation on the received signals. Commutability of the decoding and coherent integration operations is used to obtain a significant reduction in computations. The limitations of the processors are outlined. At the next level of data reduction, the measured function is parameterized by a few spectral moments that can be related to physical processes in the medium. The problems encountered in estimating the spectral moments in the presence of strong ground clutter, external interference, and noise are discussed. The graphical and statistical analysis of the inferred parameters are outlined. The requirements for special purpose processors for MST radars are discussed.

Rastogi, P. K.↗
NASA NTRSDate not supplied

Exploring Quantum State Preparation Using Tensor Networks and Sparse Wavefunction Simulations

The variational quantum eigenvalue solver is a powerful hybrid quantum-classical approach that has been suggested as a candidate method to run on near-term quantum hardware for computing ground state electronic energies of molecular systems. However, even for small molecules, the number of variational parameters and qubits required to minimize the electronic energy is beyond the reach of current quantum computers except for small basis sets. We explore a new paradigm for state preparation where we test how much of the optimization can be approximately prepared with classical computers to reduce the number of optimization steps performed using a quantum device. By adapting a recent algorithm for the factorized form of the UCC ansatz, we can study molecular electronic structure problems with up to 64 qubits. In addition, we also test a related approach of using tensor networks to optimize quantum circuits in order to benchmark various lattice models. We present results using these approaches and discuss strategies for incorporating these ideas into variational algorithms involving near-term quantum computers. Our results help demonstrate the strength of the UCC ansatz and address pressing questions about optimal initial parameterizations and circuit construction.

quantum computing↗
NASA NTRSDate not supplied

Remote Sensing of Lineage Functional Types for Modeling and Monitoring Biodiversity

Hyperspectral remote sensing has the potential to continuously scale plant function and plant diversity information from landscape to global extents. Numerous studies have indicated that VSWIR (400-2500 nm) reflectance properties of vegetation capture evolutionarily conserved biochemical, structural, and other functional attributes of plant species. Spectral properties conserved in plants provide the opportunity to both 1) aggregate species into lineages with improved classification accuracy and 2) link those lineages directly to plant traits. Full realization of this goal will enable parameterization of Land Surface Models (LSMs) with remotely sensed information, e.g., canopy nitrogen, and better representations of biodiversity and functional diversity in biogeographic studies. In this study, we use hyperspectral AVIRIS data from the 2013 HyspIRI campaign over the Southern Sierra Nevada, California flight box to investigate the potential for incorporating evolutionary thinking into landcover classification. We link the airborne hyperspectral data with vegetation plot data from roughly 1372 surveys and a phylogeny representing 1361 species. We aggregate species into lineages ranging from species level groups down to similar number of Plant Functional Types as often used in LSMs. We assessed the ability of Random Forest and Partial Least Squares Discriminant Analysis to discriminate across these different phylogenetic scales and determine the optimal number of lineages to classify. Although there are some temporal and spatial differences in our training data, our best approaches achieved moderate classification accuracy (Kappa > 0.65). Given an optimal number of lineages, we explored approaches to improve classifications including machine learning and unmixing approaches. This work suggests that lineage-based methods may be a promising way to leverage the huge amounts of data that will come from high resolution and high return interval hyperspectral data planned for the Surface Biology and Geology mission with sparsely sampled existing ground-based ecological data.

Hyperspectral↗
NASA NTRS2023

Development of High-Resolution Soil Hydraulic Parameters with Use of Earth Observations for Enhancing Root Zone Soil Moisture Product

Regional quantification of energy and water balance fluxes depends inevitably on the estimation of surface and rootzone soil moisture. The simulation of soil moisture depends on the soil retention characteristics, which are difficult to estimate at a regional scale. Thus, the present study proposes a new method to estimate high-resolution Soil Hydraulic Parameters (SHPs) which in turn help to provide high-resolution (spatial and temporal) rootzone soil moisture (RZSM) products. The study is divided into three phases—(I) involves the estimation of finer surface soil moisture (1 km) from the coarse resolution satellite soil moisture. The algorithm utilizes MODIS 1 km Land Surface Temperature (LST) and 1 km Normalized difference vegetation Index (NDVI) for downscaling 25 km C-band derived soil moisture from AMSR-2 to 1 km surface soil moisture product. At one of the test sites, soil moisture is continuously monitored at 5, 20, and 50 cm depth, while at 44 test sites data were collected randomly for validation. The temporal and spatial correlation for the downscaled product was 70% and 83%, respectively. (II) In the second phase, downscaled soil moisture product is utilized to inversely estimate the SHPs for the van Genuchten model (1980) at 1 km resolution. The numerical experiments were conducted to understand the impact of homogeneous SHPs as compared to the three-layered parameterization of the soil profile. It was seen that the SHPs estimated using the downscaled soil moisture (I-d experiment) performed with similar efficiency as compared to SHPs estimated from the in-situ soil moisture data (I-b experiment) in simulating the soil moisture. The normalized root mean square error (nRMSE) for the two treatments was 0.37 and 0.34, respectively. It was also noted that nRMSE for the treatment with the utilization of default SHPs (I-a) and AMSR-2 soil moisture (I-c) were found to be 0.50 and 0.43, respectively. (III) Finally, the derived SHPs were used to simulate both surface soil moisture and RZSM. The final product, RZSM which is the daily 1 km product also showed a nearly 80% correlation at the test site. The estimated SHPs are seen to improve the mean NSE from 0.10 (I-a experiment) to 0.50 (I-d experiment) for the surface soil moisture simulation. The mean nRMSE for the same was found to improve from 0.50 to 0.31.

soil moisture downscaling↗
NASA NTRS2007

Using Remotely-Sensed Estimates of Soil Moisture to Infer Soil Texture and Hydraulic Properties across a Semi-arid Watershed

Near-surface soil moisture is a critical component of land surface energy and water balance studies encompassing a wide range of disciplines. However, the processes of infiltration, runoff, and evapotranspiration in the vadose zone of the soil are not easy to quantify or predict because of the difficulty in accurately representing soil texture and hydraulic properties in land surface models. This study approaches the problem of parameterizing soils from a unique perspective based on components originally developed for operational estimation of soil moisture for mobility assessments. Estimates of near-surface soil moisture derived from passive (L-band) microwave remote sensing were acquired on six dates during the Monsoon '90 experiment in southeastern Arizona, and used to calibrate hydraulic properties in an offline land surface model and infer information on the soil conditions of the region. Specifically, a robust parameter estimation tool (PEST) was used to calibrate the Noah land surface model and run at very high spatial resolution across the Walnut Gulch Experimental Watershed. Errors in simulated versus observed soil moisture were minimized by adjusting the soil texture, which in turn controls the hydraulic properties through the use of pedotransfer functions. By estimating a continuous range of widely applicable soil properties such as sand, silt, and clay percentages rather than applying rigid soil texture classes, lookup tables, or large parameter sets as in previous studies, the physical accuracy and consistency of the resulting soils could then be assessed. In addition, the sensitivity of this calibration method to the number and timing of microwave retrievals is determined in relation to the temporal patterns in precipitation and soil drying. The resultant soil properties were applied to an extended time period demonstrating the improvement in simulated soil moisture over that using default or county-level soil parameters. The methodology is also applied to an independent case at Walnut Gulch using a new soil moisture product from active (C-band) radar imagery with much lower spatial and temporal resolution. Overall, results demonstrate the potential to gain physically meaningful soils information using simple parameter estimation with few but appropriately timed remote sensing retrievals.

Santanello, Joseph A.↗
NASA NTRS2018

Aerodynamic Design and Optimization of Fan Stage for Boundary Layer Ingestion Propulsion System

The present paper addresses the process of preliminary design of a low-pressure fan and outlet guide vane (OGV) of a boundary layer ingestion (BLI) propulsion system. The tail-cone thruster systems of NASA's STARC_ABL (Single-aisle Turboelectric Aircraft with an Aft Boundary-Layer propulsor) adopts an axi-symmetric BLI type inlet as opposed to other embedded engine systems. Thus, the focus of the present work is placed on maximizing the efficiency of the fan and OGV stages under a significant radial distortion. A parameterization with B-spline function for camber line angles, metal chord, thickness distribution and stacking axis of blades is presented. The flowpath lines are also parameterized by B-spline function and aggregated in the design system of blades. The design optimization with evolutionary algorithm is performed with constraints of fan pressure ratio, OGV exit swirl angle and nozzle exit properties. The inlet conditions for the turbo-machinery CFD (Computational Fluid Dynamics) domain and the design goal of the fan stage are driven by a propulsion airframe integration (PAI) model that uses a 3-D unstructured RANS (Reynolds Average Navier Stokes) solver and actuator disk model. The expected power saving of the BLI propulsor is quantified via PAI analysis and the resulting preliminary design of the fan stages is compared with a clean-inlet flow propulsor.

Turbo-machinery↗
NASA NTRS2016

Connecting Paleo and Modern Oceanographic Data to Understand Atlantic Meridional Overturning Circulation Over Decades to Centuries

Modeling is an important tool for understanding AMOC on all timescales. Mechanistic studies of modern AMOC variability have been hampered by a lack of consistency between free-running models and the sensitivity of AMOC to resolution and parameterization. Recent work within the framework of the phase two Coordinated Ocean- Reference Experiments (CORE-II) addresses this issue head on, looking at model differences of AMOC mean state and interannual variability. One consistent feature across the models is that AMOC mean transport is related to mixed layer depths and Labrador Sea salt content, whereas interannual variability is primarily associated with Labrador Sea temperature anomalies. This is consistent with the hypothesized importance of salt balance for AMOC variability on geological timescales. The simulated relationships between AMOC and subsurface temperature anomalies in fully coupled climate models reveal subsurface AMOC fingerprints that could be used to reconstruct historical AMOC variations at low frequency.With the lack of long-term AMOC observations, models of ocean state that assimilate observational data have been explored as a way to reconstruct AMOC, but comparisons between models indicate they are quite variable in their AMOC representations. Karspeck et al. (2015) found that historical reconstructions of AMOC in such models are sensitive to the details of the data assimilation procedure. The ocean data assimilation community continues to address these issues through improved models and methods for estimating and representing error information.Two objectives of paleoclimate modeling are 1) to provide mechanistic information for interpretation of paleoclimate observations, and 2) to test the ability of predictive models to simulate Earth's climate under different background forcing states. In a good example of the first objective, Schmittner and Lund (2015) and Menviel et al. (2014) provided key information about the proxy signals expected under freshwater disturbance of AMOC, which were used to support the paleoclimate observations made by Henry et al. (2016). In an example of the second objective, Muglia and Schmittner (2015) analyzed Third Paleoclimate Modeling Intercomparison Project (PMIP3) models of the Last Glacial Maximum (LGM) and found consistently more intense and deeper AMOC transports relative to preindustrial simulations, counter to the paleoclimate consensus of LGM conditions, indicating that some processes are not well represented in the PMIP3 models. One challenge is to find adequate paleo observations against which to test these models. PMIP is now in phase 4 (part of CMIP6), which includes experiments covering five periods in Earth's history: the last millennium, last glacial maximum, last interglacial, and the mid-Pliocene. Newly compiled paleoclimate datasets from the PAGES2k project, more transient simulations, and participation of isotope enabled models planned for CMIP6PMIP4 will enable richer paleo data-model comparisons in the near future.

Kilbourne, Hali↗
NASA NTRS2010

Parameterized Linear Longitudinal Airship Model

A parameterized linear mathematical model of the longitudinal dynamics of an airship is undergoing development. This model is intended to be used in designing control systems for future airships that would operate in the atmospheres of Earth and remote planets. Heretofore, the development of linearized models of the longitudinal dynamics of airships has been costly in that it has been necessary to perform extensive flight testing and to use system-identification techniques to construct models that fit the flight-test data. The present model is a generic one that can be relatively easily specialized to approximate the dynamics of specific airships at specific operating points, without need for further system identification, and with significantly less flight testing. The approach taken in the present development is to merge the linearized dynamical equations of an airship with techniques for estimation of aircraft stability derivatives, and to thereby make it possible to construct a linearized dynamical model of the longitudinal dynamics of a specific airship from geometric and aerodynamic data pertaining to that airship. (It is also planned to develop a model of the lateral dynamics by use of the same methods.) All of the aerodynamic data needed to construct the model of a specific airship can be obtained from wind-tunnel testing and computational fluid dynamics

Kulczycki, Eric↗
NASA NTRS2014

Phase Transitions in Planning Problems: Design and Analysis of Parameterized Families of Hard Planning Problems

There are two common ways to evaluate algorithms: performance on benchmark problems derived from real applications and analysis of performance on parametrized families of problems. The two approaches complement each other, each having its advantages and disadvantages. The planning community has concentrated on the first approach, with few ways of generating parametrized families of hard problems known prior to this work. Our group's main interest is in comparing approaches to solving planning problems using a novel type of computational device - a quantum annealer - to existing state-of-the-art planning algorithms. Because only small-scale quantum annealers are available, we must compare on small problem sizes. Small problems are primarily useful for comparison only if they are instances of parametrized families of problems for which scaling analysis can be done. In this technical report, we discuss our approach to the generation of hard planning problems from classes of well-studied NP-complete problems that map naturally to planning problems or to aspects of planning problems that many practical planning problems share. These problem classes exhibit a phase transition between easy-to-solve and easy-to-show-unsolvable planning problems. The parametrized families of hard planning problems lie at the phase transition. The exponential scaling of hardness with problem size is apparent in these families even at very small problem sizes, thus enabling us to characterize even very small problems as hard. The families we developed will prove generally useful to the planning community in analyzing the performance of planning algorithms, providing a complementary approach to existing evaluation methods. We illustrate the hardness of these problems and their scaling with results on four state-of-the-art planners, observing significant differences between these planners on these problem families. Finally, we describe two general, and quite different, mappings of planning problems to QUBOs, the form of input required for a quantum annealing machine such as the D-Wave II.

Problems↗
NASA NTRS2017

Remote Sensing of Multiple Cloud Layer Heights Using Multi-Angular Measurements

Cloud top height (CTH) affects the radiative properties of clouds. Improved CTH observations will allow for improved parameterizations in large-scale models and accurate information on CTH is also important when studying variations in freezing point and cloud microphysics. NASAs airborne Research Scanning Polarimeter (RSP) is able to measure cloud top height using a novel multi-angular contrast approach. For the determination of CTH, a set of consecutive nadir reflectances is selected and the cross-correlations between this set and co-located sets at other viewing angles are calculated for a range of assumed cloud top heights, yielding a correlation profile. Under the assumption that cloud reflectances are isotropic, local peaks in the correlation profile indicate cloud layers. This technique can be applied to every RSP footprint and we demonstrate that detection of multiple peaks in the correlation profile allow retrieval of heights of multiple cloud layers within single RSP footprints. This paper provides an in-depth description of the architecture and performance of the RSPs CTH retrieval technique using data obtained during the Studies of Emissions and Atmospheric Composition, Clouds and Climate Coupling by Regional Surveys (SEAC(exp. 4)RS) campaign. RSP retrieved cloud heights are evaluated using collocated data from the Cloud Physics Lidar (CPL). The method's accuracy associated with the magnitude of correlation, optical thickness, cloud thickness and cloud height are explored. The technique is applied to measurements at a wavelength of 670 nm and 1880 nm and their combination. The 1880-nm band is virtually insensitive to the lower troposphere due to strong water vapor absorption.

Sinclair, Kenneth↗
NASA NTRS1999

Practical Aerodynamic Design Optimization Based on the Navier-Stokes Equations and a Discrete Adjoint Method

Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. An efficient surface parameterization based on a free-form deformation technique has been utilized and the resulting codes have been integrated with an optimization package. Lastly, sample optimizations have been shown for inviscid and turbulent flow over an ONERA M6 wing. Drag reductions have been demonstrated by reducing shock strengths across the span of the wing. In order for large scale optimization to become routine, the benefits of parallel architectures should be exploited. Although the flow solver has been parallelized using compiler directives. The parallel efficiency is under 50 percent. Clearly, parallel versions of the codes will have an immediate impact on the ability to design realistic configurations on fine meshes, and this effort is currently underway.

Grossman, Bernard↗
NASA NTRS1999

Practical Aerodynamic Design Optimization Based on the Navier-Stokes Equations and a Discrete Adjoint Method

The technical details are summarized below: Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. . An efficient surface parameterization based on a free-form deformation technique has been utilized and the resulting codes have been integrated with an optimization package. Lastly, sample optimizations have been shown for inviscid and turbulent flow over an ONERA M6 wing. Drag reductions have been demonstrated by reducing shock strengths across the span of the wing.

Grossman, Bernard↗
NASA NTRS2009

Application of TRMM PR and TMI Measurements to Assess Cloud Microphysical Schemes in the MM5 Model for a Winter Storm

Seen by the human eye, precipitation particles are commonly drops of rain, flakes of snow, or lumps of hail that reach the ground. Remote sensors and numerical models usually deal with information about large collections of rain, snow, and hail (or graupel --also called soft hail ) in a volume of air. Therefore, the size and number of the precipitation particles and how particles interact, evolve, and fall within the volume of air need to be represented using physical laws and mathematical tools, which are often implemented as cloud and precipitation microphysical parameterizations in numerical models. To account for the complexity of the precipitation physical processes, scientists have developed various types of such schemes in models. The accuracy of numerical weather forecasting may vary dramatically when different types of these schemes are employed. Therefore, systematic evaluations of cloud and precipitation schemes are of great importance for improvement of weather forecasts. This study is one such endeavor; it pursues quantitative assessment of all the available cloud and precipitation microphysical schemes in a weather model (MM5) through comparison with the observations obtained by National Aeronautics and Space Administration (NASA) s and Japan Aerospace Exploration Agency (JAXA) s Tropical Rainfall Measuring Mission (TRMM) precipitation radar (PR) and microwave imager (TMI). When satellite sensors (like PR or TMI) detect information from precipitation particles, they cannot directly observe the microphysical quantities (e.g., water species phase, density, size, and amount etc.). Instead, they tell how much radiation is absorbed by rain, reflected away from the sensor by snow or graupel, or reflected back to the satellite. On the other hand, the microphysical quantities in the model are usually well represented in microphysical schemes and can be converted to radiative properties that can be directly compared to the corresponding PR and TMI observations. This study employs this method to evaluate the accuracy of the simulated radiative properties by the MM5 model with different microphysical schemes. It is found that the representations of particle density, size, and mass in the different schemes in the MM5 model determine the model s performance when predicting a winter storm over the eastern Pacific Ocean. Schemes lacking moderate density particles (i.e. graupel), with snow flakes that are too large, or with excessive mass of snow or graupel lead to degraded prediction of the radiative properties as observed by the TRMM satellite. This study demonstrates the uniqueness of the combination of both an active microwave sensor (PR) and passive microwave sensor (TMI) onboard TRMM on assessing the accuracy of numerical weather forecasting. It improves our understanding of the physical and radiative properties of different types of precipitation particles and provides suggestions for better representation of cloud and precipitation processes in numerical models. It would, ultimately, contribute to answering questions like "Why did it not rain when the forecast says it would?"

Han, Mei↗
NASA NTRS1994

Improvements to TOVS retrievals over sea ice and applications to estimating Arctic energy fluxes

Modeling studies suggest that polar regions play a major role in modulating the Earth's climate and that they may be more sensitive than lower latitudes to climate change. Until recently, however, data from meteorological stations poleward of 70 degs have been sparse, and consequently, our understanding of air-sea-ice interaction processes is relatively poor. Satellite-borne sensors now offer a promising opportunity to observe polar regions and ultimately to improve parameterizations of energy transfer processes in climate models. This study focuses on the application of the TIROS-N operational vertical sounder (TOVS) to sea-ice-covered regions in the nonmelt season. TOVS radiances are processed with the improved initialization inversion ('3I') algorithm, providng estimates of layer-average temperature and moisture, cloud conditions, and surface characteristics at a horizontal resolution of approximately 100 km x 100 km. Although TOVS has flown continuously on polar-orbiting satellites since 1978, its potential has not been realized in high latitudes because the quality of retrievals is often significantly lower over sea ice and snow than over the surfaces. The recent availability of three Arctic data sets has provided an opportunity to validate TOVS retrievals: the first from the Coordinated Eastern Arctic Experiment (CEAREX) in winter 1988/1989, the second from the LeadEx field program in spring 1992, and the third from Russian drifting ice stations. Comparisons with these data reveal deficiencies in TOVS retrievals over sea ice during the cold season; e.g., ice surface temperature is often 5 to 15 K too warm, microwave emissivity is approximately 15% too low at large view angles, clear/cloudy scenes are sometimes misidentified, and low-level inversions are often not captured. In this study, methods to reduce these errors are investigated. Improvements to the ice surface temperature retrieval have reduced rms errors from approximately 7 K to 3 K; correction of microwaves (50 GHz) brightness temperatures for view angle dependence now allows the surface type (open water versus sea ice) to be determined; modifications to the clear/cloud tests have improved cloud detection over sea ice, especially where it is inhomogeneous; and improved surface temperature estimates have resulted in more successful diagnoses of low-level stratification. Applications of improved retrievals to studies of the Arctic energy budget are encouraging. Preliminary calculations of air-ice stress vectors, 10-m wind speed, and horizontal advection of heat and moisture are presented.

Francis, Jennifer A.↗