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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Understanding surfaces and interfaces in nanocomposites of silicone and barium titanate through experiments and modeling

Barium titanate (BTO) is a ferroelectric perovskite used in electronics and energy storage systems because of its high dielectric constant. Decreasing the BTO particle size was shown to increase the dielectric constant of the perovskite, which is an intriguing but contested result. We investigated this result by fabricating silicone-matrix nanocomposite specimens containing BTO particles of decreasing diameter. Furthermore, density functional theory modeling was used to understand the interactions at the BTO particle surface. Combining results from experiments and modeling indicated that polymer type, particle surface interactions, and particle surface structure can influence the dielectric properties of polymer-matrix nanocomposites containing BTO.

composite↗
DOE OSTI2023

Superconducting valence bond fluid in lightly doped eight-leg $t-J$ cylinders

Superconductivity in doped quantum paramagnets has been a subject of long theoretical inquiry. In this work, we report a density matrix renormalization group study of lightly doped $t-J$ models on the finite-width square lattice (doped hole densities $δ$ = 1/12 and 1/8) with parameters for which previous studies have suggested that the undoped system in 2D is either a quantum spin liquid or a valence bond crystal. Our studies are performed on cylinders with width up to 8. Here, ground-state correlations are found to be nearly identical for the “doped quantum spin liquid” and “doped valence bond crystal.” Upon increasing the cylinder widths from 4 to 8, we observed a significant strengthening of the quasi-long-range superconducting correlations and a dramatic suppression of any “competing” charge density wave order. Extrapolating from the observed behavior of the width eight cylinders, we speculate that the system has a nodeless $d$-wave superconducting ground state in the 2D limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Directly deposited graphene nanowalls on carbon fiber for improving the interface strength in composites

Graphene nanowalls (GNWs) were grown directly on carbon fibers using a chemical vapor deposition technique which is simple and catalyst-free. We found that there is very strong π-π stacking which is a benefit for the GNWs/carbon fiber interface. This single modified filament then was embedded into an epoxy matrix to be a single-fiber composite in which was formed a “tenon-mortise” structure. Such a “tenon-mortise” model provides a simple, stable, and powerful connection between carbon fiber and the epoxy matrix. In addition, it was demonstrated that the epoxy matrix can be well embedded into GNWs through a field emission scanning electron microscope. The results of the single-fiber composite tests indicated that the interfacial strength of the composites was immensely improved by 173% compared to those specimens without GNWs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Revisiting matrix-based inversion of scanning mobility particle sizer (SMPS) and humidified tandem differential mobility analyzer (HTDMA) data

Abstract. Tikhonov regularization is a tool for reducing noise amplification during data inversion. This work introduces RegularizationTools.jl, a general-purpose software package for applying Tikhonov regularization to data. The package implements well-established numerical algorithms and is suitable for systems of up to ∼ 1000 equations. Included is an abstraction to systematically categorize specific inversion configurations and their associated hyperparameters. A generic interface translates arbitrary linear forward models defined by a computer function into the corresponding design matrix. This obviates the need to explicitly write out and discretize the Fredholm integral equation, thus facilitating fast prototyping of new regularization schemes associated with measurement techniques. Example applications include the inversion involving data from scanning mobility particle sizers (SMPSs) and humidified tandem differential mobility analyzers (HTDMAs). Inversion of SMPS size distributions reported in this work builds upon the freely available software DifferentialMobilityAnalyzers.jl. The speed of inversion is improved by a factor of ∼ 200, now requiring between 2 and 5 ms per SMPS scan when using 120 size bins. Previously reported occasional failure to converge to a valid solution is reduced by switching from the L-curve method to generalized cross-validation as the metric to search for the optimal regularization parameter. Higher-order inversions resulting in smooth, denoised reconstructions of size distributions are now included in DifferentialMobilityAnalyzers.jl. This work also demonstrates that an SMPS-style matrix-based inversion can be applied to find the growth factor frequency distribution from raw HTDMA data while also accounting for multiply charged particles. The outcome of the aerosol-related inversion methods is showcased by inverting multi-week SMPS and HTDMA datasets from ground-based observations, including SMPS data obtained at Bodega Marine Laboratory during the CalWater 2/ACAPEX campaign and co-located SMPS and HTDMA data collected at the US Department of Energy observatory located at the Southern Great Plains site in Oklahoma, USA. Results show that the proposed approaches are suitable for unsupervised, nonparametric inversion of large-scale datasets as well as inversion in real time during data acquisition on low-cost reduced-instruction-set architectures used in single-board computers. The included software implementation of Tikhonov regularization is freely available, general, and domain-independent and thus can be applied to many other inverse problems arising in atmospheric measurement techniques and beyond.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Superconductivity in a topological lattice model with strong repulsion

