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 1,315 records · Page 73

Charge Resolution of the Silicon Matrix of the ATIC Experiment

ATIC (Advanced Thin Ionization Calorimeter) is a balloon borne experiment designed to measure the cosmic ray composition for elements from hydrogen to iron and their energy spectra from approx.50 GeV to near 100 TeV. It consists of a Si-matrix detector to determine the charge of a CRT particle, a scintillator hodoscope for tracking, carbon interaction targets and a fully active BGO calorimeter. ATIC had its first flight from McMurdo, Antarctica from 28/12/2000 to 13/01/2001. The ATIC flight collected approximately 25 million events. The silicon matrix of the ATIC spectrometer is designed to resolve individual elements from proton to iron. To provide this resolution careful calibration of each pixel of the silicon matrix is required. Firstly, for each electronic channel of the matrix the pedestal value was subtracted taking into account its drift during the flight. The muon calibration made before the flight was used then to convert electric signals (in ADC channel number) to energy deposits in each pixel. However, the preflight muon calibration was not accurate enough for the purpose, because of lack of statistics in each pixel. To improve charge resolution the correction was done for the position of Helium peak in each pixel during the flight . The other way to set electric signals in electronics channels of the Si-matrix to one scale was correction for electric channel gains accurately measured in laboratory. In these measurements it was found that small different nonlinearities for different channels are present in the region of charge Z > 20. The correction for these non-linearities was not done yet. In linear approximation the method provides practically the same resolution as muon calibration plus He-peak correction. For searching a pixel with the signal of primary particle an indication from the cascade in the calorimeter was used. For this purpose a trajectory was reconstructed using weight centers of energy deposits in BGO layers. The point of intersection of this trajectory with Si-matrix and its RMS was determined. The pixel with maximal signal in 3sigma region was taken as sought. The signal in this pixel was corrected by trajectory zenith angle. The preliminary results on charge resolution of the Si-matrix in the range from protons to iron are presented.

Zatsepin, V. I.↗

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↗

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.↗

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↗

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↗

Highly parallel sparse Cholesky factorization

Several fine grained parallel algorithms were developed and compared to compute the Cholesky factorization of a sparse matrix. The experimental implementations are on the Connection Machine, a distributed memory SIMD machine whose programming model conceptually supplies one processor per data element. In contrast to special purpose algorithms in which the matrix structure conforms to the connection structure of the machine, the focus is on matrices with arbitrary sparsity structure. The most promising algorithm is one whose inner loop performs several dense factorizations simultaneously on a 2-D grid of processors. Virtually any massively parallel dense factorization algorithm can be used as the key subroutine. The sparse code attains execution rates comparable to those of the dense subroutine. Although at present architectural limitations prevent the dense factorization from realizing its potential efficiency, it is concluded that a regular data parallel architecture can be used efficiently to solve arbitrarily structured sparse problems. A performance model is also presented and it is used to analyze the algorithms.

Gilbert, John R.↗

On Dynamics of Spinning Structures

This paper provides details of developments pertaining to vibration analysis of gyroscopic systems, that involves a finite element structural discretization followed by the solution of the resulting matrix eigenvalue problem by a progressive, accelerated simultaneous iteration technique. Thus Coriolis, centrifugal and geometrical stiffness matrices are derived for shell and line elements, followed by the eigensolution details as well as solution of representative problems that demonstrates the efficacy of the currently developed numerical procedures and tools.

Rotating bodies↗

Predicting ductility of Mg/SiCp nanocomposite under multiaxial loading conditions based on unit cell modeling

This article presents an investigation on the ductility of Mg/SiCp under nine different multiaxial loading conditions based on a unit cell (or called representative volumetric element, RVE) finite element model. The constitutive relationship of magnesium matrix is characterized by a model combined Lode angle dependent plasticity model with the CPB06 plasticity model, which is able to describe the tension-compression asymmetry in magnesium matrix. The debonding of the matrix-particle interface and matrix fracture are considered by the modified Mohr-Coulomb (MMC) fracture model. Simulation results show that the interface debonding and subsequent matrix fracture are responsible for the ductility and fracture of nanocomposite. The computational results with parametric studies are compared and well correlated with experimental data, and it finally gives a stress-based full fracture locus of Mg/SiCp nanocomposite.

42 ENGINEERING↗

Transformational Challenge Reactor Safety Design and Radionuclide Retention Strategy

An integrated safety design and radionuclide (RN) retention strategy is developed to support the Transformational Challenge Reactor (TCR) demonstration. This demonstration aims to showcase viability for rapid deployment of a novel reactor by leveraging the advances in materials, manufacturing, and computational sciences through a highly integrated and agile design and development approach. This strategy provides a logical description and understanding of how RNs are contained within the facility. Rather than discussing fission product barriers individually between separate design and safety basis reports, this paper provides a consistent description and narrative to better facilitate regulatory interactions and focus safety design efforts. The principal barriers credited include the various coating layers in the tristructural isotropic (TRISO) fuel particle, the silicon carbide (SiC) matrix hosting the particles within the fuel element, the helium pressure boundary, and the confinement system. The choice and assumed performance of the credited barriers are highly conservative, which is a direct reflection of the low hazard that the TCR demonstration presents and the need to simplify and focus the safety review process accordingly. However, the strategy and the associated framework are generalized and may be adopted and tailored to support other advanced reactor demonstration efforts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Fast calculation of the tokamak vertical instability

