Search NASA⌕ Search

SEARCH · Search NASA

Results for “Generalized Linear Model”

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 523 records · Page 29

Linear dimensionality of Landsat agricultural data with implications for classification

A model for the Landsat multispectral scanner data, representing a generalization of the commonly used Gaussian model, has been formulated and analyzed. The model hypothesizes that the data for different crop types essentially lie on distinct hyperplanes in the feature space. Tests of this model reveal that: (1) the agricultural data from any single acquisition (i.e., four-channel) of Landsat are essentially two dimensional, regardless of the crop type; and (2) the data from different sites and different stages of crop development all lie on planes which are parallel. These findings have significant implications for data display, classification, feature extraction, and signature extension.

Wheeler, S. G.↗

Estimation of Unsteady Aerodynamic Models from Dynamic Wind Tunnel Data

Demanding aerodynamic modelling requirements for military and civilian aircraft have motivated researchers to improve computational and experimental techniques and to pursue closer collaboration in these areas. Model identification and validation techniques are key components for this research. This paper presents mathematical model structures and identification techniques that have been used successfully to model more general aerodynamic behaviours in single-degree-of-freedom dynamic testing. Model parameters, characterizing aerodynamic properties, are estimated using linear and nonlinear regression methods in both time and frequency domains. Steps in identification including model structure determination, parameter estimation, and model validation, are addressed in this paper with examples using data from one-degree-of-freedom dynamic wind tunnel and water tunnel experiments. These techniques offer a methodology for expanding the utility of computational methods in application to flight dynamics, stability, and control problems. Since flight test is not always an option for early model validation, time history comparisons are commonly made between computational and experimental results and model adequacy is inferred by corroborating results. An extension is offered to this conventional approach where more general model parameter estimates and their standard errors are compared.

Murphy, Patrick↗

Crustal Strain Rates in the Western United States and Their Relationship with Earthquake Rates

Abstract We present a suite of strain rate models for the western United States based on geologic and geodetic data. The geologic data consist of Quaternary fault-slip rates and the geodetic data consists of a new compilation of Global Positioning System (GPS) velocities derived from continuous, semicontinuous, and campaign measurements. We remove postseismic deformation from the GPS time series in order for our geodetic strain rate model to best capture the interseismic strain accumulation rate. We present models based on either geologic or geodetic data, but also create a hybrid model. Although there are some differences between the models, the large-scale features are the same, with the noticeable exception for the Pacific Northwest where interseismic strain is naturally more distributed than the long-term strain release. We also present a map of earthquake rate densities based on mainshocks, and the result has similar spatial features similar to the strain rate models (at least in the southwestern United States). We perform a general correlation analysis between strain rate and seismicity rate (south of Cascadia) and find a change in linearity between seismicity and strain rates from slow to faster deforming areas with seismicity rates relatively lower for the latter. The extent of that change depends a bit on assumptions made on the declustering and completeness of the catalog, but the finding of a change in slope is robust across the different strain rate models. Linearity for all areas is only expected when Gutenberg–Richter parameters and parameters involved in the conversion from strain to moment rate are uniform across the study area. We discuss these qualifications, but find no single satisfactory explanation for our observation. Moreover, when considering a rather short time and space, theoretical considerations of sampling from a power-law distribution actually predict there to be a power law instead of a linear relationship, generally consistent with our observation.

Geochemistry & Geophysics↗

Parallel transport dynamics for mixed quantum states with applications to time-dependent density functional theory

Direct simulation of the von Neumann dynamics for a general (pure or mixed) quantum state can often be expensive. One prominent example is the real-time time-dependent density functional theory (rt-TDDFT), a widely used framework for the first principle description of many-electron dynamics in chemical and materials systems. Practical rt-TDDFT calculations often avoid the direct simulation of the von Neumann equation, and solve instead a set of Schrödinger equations, of which the dynamics is equivalent to that of the von Neumann equation. However, the time step size employed by the Schrödinger dynamics is often much smaller. Here, in order to improve the time step size and the overall efficiency of the simulation, we generalize a recent work of the parallel transport (PT) dynamics for simulating pure states [An, Lin, Multiscale Model. Simul. 18, 612, 2020] to general quantum states. The PT dynamics provides the optimal gauge choice, and can employ a time step size comparable to that of the von Neumann dynamics. Going beyond the linear and near adiabatic regime in previous studies, we find that the error of the PT dynamics can be bounded by certain commutators between Hamiltonians, density matrices, and their derived quantities. Such a commutator structure is not present in the Schrödinger dynamics. We demonstrate that the parallel transport-implicit midpoint (PT-IM) method is a suitable method for simulating the PT dynamics, especially when the spectral radius of the Hamiltonian is large. The commutator structure of the error bound, and numerical results for model rt-TDDFT calculations in both linear and nonlinear regimes, confirm the advantage of the PT dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

