Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variational Analysis”

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 91 records · Page 5

Assimilation of satellite surface wind speed data using the GLA analysis/forecast system

Research at the Goddard Laboratory for Atmospheres to evaluate methods to assign directions to both real and simulated surface wind speed data is discussed. Surface wind speed measurements are obtained from satellites, including Seasat, Nimbus-7, and Geosat. The methods include the interpolation of modeled forecast winds to the wind speed datum location, the use of surface pressures with a balance relation, and a variational analysis method. It is found that the best estimates may be obtained using a multipass approach to perform a surface analysis incorporating all conventional surface data.

Bloom, S. C.↗

The influence of atmospheric stratification on scatterometer data

The effects of atmospheric stratification and the stability of the atmospheric stratification on the scatterometer data measuring surface winds over the ocean were investigated using the boundary layer model developed by Louis (1979). A variational analysis method is proposed, which allows direct assimilation of scatterometer data. It is shown that the effect of the stability of atmospheric stratification on the wind increment is relatively small. However, it is a systematic effect, and neglecting it would consistently underestimate the winds in stable regions.

Louis, Jean-Francois↗

Space-based surface wind vectors to aid understanding of air-sea interactions

A novel and unique ocean-surface wind data-set has been derived by combining the Defense Meteorological Satellite Program Special Sensor Microwave Imager data with additional conventional data. The variational analysis used generates a gridded surface wind analysis that minimizes an objective function measuring the misfit of the analysis to the background, the data, and certain a priori constraints. In the present case, the European Center for Medium-Range Weather Forecasts surface-wind analysis is used as the background.

Atlas, R.↗

Microcracking, microcrack-induced delamination, and longitudinal splitting of advanced composite structures

A combined analytical and experimental study was conducted to analyze microcracking, microcrack-induced delamination, and longitudinal splitting in polymer matrix composites. Strain energy release rates, calculated by a variational analysis, were used in a failure criterion to predict microcracking. Predictions and test results were compared for static, fatigue, and cyclic thermal loading. The longitudinal splitting analysis accounted for the effects of fiber bridging. Test data are analyzed and compared for longitudinal splitting and delamination under mixed-mode loading. This study emphasizes the importance of using fracture mechanics analyses to understand the complex failure processes that govern composite strength and life.

Nairn, John A.↗

An adaptive identification and control scheme for large space structures

A unified identification and control scheme capable of achieving space at form performance objectives under nominal or failure conditions is described. Preliminary results are also presented, showing that the methodology offers much promise for effective robust control of large space structures. The control method is a multivariable, adaptive, output predictive controller called Model Predictive Control (MPC). MPC uses a state space model and input reference trajectories of set or tracking points to adaptively generate optimum commands. For a fixed model, MPC processes commands with great efficiency, and is also highly robust. A key feature of MPC is its ability to control either nonminimum phase or open loop unstable systems. As an output controller, MPC does not explicitly require full state feedback, as do most multivariable (e.g., Linear Quadratic) methods. Its features are very useful in LSS operations, as they allow non-collocated actuators and sensors. The identification scheme is based on canonical variate analysis (CVA) of input and output data. The CVA technique is particularly suited for the measurement and identification of structural dynamic processes - that is, unsteady transient or dynamically interacting processes such as between aerodynamics and structural deformation - from short, noisy data. CVA is structured so that the identification can be done in real or near real time, using computationally stable algorithms. Modeling LSS dynamics in 1-g laboratories has always been a major impediment not only to understanding their behavior in orbit, but also to controlling it. In cases where the theoretical model is not confirmed, current methods provide few clues concerning additional dynamical relationships that are not included in the theoretical models. CVA needs no a priori model data, or structure; all statistically significant dynamical states are determined using natural, entropy-based methods. Heretofore, a major limitation in applying adaptive controllers to LSS problems has been the requirement for accurate models. We feel that MPC coupled with the CVA identification scheme will resolve this limitation. The results are based on research performed under an SBIR grant, administered by the Air Force Office of Scientific Research.

Carroll, J. V.↗

Solar Cycle Predictions Near Sunspot Minimum

