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 811 records · Page 45

Cosmology at high redshift — a probe of fundamental physics

An observational program focused on the high redshift (2<z<6) Universe has the opportunity to dramatically improve over upcoming LSS and CMB surveys on measurements of both the standard cosmological model and its extensions. Using a Fisher matrix formalism that builds upon recent advances in Lagrangian perturbation theory, we forecast constraints for future spectroscopic and 21-cm surveys on the standard cosmological model, curvature, neutrino mass, relativistic species, primordial features, primordial non-Gaussianity, dynamical dark energy, and gravitational slip. We compare these constraints with those achievable by current or near-future surveys such as DESI and Euclid, all under the same forecasting formalism, and compare our formalism with traditional linear methods. Here, our Python code FishLSS — used to calculate the Fisher information of the full shape power spectrum, CMB lensing, the cross-correlation of CMB lensing with galaxies, and combinations thereof — is publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

Emergent Superconductivity and Competing Charge Orders in Hole-Doped Square-Lattice t – J Model

The square-lattice Hubbard and closely related t-J models are considered as basic paradigms for understanding strong correlation effects and unconventional superconductivity (SC). Recent large-scale density matrix renormalization group simulations on the extended t-J model have identified d-wave SC on the electron-doped side (with the next-nearest-neighbor hopping t 2 > 0) but a dominant charge density wave (CDW) order on the hole-doped side (t 2 < 0), which is inconsistent with the SC of hole-doped cuprate compounds. We re-examine the ground-state phase diagram of the extended t-J model by employing the state-of-the-art density matrix renormalization group calculations with much enhanced bond dimensions, allowing more accurate determination of the ground state. On six-leg cylinders, while different CDW phases are identified on the hole-doped side for the doping range δ = 1/16 – 1/8, a SC phase emerges at a lower doping regime, with algebraically decaying pairing correlations and d-wave symmetry. On the wider eight-leg systems, the d-wave SC also emerges on the hole-doped side at the optimal 1=8 doping, demonstrating the winning of SC over CDW by increasing the system width. Furthermore, our results not only suggest a new path to SC in general t-J model through weakening the competing charge orders, but also provide a unified understanding on the SC of both hole- and electron-doped cuprate superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

How Bayesian methods can improve R -matrix analyses of data: The example of the d t reaction

The 3 H(d, n) 4 He reaction is of significant interest in nuclear astrophysics and nuclear applications. It is an important, early step in big-bang nucleosynthesis and a key process in nuclear fusion reactors. We use one- and two-level R-matrix approximations to analyze data on the cross section for this reaction at center-of-mass energies below 215 keV. We critically examine the data sets using a Bayesian statistical model that allows for both common-mode and additional point-to-point un- certainties. We use Markov Chain Monte Carlo sampling to evaluate this R-matrix-plus-statistical model and find two-level R-matrix results that are stable with respect to variations in the channel radii. The S factor at 40 keV evaluates to 25.36(19) MeV b (68% credibility interval). We discuss our Bayesian analysis in detail and provide guidance for future applications of Bayesian methods to R-matrix analyses. We also discuss possible paths to further reduction of the S-factor uncertainty.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling of Hydraulic Fracture Propagation at the kISMET Site Using a Fully Coupled 3D Network-Flow and Quasi- Static Discrete Element Model

Aimed at supporting the design of hydraulic fracturing experiments at the kISMET site, ~1500 m below ground in a deep mine, we performed pre-experimental hydraulic fracturing simulations in order to estimate the breakdown pressure, propagation pressure, fracture geometry, and the magnitude of induced seismicity using a newly developed fully coupled three-dimensional (3D) network flow and quasi-static discrete element model (DEM). The quasi-static DEM model, which is constructed by Delaunay tessellation of the rock volume, considers rock fabric heterogeneities by using the “disordered” DEM mesh and adding random perturbations to the stiffness and tensile/shear strengths of individual DEM elements and the elastic beams between them. A conjugate 3D flow network based on the DEM lattice is constructed to calculate the fluid flow in both the fracture and porous matrix. One distinctive advantage of the model is that fracturing is naturally described by the breakage of elastic beams between DEM elements. It is also extremely convenient to introduce mechanical anisotropy into the model by simply assigning orientation-dependent tensile/shear strengths to the elastic beams. In this paper, the 3D hydraulic fracturing model was verified against the analytic solution for a penny-shaped crack model. We applied the model to simulate fracture propagation from a vertical open borehole based on initial estimates of rock mechanical properties and in-situ stress conditions. The breakdown pressure and propagation pressure are directly obtained from the simulation. In addition, the released elastic strain energies of individual fracturing events were calculated and used as a conservative estimate for the magnitudes of the potential induced seismic activities associated with fracturing. The comparisons between model predictions and experimental results are still ongoing.

