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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Recent development of Monte Carlo shell model and its application to no-core calculations

One of the major challenges in nuclear theory is to reproduce and to predict nuclear structure from ab initio calculations with realistic nuclear forces. As the current limitation of direct diagonalization of Hamiltonian matrices by Lanczos iteration method is around the order of matrix dimensionality 10 10 in shell-model calculations, it is difficult to access heavier nuclei beyond the p shell with sufficiently large basis spaces. It is possible to overcome this difficulty by utilizing efficient approximate methods to reproduce full ab initio solutions with good precision and quantified uncertainties. Following the major success of the Monte Carlo shell model (MCSM) with an assumed inert core in the sd- and pf-shell regions and also by recent developments in the MCSM algorithm, the no-core MCSM is expected to be one of the most powerful tools to meet these conditions. We have performed benchmark calculations in the p-shell region. Results of energies are compared with those in the full configuration interaction and no-core full configuration methods. These are found to be consistent with each other within quoted uncertainties when they could be quantified. We also compare and discuss the radial density of the helium-4 ground state extracted from the MCSM and FCI many-body wave functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

ThinCurr: An open-source 3D thin-wall eddy current modeling code for the analysis of large-scale systems of conducting structures

In this paper we present a new thin-wall eddy current modeling code, ThinCurr, for studying inductively-coupled currents in 3D conducting structures -- with primary application focused on the interaction between currents flowing in coils, plasma, and conducting structures of magnetically-confined plasma devices. The code utilizes a boundary finite element method on an unstructured, triangular grid to accurately capture device structures. The new code, part of the broader Open FUSION Toolkit, is open-source and designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. Scalability to large models is enabled through use of hierarchical off-diagonal low-rank compression of the inductance matrix, which is otherwise dense. Ease of handling large models of complicated geometry is further supported by automatic determination of supplemental elements through a greedy homology approach. Here, a detailed description of the numerical methods of the code and verification of the implementation of those methods using cross-code comparisons against the VALEN code and Ansys commercial analysis software is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

A Subbasin-based framework to represent land surface processes in an Earth System Model

Realistically representing spatial heterogeneity and lateral land surface processes within and between modeling units in earth system models is important because of their implications to surface energy and water exchange. The traditional approach of using regular grids as computational units in land surface models and earth system models may lead to inadequate representation of lateral movements of water, energy and carbon fluxes, especially when the grid resolution increases. Here a new subbasin-based framework is introduced in the Community Land Model (CLM), which is the land component of the Community Earth System Model (CESM). Local processes are represented assuming each subbasin as a grid cell on a pseudo grid matrix with no significant modifications to the existing CLM modeling structure. Lateral routing of water within and between subbasins is simulated with the subbasin version of a recently-developed physically based routing model, Model for Scale Adaptive River Routing (MOSART). As an illustration, this new framework is implemented in the topographically diverse region of the U.S. Pacific Northwest. The modeling units (subbasins) are delineated from high-resolution Digital Elevation Model while atmospheric forcing and surface parameters are remapped from the corresponding high resolution datasets. The impacts of this representation on simulating hydrologic processes are explored by comparing it with the default (grid-based) CLM representation. In addition, the effects of DEM resolution on parameterizing topography and the subsequent effects on runoff processes are investigated. Limited model evaluation and comparison showed that small difference between the averaged forcing can lead to more significant difference in the simulated runoff and streamflow because of nonlinear horizontal processes. Topographic indices derived from high resolution DEM may not improve the overall water balance, but affect the partitioning between surface and subsurface runoff. More systematic analyses are needed to determine the relative merits of the subbasin representation compared to the commonly used grid-based representation, especially when land surface models are approaching higher resolutions.

Subbasin-based representation, Community Land Mode↗
DOE OSTI2018

System and method for analyzing oscillatory stability in electrical power transmission systems

