Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density 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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Effect of heat treatment on the tensile behavior of selective laser melted Ti-6Al-4V by in situ X-ray characterization

Selective Laser Melted Ti-6Al-4V (as-SLMed) exhibits decreased yield strength, increased work hardening, and increased ductility after heat treatment at 730 degrees C (HT-730) or 900 degrees C (HT-900) for 2 h. To understand the change of mechanical properties, in situ high energy X-ray diffraction (HEXRD) is used to examine the phase composition, load partitioning, slip system activity, and dislocation density evolution in all three specimens. The as-SLMed specimen is dominated by martensitic α'. After heat treatment, α' partly or fully decomposes into α+β, reducing the yield strength. In HT-730, beta precipitates with confined size show much higher lattice strain than the α'/α matrix during deformation; in HT-900, the lattice strain difference is mostly eliminated. This is a key reason for the increased ductility in HT-900. From the anisotropic lattice strain development, basal slip is identified as the easiest slip system in α'/α. Using an elasto-plastic self-consistent (EPSC) model, the critical resolved shear stress ratio between prismatic slip and basal slip (CRSS prism /CRSS basal ) is estimated to be 1.31 and 1.16 in the as -SLMed and the HT-900 specimens, respectively. Further, the alpha phase in HT-900 is able to activate multiple slip systems and accumulate more dislocations during plastic deformation. This explains why HT-900 has better ductility and higher work hardening rate than the other two specimens.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exploring the speciation of actinide salts in the presence of contaminants related to molten salt reactors

Molten salt reactors (MSR) are up-and-coming Generation IV nuclear reactors with either coolant and/or radioactive fuel in molten salt form. Due to their improved safety, efficiency, affordability and convenient waste processing system, these reactors are considered a superior alternative to conventional ones. However, prior to commercial implementation, it is necessary to have an extensive understanding of the reactions occurring in the MSR. The safety parameters of the MSR are determined by the rheological properties such as density and viscosity of the molten salt and these properties are governed by the local structure of the molten salts. Due to the high operational temperatures of MSR's, corrosion plays an important role. Because of their influence, the container corrosion products become part of the fuel salt matrix. In this work, we are investigating how these corrosion products, e.g., Mo, Ni, Cr, would impact the speciation of the actinide metal centers. For the initial studies, lanthanides are used as a surrogate for actinides. Lanthanide chloride salts (LnCl3) are mixed with transition metals and excess of alkali or alkali earth metal salts or salt mixtures (LiCl, NaCl, KCl). These are added to an alumina crucible, placed in a quartz tube, and sealed under a vacuum. This is heated up to around 1000 °C and slowly cooled to a temperature above the melting point of the salt and at that temperature, the reaction is taken out of the furnace. The excess salt is removed, and the resulting product is analyzed using X-ray diffraction techniques, scanning electron microscopy, and spectroscopic methods such as Raman, UV-Visible, and FTIR spectroscopy.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2024

Coupled Gel/Seeded Growth of Cs2AgBiBr6 Single Crystals for Radiation Detection

Abstract— The all-inorganic lead free Cs2AgBiBr6 perovskite represents a highly promising environment-friendly radiation detector material. The first attempt to grow Cs2AgBiBr6 single crystals using the unique gel growth method is reported in this study. This technique enables controlled reagent diffusion in a gel medium at room temperature, fostering the growth of visible single crystals. It offers advantages over conventional methods, leading to fewer equilibrium and non-equilibrium defects and simplified execution. Interestingly, the as-grown crystals exhibited variations in size and colors within the gel matrix, indicating the influence of Bi solutes in Bi-poor and Bi-rich Cs2AgBiBr6 crystals. Once the crystals were confirmed by powder xray diffraction, a series of comprehensive material characterizations was performed. The resistivity of Cs2AgBiBr6 single crystals was measured at 1.12x107 Ω•cm. Trap density analysis yielded a value of 2.82 x1010 cm-3. The charge carrier mobility was determined to be 32.9 cm2V-1S-1, suggesting the suitability for radiation detector applications. Additional growth optimization is going on and further characterization will be reported during the presentation. Our work helps accelerate the development of lead-free perovskite materials toward radiation detector applications.

James, Ralph B.↗
DOE OSTI2022

Diabatic Hamiltonian matrix elements made simple

With a view to applying the generator coordinate method to large configuration spaces, we propose a simple approximate formula to compute diabatic many-body matrix elements without having to evaluate two-body interaction matrix elements. Here, the method is illustrated with two analytically solvable Hamiltonians based on the harmonic oscillator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Analytical derivatives of the individual state energies in ensemble density functional theory. II. Implementation on graphical processing units (GPUs)

