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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Operations for Learning with Graphical Models

This paper is a multidisciplinary review of empirical, statistical learning from a graphical model perspective. Well-known examples of graphical models include Bayesian net- works, directed graphs representing a Markov chain, and undirected networks representing a Markov field. These graphical models are extended to model data analysis and empirical learning using the notation of plates. Graphical operations for simplifying and manipulating a problem are provided including decomposition, differentiation, and the manipulation of probability models from the exponential family. These operations adapt existing techniques from statistics and automatic differentiation to graphs. Two standard algorithm schemes for learning are reviewed in a graphical framework: Gibbs sampling and the expectation maximization algorithm. Some algorithms are developed in this graphical framework including a generalized version of linear regression, techniques for feed-forward networks, and learning Gaussian and discrete Bayesian networks from data. The paper concludes by sketching some implications for data analysis and summarizing some popular algorithms that fall within the framework presented. The main original contributions here are the decomposition techniques and the demonstration that graphical models provide a framework for understanding and developing complex learning algorithms.

Buntine, Wray L.↗

DOE OSTI2010

FLUCTUATING ENERGY STORAGE AND NONLINEAR CASCADE IN AN INHOMOGENEOUS CORONAL LOOP

The dynamics and the energy balance of large-scale fluctuations in a coronal loop are studied. The loop is represented by a simplified structure where the curvature is neglected and the background magnetic field is uniform. In a previous paper, we studied a similar model where a uniform background density was assumed. The present paper represents a generalization of the previous one and it has the purpose of investigating possible modifications to the large-scale energy balance and dynamics due to a more realistic longitudinally nonuniform density. Large-scale fluctuations are dominated by coherent eigenmodes that nonlinearly couple to produce an energy cascade to smaller scales. Eigenmodes properties are calculated by a simplified linear dissipative model, deriving an expression for the input energy flux that is not substantially modified by the presence of the density inhomogeneity and is independent of dissipation. For typical values of the parameters, the derived input energy flux is comparable with that required to heat the active region corona. Nonlinear couplings are dominated by coherence effects due to the symmetry properties of eigenmodes; the consequences are that the system is in a weakly nonlinear regime that produces fluctuating energy storage in the loop, and that the kinetic and magnetic nonlinear energy fluxes are of the same order, despite the dominance of magnetic energy at large scales. From the energy balance, an expression for the velocity fluctuation is derived, which is valid in the more general case of a nonuniform background density; this estimate is in agreement both with measures of nonthermal velocities in the solar corona and with previous numerical results.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2021

Fitting strategies of accretion column models and application to the broadband spectrum of Cen X-3

Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, e.g., a cut off power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. In particular, Becker & Wolff(2007) presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this model that is both fast and accurate enough to be fitted to observed spectra is available in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during theχ2-minimization which simplifies the application of the model considerably. We demonstrate this approach as well as a study of possible parameter degeneracies with a comprehensive Markov-Chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT dataset of Cen X-3.The derived accretion-flow structure features a small column radius of∼63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies.

X-rays: binaries↗

NASA NTRS2021

Fitting Strategies of Accretion Column Models and Application to the Broadband Spectrum of Cen X-3

Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, e.g., a cut off power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. In particular, Becker & Wolff(2007) presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this model that is both fast and accurate enough to be fitted to observed spectra is available in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during theχ2-minimization which simplifies the application of the model considerably. We demonstrate this approach as well as a study of possible parameter degeneracies with a comprehensive Markov-Chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT data-set of Cen X-3.The derived accretion-flow structure features a small column radius of∼63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies

Philipp Thalhammer↗

NASA NTRS1994

An HP Adaptive Discontinuous Galerkin Method for Hyperbolic Conservation Laws

This dissertation addresses various issues for model classes of hyperbolic conservation laws. The basic approach developed in this work employs a new family of adaptive, hp-version, finite element methods based on a special discontinuous Galerkin formulation for hyperbolic problems. The discontinuous Galerkin formulation admits high-order local approximations on domains of quite general geometry, while providing a natural framework for finite element approximations and for theoretical developments. The use of hp-versions of the finite element method makes possible exponentially convergent schemes with very high accuracies in certain cases; the use of adaptive hp-schemes allows h-refinement in regions of low regularity and p-enrichment to deliver high accuracy, while keeping problem sizes manageable and dramatically smaller than many conventional approaches. The use of discontinuous Galerkin methods is uncommon in applications, but the methods rest on a reasonable mathematical basis for low-order cases and has local approximation features that can be exploited to produce very efficient schemes, especially in a parallel, multiprocessor environment. The place of this work is to first and primarily focus on a model class of linear hyperbolic conservation laws for which concrete mathematical results, methodologies, error estimates, convergence criteria, and parallel adaptive strategies can be developed, and to then briefly explore some extensions to more general cases. Next, we provide preliminaries to the study and a review of some aspects of the theory of hyperbolic conservation laws. We also provide a review of relevant literature on this subject and on the numerical analysis of these types of problems.

Bey, Kim S.↗

