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 649 records · Page 36

BEYONDPLANCK IV. On end-to-end simulations in CMB analysis — Bayesian versus frequentist statistics

End-to-end simulations play a key role in the analysis of any high-sensitivity cosmic microwave background (CMB) experiment, providing high-fidelity systematic error propagation capabilities that are unmatched by any other means. In this paper, we address an important issue regarding such simulations, namely, how to define the inputs in terms of sky model and instrument parameters. These may either be taken as a constrained realization derived from the data or as a random realization independent from the data. We refer to these as posterior and prior simulations, respectively. We show that the two options lead to significantly different correlation structures, as prior simulations (contrary to posterior simulations) effectively include cosmic variance, but they exclude realization-specific correlations from non-linear degeneracies. Consequently, they quantify fundamentally different types of uncertainties. We argue that as a result, they also have different and complementary scientific uses, even if this dichotomy is not absolute. In particular, posterior simulations are in general more convenient for parameter estimation studies, while prior simulations are generally more convenient for model testing. Before BEYONDPLANCK, most pipelines used a mix of constrained and random inputs and applied the same hybrid simulations for all applications, even though the statistical justification for this is not always evident. BEYONDPLANCK represents the first end-to-end CMB simulation framework that is able to generate both types of simulations and these new capabilities have brought this topic to the forefront. The BEYONDPLANCK posterior simulations and their uses are described extensively in a suite of companion papers. In this work, we consider one important applications of the corresponding prior simulations, namely, code validation. Specifically, we generated a set of one-year LFI 30 GHz prior simulations with known inputs and we used these to validate the core low-level BEYONDPLANCK algorithms dealing with gain estimation, correlated noise estimation, and mapmaking.

79 ASTRONOMY AND ASTROPHYSICS↗

A novel CMB component separation method: hierarchical generalized morphological component analysis

ABSTRACT We present a novel technique for cosmic microwave background (CMB) foreground subtraction based on the framework of blind source separation. Inspired by previous work incorporating local variation to generalized morphological component analysis (GMCA), we introduce hierarchical GMCA (HGMCA), a Bayesian hierarchical graphical model for source separation. We test our method on Nside = 256 simulated sky maps that include dust, synchrotron, free–free, and anomalous microwave emission, and show that HGMCA reduces foreground contamination by $25{{\ \rm per\ cent}}$ over GMCA in both the regions included and excluded by the Planck UT78 mask, decreases the error in the measurement of the CMB temperature power spectrum to the 0.02–0.03 per cent level at ℓ > 200 (and $\lt 0.26{{\ \rm per\ cent}}$ for all ℓ), and reduces correlation to all the foregrounds. We find equivalent or improved performance when compared to state-of-the-art internal linear combination type algorithms on these simulations, suggesting that HGMCA may be a competitive alternative to foreground separation techniques previously applied to observed CMB data. Additionally, we show that our performance does not suffer when we perturb model parameters or alter the CMB realization, which suggests that our algorithm generalizes well beyond our simplified simulations. Our results open a new avenue for constructing CMB maps through Bayesian hierarchical analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

Monolithic ceramic analysis using the SCARE program

The Structural Ceramics Analysis and Reliability Evaluation (SCARE) computer program calculates the fast fracture reliability of monolithic ceramic components. The code is a post-processor to the MSC/NASTRAN general purpose finite element program. The SCARE program automatically accepts the MSC/NASTRAN output necessary to compute reliability. This includes element stresses, temperatures, volumes, and areas. The SCARE program computes two-parameter Weibull strength distributions from input fracture data for both volume and surface flaws. The distributions can then be used to calculate the reliability of geometrically complex components subjected to multiaxial stress states. Several fracture criteria and flaw types are available for selection by the user, including out-of-plane crack extension theories. The theoretical basis for the reliability calculations was proposed by Batdorf. These models combine linear elastic fracture mechanics (LEFM) with Weibull statistics to provide a mechanistic failure criterion. Other fracture theories included in SCARE are the normal stress averaging technique and the principle of independent action. The objective of this presentation is to summarize these theories, including their limitations and advantages, and to provide a general description of the SCARE program, along with example problems.

Manderscheid, Jane M.↗

Sparse measurement medical CT reconstruction using multi-fused block matching denoising priors

