Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Derived Matrices”

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 289 records · Page 16

Examining the justification for the introduction of a fermion localization function

Becke and Edgecombe suggested in 1990 a theoretical tool to describe electron localization in atoms and molecules, an idea which was borrowed by a large number of nuclear theorists since 2011 to describe nucleon localization in nuclear systems. I argue here that these arguments are highly questionable and cannot be used in interacting systems where effects beyond the naive mean field or the simple Hartree-Fock framework are important and the inclusion of correlations induced by particle interactions is necessary in order to introduce such a localization function. Furthermore, I also describe several aspects of the exchange and irreducible two-body density matrices, which depend on the character and strength of the two-particle interaction and which can be useful in justifying the derivation of an appropriate energy density functional.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

State-Space System Realization with Input- and Output-Data Correlation

This paper introduces a general version of the information matrix consisting of the autocorrelation and cross-correlation matrices of the shifted input and output data. Based on the concept of data correlation, a new system realization algorithm is developed to create a model directly from input and output data. The algorithm starts by computing a special type of correlation matrix derived from the information matrix. The special correlation matrix provides information on the system-observability matrix and the state-vector correlation. A system model is then developed from the observability matrix in conjunction with other algebraic manipulations. This approach leads to several different algorithms for computing system matrices for use in representing the system model. The relationship of the new algorithms with other realization algorithms in the time and frequency domains is established with matrix factorization of the information matrix. Several examples are given to illustrate the validity and usefulness of these new algorithms.

Juang, Jer-Nan↗

A reduced order observer based controller design for H-infinity optimization

In this paper the H-infinity control problem is investigated. It is well-known that for this problem, in general, controllers of the same dynamic order as the given system are needed. However, in the case that the standard assumptions on two direct feedthrough matrices are not satisfied, it is here shown that dynamic compensators of a lower dynamical order can be found. This result can be derived by using the standard reduced order observer based controllers in the case that one or more states are measured without noise.

Stoorvogel, Anton A.↗

Computing the projected area of a spacecraft solar array with matrices

Matrix methods for computing the projected area of a solar array as a function of the rotational position of a spacecraft and of array position on a spacecraft are presented. Formulas are derived which provide the optimum solar array pitch, cant and tracking angles for a given spacecraft configuration and orbit. These formulas are general and applicable to many spacecraft. Formulas are also provided for determining the energy output from an array for a given orbit and the resultant energy available to spacecraft loads. Results are simply obtained and realizable with a hand calculator. The methods above can be extended to the case of computer analysis of solar array shadowing. How this can be done is outlined and results are presented from a spacecraft study program.

E. M. Gaddy↗

On the symbolic manipulation and code generation for elasto-plastic material matrices

A computerized procedure for symbolic manipulations and FORTRAN code generation of an elasto-plastic material matrix for finite element applications is presented. Special emphasis is placed on expression simplifications during intermediate derivations, optimal code generation, and interface with the main program. A systematic procedure is outlined to avoid redundant algebraic manipulations. Symbolic expressions of the derived material stiffness matrix are automatically converted to RATFOR code which is then translated into FORTRAN statements through a preprocessor. To minimize the interface problem with the main program, a template file is prepared so that the translated FORTRAN statements can be merged into the file to form a subroutine (or a submodule). Three constitutive models; namely, von Mises plasticity, Drucker-Prager model, and a concrete plasticity model, are used as illustrative examples.

Chang, T. Y.↗

Polymeric precursors for fibers and matrices

Candidate polymeric precursors for ceramic fiber and matrix processing are discussed, with a view to the advantages and disadvantages of this approach relative to existing alternatives. The properties of ceramic products thus derived are noted to strongly depend on the molecular weight and structure of the starting polymer; in particular, the ceramic's composition and morphology are dependent on the character and extent of crosslinking, as well as on the path of pyrolysis. While large and complex structural ceramic components may ultimately be obtainable by these means, the polymer-precursor method is still in its developmental infancy.