Observations of solar magnetic activity and the dynamics of the solar convection zone have produced severe constraints on models of the Sun's magnetic dynamo. These constraints are so severe that, at present, we do not have numerical models that can accept the current conditions and then march forward in time to predict future activity. Given this state of solar dynamo theory we are forced to examine previous behavior to discover patterns and trends that afford us some measure of predictability. Here we examine the behavior of several indicators of solar activity near solar minimum that are well correlated with the amplitude of the following solar maximum to predict the level of solar activity over cycle 23. Sunspot numbers, areas, and positions are useful for characterizing solar cycle behavior due to the extent of the data (12 cycles or more). These data exhibit several patterns that relate future activity to past behavior. With the Odd-Even effect the odd numbered cycles have been larger than their even numbered predecessors for each of the last six cycle pairs. With the Amplitude-Period effect short period cycles have been followed by large amplitude cycles and long period cycles have been followed by small amplitude cycles for 10 of the last 13 cycles. With the Maximum-Minimum effect the sunspot number at minimum is directly correlated with the sunspot number at maximum for a given cycle. The geomagnetic indices aa and Ap are also related to solar activity by the connections between disturbances in the solar wind and variations in the Earth's magnetic field. Like the Maximum-Minimum effect for sunspots, the size of the aa and Ap indices at minimum are directly related to the amplitude of the following maximum. The number of geomagnetically disturbed days (days with Ap >= 25) over the course of a cycle is another indicator for the size of the next cycle. The aa and Ap indices can each be separated into a component in phase with the current sunspot cycle and an "extended cycle" component that is directly related to the following cycle. Individually each of these indicators provides some measure of predictive ability for the amplitudes of previous cycles. Taken in combination we find that the predictions are much more accurate. Our multi-variate analysis indicates that a maximum sunspot number of 160 +/- 20 should occur near the end of 1999 for cycle 23.

Hathaway, David H.↗

Current and Future Plans of the NASA Data Assimilation Office (DAO)

The mission of the Data Assimilation Office (DAO) is to advance the state of the art of data assimilation and produce research-quality assimilated data sets which make optimal use of space-based observations. Development efforts over the last few years have focused on delivering a production data assimilation system in support of NASA's Terra platform. That system, called the Goddard Earth Observing System - version 2 or GEOS-2, represents a major upgrade to the baseline GEOS-1 system employed in NASA's first reanalysis effort. GEOS-2 includes a physical-space three dimensional variational analysis algorithm (the Physical-space Statistical Analysis System or PSAS) and numerous improvements to the general circulation model. The latter include a Soil-Vegetation-Atmosphere Transfer (SVAT) land surface scheme, a level 2.5 moist turbulence scheme and new Short Wave (SW) and Long Wave (LW) radiation code. The system also includes an off-line ozone assimilation system, and the capability to assimilate scatterometer surface winds, and TIROS Operational Vertical Sounder (TOVS) and Special Sensor Microwave Imager (SSM/I) moisture data. GEOS-2 is currently run at 1 degree horizontal resolution and 48 levels extending to O.Olmb. Experimental versions of GEOS-2 are run with a global stretched grid allowing enhanced (e.g. 1/4 deg) regional resolution. Other capabilities being developed include, the assimilation of Tropical Rainfall Measuring Mission (TRMM) precipitation and Global Positioning System (GPS) data, an off-line land surface assimilation system, and a retrospective analysis scheme. The DAO is also engaged in a number of collaborative efforts to help accelerate the development of the next generation data assimilation system. These include, a joint modeling effort between the DAO and NCAR/CGDD to develop a new Global Circulation Model (GCM), and a Department of Energy Lawrence Livermore National Laboratory (DOE/LLNL) collaboration on model parallelization. Plans for the next reanalysis will be discussed in the context of current and near term system quality and computing capabilities, and the need for multiple reanalysis products.

Atlas, Robert↗

Surface-Wind Anomalies in North-Atlantic and North Pacific from SSM/I Observations: Influence on Temperature of Adjoining Land Regions

