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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Radionuclide Screening Analysis and Transport Parameters for Pahute Mesa Detonations, Nevada National Security Site

A “screening analysis” is implemented at Pahute Mesa (PM), Nevada National Security Site (NNSS), to determine subsets of the 43 radionuclides listed in the radionuclide inventory of Finnegan et al. (2016) that are of relevance or potential relevance to the hydrologic source term (HST) for assessment of radionuclide transport in groundwater. Consideration is added to how levels of contamination are defined and whether any other radionuclides not in the inventory are relevant to the HST. A model is developed to estimate a range of possible source concentrations in groundwater that account for uncertainty in partitioning into melt glass and sorption into surrounding rock of the exchange zone for the 82 PM underground nuclear tests detonated in vertical shafts. The model also accounts for the varied hydrogeochemical settings. Transport parameters needed for the screening model calculations are developed from databases for hydrogeologic units, chemistry, mineralogy, fracture spacing, fracture aperture, fracture openness, matrix porosity, bulk density, saturation, and alteration. Of key importance is consideration of diagenetic zonation of the mineralogy related to hydrothermal alteration. The screening analysis compares model results with available groundwater radiochemistry data. Consistency between screening model results and available data helps validate the model for application to all 43 radionuclides, most of which have no measurements or only non-detect measurements of concentration in groundwater. The screening analysis determines that ten radionuclides – tritium, Sr-90, I-129, Cs-137, U-232, U-233, U-234, U-238, and the total of Pu-239 and Pu-240 are relevant to the HST. Determination of relevance is primarily based on data and/or model results indicating source concentration exceeding a maximum contaminant level (MCL) out to 100 years since the time of detonation. U-238 is relevant to assessment of contamination relative to MCLs for uranium and gross alpha particle activity. Five other radionuclides – C-14, Cl-36, Tc-99, Np-237, and Pu-238 – are determined potentially relevant to the HST based on a 0.1 MCL threshold. The screening analysis recommends additional attention to the natural daughter radionuclides of the uranium series (U-238) and thorium series (Th-232) decay chains, particularly Ra-226 and Ra-228, that were not included the inventory.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

A Study of a Surrogate Nuclear Thermal Propulsion Fuel Element

We report a surrogate fuel, Zirconium Carbide/Titanium Nitride has been prepared using a simple method to obtain a more uniform distribution of titanium nitride in the zirconium carbide matrix. Powders were sintered using Spark Plasma Sintering at 1750oC and 1850oC using a modified sintering profile. This resulted in a fairly uniform distribution of titanium nitride in the zirconium carbide matrix with a density of 95% of theoretical at 1850oC. Scanning electron microscopy, indentation hardness and x-ray diffraction results are given.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Binder for electrolyte matrix for molten carbonate fuel cells

A binder solution for an electrolyte matrix for use with molten carbonate fuel cells is provided. The binder solution includes a first polymer with a molecular weight of less than about 150,000 and a second binder with a molecular weight of greater than about 200,000. The binder solution produces an electrolyte matrix with improved flexibility, matrix particle packing density, strength, and pore structure.

Surendranath, Arun↗
DOE OSTI2020

Methodology for Generating Covariance Data of Thermal Neutron Scattering Cross Sections

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

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Analytic techniques for solving the transport equations in electroweak baryogenesis

We develop an efficient method for solving transport equations, particularly in the context of electroweak baryogenesis. It provides fully-analytical results under mild approximations and can also test semi-analytical results, which are applicable in more general cases. Key elements of our method include the reduction of the second-order differential equations to first order, representing the set of coupled equations as a block matrix of the particle densities and their derivatives, identification of zero modes, and block decomposition of the matrix. We apply our method to calculate the baryon asymmetry of the Universe (BAU) in a Standard Model effective field theory framework of complex Yukawa couplings to determine the sensitivity of the resulting BAU to modifications of various model parameters and rates, and to estimate the effect of the commonly-used thin-wall approximation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Block Mott insulating state induced by next-nearest-neighbor hopping in the S = 3 2 zigzag chain BaCoTe 2 O 7

