Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Optimal Estimation of Clock Values and Trends from Finite Data

We show how to solve two problems of optimal linear estimation from a finite set of phase data. Clock noise is modeled as a stochastic process with stationary dth increments. The covariance properties of such a process are contained in the generalized autocovariance function (GACV). We set up two principles for optimal estimation: with the help of the GACV, these principles lead to a set of linear equations for the regression coefficients and some auxiliary parameters. The mean square errors of the estimators are easily calculated. The method can be used to check the results of other methods and to find good suboptimal estimators based on a small subset of the available data.

stationary increments↗
DOE OSTI2026

Simplified spin dependence in dark matter direct detection

The interactions of dark matter with Standard Model particles can be systematically studied in the language of effective field theories. We investigate dark matter interactions with Standard Model particles, including spin-dependent interactions, for direct detection experiments and demonstrate that, although the scattering rate generally depends on multiple types of material response functions, certain linear combinations of these material response functions vanish if the initial and final electronic states share the same Hamiltonian. We also find that several other response functions vanish in parity-symmetric materials, making these systems as simple as isotropic detectors in some respects. Finally, we present the scattering rate for an anisotropic, possibly chiral detector, for generic dark matter-electron spin interactions. These relations reduce the number of independent response functions needed, thereby simplifying the computational complexity for a broad class of dark matter models. Our results provide a complete and efficient toolkit for analyzing electron recoil signals in diverse detector materials.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A machine learning framework for accurate and robust analysis of radiation detector pulses

The microscopic properties of atomic nuclei are used to study various scientific questions. They are essential for understanding the fundamental forces of nature and the chemical evolution of the universe. Detecting decay radiation from radioactive nuclei makes it possible to probe these fundamental nuclear properties. Detector waveform traces may contain additional information about the radiation. Generally, advanced signal processing techniques are needed to extract this additional information, often involving fitting the waveform with model response functions using non-linear least-squares optimization with second-order gradient methods. While this is a powerful technique, it is also computationally expensive, leading to slow processing time, which scales with the volume of data. To address this problem, we have developed a machine learning (ML) approach that infers the characteristics of traces from a model detector response function. In particular, we are interested in classifying whether a single recorded trace consists of one or two pulse constituents and estimating the pulse parameters. Furthermore, our proposed ML method can precisely extract the pulses’ parameters, such as energy and timing information, and accurately classify the pulse multiplicity of a trace. Unlike non-learning-based approaches, our ML approach uses neural networks that are significantly faster at inference, as they do not require any optimization during this stage.

Curve fitting↗
DOE OSTI2020

Energy-Casimir, dynamically accessible, and Lagrangian stability of extended magnetohydrodynamic equilibria

The formal stability analysis of Eulerian extended magnetohydrodynamics (XMHD) equilibria is considered within the noncanonical Hamiltonian framework by means of the energy-Casimir variational principle and the dynamically accessible stability method. Specifically, we find explicit sufficient stability conditions for axisymmetric XMHD and Hall MHD (HMHD) equilibria with toroidal flow and for equilibria with arbitrary flow under constrained perturbations. The dynamically accessible, second-order variation of the Hamiltonian, which can potentially provide explicit stability criteria for generic equilibria, is also obtained. Moreover, we examine the Lagrangian stability of the general quasineutral two-fluid model written in terms of MHD-like variables, by finding the action and the Hamiltonian functionals of the linearized dynamics, working within a mixed Lagrangian-Eulerian framework. Upon neglecting electron mass, we derive a HMHD energy principle, and in addition, the perturbed induction equation arises from Hamilton’s equations of motion in view of a consistency condition for the relation between the perturbed magnetic potential and the canonical variables.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Transfer of contrast sensitivity in linear visual networks

Contrast sensitivity is a useful measure of the ability of an observer to distinguish contrast signals from noise. Although usually applied to human observers, contrast sensitivity can also be defined operationally for individual visual neurons. In a model linear neuron consisting of a filter and noise source, this operational measure is a function of filter gain, noise power spectrum, signal duration, and a performance criterion. This definition allows one to relate the sensitivities of linear neurons at different levels in the visual pathway. Mathematical formulas describing these relationships are derived, and the general model is applied to the specific problem of relating the sensitivities of parvocellular LGN neurons and cortical simple cells in the primate.

Watson, Andrew B.↗
NASA NTRS1993

Calculation of AGARD Wing 445.6 Flutter Using Navier-Stokes Aerodynamics

