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 541 records · Page 30

Research in Theoretical Particle Physics: Neutrino Masses and The Origin of CP-Violation

The main goal of this project was to investigate the theories of neutrino masses at the low scale and the origin of CP-violation in physics beyond the Standard Model. The origin of neutrino masses is one of the most pressing issues in particle physics. In this proposal, we investigated the theories for neutrino masses based on local total lepton number. In these theories anomaly cancellation predicts the existence of extra fermions with lepton number, and one of them is a good dark matter (DM) candidate. The cosmological constraints on the DM relic density implies that the lepton number symmetry breaking scale must be below the multi-TeV scale. Therefore, one can hope to test the origin of neutrino masses in current or future experiments. We investigated in great detail the anomaly cancellation in these gauge theories, study the different mechanisms for neutrino masses, study the predictions for direct and indirect dark matter experiments. Since these theories predict new sources of CP violation, we investigated the predictions for the electric dipole moments (EDM). The Higgs decays and the different signatures at the Large Hadron Collider were investigated in great detail. Finally, we studied the possible baryogenesis mechanisms in these theories in agreement with the EDM, DM and collider constraints. The origin of CP violation in the Standard Model is unknown. In this proposal, we investigated the simplest mechanisms for spontaneous CP-violation to explain the CP-violation in the CKM matrix and the value of the QCD vacuum angle. We discussed the Nelson-Barr mechanism in gauge theories predicting vector-like quarks from anomaly cancellation such as theories for local baryon number. We discussed the Bento-Branco-Parada mechanism in gauge theories to explain the CP-Violation in the CKM matrix. We will investigate models with two Higgs doublets in the context of the minimal theory for quark-lepton unification. We will show the non-decoupling effects in the Higgs sector, the main constraints coming from flavour violating processes when the Yukawa couplings are related by the gauge symmetry in theories for quark-lepton unification. We will investigate the predictions for electric dipole moments in these theories. The possibility to have successful baryogenesis was investigated. Finally, we investigated the relation between CP-violation in the quark and leptonic sectors. The future results from these studies could help us to understand two main issues in physics beyond the Standard Model: The Origin of Neutrino Masses and CP-violation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Materials structure–property factorization for identification of synergistic phase interactions in complex solar fuels photoanodes

Abstract Properties can be tailored by tuning composition in high-order composition spaces. For spaces with complex phase behavior, modeling the properties as a function of composition and phase distribution remains a formidable challenge. We present materials structure–property factorization (MSPF) as an approach to automate modeling of such data and identify synergistic phase interactions. MSPF is an interpretable machine learning algorithm that couples phase mapping via Deep Reasoning Networks (DRNets) to matrix factorization-based modeling of the representative properties of each phase in a dataset. MSPF is demonstrated for Bi–Cu–V oxide photoanodes for solar fuel generation, which contains 25 different phase combinations and correspondingly exhibits complex composition-structure-photoactivity relationships. Comparing the measured photoactivity to a learned model for non-interacting phases, synergistic phase interactions are identified to guide further photoactivity optimization and understanding. MSPF identifies synergistic interactions of a BiVO 4 -like phase with both Cu 2 V 2 O 7 -like and CuV 2 O 6 -like phases, creating avenues for understanding complex photoelectrocatalysts.

36 MATERIALS SCIENCE↗

Efficient machine-learning model for fast assessment of elastic properties of high-entropy alloys

We combined descriptor-based analytical models for stiffness-matrix and elastic-moduli with mean-field methods to accelerate assessment of technologically useful properties of high-entropy alloys, such as strength and ductility. Model training for elastic properties uses Sure-Independence Screening (SIS) and Sparsifying Operator (SO) method yielding an optimal analytical model, constructed with meaningful atomic features to predict target properties. Computationally inexpensive analytical descriptors were trained using a database of elastic properties determined from density functional theory for binary and ternary subsets of Nb-Mo-Ta-W-V refractory alloys. The optimal Elastic-SISSO models, extracted from an exponentially large feature space, give an extremely accurate prediction of target properties, similar to or better than other models, with some verified from existing experiments. Here we also show that electronegativity variance and elastic-moduli can directly predict trends in ductility and yield strength of refractory HEAs, and reveals promising alloy concentration regions.

36 MATERIALS SCIENCE↗

Dual-Channel Enzymatic Inhibition Measurement (DEIM) Coupling Isotope Substrate via Matrix-Assisted Laser Desorption/Ionization Time of Flight Mass Spectrometry

