Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix approximation”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Generalized spin σ -SCF method

We introduce a generalization of the σ-SCF method to approximate noncollinear spin ground and excited single-reference electronic states by minimizing the Hamiltonian variance. The new method is based on the σ-SCF method, originally proposed by Ye et al. [J. Chem. Phys. 147, 214104 (2017)], and provides a prescription to determine ground and excited noncollinear spin states on an equal footing. Our implementation was carried out utilizing an initial simulated annealing stage followed by a mean-field iterative self-consistent approach to simplify the cumbersome search introduced by generalizing the spin degrees of freedom. The simulated annealing stage ensures a broad exploration of the Hilbert space spanned by the generalized spin single-reference states with random complex element-wise rotations of the generalized density matrix elements in the simulated annealing stage. The mean-field iterative self-consistent stage employs an effective Fockian derived from the variance, which is utilized to converge tightly to the solutions. This process helps us to easily find complex spin structures, avoiding manipulating the initial guess. As proof-of-concept tests, we present results for Hn (n = 3–7) planar rings and polyhedral clusters with geometrical spin frustration. We show that most of these systems have noncollinear spin excited states that can be interpreted in terms of geometric spin frustration. These states are not directly targeted by energy minimization methods, which are meant to converge to the ground state. This stresses the capability of the σ-SCF methodology to find approximate noncollinear spin structures as mean-field excited states.

Chemistry↗
DOE OSTI2022

Twist-3 partonic distributions from lattice QCD

Twist-3 partonic distributions contain important information that characterizes nucleon’s structure. In this work, we show our lattice exploration of the twist-3 PDFs g_T(x) g T ( x ) , and h_L(x) h L ( x ) . We also present our preliminary results on the twist-3 GPD \tilde{G}_2(x) G ̃ 2 ( x ) . We use the quasi-distribution approach to connect the lattice-extracted matrix elements, renormalized in the RI/MOM scheme, to light-cone distributions, applying the matching procedure that we developed in parallel. We also calculate the twist-2 counterparts of g_T(x) g T ( x ) and h_L(x) h L ( x ) , i.e. g_1(x) g 1 ( x ) , and h_1(x) h 1 ( x ) , and test the Wandzura-Wilczek approximation.

Bhattacharya, Shohini↗
DOE OSTI2020

Structural Insights into Low and High Recalcitrance Natural Poplar Variants Using Neutron and X-ray Scattering

Here, we used small-angle neutron scattering (SANS) and wide-angle X-ray scattering (WAXS) to study two Populus trichocarpa variants, GW-11012 and BESC-316, with lignin contents of 17.8 and 23.2%, respectively, to understand why the low lignin variant has increased sugar release during liquid hot water (LHW) pretreatment. SANS showed differences in cell wall polymer arrangement in the variants. In GW-11012, the cellulose microfibrils were similar in size to BESC-316 but were not equally spaced and lignin was aggregated rather than homogenously interspersed with other matrix co-polymers. However, after LHW pretreatments, the SANS profiles of the pretreated variants were almost indistinguishable, indicating that the nanoscale structural features could not adequately explain the increased sugar release observed in GW-11012. WAXS analysis revealed higher cellulose crystallinity and larger crystallite size in native and pretreated GW-11012. However, the most significant difference was that the crystallite orientation was approximately 2-fold higher in GW-11012. This suggests that higher crystallite orientation may be favorable for the action of processive cellulases leading to increased sugar release from GW-11012. Our study shows that there is a complex interplay between cell wall polymers during their deposition resulting in changes in lignin distribution and cellulose structure, highlighting the importance of investigating plant cell wall architecture across multiple length scales.

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

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

Seasonal Disorder in Urban Traffic Patterns: A Low Rank Analysis