Hurwitz, Frances I.↗

Symbolic generation of elastic rotor blade equations using a FORTRAN processor and numerical study on dynamic inflow effects on the stability of helicopter rotors

The process of performing an automated stability analysis for an elastic-bladed helicopter rotor is discussed. A symbolic manipulation program, written in FORTRAN, is used to aid in the derivation of the governing equations of motion for the rotor. The blades undergo coupled bending and torsional deformations. Two-dimensional quasi-steady aerodynamics below stall are used. Although reversed flow effects are neglected, unsteady effects, modeled as dynamic inflow are included. Using a Lagrangian approach, the governing equations are derived in generalized coordinates using the symbolic program. The program generates the steady and perturbed equations and writes into subroutines to be called by numerical routines. The symbolic program can operate on both expressions and matrices. For the case of hovering flight, the blade and dynamic inflow equations are converted to equations in a multiblade coordinate system by rearranging the coefficients of the equations. For the case of forward flight, the multiblade equations are obtained through the symbolic program. The final multiblade equations are capable of accommodating any number of elastic blade modes. The computer implementation of this procedure consists of three stages: (1) the symbolic derivation of equations; (2) the coding of the equations into subroutines; and (3) the numerical study after identifying mass, damping, and stiffness coefficients. Damping results are presented in hover and in forward flight with and without dynamic inflow effects for various rotor blade models, including rigid blade lag-flap, elastic flap-lag, flap-lag-torsion, and quasi-static torsion. Results from dynamic inflow effects which are obtained from a lift deficiency function for a quasi-static inflow model in hover are also presented.

Reddy, T. S. R.↗

Assessing Gas-Phase Ion Reactivity of 50 Elements with NO and the Direct Application for 239 Pu in Complex Matrices Using ICP-MS/MS

Understanding the reactivity of metal cations with various reaction gases in ICP-MS/MS is important to determine the best gas to use for a given analyte/interference pair. In this study, nitric oxide (NO) was investigated as the reaction gas following previous experimental designs. The reactions with 50 elements were investigated to examine periodic trends in reactivity, validate theoretical modeling of reaction enthalpies as a method to screen reactant gases, and provide a baseline for potential in-line gas separation methods. ICP-MS/MS studies involving actinides are typically limited to Th, U and Pu, with analyses of Np and Am rarely reported in the literature. To date, only two previous methods have investigated the use of NO in ICP-MS/MS analyses. To showcase the utility of NO, a method was developed to measure 239 Pu in the presence of environmental matrix constituent and other actinides, like what could be expected from post-detonation debris, with no chemical separation prior to analysis. 239 Pu + was reacted to form 239 Pu 16 O + , eliminating interferences derived from the sample matrix by measuring 239 Pu + intensity at m/z = 255 ( 239 Pu 16 O + ). To validate NO for 238 U 1 H + interference removal in environmental matrices, standard reference materials were diluted to 1 mg per gram of solution and spiked to 0.05 pg/g of 239 Pu and 1 µg/g 238 U (Pu/U = 5 x 10 -8 ). Measured 239 Pu concentrations were within 6% of the spiked value. In conclusion, these results demonstrate that reliable 239 Pu measurements can be made at levels relevant to nuclear forensics without the need for extensive chemical matrix separation prior to analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Geometric conditions for saturating the data processing inequality

Abstract The data processing inequality (DPI) is a scalar inequality satisfied by distinguishability measures on density matrices. For some distinguishability measures, saturation of the scalar DPI implies an operator equation relating the arguments of the measure. These results are typically derived using functional analytic techniques. In a complementary approach, we use geometric techniques to derive a formula that gives an operator equation from DPI saturation for any distinguishability measure; moreover, for a broad class of distinguishability measures, the derived operator equation is sufficient to imply saturation as well. Our operator equation coincides with known results for the sandwiched Rényi relative entropies, and gives new results for α - z Rényi relative entropies and a family of of quantum f -divergences, which we compute explicitly.

