Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Covariance Matrix”

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 487 records · Page 27

140,142 Ce Neutron Cross Section Resolved Resonance Region Evaluation [Abstract]

A resolved resonance region evaluation of 140,142 Ce has been carried out by Oak Ridge National Laboratory. Requested by the US Nuclear Criticality Safety Program, this evaluation is based on recent high-resolution transmission and capture high-resolution measurements of nat Ce and 142 Ce conducted at JRC-GEEL at the Geel Linear Accelerator facility, as well as recently measured thermal constants available from the EXFOR database. Starting from the resonance parameters from the ENDF/B-VIII.0 library and following a preliminary R-matrix analysis, an updated set of resonance parameters and corresponding covariance information was derived by the fit of these experimental datasets using the Reich-Moore approximation of the R-matrix theory as implemented in the SAMMY code system. The resolved resonance region upper energy limit for 140 Ce was kept at 200 keV while the 142 Ce resonance region was extended from 13 to 26 keV. This new evaluation was found to be in good agreement not only with several integral quantities of interest to the reactor physics community, but also with the stellar Maxwellian-averaged cross section.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measuring attitude with a gradiometer

This paper explores using a gravity gradiometer to measure the attitude of a satellite, given that the gravity field is accurately known. Since gradiometers actually measure a combination of the gradient and attitude rate and acceleration terms, the answer is far from obvious. The paper demonstrates that it can be done and at microradian accuracy. The technique employed is dynamic estimation, based on the momentum biased Euler equations. The satellite is assumed nominally planet pointed, and subject to control, gravity gradient, and partly radom drag torques. The attitude estimator is unusual. While the standard method of feeding back measurement residuals is used, the feedback gain matrix isn't derived from Kalman theory. instead, it's chosen to minimize a measure of the terminal covariance of the error in the estimate. This depends on the gain matrix and the power spectra of all the process and measurement noises. An integration is required over multiple solutions of Lyapunov equations.

Sonnabend, David↗

Bayesian Monte Carlo Evaluation of Imperfect (n, 233 U) Data and Model

Conventional nuclear data evaluation methods using generalized linear least squares make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate to minimize the cost function (even for nonlinear models); the model (e.g., of neutron cross section) and experimental data (including covariance data) are without defect and prior PDFs of parameters and measured data are known perfectly. Neglect of covariance between model parameters and measured data in conventional evaluations contributes to imperfections. These assumptions are inherent to the generalized linear least squares minimization method commonly used for resolved resonance region neutron cross section evaluations but are often not justified due to the presence of non-normal PDFs, nonlinear models (e.g., R-matrix formalism), and inherent imperfections in data and models (e.g., imperfect covariance data). Here, these assumptions are removed in a mathematical framework of Bayes’ theorem, which is implemented using the Metropolis-Hastings Monte Carlo method. Most importantly, new parameters are introduced to parameterize discrepancies between the theoretical model and measured data to quantify judgement about discrepancies or imperfections in a reproducible manner. An evaluation of 233U in the eV region using the ENDF-B/VIII.0 library and transmission data (Guber et al.) is presented, and posterior parameters are compared to those obtained by conventional evaluation methods. This example illustrates the effects of removing the most harmful assumption: that of model-data perfection.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Progress on the R-matrix Analysis for the n+ 181 Ta Evaluation [Slides]

This presentation discusses how nuclear data evaluations can influence the accuracy of computational nuclear design calculations. It also discusses the motivation for a new evaluation of 181 Ta. The goals of the evaluations are to generate a set of resonance parameters in the Resolved Resonance Region (RRR) that extends up to 2.6 keV, to generate a set of average resonance parameters in the Unresolved Resonance Region (URR) up to 100 keV (derived from the RRR parameters), and to generate covariance information for both RRR and URR. The presentation discusses how several different experimental data sets will be fit simultaneously in the evaluation of 181 Ta. It states that the level density of 181 Ta can make R-matrix fitting a challenge as the width of the levels increase with energy (as well as decreasing experimental resolution) and the measured resonances can overlap. The presentation concludes by discussing R-matrix analysis and covariance information as the latter is important when propagated to quantify the uncertainty on output responses of nuclear design calculations such as the multiplication factor k .