Surface winds over the oceans are derived from Special Sensor Microwave Imager (SSM/I) measurements, assigning direction by Variational Analysis Method (VAM). Validations by comparison with other measurements indicate highly-satisfactory data quality. Providing global coverage from 1988, the dataset is a convenient source for surface-wind climatology. In this study, the interannual variability of zonal winds is analyzed concentrating on the westerlies in North Atlantic and North Pacific, above 30 N. Interannual differences in the westerlies exceeding 10 m sec (exp -1) are observed over large regions, often accompanied by changes of the same magnitude in the easterlies below 30 N. We concentrate on February/March, since elevated temperatures, by advancing snow-melt, can produce early spring. The extremely strong westerlies in 1997 observed in these months over North Atlantic (and also North Pacific) apparently contributed to large surface-temperature anomalies in western Europe, on the order of +3 C above the climatic monthly average for England and France. At these latitudes strong positive anomalies extended in a ring around the globe. We formulated an Index of South westerlies for the North Atlantic, which can serve as an indicator for day-by-day advection effects into Europe. In comparing 1997 and 1998 with the previous years, we establish significant correlations with the temperature anomalies (one to five days later, depending on the region, and on the season). This variability of the ocean-surface winds and of the temperature anomalies on land may be related to the El Nino/La Nina oscillations. Such large temperature fluctuations over large areas, whatever the cause, can be regarded as noise in attempts to assess long-term trends in global temperature.

Otterman, Joseph↗

The Special Sensor Microwave Imager Wind Dataset: A Source of Quantitative Information for the Ocean-to-Land Advection

It is well recognized that advection from the North Atlantic has a profound effect on the climatic conditions in central Europe. A new dataset of the ocean-surface winds, derived from the Special Sensor Microwave Imager, SSM/1, is now available. This satellite instrument measures the wind speed, but not the direction. However, variational analysis developed at the Data Assimilation Office, NASA Goddard Space Flight Center, by combining the SSM/I measurements with wind vectors measured from ships, etc., produced global maps of the ocean surface winds suitable for climate analysis. From this SSM/I dataset, a specific index I(sub na) of the North Atlantic surface winds has been developed, which pertinently quantifies the low-level advection into central Europe. For a selected time-period, the index I(sub na) reports the average of the amplitude of the wind, averaging only the speed when the direction is from the southwest (when the wind is from another direction, the contribution counts to the average as zero speed). Strong correlations were found between February I(sub na) and the surface air temperatures in Europe 50-60 deg N. In the present study, we present the correlations between I(sub na) and temperature I(sub s), and also the sensitivity of T(sub s), to an increase in I(sub na), in various seasons and various regions. We specifically analyze the flow of maritime-air from the North Atlantic that produced two extraordinary warm periods: February 1990, and early-winter 2000/2001. The very cold December 2001 was clearly due to a northerly flow. Our conclusion is that the SSM/I dataset is very useful for providing insight to the forcing of climatic fluctuations in Europe.

Otterman, J.↗

North Atlantic Basin Tropical Cyclone Activity in Relation to Temperature and Decadal- Length Oscillation Patterns

Yearly frequencies of North Atlantic basin tropical cyclones, their locations of origin, peak wind speeds, average peak wind speeds, lowest pressures, and average lowest pressures for the interval 1950-2008 are examined. The effects of El Nino and La Nina on the tropical cyclone parametric values are investigated. Yearly and 10-year moving average (10-yma) values of tropical cyclone parameters are compared against those of temperature and decadal-length oscillation, employing both linear and bi-variate analysis, and first differences in the 10-yma are determined. Discussion of the 2009 North Atlantic basin hurricane season, updating earlier results, is given.

Wilson, Robert M.↗

Using High-Resolution Forward Model Simulations of Ideal Atmospheric Tracers to Assess the Spatial Information Content of Inverse CO2 Flux Estimates