This article proposes several advances to sparse nonnegative matrix factorization (SNMF) as a way to identify large-scale patterns in urban traffic data. The input to our model is traffic counts organized by time and location. Nonnegative matrix factorization additively decomposes this information, organized as a matrix, into a linear sum of temporal signatures. Penalty terms encourage this factorization to concentrate on only a few temporal signatures, with weights which are not too large. Our interest here is to quantify and compare the regularity of traffic behavior, particularly across different broad temporal windows. In addition to the rank and error, we adapt a measure introduced by Hoyer to quantify sparsity in the representation. Combining these, we construct several curves which quantify error as a function of rank (the number of possible signatures) and sparsity; as rank goes up and sparsity goes down, the approximation can be better and the error should decreases. Plots of several such curves corresponding to different time windows leads to a way to compare disorder/order at different time scalewindows. In this paper, we apply our algorithms and procedures to study a taxi traffic dataset from New York City. In this dataset, we find weekly periodicity in the signatures, which allows us an extra framework for identifying outliers as significant deviations from weekly medians. We then apply our seasonal disorder analysis to the New York City traffic data and seasonal (spring, summer, winter, fall) time windows. We do find seasonal differences in traffic order.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

The importance of maldistribution matching for thermal performance of compact heat exchangers

Compact heat exchangers have gained increased attention in recent years, particularly in demanding applications where high temperatures, high pressures, and/or high power densities are required. For decades, the heat exchanger (HX) community believes that flow maldistribution is a key factor for HX effectiveness, that is, reducing the degree of flow maldistribution (MALD) can help increase the HX effectiveness. Therefore, significant efforts have been devoted in the past to optimizing the header geometry to minimize flow maldistribution. This work was initially motivated by this, and the original goal was to figure out a HX header design with the lowest maldistribution. However, by systematically constructing a comprehensive maldistribution matrix, the analysis revealed that the HX effectiveness is not actually determined by the MALD, but instead dominated by the degree of maldistribution mismatch (MISM). This conclusion was also theoretically generalized, which indicated that matching of the local heat capacity rate is key for achieving maximum performance. The MISM provides a local means of tracking this information, while the MALD only provides a global approximation of the maldistribution itself. With this new perspective, flow maldistribution needs not necessarily be avoided, but instead matched between two fluid streams, to improve the HX performance. We demonstrated that by carefully designing the header geometry to match the velocity profiles of the two fluids in a 2 MW PCHE with molten salt and supercritical carbon dioxide (sCO2) as the heat transfer fluids, the HX could achieve a higher effectiveness even when the maldistribution increased. Finally, a technoeconomic study using a CSP system as an example revealed that the use of this new HX design paradigm could result in CSP capital cost savings as large as 16.6%.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Real-space solution to the electronic structure problem for nearly a million electrons

We report a Kohn–Sham density functional theory calculation of a system with more than 200 000 atoms and 800 000 electrons using a real-space high-order finite-difference method to investigate the electronic structure of large spherical silicon nanoclusters. Our system of choice was a 20 nm large spherical nanocluster with 202 617 silicon atoms and 13 836 hydrogen atoms used to passivate the dangling surface bonds. To speed up the convergence of the eigenspace, we utilized Chebyshev-filtered subspace iteration, and for sparse matrix–vector multiplications, we used blockwise Hilbert space-filling curves, implemented in the PARSEC code. For this calculation, we also replaced our orthonormalization + Rayleigh–Ritz step with a generalized eigenvalue problem step. We utilized all of the 8192 nodes (458 752 processors) on the Frontera machine at the Texas Advanced Computing Center. We achieved two Chebyshev-filtered subspace iterations, yielding a good approximation of the electronic density of states. Our work pushes the limits on the capabilities of the current electronic structure solvers to nearly 106 electrons and demonstrates the potential of the real-space approach to efficiently parallelize large calculations on modern high-performance computing platforms.

Chemistry↗
DOE OSTI2024

Development of creep resistant ferritic-martensitic steels for operation at 650°C

Creep resistant ferritic-martensitic steels capable of operating at 650°C have long been under development for electricity generation. This presentation will showcase the design and development of alloys CPJ-7 and JMP with approximately 420,000 h of in-house creep testing, the equivalent of ~50 years of cumulative creep tests. Approximately fifty different alloy compositions were manufactured as part of an ICME approach. CPJ-7 present superior creep properties compared to COST alloys for steam turbine and P91/92 for boiler applications. The prolonged creep life was attributed to slowing down the destabilization of fine precipitates at 650°C. The JMP steels were designed with higher Co for solid solution strengthening and increased W for matrix strength and stability. The creep life compared to CPJ-7 increased between 118 and 150% at 650°C depending on the stress. On a Larson-Miller plot, the performance of the JMP steels surpasses that of MARBN steel. Composition-microstructure-properties relationships will be discussed.

