Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Systematic planning of moving target defence for maximising detection effectiveness against false data injection attacks in smart grid

Abstract Moving target defence (MTD) has been gaining traction to thwart false data injection attacks against state estimation (SE) in the power grid. MTD actively perturbs the reactance of transmission lines equipped with distributed flexible AC transmission system (D‐FACTS) devices to falsify the attacker's knowledge about the system configuration. However, the existing literature has not systematically studied what influences the detection effectiveness of MTD and how it can be improved based on the topology analysis. These problems are tackled here from the perspective of an MTD plan in which the D‐FACTS placement is determined. We first exploit the relation between the rank of the composite matrix and the detecting effectiveness. Then, we rigorously derive upper and lower bounds on the attack detecting probability of MTDs with a given rank of the composite matrix. Furthermore, we analyse existing planning methods and highlight the importance of bus coverage by D‐FACTS devices. To improve the detection effectiveness, we propose a novel graph theory–based planning algorithm to retain the maximum rank of the composite matrix while covering all necessary buses. Comparative results on multiple systems show the high detecting effectiveness of the proposed algorithm in both DC‐ and AC‐SE.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Dam-Break Flow Measurements of Liquid–Granular Mixtures in a Channelized Reservoir

Dam-break flows of granular materials exhibit a complex non-Newtonian rheology, which is strongly influenced by various factors, such as the characteristics of the granular material, the degree of initial compaction, and saturation level of the granular matrix. Specially adapted nonintrusive techniques are required to investigate these flows in the laboratory. This paper investigated the nonintrusive technique of image analysis. The methodologies used in the measurement of bulk flow parameters, such as flow depth, front velocity, and flow rate, and the extraction of detailed flow characteristics such as velocity profiles, were discussed in the context of a series of granular dam-break flow experiments. These experiments were performed under three types of initial conditions: dry, saturated (water filling the pore space), and submerged (a layer of saturated granular material submerged under water) using four types of granular materials—PET pellets, crushed walnut shells, urea, and sand—with grain sizes ranging from 200 μm to 2.9 mm. High-speed cameras were used to capture the evolution and local properties of the flow, such as velocity fields using particle image velocimetry (PIV). Bulk properties of the flow were obtained by tracking the surface profile of the granular matrix and the water table. The extracted profiles were used to measure the evolution of the volume of the granular matrix and the fluid phase separately upstream and downstream of the gate, which allowed measuring the dilatancy (increase in pore space). The study showed that the image analysis techniques are versatile and can be adapted to measure different flow parameters in a variety of granular flows to better understand and analyze their complex behavior. Finally, the advantages and drawbacks of each image analysis technique were discussed and best-practice recommendations were provided.

42 ENGINEERING↗
DOE OSTI2020

Modeling of resistive plasma response in toroidal geometry using an asymptotic matching approach

The method of solving the linear resistive plasma response, based on the asymptotic matching approach, is developed for full toroidal tokamaks by upgrading the resistive DCON code [A. H. Glasser, Z. R. Wang, and J.-K. Park, Phys. Plasmas 23, 112506 (2016)]. The derived matching matrix, asymptotically matching the outer and inner regions, indicates that the applied three dimension (3-D) magnetic perturbations contribute additional small solutions at each resonant surface due to the toroidal coupling of poloidal modes. In contrast, the resonant harmonic only affects the corresponding resonant surface in the cylindrical plasma. Additionally, the solution of the ideal outer region is critical to the asymptotic matching and is challenging in toroidal geometry due to the singular power series solution at the resonant surfaces. Thus, a systematic verification of the outer region Δ' matrix is made by reproducing the well-known analytical Δ' result in Furth et al. [Phys. Fluids 16, 1054–1063 (1073)] and by making a quantitative benchmark with the PEST3 code [A. Pletzer and R. L. Dewar, J. Plasma Phys. 45, 427–451 (1991)]. Finally, the reconstructed numerical solution of the resistive plasma response from the toroidal matching matrix is presented. Compared with the ideal plasma response, the global structure of the response can be affected by the small finite island at the resonant surfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Enabling complete multichannel nonadiabatic dynamics: A global representation of the two-channel coupled, 1,2 1 A and 1 3 A states of NH 3 using neural networks