The highly tunable nature of synthetic quantum materials-both in the solid-state and cold atom contexts-invites examining which microscopic ingredients aid in the realization of correlated phases of matter such as superconductors. Recent experimental advances in moiré materials suggest that unifying the features of the Fermi-Hubbard model and quantum Hall systems creates a fertile ground for the emergence of such phases. Here, we introduce the "double Hofstadter"model, a minimal 2D lattice model that incorporates exactly these features: Time-reversal symmetry, band topology, and strong repulsive interactions. By using infinite cylinder density matrix renormalization group methods (cylinder iDMRG), we investigate the ground state phase diagram of this model. We find that it hosts an interaction-induced quantum spin Hall insulator and demonstrate that weakly hole-doping this state gives rise to a superconductor at a finite circumference, with indications that this behavior persists on larger cylinders. Further, at the aforementioned circumference, the superconducting phase is surprisingly robust to perturbations including additional repulsive interactions in the pairing channel. By developing a technique to probe the superconducting gap function in iDMRG, we phenomenologically characterize the superconductor. Namely, we demonstrate that it is formed from the weak pairing of holes atop the quantum spin Hall insulator. Furthermore, we determine the pairing symmetry of the superconductor, finding it to be p-wave-reminiscent of the unconventional superconductivity reported in experiments on twisted bilayer graphene (TBG). Motivated by this, we elucidate structural similarities and differences between our model and those of TBG in its chiral limit. Finally, to provide a more direct experimental realization, we detail an implementation of our Hamiltonian in a system of cold fermionic alkaline-earth atoms in an optical lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

A patterned phase-changing vanadium dioxide film stacking with VO 2 nanoparticle matrix for high performance energy-efficient smart window applications

A vanadium dioxide (VO 2 ) based solid-to-solid phased changing material has been attracting great interest in smart window applications. However, achieving high solar modulation and high transparency simultaneously in visible light is the major challenge for the practical application of this smart material. To resolve this issue, in this paper, a smart film composed of a VO 2 nanoparticle matrix and a patterned VO 2 film is presented. Furthermore, numerical modeling and electromagnetic simulation are carried out to characterize the performance in terms of solar modulation and luminous transmittance, and a parametric study is carried out to optimize the proposed smart window film. Compared with the VO 2 nanoparticle matrix, the proposed structure can obtain 23% solar modulation and 57% luminous transmittance but with a much thinner thickness, which will significantly reduce the cost and fabrication complexity and extend the environment stability.

42 ENGINEERING↗
DOE OSTI2013

Understanding the solidification and microstructure evolution during CSC-MIG welding of Fe–Cr–B-based alloy

The present is a study of the solidification and microstructure of Fe–28.2%Cr–3.8%B–1.5%Si–1.5%Mn (wt.%) alloy deposited onto a 1020 plain carbon steel substrate using the controlled short-circuit metal inert gas welding process. The as-solidified alloy was a metal matrix composite with a hypereutectic microstructure. Thermodynamic calculation based on the Scheil–Gulliver model showed that a primary (Cr,Fe){sub 2}B phase formed first during solidification, followed by an eutectic formation of the (Cr,Fe){sub 2}B phase and a body-centered cubic Fe-based solid solution matrix, which contained Cr, Mn and Si. Microstructure analysis confirmed the formation of these phases and showed that the shape of the (Cr,Fe){sub 2}B phase was irregular plate. As the welding heat input increased, the weld dilution increased and thus the volume fraction of the (Cr,Fe){sub 2}B plates decreased while other microstructural characteristics were similar. - Highlights: • We deposit Fe–Cr–B-based alloy onto plain carbon steel using the CSC-MIG process. • We model the solidification behavior using thermodynamic calculation. • As deposited alloy consists of (Cr,Fe){sub 2}B plates embedded in Fe-based matrix. • We study the effect of the welding heat input on the microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Exciton Level Structure and Dynamics in Tubular Porphyrin Aggregates

