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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Electron-impact excitation data for W 2+ in support of tungsten spectroscopy and re-deposition measurements for magnetically-confined plasmas

Abstract To better understand plasma wall interactions involving tungsten, accurate atomic structure and electron-impact driven collisional processes for near-neutral ion stages of tungsten are required. Complementing existing work on neutral and singly ionised tungsten, atomic structure and collisional calculations for W 2+ electron-impact excitation have been completed. These excitation calculations are an important component of S/XB coefficients for near-neutral charge states, which may be used to spectroscopically infer re-deposition of tungsten at the plasma-solid boundary of fusion relevant devices. With W 2+ in particular having emission lines that can be observed at ultraviolet (UV) wavelengths, while higher charge states of tungsten are unlikely to have lines possible to observe outside of the vacuum UV range. The atomic structure was generated using the General-purpose Relativistic Atomic Structure Package (GRASP 0 ), implementing the Multi-configuration Dirac Fock approach. This structure was the basis for a subsequent Dirac R -matrix electron-impact excitation calculation to provide Maxwellian averaged rate coefficients. A synthetic spectrum was generated from this data using a collisional-radiative model to predict the strongest W III spectral lines and these lines were compared to emission from the Compact Toroidal Hybrid (CTH) plasma device. Several of the strongest W III lines are observed in CTH and agree well with the modelled line wavelengths and intensities, a table of these lines is provided that could be observed in other devices.

McCann, M. (ORCID:0000000215321240)↗
DOE OSTI2020

Measurement of the Li 7 ( γ , t ) He 4 ground-state cross section between E γ = 4.4 and 10 MeV

Here, the Li 7 ( γ , t ) He 4 ground state cross section was measured for the first time using monoenergetic γ rays with energies between 4.4 and 10 MeV at the High Intensity Gamma-ray Source. The reaction is important for the primordial Li problem and for testing our understanding of the mirror α -capture reactions H 3 ( α , γ ) Li 7 and He 3 ( α , γ ) Be 7 . Although over the last 30 years most measurements of the H 3 ( α , γ ) Li 7 reaction have concentrated in an energy range below E γ = 3.65 MeV, measurements at higher energies could potentially restrict the extrapolation to astrophysically important energies. The experimental arrangement for measuring the Li 7 ( γ , t ) He 4 reaction included a large-area silicon detector array and several beam characterization instruments. The experimental astrophysical S factor of H 3 ( α , γ ) calculated from the present data was fitted using the R -matrix formalism. The results are in disagreement with previous experimental measurements in the same energy range but the extrapolated S factor agrees with the potential model calculation and lower energy experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Hyper Spectral Anomaly Detection

The HSA is a statistics based anomaly detection model. The model performs unsupervised anomaly detection, based on a datapoint's density and similarity within a dataset. Density and similarity data are encoded into an affinity matrix. The affinity matrix is evolved to summarize the data's structure on greater topographical scales within the data's function space. The set of evolved affinity matrices and an anomaly score vector are passed to a user defined penalized objective function. The penalized objective function of anomaly scores is then minimized. Data points where the absolute value of the z-scores of anomaly scores greater than a specified threshold are predicted as anomalies. A novel multi-filter feature has also been implemented. To reduce false positive rates, the multi-filter records the indexes of the HSA predictions. A new dataset and data loader are instantiated consisting of all the initial HSA predictions and non-anomalous data points in a 10% and 90% split respectively. The HSA is then run through this data set and a count of number of times a data point is predicted is kept. In this way the initial predictions may be compared with data spanning the entire dataset. After the multi-filter is complete, all datapoints will have an associated anomaly score, as well as a multi-filter prediction count to further filter the anomalous predictions.