DOE OSTI2019

HALO: A full-orbit model of nonlinear interaction of fast particles with eigenmodes

HALO (HAgis LOcust) solves the initial value Vlasov-Maxwell problem perturbatively for application to certain nonlinear wave-particle problems in tokamak plasmas. It uses the same basic approach as the HAGIS code for wave evolution but is built on the LOCUST-GPU full-orbit code for the solution of the Hamiltonian fast particle motion in cylindrical coordinates. The wave amplitude and particle evolution include all finite Larmour radius effects. We describe and benchmark the currently implemented Alfvén eigenmode workflow, demonstrating correct particle motion, linear and nonlinear power transfer. Here, the formulation and numerical scheme are sufficiently general as to allow easy future implementation of different kinds of eigenmodes, such as modes close to the ion-cyclotron frequency. The code can model multiple eigenmodes and multiple fast ion species simultaneously, and supports the general form of the equilibrium distribution in constants of motion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2018

Orbifolds and Exact Solutions of Strongly-Coupled Matrix Models

Here we find an exact solution to strongly-coupled matrix models with a single-trace monomial potential. Our solution yields closed form expressions for the partition function as well as averages of Schur functions. The results are fully factorized into a product of terms linear in the rank of the matrix and the parameters of the model. We extend our formulas to include both logarthmic and finite-difference deformations, thereby generalizing the celebrated Selberg and Kadell integrals. We conjecture a formula for correlators of two Schur functions in these models, and explain how our results follow from a general orbifold-like procedure that can be applied to any one-matrix model with a single-trace potential.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

An Open‐Source, Physics‐Based, Tropical Cyclone Downscaling Model With Intensity‐Dependent Steering

Abstract An open‐source, physics‐based tropical cyclone (TC) downscaling model is developed, in order to generate a large climatology of TCs. The model is composed of three primary components: (a) a random seeding process that determines genesis, (b) an intensity‐dependent beta‐advection model that determines the track, and (c) a non‐linear differential equation set that determines the intensification rate. The model is entirely forced by the large‐scale environment. Downscaling ERA5 reanalysis data shows that the model is generally able to reproduce observed TC climatology, such as the global seasonal cycle, genesis locations, track density, and lifetime maximum intensity distributions. Inter‐annual variability in TC count and power‐dissipation is also well captured, on both basin‐wide and global scales. Regional TC hazard estimated by this model is also analyzed using return period maps and curves. In particular, the model is able to reasonably capture the observed return period curves of landfall intensity in various sub‐basins around the globe. The incorporation of an intensity‐dependent steering flow is shown to lead to regionally dependent changes in power dissipation and return periods. Advantages and disadvantages of this model, compared to other downscaling models, are also discussed.

Meteorology & Atmospheric Sciences↗

DOE OSTI2013

NON-ZEEMAN CIRCULAR POLARIZATION OF MOLECULAR ROTATIONAL SPECTRAL LINES

We present measurements of circular polarization from rotational spectral lines of molecular species in Orion KL, most notably {sup 12}CO (J = 2 {yields} 1), obtained at the Caltech Submillimeter Observatory with the Four-Stokes-Parameter Spectral Line Polarimeter. We find levels of polarization of up to 1%-2% in general; for {sup 12}CO (J = 2 {yields} 1) this level is comparable to that of linear polarization also measured for that line. We present a physical model based on resonant scattering in an attempt to explain our observations. We discuss how slight differences in scattering amplitudes for radiation polarized parallel and perpendicular to the ambient magnetic field, responsible for the alignment of the scattering molecules, can lead to the observed circular polarization. We also show that the effect is proportional to the square of the magnitude of the plane of the sky component of the magnetic field and therefore opens up the possibility of measuring this parameter from circular polarization measurements of Zeeman insensitive molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2010

Parameterizing Bayesian network Representations of Social-Behavioral Models by Expert Elicitation

Bayesian networks provide a general framework with which to model many natural phenomena. The mathematical nature of Bayesian networks enables a plethora of model validation and calibration techniques: e.g parameter estimation, goodness of fit tests, and diagnostic checking of the model assumptions. However, they are not free of shortcomings. Parameter estimation from relevant extant data is a common approach to calibrating the model parameters. In practice it is not uncommon to find oneself lacking adequate data to reliably estimate all model parameters. In this paper we present the early development of a novel application of conjoint analysis as a method for eliciting and modeling expert opinions and using the results in a methodology for calibrating the parameters of a Bayesian network.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2013

Uniform and Half‐Space Models for Magnetic Dipole Sources in a Conducting Medium: I. Theory and Results

