Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 325 records · Page 18

Small Signal Modelling and Control of the Hydrogen Mixer for Facility E1

We have undertaken the theoretical modelling of an existing liquid hydrogen (LH2) and gas hydrogen (GH2) mixer subsystem of the E1 Ground Test Facility at NASA John C. Stennis Space Center. The E1 test facility carries out comprehensive ground-based testing and certification of various liquid rocket engines and their components. The mixer described in this work is responsible for combining high pressure LH2 and GH2 to produce a hydrogen flow that meets certain thermodynamic properties before it is fed into a test article. The desired properties are maintained by precise control of the mixture of LH2 and GH2 flows. The mixer is modelled as a general multi-flow lumped volume for single constituent fluids using density and internal energy as states. The set of nonlinear differential equations is linearized about an equilibrium point and the resulting two-state, 3-input linear model is analyzed as a possible candidate for control design.

Barbieri, Enrique↗

Small Signal Modeling and Control of the Hydrogen Mixer for Facility E1

We have undertaken the theoretical modelling of an existing liquid hydrogen (LH2) and gas hydrogen (GH2) mixer subsystem of the E1 Ground Test Facility at NASA John C. Stennis Space Center. The E1 test facility carries out comprehensive ground-based testing and certification of various liquid rocket engines and their components. The mixer described in this work is responsible for combining high pressure before it is fed into a test article. The desired properties are maintained by precise control of the mixture of LH2 and GH2 flows. The mixer is modelled as a general multi-flow lumped volume for single constituent fluids using density and internal energy as states. The set of nonlinear differential equations is linearized about an equilibrium point and the resulting two-state, 3-input linear model is analyzed as a possible candidate for control design.

Barbieri, Enrique↗

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↗

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.↗

Generative Patch Priors for Practical Compressive Image Recovery

In this work, we propose the generative patch prior (GPP) that defines a generative prior for compressive image recovery, based on patch-manifold models. Unlike learned, image-level priors that are restricted to the range space of a pre-trained generator, GPP can recover a wide variety of natural images using a pre-trained patch generator. Additionally, GPP retains the benefits of generative priors like high reconstruction quality at extremely low sensing rates, while also being much more generally applicable. We show that GPP outperforms several unsupervised and supervised techniques on three different sensing models – linear compressive sensing with known, and unknown calibration settings, and the non-linear phase retrieval problem. Finally, we propose an alternating optimization strategy using GPP for joint calibration-and-reconstruction which performs favorably against several baselines on a real world, uncalibrated compressive sensing dataset. The code and models for GPP are available on github.

compressed sensing↗

Data for "Soil carbon dynamics during drying vs. rewetting: Importance of antecedent moisture conditions"

This dataset contains data used for the paper "Soil carbon dynamics during drying vs. rewetting: importance of antecedent moisture conditions". The Related References field will be updated with a full citation when available.Soil moisture influences soil carbon dynamics, including microbial growth and respiration. The response of such ‘soil respiration’ to moisture changes is generally assumed to be linear and reversible, i.e. to depend only on the current moisture state. Current models thus do not account for antecedent soil moisture conditions when determining soil respiration or the available substrate pool. We conducted a laboratory incubation to determine how the antecedent conditions of drought and flood influenced soil organic matter (SOM) chemistry, bioavailability, and respiration. We sampled soils from an upland coastal forest, Beaver Creek, WA USA, and subjected them to drying and rewetting treatments. For the drying treatment, field moist soils were saturated and then dried to 75, 50, 35, and 5 % saturation. In the rewetting treatment, field moist soils were air-dried and then rewet to 35, 50, 75, and 100 % saturation. We measured respiration and water extractable organic carbon (WEOC) concentrations and used 1H-NMR and FT-ICR-MS to characterize the WEOC pool across the treatments. The drying vs. wetting treatment strongly influenced SOM bioavailability, as rewet soils (with antecedent drought) had greater WEOC concentrations and respiration fluxes compared to the drying soils (with antecedent flood). In addition, air-dry soils had the highest WEOC concentrations, and the NMR-resolved peaks showed a strong contribution of protein groups in these soils. Both NMR and FT-ICR-MS analyses indicated increased contribution of complex aromatic groups/molecules in the rewet soils, compared to the drying soils. We suggest that drying introduced organic matter into the WEOC pool via desorption of aromatic molecules and/or by microbial cell lysis, and this stimulated microbial mineralization rates. Our work indicates that even short-term shifts in antecedent moisture conditions can strongly influence soil C dynamics at the core scale. The predictive uncertainties in current soil models may be reduced by a more accurate representation of soil water and C persistence that includes a mechanistic and quantitative understanding of the impact of antecedent moisture conditions.This dataset contains a compressed (.zip) archive of the data and R scripts used for this manuscript. The dataset includes files in .csv and .txt format, which can be accessed and processed using MS Excel or R. NMR data are provided as raw output data (accessed in Bruker TopSpin or MestreNova) as well as the MestreNova-processed files. This archive can also be accessed on GitHub at https://github.com/kaizadp/hysteresis_and_soil_carbon (DOI: 10.5281/zenodo.4432885).

1H-NMR↗