A novel dual-channel enzymatic inhibition measurement (DEIM) method was developed to improve the repeatability with light/heavy isotope substrates, producing reliable relative standard deviations (< 3%) by employing acetylcholinesterase (AChE) as the model enzyme. The matrix-assisted laser desorption/ionization time of flight mass spectrometry (MALDI-TOF MS) was adapted for enzyme-inhibited method due to its good salt-tolerance and high throughput; meanwhile, dual-channel enzymatic reactions were performed to improve the repeatability of each well. The acetylcholinesterase inhibition measurement was conducted by mixing the quenched enzyme reaction solution of blank group (with heavy isotope as substrate) and experimental group (with light isotope as substrate), of which the inhibition rate might be affected by isotope effects. Hence, inverse study and Km measurement were implemented to validate the method. The inverse study shows similar inhibition rate (68.9 and 70.3%) and the Km of isotope substrates are analogous (0.139 and 0.135 mM), which demonstrated that the novel method is feasible to AChE inhibition measurement. Finally, the method was applied to herb extracts, half of which exhibit inhibition to AChE. The precise dual-channel enzymatic inhibition measurement (DEIM) method could be regarded as a promising approach to potential enzyme inhibitor screening. .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Scaling up the transcorrelated density matrix renormalization group

Explicitly correlated methods, such as the transcorrelated method which shifts a Jastrow or Gutzwiller correlator from the wave function to the Hamiltonian, are designed for high-accuracy calculations of electronic structures, but their application to larger systems has been hampered by the computational cost. We develop improved techniques for the transcorrelated density-matrix renormalization group (DMRG), in which the ground state of the transcorrelated Hamiltonian is represented as a matrix product state (MPS), and demonstrate large-scale calculations of the ground-state energy of the two-dimensional Fermi-Hubbard model. Our developments stem from three technical inventions: (i) constructing matrix product operators (MPOs) of transcorrelated Hamiltonians with low bond dimension and high sparsity, (ii) exploiting the entanglement structure of the ground states to increase the accuracy of the MPS representation, and (iii) optimizing the nonlinear parameter of the Gutzwiller correlator to mitigate the nonvariational nature of the transcorrelated method. Here, we examine systems of size up to 12×12 lattice sites, four times larger than previous transcorrelated DMRG studies, and demonstrate that transcorrelated DMRG yields significant improvements over standard nontranscorrelated DMRG for equivalent computational effort. Transcorrelated DMRG reduces the error of the ground-state energy by 2.4×–14×, with the smallest improvement seen for a small system at half filling and the largest improvement in a dilute closed-shell system.

Density matrix renormalization group↗

Crystallographic orientation-dependent strain hardening in a precipitation-strengthened Al-Cu alloy

While the strengthening of Al-Cu alloys due to precipitation has been extensively studied, the effect of crystallographic orientation of the matrix and precipitates, as well as precipitate morphology, on the strain hardening behavior is not well understood. In this paper, we investigate this effect with in situ neutron diffraction during deformation of an Al-Cu alloy (206) after multiple aging treatments. Precipitate-dislocation interactions were found to change from precipitate shearing for microstructures predominantly containing GPI and θ" precipitates to Orowan looping for microstructures with primarily θ' and θ precipitates. Notably, significant anisotropy in strain hardening behavior was observed when θ' precipitates were present, which was attributed to crystallographic orientation dependent load transfer from the Al matrix to the θ' precipitates. The anisotropic load transfer is hypothesized to be caused by the extent of rotation of high aspect-ratio θ' precipitates, owing to dislocations looping around them during plastic deformation of the matrix. Predictions from an analytical model describing the anisotropic magnitude of load transfer from precipitate rotation agree well with experimental results, successfully validating the precipitate rotation hypothesis and explaining the anisotropic strain hardening behavior. This model allows for the prediction of stresses separately in the precipitate and matrix phases as a function of crystallographic orientation, only given the bulk mechanical properties.

36 MATERIALS SCIENCE↗

Analytic error analysis of cross section interpolation methods in nodal diffusion codes - I : Theory

This paper discusses two cross section interpolation methods commonly found in popular nodal codes; the partial derivatives and multiple tables models. The motivation for choosing a model, and thus a case matrix structure, is a trade off between accuracy and computational cost. Due to decades of experience, there are default structures that are sufficient for current light water reactor analysis. However, this is not necessarily the case for advanced reactor designs. Therefore, it is advantageous to understand the sources of error in cross section interpolation models so that the quality of a case matrix may be improved. A mathematical framework for these models is presented in this work that provides a more rigorous connection between the nuclear engineering field's cross section interpolation methods and the broader mathematical field of function approximation. The two cross section models examined in this paper were found to utilize Lagrange interpolation and are a subset of Lagrange tensor products. Classical results of Lagrange polynomial error analysis were then applied to the partial derivative and multiple tables models to derive expressions for the total point-wise error. The analytical results classify the total error into two parts: the model form error and interpolation error. Finally, based on our observations, a better foundation for improving the quality of a case matrix is proposed. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Graph convolutional networks applied to unstructured flow field data