07 ISOTOPE AND RADIATION SOURCES↗

On a stochastic control method for weakly coupled linear systems

The stochastic control of two weakly coupled linear systems with different controllers is considered. Each controller only makes measurements about his own system; no information about the other system is assumed to be available. Based on the noisy measurements, the controllers are to generate independently suitable control policies which minimize a quadratic cost functional. To account for the effects of weak coupling directly, an approximate model, which involves replacing the influence of one system on the other by a white noise process is proposed. Simple suboptimal control problem for calculating the covariances of these noises is solved using the matrix minimum principle. The overall system performance based on this scheme is analyzed as a function of the degree of intersystem coupling.

Kwong, R. H.↗

Nuclear Criticality Safety Integral Experiment Covariance Determination

Integral benchmarks for criticality safety and nuclear data validation require expensive uncertainty quantification studies. Commonly, the uncertainty quantification ignores correlations between experiments that share components. Experiments such as the TEX (Thermal/Epithermal eXperiments) campaigns consist of many shared parts between experiments, such as fuel, which creates a strong correlation in their errors. While these correlations are known to exist, they are often not estimated due to the complexity of such calculations. This paper describes a software package that uses an intuitive method of determining the covariance for each of the experimental components, providing a correlation matrix for each family of parts across the multiple cases examined within a benchmark. The code uses the TEX-HEU campaign as a proof of concept, and we show that the correlations can be calculated with information commonly found in ICSBEP (International Criticality Safety Benchmark Evaluation Project) benchmarks. The estimated covariances are used in χ 2 trending studies to evaluate their impact on nuclear data validation. The covariance determination code can be easily integrated into current benchmark evaluations as well as reevaluating legacy benchmark uncertainties. Uncertainty correlation calculations should become the baseline for criticality safety integral experiment benchmarks and can now be easily calculated with the described software package.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Nuclear Criticality Safety Integral Experiment Covariance Determination

Integral benchmarks for criticality safety and nuclear data validation require expensive uncertainty quantification studies. Commonly, the uncertainty quantification ignores correlations between experiments that share components. Experiments such as the TEX (Thermal/Epithermal eXperiments) campaigns consist of many shared parts, such as fuel, which create a strong correlation in their uncertainties. While these correlations are known to exist, they are often not estimated due to the complexity of such calculations. This paper describes a software package that uses an intuitive method of determining the covariance for each of the experimental components, providing a correlation matrix for each family of parts across the multiple cases examined within a benchmark. The code uses the TEX-HEU campaign as a proof of concept, and we show that the correlations can be calculated with information commonly found in ICSBEP (International Criticality Safety Benchmark Evaluation Project) benchmarks. The estimated covariances are used in χ 2 trending studies to evaluate their impact on nuclear data validation. Without covariances, χ 2 per degree of freedom was calculated as 2.203 and with covariances it was 1.179. The difference shows that omitting covariance information may cause overly pessimistic bias quantifications. The covariance determination code can be easily integrated into current benchmark evaluations as well as reevaluating legacy benchmark uncertainties. Uncertainty correlation calculations should become the baseline for criticality safety integral experiment benchmarks and can now be easily calculated with the described software package.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Design and Implementation of a Parallel Multivariate Ensemble Kalman Filter for the Poseidon Ocean General Circulation Model