Conical intersections control excited state reactivity, and thus, elucidating and predicting their geometric and energetic characteristics are crucial for understanding photochemistry. Locating these intersections requires accurate and efficient electronic structure methods. Unfortunately, the most accurate methods (e.g., multireference perturbation theories such as XMS-CASPT2) are computationally challenging for large molecules. The state-interaction state-averaged restricted ensemble referenced Kohn–Sham (SI-SA-REKS) method is a computationally efficient alternative. The application of SI-SA-REKS to photochemistry was previously hampered by a lack of analytical nuclear gradients and nonadiabatic coupling matrix elements. We have recently derived analytical energy derivatives for the SI-SA-REKS method and implemented the method effectively on graphical processing units. We demonstrate that our implementation gives the correct conical intersection topography and energetics for several examples. Furthermore, our implementation of SI-SA-REKS is computationally efficient, with observed sub-quadratic scaling as a function of molecular size. This demonstrates the promise of SI-SA-REKS for excited state dynamics of large molecular systems.

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

Fault isolation and fault-tolerant control for Takagi-Sugeno fuzzy time-varying delay stochastic distribution systems

A fault isolation, estimation, and fault-tolerant control (FTC) scheme for nonlinear time-varying delay stochastic distribution control systems was presented in this paper. The Takagi-Sugeno fuzzy model was adopted to approach the nonlinear dynamics of time-varying delay systems. According to the output equivalence principle and Laplace transformation, an augmented state vector was given to solve the time-varying delay problem. When multiple actuator faults and interference occur simultaneously, fault detection, isolation and fault estimation was designed to obtained the fault information. To decouple faults and obtain the value and location information of the fault, the system was separated into two parts through the designed multiple conversion matrices, in which one subsystem was only affected by one actuator fault. This has simplified the design of fault isolation and estimation. A adaptive observer for fault estimation was given. Then, fault information such as the time, location, and size was determined. The observer gain matrices were calculated using linear matrix inequality (LMI). When a fault was detected and diagnosed, a FTC algorithm was devised using the proportional-integral control scheme to compensate the fault as much as possible. It has been shown that even if multiple faults actuator occurred simultaneously, the FTC controller still ensured the output probability density function of the system traced the desired probability density function when a fault occurred. Finally, the expected results were obtained through the simulation example, which confirmed the effectiveness of the method.

42 ENGINEERING↗
DOE OSTI2022

Redox-homogeneous, gel electrolyte-embedded high-mass-loading cathodes for high-energy lithium metal batteries

Lithium metal batteries have higher theoretical energy than their Li-ion counterparts, where graphite is used at the anode. However, one of the main stumbling blocks in developing practical Li metal batteries is the lack of cathodes with high-mass-loading capable of delivering highly reversible redox reactions. To overcome this issue, here we report an electrode structure that incorporates a UV-cured non-aqueous gel electrolyte and a cathode where the LiNi 0.8 Co 0.1 Mn 0.1 O 2 active material is contained in an electron-conductive matrix produced via simultaneous electrospinning and electrospraying. This peculiar structure prevents the solvent-drying-triggered non-uniform distribution of electrode components and shortens the time for cell aging while improving the overall redox homogeneity. Moreover, the electron-conductive matrix eliminates the use of the metal current collector. When a cathode with a mass loading of 60 mg cm -2 is coupled with a 100 µm thick Li metal electrode using additional non-aqueous fluorinated electrolyte solution in lab-scale pouch cell configuration, a specific energy and energy density of 321 Wh kg -1 and 772 Wh L -1 (based on the total mass of the cell), respectively, can be delivered in the initial cycle at 0.1 C (i.e., 1.2 mA cm -2 ) and 25 °C.

25 ENERGY STORAGE↗
DOE OSTI2023

Subspace recursive Fermi-operator expansion strategies for large-scale DFT eigenvalue problems on HPC architectures