Attribution of observed atmospheric carbon concentrations to emissions on the country, state or city level is often inferred using "inversion" techniques. Such computations are often performed using advanced mathematical techniques, such as synthesis inversion or four-dimensional variational analysis, that invoke tracing observed atmospheric concentrations backwards through a transport model to a source region. It is, to date, not well understood how well such techniques can represent fine spatial (and temporal) structure in the inverted flux fields. This question is addressed using forward-model computations with idealized tracers emitted at the surface in a large number of grid boxes over selected regions and examining how distinctly these emitted tracers can be detected downstream. Initial results show that tracers emitted in half-degree grid boxes over a large region of the Eastern USA cannot be distinguished from each other, even at short distances over the Atlantic Ocean, when they are emitted in grid boxes separated by less than five degrees of latitude - especially when only total-column observations are available. A large number of forward model simulations, with varying meteorological conditions, are used to assess how distinctly three types observations (total column, upper tropospheric column, and surface mixing ratio) can separate emissions from different sources. Inferences inverse modeling and source attribution will be drawn.

Pawson, Steven↗

RACORO Continental Boundary Layer Cloud Investigations: 1. Case Study Development and Ensemble Large-Scale Forcings

Observation-based modeling case studies of continental boundary layer clouds have been developed to study cloudy boundary layers, aerosol influences upon them, and their representation in cloud- and global-scale models. Three 60 h case study periods span the temporal evolution of cumulus, stratiform, and drizzling boundary layer cloud systems, representing mixed and transitional states rather than idealized or canonical cases. Based on in situ measurements from the Routine AAF (Atmospheric Radiation Measurement (ARM) Aerial Facility) CLOWD (Clouds with Low Optical Water Depth) Optical Radiative Observations (RACORO) field campaign and remote sensing observations, the cases are designed with a modular configuration to simplify use in large-eddy simulations (LES) and single-column models. Aircraft measurements of aerosol number size distribution are fit to lognormal functions for concise representation in models. Values of the aerosol hygroscopicity parameter, kappa, are derived from observations to be approximately 0.10, which are lower than the 0.3 typical over continents and suggestive of a large aerosol organic fraction. Ensemble large-scale forcing data sets are derived from the ARM variational analysis, European Centre for Medium-Range Weather Forecasts, and a multiscale data assimilation system. The forcings are assessed through comparison of measured bulk atmospheric and cloud properties to those computed in "trial" large-eddy simulations, where more efficient run times are enabled through modest reductions in grid resolution and domain size compared to the full-sized LES grid. Simulations capture many of the general features observed, but the state-of-the-art forcings were limited at representing details of cloud onset, and tight gradients and high-resolution transients of importance. Methods for improving the initial conditions and forcings are discussed. The cases developed are available to the general modeling community for studying continental boundary clouds.

Boundary layers↗

Decoding the effects of synonymous variants

Synonymous single nucleotide variants (sSNVs) are common in the human genome but are often overlooked. However, sSNVs can have significant biological impact and may lead to disease. Existing computational methods for evaluating the effect of sSNVs suffer from the lack of gold-standard training/evaluation data and exhibit over-reliance on sequence conservation signals. We developed synVep (synonymous Variant effect predictor), a machine learning-based method that overcomes both of these limitations. Our training data was a combination of variants reported by gnomAD (observed) and those unreported, but possible in the human genome (generated). We used positive-unlabeled learning to purify the generated variant set of any likely unobservable variants. We then trained two sequential extreme gradient boosting models to identify subsets of the remaining variants putatively enriched and depleted in effect. Our method attained 90% precision/recall on a previously unseen set of variants. Furthermore, although synVep does not explicitly use conservation, its scores correlated with evolutionary distances between orthologs in cross-species variation analysis. synVep was also able to differentiate pathogenic vs. benign variants, as well as splice-site disrupting variants (SDV) vs. non-SDVs. Thus, synVep provides an important improvement in annotation of sSNVs, allowing users to focus on variants that most likely harbor effects.

Zishuo Zeng↗

The Gmao Hybrid 4d-Envar Observing System Simulation Experiment Framework.