Quasi-one-dimensional correlated electronic multiorbital systems with either ladder or chain geometries continue attracting considerable interest due to their complex electronic phases arising from the interplay of the hopping matrix, the crystal-field splitting, the electronic correlations (Hubbard repulsion U and Hund coupling J H ), and strong quantum fluctuations. Recently, the intriguing cobalt zigzag chain system BaCoTe 2 ⁢ O 7 , with electronic density n=7, was prepared experimentally. Here, we systematically study the electronic and magnetic properties of this quasi-one-dimensional compound from the theoretical perspective. Based on first-principles density functional theory calculations, strongly anisotropic one-dimensional electronic Co 3⁢d bands were found near the Fermi level. By evaluating the relevant hopping amplitudes, we provide the magnitude and origin of the nearest-neighbor (NN) and next-nearest-neighbor (NNN) hopping matrices in BaCoTe 2 ⁢ O 7 . With this information, we constructed a three-orbital electronic Hubbard model for this zigzag chain system, and studied two cases: with only a NN hopping matrix, and with NN plus NNN hopping matrices. Introducing the Hubbard and Hund couplings and studying the model via the density matrix renormalization group method, we constructed the ground-state phase diagram. A robust staggered ↑ - ↓ - ↑ - ↓ antiferromagnetic (AFM) region was found when only the NN hopping matrix in the chain direction was employed. However, for the realistic case where the NNN hopping matrix is also included, the dominant state becomes instead a block AFM ↑ - ↑ - ↓ - ↓ order, in agreement with experiments. The system displays Mott insulator characteristics with three half-filled orbitals, when the block AFM order is stable. In conclusion, our results for BaCoTe 2 ⁢O 7 provide guidance to experimentalists and theorists working on this zigzag one-dimensional chain and related materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Accelerating self-consistent field iterations in Kohn-Sham density functional theory using a low-rank approximation of the dielectric matrix

We present an efficient preconditioning technique for accelerating the fixed-point iteration in real-space Kohn-Sham density functional theory (DFT) calculations. The preconditioner uses a low-rank approximation of the dielectric matrix (LRDM) based on Gâteaux derivatives of the residual of fixed-point iteration along appropriately chosen direction functions. We develop a computationally efficient method to evaluate these Gâteaux derivatives in conjunction with the Chebyshev filtered subspace iteration procedure, an approach widely used in large-scale Kohn-Sham DFT calculations. Further, we propose a variant of LRDM preconditioner based on adaptive accumulation of low-rank approximations from previous self-consistent field iterations, and also extend the LRDM preconditioner to spin-polarized Kohn-Sham DFT calculations. We demonstrate the robustness and efficiency of the LRDM preconditioner against other widely used preconditioners on a range of benchmark systems with sizes ranging from ~100 to 1100 atoms (~500–20,000 electrons). The benchmark systems include various combinations of metal-insulating-semiconducting heterogeneous material systems, nanoparticles with localized d orbitals near the Fermi energy, nanofilm with metal dopants, and magnetic systems. In all benchmark systems, the LRDM preconditioner converges robustly within 20–30 iterations. In contrast, other widely used preconditioners show slow convergence in many cases, as well as divergence of the fixed-point iteration in some cases. Lastly, we demonstrate the computational efficiency afforded by the LRDM method, with up to 3.4-fold reduction in computational cost for the total ground-state calculation compared to other preconditioners.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Evaluation of a Geotechnical Compaction Model to Estimate Sediment Deformation

When surface ruts from vehicle activity occur above archeological deposits, there is a potential for an impact to the archeological record to occur. While impacts to the soil surface are relatively straightforward to describe, sub-surface disturbances are more difficult to quantify. In particular, disturbances to buried archeological materials can be difficult to characterize. Materials can be damaged by surface disturbances, but they can also be impacted even if they stay intact. Interpretation of intact materials relocated within the sediment column can be problematic. Traditional characterization methods (e.g. excavation) can also disturb archaeological materials. Here, we present a geotechnical compaction model that can be used to assist with sub-surface characterization under vehicle ruts formed during off-road travel. The model is based on a Boussinesq analysis of stress increase below a surface load. The change in bulk density resulting from compaction is assumed to be proportional to the stress increase. The model results were compared against pilot scale tests conducted with homogeneous sediment and field scale in-situ tests. The results indicate that the model generally agreed with measurements to within 20%, although only a single soil type was evaluated in these studies, with a narrow range of initial bulk densities and rut depths. While application of the model across a broader matrix of soil types, bulk densities and rut depths would increase confidence in the model, it appears that the model may be suitable for use in archeological field investigations. Potential uses include both compliance and research based field applications.