Quantum mechanical calculations for material modeling using Kohn–Sham density functional theory (DFT) involve the solution of a nonlinear eigenvalue problem for N smallest eigenvector-eigenvalue pairs, with N proportional to the number of electrons in the material system. Here, these calculations are computationally demanding and have asymptotic cubic scaling complexity with the number of electrons. Large-scale matrix eigenvalue problems arising from the discretization of the Kohn–Sham DFT equations employing a systematically convergent basis traditionally rely on iterative orthogonal projection methods, which are shown to be computationally efficient and scalable on massively parallel computing architectures. However, as the size of the material system increases, these methods are known to incur dominant computational costs through the Rayleigh–Ritz projection step of the discretized Kohn–Sham Hamiltonian matrix and the subsequent subspace diagonalization of the projected matrix. This work explores the potential of polynomial expansion approaches based on recursive Fermi-operator expansion as an alternative to the subspace diagonalization of the projected Hamiltonian matrix to reduce the computational cost. Subsequently, we perform a detailed comparison of various recursive polynomial expansion approaches to the traditional approach of explicit diagonalization on both multi-node central processing unit and graphics processing unit architectures and assess their relative performance in terms of accuracy, computational efficiency, scaling behavior, and energy efficiency.

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

Modeling graphene sheet growth and dynamical matrix calculations using molecular dynamics

Molecular dynamics (MD) has been an incredibly useful tool to model physical processes that were synthesized experimentally but not fully understood. MD, through the use of semi-empirical inter-atomic potentials, has allowed understanding of different physical processes in materials science. Yet as well as providing useful insights into materials science, molecular dynamics has a wider range of usability. In this report, I will be detailing how MD can be used to study graphene formation from a carbon liquid which requires high temperatures and pressures. Beyond this, I will describe the usefulness of MD for understanding the physics for phonon transport quantum sensors. To do this, MD was employed to determine the dynamical matrix by treating atoms as coupled oscillators. An accurate understanding of the dynamical matrix of a system is required to calculate the non-equilibrium Green’s function used to describe the phonon transport within phonon wave-guides. I found that, across multiple pressures and temperatures, randomly placed carbon atoms will show evidence of pent-first formation with semi-empirical models. Density functional theory (DFT), on the other hand, was too computationally expensive to use for full scale MD simulations, but we have the possibility of training a machine learned interatomic potential to approximate DFT for carbon in the environments being studied for pent-first graphene sheet formation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Phase-Field Modeling of Mechanical Damages in Ceramic Matrix Composites

Developed a phase-field model for mechanical damages in CMCs, which incorporates the CMC microstructures, matrix cracking, fiber breakage, and interfacial sliding. Two types of mechanisms, fiber bridging and fiber pull-out, are considered. The obtained simulation results agree with experimental observations and an analytical solution. Simulation results suggest that the performance of CMCs would be enhanced with thicker fibers, longer fibers, and higher fiber density. Opposite trends of interfacial sliding resistances are suggested for the two types of situations. In reality, a mixture of the two situations may exist, and then an intermediate interfacial sliding resistance may be optimal.

Xue, Fei↗
DOE OSTI2022

Inverse Calculation of Burden Distribution Matrix Using B-spline Model Based PDF control in Blast Furnace Burden Charging Process

The inverse calculation of burden distribution matrix (BDM) is one of the most important challenges in the blast furnace operation in iron-making processes. In general, blast furnace consumes 65% of the total energy for the whole steel-making. Focusing on this practical challenge, this article proposes a new burden distribution spatial model in calculating burden charging process, and develops a B-spline approximation-based probability density function (PDF) control algorithm to assign the expected thickness distribution of burden layer and, thus, develops a new method for the required inverse calculation of BDM. First, a novel method for the thickness distribution of burden layer is given using B-spline model to produce an expected distribution shape subjected to a desired tracking within a specific spatial constraint. Then, according to the coexistence of continuous and bounded discrete variables in BDM, a novel hybrid optimization control method by combining integer programming and PDF tracking is further established for the effective inverse calculation of BDM. Finally, the proposed PDF-based iterative inverse calculation of BDM using B-spline models are tested using various data from industrial examples. Furthermore, the simulation results show that the proposed method is well suited to solve the BDM inverse calculation problem in practice.

42 ENGINEERING↗
DOE OSTI2022

Developing Reaction Chemistry Models from Reactive Molecular Dynamics: TATB

Reactive Molecular Dynamic (RMD) are used to simulate the cook-off chemistry of TATB at a variety of fixed density and fixed temperature conditions. The chemical transformations are monitored using a Coordination Geometry Analysis (CGA) approach which tracks which atom types are bonded to each specific atom. This particularly identifies oxidation state changes that occur during the transformations. Correlations between these different chemical changes are identified using a Non-negative Matrix Factorization (NMF) approach. These identify reduced order chemistry models for the TATB system which contains six components whose concentration profiles are a function of both the temperature and density/pressure. The time histories of these transformations appear to show exponential growth/decay properties that could be fit with Arrhenius rates. These components should form the basis of deflagration rate models for these materials which could then be used in mesoscale simulations to analyze accidental initiation, shock-to-detonation and detonation propagation.

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