This work describes the extension of the Global Modeling and Assimilation Office (GMAO) Observing System Simulation Experiment (OSSE) framework to use a hybrid 4D-EnVar scheme instead of 3D-Var. The original 3D-Var and hybrid 4D-EnVar OSSEs use the same version of the data assimilation system (DAS) so that a direct comparison is possible in terms of the validation with respect to their corresponding real cases. Rather than quantifying the differences between the two data assimilation methodologies, a short inter-comparison of upgrading from a 3D- to a 4D-OSSE is provided to highlight aspects where this change matters to the OSSE community and to identify features of data assimilation that can only be explored in a four-dimensional OSSE framework. A short validation of the hybrid 4D-EnVar OSSE shows that conclusions from previous assessments of the 3D-Var OSSE in its ability to mimic the behavior of the real system still hold with the same caveats. Furthermore, some aspects of the ensemble configuration and behavior are discussed along with forecast sensitivity to observation impacts (FSOI). Estimates of error standard deviations are shown to be smaller in the hybrid 4D-EnVar OSSE but with little impact on the character of the error. A discussion on future work directions focuses on exploring the four-dimensional aspect such as the error distribution within the assimilation window or four-dimensional handling of high-temporal density observations.

Data assimilation↗

A multivariate variational objective analysis-assimilation method. Part 1: Development of the basic model

The variational method of undetermined multipliers is used to derive a multivariate model for objective analysis. The model is intended for the assimilation of 3-D fields of rawinsonde height, temperature and wind, and mean level temperature observed by satellite into a dynamically consistent data set. Relative measurement errors are taken into account. The dynamic equations are the two nonlinear horizontal momentum equations, the hydrostatic equation, and an integrated continuity equation. The model Euler-Lagrange equations are eleven linear and/or nonlinear partial differential and/or algebraic equations. A cyclical solution sequence is described. Other model features include a nonlinear terrain-following vertical coordinate that eliminates truncation error in the pressure gradient terms of the horizontal momentum equations and easily accommodates satellite observed mean layer temperatures in the middle and upper troposphere. A projection of the pressure gradient onto equivalent pressure surfaces removes most of the adverse impacts of the lower coordinate surface on the variational adjustment.

Achtemeier, Gary L.↗

Variational Methods in Sensitivity Analysis and Optimization for Aerodynamic Applications

Variational methods (VM) sensitivity analysis, which is the continuous alternative to the discrete sensitivity analysis, is employed to derive the costate (adjoint) equations, the transversality conditions, and the functional sensitivity derivatives. In the derivation of the sensitivity equations, the variational methods use the generalized calculus of variations, in which the variable boundary is considered as the design function. The converged solution of the state equations together with the converged solution of the costate equations are integrated along the domain boundary to uniquely determine the functional sensitivity derivatives with respect to the design function. The determination of the sensitivity derivatives of the performance index or functional entails the coupled solutions of the state and costate equations. As the stable and converged numerical solution of the costate equations with their boundary conditions are a priori unknown, numerical stability analysis is performed on both the state and costate equations. Thereafter, based on the amplification factors obtained by solving the generalized eigenvalue equations, the stability behavior of the costate equations is discussed and compared with the state (Euler) equations. The stability analysis of the costate equations suggests that the converged and stable solution of the costate equation is possible only if the computational domain of the costate equations is transformed to take into account the reverse flow nature of the costate equations. The application of the variational methods to aerodynamic shape optimization problems is demonstrated for internal flow problems at supersonic Mach number range. The study shows, that while maintaining the accuracy of the functional sensitivity derivatives within the reasonable range for engineering prediction purposes, the variational methods show a substantial gain in computational efficiency, i.e., computer time and memory, when compared with the finite difference sensitivity analysis.

Ibrahim, A. H.↗

A multivariate variational objective analysis-assimilation method. Part 2: Case study results with and without satellite data

The variational multivariate assimilation method described in a companion paper by Achtemeier and Ochs is applied to conventional and conventional plus satellite data. Ground-based and space-based meteorological data are weighted according to the respective measurement errors and blended into a data set that is a solution of numerical forms of the two nonlinear horizontal momentum equations, the hydrostatic equation, and an integrated continuity equation for a dry atmosphere. The analyses serve first, to evaluate the accuracy of the model, and second to contrast the analyses with and without satellite data. Evaluation criteria measure the extent to which: (1) the assimilated fields satisfy the dynamical constraints, (2) the assimilated fields depart from the observations, and (3) the assimilated fields are judged to be realistic through pattern analysis. The last criterion requires that the signs, magnitudes, and patterns of the hypersensitive vertical velocity and local tendencies of the horizontal velocity components be physically consistent with respect to the larger scale weather systems.

Achtemeier, Gary L.↗