Global coupled three-state two-channel potential energy and property/interaction (dipole and spin–orbit coupling) surfaces for the dissociation of NH 3 (Ã) into NH + H 2 and NH 2 + H are reported. The permutational invariant polynomial-neural network approach is used to simultaneously fit and diabatize the electronic Hamiltonian by fitting the energies, energy gradients, and derivative couplings of the two coupled lowest-lying singlet states as well as fitting the energy and energy gradients of the lowest-lying triplet state. The key issue in fitting property matrix elements in the diabatic basis is that the diabatic surfaces must be smooth, that is, the diabatization must remove spikes in the original adiabatic property surfaces attributable to the switch of electronic wavefunctions at the conical intersection seam. Here, we employ the fit potential energy matrix to transform properties in the adiabatic representation to a quasi-diabatic representation and remove the discontinuity near the conical intersection seam. The property matrix elements can then be fit with smooth neural network functions. The coupled potential energy surfaces along with the dipole and spin–orbit coupling surfaces will enable more accurate and complete treatment of optical transitions, as well as nonadiabatic internal conversion and intersystem crossing.

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

Novel method to control explosive shock sensitivity: A mesoscale study to understand the effect of thermally expandable microsphere (TEM) inclusions in high explosives (HE) microstructure

When the void content and/or void structure of a high explosive (HE) is altered by some means (i.e., bulk heating or mechanical damage), the shock initiation behavior of the material changes. The ability to precisely predict the change in shock sensitivity after an HE has undergone microstructural changes is a crucial capability in multi-scale reactive flow models. Here, we utilize thermally expandable microspheres (TEMs) as a dopant in a polymer bonded explosive (PBX) matrix to alter the shock initiation properties in a controlled fashion. Here, using a mesoscale modeling approach, we evaluated how a single TEM (before and after thermal expansion) behaves under shock compression, as well as how the matrix PBX in the direct vicinity of the TEM is affected. We first examined the effect of an unexpanded TEM in the explosive matrix and found that its presence does not significantly perturb the bulk flow and by extension will not affect bulk sensitivity. Next, we examined the effect of an expanded TEM and found that its presence significantly perturbs the flow via hydrodynamic jetting, which causes a secondary shock wave with a strength that exceeds that of the incident wave. Finally, we showed that this secondary shock interacts with the downstream porosity to ignite a larger fraction of the overall pore volume, commensurate with the secondary shock strength and the affected volume, increasing the global (bulk) shock sensitivity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Dynamic mode decomposition with core sketch

With the increase in collected data volumes, either from experimental measurements or high fidelity simulations, there is an ever-growing need to develop computationally efficient tools to process, analyze, and interpret these datasets. Modal analysis techniques have gained great interest due to their ability to identify patterns in the data and extract valuable information about the system being considered. Dynamic mode decomposition (DMD) relies on elements of the Koopman approximation theory to compute a set of modes, each associated with a fixed oscillation frequency and a decay/growth rate. Extracting these details from large datasets can be computationally expensive due to the need to implement singular value decomposition of the input data matrix. Sketching algorithms have become popular in numerical linear algebra where statistical theoretic approaches are utilized to reduce the cost of major operations. A sketch of a matrix is another matrix, which is significantly smaller, but still sufficiently approximates the original system. We put forth an efficient DMD framework, SketchyDMD, based on a core sketching algorithm that captures information about the range and corange (their mutual relationship) of input data. The proposed sketching-based framework can accelerate various portions of the DMD routines, compared to classical methods that operate directly on the raw input data. We conduct numerical experiments using the spherical shallow water equations as a prototypical model in the context of geophysical flows. In conclusion, we show that the proposed SketchyDMD is superior to existing randomized DMD methods that are based on capturing only the range of the input data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Thermal expansion of plasma-exposed tungsten