Effect of thermal treatments on electrochemical behavior of binder jetted 17-4 PH stainless steel

Binder jetted 17–4 PH stainless steel was post-processed to relate heat treatment, microstructure, and corrosion in 3.5 wt% NaCl. Specimens were sintered at 1380 or 1400 °C, solution-annealed at 1055 °C for 1 h, and aged at 482 °C for 1 h. Here, as-sintered parts showed α′-martensitic matrix with a δ-ferrite network and Cu-rich precipitates in ferrite; inclusions (MnS, NbC) promoted localized attack. solutionizing redistributed elements and reduced ferrite, while aging generated coherent Cu nano-precipitates. Corrosion resistance was highly sensitive to post-processing in which aged specimens exhibited the lowest corrosion current density and formed a thicker, stable Cr 2 O 3 -rich passive film, whereas sintered specimens degraded most. Pitting potential depended on sintering temperature and microstructure, with 1400 °C sintering yielding more positive pitting potentials and the best overall performance after aging. Although the 1400-solutionized condition showed a relatively noble pitting response versus 1380-solutionized, it displayed unstable corrosion kinetics attributed to an imperfect passive film linked to higher NbC density. XPS depth profiles corroborated these trends, showing thicker, more continuous Cr 2 O 3 in aged states and discontinuous/thinning oxides in less resistant conditions. Practically, high-temperature sintering (∼1400 °C) followed by solutionizing and aging is recommended, with further gains expected from reducing NbC/MnS populations and porosity via powder and process control.

36 MATERIALS SCIENCE↗
DOE OSTI2024

ElasTool v3.0: Efficient computational and visualization toolkit for elastic and mechanical properties of materials

Efficient computation and visualization of elastic and mechanical properties are crucial in the selection of materials and the design of new materials. Here, the ElasTool v3.0 toolkit marks a significant advancement in the computational analysis and visualization of elastic and mechanical properties of materials, essential in material selection and design. This enhanced version extends beyond standard calculations like elastic tensor, Young's modulus, bulk modulus, and Poisson's ratio. It introduces capabilities for computing minimum thermal conductivity, linear compressibility, rendering the Christoffel equation, and elastic energy density. Notably, it integrates advanced visualization tools, including compatibility with Plotly and Elate web platforms for interactive web-based property exploration. A key feature of ElasTool v3.0 is the implementation of second-order elastic constants (SOECs) for tubular 2D-based nanostructures and nanotubes. Leveraging high-efficiency strain-matrix sets (OHESS), the toolkit now facilitates efficient computation of elastic constants and mechanical properties at both zero and finite temperatures for 1D, 2D, and 3D dimensions. ElasTool is openly accessible on GitHub: https://github.com/gmp007/elastool.

1D, 2D, 3D, and tubular 2D nanostructure and nanot↗
DOE OSTI2024

Nanoscale Ion Transport Enhances Conductivity in Solid Polymer-Ceramic Lithium Electrolytes

The predictive design of flexible and solvent-free polymer electrolytes for solid-state batteries requires an understanding of the fundamental principles governing the ion transport. In this work, we establish a correlation among the composite structures, polymer segmental dynamics, and lithium ion (Li + ) transport in a ceramic-polymer composite. Elucidating this structure–property relationship will allow tailoring of the Li + conductivity by optimizing the macroscopic electrochemical stability of the electrolyte. The ion dissociation from the slow polymer segmental dynamics was found to be enhanced by controlling the morphology and functionality of the polymer/ceramic interface. The chemical structure of the Li + salt in the composite electrolyte was correlated with the size of the ionic cluster domains, the conductivity mechanism, and the electrochemical stability of the electrolyte. Polyethylene oxide (PEO) filled with lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) or lithium bis(fluorosulfonyl) imide (LiFSI) salts was used as a matrix. A garnet electrolyte, aluminum substituted lithium lanthanum zirconium oxide (Al-LLZO) with a planar geometry, was used for the ceramic nanoparticle moieties. Further, the dynamics of the strongly bound and highly mobile Li + were investigated using dielectric relaxation spectroscopy. The incorporation of the Al-LLZO platelets increased the number density of more mobile Li + . The structure of the nanoscale ion-agglomeration was investigated by small-angle X-ray scattering, while molecular dynamics (MD) simulation studies were conducted to obtain the fundamental mechanism of the decorrelation of the Li + in the LiTFSI and LiFSI salts from the long PEO chain.

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