Rogers, DempseyD [Idaho National Laboratory (INL),↗
DOE OSTI2024

Geometry-aware training of factorized layers in tensor Tucker format

Reducing parameter redundancies in neural network architectures is crucial for achieving feasible computational and memory requirements during train and inference of large networks. Given its easy implementation and flexibility, one promising approach is layer factorization, which reshapes weight tensors into a matrix format and parameterizes it as the product of two rank-r matrices. However, this family of approaches often requires an initial full-model warm-up phase, prior knowledge of a feasible rank, and it is sensitive to parameter initialization.In this work, we introduce a novel approach to train the factors of a Tucker decomposition of the weight tensors. Our training proposal proves to be optimal in locally approximating the original unfactorized dynamics and stable for the initialization. Furthermore, the rank of each mode is dynamically updated during training.We provide a theoretical analysis of the algorithm, showing convergence, approximation and local descent guarantees. The method's performance is further illustrated through a variety of experiments, showing remarkable training compression rates and comparable or even better performance than the full baseline and alternative layer factorization strategies.

Zangrando, Emanuele [Gran Sasso Science Institute ↗
DOE OSTI2025

Exploring the energy landscape of RBMs: reciprocal space insights into bosons, hierarchical learning and symmetry breaking

Deep generative models have become ubiquitous due to their ability to learn and sample from complex distributions. Despite the proliferation of various frameworks, the relationships among these models remain largely unexplored, a gap that hinders the development of a unified theory of AI learning. In this work, we address two central challenges: clarifying the connections between different deep generative models and deepening our understanding of their learning mechanisms. We focus on Restricted Boltzmann Machines (RBMs), a class of generative models known for their universal approximation capabilities for discrete distributions. By introducing a reciprocal space formulation for RBMs, we reveal a connection between these models, diffusion processes, and systems of coupled bosons. Our analysis shows that at initialization, the RBM operates at a saddle point, where the local curvature is determined by the singular values of the weight matrix, whose distribution follows the Marc̆enko-Pastur law and exhibits rotational symmetry. During training, this rotational symmetry is broken due to hierarchical learning, where different degrees of freedom progressively capture features at multiple levels of abstraction. This leads to a symmetry breaking in the energy landscape, reminiscent of Landau’s theory. This symmetry breaking in the energy landscape is characterized by the singular values and the weight matrix eigenvector matrix. We derive the corresponding free energy in a mean-field approximation. We show that in the limit of infinite size RBM, the reciprocal variables are Gaussian distributed. Our findings indicate that in this regime, there will be some modes for which the diffusion process will not converge to the Boltzmann distribution. To illustrate our results, we trained replicas of RBMs with different hidden layer sizes using the MNIST dataset. Our findings not only bridge the gap between disparate generative frameworks but also shed light on the fundamental processes underpinning learning in deep generative models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A novel in vivo mouse intervertebral disc degeneration model induced by compressive suture

Intervertebral disc degeneration (IDD) is the root cause of many musculoskeletal disorders of the spine. However, the etiology of IDD is complex and still not well understood. Animal models of IDD would be useful in deciphering the underlying mechanisms. But the existing animal models have their limitations. Therefore, to establish a novel mouse model that can simulate the human IDD process in vivo, we proposed to carefully circumcise the 2 mm-wide tail skin and then compressively sutured the defect with a simple end-to-end suture to exert excessive pressure on the disc. After 1-week, 2-week, and 4-week compression, the mice were sacrificed and the intervertebral discs were harvested for tissue analysis. The radiological, morphological, and molecular modifications of intervertebral discs were measured to characterize this model. Radiologically, the water content of the intervertebral disc decreased significantly after 2-week compression. Morphologically, the nucleus pulposus showed a decrease in volume and the number of notochordal cells. The compressive suture also broke the balance between anabolic and catabolic enzymes in nucleus pulposus, which led to the remodeling of the extracellular matrix in nucleus pulposus as the content of aggrecan and collagen II decreased. The compressive suture could induce intervertebral discs degeneration in a more reasonable way, which was solely influenced by mechanical loading, as the mice caudal vertebrae still moved freely after the operation. This kind of animal model could be adapted as a reliable in vivo mouse IDD model for the research regarding the etiology and treatments of IDD.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Design and Evaluation of Environmental Barrier Coatings for Protection of Ceramic Matrix Composites in Hydrogen-Based Turbines

This presentation at the 2024 FECM Spring Project Review Meeting summarizes ongoing work at NETL aimed at enabling the use of ceramic matrix composite (CMC) materials in next-generation high-efficiency power generation turbines fueled by hydrogen, natural gas, or other fuels. The work broadly involves three aspects (1) phase-field damage modeling of CMCs and the environmental barrier coatings (EBCs) used to protect them (2) computational design of new and improved EBCs and (3) performance testing of EBCs using a novel gas turbine combustion simulation rig developed at NETL. An overview on current status and next steps for each aspect of the project are presented.

ceramic matrix composite (CMC)↗
DOE OSTI2020

Construction of a Code Verification Matrix for Heat Conduction With Finite Element Code Applications

When establishing the pedigree of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model. During this process, numerical results on various meshes are systematically compared to a reference analytic solution. The selection of analytic solutions can be a laborious process, as it is difficult to establish adequate code confidence without performing redundant work. In this report we address this issue by applying a physics-based process that establishes a set of reference problems. In this process, code simulation options are categorized and systematically tested, which ensures that gaps in testing are easily identified and addressed. The resulting problems are primarily intended for code verification analysis but may also be useful for comparison to other simulation codes, troubleshooting activities, or training exercises. The process is used to select fifteen code verification problems relevant for the one-dimensional steady-state heat conduction equation. These problems are applicable to a wide variety of simulation tools, but, in this work, a demonstration is performed using the finite element-based nuclear fuel performance code BISON. Convergence to the analytic solution at the theoretical rate is quantified for a selection of the problems, which establishes a baseline pedigree for the code. Not only can this standard set of conduction solutions be used for verification of other codes, but also the physics-based process for selecting problems can be utilized to quantify and expand testing for any simulation tool.

42 ENGINEERING↗
DOE OSTI2024

Physics-informed Estimation of the Covariance Matrix for Various Neutron Spectra

A method for estimating covariance matrices which capture the uncertainties in calculated reactor spectra has been developed. This method is based on perturbing the parameters of a physics-based analytic model fitted to a calculated spectrum. The covariance of the perturbed analytic spectra imposes energy-dependent correlations due to the physics of the neutron processes in the reactor, i.e., a fission component, a 1/E down-scatting component, and a thermal Maxwellian component. An analytic model is developed which is shown to produce good fits to several reactor environments. The covariance matrices produced via this method are then used as the prior spectrum in STAYSL least squares spectrum adjustment where it is combined with integral metrics, such as activation measurements, to produce a high-fidelity neutron spectrum characterization. It was concluded that the methodology showed agreeable results for the ACRR free-field spectrum adjustment in STAYSL resulting in a 𝜒 2 value of 2.21 (per degree of freedom), but further work is needed to describe scattering and interface regions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Hydrogen Uptake Kinetics of 1,4-Bis(phenylethynyl)benzene Rubberized O-Rings: Measurements, Modeling, and Comparison with Additional Forms of Organic Getters

A class of molecules called “hydrogen getters” can react with, or scavenge, H 2 in applications where the hydrogen presence and/or buildup are not desirable. One such “getter”, 1,4-bis(phenylethynyl)benzene or DEB, can be incorporated into a silicone matrix in the form of an O-ring for added flexibility, environmental resiliency, and convenient use as a gasket in sealed applications. However, the performance and kinetics of this DEB-loaded rubberized O-ring have not yet been characterized. Here in this work, the hydrogen uptake kinetics of the rubberized DEB O-ring were extracted by the isoconversional analysis from isothermal isobaric data under conditions of 13,332 Pa H 2 and 305–325 K. The isoconversional and cylindrical diffusion approximations were then used to predict and model the hydrogen uptake of rubberized DEB O-rings under any arbitrary condition, as illustrated in this work for the simple case of a constant rate of hydrogen generation/input. In comparison with other Pd/carbon-based/organic getter systems, these rubberized DEB O-rings provide a type of hydrogen getter material with enhanced flexibility and catalyst protection in applications where an O-ring seal is needed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Uncertainty quantification in elastic constants of SiC f /SiC m tubular composites using global sensitivity analysis

Silicon carbide fiber and silicon carbide matrix (SiC f /SiC m ) tubes produced through the chemical vapor infiltration process have become a candidate cladding material in nuclear applications. The performance of this composite is influenced by many variables such as braiding angle, porosity, material properties, etc., which vary over a range of values due to the inherent fluctuations in the manufacturing process. In this study, the variability in elastic constants of SiC f /SiC m composite has been quantified through multiscale finite element (FE) simulations, variable screening, and high-fidelity surrogate modeling. The key variables dominantly affecting the elastic constants of SiC f /SiC m tubes were identified using global sensitivity analysis. A surrogate to the high-fidelity FE-based model was used in Monte Carlo simulations to generate a hundred thousand samples from which the uncertainty in elastic constants was assessed. It turned out that the coefficient of variation was less than 10%.

Materials Science↗
DOE OSTI2024

A principle of maximum ignorance for semiclassical gravity

The principle of maximum ignorance posits that the coarse-grained description of a system is maximally agnostic about its underlying microscopic structure. We briefly review this principle for random matrix theory and for the eigenstate thermalization hypothesis. We then apply this principle in holography to construct ensembles of random mixed states. This leads to an ensemble of microstates which models our microscopic ignorance, and which on average reproduces the effective semiclassical physics of a given bulk state. We call this ensemble the state-averaging ansatz. The output of our model is a prediction for semiclassical contributions to variances and higher statistical moments over the ensemble of microstates. The statistical moments provide coarse-grained — yet gravitationally non-perturbative — information about the microstructure of the individual states of the ensemble. We show that these contributions exactly match the on-shell action of known wormhole configurations of the gravitational path integral. These results strengthen the view that wormholes simply parametrize the ignorance of the microstructure of a fundamental state, given a fixed semiclassical bulk description.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Non-standard neutrino interactions mediated by a light scalar at DUNE

Abstract We investigate the effect on neutrino oscillations generated by beyond-the-standard-model interactions between neutrinos and matter. Specifically, we focus on scalar-mediated non-standard interactions (NSI) whose impact fundamentally differs from that of vector-mediated NSI. Scalar NSI contribute as corrections to the neutrino mass matrix rather than the matter potential and thereby predict distinct phenomenology from the vector-mediated ones. Similar to vector-type NSI, the presence of scalar-mediated neutrino NSI can influence measurements of oscillation parameters in long-baseline neutrino oscillation experiments, with a notable impact on CP measurement in the case of DUNE. Our study focuses on the effect of scalar NSI on neutrino oscillations, using DUNE as an example. We introduce a model-independent parameterization procedure that enables the examination of the impact of all non-zero scalar NSI parameters simultaneously. Subsequently, we convert DUNE’s sensitivity to the NSI parameters into projected sensitivity concerning the parameters of a light scalar model. We compare these results with existing non-oscillation probes. Our findings reveal that the region of the light scalar parameter space sensitive to DUNE is predominantly excluded by non-oscillation probes, especially when considering all nonzero parameters simultaneously for DUNE.

Physics↗
DOE OSTI2026

Enhanced HLW glass property-composition models - phase 3

During the present phase (Phase 3) of work to enhance and expand the HLW glass property-composition models, test data for 137 glasses were collected and incorporated into the combined WTP/ORP database. The new data include those collected on glasses from two statistically designed matrices to augment the moderate alumina region (the HLW16-MA matrix covering the region of 9 wt% < Al2O3 < 12 wt% with 16 glasses) and the ultra-high alumina region (the HLW16-UHA matrix covering the region of 26 wt% < Al2O3 < 29 wt% with 29 glasses). The addition of ORP glasses over the different phases of testing has significantly expanded the combined WTP/ORP database. The compiled data include glass compositions, PCT (B, Li and Na) releases, spinel T1%, electrical conductivity, viscosity, TCLP-Cd releases, and nepheline formation upon CCC. With the exception of the CCC spinel formation data – the CCC data will be used to support the further development of nepheline models for Hanford HLW glasses.

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

Anyon superconductivity from topological criticality in a Hofstadter–Hubbard model

We argue that the combination of strong repulsive interactions and high magnetic fields can generate electron pairing and superconductivity. Inspired by the large lattice constants of moiré materials, which make large flux per unit cell accessible at laboratory fields, we study the triangular lattice Hofstadter-Hubbard model at one-quarter flux quantum per plaquette, where previous literature has argued that a chiral spin liquid separates a weak-coupling integer quantum Hall phase and a strong-coupling topologically trivial antiferromagnetic insulator at a density of one electron per site. We argue that topological superconductivity emerges upon doping in the vicinity of the integer quantum Hall to chiral spin liquid transition. We employ exact diagonalization and density matrix renormalization group methods to examine this theoretical scenario and find that electronic pairing indeed occurs on both sides of criticality over a remarkably broad range of interaction strengths. On the chiral spin liquid side, our results provide a concrete model realization of the long-hypothesized mechanism of anyon superconductivity. Our study thus establishes a beyond-Bardeen-Cooper-Schrieffer route to electron pairing in a well-controlled limit, relying crucially on the interplay between electron correlations and band topology.

Hubbard model↗
DOE OSTI2021

Learning constitutive relations using symmetric positive definite neural networks

In this work, we present a new neural-network architecture, called the Cholesky-factored symmetric positive definite neural network (SPD-NN), for modeling constitutive relations in computational mechanics. Instead of directly predicting the stress of the material, the SPD-NN trains a neural network to predict the Cholesky factor of the tangent stiffness matrix, based on which the stress is calculated in incremental form. As a result of this special structure, SPD-NN weakly imposes convexity on the strain energy function, satisfies the second order work criterion (Hill's criterion) and time consistency for path-dependent materials, and therefore improves numerical stability, especially when the SPD-NN is used in finite element simulations. Depending on the types of available data, we propose two training methods, namely direct training for strain and stress pairs and indirect training for loads and displacement pairs. We demonstrate the effectiveness of SPD-NN on hyperelastic, elasto-plastic, and multiscale fiber-reinforced plate problems from solid mechanics. The generality and robustness of SPD-NN make it a promising tool for a wide range of constitutive modeling applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Colloid-facilitated transport of 238 Pu, 233 U and 137 Cs through fractured chalk: Laboratory experiments, modelling, and implications for nuclear waste disposal

The influence of montmorillonite colloids on the mobility of 238 Pu, 233 U and 137 Cs through a chalk fracture was investigated to assess the transport potential for radioactive waste. Radioisotopes of each element, along with the conservative tracer tritium, were injected in the presence and absence of montmorillonite colloids into a naturally fractured chalk core. In parallel, batch experiments were conducted to obtain experimental sorption coefficients (K d , mL/g) for both montmorillonite colloids and the chalk fracture material. Breakthrough curves were modelled to determine diffusivity and sorption of each radionuclide to the chalk and the colloids under advective conditions. Uranium sorbed sparingly to chalk (log K d = 0.7 ± 0.2) in batch sorption experiments. 233 U(VI) breakthrough was controlled primarily by the matrix diffusion and sorption to chalk (15 and 25% recovery with and without colloids, respectively). Cesium, in contrast, sorbed strongly to both the montmorillonite colloids and chalk (batch log K d = 3.2 ± 0.01 and 3.9 ± 0.01, respectively). The high affinity to chalk and low colloid concentrations overwhelmed any colloidal Cs transport, resulting in very low 137 Cs breakthrough (1.1–5.5% mass recovery). Batch and fracture transport results, and the associated modelling revealed that Pu migrates both as Pu (IV) sorbed to montmorillonite colloids and as dissolved Pu(V) (7% recovery). Transport experiments revealed differences in Pu(IV) and Pu(V) transport behavior that could not be quantified in simple batch experiments but are critical to effectively predict transport behavior of redox-sensitive radionuclides. Finally, a brackish groundwater solution was injected after completion of the fracture flow experiments and resulted in remobilization and recovery of 2.2% of the total sorbed radionuclides which remained in the core from previous experiments. In general, our study demonstrates consistency in sorption behavior between batch and advective fracture transport. Finally, the results suggest that colloid-facilitated radionuclide transport will enhance radionuclide migration in fractured chalk for those radionuclides with exceedingly high affinity for colloids.

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

Exploring Hilbert space on a budget: Novel benchmark set and performance metric for testing electronic structure methods in the regime of strong correlation

This work explores the ability of classical electronic structure methods to efficiently represent (compress) the information content of full configuration interaction (FCI) wave functions. We introduce a benchmark set of four hydrogen model systems of different dimensionalities and distinctive electronic structures: a 1D chain, a 1D ring, a 2D triangular lattice, and a 3D close-packed pyramid. To assess the ability of a computational method to produce accurate and compact wave functions, we introduce the accuracy volume, a metric that measures the number of variational parameters necessary to achieve a target energy error. Using this metric and the hydrogen models, we examine the performance of three classical deterministic methods: (i) selected configuration interaction (sCI) realized both via an a posteriori (ap-sCI) and variational selection of the most important determinants, (ii) an a posteriori singular value decomposition (SVD) of the FCI tensor (SVD-FCI), and (iii) the matrix product state representation obtained via the density matrix renormalization group (DMRG). We find that the DMRG generally gives the most efficient wave function representation for all systems, particularly in the 1D chain with a localized basis. For the 2D and 3D systems, all methods (except DMRG) perform best with a delocalized basis, and the efficiency of sCI and SVD-FCI is closer to that of DMRG. For larger analogs of the models, the DMRG consistently requires the fewest parameters but still scales exponentially in 2D and 3D systems, and the performance of SVD-FCI is essentially equivalent to that of ap-sCI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