The flutter characteristics of the first AGARD standard aeroelastic configuration for dynamic response, Wing 445.6, are studied using an unsteady Navier-Stokes algorithm in order to investigate a previously noted discrepancy between Euler flutter characteristics and the experimental data. The algorithm, which is a three-dimensional, implicit, upwind Euler/Navier-Stokes code (CFL3D Version 2.1), was previously modified for the time-marching, aeroelastic analysis of wings using the unsteady Euler equations. These modifications include the incorporation of a deforming mesh algorithm and the addition of the structural equations of motion for their simultaneous time integration with the governing flow equations. In this paper, the aeroelastic method is extended and evaluated for applications that use the Navier- Stokes aerodynamics. The paper presents a brief description of the aeroelastic method and presents unsteady calculations which verify this method for Navier-Stokes calculations. A linear stability analysis and a time-marching aeroelastic analysis are used to determine the flutter characteristics of the isolated 45 deg. swept-back wing. Effects of fluid viscosity, structural damping, and number of modes in the structural model are investigated. For the linear stability analysis, the unsteady generalized aerodynamic forces of the wing are computed for a range of reduced frequencies using the pulse transfer-function approach. The flutter characteristics of the wing are determined using these unsteady generalized aerodynamic forces in a traditional V-g analysis. This stability analysis is used to determine the flutter characteristics of the wing at free-stream Mach numbers of 0.96 and 1.141 using the generalized aerodynamic forces generated by solving the Euler equations and the Navier-Stokes equations. Time-marching aeroelastic calculations are performed at a free-stream Mach number of 1.141 using the Euler and Navier-Stokes equations to compare with the linear V-g flutter analysis method. The V-g analysis, which is used in conjunction with the time-marching analysis, indicates that the fluid viscosity has a significant effect on the supersonic flutter boundary for this wing while the structural damping and number of modes in the structural model have a lesser effect.

Lee-Rausch, Elizabeth M.↗
DOE OSTI2020

The Generality of the GUGA MRCI Approach in COLUMBUS for Treating Complex Quantum Chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

hi_class background evolution, initial conditions and approximation schemes