MULTI-DIMENSIONAL RADIATIVE TRANSFER TO ANALYZE HANLE EFFECT IN Ca II K LINE AT 3933 A

Radiative transfer (RT) studies of the linearly polarized spectrum of the Sun (the second solar spectrum) have generally focused on line formation, with an aim to understand the vertical structure of the solar atmosphere using one-dimensional (1D) model atmospheres. Modeling spatial structuring in the observations of the linearly polarized line profiles requires the solution of multi-dimensional (multi-D) polarized RT equation and a model solar atmosphere obtained by magnetohydrodynamical (MHD) simulations of the solar atmosphere. Our aim in this paper is to analyze the chromospheric resonance line Ca II K at 3933 A using multi-D polarized RT with the Hanle effect and partial frequency redistribution (PRD) in line scattering. We use an atmosphere that is constructed by a two-dimensional snapshot of the three-dimensional MHD simulations of the solar photosphere, combined with columns of a 1D atmosphere in the chromosphere. This paper represents the first application of polarized multi-D RT to explore the chromospheric lines using multi-D MHD atmospheres, with PRD as the line scattering mechanism. We find that the horizontal inhomogeneities caused by MHD in the lower layers of the atmosphere are responsible for strong spatial inhomogeneities in the wings of the linear polarization profiles, while the use of horizontally homogeneous chromosphere (FALC) produces spatially homogeneous linear polarization in the line core. The introduction of different magnetic field configurations modifies the line core polarization through the Hanle effect and can cause spatial inhomogeneities in the line core. A comparison of our theoretical profiles with the observations of this line shows that the MHD structuring in the photosphere is sufficient to reproduce the line wings and in the line core, but only line center polarization can be reproduced using the Hanle effect. For a simultaneous modeling of the line wings and the line core (including the line center), MHD atmospheres with inhomogeneities in the chromosphere are required.

79 ASTRONOMY AND ASTROPHYSICS↗

A Formally Verified Conflict Detection Algorithm for Polynomial Trajectories

In air traffic management, conflict detection algorithms are used to determine whether or not aircraft are predicted to lose horizontal and vertical separation minima within a time interval assuming a trajectory model. In the case of linear trajectories, conflict detection algorithms have been proposed that are both sound, i.e., they detect all conflicts, and complete, i.e., they do not present false alarms. In general, for arbitrary nonlinear trajectory models, it is possible to define detection algorithms that are either sound or complete, but not both. This paper considers the case of nonlinear aircraft trajectory models based on polynomial functions. In particular, it proposes a conflict detection algorithm that precisely determines whether, given a lookahead time, two aircraft flying polynomial trajectories are in conflict. That is, it has been formally verified that, assuming that the aircraft trajectories are modeled as polynomial functions, the proposed algorithm is both sound and complete.

Narkawicz, Anthony↗

Sequence-to-sequence neural networks for short-term electrical load forecasting in commercial office buildings

The U.S. power grid is transforming to become smarter, cleaner, and more effi- cient. This is leading to the addition of significant distributed variable renew- able generation. Due to the variable nature of renewable generation, the short- and long-term supply-demand imbalances are less predictable, and conventional approaches to mitigating the imbalance will not be efficient or cost-effective. To address this challenge, transactive control technologies have been proposed which balance energy generation and consumption with market activity and in- frastructural limitations. Transactive control requires the ability of individual end-use loads to express flexibility as a function of a transactive signal (e.g., price). Empirical gray- and black-box models have been widely used to express flexibility, and although these approaches are generally easy to construct and simple to use, they do not capture the non-linear behavior that some end-use loads represent . Machine learning approaches have been proposed to address this limitation. Although deep learning approaches for forecasting end-use loads have been explored, certain aspects of the application of deep models to load forecasting are not well understood. These aspects include how much training data is required, and how models should be structured and trained. To that end, this work explores how to approach applying deep recurrent neural networks to short-term electrical load forecasting with a case study of four commercial office buildings. We identify data requirements for training accurate models of whole building electricity use conditioned on outdoor temperature, provide insight into model hyperparameter sensitivity, and demonstrate how readily models can be generalized to unseen buildings.