Abstract Many scientific and engineering processes produce spatially unstructured data. However, most data-driven models require a feature matrix that enforces both a set number and order of features for each sample. They thus cannot be easily constructed for an unstructured dataset. Therefore, a graph based data-driven model to perform inference on fields defined on an unstructured mesh, using a graph convolutional neural network (GCNN) is presented. The ability of the method to predict global properties from spatially irregular measurements with high accuracy is demonstrated by predicting the drag force associated with laminar flow around airfoils from scattered velocity measurements. The network can infer from field samples at different resolutions, and is invariant to the order in which the measurements within each sample are presented. The GCNN method, using inductive convolutional layers and adaptive pooling, is able to predict this quantity with a validation R 2 above 0.98, and a Normalized Mean Squared Error below 0.01, without relying on spatial structure.

Ogoke, Francis (ORCID:0000000224327783)↗

Nuclear Matrix Elements for Neutrinoless Double-Beta Decay

Neutrinoless double-beta decay ($0\nu\beta\beta$) is a rare hypothesised process that, if discovered, would establish that the neutrino is Majorana, that is, it is its own antiparticle. Interpretation of experimental results relies on knowledge of nuclear matrix elements, whose large model uncertainty is the limiting factor in comparing measured (bounds on) half-lives to the neutrino mass. Nuclear effective field theory and lattice QCD have the potential to compute these matrix elements with better control over uncertainties, enhancing the discovery potential of next-generation $0\nu\beta\beta$ experiments. This work will survey various lattice QCD double-beta decay calculations and discuss their implications.