Physics↗

BICEP/Keck. XX. Component-separated Maps of the Polarized Cosmic Microwave Background and Thermal Dust Emission Using Planck and BICEP/Keck Observations through the 2018 Observing Season

We present component-separated polarization maps of the cosmic microwave background (CMB) and Galactic thermal dust emission, derived using data from the BICEP/Keck experiments through the 2018 observing season and Planck. By employing a maximum-likelihood method that utilizes observing matrices, we produce unbiased maps of the CMB and dust signals. We outline the computational challenges and demonstrate an efficient implementation of the component map estimator. We show methods to compute and characterize power spectra of these maps, opening up an alternative way to infer the tensor-to-scalar ratio from our data. We compare the results of this map-based separation method with the baseline BICEP/Keck analysis. Our analysis demonstrates consistency between the two methods, finding an 84% correlation between the pipelines.

cosmic inflation↗

Multivariable control theory applied to hierarchial attitude control for planetary spacecraft

Multivariable control theory is applied to the design of a hierarchial attitude control system for the CARD space vehicle. The system selected uses reaction control jets (RCJ) and control moment gyros (CMG). The RCJ system uses linear signal mixing and a no-fire region similar to that used on the Skylab program; the y-axis and z-axis systems which are coupled use a sum and difference feedback scheme. The CMG system uses the optimum steering law and the same feedback signals as the RCJ system. When both systems are active the design is such that the torques from each system are never in opposition. A state-space analysis was made of the CMG system to determine the general structure of the input matrices (steering law) and feedback matrices that will decouple the axes. It is shown that the optimum steering law and proportional-plus-rate feedback are special cases. A derivation of the disturbing torques on the space vehicle due to the motion of the on-board television camera is presented. A procedure for computing an upper bound on these torques (given the system parameters) is included.

Boland, J. S., III↗

Necessary and sufficient conditions for the complete controllability and observability of systems in series using the coprime factorization of a rational matrix

The series connection of two linear time-invariant systems that have minimal state space system descriptions is considered. From these descriptions, strict-system-equivalent polynomial matrix system descriptions in the manner of Rosenbrock are derived. They are based on the factorization of the transfer matrix of the subsystems as a ratio of two right or left coprime polynomial matrices. They give rise to a simple polynomial matrix system description of the tandem connection. Theorem 1 states that for the complete controllability and observability of the state space system description of the series connection, it is necessary and sufficient that certain 'denominator' and 'numerator' groups are coprime. Consequences for feedback systems are drawn in Corollary 1. The role of pole-zero cancellations is explained by Lemma 3 and Corollaires 2 and 3.

Callier, F. M.↗

A higher-order finite-element implementation of the nonlinear Fokker–Planck collision operator for charged particle collisions in a low density plasma

Collisions between particles in a low density plasma are described by the Fokker–Planck collision operator. In applications, this nonlinear integro-differential operator is often approximated by linearised or ad-hoc model operators due to computational cost and complexity. In this work, we present an implementation of the nonlinear Fokker–Planck collision operator written in terms of Rosenbluth potentials in the Rosenbluth–MacDonald–Judd (RMJ) form. The Rosenbluth potentials may be obtained either by direct integration or by solving partial differential equations (PDEs) similar to Poisson's equation: we optimise for performance and scalability by using sparse matrices to solve the relevant PDEs. We represent the distribution function using a tensor-product continuous-Galerkin finite-element representation and we derive and describe the implementation of the weak form of the collision operator. We present tests demonstrating a successful implementation using an explicit time integrator and we comment on the speed and accuracy of the operator. Finally, we speculate on the potential for applications in the current and next generation of kinetic plasma models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Petrology of the Cangas de Onis and nulles regolith breccias Implications for parent body history

