Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix models”

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 847 records · Page 47

Status on lattice calculations of the proton spin decomposition

Abstract Lattice calculations of the proton spin components is reviewed. The lattice results of the quark spin from the axial-vector current matrix element at ∼ 0.3−0.4 is smaller than those from the constituent quark models. This is largely due to the fact that the vacuum polarization contribution from the disconnected insertion is negative. Its connection with the anomalous Ward identity is clarified and verified numerically. This resolves the contentious issue in the “proton spin crisis.” The glue spin and angular momentum are found to be large and there is notable contribution from the quark orbital angular momentum. Renormalization, mixing, and normalization of the quark and glue angular momenta are discussed. With sufficient precision, they can be compared with more precise experimental measurements when the electron-ion collider facility is available.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Extending Orbital-Optimized Density Functional Theory to L-Edge XPS and Beyond: Spin–Orbit Coupling via Nonorthogonal Quasi-Degenerate Perturbation Theory

Quantum mechanical calculations of core electron binding energies (CEBEs) are relevant to interpreting X-ray photoelectron spectroscopy (XPS). Orbital-optimized density functional theory (OO-DFT) accurately predicts K-edge CEBEs but is challenged by the presence of significant spin–orbit coupling (SOC) at L- and higher edges involving inner-shell orbitals with nonzero angular momentum. Here, to extend OO-DFT to L-edges and higher, our method utilizes scalar-relativistic, spin-restricted open-shell OO-DFT to construct a minimal, quasi-degenerate basis of core-hole states corresponding to a chosen inner-shell (e.g., ionizing all six possible 2p spin orbitals). Nonorthogonal configuration interaction (NOCI) is then used to obtain the matrix elements of the full Hamiltonian including SOC in this quasi-degenerate model space of determinants. Using a screened 1-electron SOC operator parametrized with the Dirac-Coulomb-Breit (DCB) Hamiltonian results in doublet splitting (DS) values for third row elements that are nearly in quantitative agreement with experiment. The resulting NOCI eigenvalues are shifted by the average of the (scalar) OO-DFT CEBEs to yield CEBEs (split by SOC) corrected for dynamic correlation. Comparing calculations on gas phase molecules with experimental results establishes that NO-QDPT with the SCAN functional (NO-QDPT/SCAN), using the DCB screened 1-electron SOC operator is accurate to about 0.2 eV for L-edge CEBEs of molecules containing third row atoms. However, this NO-QDPT approach becomes less accurate for fourth-row elements starting in the middle of the 3d transition metal series, with errors increasing as atomic number increases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Complex Fluid‐Driven Fractures Caused by Crack‐Parallel Stress

Abstract Managing fluid‐driven fracture networks is crucial for subsurface resource utilization, yet the current understanding of the key controlling factors remains insufficient. While geologic discontinuities have been shown to significantly influence fracture network complexity, this study identifies another major contributor. We conducted a new set of experiments using a transparent true triaxial cell, which enabled video recording of the temporal evolution of fluid‐driven fracture paths. Using pseudo‐2D samples without macroscale structural discontinuities, we observed multiple occurrences of hydraulic fracture curving and branching under anisotropic boundary stresses. We proposed a theoretical model demonstrating that the stress parallel to the crack line in the solid matrix near the crack tip (i.e., the T ‐stress) accounts for the observed fracture curving behavior. This finding suggests that T ‐stress is an additional mechanism contributing to the complexity of fluid‐driven fracture networks in the subsurface, besides the geologic discontinuities.

58 GEOSCIENCES↗

Compact spectroscopy of keV to MeV X-rays from a laser wakefield accelerator