There has been recent interest in fast calculations of the tokamak axisymmetric vertical instability for real time feedback control purposes. It is shown that the maximum eigenvalue for the basic rigid version of this stability problem can be obtained by finding the positive root to a simple scalar function. This function can be generalized to include plasma mass and has complexity linear in the number of conductive elements. The formulation is based on standard matrix decompositions of the fixed-geometry part of the eigenproblem. The calculation bottleneck is the summary of mutual inductances from the reconstructed equilibrium current density. The with-mass spectrum can be made fully real-valued by the addition of a critical amount of damping with negligible effect on the vertical growth rate. Furthermore, the calculation has been implemented in the plasma control system at the DIII-D tokamak and used in experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Precise β branching-ratio measurement for the 0 + → 0 + superallowed decay of 34 Ar

We have measured the branching ratio for the superallowed 0 + → 0 + β transition from 34 Ar to be 0.9448(8), and determined its ft value to be 3058.1(28) s, a result with ±0.09% precision, which is a factor of three improvement over the previous result based on current world data. The ft-value ratio for the mirror pair of superallowed transitions 34 Ar → 34 Cl and 34 Cl → 34 S, becomes the most precise yet measured and, in a sensitive test of the method used to calculate the isospin-symmetry-breaking correction, δ C , it agrees well with the ratio as calculated with Woods-Saxon radial wave functions. This confirms the method used in the most recent survey of superallowed decays to extract V ud , the up-down quark-mixing element of the Cabibbo-Kobayashi-Maskawa matrix. Additionally, our branching-ratio results for the four observed Gamow-Teller branches to 1 + states in 34 Cl are shown to agree well with shell-model calculations based on the same effective interactions that were used in the determination of δ C .

39 ≤ A ≤ 58↗

Superallowed 0 + → 0 + nuclear β decays: 2020 critical survey, with implications for V ud and CKM unitarity

A new critical survey of all half-life, decay-energy and branching-ratio measurements related to 23 superallowed 0 + → 0 + β decays is presented. Included are 222 individual measurements of comparable precision obtained from 174 published references. Compared with our last survey in 2015, we have added results from 28 new publications and eliminated an approximately equal number whose results have been superseded by much more precise modern data. We obtain world-average ft values for each of the 21 transitions that have a complete set of data, then apply radiative and isospin-symmetry-breaking corrections to extract “corrected” Ft values. Fifteen of these Ft values now have a precision of 0.3% or better and all take the same value within statistics, as expected from conservation of the vector current. Their average, Ft¯, when combined with the muon lifetime, yields the up-down quark-mixing element of the Cabibbo-Kobayashi-Maskawa matrix, V ud = 0.97373 ± 0.00031. This is lower than our 2015 result by one standard deviation and its uncertainty is increased by 50%. This is a consequence, not of any shifts in the experimental data, but of new calculations for the radiative corrections. The lower V ud value now leads to greater tension in the top-row test of unitarity in the CKM matrix. Updates in experimental data have independently led to a factor-of-two tighter limit being set on the possible existence of a scalar interaction. In conclusion, the new limit on Fierz interference is b F ≤ 0.0033 at the 90% confidence level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Novel approaches to determine B ± and B 0 meson production fractions

We propose novel methods to determine the ϒ ( 4 S ) → B + B − and ϒ ( 4 S ) → B 0 B ¯ 0 decay rates. The precision to which they and their ratio are known yields at present a limiting uncertainty around 2% in measurements of absolute B decay rates, and thus in a variety of applications, such as precision determinations of elements of the Cabibbo-Kobayashi-Maskawa matrix and flavor-symmetry relations. The new methods we propose are based in one case on exploiting the ϒ ( 5 S ) datasets, in the other case on the different average number of charged tracks in B ± and B 0 decays. We estimate future sensitivities using these methods and discuss possible measurements of f d / f u at the (HL-)LHC. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Tacho

SAND2022-7470 O Tacho is a performance portable sparse direct solver for use with various sparse matrix problems which typically arise from finite element methods. The code relies on the Kokkos parallel programming model porting to both CPUs and GPUs. It is also part of the Trilinos library. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Kim, Kyungjoo↗

Sparse matrix methods based on orthogonality and conjugacy

A matrix having a high percentage of zero elements is called spares. In the solution of systems of linear equations or linear least squares problems involving large sparse matrices, significant saving of computer cost can be achieved by taking advantage of the sparsity. The conjugate gradient algorithm and a set of related algorithms are described.

Lawson, C. L.↗

An algorithm for reducing the bandwidth and profile of a sparse matrix

A new algorithm for reducing the bandwidth and profile of a sparse matrix is described. Extensive testing on finite element matrices indicates that the algorithm typically produces bandwidth and profile which are comparable to those of the commonly-used reverse Cuthill-McKee algorithm, yet requires significantly less computation time.

Gibbs, N. E.↗

Time dependent difference theory for sound propagation in axisymmetric ducts with plug flow

The time dependent governing acoustic-difference equations and boundary conditions are developed and solved for sound propagation in an axisymmetric (cylindrical) hard wall duct with a plug mean flow and spinning acoustic modes. The analysis begins with a harmonic sound source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the transient as well as the 'steady' state solutions of the acoustic field. The time dependent finite difference analysis has two advantages over the steady state finite difference and finite element techniques: (1)the elimination of large matrix storage requirements, and (2)shorter solution times under most conditions.

Baumeister, K. J.↗