Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix elements”

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

Anomalous Hall effect in conical helimagnetic crystals

Spin-spiral texture can substantially change charge transport properties in helimagnets. Here we find the anomalous Hall effect (AHE) exhibiting the dramatic behavior with respect to chemical potential μ in conical magnetic structures. The direct conductivity demonstrates kinks, and the anomalous Hall current exhibits minima and maxima changing the sign. We analytically derive the expression for energy bands and eigenstates in the most general case. Because of the conical potential, the energy bands are split into two nonparabolic bands where the lower band can have one- or two-minima shapes in the k z direction ($\hat{z}$ is a direction of the spiral axis). We prove that the origin of the anomalous Hall effect is not topological and is due to the interplay between the asymmetry of energy bands in the x and z directions and spin restrictions in the phase space due to the conical potential. We also investigate the dependence of transport properties on cone half- angle θ, and find that the effects are most pronounced at θ=π/2 (a helical state). Electric current is calculated using the Boltzmann equation where the relaxation is caused by electron-acoustic phonon interaction. The transition probability is found to be a 2×2 matrix with nonvanishing off-diagonal elements indicating the strong interband transitions. The origin of interband transitions is because of the nature of the conical potential where conduction electron spins interact with localized magnetic moments. To verify the proposed theory, we calculate the temperature dependence of resistivity for MnSi crystals and find the discontinuity at the phase transition between conical and paramagnetic phases. The calculations are in the excellent agreement with the experimental data. In addition, we predict the discontinuity behavior for the anomalous Hall resistivity at the phase transition where the resistivity exhibits the abrupt change at T=T C , (a) to zero if the relativistic effects for the conduction electrons are small or (b) to a nonzero value if Rashba/Dresselhaus effects are taken into account.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Simulation and Reynolds-averaged Navier-Stokes modeling of a three-component Rayleigh-Taylor mixing problem with thermonuclear burn

Rayleigh-Taylor mixing in the presence of a third component with intermediate density is investigated through three-dimensional large-eddy simulation (LES) with a high-order compact finite-difference code. Two configurations are considered: (1) a symmetric configuration in which the Atwood number between the heavy and intermediate components matches the Atwood number between the intermediate and light components and (2) an asymmetric configuration in which the Atwood number between the heavy and intermediate components is an order of magnitude greater than the Atwood number between the intermediate and light components. Mass fraction covariances are extracted, and proposed Reynolds-averaged Navier-Stokes (RANS) closures for density-specific-volume and density-mass-fraction covariances are evaluated in an a priori fashion. Additionally, a multicomponent extension of the k - Φ - L - a - V RANS model [Morgan, Phys. Rev. E 104, 015107 (2021)] is presented which includes model equations for the upper-triangular elements of the mass fraction covariance matrix. This model, referred to as the k - Φ - L - a - C model, is compared against results from LES and against other RANS models. Profiles of average mass fraction, mass-fraction covariance, and density-specific-volume covariance obtained with the k - Φ - L - a - C model are found to agree well with LES data. Finally, the impact of three-component turbulent mixing on average reaction rate is investigated in both premixed and nonpremixed cases by heating the mixing layer and allowing it to undergo thermonuclear (TN) burn. A closure model for average reaction rate is proposed for use with the k - Φ - L - a - C model, and when this model is applied, improved agreement is obtained between LES and RANS in total TN neutron production.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Ultrafast Molecular Frame Quantum Tomography

Here, we develop and experimentally demonstrate a methodology for a full molecular frame quantum tomography (MFQT) of dynamical polyatomic systems. We exemplify this approach through the complete characterization of an electronically nonadiabatic wave packet in ammonia (NH 3 ). The method exploits both energy and time-domain spectroscopic data, and yields the lab frame density matrix (LFDM) for the system, the elements of which are populations and coherences. The LFDM fully characterizes electronic and nuclear dynamics in the molecular frame, yielding the time- and orientation-angle dependent expectation values of any relevant operator. For example, the time-dependent molecular frame electronic probability density may be constructed, yielding information on electronic dynamics in the molecular frame. In NH 3 , we observe that electronic coherences are induced by nuclear dynamics which nonadiabatically drive electronic motions (charge migration) in the molecular frame. Here, the nuclear dynamics are rotational and it is nonadiabatic Coriolis coupling which drives the coherences. Interestingly, the nuclear-driven electronic coherence is preserved over longer timescales. In general, MFQT can help quantify entanglement between electronic and nuclear degrees of freedom, and provide new routes to the study of ultrafast molecular dynamics, charge migration, quantum information processing, and optimal control schemes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Performance of High Stopping Power Bismuth-Loaded Plastic Scintillators for Radiation Portal Monitors