We reconstruct spectra of secondary X-rays from a tunable 250–350 MeV laser wakefield electron accelerator from single-shot X-ray depth-energy measurements in a compact (7.5 × 7.5 × 15 cm), modular X-ray calorimeter made of alternating layers of absorbing materials and imaging plates. X-rays range from few-keV betatron to few-MeV inverse Compton to > 100 MeV bremsstrahlung emission, and are characterized both individually and in mixtures. Geant4 simulations of energy deposition of single-energy X-rays in the stack generate an energy-vs-depth response matrix for a given stack configuration. An iterative reconstruction algorithm based on analytic models of betatron, inverse Compton and bremsstrahlung photon energy distributions then unfolds X-ray spectra, typically within a minute. We discuss uncertainties, limitations and extensions of both measurement and reconstruction methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bond order predictions using deep neural networks

Machine learning is an extremely powerful tool for the modern theoretical chemist since it provides a method for bypassing costly algorithms for solving the Schrödinger equation. Already, it has proven able to infer molecular and atomic properties such as charges, enthalpies, dipoles, excited state energies, and others. Most of these machine learning algorithms proceed by inferring properties of individual atoms, even breaking down total molecular energy into individual atomic contributions. In this paper, we introduce a modified version of the Hierarchically Interacting Particle Neural Network (HIP-NN) capable of making predictions on the bonds between atoms rather than on the atoms themselves. We train the modified HIP-NN to infer bond orders for a large number of small organic molecules as computed via the Natural Bond Orbital package. We demonstrate that the trained model is extensible to molecules much larger than those in the training set by studying its performance on the COMP6 dataset. This method has applications in cheminformatics and force field parameterization and opens a promising future for machine learning models to predict other quantities that are defined between atoms such as density matrix elements, Hamiltonian parameters, and molecular reactivities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Gapless spin liquid and pair density wave of the Hubbard model on three-leg triangular cylinders

We study the ground state properties of the Hubbard model on three-leg triangular cylinders using large-scale density-matrix renormalization group simulations. At half-filling, we identify an intermediate gapless spin liquid phase, which has one gapless spin mode and algebraic spin–spin correlations but exponential decay scalar chiral–chiral correlations, between a metallic phase at weak coupling and Mott insulating dimer phase at strong interaction. Upon light doping the gapless spin liquid, the system exhibits power-law charge-density-wave (CDW) correlations but short-range single-particle, spin–spin, and chiral–chiral correlations. Similar to CDW correlations, the superconducting correlations also decay in power-law but oscillate in sign as a function of distance, which is consistent with the striped pair-density wave. When further doping the gapless spin liquid phase or doping the dimer order phase, another phase takes over, which has similar CDW correlations but all other correlations decay exponentially.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Baryon acoustic oscillations in thin redshift shells from BOSS DR12 and eBOSS DR16 galaxies

ABSTRACT In an age of large astronomical data sets and severe cosmological tensions, the case for model independent analyses is compelling. We present a set of 14 baryon acoustic oscillations measurements in thin redshift shells with $3\,\mathrm{ per} \,\mathrm{ cent}$ precision that were obtained by analysing BOSS DR12 and eBOSS DR16 galaxies in the redshift range 0.32 < z < 0.66. Thanks to the use of thin shells, the analysis is carried out using just redshifts and angles so that the fiducial model is only introduced when considering the mock catalogues, necessary for the covariance matrix estimation and the pipeline validation. We compare our measurements, with and without supernova data, to the corresponding constraints from Planck 2018, finding good compatibility. A Monte Python module for this likelihood is available at github.com/ranier137/angularBAO.

79 ASTRONOMY AND ASTROPHYSICS↗

Niobium hydride formation in superconducting qubit thin films

The formation of nonsuperconducting hydrides in 160–170-nm-thick films of niobium is examined. We identify six elastically distinct orientation relationships ɛ−Nb 4 ⁢H 3 takes within the matrix, solid-solution 𝛼−NbH 𝑥 phase. We employ a phase field model to assess the impact of elastic energy induced by the strain of phase transformation on the morphology and transformation dynamics of ɛ precipitates within a thin film. We consider the dimensions of the thin film, crystallographic growth direction, and diffusion rates to predict the timescale of hydride evolution. Leveraging the finite element method, we predict two-dimensional and three-dimensional equilibrium shapes of ɛ-hydrides within a bulk sample and in a thin film that has a traction-free surface. Our results suggest that niobium hydrides migrate to the free surface of the film. Precipitates which reach the free surface coarsen, while precipitates within the film dissolve. Precipitates in both two dimensions and three dimensions experience a repulsive interaction force at the free surface, that is attractive in the bulk, shown in experiment and theory of previous studies.