We present an account of the optical properties of the Frenkel excitons in self-assembled porphyrin tubular aggregates that represent an analog to natural photosynthetic antennae. Using a combination of ultrafast optical spectroscopy and stochastic exciton modeling, we address both linear and nonlinear exciton absorption, relaxation pathways, and the role of disorder. The static disorder-dominated absorption and fluorescence line widths show little temperature dependence for the lowest excitons (Q band), which we successfully simulate using a model of exciton scattering on acoustic phonons in the host matrix. Temperature-dependent transient absorption of and fluorescence from the excitons in the tubular aggregates are marked by nonexponential decays with time scales ranging from a few picoseconds to a few nanoseconds, reflecting complex relaxation mechanisms. Combined experimental and theoretical investigations indicate that nonradiative pathways induced by traps and defects dominate the relaxation of excitons in the tubular aggregates. We model the pumpprobe spectra and ascribe the excited-state absorption to transitions from one-exciton states to a manifold of mixed one- and two-exciton states. Our results demonstrate that while the delocalized Frenkel excitons (over 208 (1036) molecules for the optically dominant excitons in the Q (B) band) resulting from strong intermolecular coupling in these aggregates could potentially facilitate efficient energy transfer, fast relaxation due to defects and disorder probably present a major limitation for exciton transport over large distances.

biomimetic photosynthetic antennae↗
DOE OSTI2017

Modeling Dynamic Helium Release as a Tracer of Rock Deformation

Here, we use helium released during mechanical deformation of shales as a signal to explore the effects of deformation and failure on material transport properties. A dynamic dual-permeability model with evolving pore and fracture networks is used to simulate gases released from shale during deformation and failure. Changes in material properties required to reproduce experimentally observed gas signals are explored. We model two different experiments of 4He flow rate measured from shale undergoing mechanical deformation, a core parallel to bedding and a core perpendicular to bedding. We also found that the helium signal is sensitive to fracture development and evolution as well as changes in the matrix transport properties. We constrain the timing and effective fracture aperture, as well as the increase in matrix porosity and permeability. Increases in matrix permeability are required to explain gas flow prior to macroscopic failure, and the short-term gas flow postfailure. Increased matrix porosity is required to match the long-term, postfailure gas flow. This model provides the first quantitative interpretation of helium release as a result of mechanical deformation. The sensitivity of this model to changes in the fracture network, as well as to matrix properties during deformation, indicates that helium release can be used as a quantitative tool to evaluate the state of stress and strain in earth materials.

58 GEOSCIENCES↗
DOE OSTI2020

Quantitative and three-dimensional assessment of holdup material

Nuclear material deposited in equipment, transfer lines, and ventilation systems of a processing facility is usually referred to as holdup. Operators need to know the location and amount of holdup in their facility for radiation safety and material accountability. In this work, we propose to use an array of detectors co-axial to the inspected pipe to measure the holdup material. This method is implementable into an automated system capable of crawling on surfaces and pipes of various curvatures, which would enable faster, easier, and more accurate holdup safeguards measurements. We first demonstrated that the current holdup assay procedure could lead to a non-negligible bias in the estimate of special nuclear material mass, due to the simplified assumption of deposited geometry introduced by the Generalized Geometry Holdup (GGH) model. The new approach consists of imaging the inner holdup material by characterizing the detector array’s response and unfolding it from the measured light output. Our experimental proof of principle consists of three NaI(Tl) detectors surrounding an aluminum pipe containing two cesium-137 ( 137 Cs) sources. We derived the source distribution inside the pipe by first calculating the detector response matrix using a method adaptive to the surface geometry of the object containing the measured holdup material. Creating a matrix of the detector array’s measured counts, we then proceed to solve an inverse problem, resulting in an accurately located source position and activity distribution within the response matrix’s spatial resolution. Furthermore, we then developed a simulated model of the envisioned experimental setup, which accurately described both the activity and position of the source in 2D. Finally, we extended our model onto a discretized three-dimensional model of the system, encompassing 36 detectors. The model was simulated to measure the source distribution and activity of multiple sources in pipes of varied geometries, which could then be conveniently translated to a real-time holdup measurement protocol. For the 3D simulation of four different source geometries, the model accurately localized the source position in 3D, while the activity retained a maximum relative error of ±5.32%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Exciton dissociation at donor-acceptor heterojunctions: Dynamics using the collective effective mode representation of the spin-boson model