Underground nuclear events result in electromagnetic signals near the event location. Of particular interest, because its signal can be potentially enhanced by the use of large cavities to decouple the seismic signal, is the generation of an EMP signal within the nuclear detonation (NUDET) cavity, from the formation of a bubble in the local magnetic field. This mechanism is equivalent to a buried time varying magnetic dipole. In order to model such signals Los Alamos National Laboratory is developing a detailed model for the propagation of RF signals from buried sources to nearby sensors, with a focus on sources that are magnetic dipole in character. In order to validate this code, models for simple geometries involving magnetic dipole sources have been implemented. These models represent the fields from vertical and horizontal dipoles in cylindrical coordinates, either in a uniform linear isotropic electrical medium, or in a space consisting of two half spaces interfacing at the plane z=0, with the upper half space having the electrical properties of free space and the lower half space having arbitrary linear conductivity and electrical permittivity. While primarily intended for the validation of the more detailed model, the shorter turnaround time of these models facilitates the study of the general properties of geomagnetic bubble signals. This report summarizes the theory and example results for those models. The use of the model for model validation is illustrated with the near and surface fields for dipoles with Gaussian time dependence. Also examined are the surface fields expected for the geomagnetic bubble predicted by the numerical model of Kovalenko et al for a 1 kt nuclear detonation. The model of Kovalenko et al was also rescaled to allow comparison with Zablocki’s electric field observations for Bilby, Event III, and Event V. While the rescaled Kovalenko et al model had approximately the right time scale for the duration of the main peak, and a slower fall than rise also in qualitative agreement with the observations, the predicted signals had very different polarizations, were too low by more than an order of magnitude, and had less temporal structure. A simple comparison with the vertical B field data of Sweeney, where the detailed event geometries and yields are classified, but the sensor range and bounds on the yields are given, suggest that the scaled Kovalenko et al model gives too small a peak field compared to observations for that data set as well, except perhaps for the Borate event. A companion report will give a detailed discussion of the implementation and code.

58 GEOSCIENCES↗

NASA NTRS1979

A two-dimensional thermal analysis of a new high-performance tubular solar collector

The first of two articles are presented which describe and analyze the thermal performance of a vacuum tube solar collector. The assumptions and mathematical modeling are presented. The problem is reduced to the formulation of two simultaneous linear differential equations characterizing the collector thermal behavior. After applying the boundary conditions, a general solution is obtained which is found similar to the general Hottel, Whillier and Bliss form but with a complex flow factor.

Lansing, F. L.↗

NASA NTRS1981

An optimization model for energy generation and distribution in a dynamic facility

An analytical model is described using linear programming for the optimum generation and distribution of energy demands among competing energy resources and different economic criteria. The model, which will be used as a general engineering tool in the analysis of the Deep Space Network ground facility, considers several essential decisions for better design and operation. The decisions sought for the particular energy application include: the optimum time to build an assembly of elements, inclusion of a storage medium of some type, and the size or capacity of the elements that will minimize the total life-cycle cost over a given number of years. The model, which is structured in multiple time divisions, employ the decomposition principle for large-size matrices, the branch-and-bound method in mixed-integer programming, and the revised simplex technique for efficient and economic computer use.

Lansing, F. L.↗

NASA NTRS1983

A compressible model for the simulation of moist mountain waves

A two-dimensional, nonlinear, nonhydrostatic model is described which allows the calculation of moist airflow in mountainous terrain. The model is compressible, uses a terrain-following coordinate system, and employs lateral and upper boundary conditions which minimize wave reflections. The model's accuracy and sensitivity are examined. These tests suggest that in numerical simulations of vertically propagating, highly nonlinear mountain waves, a wave absorbing layer does not accurately mimic the effects of wave breakdown and dissipation at high levels in the atmosphere. In order to obtain a correct simulation, the region in which the waves are physically absorbed must generally be included in the computational domain (a nonreflective upper boundary condition should be used as well). The utility of the model is demonstrated in two examples (linear waves in a uniform atmosphere and the Jan. 11, 1972 Boulder windstorm) which illustrate how the presence of moisture can influence propagating waves. In both cases, the addition of moisture to the upstream flow greatly reduces the wave response.

Durran, D. R.↗

NASA NTRS2014

In-Flight Aeroelastic Stability of the Thermal Protection System on the NASA HIAD, Part I: Linear Theory

Conical shell theory and piston theory aerodynamics are used to study the aeroelastic stability of the thermal protection system (TPS) on the NASA Hypersonic Inflatable Aerodynamic Decelerator (HIAD). Structural models of the TPS consist of single or multiple orthotropic conical shell systems resting on several circumferential linear elastic supports. The shells in each model may have pinned (simply-supported) or elastically-supported edges. The Lagrangian is formulated in terms of the generalized coordinates for all displacements and the Rayleigh-Ritz method is used to derive the equations of motion. The natural modes of vibration and aeroelastic stability boundaries are found by calculating the eigenvalues and eigenvectors of a large coefficient matrix. When the in-flight configuration of the TPS is approximated as a single shell without elastic supports, asymmetric flutter in many circumferential waves is observed. When the elastic supports are included, the shell flutters symmetrically in zero circumferential waves. Structural damping is found to be important in this case. Aeroelastic models that consider the individual TPS layers as separate shells tend to flutter asymmetrically at high dynamic pressures relative to the single shell models. Several parameter studies also examine the effects of tension, orthotropicity, and elastic support stiffness.

Goldman, Benjamin D.↗

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

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