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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Unveiling the hidden structure of complex stochastic biochemical networks

Complex Markov models are widely used and powerful predictive tools to analyze stochastic biochemical processes. However, when the network of states is unknown, it is necessary to extract information from the data to partially build the network and estimate the values of the rates. The short-time behavior of the first-passage time distributions between two states in linear chains has been shown recently to behave as a power of time with an exponent equal to the number of intermediate states. For a general Markov model we derive the complete Taylor expansion of the first-passage time distribution between two arbitrary states. By combining algebraic methods and graph theory approaches it is shown that the first term of the Taylor expansion is determined by the shortest path from the initial state to the final state. When this path is unique, we prove that the coefficient of the first term can be written in terms of the product of the transition rates along the path. It is argued that the application of our results to first-return times may be used to estimate the dependence of rates on external parameters in experimentally measured time distributions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

More about a successful vector-tensor theory of gravitation

The vector-tensor (VT) theory of gravitation revisited in this article was studied in previous papers, where it was proved that VT works and deserves attention. New observational data and numerical codes have motivated further development which is presented here. New research has been planed with the essential aim of proving that current cosmological observations, including Planck data, baryon acoustic oscillations (BAO), and so on, may be explained with VT, a theory which accounts for a kind of dark energy which has the same equation of state as vacuum. New versions of the codes CAMB and COSMOMC have been designed for applications to VT, and the resulting versions have been used to get the cosmological parameters of the VT model at suitable confidence levels. The parameters to be estimated are the same as in general relativity (GR), plus a new parameter D . For D = 0, VT linear cosmological perturbations reduces to those of GR, but the VT background may explain dark energy. The fits between observations and VT predictions lead to non vanishing | D | upper limits at the 1σ confidence level. The value D = 0 is admissible at this level, but this value is not that of the best fit in any case. Results strongly suggest that VT may explain current observations, at least, as well as GR; with the advantage that, as it is proved in this paper, VT has an additional parameter which facilitates adjustments to current observational data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Space–time reduced order model for large-scale linear dynamical systems with application to Boltzmann transport problems

A classical reduced order model for dynamical problems involves spatial reduction of the problem size. However, temporal reduction accompanied by the spatial reduction can further reduce the problem size without losing much accuracy, which results in a considerably more speed-up than the spatial reduction only. Recently, a novel space–time reduced order model for dynamical problems has been developed [17], where the space–time reduced order model shows an order of a hundred speed-up with a relative error of 10 –4 for small academic problems. However, in order for the method to be applicable to a large-scale problem, an efficient space–time reduced basis construction algorithm needs to be developed. Here we present the incremental space–time reduced basis construction algorithm. The incremental algorithm is fully parallel and scalable. Additionally, the block structure in the space–time reduced basis is exploited, which enables the avoidance of constructing the reduced space–time basis. These novel techniques are applied to a large-scale particle transport simulation with million and billion degrees of freedom. The numerical example shows that the algorithm is scalable and practical. Also, it achieves a tremendous speed-up, maintaining a good accuracy. Finally, error bounds for space-only and space–time reduced order models are derived.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

On effective temperature in network models of collective behavior

Collective behavior of self-propelled units is studied analytically within the Vectorial Network Model (VNM), a mean-field approximation of the well-known Vicsek model. We propose a dynamical systems framework to study the stochastic dynamics of the VNM in the presence of general additive noise. We establish that a single parameter, which is a linear function of the circular mean of the noise, controls the macroscopic phase of the system—ordered or disordered. By establishing a fluctuation–dissipation relation, we posit that this parameter can be regarded as an effective temperature of collective behavior. The exact critical temperature is obtained analytically for systems with small connectivity, equivalent to low-density ensembles of self-propelled units. Numerical simulations are conducted to demonstrate the applicability of this new notion of effective temperature to the Vicsek model. The identification of an effective temperature of collective behavior is an important step toward understanding order–disorder phase transitions, informing consistent coarse-graining techniques and explaining the physics underlying the emergence of collective phenomena.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Modeling and characterizing anisotropic inclusion orientation in heterogeneous material via directional cluster functions and stochastic microstructure reconstruction