36 - MATERIALS SCIENCE↗

FeederGAN: Synthetic Feeder Generation via Deep Graph Adversarial Nets

This paper presents a novel, automated, generative adversarial networks (GAN) based synthetic feeder generation mechanism, abbreviated as FeederGAN. FeederGAN digests real feeder models represented by directed graphs via a deep learning framework powered by GAN and graph convolutional networks (GCN). Information of a distribution feeder circuit is extracted from its model input files so that the device connectivity is mapped onto the adjacency matrix and the device characteristics, such as circuit types (i.e., 3-phase, 2-phase, and 1-phase) and component attributes (e.g., length and current ratings), are mapped onto the attribute matrix. Then, Wasserstein distance is used to optimize the GAN and GCN is used to discriminate the generated graphs from the actual ones. A greedy method based on graph theory is developed to reconstruct the feeder using the generated adjacency and attribute matrices. Furthermore, our results show that the GAN generated feeders resemble the actual feeder in both topology and attributes verified by visual inspection and by empirical statistics obtained from actual distribution feeders.

42 ENGINEERING↗

Optimized Carbon Fiber Intermediate Development to Enable High-Volume Manufacturing of Lightweight Automotive Composites

The ongoing pursuit of improved fuel economy and reduced greenhouse gas emissions has resulted in sustained interest for lightweight materials technologies. In this context, carbon fiber composites have captured the imagination of automotive engineers due to the potential to achieve substantial mass reduction when compared to traditional steel construction. That stated, the use of carbon fiber composites in automotive has been limited for the most part to premium supercars and other derivative platforms. In these cases, manufacturing costs are less of an obstacle to implementation, and the performance benefits of carbon fiber have enabled production of structures offering more than 50% weight savings. In practice, translating these low-volume demonstrations onto high-volume vehicle platforms has remained challenging. This can be attributed to several factors, with the absence of suitable high throughput production methodologies being a key impediment. To date, structural, crash critical components have relied upon manufacturing techniques born out of the aerospace industry. This has created a disconnect between automotive production systems that are accustomed to manufacturing multiple parts per minute and the aerospace technologies that have cycle times in the order of hours. Consequently, the focus of this project is the development of manufacturing process technology for carbon fiber composites that can support a mainstream vehicle program at an assumed throughput of 100,000 vehicles per year. In practice, this translates to a part-to-part cycle time of less than 3 minutes. Project participation included contributions from a broad range of academic, industrial and national lab partners. The primary scope of work, being the development of new carbon fiber epoxy compounds that are stable at room temperature and suited to high throughput automated processing. For project management, the work streams were divided into six key areas, with the lead organization in parentheses. • Carbon fiber/epoxy materials formulation development and scale up (Dow) • Simulation of discontinuous near isotropic meso-structure intermediates (Purdue) • Simulation of mechanical performance of compression molded components (Purdue) • Meso-Scale morphological analysis and correlation with structural performance (UTK) • Paint and adhesion durability analysis (MSU) • Demonstrator part design, prototype production, and validation testing (Ford). The primary goal at the commencement of the project was development of a chopped carbon fiber sheet molding compound (SMC) that offered a three times improvement in tensile modulus over a comparable glass-based SMC. In addition to meeting mechanical performance targets, the resin kinetics were modified to achieve a processing cycle time of less than 3 minutes. Other critical-to-quality (CTQ) specifications were also stipulated to account for a broad range of materials and processing characteristics. To achieve the above, staff scientists at Dow Chemical created an extensive series of new epoxy blends for testing and validation. Throughout this development, a key challenge was attaining material performance goals without impacting processing behavior and paintability of finished components. The latter required a new internal mold release system being developed by Dow that was designed to complement the kinetics of the rapid cure epoxy. As a complement to work studies at the industrial partners, the teams from academia executed a series of analytical and experimental studies to investigate potential factors influencing CF-SMC performance. Unit cell models were developed to capture the meso-scale representations of the fiber matrix architecture. Results of this analysis and subsequent morphological investigations led to the design of a novel composite derivative comprising carbon fiber platelets embedded in an epoxy matrix; the platelet size and aspect ratio playing significant role in final composite properties. This approach was a departure from previous research in CF-SMC development whereby bulk filamentization or disassembly of the carbon fiber rovings had been considered the most effective means of achieving both fiber wet through and wet-out. As the course of the academia studies progressed, the aspect ratio of the fiber constituents was further optimized before finalizing material attributes and processing conditions. For the purposes of technology validation, the Ford team led a work stream devoted to the design, fabrication and testing of demonstration components. The carbon fiber SMC material has the potential to displace numerous stampings and castings on an automotive structure but ultimately vehicle closure applications were selected to showcase the abilities of the CF-SMC to achieve both mass reduction and business case for large complex structures. Using target properties established by the Dow staff scientists, the complete closure system for a full-size sedan decklid and a mid-size wagon liftgate were engineered. Prototypes for both applications were fabricated using production representative processing methods to allow for physical testing and performance validation of the CF-SMC structures. Following completion of a testing program that concluded with a FMVSS301 55 mph offset rear crash, the CF-SMC formulation was declared by the Ford team to have met all engineering requirements. To summarize, the joint development activities during this project led to significant technical breakthroughs and achievement of all milestones. The result was the development of a novel, tack-free carbon fiber molding compound that is suited to automated processing. This combined room temperature stability, fast cure kinetics, and internal mold release system facilitates cycle times that are conducive to high-volume production. The VORAFUSE M6400 successfully passed technology validation at Ford and is now eligible for consideration on future production commercial vehicle programs.