Plastic scintillators are widely used in radiation portal monitors because of their low cost and availability in large sizes. However, due to their low density and low effective atomic number (Z), they offer low intrinsic efficiency and little spectroscopic information. The addition of high-Z constituents to these plastics can greatly increase both their total stopping power and the amount of photoelectric absorption, leading to full-energy deposition and thus useful gamma spectra. Herein, we present the performance of the latest formulation of Bi-loaded plastic scintillators showing their useful spectroscopic information up to relatively high energy (~1 MeV) due to their high stopping power compared to the current commercially available plastics. These Bi-loaded plastics use 20 weight percent (wt%) Bi-pivalate (8 wt% elemental Bi) dissolved in polyvinyltoluene (PVT) matrix and conventional fast fluors (~10 ns decay time). These Bi-loaded plastics achieve up to approximately 6000 photons/MeV and have been produced in sizes up to 17 in 3 . The performance of these Bi-loaded plastics is also demonstrated in the existing portal monitor hardware (Rapiscan Model TSA Trainer 770) showing the possibility to provide improved sensitivity as a drop-in replacement with continued scale-up.

42 ENGINEERING↗
DOE OSTI2021

New Modulation and Impact of Transformer Leakage Inductance on Current-Source Solid-State Transformer

This article presents a novel modulation scheme for device voltage stress mitigation and comprehensive analysis of the impact of transformer leakage inductance in a current-source solid-state transformer (SST). Different from dual active bridge (DAB) SST, the operation of the current-source SST is similar to that of a flyback converter. The device bridge on only one side of the transformer is active to store energy into or release energy from the magnetizing inductance which acts as a current-source dc link. Such flyback operation with reverse-blocking switches can lead to additional device voltage stress and incomplete zero-voltage switching (ZVS) on the current-source soft-switching solid-state transformer (S4T) under conventional modulation. A new modulation scheme is proposed to address this issue. Moreover, different from the DAB, the leakage inductance of the medium-frequency transformer (MFT) in the S4T is a parasitic element similar to that in a matrix SST and can cause additional device voltage stress. Though the resonant capacitors, originally added to achieve ZVS, can absorb and recycle the leakage energy in the S4T, these capacitors need to be increased with larger leakage inductances to limit the voltage stress. However, large resonant capacitors can result in more lost duty cycles and reduced efficiency. The impact of such leakage inductance on device voltage stress is analyzed comprehensively, which is critical to guide future research and design of the S4T. Experimental results from S4T prototypes for DC-DC, multiport AC-DC, and AC-AC conversion with 1:1 and 4:1 MFT comprehensively verify the proposed concepts. Lastly, a case study of a three-phase AC-AC S4T over a power range from 1 kVA to 100 kVA each module reveals that the MFT leakage inductance should be less than 1% of the magnetizing inductance for safe operation.

14 SOLAR ENERGY↗
DOE OSTI2021

Molten salt compositions with enhanced heat transfer and reduced corrosion properties

A heat transfer (exchange) composition comprising a halide salt matrix having dispersed therein nanoparticles comprising elemental carbon in the absence of water and surfactants, wherein said halide is fluoride or chloride, wherein the halide salt may be an alkali halide salt (e.g., lithium fluoride, sodium fluoride, potassium fluoride, rubidium fluoride, sodium chloride, potassium chloride, rubidium chloride, and eutectic mixtures thereof) or an alkaline earth halide salt (e.g., fluoride or chloride salt of beryllium, magnesium, calcium, strontium, or barium), and wherein the nanoparticles comprising elemental carbon may be solid or hollow, and wherein the composition may further include nanoparticles comprising a fissile material (e.g., U, Th, or Pu) dispersed within the composition. Molten salt reactors (MSRs) containing these heat transfer compositions in coolant loops in thermal exchange with a reactor core, as well operation of such MSRs, are also described.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Formation and Consolidation of UN Nanostructures in U-Mo Fuel via Mechanical Alloying and Spark Plasma Sintering

