Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic matrix completion”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Development of an engineering analysis of progressive damage in composites during low velocity impact

A computerized, analytical methodology was developed to study damage accumulation during low velocity lateral impact of layered composite plates. The impact event was modeled as perfectly plastic with complete momentum transfer to the plate structure. A transient dynamic finite element approach was selected to predict the displacement time response of the plate structure. Composite ply and interlaminar stresses were computed at selected time intervals and subsequently evaluated to predict layer and interlaminar damage. The effects of damage on elemental stiffness were then incorporated back into the analysis for subsequent time steps. Damage predicted included fiber failure, matrix ply failure and interlaminar delamination.

Humphreys, E. A.↗
DOE OSTI2022

Hybrid quantum-classical simulations of magic angle spinning dynamic nuclear polarization in very large spin systems

Solid-state nuclear magnetic resonance can be enhanced using unpaired electron spins with a method known as dynamic nuclear polarization (DNP). Fundamentally, DNP involves ensembles of thousands of spins, a scale that is difficult to match computationally. This scale prevents us from gaining a complete understanding of the spin dynamics and applying simulations to design sample formulations. We recently developed an ab initio model capable of calculating DNP enhancements in systems of up to ~1000 nuclei; however, this scale is insufficient to accurately simulate the dependence of DNP enhancements on radical concentration or magic angle spinning (MAS) frequency. We build on this work by using ab initio simulations to train a hybrid model that makes use of a rate matrix to treat nuclear spin diffusion. We show that this model can reproduce the MAS rate and concentration dependence of DNP enhancements and build-up time constants. Finally, we then apply it to predict the DNP enhancements in core–shell metal-organic-framework nanoparticles and reveal new insights into the composition of the particles’ shells.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Gauge-fixing quantum density operators at scale

We provide a theory, algorithms, and simulations of nonequilibrium quantum systems using a one-dimensional (1D) completely positive (CP), matrix-product (MP) density-operator (𝜌) representation. By generalizing the matrix product state's orthogonality center, to additionally store positive classical mixture correlations, the MP⁢𝜌 factorization naturally emerges. In this setting, we analytically and numerically examine the virtual gauge freedoms associated with the representation of quantum density operators. Based on this perspective, we simplify algorithms in certain limits to speed up the integration of the canonical-form master-equation dynamics. This enables us to quickly evolve under the dynamics of two-body quantum channels without resorting to optimization-based methods. In addition to this technical advance, we also scale up numerical examples and discuss implications for accurately modeling hardware architectures and predicting their performance in the near term. This includes an example of the quantum to classical transition of informationally leaky, i.e., decohering, qubits. In this setting, because of loss from environmental interactions, nonlocal complex coherence correlations are converted into global incoherent classical statistical mixture correlations. Lastly, the representation of both global and local correlations is discussed. We expect this work to have applications in additional nonequilibrium settings, beyond qubit engineering.