A computer-based method for contingency analysis of oscillatory stability in an electrical power transmission system is provided. The method uses at least one processor. The method includes receiving, by the at least one processor, a plurality of component inputs from a plurality of system components within the electrical power transmission system. The method also includes generating a nominal matrix for the electrical power transmission system. The nominal matrix includes a set of equations at least partially modeling the electrical power transmission system. The method further includes calculating eigenvalues and eigenvectors of the nominal matrix. The method also includes identifying a contingency representing a postulated disturbance of the electrical power transmission system. The method further includes estimating a contingency eigenvalue for the contingency using the eigenvalues and eigenvectors of the nominal matrix.

Baone, Chaitanya Ashok↗
DOE OSTI2014

Calculations of single crystal elastic constants for yttria partially stabilised zirconia from powder diffraction data

Yttria Stabilised Zirconia (YSZ) is a tough, phase-transforming ceramic that finds use in a wide range of commercial applications from dental prostheses to thermal barrier coatings. Micromechanical modelling of phase transformation can deliver reliable predictions in terms of the influence of temperature and stress. However, models must rely on the accurate knowledge of single crystal elastic stiffness constants. Some techniques for elastic stiffness determination are well-established. The most popular of these involve exploiting frequency shifts and phase velocities of acoustic waves. However, the application of these techniques to YSZ can be problematic due to the micro-twinning observed in larger crystals. Here, we propose an alternative approach based on selective elastic strain sampling (e.g., by diffraction) of grain ensembles sharing certain orientation, and the prediction of the same quantities by polycrystalline modelling, for example, the Reuss or Voigt average. The inverse problem arises consisting of adjusting the single crystal stiffness matrix to match the polycrystal predictions to observations. In the present model-matching study, we sought to determine the single crystal stiffness matrix of tetragonal YSZ using the results of time-of-flight neutron diffraction obtained from an in situ compression experiment and Finite Element modelling of the deformation of polycrystalline tetragonal YSZ. The best match between the model predictions and observations was obtained for the optimized stiffness values of C11 = 451, C33 = 302, C44 = 39, C66 = 82, C12 = 240, and C13 = 50 (units: GPa). Considering the significant amount of scatter in the published literature data, our result appears reasonably consistent.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A new woven composite constitutive model validated by shock wave experiments

In this paper, we present results of plate impact simulations of shock compressed woven glass fiber-reinforced plastic (GRP) performed using the Arbitrary Lagrangian–Eulerian three-dimensional finite element code. A hyperelastic large-strain-based empirical Continuum Damage Mechanics (CDM) formulation is employed to describe damage initiation and growth in the shock-compressed GRP. The model parameters calibration scheme utilizes the Velocity Interferometer System for Any Reflector normal particle velocity measurements at the free surface of the GRP target plates. The impact velocity in the experiments ranged from 8.5 to 418 m/s. The finite element model considered planar 0°/90° bidirectional plies with an individual ply thickness of 0.68 mm, stacked to reach a total laminate thickness of 6.8 mm. The anisotropic elastic strains were estimated from the experimentally determined tetragonal symmetry stiffness matrix for the GRP. The strain-based damage model captures several salient features observed in the measured free surface particle wave profiles, including the shock rise time, onset of Elastic—Elastic Cracking, and the shape of the nonlinear portion of the experimental particle velocity profiles. Furthermore, the CDM model predicts the dominant damage mode to be matrix microcracking due to shear and the associated bulk expansion (bulking) under the global compressive loading in the plate impact configuration.

42 ENGINEERING↗
DOE OSTI2021

Isospin mixing and the cubic isobaric multiplet mass equation in the lowest T=2, A=32 quintet