Fritz, Brad G.↗
DOE OSTI2023

Mechanical and Thermal Characterization of Phase-Change Material and High-Density Polyethylene Functional Composites for Thermal Energy Storage

Abstract Phase-change materials (PCMs) can be used to develop thermal energy storage systems as they absorb large amount of latent heat nearly at a constant temperature when changing phase from a solid to a liquid. To prevent leakage when in a liquid state, PCM is shape stabilized in a polymer matrix of high-density polyethylene (HDPE). The present research explores the injection-molded mechanical and thermal properties of different PCM/HDPE composite ratios. The tensile strength and modulus of elasticity at room temperature and with the PCM fully melted within the composite are measured. Additionally, the hardness, latent heat of fusion, phase-change temperature, and thermal conductivity are investigated. An analysis of microstructures of the composite is used to support the findings. The PCM within the PCM/HDPE composite gives it the benefit of thermal storage but causes a decrease in mechanical properties.

Energy & Fuels↗
DOE OSTI2020

Charge-density-wave melting in the one-dimensional Holstein model

We study the Holstein model of spinless fermions, which at half filling exhibits a quantum phase transition from a metallic Tomonaga-Luttinger liquid phase to an insulating charge-density-wave (CDW) phase at a critical electron-phonon coupling strength. In our work, we focus on the real-time evolution starting from two different types of initial states that are CDW ordered: (i) ideal CDW states with and without additional phonons in the system and (ii) correlated ground states in the CDW phase. We identify the mechanism for CDW melting in the ensuing real-time dynamics and show that it strongly depends on the type of initial state. We focus on the far-from-equilibrium regime and emphasize the role of electron-phonon coupling rather than dominant electronic correlations, thus complementing a previous study of photoinduced CDW melting [H. Hashimoto and S. Ishihara, Phys. Rev. B 96, 035154 (2017)]. The numerical simulations are performed by means of matrix-product-state based methods with a local basis optimization (LBO). Within these techniques, one rotates the local (bosonic) Hilbert spaces adaptively into an optimized basis that can then be truncated while still maintaining a high precision. In this work, we extend the time-evolving block decimation (TEBD) algorithm with LBO, previously applied to single-polaron dynamics, to a half-filled system. We demonstrate that in some parameter regimes, a conventional TEBD method without LBO would fail. Furthermore, we introduce and use a ground-state density-matrix renormalization group method for electron-phonon systems using local basis optimization. In our examples, we account for up to M ph =40 bare phonons per site by working with O(10) optimal phonon modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Elastic proton scattering off nonzero spin nuclei

In recent years, we constructed a microscopic optical potential (OP) for elastic nucleon-nucleus ( NA ) scattering using modern approaches based on chiral theories for the nucleon-nucleon ( NN ) interaction. The OP was derived at first order of the spectator expansion in Watson multiple scattering theory and its final expression was a folding integral between the NN t matrix and the nuclear density of the target. Two- and three-body forces are consistently included both in the target and in the projectile description. The purpose of this work is to apply our microscopic OP to nuclei characterized by a ground state of spin-parity quantum numbers J π ≠ 0 + . We extended our formalism to include the spin of the target nucleus. The full amplitudes of the NN reaction matrix are retained in the calculations starting from two- and three-body chiral forces. The microscopic OP can be applied in the energy range 100 ≤ E ≤ 350 MeV. Finally, we show a remarkable agreement with experimental data for the available observables and, simultaneously, provide reliable estimates for the theoretical uncertainties. This work paves the way toward a full microscopic approach to inelastic NA scattering, showing that the derivation of optical potentials between states with J π ≠ 0 + is completely under control.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Understanding the setting and hardening process of wollastonite-based brushite cement. Part 2: Influence of the boron and aluminum concentrations in the mixing solution