Uranium-molybdenum (U-Mo) alloys show promise as a nuclear fuel system due to their high thermal conductivity and high fuel loading capability. However, U-Mo systems suffer from irradiation induced swelling ultimately affecting the cladding via mechanical and chemical interaction. To address these shortcomings, this thesis investigated the formation of uranium mononitride (UN) nanoparticles within a U-10Mo matrix. To promote the formation of UN, U-10Mo powders were mechanically alloyed under a high purity nitrogen cover gas (99.9995%) atmosphere. Variations were made in stainless steel milling media size, duration of milling, and number of times the milling jar was re-aerated with nitrogen gas. UN nanoparticulates were successfully formed within the U-Mo matrix and was characterized utilizing light element analysis, X-ray diffractometry, scanning and transmission-electron microscopy, electron energy loss spectroscopy, and atom probe tomography. Presence of the UN nanoparticles was found as early as 1-hour into the mechanical alloying process, and significant iron contamination found after 10-hours of milling. The U-Mo/UN powders were compacted using spark plasma sintering techniques in order to stabilize the nanostructure. U-Mo/UN compacts containing small amounts of neodymium and helium were analyzed to observe the defect migration tendencies of fission products within the fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

End-to-end GPU acceleration of low-order-refined preconditioning for high-order finite element discretizations

In this article, we present algorithms and implementations for the end-to-end GPU acceleration of matrix-free low-order-refined preconditioning of high-order finite element problems. The methods described here allow for the construction of effective preconditioners for high-order problems with optimal memory usage and computational complexity. The preconditioners are based on the construction of a spectrally equivalent low-order discretization on a refined mesh, which is then amenable to, for example, algebraic multigrid preconditioning. The constants of equivalence are independent of mesh size and polynomial degree. For vector finite element problems in H(curl) and H(div) (e.g., for electromagnetic or radiation diffusion problems), a specially constructed interpolation–histopolation basis is used to ensure fast convergence. Detailed performance studies are carried out to analyze the efficiency of the GPU algorithms. The kernel throughput of each of the main algorithmic components is measured, and the strong and weak parallel scalability of the methods is demonstrated. The different relative weighting and significance of the algorithmic components on GPUs and CPUs is discussed. Results on problems involving adaptively refined nonconforming meshes are shown, and the use of the preconditioners on a large-scale magnetic diffusion problem using all spaces of the finite element de Rham complex is illustrated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Brittle failure analysis and modeling of high-burnup PWR fuel cladding alloys

The aim of this research is the development of methods for predicting mechanical behavior and identification of limiting conditions to prevent brittle failure of high-burnup (HBU) pressure water reactor (PWR) fuel cladding alloys. A finite element (FE) model of the ring compression test (RCT) was created to analyze the failure behavior of zirconium-based alloys with radial hydrides during the RCT. An elastic-plastic material model describes the zirconium alloy. The stress-strain curve needed for the elastic-plastic material model was derived by inverse finite element analyses. Cohesive zone modeling is used to reproduce sudden load drops during RCT loading. Based on the failure mechanism in non-irradiated ZIRLO (R) claddings, a micro-mechanical model was developed that distinguishes between brittle failure along hydrides and ductile failure of the zirconium matrix. Two different cohesive laws representing these types of failure are present in the same cohesive interface. The key differences between these constitutive laws are the cohesive strength, the stress at which damage initiates, and the cohesive energy, which is the damage energy dissipated by the cohesive zone. Statistically generated matrix-hydride distributions were mapped onto the cohesive elements and simulations with focus on the first load drop were performed. Computational results are in good agreement with the RCT results conducted on high-burnup M5 (R) samples. It could be shown that crack initiation and propagation strongly depend on the specific configuration of hydrides and matrix material in the fracture area.

Simbruner, Kai↗
DOE OSTI2021

Analysis of liquid petroleum using a laser-induced breakdown spectroscopy instrument

Here, a prototype analyzer for the direct LIBS analysis of nebulized liquid samples was developed and tested, particularly for the analysis of petroleum, organic solvents and aqueous solutions. The LIBS analyzer requires about 750 µl of liquid sample, 1 liter of N 2 gas, and 10 s of the analysis time to record 100 spectra. The limits of detection in oil and solvent are as low as 0.01–0.04 ppm for Li, Mg, and Cu. They increase for the difficult elements, such as Pb and Hg (7-10 ppm), Cl (250 ppm), and S (~0.7%). The relative standard deviation of measuring 100 ppm vanadium in oil and solvent was 1.5%. The LIBS detection limits and repeatability are better than required by the standard method ASTM D5185 for the analysis of lubricating oils in ICP-OES. Several petroleum samples were analyzed by LIBS and the quantitative results for V, Ni, and Fe compared to the ICP-OES data. Light crude oils can be nebulized and analyzed directly. Medium crude oils require minimal dilution at least 1:1, otherwise errors of determination become large. Presumably, utilization of the internal standard and chemometrics can be useful to correct for the matrix effects. In addition to the trace element analysis, the LIBS prototype demonstrated ability to measure the hydrogen-to-carbon ratio in organic liquid samples. Several unidentified features were observed in the carbon spectrum. Their possible origin is discussed.

