A Robust Control Approach for the Formation Flying of Multiple Spacecraft
We present a disturbance rejection mechanism for the formation flying of multiple spacecraft based on a robust control approach in terms of an H control problem.
formation↗
SEARCH · Engineering Papers
Results for “matrix multiplication”
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.
Direct quarkonium production in DIS from a joint CGC and NRQCD framework
We compute the differential cross section for direct quarkonium production in high-energy electron-nucleus collisions at small 𝑥. Our computation is performed within the nonrelativistic QCD factorization formalism that separates the calculation into short distance coefficients and long distance matrix elements that depend on the color and spin of the state. We obtain the short distance coefficients of the production of the heavy quark pair within the framework of the color glass condensate effective field theory, which resums coherent multiple interactions of the heavy quark pair with the nucleus to all orders. Our results are expressed as the convolution of perturbatively calculable functions with multipoint lightlike Wilson line correlators. In the correlation limit, we establish the correspondence between our color glass condensate formulation with calculations employing the transverse momentum dependent (TMD) framework. We extend this correspondence by resumming kinematic power corrections within the improved TMD framework, which interpolates between the TMD formalism and 𝑘 ⊥ -factorization formalism. We present a detailed numerical analysis, focusing on 𝐽/𝜓 production in the kinematics accessible at the future Electron-Ion Collider, highlighting the importance of genuine higher-order saturation contributions when the electron collides with a large nucleus. Our results are also valid in the photoproduction limit where we expect the largest contribution from genuine higher-order saturation contributions which could be accessed in ultraperipheral collisions of relativistic heavy ions.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Simulation of Deuterium Oxidation in a Room with a Fire
A series of experiments where deuterium was released in trace amounts into a room with a fire were performed at Sandia National Laboratories. This report describes the corresponding effort to model the test series using SIERRA/FUEGO. The objective of this modeling effort was to produce a simulation test matrix that can be utilized to help interpret the corresponding experiments and be used to assess the credibility of using the SIERRA/FUEGO simulations as a surrogate for real tritium reaction data for mock fire and release studies. The simulations assessed multiple configurations of fire size and location and overall found that there was low conversion of deuterium to oxidized deuterium. There are some notable differences between the simulation predictions and the experimental results, such as the simulations generally predicting higher amounts of deuterium conversion compared to the experiment. A strategy to address modeling uncertainties and unknowns is also given in this report.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Ultra high tip speed (670.6 m/sec) fan stage with composite rotor: Aerodynamic and mechanical design
A highly loaded, single-stage compressor having a tip speed of 670.6 m/sec was designed for the purpose of investigating very high tip speeds and high aerodynamic loadings to obtain high stage pressure ratios at acceptable levels of efficiency. The design pressure ratio is 2.8 at an adiabatic efficiency of 84.4%. Corrected design flow is 83.4 kg/sec; corrected design speed is 15,200 rpm; and rotor inlet tip diameter is 0.853 m. The rotor uses multiple-circular-arc airfoils from 0 to 15% span, precompression airfoils assuming single, strong oblique shocks from 21 to 43% span, and precompression airfoils assuming multiple oblique shocks from 52% span to the tip. Because of the high tip speeds, the rotor blades are designed to be fabricated of composite materials. Two composite materials were investigated: Courtaulds HTS graphite fiber in a Kerimid 601 polyimide matrix and the same fibers in a PMR polyimide matrix. In addition to providing a description of the aerodynamic and mechanical design of the 670.0 m/sec fan, discussion is presented of the results of structural tests of blades fabricated with both types of matrices.
Halle, J. E.↗
Comparison of Two Algebraic Methods for Curve/curve Intersection
Most geometric modeling systems use either polynomial or rational functions to represent geometry. In such systems most computational problems can be formulated as systems of polynomials in one or more variables. Classical elimination theory can be used to solve such systems. Here Cayley's method of elimination is summarized and it is shown how it can best be used to solve the curve/curve intersection problem. Cayley's method was found to be a more straightforward approach. Furthermore, it is computationally simpler, since the elements of the Cayley matrix are one variable instead of two variable polynomials. Researchers implemented and tested both methods and found Cayley's to be more efficient. Six pairs of curves, representing mixtures of lines, circles, and cubic arcs were used. Several examples had multiple intersection points. For all six cases Cayley's required less CPU time than the other method. The average time ratio of method 1 to method 2 was 3.13:1, the least difference was 2.33:1, and the most dramatic was 6.25:1. Both of the above methods can be extended to solve the surface/surface intersection problem.
Demontaudouin, Y.↗
Landsat D Thematic Mapper image dimensionality reduction and geometric correction accuracy
To characterize and quantify the performance of the Landsat thematic mapper (TM), techniques for dimensionality reduction by linear transformation have been studied and evaluated and the accuracy of the correction of geometric errors in TM images analyzed. Theoretical evaluations and comparisons for existing methods for the design of linear transformation for dimensionality reduction are presented. These methods include the discrete Karhunen Loeve (KL) expansion, Multiple Discriminant Analysis (MDA), Thematic Mapper (TM)-Tasseled Cap Linear Transformation and Singular Value Decomposition (SVD). A unified approach to these design problems is presented in which each method involves optimizing an objective function with respect to the linear transformation matrix. From these studies, four modified methods are proposed. They are referred to as the Space Variant Linear Transformation, the KL Transform-MDA hybrid method, and the First and Second Version of the Weighted MDA method. The modifications involve the assignment of weights to classes to achieve improvements in the class conditional probability of error for classes with high weights. Experimental evaluations of the existing and proposed methods have been performed using the six reflective bands of the TM data. It is shown that in terms of probability of classification error and the percentage of the cumulative eigenvalues, the six reflective bands of the TM data require only a three dimensional feature space. It is shown experimentally as well that for the proposed methods, the classes with high weights have improvements in class conditional probability of error estimates as expected.
Ford, G. E.↗
A new approach to mixed H2/H infinity controller synthesis using gradient-based parameter optimization methods
In the past few years, the mixed H(sub 2)/H-infinity control problem has been the object of much research interest since it allows the incorporation of robust stability into the LQG framework. The general mixed H(sub 2)/H-infinity design problem has yet to be solved analytically. Numerous schemes have considered upper bounds for the H(sub 2)-performance criterion and/or imposed restrictive constraints on the class of systems under investigation. Furthermore, many modern control applications rely on dynamic models obtained from finite-element analysis and thus involve high-order plant models. Hence the capability to design low-order (fixed-order) controllers is of great importance. In this research a new design method was developed that optimizes the exact H(sub 2)-norm of a certain subsystem subject to robust stability in terms of H-infinity constraints and a minimal number of system assumptions. The derived algorithm is based on a differentiable scalar time-domain penalty function to represent the H-infinity constraints in the overall optimization. The scheme is capable of handling multiple plant conditions and hence multiple performance criteria and H-infinity constraints and incorporates additional constraints such as fixed-order and/or fixed structure controllers. The defined penalty function is applicable to any constraint that is expressible in form of a real symmetric matrix-inequity.
Ly, Uy-Loi↗
Expendable launch vehicle studies
Analytical support studies of expendable launch vehicles concentrate on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed and linearized. The information from LeRC about the grids, masses, and modes is incorporated into the model. The eigenvalues of the plant are analyzed for several modeling factors: utilizing diagonal mass matrix, uniform beam assumption, inclusion of aerodynamics, and the interaction between the aerodynamics and the flexible bending motion. Preliminary PID, LQR, and LQG control designs with sensor and actuator dynamics for this system and simulations are also conducted. The initial analysis for comparison of PD (proportional-derivative) and full state feedback LQR Linear quadratic regulator) shows that the split weighted LQR controller has better performance than that of the PD. In order to meet both the performance and robustness requirements, the H(sub infinity) robust controller for the expendable launch vehicle is developed. The simulation indicates that both the performance and robustness of the H(sub infinity) controller are better than that for the PID and LQG controllers. The modelling and analysis support studies team has continued development of methodology, using eigensensitivity analysis, to solve three classes of discrete eigenvalue equations. In the first class, the matrix elements are non-linear functions of the eigenvector. All non-linear periodic motion can be cast in this form. Here the eigenvector is comprised of the coefficients of complete basis functions spanning the response space and the eigenvalue is the frequency. The second class of eigenvalue problems studied is the quadratic eigenvalue problem. Solutions for linear viscously damped structures or viscoelastic structures can be reduced to this form. Particular attention is paid to Maxwell and Kelvin models. The third class of problems consists of linear eigenvalue problems in which the elements of the mass and stiffness matrices are stochastic. dynamic structural response for which the parameters are given by probabilistic distribution functions, rather than deterministic values, can be cast in this form. Solutions for several problems in each class will be presented.
Bainum, Peter M.↗
Investigation of Photoluminescence and Photocurrent in InGaAsP/InP Strained Multiple Quantum Well Heterostructures
Multiple quantum well InGaAsP/InP p-i-n laser heterostructures with different barrier thicknesses have been investigated using photoluminescence (PL) and photocurrent (PC) measurements. The observed PL spectrum and peak positions are in good agreement with those obtained from transfer matrix calculations. Comparing the measured quantum well PC with calculated carrier escape rates, the photocurrent changes are found to be governed by the temperature dependence of the electron escape time.
Raisky, O. Y.↗
FIB-TEM Investigations of Fe-NI-Sulfides in the CI Chondrites Alais and Orgueil
The CI chondrites are primitive meteorites with bulk compositions matching the solar photosphere for all but the lightest elements. They have been extensively aqueously altered, and are composed primarily of fine-grained phyllosilicate matrix material which is host to carbonates, sulfates, sulfides, and minor amounts of olivine and pyroxene. The alteration, while extensive, is heterogeneous. For example, CI-chondrite cubanite and carbonate grains differ on mm to sub-mm scales, demonstrating multiple aqueous episodes. CI-chondrite variability is also evidenced by degree of brecciation, abundance and size of coarse-grained phyllosilicates, olivine and pyroxene abundance, as well as Ni-content and size of sulfide grains. Our previous work revealed Orgueil sulfide grains with variable Ni-contents, metal:S ratios, crystal structures and textures. We continue to explore the variability of CI-chondrite pyrrhotite (Po, (FeNi)1-xS) and pentlandite (Pn, (Fe,Ni)9S8) grains. We investigate the microstructure of sulfides within and among CI-chondrite meteorites in order to place constraints on the conditions under which they formed.
Berger, Eve L.↗
A method for characterization of multiple dynamic constitutive parameters of FRCs
We propose a method to measure multiple dynamic material constitutive parameters of unidirectional fiber reinforced composites (FRCs) in a single experiment. Dynamic short-beam shear (DSBS) experiments were performed on a modified Kolsky compression bar, with integration of high-speed imaging and digital image correlation (DIC). The unidirectional FRCs investigated were S-2 glass fiber reinforced matrix of TGDDM-Jeffamine® D230 with monoamine functionalized partially reacted substructures (mPRS) and commercially available SC-15. Analytical solutions of normal and shear strains of a composite beam were derived based on Timoshenko beam theory, assuming material to be transversely isotropic and have different moduli in tension and compression in each principle material orientation. Tensile and compressive moduli were inversely computed through monitoring normal strain slope when specimen was constantly loaded at a speed of ~7.3 m/s within a small deflection. Non-linear shear stress-strain behavior of the composite was described via Ramberg–Osgood equation. Finite element (FE) analysis was conducted in ABAQUS, simultaneously defining via user subroutine UMAT the transverse isotropy of material, bi-modulus constitutive model, and non-linear shear stress-strain relation. The method proposed in this work was validated by comparing strain distributions computed by FE model and DIC measurements. Comparing with traditional dynamic tensile, compressive, and shear experiments on FRCs, this method significantly simplifies the specimen preparation and design of complicated gripping fixtures for multiple experiments. Here, systematic errors resulting from variations of specimen geometry and dimension, loading direction, and instrumentation are reduced, thereby providing compatible data for numerical studies on impact behavior of composites.
36 MATERIALS SCIENCE↗
Multiple input/output random vibration control system
A multi-input/output random vibration control algorithm was developed based on system identification concepts derived from random vibration spectral analysis theory. The unique features of the algorithm are: (1) the number of input excitors and the number of output control responses need not be identical; (2) the system inverse response matrix is obtained directly from the input/output spectral matrix; and (3) the system inverse response matrix is updated every control loop cycle to accommodate system amplitude nonlinearities. A laboratory demonstration case of two imputs with three outputs is presented to demonstrate the system capabilities.
Unruh, James F.↗
Evading strength and ductility trade-off in an inverse nacre structured magnesium matrix nanocomposite
Mimicking nacre structure is an often-used strategy for developing materials which combine high strength and high toughness. Here, a novel strategy, named inverse nacre structure, is proposed where elongated and curled soft constituents are embedded in a matrix of hard constituent to form an ordered brick-and-mortar arrangement. This strategy is demonstrated in a particle reinforced metal matrix composite, using pure Mg as the soft phase and SiC nanoparticle reinforced Mg as the hard phase. In this work, the resulting nanocomposite yields high strength, and especially high tensile elongation which is five times that of the homogeneous counterpart. The superior comprehensive strengthening-toughening effect originates from the tailored inverse nacre structure which enables a stable microcrack-multiplication process by sufficiently promoting strain hardening so as to avoid any significant decrease in post-ultimate tensile strength. This strategy, while demonstrated here for the specific case of particle reinforced Mg, is not restricted to any particular material system but constitutes a generic pathway for the development of high-performance materials for engineering applications.
36 MATERIALS SCIENCE↗
Skin factor and potential formation damage from chemical and mechanical processes in a naturally fractured carbonate aquifer with implications to CO 2 sequestration
Here, in this study, we investigate formation damage due to acidization and water injection tests into the naturally fractured carbonate Middle Duperow Formation at Kevin Dome, Montana, potentially diminishing the chance of a future successful Geological Carbon Sequestration (GCS) project. Multiple well-test analytical models, correlated with core description and lithology data, are used to determine flow behavior and communication between the water injection interval and surrounding formations. An improved three-dimensional (3D) geologic model with dual-continuum matrix and fracture properties is constructed based on most recent seismic, core, and water sample measurements. Brine injection is simulated to verify the interpretation from the analytical models, followed by CO 2 injection simulation. Geochemical calculations are performed to understand the in-situ processes that led to formation damage. Our findings suggest: (1) there are two possible scenarios that could lead to a positive total effective skin factor and permeability decline: partial penetration and formation damage; (2) analytical models indicate a positive total skin factor, contradicting results of a previous study suggesting that the well was mildly stimulated; (3) numerical simulation supports the formation damage hypothesis by matching the pressure buildup observed during the latter two brine injection tests; (4) several mechanical and chemical processes may have occurred during injection to clog the matrix/fracture system: anhydrite fines migration and/or calcite precipitation. We then make preventative suggestions for future GCS projects into carbonate reservoirs and remediation recommendations for GCS operation at the Kevin Dome.
54 ENVIRONMENTAL SCIENCES↗
Tailoring light-induced charge transfer and intersystem crossing in FeCO using time-dependent spin–orbit configuration interaction
Real-time (RT) electronic structure methods provide a natural framework for describing light–matter interactions in arbitrary time-dependent electromagnetic fields (EMF). Optically induced excited state transitions are of particular interest, which require tuned EMF to drive population transfer to and from the specific state(s) of interest. Intersystem crossing, or spin-flip, may be driven through shaped EMF or laser pulses. These transitions can result in long-lived “spin-trapped” excited states, which are especially useful for materials requiring charge separation or protracted excited state lifetimes. Time-dependent configuration interaction (TDCI) is unique among RT methods in that it may be implemented in a basis of eigenstates, allowing for rapid propagation of the time-dependent Schrödinger equation. The recent spin–orbit TDCI (TD-SOCI) enables a real-time description of spin-flip dynamics in an arbitrary EMF and, therefore, provides an ideal framework for rational pulse design. The present study explores the mechanism of multiple spin-flip pathways for a model transition metal complex, FeCO, using shaped pulses designed to drive controlled intersystem crossing and charge transfer. These results show that extremely tunable excited state dynamics can be achieved by considering the dipole transition matrix elements between the states of interest.
Chemistry↗
Efficient simulation of open quantum systems coupled to a reservoir through multiple channels
It is challenging to simulate open quantum systems that are connected to a reservoir through multiple channels. For example, vibrations may induce fluctuations in both energy gaps and electronic couplings, which represent two independent channels of system–bath couplings. Systems of this kind are ubiquitous in the processes of excited state radiationless decay. Combined with density matrix renormalization group (DMRG) and matrix product states (MPS) methods, we develop an interaction-picture chain mapping strategy for vibrational reservoirs to simulate the dynamics of these open systems, resulting in time-dependent spatially local system–bath couplings in the chain-mapped Hamiltonian. This transformation causes the entanglement generated by the system–bath interactions to be restricted within a narrow frequency window of vibrational modes, enabling efficient DMRG/MPS dynamical simulations. We demonstrate the utility of this approach by simulating singlet fission dynamics using a generalized spin-boson Hamiltonian with both diagonal and off-diagonal system–bath couplings. This approach generalizes an earlier interaction-picture chain mapping scheme, allowing for efficient and exact simulation of systems with multi-channel system–bath couplings using matrix product states, which may further our understanding of nonlocal exciton–phonon couplings in exciton transport and the non-Condon effect in energy and electron transfer.
Chemistry↗
Quantum simulations of hadron dynamics in the Schwinger model using 112 qubits
Hadron wave packets are prepared and time evolved in the Schwinger model using 112 qubits of IBM’s 133-qubit Heron quantum computer ibm_torino. The initialization of the hadron wave packet is performed in two steps. First, the vacuum is prepared across the whole lattice using the recently developed SC-ADAPT-VQE algorithm and workflow. SC-ADAPT-VQE is then extended to the preparation of localized states, and used to establish a hadron wave packet on top of the vacuum. This is done by adaptively constructing low-depth circuits that maximize the overlap with an adiabatically prepared hadron wave packet. Due to the localized nature of the wavepacket, these circuits can be determined on a sequence of small lattices using classical computers, and then robustly scaled to prepare wave packets on large lattices for simulations using quantum computers. Time evolution is implemented with a second-order Trotterization. To reduce both the required qubit connectivity and circuit depth, an approximate quasilocal interaction is introduced. This approximation is made possible by the emergence of confinement at long distances, and converges exponentially with increasing distance of the interactions. Using multiple error-mitigation strategies, up to 14 Trotter steps of time evolution are performed, employing 13,858 two-qubit gates (with a CNOT depth of 370). The propagation of hadrons is clearly identified, with results that compare favorably with Matrix Product State simulations. Finally, prospects for a near-term quantum advantage in simulations of hadron scattering are discussed.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
On obtaining the forward phase functions of Saturn ring features from radio occultation observations
It is noted that the near-forward scattering functions of particles in Saturn ring features are related to 3.6 cm radio occultation power spectra by a Fredholm integral equation of the first kind. The equation reduces to an algebraic system of equation whose solution by usual inversion techniques (that is, least mean squares) is ruled out by the near singularity of the forward transformation matrix. A combination of constrained linear inversion and a filtering algorithm based on eigenvector decomposition of the matrix reduces the instabilities; this yields derived phase functions valid over the range of zero to about 12 mrad. It is noted that these functions represent the collective forward diffraction lobe of particles greater than about 1 m in radius. Since multiple scattering of the signal is a significant effect, the measured phase functions must be adjusted to obtain the singly scattered component. This single-scattering correction is examined for two physical models, namely the monolayer and the classical discrete random slab; in addition, the fraction of opacity in submeter particles for each model for particular ring features is estimated.
Zebker, H. A.↗