Following the recent quantum dynamics investigation of the charge transfer at an oligothiophene-fullerene heterojunction by the multi-configuration time dependent Hartree method [H. Tamura, R. Martinazzo, M. Ruckenbauer and I. Burghardt, J. Chem. Phys. 137, 22A540 (2012)], we revisit the transfer process by a perturbative non-Markovian master equation treated by the time local auxiliary density matrix approach. We compare the efficiency of the spin-boson model calibrated by quantum chemistry with the effective mode representation. A collective mode is extracted from the spin-boson spectral density. It is weakly coupled to a residual bath of vibrational modes, allowing second-order dynamics. The electron transfer is analyzed for a sampling of inter-fragment distances showing the fine interplay of the electronic coupling and energy gap on the relaxation. The electronic coherence, expected to play a role in the process, is preserved during about 200 fs.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Permittivity effects of particle agglomeration in ferroelectric ceramic-epoxy composites using finite element modeling

Here, the size dependence of the dielectric constants of barium titanate or other ferroelectric particles can be explored by embedding particles into an epoxy matrix whose dielectric constant can be measured directly. However, to extract the particle dielectric constant requires a model of the composite medium. We compare a finite element model for various volume fractions and particle arrangements to several effective medium approximations, which do not consider particle arrangement explicitly. For a fixed number of particles, the composite dielectric constant increases with the degree of agglomeration, and we relate this increase to the number of regions of enhanced electric field along the applied field between particles in an agglomerate. Additionally, even for dispersed particles, we find that the composite method of assessing the particle dielectric constant may not be effective if the particle dielectric constant is too high compared to the background medium dielectric constant.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A numerically exact description of ultrafast vibrational decoherence in vibration-coupled electron transfer

Broadband pump–probe spectroscopy has been widely used to measure vibrational decoherence associated with the reaction coordinate in photoinduced ultrafast vibration-coupled electron transfer (VCET) reactions. These experiments provide insight into the interplay of intramolecular coordinates along the reaction coordinate. However, a general theoretical foundation for analyzing, and even for explaining rigorously, these data is lacking. In this work, we study vibrational decoherence in a model VCET reaction using the nearly exact time-dependent density matrix renormalization group simulation method. We explore how analyzing the density matrix with quantum information measures can help elucidate the evolution of vibrational coherence in simulations of dynamics. We examine how vibrational coherence is affected by electron transfer on the timescale of approximately 100 femtoseconds. Our results suggest that electron transfer, in the nonadiabatic model, changes the vibrational equilibrium position abruptly—an example of a “quantum quench” event. This explains the concomitant vibrational decoherence. We find that abrupt vibrational decoherence can be mitigated by wavepacket motion occurring on the timescale of the electron transfer.

Science & Technology - Other Topics↗
DOE OSTI2016

Symmetric quadratic Hamiltonians with pseudo-Hermitian matrix representation

We prove that any symmetric Hamiltonian that is a quadratic function of the coordinates and momenta has a pseudo-Hermitian adjoint or regular matrix representation. The eigenvalues of the latter matrix are the natural frequencies of the Hamiltonian operator. When all the eigenvalues of the matrix are real, then the spectrum of the symmetric Hamiltonian is real and the operator is Hermitian. As illustrative examples we choose the quadratic Hamiltonians that model a pair of coupled resonators with balanced gain and loss, the electromagnetic self-force on an oscillating charged particle and an active LRC circuit. -- Highlights: •Symmetric quadratic operators are useful models for many physical applications. •Any such operator exhibits a pseudo-Hermitian matrix representation. •Its eigenvalues are the natural frequencies of the Hamiltonian operator. •The eigenvalues may be real or complex and describe a phase transition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Fuel Performance Analysis of Fast Flux Test Facility MFF-3 and -5 Fuel Pins Using BISON with Post Irradiation Examination Data