The mixing solution of wollastonite-based brushite cement is a phosphoric acid solution containing metallic cations and borax. This work complements a previous study devoted to the influence of the H{sub 3}PO{sub 4} concentration, Ca/P and liquid-to-solid (l/s) ratios on the setting and hardening process of the binder by providing new insight into the role of aluminum and boron. Boron retards the setting and decreases the heat released during the process. It also contributes to reduce the macroporosity of the hardened material but yields to poor compressive strength. With aluminum in the mixing solution, the mechanical properties are greatly improved thanks to the precipitation of an amorphous aluminophosphate which increases the density of the cement matrix. But aluminum alone leads to fast setting. A joint addition of boron and aluminum to the mixing solution makes it possible to get a material with optimized properties both in fresh and hardened states.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Noise and error analysis and optimization in particle-based kinetic plasma simulations

In this paper we analyze the noise in macro-particle methods used in plasma physics and fluid dynamics, leading to approaches for minimizing the total error, focusing on electrostatic models in one dimension. We begin by describing kernel density estimation for continuous values of the spatial variable x, expressing the kernel in a form in which its shape and width are represented separately. The covariance matrix of the noise in the density is computed, first for uniform true density. The bandwidth of the covariance matrix C(x,y) is related to the width of the kernel. A feature that stands out is the presence of constant negative terms in the elements of the covariance matrix both on and off-diagonal. These negative correlations are related to the fact that the total number of particles is fixed at each time step; they also lead to the property ∫C(x,y)dy = 0. We investigate the effect of these negative correlations on the electric field computed by Gauss's law, finding that the noise in the electric field is related to a process called the Ornstein-Uhlenbeck bridge, leading to a covariance matrix of the electric field with variance significantly reduced relative to that of a Brownian process. For non-constant density, p(x), still with continuous x, we analyze the total error in the density estimation and discuss it in terms of bias-variance optimization (BVO). For some characteristic length l, determined by the density and its second derivative, and kernel width h, having too few particles within h leads to too much variance; for h that is large relative to l, there is too much smoothing of the density. The optimum between these two limits is found by BVO. For kernels of the same width, it is shown that this optimum (minimum) is weakly sensitive to the kernel shape. Next, we repeat the analysis for x discretized on a grid. In this case the charge deposition rule is determined by a particle shape. An important property to be respected in the discrete system is the exact preservation of total charge on the grid; this property is necessary to ensure that the electric field is equal at both ends, consistent with periodic boundary conditions. We find that if the particle shapes satisfy a partition of unity property, the particle charge deposited on the grid is conserved exactly. Further, if the particle shape is expressed as the convolution of a kernel with another kernel that satisfies the partition of unity, then the particle shape obeys the partition of unity. This property holds for kernels of arbitrary width, including widths that are not integer multiples of the grid spacing. Furthermore, we show results relaxing the approximations used to do BVO optimization analytically, by doing numerical computations of the total error as a function of the kernel width, on a grid in x. The comparison between numerical and analytical results shows good agreement over a range of particle shapes. We discuss the practical implications of our results, including the criteria for design and implementation of computationally efficient particle shapes that take advantage of the developed theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Elastic Antiproton-Nucleus Scattering from Chiral Forces

Elastic scattering of antiprotons off He 4 , C 12 , and O 16 , 18 is described for the first time with a consistent microscopic approach based on the calculation of an optical potential (OP) describing the antiproton-target interaction. The OP is derived using the recent antiproton-nucleon ( p ¯ N ) chiral interaction to calculate the p ¯ N t matrix, while the target densities are computed with the ab initio no-core shell model using chiral interactions as well. Here, our results are in good agreement with the existing experimental data and the results computed at different chiral orders of the p ¯ N interaction display a well-defined convergence pattern.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Preliminary Thermal-hydraulics and Selected Safety Analysis for Holos-Quad Reactor Design