Cosmological datasets have great potential to elucidate the nature of dark energy and test gravity on the largest scales available to observation. Theoretical predictions can be computed with hi_class (www.hiclass-code.net), an accurate, fast and flexible code for linear cosmology, incorporating a wide range of dark energy theories and modifications to general relativity. We introduce three new functionalities into hi_class: (1) Support for models based on covariant Lagrangians, including a constraint-preserving integration scheme for the background evolution and a series of worked-out examples: Galileon, nKGB, quintessence (monomial, tracker) and Brans-Dicke. (2) Consistent initial conditions for the scalar-field perturbations in the deep radiation era, identifying the conditions under which modified-gravity isocurvature perturbations may grow faster than adiabatic modes leading to a loss of predictivity. (3) An automated quasistatic approximation scheme allowing order-of-magnitude improvement in computing performance without sacrificing accuracy for wide classes of models. These enhancements bring the treatment of dark energy and modified gravity models to the level of detail comparable to software tools restricted to standard $\Lambda$CDM cosmologies. Finally, the hi_class code is publicly available (https://github.com/miguelzuma/hi_class_public), ready to explore current data and prepare for next-generation experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Towards enhancing and delaying disturbances in free shear flows

The family of shear flows comprising the jet, wake, and the mixing layer are subjected to perturbations in an inviscid incompressible fluid. By modeling the basic mean flows as parallel with piecewise linear variations for the velocities, complete and general solutions to the linearized equations of motion can be obtained in closed form as functions of all space variables and time when posed as an initial value problem. The results show that there is a continuous as well as the discrete spectrum that is more familiar in stability theory and therefore there can be both algebraic and exponential growth of disturbances in time. These bases make it feasible to consider control of such flows. To this end, the possibility of enhancing the disturbances in the mixing layer and delaying the onset in the jet and wake is investigated. It is found that growth of perturbations can be delayed to a considerable degree for the jet and the wake but, by comparison, cannot be enhanced in the mixing layer. By using moving coordinates, a method for demonstrating the predominant early and long time behavior of disturbances in these flows is given for continuous velocity profiles. It is shown that the early time transients are always algebraic whereas the asymptotic limit is that of an exponential normal mode. Numerical treatment of the new governing equations confirm the conclusions reached by use of the piecewise linear basic models. Although not pursued here, feedback mechanisms designed for control of the flow could be devised using the results of this work.

Criminale, W. O.↗
NASA NTRS1993

The analytical representation of viscoelastic material properties using optimization techniques

This report presents a technique to model viscoelastic material properties with a function of the form of the Prony series. Generally, the method employed to determine the function constants requires assuming values for the exponential constants of the function and then resolving the remaining constants through linear least-squares techniques. The technique presented here allows all the constants to be analytically determined through optimization techniques. This technique is employed in a computer program named PRONY and makes use of commercially available optimization tool developed by VMA Engineering, Inc. The PRONY program was utilized to compare the technique against previously determined models for solid rocket motor TP-H1148 propellant and V747-75 Viton fluoroelastomer. In both cases, the optimization technique generated functions that modeled the test data with at least an order of magnitude better correlation. This technique has demonstrated the capability to use small or large data sets and to use data sets that have uniformly or nonuniformly spaced data pairs. The reduction of experimental data to accurate mathematical models is a vital part of most scientific and engineering research. This technique of regression through optimization can be applied to other mathematical models that are difficult to fit to experimental data through traditional regression techniques.

Hill, S. A.↗
DOE OSTI2020

Perturbation theory for modeling galaxy bias: validation with simulations of the Dark Energy Survey

We describe perturbation theory (PT) models of galaxy bias for applications to photometric galaxy surveys. We model the galaxy-galaxy and galaxy-matter correlation functions in configuration space and validate against measurements from mock catalogs designed for the Dark Energy Survey (DES). We find that an effective PT model with five galaxy bias parameters provides a good description of the 3D correlation functions above scales of 4 Mpc/$h$ and $z < 1$. Our tests show that at the projected precision of the DES-Year 3 analysis, two of the non-linear bias parameters can be fixed to their co-evolution values, and a third (the $k^2$ term for higher derivative bias) set to zero. The agreement is typically at the 2 percent level over scales of interest, which is the statistical uncertainty of our simulation measurements. To achieve this level of agreement, our {\it fiducial} model requires using the full non-linear matter power spectrum (rather than the 1-loop PT one). We also measure the relationship between the non-linear and linear bias parameters and compare them to their expected co-evolution values. We use these tests to motivate the galaxy bias model and scale cuts for the cosmological analysis of the Dark Energy Survey; our conclusions are generally applicable to all photometric surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2003

Stochastic Mixing Model with Power Law Decay of Variance

Here we present a simple stochastic mixing model based on the law of large numbers (LLN). The reason why the LLN is involved in our formulation of the mixing problem is that the random conserved scalar c = c(t,x(t)) appears to behave as a sample mean. It converges to the mean value mu, while the variance sigma(sup 2)(sub c) (t) decays approximately as t(exp -1). Since the variance of the scalar decays faster than a sample mean (typically is greater than unity), we will introduce some non-linear modifications into the corresponding pdf-equation. The main idea is to develop a robust model which is independent from restrictive assumptions about the shape of the pdf. The remainder of this paper is organized as follows. In Section 2 we derive the integral equation from a stochastic difference equation describing the evolution of the pdf of a passive scalar in time. The stochastic difference equation introduces an exchange rate gamma(sub n) which we model in a first step as a deterministic function. In a second step, we generalize gamma(sub n) as a stochastic variable taking fluctuations in the inhomogeneous environment into account. In Section 3 we solve the non-linear integral equation numerically and analyze the influence of the different parameters on the decay rate. The paper finishes with a conclusion.

Fedotov, S.↗
DOE OSTI2023

Linear model decision trees as surrogates in optimization of engineering applications

Machine learning models are promising as surrogates in optimization when replacing difficult to solve equations or black-box type models. This work demonstrates the viability of linear model decision trees as piecewise-linear surrogates in decision-making problems. Linear model decision trees can be represented exactly in mixed-integer linear programming (MILP) and mixed-integer quadratic constrained programming (MIQCP) formulations. Furthermore, they can represent discontinuous functions, bringing advantages over neural networks in some cases. We present several formulations using transformations from Generalized Disjunctive Programming (GDP) formulations and modifications of MILP formulations for gradient boosted decision trees (GBDT). We then compare the computational performance of these different MILP and MIQCP representations in an optimization problem and illustrate their use on engineering applications. Importantly, we observe faster solution times for optimization problems with linear model decision tree surrogates when compared with GBDT surrogates using the Optimization and Machine Learning Toolkit (OMLT).

42 ENGINEERING↗
NASA NTRS1992

Transition operators in electromagnetic-wave diffraction theory - General theory

A formal theory is developed for the scattering of time-harmonic electromagnetic waves from impenetrable immobile obstacles with given linear, homogeneous, and generally nonlocal boundary conditions of Leontovich (impedance) type for the wave of the obstacle's surface. The theory is modeled on the complete Green's function and the transition (T) operator in time-independent formal scattering theory of nonrelativistic quantum mechanics. An expression for the differential scattering cross section for plane electromagnetic waves is derived in terms of certain matrix elements of the T operator for the obstacle.

Hahne, G. E.↗
NASA NTRS1989

A methodology for formulating a minimal uncertainty model for robust control system design and analysis

In the design and analysis of robust control systems for uncertain plants, the technique of formulating what is termed an M-delta model has become widely accepted and applied in the robust control literature. The M represents the transfer function matrix M(s) of the nominal system, and delta represents an uncertainty matrix acting on M(s). The uncertainty can arise from various sources, such as structured uncertainty from parameter variations or multiple unstructured uncertainties from unmodeled dynamics and other neglected phenomena. In general, delta is a block diagonal matrix, and for real parameter variations the diagonal elements are real. As stated in the literature, this structure can always be formed for any linear interconnection of inputs, outputs, transfer functions, parameter variations, and perturbations. However, very little of the literature addresses methods for obtaining this structure, and none of this literature addresses a general methodology for obtaining a minimal M-delta model for a wide class of uncertainty. Since have a delta matrix of minimum order would improve the efficiency of structured singular value (or multivariable stability margin) computations, a method of obtaining a minimal M-delta model would be useful. A generalized method of obtaining a minimal M-delta structure for systems with real parameter variations is given.

Belcastro, Christine M.↗
NASA NTRS1985

Magnetospheric impulse response for many levels of geomagnetic activity

The temporal relationship between the solar wind and magnetospheric activity has been studied using 34 intervals of high time resolution IMP 8 solar wind data and the corresponding AL auroral activity index. The median values of the AL index for each interval were utilized to rank the intervals according to geomagnetic activity level. The linear prediction filtering technique was then applied to model magnetospheric response as measured by the AL index to the solar wind input function VB(s). The linear prediction filtering routine produces a filter of time-lagged response coefficients which estimates the most general linear relationship between the chosen input and output parameters of the magnetospheric system. It is found that the filters are composed of two response pulses speaking at time lags of 20 and 60 min. The amplitude of the 60-min pulse is the larger for moderate activity levels, while the 20-min pulse is the larger for strong activity levels. A possible interpretation is that the 20-min pulse represents magnetospheric activity driven directly by solar wind coupling and that the 60-min pulse represents magnetospheric activity driven by the release of energy previously stored in the magnetotail. If this interpretation is correct, the linear filtering results suggest that both the driven and the unloading models of magnetospheric response are important facets of a more comprehensive response model.

Bargatze, L. F.↗
NASA NTRS2005

Constraining Oceanic Sources and Sinks of CO2 by Inverse Modeling: A First Step Toward Assimilation of In Situ and Remote Data

The primary goal of this collaborative project, lead by Nicholas Gruber at UCLA, is to better quantify the air-sea fluxes of pre-industrial and anthropogenic carbon dioxide CO2 for the global ocean with a newly developed inverse method. The method involves using ocean general circulation models (OGCMs) to characterize the subsurface distribution of tracers injected for different regions of the ocean surface. By linearly combining a series of these basis (or Greens) functions, one should be able to recreate any arbitrary ocean tracer field. In our case, we are using the global database of ocean dissolved inorganic carbon (DIC) from the WOCE/JGOFS global CO2 survey from the 1990s. The observed DIC field has been partitioned into pre-industrial and anthropogenic components using the empirical C* technique. We then solve for an optimal combination of the regional basis functions that best match the two sets observations in a least-squares sense with singular value decomposition. This in turn provides solutions for the regional air-sea CO2 fluxes (pre-industrial and anthropogenic) that are independent of more traditional estimates derived from surface water disequilibria in the partial pressure of CO2 and wind speed gas exchange parameterizations. We are apply the method to a suite of nine different OGCMs to test the robustness of the derived surface fluxes fields to differences in ocean circulation. Together with a collegue Keith Lindsay at NCAR, S. Doney's role in this project was to create the appropriate steady-state and temporally evolving basis functions from the NCAR CSM-1 ocean circulation model. Because of a slow start across the project in defining a common set of experimental protocols, we requested and were given a one year, no cost extension. Both the steady-state and temporally evolving basis functions were delivered to the UCLA group for further analysis. Additional information is included in the original extended abstract.

Doney, Scott↗
DOE OSTI2024

Linear Time-Invariant Models of a Large Cumulus Ensemble

Abstract Methods in system identification are used to obtain linear time-invariant state-space models that describe how horizontal averages of temperature and humidity of a large cumulus ensemble evolve with time under small forcing. The cumulus ensemble studied here is simulated with cloud-system-resolving models in radiative–convective equilibrium. The identified models extend steady-state linear response functions used in past studies and provide accurate descriptions of the transfer function, the noise model, and the behavior of cumulus convection when coupled with two-dimensional gravity waves. A novel procedure is developed to convert the state-space models into an interpretable form, which is used to elucidate and quantify memory in cumulus convection. The linear problem studied here serves as a useful reference point for more general efforts to obtain data-driven and interpretable parameterizations of cumulus convection.

Meteorology & Atmospheric Sciences↗