36 MATERIALS SCIENCE↗

Transport Simulations on Scanning Transmission Electron Microscope Images of Nanoporous Shale

Digital rock physics is an often-mentioned approach to better understand and model transport processes occurring in tight nanoporous media including the organic and inorganic matrix of shale. Workflows integrating nanometer-scale image data and pore-scale simulations are relatively undeveloped, however. In this paper, a workflow is demonstrated progressing from sample acquisition and preparation, to image acquisition by Scanning Transmission Electron Microscopy (STEM) tomography, to volumetric reconstruction to pore-space discretization to numerical simulation of pore-scale transport. Key aspects of the workflow include (i) STEM tomography in high angle annular dark field (HAADF) mode to image three-dimensional pore networks in µm-sized samples with nanometer resolution and (ii) lattice Boltzmann method (LBM) simulations to describe gas flow in slip, transitional, and Knudsen diffusion regimes. It is shown that STEM tomography with nanoscale resolution yields excellent representation of the size and connectivity of organic nanopore networks. In turn, pore-scale simulation on such networks contributes to understanding of transport and storage properties of nanoporous shale. Interestingly, flow occurs primarily along pore networks with pore dimensions on the order of tens of nanometers. Smaller pores do not form percolating pathways in the sample volume imaged. Apparent gas permeability in the range of 10−19 to 10−16 m2 is computed.

42 ENGINEERING↗

A stable potential-based time-domain method for wideband elec- tromagnetic analysis

In previous research, the frequency-domain A-ϕ formulation has been validated using the finite element method for electromagnetic simulations of low-frequency and multi- scale problems, demonstrating excellent numerical accuracy, good matrix condition, and high computational efficiency. Time- domain simulations provide significant advantages for modeling wideband problems and are crucial for multiphysics applications. In this paper, the frequency-domain A-ϕ formulation is extended to the time domain. The central difference scheme is employed for temporal discretization to ensure both accuracy and stability. A numerical example is presented to demonstrate the capability of the proposed time-domain method in wideband electromagnetic analysis.