We report results from a systematic analysis of thermal expansion of plasma-exposed tungsten based on molecular-dynamics simulations using models of tungsten with distributions of helium (He) bubbles in the tungsten matrix. We distinguish between two approaches of filling the bubbles with He, where the amount of He in the bubble can or cannot vary with temperature. In the former case, the thermal expansion coefficient decreases monotonically with the porosity and He content of the tungsten matrix, while in the latter case, the thermal expansivity increases monotonically with increasing porosity and He content. The latter condition, where the He content in the bubble is determined at the implantation temperature and remains constant with varying temperature in the tungsten matrix, is consistent with He species transport in tungsten used as a plasma-facing component (PFC) in nuclear fusion reactors and implies the development of biaxial compressive thermal strains in the PFC material that contribute to accelerating the growth of a nanostructure on PFC tungsten surfaces. Furthermore our analysis advances the fundamental understanding of thermal expansion in PFC tungsten and contributes to the development of a thermophysical property database for properly incorporating effects of realistic heat loads into modeling the dynamical response of PFC tungsten under fusion reactor operating conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

SQMBox: Interfacing a semiempirical integral library to modular ab initio electronic structure enables new semiempirical methods

Ab initio and semiempirical electronic structure methods are usually implemented in separate software packages or use entirely different code paths. As a result, it can be time-consuming to transfer an established ab initio electronic structure scheme to a semiempirical Hamiltonian. Here we present an approach to unify ab initio and semiempirical electronic structure code paths based on a separation of the wavefunction ansatz and the needed matrix representations of operators. With this separation, the Hamiltonian can refer to either an ab initio or semiempirical treatment of the resulting integrals. We built a semiempirical integral library and interfaced it to the GPU-accelerated electronic structure code TeraChem. Equivalency between ab initio and semiempirical tight-binding Hamiltonian terms is assigned according to their dependence on the one-electron density matrix. The new library provides semiempirical equivalents of the Hamiltonian matrix and gradient intermediates, corresponding to those provided by the ab initio integral library. This enables the straightforward combination of semiempirical Hamiltonians with the full pre-existing ground and excited state functionality of the ab initio electronic structure code. We demonstrate the capability of this approach by combining the extended tight-binding method GFN1-xTB with both spin-restricted ensemble-referenced Kohn–Sham and complete active space methods. We also present a highly efficient GPU implementation of the semiempirical Mulliken-approximated Fock exchange. The additional computational cost for this term becomes negligible even on consumer-grade GPUs, enabling Mulliken-approximated exchange in tight-binding methods for essentially no additional cost.

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

Generating accurate density matrices on the tangent space of a Grassmann manifold

Interpolating a density matrix from a set of known density matrices is not a trivial task. This is because a linear combination of density matrices does not necessarily correspond to another density matrix. In this Communication, density matrices are examined as objects of a Grassmann manifold. Although this manifold is not a vector space, its tangent space is a vector space. As a result, one can map the density matrices on this manifold to their corresponding vectors in the tangent space and then perform interpolations on that tangent space. The resulting interpolated vector can be mapped back to the Grassmann manifold, which can then be utilized (1) as an optimal initial guess for a self-consistent field (SCF) calculation or (2) to derive energy directly without time-consuming SCF iterations. Such a promising approach is denoted as Grassmann interpolation (G-Int). The hydrogen molecule has been used to illustrate that the described interpolated method in this work preserves the essential attributes of a density matrix. For phosphorus mononitride and ferrocene, it was demonstrated numerically that reference points for the definition of the corresponding tangent spaces can be chosen arbitrarily. In addition, the interpolated density matrices provide a superior and essentially converged initial guess for an SCF calculation to make the SCF procedure itself unnecessary. Finally, this accurate, efficient, robust, and systematically improved G-Int strategy has been used for the first time to generate highly accurate potential energy surfaces with fine details for the difficult case, ferrocene.

Chemistry↗
DOE OSTI2023

Generalized Mercier stability criterion for stellarators