Detrois, Martin↗
DOE OSTI2020

Thermal shock resistance of multilayer silicon carbide receiver tubes for 800°C molten salt concentrating solar power application

CSP power tower receiver systems during rapid transient weather periods can be vulnerable to thermal shock conditions from rain that which can facilitate the onset of leaks and failures that can have catastrophic consequences. Silicon carbide (SiC) materials have attractive receiver application characteristics for being light weight, having high-strength and excellent thermal shock resistance performance which make them a particularly good fit for receiver absorber materials in CSP. In this investigation, the performance characteristics of Ceramic Tubular Products (CTP) SiC ceramic matrix composite (CMC), multilayered tubes were explored with respect to thermal shock performance for solar receiver applications in next generation CSP plants. Here, thermal shock testing was performed at the Sandia National Laboratories (SNL) Solar Furnace facility using a dynamic stage and thermal shock tube test setup. The tubes tested under incident solar heat flux of 100 W/cm2 were heated with inner tube temperatures reaching approximately 800 °C, with outer temperatures exceeding or just reaching 1000 El for the multilayer and monolithic SiC tubes respectively. The tubes were then quenched with simulated rain. The tubes were then cooled and subjected to hoop stress analysis using an Instron device to assess their subsequent mechanical strength. Furthermore, the on-sun study experimental results indicate an average of 24.2% and 97% higher hoop strength for the CMC tubes than those composed of monolithic SiC and aluminum oxide (A1203) respectively.

14 SOLAR ENERGY↗
DOE OSTI2020

Stochastically Realized Observables for Excitonic Molecular Aggregates

In this work, we show that a stochastic approach enables calculations of the optical properties of large 2-dimensional and nanotubular excitonic molecular aggregates. Previous studies of such systems relied on numerically diagonalizing the dense and disordered Frenkel Hamiltonian, which scales approximately as $\mathcal{O}$($N$ 3 ) for $N$ dye molecules. Our approach scales much more efficiently as $\mathcal{O}$($N$ log($N$)), enabling quick study of systems with a million of coupled molecules on the micrometer size scale. We calculate several important experimental observables, including the optical absorption spectrum and density of states, and develop a stochastic formalism for the participation ratio. Quantitative agreement with traditional matrix diagonalization methods is demonstrated for both small- and intermediate-size systems. The stochastic methodology enables the study of the effects of spatial-correlation in site energies on the optical signatures of large 2D aggregates. Our results demonstrate that stochastic methods present a path forward for screening structural parameters and validating experiments and theoretical predictions in large excitonic aggregates.

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

Gravitational wave modes in matter

Here, a general linear gauge-invariant equation for dispersive gravitational waves (GWs) propagating in matter is derived. This equation describes, on the same footing, both the usual tensor modes and the gravitational modes strongly coupled with matter. It is shown that the effect of matter on the former is comparable to diffraction and therefore negligible within the geometrical-optics approximation. However, this approximation is applicable to modes strongly coupled with matter due to their large refractive index. GWs in ideal gas are studied using the kinetic average-Lagrangian approach and the gravitational polarizability of matter that we have introduced earlier. In particular, we show that this formulation subsumes the kinetic Jeans instability as a collective GW mode with a peculiar polarization, which is derived from the dispersion matrix rather than assumed a priori. This forms a foundation for systematically extending GW theory to GW interactions with plasmas, where symmetry considerations alone are insufficient to predict the wave polarization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Advancing attenuation estimation through integration of the Hessian in multiparameter viscoacoustic full-waveform inversion