This report presents thermal hydraulics and safety analyses of the 10 MWe Holos-Quad micro-reactor design developed by HolosGen LLC. These analyses were executed under steady-state and station black-out (SBO) transient conditions using the System Analysis Module (SAM), with a focus on ensuring peak fuel temperatures do not exceed the safety thresholds of 1,250°C for steady state normal operating conditions and 1,600°C for transient conditions. For steady-state temperature predictions, a coupled 1D fluid-to-3D solid heat conduction methodology was used. The steady-state simulations included a model of the central core region where power peaking occurs to predict the maximum steady-state fuel temperature, and a coarse-mesh full-core model to provide the initial condition for the consequent transient SBO decay heat removal simulation. The maximum fuel temperature predicted from the central core region simulations under the normal operating condition was 1,222°C, which is below the design limit 1,250°C. To simulate a SBO transient, a full core model was necessary because of the non-axial symmetric core design and heat loss from the outer core structures to the environment. To reduce the required computational costs, a novel scaled subassembly approach was used to model the full core domain. With the full core steady-state temperature distribution provided as the initial condition, the passive decay heat removal of the core was simulated with the scaled subassembly full core model under SBO transient conditions. This analysis showed that the peak fuel temperatures of the Holos-Quad core remained below their steady-state values during the transient, therefore much lower than the safety limit 1,600°C. This is mainly due to the large thermal inertia of the graphite matrix, relatively low power density, and the large surface-to-volume-ratio of the core. The preliminary analyses presented in this report confirm the inherent safety of the Holos-Quad micro-reactor thermal-hydraulics design, as the peak fuel temperatures under both normal and off-normal operational conditions remain within design and safety limits.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Graphene oxide/metal nanocrystal multilaminates the atomic limit for safe, selective hydrogen storage

Various embodiments of the invention describe an environmentally stable, and exceptionally dense hydrogen storage (6.5 wt % and 0.105 kg H2/L in the total composite, 7.56 wt % in Mg) using atomically thin and gas-selective reduced graphene oxide sheets as encapsulants. Other approaches to protecting reactive materials involve energy intensive introduction of considerable amounts of inactive, protective matrix which compromises energy density. However, these multilaminates are able to deliver exceptionally dense hydrogen storage far-exceeding 2020 DOE target metrics for gravimetric capacity (5.5 wt %), and ultimate full-fleet volumetric targets (0.070 kg H 2 /L) for fuel cell electric vehicles. Methods of stabilizing reactive nanocrystalline metals in zero-valency also has wide-ranging applications for batteries, catalysis, encapsulants, and energetic materials.

Urban, Jeffrey J.↗
DOE OSTI2024

Controlling interactions between high-frequency phonons and single quantum systems using phononic crystals

The ability to control phonons in solids is key in many fields of quantum science, ranging from quantum information processing to sensing. Phonons often act as a source of noise and decoherence when solid-state quantum systems interact with the phonon bath of their host matrix. In this study, we demonstrate the ability to control the phononic local density of states of the host matrix using phononic crystals and measure its positive impact on single quantum systems. We design and fabricate diamond phononic crystals with features down to around 20 nm, resulting in a high-frequency complete phononic bandgap from 50 to 70 GHz. The engineered local density of states is probed using single silicon-vacancy colour centres embedded in the phononic crystals. We observe an 18-fold reduction in the phonon-induced orbital relaxation rate of the emitters compared to bulk, thereby demonstrating that the phononic crystal suppresses spontaneous single-phonon processes. Furthermore, we show that our approach can efficiently suppress single-phonon-emitter interactions up to 20 K, allowing the investigation of multi-phonon processes in the emitters. Our results represent an important step towards the realization of efficient phonon-emitter interfaces that can be used for quantum acoustodynamics and quantum phononic networks.

quantum information↗
DOE OSTI2023

Nuclear-Level Effective Theory of 𝜇 → 𝑒 Conversion

The Mu2e and COMET 𝜇 → 𝑒 conversion experiments are expected to significantly advance limits on new sources of charged lepton flavor violation. Almost all theoretical work in the field has focused on just two operators. However, general symmetry arguments lead to a 𝜇 → 𝑒 conversion rate with six response functions, each of which, in principle, is observable by varying nuclear properties of targets. We construct a nucleon-level nonrelativistic effective theory (NRET) to clarify the microscopic origin of these response functions and to relate rate measurements in different targets. This exercise identifies three operators and their small parameters that control the NRET operator expansion. We note inconsistencies in past treatments of these parameters. The NRET is technically challenging, involving 16 operators, several distorted electron partial waves, bound muon upper and lower components, and an exclusive nuclear matrix element. We introduce a trick for treating the electron Coulomb effects accurately, which enables us to include all of these effects while producing transition densities whose one-body matrix elements can be evaluated analytically, greatly simplifying the nuclear physics. We derive bounds on operator coefficients from existing and anticipated 𝜇 → 𝑒 conversion experiments. We discuss how similar NRET formulations have impacted dark matter phenomenology, noting that the tools this community has developed could be adapted for charged lepton flavor violation studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