15 GEOTHERMAL ENERGY↗

Reflectance properties of one-dimensional metal-dielectric ternary photonic crystal

Metallic photonic crystal has a very important application in absorption enhancement in solar cells. It has been found that an ultra-thin metallic layer becomes transparent due to internal scattering of light through the each interface of the dielectric and metal surfaces. The metal has absorption due to their surface plasmon and the plasmon has important parameters for changing optical properties of the metal. We consider ternary metallic-dielectric photonic crystal (MDPC) for having large probabilities to change the optical properties of the MDPC and the photonic crystals may be changed by changing dimensionality, symmetry, lattice parameters, Filling fraction and effective refractive index refractive index contrast. In this present communication, we try to show that the photonic band gap in ternary metal-dielectric photonic crystal can be significantly enlarged when air dielectric constant is considered. All the theoretical analyses are made based on the transfer matrix method together with the Drude model of metal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Draft Experimental Plan for Commercial SNF Degradation in Repository Environments with a Focus on Fracture Matrix Degradation

This report documents work performed supporting US Department of Energy (DOE) Nuclear Energy Spent Fuel and Waste Disposition, Spent Fuel and Waste Science and Technology, under work breakdown structure element 1.08.01.03.04, “Inventory and Waste Form Characteristics and Performance.” In particular, this report fulfills the M3 milestone, M3SF-21OR010309072 "Draft Experimental Plan for Commercial SNF Degradation in Repository Environments” as described in work package SF-21OR01030907 "SNF Degradation Testing-ORNL".This experimental plan is a collaboration between ORNL, SNL, ANL, and PNNL. It is intended to support source term model development and spent fuel degradation testing in a range of generic repository conditions. Potential testing and experiments to support validation of the Fuel Matrix Dissolution Model (FMDM) and the overall development/validation of repository source term modeling are discussed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Matrix Multiply Performance of GPUs on Exascale-class HPE/Cray Systems

The computation of dense matrix-matrix products (GEMMs) is central to many modeling and simulation workloads as well as AI/ML deep learning campaigns. In fact, millions of dollars are spent annually on computing GEMMs, and large model training demands are increasing exponentially. Specialized processors such as GPUs are designed to perform well for these operations. However, the performance of GEMMs on GPUs can exhibit complex behaviors depending on many factors, making it challenging to optimize the performance of GEMMs on these processors. In this study we undertake an examination of GEMM performance on several leading GPU models taken from product lines of GPUs to be deployed in forthcoming exascale computing systems. We show results to illustrate the many factors that can affect performance of GEMMs on GPUs. We then present data collected from a large number of test runs for an example GEMM operation to show the dependence behaviors of GEMM rate on matrix dimensions. Finally, we show results from machine learning-based performance models using novel feature engineering methods to fit the measured performance, providing a potential basis for GEMM performance tuning and autotuning methods for GPUs. Recommendations are also given for how to achieve high GEMM performance on modern GPUs.

Melesse Vergara, Veronica↗

Atomistic model of xenon gas bubble re-solution rate due to thermal spike in uranium oxide

Atomistic simulations are performed to study the response of Xe gas bubbles in UO 2 to ionizing fission products through the thermal spike approximation. A portion of the total electronic stopping power (S e ) is taken as the thermal spike energy through a ratio variable ζ. The thermal spike energy causes extreme melting within the fission track cylindrical region. Molecular dynamics is employed to quantify the probability of a Xe gas atom to be re-solved (re-dissolved) back into the UO 2 matrix. Subsequently, a re-solution model is developed and parametrized as a function of bubble radius (R), off-centered distance (r), and thermal spike energy (ζS e ). The off-centered distance measures the shift of the thermal spike axis from the bubble center. To evaluate the re-solution model, independent fission product yield of U-235 fission due to thermal neutrons (0.0253 eV), taken from the JEFF-3.3 database, is used. The kinetic energy of the fission products is taken from the EXFOR database. Subsequently, the decay of Se over distance for each fission product is simulated. Finally, the evaluated re-solution rate (re-solution probability per second) is presented as a function of bubble radius for a range of ζ.