The Mercier criterion is a well-known stability criterion for tokamaks. It is derived from a 2 × 2 matrix eigenvalue problem arising from the expansion of resonant solutions about a singular surface where m−nq=0, with m and n being the poloidal and toroidal mode numbers, respectively, and q being the safety factor. The stability criterion is that the eigenvalues must be real, otherwise, the solution oscillates, violating the Newcomb crossing criterion. Because of the non-axisymmetry of stellarators, different toroidal as well as poloidal harmonics couple to each other. It follows that each singular surface can have multiple resonant harmonics, with multiplicity M≥1. The corresponding matrix eigenvalue problem involves a 2M×2M matrix, resulting in M pairs of positive and negative eigenvalues. The generalized stability criterion is that all eigenvalues must be real. While the original Mercier criterion can be expressed in terms of quadratures of equilibrium quantities over the singular surface, which can be evaluated anywhere, the generalized Mercier criterion can only be evaluated on rational q surfaces with a given set of resonant harmonics.

Physics↗
DOE OSTI2023

Symmetry and scaling in one-dimensional compressible two-phase flow

Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Block2 : A comprehensive open source framework to develop and apply state-of-the-art DMRG algorithms in electronic structure and beyond

block2 is an open source framework to implement and perform density matrix renormalization group and matrix product state algorithms. Out-of-the-box it supports the eigenstate, time-dependent, response, and finite-temperature algorithms. In addition, it carries special optimizations for ab initio electronic structure Hamiltonians and implements many quantum chemistry extensions to the density matrix renormalization group, such as dynamical correlation theories. The code is designed with an emphasis on flexibility, extensibility, and efficiency and to support integration with external numerical packages. Here, we explain the design principles and currently supported features and present numerical examples in a range of applications.

Chemistry↗
DOE OSTI2024

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↗
DOE OSTI2025

Effects of a refractory metal protective barrier on uranium carbonitride ceramic–metallic composites exposed to a high temperature hydrogen environment

Uranium nitride has been under consideration at NASA and the Department of Energy for nuclear thermal propulsion applications intended for deep space exploration vehicles. It is critical to the success of the technology development that the behavior of the fuel system is properly characterized in operational environments. This body of work fabricated ceramic metal composites consisting of uranium ceramic fuel in a molybdenum and molybdenum–tungsten metal matrix, both with and without an exterior molybdenum barrier. Specimens were introduced to a high temperature hydrogen environment and characterized through x-ray diffraction, scanning electron microscopy, and energy dispersive spectroscopy to identify characteristic changes. A portion of the specimens was fabricated with a 1.27 mm thick molybdenum barrier in order to isolate the material from the hydrogen environment. Instability of the uranium nitride occurred as low as 2073 K, leading to aggressive instability approaching 2658 K. Carbon reactions with the metal matrix were observed as well as alloying between the uranium and metal matrix constituents. The liquid uranium resulting from the instability of the uranium nitride was observed to attack the grain boundaries of the refractory metal, resulting in the formation of a liquid uranium molybdenum eutectic melt. Specimens encased in a barrier experienced similar results. The authors conclude this is a significant challenge to applications of this cermet at these temperatures. Alternate configurations with resilient hermetic seals or lower temperature applications such as power reactors may have more success.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidation of the physical factors that control activated transport of penetrants in chemically complex glass-forming liquids

Understanding the activated transport of penetrant or tracer atoms and molecules in condensed phases is a challenging problem in chemistry, materials science, physics, and biophysics. Many angstrom- and nanometer-scale features enter due to the highly variable shape, size, interaction, and conformational flexibility of the penetrant and matrix species, leading to a dramatic diversity of penetrant dynamics. Based on a minimalist model of a spherical penetrant in equilibrated dense matrices of hard spheres, a recent microscopic theory that relates hopping transport to local structure has predicted a novel correlation between penetrant diffusivity and the matrix thermodynamic dimensionless compressibility, S 0 ( T ) (which also quantifies the amplitude of long wavelength density fluctuations), as a consequence of a fundamental statistical mechanical relationship between structure and thermodynamics. Moreover, the penetrant activation barrier is predicted to have a factorized/multiplicative form, scaling as the product of an inverse power law of S 0 ( T ) and a linear/logarithmic function of the penetrant-to-matrix size ratio. This implies an enormous reduction in chemical complexity that is verified based solely on experimental data for diverse classes of chemically complex penetrants dissolved in molecular and polymeric liquids over a wide range of temperatures down to the kinetic glass transition. The predicted corollary that the penetrant diffusion constant decreases exponentially with inverse temperature raised to an exponent determined solely by how S 0 ( T ) decreases with cooling is also verified experimentally. Our findings are relevant to fundamental questions in glassy dynamics, self-averaging of angstrom-scale chemical features, and applications such as membrane separations, barrier coatings, drug delivery, and self-healing.