A major challenge for medical X-ray CT imaging is reducing the number of X-ray projections to lower radiation dosage and reduce scan times without compromising image quality. However these under-determined inverse imaging problems rely on the formulation of an expressive prior model to constrain the solution space while remaining computationally tractable. Traditional analytical reconstruction methods like Filtered Back Projection (FBP) often fail with sparse measurements, producing artifacts due to their reliance on the Shannon-Nyquist Sampling Theorem. Consensus Equilibrium, which is a generalization of Plug and Play, is a recent advancement in Model-Based Iterative Reconstruction (MBIR), has facilitated the use of multiple denoisers are prior models in an optimization free framework to capture complex, non-linear prior information. However, 3D prior modelling in a Plug and Play approach for volumetric image reconstruction requires long processing time due to high computing requirement. Instead of directly using a 3D prior, this work proposes a BM3D Multi Slice Fusion (BM3D-MSF) prior that uses multiple 2D image denoisers fused to act as a fully 3D prior model in Plug and Play reconstruction approach. Our approach does not require training and are thus able to circumvent ethical issues related with patient training data and are readily deployable in varying noise and measurement sparsity levels. In addition, reconstruction with the BM3D-MSF prior achieves similar reconstruction image quality as fully 3D image priors, but with significantly reduced computational complexity. We test our method on clinical CT data and demonstrate that our approach improves reconstructed image quality.

Hossain, Maliha [ORNL]↗

A manifold learning perspective on surrogate modeling of nitrate concentration in the Kansas River

Abstract A non-linear surrogate model of nitrate concentration in the Kansas River (USA) is described. The model is an (almost) Piece-wise Linear response surface that provides a mean field approximation to the dynamics of the measured data for nitrate plus nitrite (target product) correlations to turbidity and chlorophyll-a concentrations (input variables). The method extends the United States Geological Survey’s linear procedures for surrogate data modeling allowing for better approximations for river systems exhibiting algal blooms due to nutrient-rich source waters. The model and visualization procedures illustrated in the Kansas River example should be generally applicable to many medium-size rivers in agricultural regions.

Tufillaro, Nicholas (ORCID:0009000628968832)↗

Quantum Non-Equilibrium Dynamics Prediction of Electronic Transport Coefficients in Non-Linear Regimes

The goal of this work was to develop unique LLNL capabilities for the computational prediction of electrical, optical and thermal transport properties from quantum simulations. Using a recently developed theory of microscopic quantum transport, we were able to model compressed metallic hydrogen and it’s non-linear and out-of-equilibrium contributions in electrical and thermal transport. This theory was implemented in the open-source, real-space time-dependent density functional theory code Octopus. This model is general and widely applicable to high energy density systems that are important to LLNL and NNSA mission areas.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A statistical study of Pc 3 magnetic pulsations at synchronous orbit, ATS 6

Results are reported for a study designed to determine statistically the nature and variation of the polarization parameters of Pc3 magnetic pulsations at synchronous orbit. Data obtained with a triaxial flux gate magnetometer on the geosynchronous satellite ATS 6 are subjected to routine power spectral and polarization analyses, and the events examined are identified as peaks in the spectra. Coherence analysis of spectral matrices for these events indicates that Pc3 are very common at synchronous orbit and are primarily a local morning phenomenon with an occurrence peak around 1000 LT; that the pulsations occur at all frequencies in the Pc3 band, with the most common frequency at about 0.35 Hz; that the Pc3 at synchronous orbit are most often linear and generally transverse, but some have significant compressional components; and that based on azimuth angle, there are two distinct classes of Pc3 (azimuthal and radial) with different characteristics. These statistical results are compared with model predictions and ground-based observations.

Arthur, C. W.↗

Alternatives for jet engine control

Alternatives to linear quadratic regulator theory in the linear case are examined along with nonlinear modelling and optimization approaches for global control. Context for the studies has been set by the DYNGEN digital simulator and by models generated for various phases of the F100 Multivariable Control Synthesis Program. With respect to the linear alternatives, the multivariable frequency domain is stressed. Progress is reported in both the direct algebraic approach to exact model matching, by means of stimulating work on the basic computational issues, and in the indirect generalized Nyquist approach. With respect to nonlinear modelling and optimization, the emphasis is twofold: the development of analytical nonlinear models of the jet engine and the use of these models in conjunction with techniques of mathematical programming in order to study global control over nonincremental portions of the flight envelope. The possibility of using tensor methods is explored.

Sain, M. K.↗

Measured and predicted structural behavior of the HiMAT tailored composite wing