A multivariate ensemble Kalman filter (MvEnKF) implemented on a massively parallel computer architecture has been implemented for the Poseidon ocean circulation model and tested with a Pacific Basin model configuration. There are about two million prognostic state-vector variables. Parallelism for the data assimilation step is achieved by regionalization of the background-error covariances that are calculated from the phase-space distribution of the ensemble. Each processing element (PE) collects elements of a matrix measurement functional from nearby PEs. To avoid the introduction of spurious long-range covariances associated with finite ensemble sizes, the background-error covariances are given compact support by means of a Hadamard (element by element) product with a three-dimensional canonical correlation function. The methodology and the MvEnKF configuration are discussed. It is shown that the regionalization of the background covariances; has a negligible impact on the quality of the analyses. The parallel algorithm is very efficient for large numbers of observations but does not scale well beyond 100 PEs at the current model resolution. On a platform with distributed memory, memory rather than speed is the limiting factor.

Keppenne, Christian L.↗

Resolved Resonance Evaluation for Neutron Interactions with 103 Rh up to 8 keV

A neutron cross-section evaluation for the n + 103 Rh reaction in the resolved resonance region was carried out in the energy range 10−5 eV to 8 keV encompassing thermal energy at 0.0253 eV. The scope of this work is to generate resonance parameters and resonance parameter covariances based on the Reich-Moore reduced R-matrix formalism using the code SAMMY. Some features of the new evaluation are the inclusion of high-resolution capture data in the SAMMY evaluation process and the extension of the resolved resonance range from 4 to 8 keV. Furthermore, the evaluation employs more accurate resonance parameter representation by exploring the use of the LRF = 7 ENDF feature and also the use of the LCOMP = 2 compact format for resonance parameter covariance representation. Included in the SAMMY evaluation are transmission data, capture cross-section data, and neutron scattering length information. Thermal cross-section values listed in the literature, as well as capture resonance integrals, were also incorporated into the evaluation process.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A covariant multiple scattering series for elastic projectile-target scattering

A covariant formulation of the multiple scattering series for the optical potential is presented. The case of a scalar nucleon interacting with a spin zero isospin zero A-body target through meson exchange, is considered. It is shown that a covariant equation for the projectile-target t-matrix can be obtained which sums the ladder and crossed ladder diagrams efficiently. From this equation, a multiple scattering series for the optical potential is derived, and it is shown that in the impulse approximation, the two-body t-matrix associated with the first order optical potential is the one in which one particle is kept on mass-shell. The meaning of various terms in the multiple scattering series is given. The construction of the first-order optical potential for elastic scattering calculations is described.

Gross, Franz↗

Integrated Model Reduction and Control of Aircraft with Flexible Wings

This paper presents an integrated approach to the modeling and control of aircraft with exible wings. The coupled aircraft rigid body dynamics with a high-order elastic wing model can be represented in a nite dimensional state-space form. Given a set of desired output covariance, a model reduction process is performed by using the weighted Modal Cost Analysis (MCA). A dynamic output feedback controller, which is designed based on the reduced-order model, is developed by utilizing output covariance constraint (OCC) algorithm, and the resulting OCC design weighting matrix is used for the next iteration of the weighted cost analysis. This controller is then validated for full-order evaluation model to ensure that the aircraft's handling qualities are met and the uttering motion of the wings suppressed. An iterative algorithm is developed in CONDUIT environment to realize the integration of model reduction and controller design. The proposed integrated approach is applied to NASA Generic Transport Model (GTM) for demonstration.

AEROELASTIC↗

Results on the two population feature selection problem using probability of correct classification as a criterion

Variational equations are presented for maximizing the probability of correct classification as a function of a 1xn feature selection matrix B for the two-population problem. For the special case of equal covariance matrices the optimal B is unique up to scalar multiples and rank one sufficient. For equal population means, the best 1xn B is an eigenvector corresponding either to the largest or smallest eigenvalue of sigma sub 2 to the minus 1 power sigma sub 1 where sigma sub 1 and sigma sub 2 are the nxn covariance matrices of the two populations. The transformed probability of correct classification depends only on the eigenvalue. Finally, a procedure is proposed for constructing an optimal or nearly optimal kxn matrix of rank k without solving the k-dimensional variational equation.