Vertical Architecture Solution-Processed Quantum Dot Photodetectors with Amorphous Selenium Hole Transport Layer

Colloidal quantum dots (CQDs) provide wide spectral tunability and high absorption coefficients owing to quantum confinement and large oscillator strengths, which along with solution processability, allow a facile, low-cost, and room-temperature deposition technique for the fabrication of photonic devices. However, many solution-processed CQD photodetector devices demonstrate low specific-detectivity and slow temporal response. To achieve improved photodetector characteristics, limiting carrier recombination and enhancing photogenerated carrier separation are crucial. In this study, we develop and present an alternate vertical-stack photodetector wherein we use a solution-processed quantum dot photoconversion layer coupled to an amorphous selenium (a-Se) wide-bandgap charge transport layer that is capable of exhibiting single-carrier hole impact ionization and is compatible with active-matrix readout circuitry. This a-Se chalcogenide transport layer enables the fabrication of high-performance and reliable solution-processed quantum dot photodetectors, with enhanced charge extraction capabilities, high specific detectivity (D* ~ 0.5–5 × 10 12 Jones), fast 3 dB electrical bandwidth (3 dB BW ~ 22 MHz), low dark current density (J D ~ 5–10 pA/cm 2 ), low noise current (i n ~ 20–25 fW/Hz 1/2 ), and high linear dynamic range (LDR ~130–150 dB) across the measured visible electromagnetic spectrum (~405–656 nm).

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Photo-actuators via epitaxial growth of microcrystal arrays in polymer membranes

Photomechanical crystals composed of three-dimensionally ordered and densely packed photochromes hold promise for high-performance photochemical actuators. However, bulk crystals with high structural ordering are severely limited in their flexibility, resulting in poor processibility and a tendency to fragment upon light exposure, while previous nano- or microcrystalline composites have lacked global alignment. Here we demonstrate a photon-fuelled macroscopic actuator consisting of diarylethene microcrystals in a polyethylene terephthalate host matrix. These microcrystals survive large deformations and show a high degree of three-dimensional ordering dictated by the anisotropic polyethylene terephthalate, which critically also has a similar stiffness. Overall, these ordered and compliant composites exhibit rapid response times, sustain a performance of over at least hundreds of cycles and generate work densities exceeding those of single crystals. In conclusion, our composites represent the state-of-the-art for photochemical actuators and enable properties unattainable by single crystals, such as controllable, reversible and abrupt jumping (photosalient behaviour).

36 MATERIALS SCIENCE↗
DOE OSTI2022

Impact of CH 4 addition on the electron properties and electric field dynamics in a Ar nanosecond-pulsed dielectric barrier discharge

Non-equilibrium plasmas derive their low temperature reactivity from producing and driving energetic electrons and active species under large electric fields. Therefore, the impact of reactants on the plasma properties including electron number density, electric field, and electron temperature is critical for applications such as plasma CH 4 reforming. Due to experimental complexity, electron properties and the electric field are rarely measured together in the same discharge. Here, we combine time-resolved Thomson scattering and electric field induced second harmonic generation (EFISH) to probe electron temperature, electron density, and electric field strength in a 60 Torr CH 4 /Ar nanosecond-pulsed dielectric barrier discharge (ns-DBD) while varying the CH 4 mole fraction from 0 to 8%. These measurements are compared to a 1-D numerical model to benchmark its predictions and identify areas of uncertainty. Nonlinear coupling between CH 4 addition, electron temperature, electron density, and the electric field was directly observed. Contrary to previous measurements in He, the electron temperature increased with CH 4 mole fraction. This rise in electron temperature is identified as electron heating by residual electric fields that increased with larger CH 4 mole fraction. Moreover, the electron number density has been found to decrease rapidly with the increase of methane mole fraction. Comparison of these measurements with the model yielded better agreement at higher CH 4 mole fractions and with the usage of ab initio calculated Ar electron-impact cross-sections from the B-spline R-matrix (BSR) database. Furthermore, the calculated plasma properties are shown to be sensitive to the residual surface charge implanted on the quartz dielectric surfaces. Without considering surface charge in the simulations, the calculated electric field profiles agreed well with the measurements, but the electron properties were underpredicted by more than a factor of three. Therefore, measurements of either the electric field or electron properties measurements alone are insufficient to fully validate modelling predictions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