02 PETROLEUM↗
DOE OSTI2021

Data compression and covariance matrix inspection: Cosmic shear

Covariance matrices are among the most difficult pieces of end-to-end cosmological analyses. In principle, for two-point functions, each component involves a four-point function, and the resulting covariance often has hundreds of thousands of elements. We investigate various compression mechanisms capable of vastly reducing the size of the covariance matrix in the context of cosmic shear statistics. This helps identify which of its parts are most crucial to parameter estimation. We start with simple compression methods, by isolating and “removing” 200 modes associated with the lowest eigenvalues, then those with the lowest signal-to-noise ratio, before moving on to more sophisticated schemes like compression at the tomographic level and, finally, with the massively optimized parameter estimation and data compression (MOPED). We find that, while most of these approaches prove useful for a few parameters of interest, like Ω m , the simplest yield a loss of constraining power on the intrinsic alignment (IA) parameters as well as S 8 . For the case considered—cosmic shear from the first year of data from the Dark Energy Survey—only MOPED was able to replicate the original constraints in the 16-parameter space. Finally, we apply a tolerance test to the elements of the compressed covariance matrix obtained with MOPED and confirm that the IA parameter A IA is the most susceptible to inaccuracies in the covariance matrix.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Low-Order Preconditioning for the High-Order Finite Element de Rham Complex

Here, we present a unified framework for constructing spectrally equivalent low-order-refined discretizations for the high-order finite element de Rham complex. This theory covers diffusion problems in H 1 , H(curl), and H(div) and is based on combining a low-order discretization posed on a refined mesh with a high-order basis for Nédélec and Raviart–Thomas elements that makes use of the concept of polynomial histopolation (polynomial fitting using prescribed mean values over certain regions). This spectral equivalence, coupled with algebraic multigrid methods constructed using the low-order discretization, results in highly scalable matrix-free preconditioners for high-order finite element problems in the full de Rham complex. Additionally, a new lowest-order (piecewise constant) preconditioner is developed for high-order interior penalty discontinuous Galerkin (DG) discretizations, for which spectral equivalence results and convergence proofs for algebraic multigrid methods are provided. In all cases, the spectral equivalence results are independent of polynomial degree and mesh size; for DG methods, they are also independent of the penalty parameter. These new solvers are flexible and easy to use; any “black-box” preconditioner for low-order problems can be used to create an effective and efficient preconditioner for the corresponding high-order problem. A number of numerical experiments are presented, based on an implementation in the finite element library MFEM. A range of challenging three-dimensional problems are used to corroborate the theoretical properties and demonstrate the flexibility and scalability of the method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Efficient exascale discretizations: High-order finite element methods

Efficient exploitation of exascale architectures requires rethinking of the numerical algorithms used in many large-scale applications. These architectures favor algorithms that expose ultra fine-grain parallelism and maximize the ratio of floating point operations to energy intensive data movement. One of the few viable approaches to achieve high efficiency in the area of PDE discretizations on unstructured grids is to use matrix-free/partially assembled high-order finite element methods, since these methods can increase the accuracy and/or lower the computational time due to reduced data motion. In this paper we provide an overview of the research and development activities in the Center for Efficient Exascale Discretizations (CEED), a co-design center in the Exascale Computing Project that is focused on the development of next-generation discretization software and algorithms to enable a wide range of finite element applications to run efficiently on future hardware. CEED is a research partnership involving more than 30 computational scientists from two US national labs and five universities, including members of the Nek5000, MFEM, MAGMA and PETSc projects. We discuss the CEED co-design activities based on targeted benchmarks, miniapps and discretization libraries and our work on performance optimizations for large-scale GPU architectures. We also provide a broad overview of research and development activities in areas such as unstructured adaptive mesh refinement algorithms, matrix-free linear solvers, high-order data visualization, and list examples of collaborations with several ECP and external applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Multi deep learning-based stochastic microstructure reconstruction and high-fidelity micromechanics simulation of time-dependent ceramic matrix composite response