The isobaric multiplet mass equation (IMME) is known to break down in the first T=2,A=32 isospin quintet. In this work we combine high-resolution experimental data with state-of-the-art shell-model calculations to investigate isospin mixing as a possible cause for this violation. The experimental data are used to validate isospin-mixing matrix elements calculated with newly developed shell-model Hamiltonians. Our analysis shows that isospin mixing with nonanalog T=1 states contributes to the IMME breakdown, making the requirement of an anomalous cubic term inevitable for the multiplet.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Connecting dark matter annihilation to the vertex functions of Standard Model fermions

We consider here scenarios in which dark matter is a Majorana fermion which couples to Standard Model fermions through the exchange of charged mediating particles. The matrix elements for various dark matter annihilation processes are then related to one-loop corrections to the fermion-photon vertex, where dark matter and the charged mediators run in the loop. In particular, in the limit where Standard Model fermion helicity mixing is suppressed, the cross section for dark matter annihilation to various final states is related to corrections to the Standard Model fermion charge form factor. These corrections can be extracted in a gauge-invariant manner from collider cross sections. While current measurements from colliders are not precise enough to provide useful constraints on dark matter annihilation, improved measurements at future experiments, such as the International Linear Collider, could improve these constraints by several orders of magnitude, allowing them to surpass the limits obtainable by direct observation

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

TAUOLA update for decay channels with e + e – pairs in the final state

With the arrival of high luminosity B-factories like the Belle II experiment, τ decay measurements have become more precise than ever, allowing rarer processes to be explored, and finer details of τ decays to be studied. These are important to understand the spectrum of intermediate particles produced in τ decays. Therefore Monte Carlo generators, like the TAUOLA program, have to facilitate precision analysis as well as confront new models that constantly emerge with the availability of high statistics experimental data. New decay channels and models may lead to large variation of matrix elements size within the available phase space, as in the case when e + e - pairs are present in final state. It requires appropriate presampler of the phase space generator and a proper documentation to help users introduce their own models. While releasing a new update, it is important for the TAUOLA Monte Carlo library to maintain the general structure of the previous versions to preserve backward compatibility. The aim is to minimize changes from the user perspective. Finally, this paper presents a demonstration of new models implementation facilitated by the current update.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

SAM Enhancements and Model Developments for Molten-Salt-Fueled Reactors

To support the development and utilization of the SAM code for molten-salt-fueled reactor (MSR) safety analysis and licensing, an effort was devoted to enhancing code capabilities and developing reference models for the MSR primary loop. A reference standard problem of a prototypical reactor design is foundational to NRC to verify the adequacy of computer codes and evaluation models for a specific reactor type. A thermo-fluid model of the Molten Salt Reactor Experiment (MSRE) has been developed. The MSRE primary loop model consists of a 2-D core region and external core components in 1-D or 0-D. Both the steady-state and a transient scenario are simulated. The porous medium model is utilized for both the molten salt fluid channels and the moderator matrix in the MSRE core. The use of the porous medium model for the MSRE core is verified with the higher-fidelity simulation results using 1-D representations of the fluid channels and 3-D modeling of core structures. To further enhance SAM multi-scale simulation capabilities for MSRs, a mass transport model is developed and implemented in SAM multi-dimension flow module for modeling of species transport such as delayed neutron precursors in MSRs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Is the Matrix Completion of Reduced Density Matrices Unique?

Reduced density matrices are central to describing observables in many-body quantum systems. In electronic structure theory, the two-particle reduced density matrix (2-RDM) suffices to determine the energy and other key properties. Recent work has used matrix completion, leveraging the low-rank structure of RDMs and approximate theoretical models, to reconstruct the 2-RDM from partial data and thus reduce the computational cost. However, matrix completion is, in general, an under-determined problem. Revisiting Rosina’s theorem (Rosina, M. Queen’s Papers on Pure and Applied Mathematics , 1968, No. 11, 369), we here show that the matrix completion is unique under certain conditions, identifying the subset of 2-RDM elements that enables its exact reconstruction from incomplete information. Building on this, we introduce a hybrid quantum–stochastic algorithm that achieves exact matrix completion, demonstrated through applications to the Fermi–Hubbard model.

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