Gangapuram, Amit Jamadagni [Oak Ridge National Lab↗
DOE OSTI2024

Flax fiber-reinforced fatty acid vitrimer biocomposite with enhanced chemical recyclability

Here, this study unveils a sustainable, easily recyclable biocomposite, leveraging the dynamic nature of covalently adaptive bonds in a vitrimer matrix. The fabrication involved a fatty acid-derived vitrimer as the polymer matrix and multi-layered, nonwoven flax mat as reinforcing scaffold. The incorporation of these fibers significantly improved the mechanical performance of the vitrimer matrix uniformly. The ester-based covalently adaptive network plays a crucial role in enabling exceptional fiber-matrix bonding, as well as recyclability. The vitrimer matrix dissolves in ethylene glycol through transesterification, facilitating complete material recovery and biocomposite recycling without compromising the original properties of the matrix and reinforcing fibers.

36 MATERIALS SCIENCE↗
NASA NTRS1984

A variational method for the calculation of dynamic polarizabilities and two-photon transition moments - The dressed molecule approach

The solution of the time-dependent Schroedinger equation of the molecule + radiation field system is analyzed. A quantized radiation field is used. The relationship between the oscillatory wave function and the dressed state wave function is established, the oscillatory wave function being the solution if the radiation field is turned on at t = 0 and the dressed state being the stationary solution satisfying the boundary condition only if the field is present at t = negative infinity. In general, the oscillatory wave function can be expressed using a complete set of dressed states. However, in the presence of a nonresonant radiation field, the system is well represented by a single dressed state. Molecular properties such as dynamic polarizabilities and two-photon transition moments can be deduced from the dressed state wave function instead of the oscillatory wave function as in previous methods. Because of its stationary character, the dressed state is more amenable to approximate calculations. A CI method is developed for this purpose. The CI matrix is simple in form and can be readily constructed using existing computer codes. The present method can also be adapted to calculate other optical properties.

Huo, W. M.↗
NASA NTRS1977

Error rate information in attention allocation pilot models

The Northrop urgency decision pilot model was used in a command tracking task to compare the optimized performance of multiaxis attention allocation pilot models whose urgency functions were (1) based on tracking error alone, and (2) based on both tracking error and error rate. A matrix of system dynamics and command inputs was employed, to create both symmetric and asymmetric two axis compensatory tracking tasks. All tasks were single loop on each axis. Analysis showed that a model that allocates control attention through nonlinear urgency functions using only error information could not achieve performance of the full model whose attention shifting algorithm included both error and error rate terms. Subsequent to this analysis, tracking performance predictions for the full model were verified by piloted flight simulation. Complete model and simulation data are presented.

Faulkner, W. H.↗
DOE OSTI2025

Friction surface layer deposition of triple-phase Al 10 Cr 12 Fe 35 Mn 23 Ni 20 high entropy alloy: Process optimization and microstructural evolution

A high-strength Co-free triple-phase Al 10 Cr 12 Fe 35 Mn 23 Ni 20 high-entropy alloy (HEA) was successfully fabricated using Friction Surface Layer Deposition (FSLD), a bulk manufacturing method. Multiple single-layer deposits were produced by varying forging force (F) and traverse speeds (T r ) to optimize the process parameters. The optimized conditions (F = 40 kN & T r = 200 mm/min) were then applied to manufacture a scaled-up multi-layer specimen. The initial microstructure of the HEA consisted of coarse grains of the soft FCC-phase, long columnar dendrites of the hard BCC-phase, and small precipitates of the harder B2-phase within the BCC-dendrites. During FSLD, the FCC-matrix underwent continuous dynamic recrystallization due to high-temperature severe plastic deformation, forming finer equiaxed grains. Simultaneously, the BCC-dendrites fractured into smaller fragments, some of which experienced partial growth and coarsening under applied stress, resulting in an hourglass morphology. In contrast, the small B2-precipitates within the BCC-fragments dissolved during the elevated temperatures of FSLD and reprecipitated as substantially finer precipitates during continuous cooling post-FSLD. Additionally, the orientation relationships between the FCC and BCC/B2 phases were completely destroyed by the severe thermoplastic deformation inherent to FSLD. The microstructural refinements led to a substantial improvement in hardness from 177 HV to 283 HV, driven by Hall-Petch strengthening. The increased number of interfaces, including coherent BCC-B2 interfaces, potentially enhances the sink strength and radiation tolerance of the HEA, making it a promising candidate for nuclear applications. In conclusion, this study also highlights FSLD as a versatile technique for achieving tunable properties in HEAs, with detailed schematics illustrating the complex mechanisms of phase transformations during processing.

Additive Manufacturing↗
DOE OSTI2022

Solving quasi-free and quadratic Lindblad master equations for open fermionic and bosonic systems

The dynamics of Markovian open quantum systems are described by Lindblad master equations. For fermionic and bosonic systems that are quasi-free, i.e. with Hamiltonians that are quadratic in the ladder operators and Lindblad operators that are linear in the ladder operators, we derive the equation of motion for the covariance matrix. This determines the evolution of Gaussian initial states and the steady states, which are also Gaussian. Using ladder super-operators (a.k.a. third quantization), here we show how the Liouvillian can be transformed to a many-body Jordan normal form which also reveals the full many-body spectrum. Extending previous work by Prosen and Seligman, we treat fermionic and bosonic systems on equal footing with Majorana operators, shorten and complete some derivations, also address the odd-parity sector for fermions, give a criterion for the existence of bosonic steady states, cover non-diagonalizable Liouvillians also for bosons, and include quadratic systems. In extension of the quasi-free open systems, quadratic open systems comprise additional Hermitian Lindblad operators that are quadratic in the ladder operators. While Gaussian states may then evolve into non-Gaussian states, the Liouvillian can still be transformed to a useful block-triangular form, and the equations of motion for k-point Green’s functions form a closed hierarchy. Based on this formalism, results on criticality and dissipative phase transitions in such models are discussed in a companion paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

In-Silico Design of Next Generation Cellulose-Derived Packaging Materials (CRADA Final Report)

Developing sustainable solutions for single-use packaging is an important objective to combat the environmental crisis of plastics pollution. Most embodiments of cellulose-based packaging materials, including CellophaneTM, are completely biodegradable in both terrestrial and marine environments. However, petroleum-derived alternatives offer some performance advantages for metrics such as moisture barriers and mechanical properties. This project leverages molecular dynamics simulation to investigate how molecular modifications to cellulose-based polymer assemblies impact their material properties. An important performance criterion for the modified materials was to retain biodegradability; thus, modifications by naturally occurring, biodegradable additives were the focus of this study. Specifically, we developed models with xylan and lignin of varying monomeric compositions into the cellulose matrix. The mechanical properties were investigated by performing stress-strain simulations, and the water barrier and hydrophobicity were investigated by simulating the water contact angle. Our findings indicate that the incorporation of xylan into the cellulose matrix tends to increase the mechanical properties with an optimal loading of ~27 wt%. We also predict that orienting the nanoscale directionality of the xylan chains such that they are perpendicular to the cellulose fibrils will dramatically increase mechanical strength. In contrast, the incorporation of lignin tends to weaken the composite at all loadings investigated. Simulations of water contact angle predicted that coating polymers on the surface of the cellulose assembly creates a more hydrophobic surface than incorporating them throughout the matrix. Of the coatings investigated, lignin resulted in the most hydrophobic surface, followed by pectin and keratin, which both imparted modest increases in hydrophobicity. Future experimental work done by Futamura will focus on designing material prototypes to capitalize on the predictions of performance enhancement obtained from molecular modeling. While substantial progress was made by the simulations performed in this project, there still exists a vast parameter space that we were unable to investigate, including branching, functional group decoration, and degree of polymerization of polymer additives. However, the methods developed in this initial investigation will facilitate more rapid evaluation of the impact of molecular characteristics on the performance of biopolymer composite materials and thereby accelerate future materials discovery efforts in this area.

36 MATERIALS SCIENCE↗
NASA NTRS1977

Dynamic crystallization experiments bearing on the origin of textures in impact-generated liquids

Dynamic crystallization experiments on a synthetic glass of 14310 composition show that the density and distribution of crystals present at the initiation of cooling can drastically affect the texture of the crystallized product. Experiments were performed in one-atmosphere gas-mixing furnaces at an oxygen fugacity slightly below iron-wustite, and precooling crystallization, cooling history, and melt history were systematically varied. When cooled from a complete liquid, the experimental charges were porphyritic with 1-2 mm phenocrysts set in a very fine-grained matrix (less than 0.1 mm). This result suggests that 14310 did not cool from a complete melt. When cooled from a crystal-liquid mixture, the experimental charges had widely varying textures depending on the density and distribution of crystals that act as nucleation sites. It is suggested that the close spatial association of widely divergent textures in the matrices of some breccias and rocks might arise in this manner; local variations in cooling rates are not required.

Lofgren, G. E.↗
DOE OSTI2025

Spectral anomalies and broken symmetries in maximally chaotic quantum maps

Spectral statistics such as the level spacing statistics and spectral form factor (SFF) are widely expected to accurately identify “ergodicity,” including the presence of underlying macroscopic symmetries, in generic quantum systems ranging from quantized chaotic maps to interacting many-body systems. By studying various quantizations of maximally chaotic maps that break a discrete classical symmetry upon quantization, we demonstrate that this approach can be misleading and fail to detect macroscopic symmetries. Notably, the same classical map can exhibit signatures of different random matrix symmetry classes in short-range spectral statistics depending on the quantization. While the long-range spectral statistics encoded in the early time ramp of the SFF are more robust and correctly identify macroscopic symmetries in several common quantizations, we also demonstrate analytically and numerically that the presence of Berry-like phases in the quantization leads to spectral anomalies, which break this correspondence. Finally, we provide numerical evidence that long-range spectral rigidity remains directly correlated with ergodicity in the quantum dynamical sense of visiting a complete orthonormal basis.

Shou, Laura [Univ. of Maryland, College Park, MD (↗
DOE OSTI2022

Quantification of electron correlation for approximate quantum calculations

State-of-the-art many-body wave function techniques rely on heuristics to achieve high accuracy at an attainable computational cost to solve the many-body Schrödinger equation. By far, the most common property used to assess accuracy has been the total energy; however, total energies do not give a complete picture of electron correlation. In this work, we assess the von Neumann entropy of the one-particle reduced density matrix (1-RDM) to compare selected configuration interaction (CI), coupled cluster, variational Monte Carlo, and fixed-node diffusion Monte Carlo for benchmark hydrogen chains. A new algorithm, the circle reject method, is presented, which improves the efficiency of evaluating the von Neumann entropy using quantum Monte Carlo by several orders of magnitude. The von Neumann entropy of the 1-RDM and the eigenvalues of the 1-RDM are shown to distinguish between the dynamic correlation introduced by the Jastrow and the static correlation introduced by determinants with large weights, confirming some of the lore in the field concerning the difference between the selected CI and Slater–Jastrow wave functions.

Chemistry↗
NASA NTRS1973

Free vibration analysis of spinning structural systems.

This article presents an efficient digital computer procedure, along with the complete listing of the associated computer program, which may be conveniently utilized for the accurate solution of a wide range of practical eigenvalue problems. Important applications of the present work are envisaged in the natural frequency analysis of spinning structures discretized by the finite element technique, and in the determination of transfer functions associated with the dynamic blocks of control systems of spacecraft utilizing gas jets or reaction wheels for attitude control, as well as of spin-stabilized and dual-spin-stabilized satellites. The validity of the Sturm sequence property is first established for the related matrix formulation involving Hermitian and real symmetric, positive-definite matrices, both being usually of highly banded configuration. A numerically stable algorithm based on the Sturm sequence method is then developed which fully exploits the banded form of the associated matrices.

Gupta, K. K.↗
NASA NTRS1986

Optimally robust redundancy relations for failure detection in uncertain systems

All failure detection methods are based, either explicitly or implicitly, on the use of redundancy, i.e. on (possibly dynamic) relations among the measured variables. The robustness of the failure detection process consequently depends to a great degree on the reliability of the redundancy relations, which in turn is affected by the inevitable presence of model uncertainties. In this paper the problem of determining redundancy relations that are optimally robust is addressed in a sense that includes several major issues of importance in practical failure detection and that provides a significant amount of intuition concerning the geometry of robust failure detection. A procedure is given involving the construction of a single matrix and its singular value decomposition for the determination of a complete sequence of redundancy relations, ordered in terms of their level of robustness. This procedure also provides the basis for comparing levels of robustness in redundancy provided by different sets of sensors.

Lou, X.-C.↗
DOE OSTI2022

UAE6 - Wind Tunnel Tests Data - UAE6 - Sequence O - Raw Data

Sequence O: Sin AOA, Parked (P) This sequence was designed to quantify the blade 3-D unsteady aerodynamic response in the absence of rotational influences by varying blade pitch angle. This test sequence used an upwind, rigid turbine with a 0° cone angle. The wind speed ranged from 18.9 m/s to 39.3 m/s, and the yaw angle was 0° throughout the sequence. The rotor was parked with the instrumented blade fixed at 0° azimuth, and the rotor lock was installed (see Appendix A). Blade and probe pressure measurements were collected. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to -99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. Sinusoidal angle-of-attack histories were designed to emulate 2-D S809 dynamic data obtained at the OSU wind tunnel (Reuss Ramsey et al. 1995) as well as 1P yawed flow angle-of-attack variations the turbine could be expected to encounter in routine operation. The reduced frequency (K), mean angle of attack (αm), and oscillation amplitude (αω) were computed for each of the five primary span locations. The file lengths varied according to the time required to complete 40 oscillations at the specified frequency. The file name convention used the sequence designation O, followed by two digits for the span location (30, 47, 80, or 95), followed by two digits representing the condition specified in the test matrix, followed by the repetition digit. Due to time constraints, the test conditions corresponding to the 63% span location were eliminated. Test points with a static pitch angle of 3° were collected throughout the sequence to ascertain instrumentation fidelity.

17 WIND ENERGY↗
NASA NTRS1983

Development of the functional simulator for the Galileo attitude and articulation control system

A simulation program for verifying and checking the performance of the Galileo Spacecraft's Attitude and Articulation Control Subsystem's (AACS) flight software is discussed. The program, which is called Functional Simulator (FUNSIM), provides a simple method of interfacing user-supplied mathematical models coded in FORTRAN which describes spacecraft dynamics, sensors, and actuators; this is done with the AACS flight software, coded in HAL/S (High-level Advanced Language/Shuttle). It is thus able to simulate the AACS flight software accurately to the HAL/S statement level in the environment of a mainframe computer system. FUNSIM also has a command and data subsystem (CDS) simulator. It is noted that the input/output data and timing are simulated with the same precision as the flight microprocessor. FUNSIM uses a variable stepsize numerical integration algorithm complete with individual error bound control on the state variable to solve the equations of motion. The program has been designed to provide both line printer and matrix dot plotting of the variables requested in the run section and to provide error diagnostics.

Namiri, M. K.↗
DOE OSTI2021

Quantum-classical hybrid algorithm for the simulation of all-electron correlation

While chemical systems containing hundreds to thousands of electrons remain beyond the reach of quantum devices, hybrid quantum-classical algorithms present a promising pathway toward a quantum advantage. Hybrid algorithms treat the exponentially scaling part of the calculation-the static correlation-on the quantum computer and the non-exponentially scaling part-the dynamic correlation-on the classical computer. While a variety of algorithms have been proposed, the dependence of many methods on the total wave function limits the development of easy-to-use classical post-processing implementations. Here, we present a novel combination of quantum and classical algorithms, which computes the all-electron energy of a strongly correlated molecular system on the classical computer from the 2-electron reduced density matrix (2-RDM) evaluated on the quantum device. Significantly, we circumvent the wave function in the all-electron calculations by using density matrix methods that only require input of the statically correlated 2-RDM. Although the algorithm is completely general, we test it with two classical density matrix methods, the anti-Hermitian contracted Schrödinger equation (ACSE) and multiconfiguration pair-density functional theories, using the recently developed quantum ACSE method for simulating the statically correlated 2-RDM. Furthermore, we obtain experimental accuracy for the relative energies of all three benzyne isomers and thereby demonstrate the ability of the developed algorithm to achieve chemically relevant and accurate results on noisy intermediate-scale quantum devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

UAE6 - Wind Tunnel Tests Data - UAE6 - Sequence N - Raw Data

Sequence N: Sin AOA, Rotating (P) This sequence was designed to quantify the blade 3-D unsteady aerodynamic response in the presence of rotational influences by varying blade pitch angle. Test sequence N used an upwind, rigid turbine with a 0° cone angle. The wind speed was 15 m/s, and the yaw angle was 0° throughout the sequence. The rotor rotated at 72 RPM. Blade and probe pressure measurements were collected. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to –99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. Sinusoidal angle-of-attack histories were designed to emulate 2-D S809 dynamic data obtained at the Ohio State University (OSU) wind tunnel (Reuss Ramsey et al. 1995) and the 1P yawed flow angle-of-attack variations the turbine could be expected to encounter in routine operation. The reduced frequency (K), mean angle of attack (αm), and oscillation amplitude (αω) were computed for each of the five primary span locations. The file lengths varied according to the time required to complete 40 oscillations at the specified frequency. The file name convention used the sequence designation N, followed by two digits for the span location (30, 47, 63, 80, or 95), followed by two digits representing the condition specified in the test matrix, followed by the repetition digit. Due to time constraints, most of the test conditions corresponding to the 63% span location were eliminated. Test points with a static pitch angle of 3° were collected throughout the sequence to ascertain instrumentation fidelity.

17 WIND ENERGY↗