We present a framework to model and characterize the microstructure of heterogeneous materials with anisotropic inclusions of secondary phases based on the directional correlation functions of the inclusions. Specifically, we have devised an efficient method to incorporate both directional two-point correlation functions S{sub 2} and directional two-point cluster functions C{sub 2} that contain non-trivial topological connectedness information into the simulated annealing microstructure reconstruction procedure. Our framework is applied to model an anisotropic aluminum alloy and the accuracy of the reconstructed structural models is assessed by quantitative comparison with the actual microstructure obtained via x-ray tomography. We show that incorporation of directional clustering information via C{sub 2} significantly improves the accuracy of the reconstruction. In addition, a set of analytical “basis” correlation functions are introduced to approximate the actual S{sub 2} and C{sub 2} of the material. With the proper choice of basis functions, the anisotropic microstructure can be represented by a handful of parameters including the effective linear sizes of the iron-rich and silicon-rich inclusions along three orthogonal directions. This provides a general and efficient means for heterogeneous material modeling that enables one to significantly reduce the data set required to characterize the anisotropic microstructure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Sensitivity of a global climate model to the critical Richardson number in the boundary layer parameterization

The critical bulk Richardson number (Ri cr ) is an important parameter in planetary boundary layer (PBL) parameterization schemes used in many climate models. This paper examines the sensitivity of a Global Climate Model, the Beijing Climate Center Atmospheric General Circulation Model, BCC_AGCM to Ri cr . The results show that the simulated global average of PBL height increases nearly linearly with Ri cr , with a change of about 114 m for a change of 0.5 in Ri cr . The surface sensible (latent) heat flux decreases (increases) as Ri cr increases. The influence of Ri cr on surface air temperature and specific humidity is not significant. The increasing Ri cr may affect the location of the Westerly Belt in the Southern Hemisphere. Further diagnosis reveals that changes in Ri cr affect stratiform and convective precipitations differently. Increasing Ri cr leads to an increase in the stratiform precipitation but a decrease in the convective precipitation. Significant changes of convective precipitation occur over the inter-tropical convergence zone, while changes of stratiform precipitation mostly appear over arid land such as North Africa and Middle East.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Contrasting interannual and multidecadal NAO variability

Decadal and longer timescale variability in the winter North Atlantic Oscillation (NAO) has considerable impact on regional climate, yet it remains unclear what fraction of this variability is potentially predictable. This study takes a new approach to this question by demonstrating clear physical differences between NAO variability on interannual-decadal (<30 year) and multidecadal (>30 year) timescales. It is shown that on the shorter timescale the NAO is dominated by variations in the latitude of the North Atlantic jet and storm track, whereas on the longer timescale it represents changes in their strength instead. NAO variability on the two timescales is associated with different dynamical behaviour in terms of eddy-mean flow interaction, Rossby wave breaking and blocking. The two timescales also exhibit different regional impacts on temperature and precipitation and different relationships to sea surface temperatures. These results are derived from linear regression analysis of the Twentieth Century and NCEP-NCAR reanalyses and of a high-resolution HiGEM General Circulation Model control simulation, with additional analysis of a long sea level pressure reconstruction. Evidence is presented for an influence of the ocean circulation on the longer timescale variability of the NAO, which is particularly clear in the model data. Finally, as well as providing new evidence of potential predictability, these findings are shown to have implications for the reconstruction and interpretation of long climate records

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Time-varying coefficient vector autoregressions model based on dynamic correlation with an application to crude oil and stock markets

This paper proposes a new time-varying coefficient vector autoregressions (VAR) model, in which the coefficient is a linear function of dynamic lagged correlation. The proposed model allows for flexibility in choices of dynamic correlation models (e.g. dynamic conditional correlation generalized autoregressive conditional heteroskedasticity (GARCH) models, Markov-switching GARCH models and multivariate stochastic volatility models), which indicates that it can describe many types of time-varying causal effects. Time-varying causal relations between West Texas Intermediate (WTI) crude oil and the US Standard and Poor’s 500 (S&P 500) stock markets are examined by the proposed model. The empirical results show that their causal relations evolve with time and display complex characters. Both positive and negative causal effects of the WTI on the S&P 500 in the subperiods have been found and confirmed by the traditional VAR models. Similar results have been obtained in the causal effects of S&P 500 on WTI. In addition, the proposed model outperforms the traditional VAR model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