A series of load tests was conducted on the HiMAT tailored composite wing. Coupon tests were also run on a series of unbalanced laminates, including the ply configuration of the wing, the purpose of which was to compare the measured and predicted behavior of unbalanced laminates, including - in the case of the wing - a comparison between the behavior of the full scale structure and coupon tests. Both linear and nonlinear finite element (NASTRAN) analyses were carried out on the wing. Both linear and nonlinear point-stress analyses were performed on the coupons. All test articles were instrumented with strain gages, and wing deflections measured. The leading and trailing edges were found to have no effect on the response of the wing to applied loads. A decrease in the stiffness of the wing box was evident over the 27-test program. The measured load-strain behavior of the wing was found to be linear, in contrast to coupon tests of the same laminate, which were nonlinear. A linear NASTRAN analysis of the wing generally correlated more favorably with measurements than did a nonlinear analysis. An examination of the predicted deflections in the wing root region revealed an anomalous behavior of the structural model that cannot be explained. Both hysteresis and creep appear to be less significant in the wing tests than in the corresponding laminate coupon tests.

Nelson, Lawrence H.↗

On neglect of nonlinear momentum terms in solar nebula accretion disk models

The standard model for accretion disks (e.g., Pringle, 1981; Lin and Papaloizou, 1985) parameterizes the radial eddy angular momentum flux which is quadratically nonlinear in velocity, solely by a viscous diffusion process which is linear in velocity. This paper shows that the eddy angular momentum flux, in addition to a gradient dissipative (viscous) flux, involves a countergradient accelerating flux. Therefore, the momentum parameterization in the standard model is in general unjustified. The formalism of the parameterization can be justified only by recognizing that the eddy viscosity used in the standard model is in fact a pseudoviscosity which includes the effects of countergradient fluxes, and which therefore may have a magnitude less than the eddy diffusivity K for mixing trace constituents in the disk.

Prinn, Ronald G.↗

A study of control approaches for large-angle beam slewing with an adaptive truss

This paper presents an analytical study of large-angle beam slewing with a double-octahedral adaptive truss and documents the performance of different control laws for simultaneous vibration control and beam slewing. A general nonlinear model of an adaptive truss and flexible beam is derived which includes a degree-of-freedom previously considered negligible for a double-octahedral adaptive truss actuated at its mid-plane. Linear and nonlinear, gain-scheduled, control laws are developed and applied to the nonlinear model for various large-angle slewing maneuvers. The nonlinear, gain-scheduled, control law is found to perform better than the linear control laws.

Warrington, Thomas J.↗

The embedded young stars in the Taurus-Auriga molecular cloud. II - Models for scattered light images

We describe NIR imaging observations of embedded young stars in the Taurus-Auriga molecular cloud. We find a large range in J-K and H-K colors for these class I sources. The bluest objects have colors similar to the reddest T Tauri stars in the cloud; redder objects lie slightly above the reddening line for standard ISM dust and have apparent K extinctions of up to 5 mag. Most of these sources also show extended NIR emission on scales of 10-20 arcsec which corresponds to linear sizes of 1500-3000 AU. The NIR colors and nebular morphologies for this sample and the magnitude of linear polarization in several sources suggest scattered light produces most of the NIR emission in these objects. We present modeling results that suggest mass infall rates that agree with predictions for cold clouds and are generally consistent with rates estimated from radiative equilibrium models. For reasonable dust grain parameters, the range of colors and extinctions require flattened density distributions with polar cavities evacuated by bipolar outflows. These results support the idea that infall and outflow occur simultaneously in deeply embedded bipolar outflow sources. The data also indicate fairly large centrifugal radii and large inclinations to the rotational axis for a typical source.

Kenyon, Scott J.↗

The effects of physical aging at elevated temperatures on the viscoelastic creep on IM7/K3B

Physical aging at elevated temperature of the advanced composite IM7/K3B was investigated through the use of creep compliance tests. Testing consisted of short term isothermal, creep/recovery with the creep segments performed at constant load. The matrix dominated transverse tensile and in-plane shear behavior were measured at temperatures ranging from 200 to 230 C. Through the use of time based shifting procedures, the aging shift factors, shift rates and momentary master curve parameters were found at each temperature. These material parameters were used as input to a predictive methodology, which was based upon effective time theory and linear viscoelasticity combined with classical lamination theory. Long term creep compliance test data was compared to predictions to verify the method. The model was then used to predict the long term creep behavior for several general laminates.

Gates, Thomas S.↗

Localized tearing modes in the magnetotail driven by curvature effects