Massive and massless spin-2 scattering and asymptotic superluminality

We constrain theories of a massive spin-2 particle coupled to a massless spin-2 particle by demanding the absence of a time advance in eikonal scattering. This is an S-matrix consideration that leads to model-independent constraints on the cubic vertices present in the theory. Of the possible cubic vertices for the two spin-2 particles, the requirement of subluminality leaves a particular linear combination of cubic vertices of the Einstein-Hilbert type. Either the cubic vertices must appear in this combination or new physics must enter at a scale parametrically the same as the mass of the massive spin-2 field, modulo some standard caveats. These conclusions imply that there is a one-parameter family of ghost-free bimetric theories of gravity that are consistent with subluminal scattering. When both particles couple to additional matter, subluminality places additional constraints on the matter couplings. We lastly reproduce these constraints by considering classical scattering off of a shockwave background in the ghost-free bimetric theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Machine-Learning-based Algorithms for Automated Image Segmentation Techniques of Transmission X-ray Microscopy (TXM)

Four state-of-the-art Deep Learning-based Convolutional Neural Networks (DCNN) were applied to automate the semantic segmentation of a 3D Transmission x-ray Microscopy (TXM) nanotomography image data. The standard U-Net architecture as baseline along with UNet++, PSPNet, and DeepLab v3+ networks were trained to segment the microstructural features of an AA7075 micropillar. A workflow was established to evaluate and compare the DCNN prediction dataset with the manually segmented features using the Intersection of Union (IoU) scores, time of training, confusion matrix, and visual assessment. Comparing all model segmentation accuracy metrics, it was found that using pre-trained models as a backbone along with appropriate training encoder-decoder architecture of the Unet++ can robustly handle large volumes of x-ray radiographic images in a reasonable amount of time. This opens a new window for handling accurate and efficient image segmentation of in situ time-dependent 4D x-ray microscopy experimental datasets.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A comparative study of discrete fracture network and equivalent continuum models for simulating flow and transport in the far field of a hypothetical nuclear waste repository in crystalline host rock

One of the major challenges of simulating flow and transport in the far field of a geologic repository in crystalline host rock is related to reproducing the properties of the fracture network over the large volume of rock with sparse fracture characterization data. Various approaches have been developed to simulate flow and transport through the fractured rock. The approaches can be broadly divided into Discrete Fracture Network (DFN) and Equivalent Continuum Model (ECM). The DFN explicitly represents individual fractures, while the ECM uses fracture properties to determine equivalent continuum parameters. In this paper, we compare DFN and ECM in terms of upscaled observed transport properties through generic fracture networks. The major effort was directed on making the DFN and ECM approaches similar in their conceptual representations. This allows for separating differences related to the interpretation of the test conditions and parameters from the differences between the DFN and ECM approaches. The two models are compared using a benchmark test problem that is constructed to represent the far field (1 × 1 × 1 km 3 ) of a hypothetical repository in fractured crystalline rock. The test problem setting uses generic fracture properties that can be expected in crystalline rocks. The models are compared in terms of the: 1) effective permeability of the domain, and 2) nonreactive solute breakthrough curves through the domain. The principal differences between the models are mesh size, network connectivity, matrix diffusion and anisotropy. We demonstrate how these differences affect the flow and transport. Finally, we identify the factors that should be taken in consideration when selecting an approach most suitable for the site-specific conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2017

Stability of nanosized oxides in ferrite under extremely high dose self ion irradiations