An Extended MHD Study of the 16 October 2015 MMS Diffusion Region Crossing

The Magnetospheric Multiscale (MMS) mission has given us unprecedented access to high cadence particle and field data of magnetic reconnection at Earth's magnetopause. MMS first passed very near an X-line on 16 October 2015, the Burch event, and has since observed multiple X-line crossings. Subsequent 3-D particle-in-cell (PIC) modeling efforts of and comparison with the Burch event have revealed a host of novel physical insights concerning magnetic reconnection, turbulence-induced particle mixing, and secondary instabilities. In this study, we employ the Gkeyll simulation framework to study the Burch event with different classes of extended, multifluid magnetohydrodynamics (MHD), including models that incorporate important kinetic effects, such as the electron pressure tensor, with physics-based closure relations designed to capture linear Landau damping. Such fluid modeling approaches are able to capture different levels of kinetic physics in global simulations and are generally less costly than fully kinetic PIC. Herein we focus on the additional physics one can capture with increasing levels of fluid closure refinement via comparison with MMS data and existing PIC simulations. In particular, we find that the ten-moment model well captures the agyrotropic structure of the pressure tensor in the vicinity of the X-line and the magnitude of anisotropic electron heating observed in MMS and PIC simulations. However, the ten-moment model is found to have difficulty resolving the lower hybrid drift instability, which plays a fundamental role in heating and mixing electrons in the current layer.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Dynamically reconfigurable topological routing in nonlinear photonic systems

The propagation path of topologically protected states is bound to the interface between regions with different topology, and as such, the functionality of linear photonic devices leveraging these states is fixed during fabrication. Here, we propose a mechanism for dynamic control over a driven dissipative system’s local topology, yielding reconfigurable topological interfaces and thus tunable paths for protected routing. We illustrate our approach in non-resonantly pumped polariton lattices, where the nonlinear interaction between the polaritons and the exciton reservoir due to non-resonant pumping can yield picosecond-scale changes in the propagation paths of the chiral edge states. To analytically confirm the numerically observed topological dynamics, we generalize the spectral localizer framework to non-linear non-Hermitian Chern materials and apply this framework to a continuous model of the polariton system based on a driven-dissipative Gross-Pitaevskii equation. In doing so, we show that the local changes in the polariton lattice’s topology are captured by a local Chern marker. Looking forward, we anticipate such reconfigurable topological routing will enable the realization of novel classes of topological photonic devices.