74 ATOMIC AND MOLECULAR PHYSICS↗

Thermal power spectrum in the CFT driven cosmology

We present an overview of the recently suggested cosmological model driven by conformal field theory (CFT) with the initial conditions in the form of the microcanonical density matrix. In particular, we discuss the origin of inflationary stage in this model and a novel feature — the thermal nature of the primordial power spectrum of the CMB anisotropy. The relevant effect of ''temperature of the relic temperature anisotropy'' can be responsible for a thermal contribution to the red tilt of this spectrum, additional to its conventional vacuum component. The amplification of this effect due to recently established a-theorem in CFT is briefly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

How quantum phases on cylinders approach the two-dimensional limit

We consider the properties of $T$ = $0$ quantum phases of matter - especially superconducting and analogous spin-liquid phases - on infinite cylinders of width $L$ ⊥ and analyze the ways in which the $L$ ⊥ → ∞ (two-dimensional, 2d) limit is approached. This problem is interesting in its own right but is particularly important in the context of extrapolating accessible density matrix renormalization group (DMRG) results on model strongly interacting problems to the desired 2d limit. In conclusion, various methods for drawing firm conclusions about the quantum phases in 2d from relatively small $L$ ⊥ results are illustrated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Accelerating Transformer-based Deep Learning Models on FPGAs using Column Balanced Block Pruning

Although Transformer-based language representations achieve state-of-the-art accuracy on various natural language processing (NLP) tasks, the large model size has been challenging the resource constrained computing platforms. Weight pruning, as a popular and effective technique in reducing the number of weight parameters and accelerating the Transformer, has been investigated on GPUs. However, the Transformer acceleration using weight pruning on field-programmable gate array (FPGAs) remains unexplored. This paper investigates the column balanced block-wise pruning on Transformer and designs an FPGA acceleration engine to customize the balanced blockwise matrix multiplication. We implement the Transformer model with proper hardware scheduling, and the experiments show that the Transformer inference on FPGA achieves 10.35 ms latency with the batch size of 32, which is 10.96 × speed up comparing to CPU platform and 2.08 × speed up comparing to GPU platform.

Peng, Hongwu↗

Transport properties in dilute UN ( X ) solid solutions ( X = Xe , Kr )

Uranium nitride (UN) is a candidate fuel for current GEN III fission reactors, for which it is investigated as an accident-tolerant fuel, as well as for future GEN IV reactors. In this study, we investigate the kinetic properties of gas fission products (Xe and Kr) in UN. Binding and migration energies are obtained using density functional theory, with an added Hubbard correlation to model f electrons, and the occupation matrix control scheme to avoid metastable states. These energies are then used as input for the self-consistent mean field method which enables to determine transport coefficients for vacancy-mediated diffusion of Xe and Kr on the U sublattice. The magnetic ordering of the UN structure is explicitly taken into account, for both energetic and transport properties. Solute diffusivities are compared with experimental measurements and the effect of various parameters on the theoretical model is carefully investigated. We find that kinetic correlations are very strong in this system, and that despite atomic migration anisotropy, macroscopic solute diffusivities show limited anisotropy. Our model indicates that the discrepancy between experimental measurements probably results from different irradiation conditions, and hence different defect concentrations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Efficacy of the SU(3) scheme for ab initio large-scale calculations beyond the lightest nuclei

We report on the computational characteristics of ab initio nuclear structure calculations in a symmetry-adapted no-core shell model (SA-NCSM) framework. We examine the computational complexity of the current implementation of the SA-NCSM approach, dubbed LSU3shell, by analyzing ab initio results for 6 Li and 12 C in large harmonic oscillator model spaces and SU3-selected subspaces. We demonstrate LSU3shell’s strong-scaling properties achieved with highly-parallel methods for computing the many-body matrix elements. Results compare favorably with complete model space calculations and significant memory savings are achieved in physically important applications. In particular, a well-chosen symmetry-adapted basis affords memory savings in calculations of states with a fixed total angular momentum in large model spaces while exactly preserving translational invariance.