The stability of collisionless tearing modes is examined in the presence of curvature drift resonances and the trapped particle effects. A kinetic description for both electrons and ions is employed to investigate the stability of a two-dimensional equilibrium model. The main features of the study are to treat the ion dynamics properly by incorporating effects associated with particle trajectories in the tail fields and to include the linear coupling of trapped particle modes. Generalized dispersion relations are derived in several parameter regimes by considering two important sublayers of the reconnecting region. For a typical choice of parameters appropriate to the current sheet region, we demonstrate that localized tearing modes driven by ion curvature drift resonance effects are excited in the current sheet region with growth time of the order of a few seconds. Also, we examine nonlocal characteristics of tearing modes driven by curvature effects and show that modes growing in a fraction of a second arise when mode widths are larger than the current sheet width. Further, we show that trapped particle effects, in an interesting frequency regime, significantly enhance the growth rate of the tearing mode. The relevance of this theory for substorm onset phase and other features of the substorms is briefly discussed.

Sundaram, A. K.↗

An interleave retrieval method for oversampled limb scans

Rodgers presents a general approach for error analysis of retrievals based on a linearization of the retrieval problem. This approach allows calculation of retrieval error from uncertainties in the input radiances and in each of the forward radiance model parameters. These terms are weighted according to the sensitivity of the retrieval to each of them. We use this framework to describe the results of an error analysis performed on HALOE data. This paper will show the effect of using over sampled profile data and an interleave retrieval technique by application of formal error analysis.

Gordley, Larry L.↗

Energy spectrum transfer equations of solar wind turbulence

The recent studies of transfer equations for solar wind magnetohydrodynamic (MHD) turbulence are reviewed with emphasis on the comparison with the statistical observational results. Helios and Voyager missions provide an opportunity to study the the radial evolution of the power spectrum. the cross-helicity the Alfven ratio and the minimum variance direction. Spectrum transfer equations are considered as a tool to explore the nature of this radial evolution of the fluctuations. The transfer equations are derived from incompressible MHD equations. Generally one needs to make assumptions about the nature of the fluctuations and the nature of the turbulent non-linear interactions to obtain numerical results which can be compared with the observations. Some special model results for several simple cases SUCH as for structures or strong mixing. for Alfven waves with weak turbulent interactions. and for a superposition of structures and Alfven waves. are discussed. The difference between the various approaches to derive and handle the transfer equations are also addressed. Finally some theoretical description of the compressible fluctuations are also briefly reviewed.

Tu, C.-Y.↗

Model Checking a Byzantine-Fault-Tolerant Self-Stabilizing Protocol for Distributed Clock Synchronization Systems

This report presents the mechanical verification of a simplified model of a rapid Byzantine-fault-tolerant self-stabilizing protocol for distributed clock synchronization systems. This protocol does not rely on any assumptions about the initial state of the system. This protocol tolerates bursts of transient failures, and deterministically converges within a time bound that is a linear function of the self-stabilization period. A simplified model of the protocol is verified using the Symbolic Model Verifier (SMV) [SMV]. The system under study consists of 4 nodes, where at most one of the nodes is assumed to be Byzantine faulty. The model checking effort is focused on verifying correctness of the simplified model of the protocol in the presence of a permanent Byzantine fault as well as confirmation of claims of determinism and linear convergence with respect to the self-stabilization period. Although model checking results of the simplified model of the protocol confirm the theoretical predictions, these results do not necessarily confirm that the protocol solves the general case of this problem. Modeling challenges of the protocol and the system are addressed. A number of abstractions are utilized in order to reduce the state space. Also, additional innovative state space reduction techniques are introduced that can be used in future verification efforts applied to this and other protocols.

Malekpour, Mahyar R.↗

Verification of a Byzantine-Fault-Tolerant Self-stabilizing Protocol for Clock Synchronization

This paper presents the mechanical verification of a simplified model of a rapid Byzantine-fault-tolerant self-stabilizing protocol for distributed clock synchronization systems. This protocol does not rely on any assumptions about the initial state of the system except for the presence of sufficient good nodes, thus making the weakest possible assumptions and producing the strongest results. This protocol tolerates bursts of transient failures, and deterministically converges within a time bound that is a linear function of the self-stabilization period. A simplified model of the protocol is verified using the Symbolic Model Verifier (SMV). The system under study consists of 4 nodes, where at most one of the nodes is assumed to be Byzantine faulty. The model checking effort is focused on verifying correctness of the simplified model of the protocol in the presence of a permanent Byzantine fault as well as confirmation of claims of determinism and linear convergence with respect to the self-stabilization period. Although model checking results of the simplified model of the protocol confirm the theoretical predictions, these results do not necessarily confirm that the protocol solves the general case of this problem. Modeling challenges of the protocol and the system are addressed. A number of abstractions are utilized in order to reduce the state space.

Malekpour, Mahyar R.↗