A multi deep learning-based framework is developed for efficient, automated microstructure reconstruction and generation of stochastic representative volume elements (SRVEs) with periodic boundary conditions (PBCs) for accurate modeling of ceramic matrix composite (CMC) response. The methodology comprises a convolutional neural network coupled with regression layers to act as a vanilla regression network for semantic segmentation of the microstructure, allowing accurate characterization of the phases and their distributions at the microscale. Scanning electron microscope and confocal microscope are used to obtain C/SiNC and SiC/SiNC CMCs micrographs for vanilla regression testing. Microstructure variability in terms of fiber volume fraction and porosity are quantified through the output regression layer, ensuring accurate representation of material variability in SRVE construction. Generative adversarial network (GAN) and its variants are designed to produce high-fidelity SRVE, spanning CMCs microstructure variability space. A circular padding algorithm is developed to generate SRVEs with PBCs during training of GANs. The accuracy of the generated SRVEs is established through micromechanics simulations, where an efficient formulation of the high-fidelity generalized methods of cells (HFGMC) approach is used to compute the effective mechanical properties. Furthermore, an iterative algorithm is implemented in the HFGMC solver to simulate time-dependent deformation of SiC/SiNC subjected to creep loading conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effects of as-received defects on ceramic matrix composites properties using high-fidelity microstructures with periodic boundary conditions

The presence of microstructural defects in as-received specimens of ceramic matrix composites (CMCs) significantly influences their constitutive response and damage, highlighting the importance of characterization and quantification of these defects for accurate assessment of damage and failure in the service environment. In a recent effort, the authors developed an algorithm to generate stochastic representative volume elements (SRVEs) of Carbon fiber Silicon-Carbide-Nitride matrix (C/SiNC) CMCs based on extensive multiscale material and defect characterization data. This paper implements this algorithm within a commercial finite element solver with periodic boundary conditions (PBCs) for high-fidelity micromechanics analysis and investigation of macroscopic material behavior of C/SiNC composites. Different loading directions are used to predict the global mechanical properties, and the results are in excellent agreement with theoretical (rule of mixture) predictions. Subsequently, the effects of as-received defects on the global and local responses are investigated. The results show that intratow porosity has pronounced degradation effects on the global elastic properties and results in complex stress localization patterns, which can be attributed to potential damage initiation sites.

Khafagy, Khaled H.↗
DOE OSTI2025

Novel Synthesis and Characterization of Uranium-Zirconium Carbonitride by Direct Casting

UZrCN has exhibited thermophysical properties beneficial to high temperature reactor applications such as nuclear thermal rockets. The present work investigates a novel liquid phase synthesis method involving the admission of nitrogen gas during arc melting of uranium, zirconium, and carbon. Initial microstructural examinations using SEM indicated that the samples remain heterogenous with zirconium-rich cores in a uranium-rich matrix. Heterogeneity resulted from large differences in melting temperature of the major constituents and rapid solidification. Additional EDS, combustion analysis, and inert gas fusion analysis proved that the core regions have both uranium and zirconium and establish that carbon and nitrogen are retained during the fabrication process. XRD clarifies that the light elements did not form compounds with the uranium matrix but rather incorporated into the core region forming a sub-stoichiometric UZrCN. Heat treatment performed on equimolar U-Zr-C resulted in an increase in the homogenous phase present with an XRD pattern reflecting a ternary UZrCN despite small amounts of segregation in the final button.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

On Uranium-Zirconium Carbonitride Formation by Direct Casting - A Novel Synthesis Study

UZrCN has exhibited thermophysical properties beneficial to high temperature reactor applications such as nuclear thermal rockets. The present work investigates a novel liquid phase synthesis method involving the admission of nitrogen gas during arc melting of uranium, zirconium, and carbon. Initial microstructural examinations using SEM indicated that the samples remain heterogenous with zirconium-rich cores in a uranium-rich matrix. Heterogeneity resulted from large differences in melting temperature of the major constituents and rapid solidification. Additional EDS, combustion analysis, and inert gas fusion analysis proved that the core regions have both uranium and zirconium and establish that carbon and nitrogen are retained during the fabrication process. XRD clarifies that the light elements did not form compounds with the uranium matrix but rather incorporated into the core region forming a sub-stoichiometric UZrCN. Heat treatment performed on equimolar U-Zr-C resulted in an increase in the homogenous phase present with an XRD pattern reflecting a ternary UZrCN despite small amounts of segregation in the final button.

36 - MATERIALS SCIENCE↗