Accurate seismic attenuation models of subsurface structures not only enhance subsequent migration processes by improving fidelity, resolution, and facilitating amplitude-compliant angle gather generation but also provide valuable constraints on subsurface physical properties. Leveraging full-wavefield information, multiparameter viscoacoustic full-waveform inversion ( Q-FWI) simultaneously estimates seismic velocity and attenuation ( Q) models. However, a major challenge in Q-FWI is the contamination of crosstalk artifacts, where inaccuracies in the velocity model are mistakenly mapped to the inverted attenuation model. While incorporating the Hessian is expected to mitigate these artifacts, the explicit implementation is prohibitively expensive due to its formidable computational cost. In this study, we formulate and develop a Q-FWI algorithm via the Newton-conjugate gradient (CG) framework, where the search direction at each iteration is determined through an internal CG loop. In particular, the Hessian is integrated into each CG step in a matrix-free fashion using the second-order adjoint-state method. We find through synthetic experiments that our Newton-CG Q-FWI significantly mitigates crosstalk artifacts compared with the limited-memory Broyden-Fletcher-Goldfarb-Shanno method and the CG method, albeit with a notable computational cost. In the discussion of several key implementation details, we also determine the significance of the approximate Gauss-Newton Hessian, the second-order adjoint-state method, and the two-stage inversion strategy.

Geochemistry & Geophysics↗
DOE OSTI2026

Simulation and Experimental Validation of an Integrated Heat Pump – Thermal Energy Storage Using a Room-Temperature Phase Change Material

As the dependence on electrical heat pumps (HPs) and intermittent renewables increases, grid strains are expected to grow. This necessitates an energy storage system to reduce the mismatch between energy supply and demand. Thus, a proposed dual-mode commercially available 14.1 kW HP was integrated with a single 22°C phase change material (PCM) thermal storage system (TES) to load-shift both cooling and heating loads. The HP-TES system was manufactured and experimentally tested using a novel test matrix based on AHRI 210/240 psychrometric conditions. Furthermore, transient dual-mode system-level HP-TES models were developed in Modelica and validated using the experimental test conditions. Base HP cooling and heating experimental tests at ambient temperatures of 35°C and −8.3°C show that the modified HP-TES maintained the rated system capacity and performance. The HP-TES discharge provided approximately 30% and 50% reductions in cooling and heating demand, respectively. The transient HP-TES models predicted system capacity and total power input for discharge and recharge operating modes within ±4% mean percentage error, and recharge power input within ±2%, with maximum errors occurring at the equipment startup. During system operation, the sources of model deviations are first-order polynomial fits of the PCM digital scanning calorimetry (DSC) data and unaccounted supercooling in the PCM during solidification. Nonetheless, the model predictions agree with the experimental tests, demonstrating the availability of robust, accurate, and validated transient models that can be used for further validation and the development of system controls.

25 ENERGY STORAGE↗
DOE OSTI2023

Development of Novel Ferritic-Martensitic Steels with Superior Creep Properties for Power Plant Applications

A creep resistant martensitic steel, CPJ-7, was developed with an operating temperature approaching 650°C. Subsequently, another creep resistant martensitic steel, JMP, was designed with the potential to operate at, or slightly above, 650°C. This report describes the development of these alloys from early iterations of CPJ steels to the final CPJ-7 formulation as well as the more heat resistant JMP steel formulation. The design originated from computational modeling for phase stability and precipitate strengthening using fifteen constituent elements. Approximately forty heats of CPJ and ten of JMP, each weighing ~7 kg, were vacuum induction melted. A computationally optimized heat treatment schedule was developed to homogenize the ingots prior to hot forging and rolling prior to final normalization and tempering. Overall, wrought and cast versions of CPJ-7 present superior creep properties when compared to wrought and cast versions of COST alloys for steam turbine and wrought and cast versions of P91/92 for boiler applications. For instance, the Larson Miller Parameter curve for CPJ-7 at 650°C almost coincides with that of COST E at 620°C. The prolonged creep life was attributed to slowing down the process of the destabilization of the MX and M 23 C 6 precipitates at 650°C. On the other hand, the cast version of CPJ-7 also revealed superior mechanical performance, well above commercially available cast 9% Cr martensitic steel or derivatives, especially for creep life at 650°C. The casting process employed slow cooling to simulate the conditions of a thick wall full-size steam turbine casing but utilized a separate homogenization step prior to final normalization and tempering. To advance the development of CPJ-7 for commercial applications, a process was used to scale up the production of the alloy using vacuum induction melting (VIM) and electroslag remelting (ESR), which underlined the importance of melt processing control of intentionally designed minor and trace elements in these advanced alloys. Following the work on CPJ-7, the JMP steels were designed with higher Co for increased solid solution strengthening, Si for oxidation resistance and increased W (with low Mo content) for matrix strength and stability as well as solid solution strengthening. The JMP steels showed increases in creep life compared to CPJ-7 between 118 to 150% at 650°C for testing at various stresses between 138 MPa and 207 MPa. On a Larson-Miller plot, the performance of the JMP steels surpasses that of state-of-the-art MARBN and other MARBN-type steels. The influence of various elements within the composition of the alloys on the microstructure and mechanical properties are discussed. This report presents approximately 420,000 h of in-house creep testing, the equivalent of almost 50 years of cumulative creep tests.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Eigenmode analysis of the sheared-flow Z-pinch