Skomski, Elliott↗

An interpretation and solution of ill-conditioned linear equations

Data insufficiency, poorly conditioned matrices and singularities in equations occur regularly in complex optimization, correlation, and interdisciplinary model studies. This work concerns itself with two methods of obtaining certain physically realistic solutions to ill-conditioned or singular algebraic systems of linear equations arising from such studies. Two efficient computational solution procedures that generally lead to locally unique solutions are presented when there is insufficient data to completely define the model, or a least-squares error formulation of this system results in an ill-conditioned system of equations. If it is assumed that a reasonable estimate of the uncertain data is available in both cases cited above, then we shall show how to obtain realistic solutions efficiently, in spite of the insufficiency of independent data. The proposed methods of solution are more efficient than singular-value decomposition for dealing with such systems, since they do not require solutions for all the non-zero eigenvalues of the coefficient matrix.

Ojalvo, I. U.↗

Sensitivity analysis of radiative heating and cooling rates in planetary atmospheres: general linearization and adjoint approaches

Radiative heating and cooling provide primary source and ultimate sink of energy driving lower planetary atmospheres. Evaluating the sensitivities of atmospheric dynamics models on these primary atmospheric parameters requires knowing how heating and cooling rates depend on these same parameters. We discuss two approaches that make it possible to directly compute the sensitivities of heating and cooling rates in parallel with evaluation of heating and cooling rates themselves.

radiative transfer sensitivity analysis↗

Generalized Relationship Linking Water Balance and Vegetation Productivity across Site-to-Regional Scales

Evapotranspiration (ET) is a pivotal component in catchment-scale water balance and is essential for informed watershed management. Nevertheless, uncertainties in ET observation or modeling have been hindering effective water resources management. This study addresses this gap by establishing a robust, generalized linear relationship between ET and gross primary productivity (GPP) at the catchment scale. We test the linearity of the relationships between monthly GPP and ET data at 380 near-natural catchments across various climatic and landscape conditions in the contiguous U.S., yielding Pearson’s r ≥ 0.6 for 97% of the 380 catchments. We then develop a regionalization strategy to parameterize this GPP-ET relationship at the catchment scale by identifying and utilizing the linkages between the parameter values and extensively available hydroclimatic and landscape data. We demonstrate the efficacy of the proposed GPP-ET relationship and parameter regionalization strategy by their combined predictive capacity, where the predicted monthly GPP matches well with remote-sensing-based GPP product, achieving Kling-Gupta Efficient (KGE) values ≥ 0.5 for 92% of the catchments. In addition, we verify the relationship and its parameter regionalization at 35 AmeriFlux sites with KGE ≥ 0.5 for 25 sites, suggesting that the new relationship is transferable across the site, catchment, and regional scales. Furthermore, our findings are valuable for improving remote-sensing-based estimation of monthly ET and diagnosing coupled water–carbon simulations in land surface and Earth system models.

54 ENVIRONMENTAL SCIENCES↗

A parameterization technique for nonlinear spectral models

An objective parameterization technique is developed for general nonlinear hydrodynamical systems. The typical structure of hydrodynamical systems, regardless of their complexity, is one in which the rates of change of the dependent variables depend on homogeneous quadratic and linear forms, as well as on inhomogeneous forcing terms. As a prototype of the generic problem containing this typical structure, the parameterization technique is applied to various three component subsets of a five component nonlinear spectral model of forced, dissipative quasi-geostrophic flow in a channel. The results obtained lead to specification of the necessary data coverage requirements for applying the technique in general.

Gelaro, R.↗

A parameterization technique for nonlinear spectral models

An objective parameterization technique is developed for general nonlinear hydrodynamical systems. The typical structure of hydrodynamical systems, regardless of their complexity, is one in which the rates of change of the dependent variables depend on homogeneous quadratic and linear forms, as well as on inhomogeneous forcing terms. As a prototype of the generic problem containing this typical structure, the parameterization technique is applied to various three component subsets of a five component nonlinear spectral model of forced, dissipative quasi-geostrophic flow in a channel. The results obtained lead to specification of the necessary data coverage requirements for applying the technique in general.