Mekonnen, Minyichil↗

Integrable multistate Landau–Zener models with parallel energy levels

In this paper, we discuss solvable multistate Landau–Zener (MLZ) models whose Hamiltonians have commuting partner operators with ~1/τ-time-dependent parameters. Many already known solvable MLZ models belong precisely to this class. We derive the integrability conditions on the parameters of such commuting operators, and demonstrate how to use such conditions in order to derive new solvable cases. We show that MLZ models from this class must contain bands of parallel diabatic energy levels. The structure of the scattering matrix and other properties are found to be the same as in the previously discussed completely solvable MLZ Hamiltonians.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Characterization of shale using Helium and Argon at high pressures

In order to estimate the shale gas in place and the eventual recovery during shale gas operations, one of the key requirements is to accurately characterize the shale’s petrophysical and transport properties such as porosity, permeability, diffusivity, and storage capacity. Despite the many efforts reported in the technical literature aiming to provide an improved understanding of the complex pore structures and the associated fluid flow in gas shales, complete characterization of organic-rich shale samples still poses a big challenge. Here, we have characterized mass transfer and sorption in shale at different length scales using Helium (He) and Argon (Ar) as probe gases. Thermogravimetric analysis (TGA) with a shale cube of ~1 cm 3 in size and gas expansion experiments with a full-diameter core (3.5” in diameter) were used to measure sorption kinetics/isotherms and mass transfer, respectively. Both samples are from the same depth/location in the Marcellus shale formation. The TGA steady-state technique was initially used to generate excess sorption isotherms for Ar, while dynamic TGA experiments were used to study its sorption kinetics. The TGA experiments demonstrate that Ar, which has a similar sorption potential as Methane, but is generally assumed to be inert, adsorbs onto the surfaces of the mesoporous and microporous regions of the shale samples according to a Langmuir-type behavior. Helium expansion experiments, on the full-diameter core, were used to measure the overall porosity, on the basis that He is a non-sorbing and inert gas as compared to Ar. The He expansion experiments, furthermore, allow us to delineate the mass transfer of gas across the inherent hierarchy of pore sizes, including macropores (macro- and microcracks), mesopores and micropores. Similar expansion experiments were also performed with Ar to study the combined impact of mass transfer and sorption. A triple-porosity model (TPM) was utilized to interpret the He expansion experiments with the shale core and to extract (estimate) relevant transport parameters. We report and compare here the diffusivities and permeabilities of the whole core for both He and Ar, as calculated from the modeling and fitting of the experimental data. On the premise that the shale cube is representative of the matrix region of the core, the Ar sorption kinetics from the cube experiments were subsequently combined with the transport parameters extracted from the He experiments to predict the behavior of the Ar expansion test with the full-diameter core. An excellent agreement is observed between the model predictions and the experimental data. The experimental observations and their interpretation indicate that one must be cautious when using Ar to estimate the true porosity and permeability of shales. In addition, we demonstrate that He and Ar probe gases, when used in tandem, can be employed effectively as a tool to characterize shales in terms of mass transfer and sorption dynamics across scales.

04 OIL SHALES AND TAR SANDS↗

Characterization of elastic mechanical properties of Tuscaloosa Marine Shale from well logs using the vertical transversely isotropic model