Wong, Stephan [Sandia National Laboratories (SNL-N↗
DOE OSTI2012

Magnetised Kelvin-Helmholtz instability in the intermediate regime between subsonic and supersonic regimes

The understanding of the dynamics at play at the Earth's Magnetopause, the boundary separating the Earth's magnetosphere and the solar wind plasmas, is of primary importance for space plasma modeling. We focus our attention on the low latitude flank of the magnetosphere where the velocity shear between the magnetosheath and the magnetospheric plasmas is the energetic source of Kelvin-Helmholtz instability. On the shoulder of the resulting vortex chain, different secondary instabilities are at play depending on the local plasma parameters and compete with the vortex pairing process. Most important, secondary instabilities, among other magnetic reconnection, control the plasma mixing as well as the entry of solar wind plasma in the magnetosphere. We make use of a two-fluid model, including the Hall term and the electron mass in the generalized Ohm's law, to study the 2D non-linear evolution of the Kelvin-Helmholtz instability at the magnetosheath-magnetosphere interface, in the intermediate regime between subsonic and supersonic regimes. We study the saturation mechanisms, depending on the density jump across the shear layer and the magnetic field strength in the plane. In the presence of a weak in-plane magnetic field, the dynamics of the Kelvin-Helmholtz rolled-up vortices self-consistently generates thin current sheets where reconnection instability eventually enables fast reconnection to develop. Such a system enables to study guide field multiple-island collisionless magnetic reconnection as embedded in a large-scale dynamic system, unlike the classical static, ad hoc reconnection setups. In this regime, reconnection is shown to inhibit the vortex pairing process. This study provides a clear example of nonlinear, cross-scale, collisionless plasma dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Effect of rotation zero-crossing on single-fluid plasma response to three-dimensional magnetic perturbations

In order to understand the effect of rotation on the plasma's response to three-dimensional magnetic perturbations, we perform a systematic scan of the zero-crossing of the rotation profile in a DIII-D ITER-similar shape equilibrium using linear, time-independent modeling with the M3D-C1 extended magnetohydrodynamics code. We confirm that the local resonant magnetic field generally increases as the rotation decreases at a rational surface. Multiple peaks in the resonant field are observed near rational surfaces, however, and the maximum resonant field does not always correspond to zero rotation at the surface. Furthermore, we show that non-resonant current can be driven at zero- crossings not aligned with rational surfaces if there is sufficient shear in the rotation profile there, leading to an amplification of near-resonant Fourier harmonics of the perturbed magnetic field and a decrease in the far-off -resonant harmonics. The quasilinear electromagnetic torque induced by this non-resonant plasma response provides drive to flatten the rotation, possibly allowing for increased transport in the pedestal by the destabilization of turbulent modes. In addition, this torque acts to drive the rotation zero-crossing to dynamically stable points near rational surfaces, which would allow for increased resonant penetration. By one or both of these mechanisms, this torque may play an important role in bifurcations into ELM suppression. Finally, we discuss how these changes to the plasma response could be detected by tokamak diagnostics. In particular, we show that the changes to the resonant field discussed here have a significant impact on the external perturbed magnetic field, which should be observable by magnetic sensors on the high-field side of tokamaks, but not on the low-field side. In addition, TRIP3D-MAFOT simulations show that none of the changes to the plasma response described here substantially affects the divertor footprint structure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2010

Duality of parity doublets of helicity {+-}2 in D=2+1 dimensions

In D=2+1 dimensions there are two dual descriptions of parity singlets of helicity {+-}1, namely, the self-dual model of first order (in derivatives) and the Maxwell-Chern-Simons theory of second order. Correspondingly, for helicity {+-}2 there are four models S{sub SD{+-}}{sup (r)} describing parity singlets of helicities {+-}2. They are of first, second, third, and fourth order (r=1, 2, 3, 4), respectively. Here we show that the generalized soldering of the opposite helicity models S{sub SD+}{sup (4)} and S{sub SD-}{sup (4)} leads to the linearized form of the new massive gravity suggested by Bergshoeff, Hohm, and Townsend (BHT) similarly to the soldering of S{sub SD+}{sup (3)} and S{sub SD-}{sup (3)}. We argue why in both cases we have the same result. We also find out a triple master action which interpolates between the three dual models: linearized BHT theory, S{sub SD+}{sup (3)}+S{sub SD-}{sup (3)}, and S{sub SD+}{sup (4)}+S{sub SD-}{sup (4)}. By comparing gauge invariant correlation functions we deduce dual maps between those models. In particular, we learn how to decompose the field of the linearized BHT theory in helicity eigenstates of the dual models up to gauge transformations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Signal yields, energy resolution, and recombination fluctuations in liquid xenon

This study presents an analysis of monoenergetic electronic recoil peaks in the dark-matter-search and calibration data from the first underground science run of the Large Underground Xenon (LUX) detector. Liquid xenon charge and light yields for electronic recoil energies between 5.2 and 661.7 keV are measured, as well as the energy resolution for the LUX detector at those same energies. Additionally, there is an interpretation of existing measurements and descriptions of electron-ion recombination fluctuations in liquid xenon as limiting cases of a more general liquid xenon recombination fluctuation model. Measurements of the standard deviation of these fluctuations at monoenergetic electronic recoil peaks exhibit a linear dependence on the number of ions for energy deposits up to 661.7 keV, consistent with previous LUX measurements between 2 and 16 keV with 3 H. We highlight similarities in liquid xenon recombination for electronic and nuclear recoils with a comparison of recombination fluctuations measured with low-energy calibration data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1985

Calculation of steady and unsteady airfoil flow fields via the Navier-Stokes equations

A compressible time-dependent procedure for the two-dimensional ensemble averaged Navier-Stokes equations has been applied to the isolated airfoil problem in steady and unsteady flows. The procedure solves the governing equations via the linearized block implicit technique. Turbulence is modeled either via a mixing length or turbulence energy approach. The equations are solved in general non-orthogonal form with no-slip boundary conditions applied at the airfoil surface. Results are presented for airfoils at constant incidence, an airfoil in ramp motion and an airfoil oscillating through a dynamic stall loop. In general, steady converged solutions are obtained within 70 time steps over the range of Mach numbers considered. Comparisons with measured data show good agreement between computation and measurement.

Shamroth, S. J.↗
NASA NTRS2012

Optimal Initial Perturbations for Ensemble Prediction of the Madden-Julian Oscillation during Boreal Winter

An initialization strategy, tailored to the prediction of the Madden-Julian oscillation (MJO), is evaluated using the Goddard Earth Observing System Model, version 5 (GEOS-5), coupled general circulation model (CGCM). The approach is based on the empirical singular vectors (ESVs) of a reduced-space statistically determined linear approximation of the full nonlinear CGCM. The initial ESV, extracted using 10 years (1990-99) of boreal winter hindcast data, has zonal wind anomalies over the western Indian Ocean, while the final ESV (at a forecast lead time of 10 days) reflects a propagation of the zonal wind anomalies to the east over the Maritime Continent an evolution that is characteristic of the MJO. A new set of ensemble hindcasts are produced for the boreal winter season from 1990 to 1999 in which the leading ESV provides the initial perturbations. The results are compared with those from a set of control hindcasts generated using random perturbations. It is shown that the ESV-based predictions have a systematically higher bivariate correlation skill in predicting the MJO compared to those using the random perturbations. Furthermore, the improvement in the skill depends on the phase of the MJO. The ESV is particularly effective in increasing the forecast skill during those phases of the MJO in which the control has low skill (with correlations increasing by as much as 0.2 at 20 25-day lead times), as well as during those times in which the MJO is weak.

Ham, Yoo-Geun↗
DOE OSTI2017

Reactor operator-in-the-loop dynamic modelling - 122

The description of human operator dynamic characteristics in mathematical terms compatible with control engineering practice is an essential prerequisite to the analytical treatment of manual reactor control systems. Safe reactor operation requires effective operator control through interaction with plant dynamics, manipulators and displays. Traditional static analysis methods consider only specific situations; they fail to adequately explain the mutual interactions between the operator and the reactor plant characteristics. In this paper we investigate the theory for describing operator-reactor characteristics based on the methods of conventional control engineering techniques. The primary purpose of the experiments reported is the validation of the quasi-linear operator model. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2012

An Empirical Temperature Variance Source Model in Heated Jets

An acoustic analogy approach is implemented that models the sources of jet noise in heated jets. The equivalent sources of turbulent mixing noise are recognized as the differences between the fluctuating and Favre-averaged Reynolds stresses and enthalpy fluxes. While in a conventional acoustic analogy only Reynolds stress components are scrutinized for their noise generation properties, it is now accepted that a comprehensive source model should include the additional entropy source term. Following Goldstein s generalized acoustic analogy, the set of Euler equations are divided into two sets of equations that govern a non-radiating base flow plus its residual components. When the base flow is considered as a locally parallel mean flow, the residual equations may be rearranged to form an inhomogeneous third-order wave equation. A general solution is written subsequently using a Green s function method while all non-linear terms are treated as the equivalent sources of aerodynamic sound and are modeled accordingly. In a previous study, a specialized Reynolds-averaged Navier-Stokes (RANS) solver was implemented to compute the variance of thermal fluctuations that determine the enthalpy flux source strength. The main objective here is to present an empirical model capable of providing a reasonable estimate of the stagnation temperature variance in a jet. Such a model is parameterized as a function of the mean stagnation temperature gradient in the jet, and is evaluated using commonly available RANS solvers. The ensuing thermal source distribution is compared with measurements as well as computational result from a dedicated RANS solver that employs an enthalpy variance and dissipation rate model. Turbulent mixing noise predictions are presented for a wide range of jet temperature ratios from 1.0 to 3.20.

Khavaran, Abbas↗