Gelaro, R.↗

Two dimensional nonlinear analysis of sound transmission through a near-sonic throat flow

A two-dimensional nonlinear theory of sound transmission through a nonuniform duct carrying a near-sonic throat flow is described. The mean flow in the duct is treated according to a generalized quasi-one dimensional model, and the analysis of the unsteady perturbations is carried out using the Method of Matched Asymptotic Expansions. The linearized acoustic field in the subsonic regions of the duct is approximated by a Wave Envelope model which is matched asymptotically with a nonlinear inner solution valid in the near-sonic throat region. Numerical results are presented to illustrate the predictions of the theory. It is found that, in general, shock waves develop in the acoustic field during transmission through the near-sonic flow. Unlike previous one-dimensional theories, the current study shows that dispersion can play a major role in the propagation process.

Myers, M. K.↗

A diagnostic for denitrification in the winter polar stratospheres

The pairwise correlation of NO(y) and N2O data from the Southern and Northern Hemispheres is presented. Both data sets show a linear correlation region, defined as a reference state, and regions of denitrification where the correlation breaks down. Using two-dimensional photochemical model simulations of the atmosphere, a similar linear correlation is found between NO(y) and N2O, thereby establishing a theoretical framework for the reference state. This general approach, which can be extended to other pairs of molecules, should prove to be powerful in further comparisons of aircraft data with numerical models.

Fahey, D. W.↗

Estimating anisotropic diffusion of neutrons near the boundary of a pebble bed random system

Due to the arrangement of the pebbles in a Pebble Bed Reactor (PBR) core, if a neutron is located close to a boundary wall, its path length probability distribution function in directions of flight parallel to the wall is significantly different than in other directions. Hence, anisotropic diffusion of neutrons near the boundaries arises. We describe an analysis of neutron transport in a simplified 3-D pebble bed random system, in which we investigate the anisotropic diffusion of neutrons born near one of the system's boundary walls. While this simplified system does not model the actual physical process that takes place near the boundaries of a PBR core, the present work paves the road to a formulation that may enable more accurate diffusion simulations of such problems to be performed in the future. Monte Carlo codes have been developed for (i) deriving realizations of the 3-D random system, and (ii) performing 3-D neutron transport inside the heterogeneous model; numerical results are presented for three different choices of parameters. These numerical results are used to assess the accuracy of estimates for the mean-squared displacement of neutrons obtained with the diffusion approximations of the Atomic Mix Model and of the recently introduced [1] Non-Classical Theory with angular-dependent path length distribution. The Non-Classical Theory makes use of a Generalized Linear Boltzmann Equation in which the locations of the scattering centers in the system are correlated and the distance to collision is not exponentially distributed. We show that the results predicted using the Non-Classical Theory successfully model the anisotropic behavior of the neutrons in the random system, and more closely agree with experiment than the results predicted by the Atomic Mix Model. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Effective medium theory for a system of C{sub 60} molecules

An effective medium theory is developed to study the effective permittivity of a system of C{sub 60} molecules. We use a two-dimensional, spherical, two-fluid hydrodynamic model to describe the linear response of the π and σ electrons over the each C{sub 60} molecule. A general expression for the electromagnetic wave attenuation coefficient of the system is then deduced, and its functional dependence on the filling factor is presented. Furthermore, the dispersion characteristics of electromagnetic oscillations of the system are studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simple spatial scaling rules behind complex cities

Although most of wealth and innovation have been the result of human interaction and cooperation, we are not yet able to quantitatively predict the spatial distributions of three main elements of cities: population, roads, and socioeconomic interactions. By a simple model mainly based on spatial attraction and matching growth mechanisms, we reveal that the spatial scaling rules of these three elements are in a consistent framework, which allows us to use any single observation to infer the others. All numerical and theoretical results are consistent with empirical data from ten representative cities. In addition, our model can also provide a general explanation of the origins of the universal super- and sub-linear aggregate scaling laws and accurately predict kilometre-level socioeconomic activity. Our work opens a new avenue for uncovering the evolution of cities in terms of the interplay among urban elements, and it has a broad range of applications.

97 MATHEMATICS AND COMPUTING↗