To avoid steep declines in the Tuscaloosa Marine Shale (TMS) production, wells are fracture-stimulated to release the hydrocarbons trapped in the matrix of the formation. An accurate estimation of Young’s modulus and Poisson’s ratio is essential for hydraulic fracture propagation. In addition, ignoring the highly heterogeneous and anisotropic character of TMS can lead to erroneous stress values, which subsequently affect hydraulic fracture width estimates and the overall hydraulic fracturing process. We have developed an empirical 1D geomechanical model that takes into account VTI anisotropy, and it is used to characterize the elastic mechanical properties of TMS in two wells. In the analyzed formation, the vertical Poisson’s ratio is less than the horizontal Poisson’s ratio, which suggests the necessity of an alternative to the ANNIE equations. The stiffness coefficients [Formula: see text] and [Formula: see text] were estimated using the relationships developed from the ultrasonic core data available for the two TMS. Further, correlations between the static and dynamic properties from laboratory tests were used to improve the minimum horizontal stress calculation. We compare VTI Young’s moduli, Poisson’s ratios, and minimum horizontal stress with the isotropic solution. VTI modeling improves the estimation of the elastic mechanical properties. The isotropic solution underestimates the minimum horizontal stress in the formation. Moreover, it was shown that the 20 ft shale interval below the TMS base is characterized by a low Young’s modulus (the vertical Young’s modulus is equal to 20 GPa, whereas the horizontal Young’s modulus is equal to 40 GPa) and may be a frac barrier.

Geochemistry & Geophysics↗

Magnon interactions in the quantum paramagnetic phase of CoNb 2 O 6

In this work, we study effects of magnon interactions in the excitation spectrum of CoNb 2 O 6 in the quantum paramagnetic phase in transverse field, where the 1/S spin-wave theory exhibits unphysical divergences at the critical field. We propose a self-consistent Hartree-Fock approach that eliminates such unphysical singularities while preserving the integrity of the singular threshold phenomena of magnon decay and spectrum renormalization that are present in both theory and experiment. With the microscopic parameters adopted from previous studies, this method yields a close quantitative agreement with the available experimental data for CoNb 2 O 6 in the relevant regime. Furthermore, insights into the general structure of the spin-anisotropic model of CoNb 2 O 6 and related zigzag chain materials are also provided and a discussion of the effects of additional longitudinal field on the spectrum is given.

1-dimensional spin chains↗

Intrusive Uncertainty Quantification and Optimal Experiment Design in the Open-Source Pyomo Ecosystem

This contribution describes ParmEst and Pyomo.DoE, two pillars of the open-source Python-based Pyomo ecosystem for computational optimization with (partial differential) algebraic equation mathematical models. Specifically, ParmEst facilitates intrusive frequentist parameter estimation (PE) and uncertainty quantification (UQ) through built-in features, such as covariance matrix estimation, bootstrapping, and likelihood ratio tests. Complementary, Pyomo.DoE enables optimal experiment design by maximizing various metrics of the Fisher information matrix, such as A-optimality (trace), D-optimality (determinant), E-optimality (minimum eigenvalue), and ME-optimality (condition number). ParmEst and Pyomo.DoE can solve high-dimensional optimization problems by leveraging the model structure and exact derivative information. Finally, we will discuss future opportunities to integrate PE and UQ capabilities with optimization under uncertainty, including robust optimization with non-convex models via PyROS.

97 MATHEMATICS AND COMPUTING↗

Quantum Monte Carlo calculation of 𝛿 NS in 10 C using an effective field theory approach

Here, we compute radiative corrections to the superallowed 𝛽 decay of 10 C in an effective field theory approach using nuclear matrix elements obtained from quantum Monte Carlo calculations. These corrections are an important ingredient in the extraction of the Cabibbo-Kobayashi-Masakawa quark mixing matrix element 𝑉 𝑢⁢𝑑 , and the role of this work is to illuminate the uncertainties arising from nuclear structure. Our results provide good agreement with both the traditional extraction of 𝑉 𝑢⁢𝑑 , as well as with a more recent evaluation performed using the no-core shell model and a dispersion formalism. The dominant uncertainty in this approach is the presence of two unknown low-energy constants that enter into the relevant nuclear matrix elements. Future determinations of these low-energy constants—either from QCD or modeling them with two nucleon amplitudes—would improve the precision of the extraction in this formalism.

beta decay↗