Using the BISON fuel-performance code, simulations were conducted of an automated process to read initial and operating conditions from the Pacific Northwest National Laboratory (PNNL) database and reports, which contain metallic-fuel data from the Fast Flux Test Facility (FFTF) MFF Experiments. This work builds on previous modeling efforts involving 1977 EBR-II metallic fuel pins from experiments. Coupling the FFTF PNNL reports to BISON allowed for all 338 pins from MFF-3 and MFF-5 campaigns to be simulated. Each BISON simulation contains unique power and flux histories, axial power and flux profiles, and coolant-channel flow rates. Fission-gas release (FGR), fuel axial swelling, cladding profilometry, and burnup were all simulated in BISON and compared to available post-irradiation examination (PIE) data. Cladding profilometry, FGR, and fuel axial swelling simulation results for full-length MFF metallic pins were found to be in agreement with PIE measurements using FFTF physics and models used previously for EBR-II simulations. The main two peaks observed within the cladding profilometry were able to be simulated, with fuel-cladding mechanical interaction (FCMI), fuel-cladding chemical interaction (FCCI), and thermal and irradiation-induced creep being the cause. A U-Pu-Zr hot-pressing model was included in this work to allow pore collapse within the fuel matrix. This allowed better agreement between BISON-simulated cladding profilometry and PIE measurements for the peak caused by FCMI. This work shows that metallic fuel models used to accurately represent fuel performance for smaller EBR-II pins may be used for full-length metallic fuel, such as FFTF MFF assemblies and the Versatile Test Reactor (VTR). As new material models and PIE measurements become available, FFTF MFF assessment cases will be reassessed to further BISON model development.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Modeling the sequential dissociative double ionization of O 2 by ultrashort intense infrared laser pulses

A density matrix approach for sequential double ionization (DM-SDI) of molecules has been developed recently and was applied to the N 2 molecule. In this article, we extended the DM-SDI model to O 2 , which is a more complicated system to model than N 2 , due to its electronic structures and spin-orbit and laser couplings in the manifold of doubly charged states. We obtained a good agreement on the kinetic energy release spectrum of O + + O + from previous experiments. Thanks to the low computational cost of the model, we explored the mechanism behind the ionization and dissociation dynamics as well as the effects of lasers on the spectrum. Furthermore, this work will pave the way to model sequential dissociative double ionization of larger molecules and to probe molecular dynamics by measuring kinetic energy release spectra from this process.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

The polarized-signal density matrix: A practical way to recover molecular frame information from isotropic samples

We present a novel approach to model ultrafast time-dependent nonlinear optical polarization sensitive signals emitted from randomly oriented molecules. By projecting the laboratory-frame analyzer polarization axis into the molecular-frame and linking that axis with the density matrix through a tensor product, we demonstrate an approach to find a specific molecular orientation that yields a good approximation to simulated four-wave mixing signals produced by the same model but with averaging over molecular orientation.

Thurston, Richard L↗
DOE OSTI2026

Stochastic fracture generation and thermo-hydro-mechanical modeling in an equivalent continuum framework for enhanced geothermal systems

Enhanced geothermal systems (EGS) involve fracturing low permeability material to establish well connectivity and then injecting and circulating fluid into the fractured subsurface for geothermal power production. Changes in fracture aperture from contraction of the cooling matrix rock may alter network connectivity and risk thermal short-circuiting. Thermo-hydro-mechanical (THM) models are a useful tool to study these processes. However, as fracture networks are complex, and data may be limited, fracture networks in THM models are often stochastically generated. Given reliance on stochastic fracture networks and THM modeling to represent the subsurface and assess productivity of EGS, increased understanding of the influence of such statistically derived fracture networks on flow and heat transport in THM models is needed. Here, a new fracture process model is developed in the reactive transport code PFLOTRAN to stochastically generate fracture families and simulate changes in fracture aperture over time due to temperature changes of the rock matrix. Sixty-four different fracture networks ranging from well to poorly-connected, are modeled in PFLOTRAN with and without mechanical processes (THM vs TH). Results indicate that for well-connected fracture networks, thermal short-circuiting is less of a concern due to the abundance of available alternative flowpaths. For poorly-connected fracture networks, inclusion of mechanical processes showed steep thermal drawdown coincident with increase in fracture aperture along developing colder flowpaths, demonstrating the risk of thermal short-circuiting. Simulations with additional, larger fractures engineered to establish connectivity in a poorly-fractured subsurface, indicate that while stochastic variation of fracture orientation of the background network had limited influence, such variation in the engineered fractures significantly affected flow and heat transport.

Discrete fracture networks (DFN)↗