Grebe, Anthony V. [Fermilab] (ORCID:00000003103201↗

Matrix Diffusion in Fractured Media: New Insights Into Power Law Scaling of Breakthrough Curves

Here, we develop a theoretical model for power law tailing behavior of transport in fractured rock based on the relative dominance of the decay rate of the advective travel time distribution, modeled using a Pareto distribution (with tail decaying as ~ time -(1+α) ), versus matrix diffusion, modeled using a Lévy distribution. The theory predicts that when the advective travel time distribution decays sufficiently slowly (α<1), the late-time decay rate of the breakthrough curve is -(1+α/2) rather than the classical -3/2. Yet, if α>1, the -3/2 decay rate is recovered. For weak matrix diffusion or short advective first breakthrough times, we identify an early-time regime where the breakthrough curve follows the Pareto distribution, before transitioning to the late-time decay rate. The theoretical predictions are validated against particle tracking simulations in the three-dimensional discrete fracture network simulator dfnWorks, where matrix diffusion is incorporated using a time domain random walk.

58 GEOSCIENCES↗

Ground-state-based model reduction with unitary circuits

Here, we present a method to numerically obtain low-energy effective models based on a unitary transformation of the ground state. The algorithm finds a unitary circuit that transforms the ground state of the original model to a projected wavefunction with only the low-energy degrees of freedom. The effective model can then be derived using the unitary transformation encoded in the circuit. We test our method on the one-dimensional and two-dimensional square-lattice Hubbard model at half-filling, and obtain more accurate effective spin models than the standard perturbative approach.

Hubbard model↗

Corrigendum to “High-fidelity micro-scale modeling of the thermo-visco-plastic behavior of carbon fiber polymer matrix composites” [Compos Struct 134 (2015) 132–141]

The authors would like to inform that one of the modifications proposed in the article “High-fidelity micro-scale modeling of the thermo-visco-plastic behavior of carbon fiber polymer matrix composites” [1] was found to be unnecessary: the paraboloid yield criterion is sufficient to describe the shear behavior of the epoxy matrix considered (Epoxy 3501-6). The authors recently noted that the experimental work [2] used to validate the pure matrix response considered engineering shear strain instead of its tensorial counter-part, which caused the apparent inconsistency with the paraboloid yield criterion. A recently proposed temperature dependency law for glassy polymers is evaluated herein, thus better agreement with the experimental results for this epoxy is observed.

36 MATERIALS SCIENCE↗

Modelling dislocation-graphene interactions in a BCC Fe matrix by molecular dynamics simulations and gradient plasticity theory

We report graphene nanosheets (GNS) can enhance the strength and ductility of metal-based composites as they can obstruct the propagation of dislocations. The present article employs Molecular Dynamics (MD) simulations to investigate dislocation-GNS interaction mechanisms and possible influencing factors, including the number of GNS layers, the thickness of the metallic amorphous layer and the C - C bond strength. The results indicated that the shear strength of the metal/GNS interface and the bending stiffness of GNS determined the ability of GNS to block dislocation transmission. A physically based phenomenological parameter that can capture such dislocation-GNS interactions is the mechanical interface energy that has been put forth within gradient plasticity. By fitting the theoretical expressions to the simulation data, it was possible to obtain estimates for the mechanical interface energy for the GNS. It was found that increasing the GNS layers and adding an amorphous layer resulted in a strengthening in the stress–strain response and increased the value of this interfacial parameter. This indicates that the mechanical interfacial energy can be a unified measure for capturing and tuning the strength of various interfaces such as grain boundaries, GNS, amorphous-crystalline interface and bimetal interfaces.

36 MATERIALS SCIENCE↗

Analysis of double stub tuner control stability in a phased array antenna with strong cross-coupling

Active stub tuning with a fast ferrite tuner (FFT) has greatly increased the ability of fusion ion cyclotron range of frequency (ICRF) systems (50–100 MHz) by allowing for the antenna system to respond dynamically to changes in the plasma load impedance such as during the L–H transition or edge localized modes (ELMs). A high power waveguide double-stub tuner is under development for use with the Alcator C-Mod lower hybrid current drive (LHCD) system at 4.6 GHz. The amplitude and relative phase shift between adjacent columns of an LHCD antenna are critical for control of the launched n || spectrum. Adding a double-stub tuning network will perturb the phase and amplitude of the forward wave particularly if the unmatched reflection coefficient is high. This effect can be compensated by adjusting the phase of the low power microwave drive for each klystron amplifier. Cross-coupling of the reflected power between columns of the launcher must also be considered. The problem is simulated by cascading a scattering matrix for the plasma provided by a linear coupling model with the measured launcher scattering matrix and that of the FFTs. The solution is advanced in an iterative manner similar to the time-dependent behavior of the real system. System performance is presented under a range of edge density conditions from under-dense to over-dense and a range of launched n || . Simulations suggest power reflection coefficients (Γ 2 ) of less than 1% with no contamination of the n || spectrum. Instability of the FFT tuning network can be difficult for certain plasma conditions and relative phasings, but reducing the control gain of the FFT network stabilizes the system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

TDCOSMO - XVI. Measurement of the Hubble constant from the lensed quasar WGD 2038–4008

Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant (H0). The TDCOSMO Collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar WGD 2038−4008 using new time-delay measurements and previous mass models by TDCOSMO. This is the first TDCOSMO lens to incorporate multiple lens modeling codes and the full time-delay covariance matrix into the cosmological inference. The models are fixed before the time delay is measured, and the analysis is performed blinded with respect to the cosmological parameters to prevent unconscious experimenter bias. We obtain DΔ t = 1.68−0.38+0.40 Gpc using two families of mass models, a power-law describing the total mass distribution, and a composite model of baryons and dark matter, although the composite model is disfavored due to kinematics constraints. In a flat ΛCDM cosmology, we constrain the Hubble constant to be H0 = 65−14+23 km s−1 Mpc−1. The dominant source of uncertainty comes from the time delays, due to the low variability of the quasar. Future long-term monitoring, especially in the era of the Vera C. Rubin Observatory’s Legacy Survey of Space and Time, could catch stronger quasar variability and further reduce the uncertainties. This system will be incorporated into an upcoming hierarchical analysis of the entire TDCOSMO sample, and improved time delays and spatially-resolved stellar kinematics could strengthen the constraints from this system in the future.Key words: gravitational lensing: strong / cosmological parameters / distance scale⋆ Corresponding author; kcwong19@gmail.com.⋆⋆ NHFP Einstein fellow.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of double stub tuner control stability in a many element phased array antenna with strong cross-coupling

Active stub tuning with a fast ferrite tuner (FFT) allows for the system to respond dynamically to changes in the plasma impedance such as during the L-H transition or edge localized modes (ELMs), and has greatly increased the effectiveness of fusion ion cyclotron range of frequency systems. A high power waveguide double-stub tuner is under development for use with the Alcator C-Mod lower hybrid current drive (LHCD) system. Exact impedance matching with a double-stub is possible for a single radiating element under most load conditions, with the reflection coefficient reduced from Γ to Γ{sup 2} in the “forbidden region.” The relative phase shift between adjacent columns of a LHCD antenna is critical for control of the launched n{sub ∥} spectrum. Adding a double-stub tuning network will perturb the phase of the forward wave particularly if the unmatched reflection coefficient is high. This effect can be compensated by adjusting the phase of the low power microwave drive for each klystron amplifier. Cross-coupling of the reflected power between columns of the launcher must also be considered. The problem is simulated by cascading a scattering matrix for the plasma provided by a linear coupling model with the measured launcher scattering matrix and that of the FFTs. The solution is advanced in an iterative manner similar to the time-dependent behavior of the real system. System performance is presented under a range of edge density conditions from under-dense to over-dense and a range of launched n{sub ∥}.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