Science & Technology - Other Topics↗
DOE OSTI2020

Methodology for Generating Covariance Data of Thermal Neutron Scattering Cross Sections

This paper details and implements a framework for evaluating thermal neutron scattering cross sections that provide S(α,β) data and covariance data for hydrogen in light water. This methodology involves perturbing model parameters of molecular dynamics potentials and fitting the simulation results to experimental data. The framework is general and can be applied to any material or simulation method. The fit is made using the Unified Monte Carlo method to experimentally measure double-differential scattering cross sections of light water at the Spallation Neutron Source at Oak Ridge National Laboratory. Mean values and covariance data were generated for model parameters, phonon density of states, double-differential cross sections, and total scattering cross sections. These posterior parameter values were very similar to their prior values with a maximum relative error of 0.54%. This falls within in the Unified Monte Carlo–calculated uncertainties on the order of 2.7%. Additionally, posterior double-differential cross sections agree favorably with ENDF/B-VIII.0 cross sections. The new thermal scattering law was tested by comparing it against benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook, which showed a slight improvement over the ENDF/B-VIII.0 library. Additionally, the covariance matrix of the phonon density of states was validated to confirm that the spread of k eff from the density of states used to generate the covariance matrix was similar to the spread of k eff from the density of states of the sampled covariance matrix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Stress Profile in Coating Layers of TRISO Fuel Particles in Contact with One Another

This work presents a discussion on a series of finite element analyses that assess stress evolution in the coating layers of tristructural isotropic (TRISO) particles in contact with each other while embedded in a matrix. The initial simulations were of applied uniaxial pressure versus matrix elastic modulus. These simulations predicted increasing stress in the silicon carbide coating layers of the TRISO particles with decreasing matrix elastic modulus. The second set of simulations focused on the effects of heating and cooling and the associated dimensional change on the state of stress in the coating layers. The general finding was that there was no significant difference below the coating layer’s deposition temperature. Although, above the deposition temperature, the contacting particles had higher stress compared with those that were separated. The third set of simulations focused on the effects of irradiation, specifically, creep, dimensional change, and swelling. An interface debonding model was introduced since these potential effects have a significant bearing on predicted stresses.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Toward the first gluon parton distribution from the LaMET

Abstract We present progress towards the first unpolarized gluon quasi-parton distribution function (PDF) from lattice quantum chromodynamics using high-statistics measurements for hadrons at two valence pion masses M π ≈ 310 and 690 MeV computed on an a ≈ 0.12 fm ensemble with 2 + 1 + 1-flavors of highly improved staggered quark generated by the MILC collaboration. In this study, we consider two gluon operators for which the hybrid-ratio renormalization matching kernels have been recently derived and a third operator that has been used in prior pseudo-PDF studies of the gluon PDFs. We compare the matrix elements for each operator for both the nucleon and pion, at both pion masses, and using two gauge-smearing techniques. Focusing on the more phenomenologically studied nucleon gluon PDF, we compare the ratio and hybrid-ratio renormalized matrix elements at both pion masses and both smearings to those reconstructed from the nucleon gluon PDF from the CT18 global analysis. We identify the best choice of operator to study the gluon PDF and present the first gluon quasi-PDF under some caveats. Additionally, we explore the recent idea of Coulomb gauge fixing to improve signal at large Wilson-line displacement and find it could be a major help in improving the signal in the gluon matrix elements. This work helps identify the best operator for studying the gluon quasi-PDF, shows higher hadron boost momentum is needed to implement hybrid-ratio renormalization reliably, and suggests the need to study more diverse set of operators with their corresponding perturbative calculations for hybrid-ratio renormalization to further gluon quasi-PDF study.

Good, William (ORCID:0000000184281003)↗