74 ATOMIC AND MOLECULAR PHYSICS↗

Peculiarities of the formation of photonic nanojets by a matrix of dielectric microtoroids

We report the results of theoretical modelling of photonic nanojets (PNJs) formed by laser radiation scattering on a single-layer ordered assembly of dielectric microtoroids placed on the surface of a transparent matrix (silicone film). Using the method of computational electrodynamics (FDTD), the main PNJ parameters (length, width, peak intensity) are analysed under the conditions of mutual influence of the light fields of neighbouring microparticles. It is shown that the main factor affecting the PNJ characteristics under study is the spatial configuration of the radiation-scattering particle, namely, the internal diameter of its cross section. It is found that for certain configurations of toroid placement in a cluster, a PNJ ensemble is implemented with parameters significantly better than those for a single toroid. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Efficacy of the SU(3) scheme for ab initio large-scale calculations beyond the lightest nuclei

We report on the computational characteristics of ab initio nuclear structure calculations in a symmetry-adapted no-core shell model (SA-NCSM) framework. We examine the computational complexity of the current implementation of the SA-NCSM approach, dubbed LSU3shell, by analyzing ab initio results for 6 Li and 12 C in large harmonic oscillator model spaces and SU(3)-selected subspaces. We demonstrate LSU3shell's strong-scaling properties achieved with highly-parallel methods for computing the many-body matrix elements. Results compare favorably with complete model space calculations and signi cant memory savings are achieved in physically important applications. In particular, a well-chosen symmetry-adapted basis a ords memory savings in calculations of states with a fixed total angular momentum in large model spaces while exactly preserving translational invariance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Detection Of Special Nuclear Materials Tagged Neutrons

In order to detect Special Nuclear Materials (SNM) in unattended luggage or cargo containers in the field of homeland security, fissions are induced by 14 MeV neutrons produced by an associated particle DT neutron generator, and prompt fission particles correlated with tagged neutron are detected by plastic scintillators. SMN produce high multiplicity events due to induced fissions, whereas nonnuclear materials produce low multiplicity events due to cross-talk, (n,2n) or (n,n'γ) reactions. The data acquisition electronics is made of compact FPGA boards. The coincidence window is triggered by the alpha particle detection, allowing to tag the emission date and direction of the 14 MeV interrogating neutron. The first part of the paper presents experiment vs. calculation comparisons to validate MCNP-PoliMi simulations and the post-processing tools developed with the data analysis framework ROOT. Measurements have been performed using different targets (iron, lead, graphite), first with small plastic scintillators (10 x 10 x 10 cm{sup 3}) and then with large detectors (10 x 10 x 100 cm{sup 3}) to demonstrate that nuclear materials can be differentiated from nonnuclear dense materials (iron, lead) in iron and wood matrixes. A special attention is paid on SNM detection in abandoned luggage. In the second part of the paper, the performances of a cargo container inspection system are studied by numerical simulation, following previous work reported in. Detectors dimensions and shielding against the neutron generator background are optimized for container inspection. Events not correlated to an alpha particle (uncorrelated background), counting statistics, time and energy resolutions of the data acquisition system are all taken into account in a realistic numerical model. The impact of the container matrix (iron, ceramic, wood) has been investigated by studying the system capability to detect a few kilograms of SNM in different positions in the cargo container, within 10 min acquisitions. (authors)

36 MATERIALS SCIENCE↗

Uncertainties in nuclear transition matrix elements of neutrinoless ββ decay

To estimate the uncertainties associated with the nuclear transition matrix elements M{sup (K)} (K=0ν/0N) for the 0{sup +} → 0{sup +} transitions of electron and positron emitting modes of the neutrinoless ββ decay, a statistical analysis has been performed by calculating sets of eight (twelve) different nuclear transition matrix elements M{sup (K)} in the PHFB model by employing four different parameterizations of a Hamiltonian with pairing plus multipolar effective two-body interaction and two (three) different parameterizations of Jastrow short range correlations. The averages in conjunction with their standard deviations provide an estimate of the uncertainties associated the nuclear transition matrix elements M{sup (K)} calculated within the PHFB model, the maximum of which turn out to be 13% and 19% owing to the exchange of light and heavy Majorana neutrinos, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