The present study of the Cangas de Onis and Nulles H chondrite regolith breccias indicates that the minerals in the matrices and equilibrated clasts have essentially the same compositional distributions, so that much of the material in the castic matrix would have to have been derived from the impact comminution of clats. The apparently exclusive occurrence of H6 clasts in Cangas de Onis, and H4 clasts in Nulles, suggests that, at the locations where these breccias formed, the regolith predominantly consisted of H6 and H4 material, respectively.

Williams, C. V.↗

On the Derivation of Quasi-Newton Formulas for Optimization in Function Spaces

Newton’s method is usually preferred when solving optimization problems due to its superior convergence properties compared to gradient-based or derivative-free optimization algorithms. However, deriving and computing second-order derivatives needed by Newton’s method often is not trivial and, in some cases, not possible. In such cases quasi-Newton algorithms are a great alternative. In this paper, we provide a new derivation of well-known quasi-Newton formulas in an infinite-dimensional Hilbert space setting. Furthermore, it is known that quasi-Newton update formulas are solutions to certain variational problems over the space of symmetric matrices. In this paper, we formulate similar variational problems over the space of bounded symmetric operators in Hilbert spaces. By changing the constraints of the variational problem we obtain updates (for the Hessian and Hessian inverse) not only for the Broyden-Fletcher-Goldfarb-Shanno (BFGS) quasi-Newton method but also for Davidon–Fletcher–Powell (DFP), Symmetric Rank One (SR1), and Powell-Symmetric-Broyden (PSB). In addition, for an inverse problem governed by a partial differential equation (PDE), we derive DFP and BFGS “structured” secant formulas that explicitly use the derivative of the regularization and only approximates the second derivative of the misfit term. We show numerical results that demonstrate the desired mesh-independence property and superior performance of the resulting quasi-Newton methods.

97 MATHEMATICS AND COMPUTING↗

Preliminary data on boulders at station 6, Apollo 17 landing site

A cluster of boulders at Station 6 (Apollo 17 landing site) consists of breccias derived from the North Massif. Three preliminary lithologic units were established, on the basis of photogeologic interpretations; all lithologies identified photogeologically were sampled. Breccia clasts and matrices studied petrographically and chemically fall into two groups by modal mineralogy: (1) low-K Fra Mauro or high basalt composition, consisting of 50-60% modal feldspar, approximately 45% orthopyroxene and 1-7% Fe-Ti oxide; (2) clasts consisting of highland basalt composition, consisting of 70% feldspar, 30% orthopyroxene and olivine and a trace of Fe-Ti oxide.

Heiken, G. H.↗

Infinite order sudden approximation for rotational energy transfer in gaseous mixtures

Rotational energy transfer in gaseous mixtures is analyzed within the framework of the infinite order sudden (IOS) approximation, and a new derivation of the IOS from the coupled states Lippman-Schwinger equation is presented. This approach shows the relation between the IOS and coupled state T matrices. The general IOS effective cross section can be factored into a finite sum of 'spectroscopic coefficients' and 'dynamical coefficients'. The evaluation of these coefficients is considered. Pressure broadening for the systems HD-He, HCl-He, CO-He, HCl-Ar, and CO2-Ar is calculated, and results based on the IOS approximation are compared with coupled state results. The IOS approximation is found to be very accurate whenever the rotor spacings are small compared to the kinetic energy, provided closed channels do not play too great a role.

Goldflam, R.↗

On the Use of Kronecker Operators for the Solution of Generalized Stochastic Petri Nets

We discuss how to describe the Markov chain underlying a generalized stochastic Petri net using Kronecker operators on smaller matrices. We extend previous approaches by allowing both an extensive type of marking-dependent behavior for the transitions and the presence of immediate synchronizations. The derivation of the results is thoroughly formalized, including the use of Kronecker operators in the treatment of the vanishing markings and the computation of impulse-based reward measures. We use our techniques to analyze a model whose solution using conventional methods would fail because of the state-space explosion. In the conclusion, we point out ideas to parallelize our approach.

Ciardo, Gianfranco↗