We produced a nanostructured ferritic alloy (NFA), 14YWT, in the form of thin walled tubing. The stability of the nano-oxides (NOs) was determined under 3.5 MeV Fe +2 irradiations up to a dose of ~585 dpa at 450 °C. Transmission electron microscopy (TEM) and atom probe tomography (APT) show that severe ion irradiation results in a ~25% reduction in size between the unirradiated and irradiated case at 270 dpa while no further reduction within the experimental error was seen at higher doses. Conversely, number density increased by ~30% after irradiation. Moreover, this ‘inverse coarsening’ can be rationalized by the competition between radiation driven ballistic dissolution and diffusional NO reformation. There were no significant changes in the composition of the matrix or NOs observed after irradiation. Modeling the experimental results also indicated a dissolution of the particles.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Isothermal annealing behaviour of nuclear grade 20Cr-25Ni austenitic stainless steel

Highlights: • At 930 °C recrystallisation of deformed 20/25 Nb steel is completed after 120 s. • Recrystallisation proceeds by strain induced boundary migration. • An increase of cube oriented grains at the expense of S oriented is observed. • The boundary migration induced the precipitation of niobium carbides. • The grain size is limited by the particle dispersion according to Zener pinning. - Abstract: We have performed an in-depth characterisation of the microstructure evolution of 20Cr-25Ni Nb-stabilised austenitic stainless steel during isothermal annealing at 930 °C using scanning and transmission electron microscopy. This steel grade is used as cladding material in advanced gas-cooled fission reactors, due to its resistance to thermal creep and water corrosion. The initial deformed microstructure undergoes recrystallisation via a strain-induced boundary migration mechanism, attaining a fully recrystallised microstructure after 120 s of annealing. The transition from low-to-high grain boundaries has already occurred after 15 s, together with an increase in the cube grain orientation at the expense of the S texture component. After 120 s, the grain boundary migration induces the formation of new fine Nb(C,N) particles, whereas the pre-existing particles become enriched in Ni and Si. The resulting particle population limits the grain growth in the austenitic matrix, based on the Zener pinning model, resulting in relatively small recrystallised austenite grains and a high density of high-angle and special coincidence-lattice-site grain boundaries, together with a large number of particle/matrix interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multiscale and multiphysics FEA simulation and materials optimization for laser ultrasound transducers

In this study, the relationship between the nanocomposite design and the laser ultrasound transducer (LUT) characteristics was investigated through simulations in multiple scale levels for material behavior, device response, and acoustic wave propagation in media. First, the effects of the nanoparticle size and concentration on the effective properties of composites were quantitatively investigated with the finite element analysis (FEA) method. Second, the effective properties of the nanocomposite were assigned to the layer, which is modeled as a homogeneous material, in the FEA for the LUT simulating the energy conversion from the incident laser to the acoustic wave. Finally, the ultrasound propagation in the water was calculated by a theoretical wave propagation model. The FEA-based prediction was compared with the experimental data in the literature and a theoretical analysis for LUT based on Thermal-Acoustic coupling. As a result, the ultrasound waves on the transducer surface and at a distance in the water could be predicted. Based on the hierarchically integrated prediction procedure, the optimal conditions of the photoacoustic nanocomposites were investigated through the parametric study with the particle size and concentration as variables. The results guide the material designs optimized for different device characteristics, such as high pressure and broad bandwidth.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Experimental and theoretical studies of laser-induced breakdown spectroscopy emission from iron oxide: Studies of atmospheric effects

In this paper, we report on a comprehensive study of the emission spectra from laser-induced breakdown spectroscopy (LIBS) measurements on iron oxide. Measurements have been made of the emission from Fe 2 O 3 under atmospheres of air, He, and Ar, and at different atmospheric pressures. The effect of varying the time delay of the measurement is also explored. Theoretical calculations were performed to analyze the plasma conditions and find that a reasonably consistent picture of the change in plasma temperature and density for different atmospheric conditions can be reached. We also investigate the sensitivity of the OI 777 nm emission lines to the plasma conditions, something that has not been explored in detail in the previous work. In conclusion, we also show that LIBS can be used to differentiate between FeO and Fe 2 O 3 by examining the ratio of the intensities of selected Fe emission to O emission lines.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