Peters, B. C.↗

Computation of optimal output-feedback compensators for linear time-invariant systems

The control of linear time-invariant systems with respect to a quadratic performance criterion was considered, subject to the constraint that the control vector be a constant linear transformation of the output vector. The optimal feedback matrix, f*, was selected to optimize the expected performance, given the covariance of the initial state. It is first shown that the expected performance criterion can be expressed as the ratio of two multinomials in the element of f. This expression provides the basis for a feasible method of determining f* in the case of single-input single-output systems. A number of iterative algorithms are then proposed for the calculation of f* for multiple input-output systems. For two of these, monotone convergence is proved, but they involve the solution of nonlinear matrix equations at each iteration. Another is proposed involving the solution of Lyapunov equations at each iteration, and the gradual increase of the magnitude of a penalty function. Experience with this algorithm will be needed to determine whether or not it does, indeed, possess desirable convergence properties, and whether it can be used to determine the globally optimal f*.

Platzman, L. K.↗

Optimal member damper controller design for large space structures

Consideration is given to the selection of velocity feedback gains for individual dampers for the members of a structurally controlled large flexible space structure. The problem is formulated as an optimal output feedback regulator problem, and necessary conditions are derived for minimizing a quadratic performance function. The diagonal nature of the gain matrix is taken into account, along with knowledge of noise covariances. It is pointed out that the method presented offers a systematic approach to the design of a class of controllers for enhancing structural damping, which have significant potential if used in conjunction with a reduced-order optimal controller for rigid-body modes and selected structural modes.

Joshi, S. M.↗

Optoelectronic signal processing for phased-array antennas; Proceedings of the Meeting, Los Angeles, CA, Jan. 12, 13, 1988

Papers are presented on fiber optic links for airborne satellite applications, optoelectronic techniques for broadband switching, and GaAs circuits for a monolithic optical controller. Other topics include the optical processing of covariance matrices for adaptive processors, an optical linear heterodyne matrix-vector processor, and an EHF fiber optic-based array. An adaptive optical signal processing architecture using a signed-digit number system is considered along with microwave fiber optic links for phased arrays.

Bhasin, Kul B.↗

LMI-Based Multiobjective Optimization and Control of Flexible Aircraft Using VCCTEF

This paper considers the control of coupled aeroelastic aircraft model with Variable Camber Continuous Trailing Edge Flap (VCCTEF) system. The relative motion between two adjacent flaps is constrained and this actuation constraint problem is converted into an output covariance constraint problem, and therefore can be formulated using linear matrix inequalities (LMIs). A set of LMI conditions is derived for the design of an observer-based dynamic output feedback controller for VCCTEF configured aeroelastic aircraft model. The proposed controller is then applied to the NASA Generic Transport Model (GTM) for simulation, and the results demonstrate the efficacy of the proposed approach.

LMI control↗

Modeling and Control of Robot-Structure Coupling During In-Space Structure Assembly

This paper considers the problem of robot-structure coupling dynamics during in-space robotic assembly of large flexible structures. A two-legged walking robot is used as a construction agent, whose primary goal is to stably walking on the flexible structure while carrying a substructure component to a designated location. The reaction forces inserted by the structure to the walking robot are treated as bounded disturbance inputs, and a trajectory tracking robotic controller is proposed that combines the standard full state feedback motion controller and an adaptive controller to account for the disturbance inputs. In this study, a reduced-order Euler-Bernoulli beam structure model is adapted, and a finite number of co-located sensors and actuators are distributed along the span of the beam structure. The robot-structure coupling forces are treated as a bounded external forcing function to the structure, and hence an output covariance constraint problem can be formulated, in terms of linear matrix inequality, for optimal structure control by utilizing the direct output feedback controllers. The numerical simulations show the effectiveness of the proposed robot-structure modeling and control methodology.

Swei, Sean↗