Experiments have demonstrated that a Z-pinch can persist for thousands of times longer than the growth time of global magnetohydrodynamic (MHD) instabilities such as the m=0 sausage and m=1 kink modes. These modes have growth times on the order of ta=a/vi, where vi is the ion thermal speed and a is the pinch radius. Axial flows with duz/dr ≲ vi/a have been measured during the stable period, and the commonly accepted theory is that this amount of shear is sufficient to stabilize these modes as predicted by numerical studies using the ideal MHD equations. However, these studies only consider specific equilibrium profiles that typically have a modest magnitude for the logarithmic pressure gradient, qP≡d ln P/d ln r, and may not represent experimental conditions. Linear stability of the sheared-flow Z-pinch is studied here via a direct eigen-decomposition of the matrix operator obtained from the linear ideal MHD equations. Several equilibrium profiles with a large variation of qP are examined. Considering a practical range of k, 1/3 ≲ ka ≲ 10, it is shown that the shear required to stabilize m=0 modes can be expressed as duz/dr≥Cγ0/(ka)α. Here, γ0=γ0(ka) is the profile-specific growth rate in the absence of shear, which scales approximately with |qP|. Both C and α are profile-specific constants, but C is order unity and α≈1. It is further demonstrated that even a large value of shear, duz/dr=3vi/a, is not sufficient to provide linear stabilization of the m=1 kink mode for all profiles considered. This result is in contrast to the currently accepted theory predicting stabilization at much lower shear, duz/dr=0.1vi/a, and suggests that the experimentally observed stability cannot be explained within the linear ideal-MHD model.

Angus, J. R. (ORCID:0000000314740002)↗
DOE OSTI2021

Comparing different density-matrix expansions for long-range pion exchange

Empirical energy density functionals (EDFs) are generally successful in describing nuclear properties across the table of nuclides. But their limitations motivate using the density-matrix expansion (DME) to embed long-range pion interactions into a Skyrme functional. Recent results on the impact of the pion were both encouraging and puzzling, necessitating a careful re-examination of the DME implementation. In this study, we take the first steps, focusing on two-body scalar terms in the DME. Exchange energies with long-range one-pion contributions are well approximated by all DME implementations considered, with preference for variants that do not truncate at two derivatives in every EDF term. The use of the DME for chiral pion contributions is therefore supported by this investigation. For scalar-isovector energies it is important to treat neutrons and protons separately. The results are found to apply under broad conditions, although self-consistency is not yet tested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Dynamics of Polymerization and Gelation in Epoxy Nanocomposites via X-ray Photon Correlation Spectroscopy

The details of the curing process of epoxy resins are notoriously difficult to ascertain, hampering polymer matrix composite manufacturing. Here, the curing of a series of six epoxy resins containing dilute spherical nanoparticles was investigated via X-ray photon correlation spectroscopy (XPCS). The resin formulation was varied to achieve a range of crosslink densities and topologies. Prior to gelation, the viscosities obtained from XPCS agreed with bulk rheology. A clear dynamic transition was detected in all samples, where the q-scaling changes from approximately q2 .2 to q 1 , indicating a transition from sub-diffusive to ballistic-like nanoparticle motion. For the samples with a high crosslink density, this change was abrupt and coincided with gelation, and there was no evidence of crosslink heterogeneity. Chain entanglement played a role in the transition for samples with lower crosslink density. Finally, this work reveals the potential of XPCS to deliver new quantitative insights into the curing process of thermosetting resins.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