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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Multiscale modeling of solute diffusion in triblock copolymer membranes

We develop a multiscale simulation model for diffusion of solutes through porous triblock copolymer membranes. The approach combines two techniques: self-consistent field theory (SCFT) to predict the structure of the self-assembled, solvated membrane and on-lattice kinetic Monte Carlo (kMC) simulations to model diffusion of solutes. Solvation is simulated in SCFT by constraining the glassy membrane matrix while relaxing the brush-like membrane pore coating against the solvent. The kMC simulations capture the resulting solute spatial distribution and concentration-dependent local diffusivity in the polymer-coated pores; we parameterize the latter using particle-based simulations. We apply our approach to simulate solute diffusion through nonequilibrium morphologies of a model triblock copolymer, and we correlate diffusivity with structural descriptors of the morphologies. We also compare the model’s predictions to alternative approaches based on simple lattice random walks and find our multiscale model to be more robust and systematic to parameterize. Furthermore, our multiscale modeling approach is general and can be readily extended in the future to other chemistries, morphologies, and models for the local solute diffusivity and interactions with the membrane.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Semi-Annual Report for Horizontal Compact High Temperature Gas Reactor (HC-HTGR) Development during Performance Period April 2023 – September 2023

Horizontal Compact High Temperature Gas Reactor (HC-HTGR) is being designed by a multi-disciplinary team of nuclear, mechanical, and structural engineers under the support of a DOE-NE Advanced Reactor Demonstration Program’s Advanced Reactor Concepts-20 (ARC-20) award. The objective of this ARC-20 project is to deliver a conceptual design for the proposed HC-HTGR in 3 years and support its commercialization as a safe, low-cost HTGR. Argonne National Laboratory (Argonne) is responsible for the design and analysis of the reactor cavity cooling system (RCCS) as a safety system for passive decay heat removal of the reactor concept. Additionally, Argonne is providing analysis of the primary coolant system to ensure temperatures within the core remain below safety margins during steady-state and potential accident scenarios. This fourth semi-annual report summarized the progress made at Argonne on the two tasks during the second half of FY23. As a part of the RCCS design task, recent efforts have been made to complete a conceptual design of the RCCS for the HC-HTGR, including the design update of the water panel and system configuration favorable in point of view of fabrication and system operation. Design calculations were conducted under various heat load conditions to validate the system design. Transient simulations using RELAP5-3D were conducted to investigate system dynamics under transients of interest and to evaluate the system performance in the design condition. The results demonstrated the overall system feasibility that the RCCS design maintains structures temperatures lower than maximum allowable temperature with sufficient system inventory without any active heat sink. In the primary system thermal hydraulics task, the preliminary analysis was performed for a long term pressurized conduction cooldown (PCC) transient. This analysis used a combination of a fully resolved and coarse homogenized mesh to predict the temperature distribution for the steady-state initial condition and the PCC transient. The steady-state initial condition was determined using a fully resolved full core model with 3D solid to 1D fluid coupling. The fully resolved mesh was also used to model the first 20 seconds of the PCC. The temperature difference between fuel pins and the graphite matrix becomes minimal and the dominant heat transfer shifts to a larger scale radially towards the RCCS. After 20 seconds, a homogenized coarse mesh is used, greatly reducing the computational costs of the model. These results demonstrate that the core is designed to passively remove enough decay heat in a protected loss of primary coolant flow to prevent an unsafe rise of core temperatures.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Modeling and Harmonic Stability of MMC-HVDC With Passive Circulating Current Filters

A modular multilevel converter (MMC) is an emerging converter technology widely used in wind power applications using high-voltage direct current (HVDC) transmission, and it has been largely investigated in medium-voltage solar harvesting as well as electric motor drives. Different from studies dealing with active circulating current control loops, this paper focuses on impedance modeling and harmonic stability studies for MMCs with passive circulating current filters (PCCFs). A power stage circuit including a PCCF is transformed and remodeled to obtain small-signal formulations. This transformation leads to additional high-order matrix computing complexity due to the added impedance subnetwork matrix from the PCCF as well as relevant crossed-frequency impacts. To simplify the computational burden, the proposed model aims to solve one arm equation instead of all six MMC arm equations. To achieve this result, additional challenges occur when the MMC is connected to a renewable energy source instead of the grid voltage, where low-level zero-sequence components might exist in the neutral point of three phases, especially in the case of no integrated advanced modular voltage balancing control, which prevents computational simplification. To overcome this challenge, further impedance matrix adjustments are conducted in this paper to theoretically suppress the impact of zero-sequence harmonics when obtaining small-signal impedance, taking into account the frequency coupling effect. Finally, simulations are carried out to validate the developed impedance model under different operating scenarios. Harmonic stability case studies of MMCs with PCCFs connected to renewable energy current sources are also presented, where frequency-domain stability analyses based on Bode diagrams are compared to time-domain simulation waveforms and FFTs, which validate the effectiveness of the proposed impedance model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Analog B ( M 1 ) strengths in the T z = ± 3 2 mirror nuclei Mn 47 and Ti 47

The lifetimes of the first excited 7 2 − states in the T z = ± 3 2 mirror nuclei Mn 47 and Ti 47 have been extracted utilizing the γ -ray line shape method, giving τ = 687 ( 36 ) ps and τ = 331 ( 15 ) ps respectively. Since these transitions are essentially pure M 1 transitions, these results allow for a high-precision comparison of analog M 1 strengths in mirror nuclei. The two analog B ( M 1 ) s are observed to be identical to a precision of about 10 % . The expected dependence of the transition matrix element with T z has been used to extract the separate isoscalar and isovector components of the transition strength, and the results are discussed in the context of predictions, based on the isospin formalism, regarding analog B ( M 1 ) strengths. Published by the American Physical Society 2024

39 ≤ A ≤ 58↗
DOE OSTI2026

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Particle-Based Conditional Sampling Scheme for the Simulation of Transport in Fractured Rock With Diffusion Into Stagnant Water and Rock Matrix

In situ experiments and field-scale characterization studies have pointed out that, in fractured crystalline media, groundwater flow is highly channelized. This implies that, at the scale of a single fracture, only part of the fracture surface area is in contact with flowing water, while the rest of in-plane water is essentially stagnant and can be accessed by solutes via molecular diffusion. Despite their importance for contaminant retention, to date, there are no numerical or analytical approaches that could be used to assess the implication of stagnant water zones on solute transport in realistic large-scale Discrete Fracture Network-based models. In this study, we present an efficient and flexible algorithm for the simulation of transport in fractured rock with diffusion into stagnant water and rock matrix. The algorithm is a generalization of a previously developed numerical framework for time domain particle tracking in sparsely fractured rock. The key of the generalization is that total time in fracture (τ f ) is first evaluated using a Monte Carlo sampling and then a second sampling is performed conditioned on τ f . The algorithm has been successfully validated against existing independent solutions and the implication of diffusion into stagnant water and secondary diffusion into the matrix has been assessed for a realistic modeling scenario. The results show that, due to diffusion into stagnant water, contaminants are more strongly retarded. This increased retention is more significant for sorbing species, as a larger number of sorption sites is accessible. A high sensitivity to the flowing channel/stagnant water zone geometry has also been observed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The effects of continuum fitting on Lyman-α forest correlations

Correlations of fluctuations of the flux in Lyman-α forests of high-redshift quasars have been observed by the Baryonic Acoustic Oscillation Spectroscopy Survey (BOSS) and the Dark Energy Spectroscopy Instrument (DESI) survey where they have revealed the effects of baryon acoustic oscillations (BAO). In order to fit the correlation functions to a physical model and thereby constrain cosmological parameters, it is necessary to take into account the effects of fitting the observed spectra to a template about which the fluctuations are measured. In this paper we use mock spectra to test the distortion matrix technique that has been used since the final BOSS data release to appropriately distort the models. We show that while percent-level effects on the derived forest bias parameters may be present, the technique works sufficiently well that the determination of the BAO peak position is not affected at the percent level. We introduce modifications in the technique used by DESI that were not in the original applications and suggest further possibilities for improvements.

Lyman alpha forest↗
DOE OSTI2023

Quantum Phase Diagram and Spontaneously Emergent Topological Chiral Superconductivity in Doped Triangular-Lattice Mott Insulators

The topological superconducting state is a highly sought-after quantum state hosting topological order and Majorana excitations. Here, in this Letter, we explore the mechanism to realize the topological superconductivity (TSC) in the doped Mott insulators with time-reversal symmetry (TRS). Through large-scale density matrix renormalization group study of an extended triangular-lattice t–J model on the six- and eight-leg cylinders, we identify a d + id-wave chiral TSC with spontaneous TRS breaking, which is characterized by a Chern number C=2 and quasi-long-range superconducting order. We map out the quantum phase diagram with by tuning the next-nearest-neighbor (NNN) electron hopping and spin interaction. In the weaker NNN-coupling regime, we identify a pseudogaplike phase with a charge stripe order coexisting with fluctuating superconductivity, which can be tuned into d -wave superconductivity by increasing the doping level and system width. The TSC emerges in the intermediate-coupling regime, which has a transition to a d-wave superconducting phase with larger NNN couplings. The emergence of the TSC is driven by geometrical frustrations and hole dynamics which suppress spin correlation and charge order, leading to a topological quantum phase transition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Artificial extracellular matrix scaffolds of mobile molecules enhance maturation of human stem cell-derived neurons

Human induced pluripotent stem cell (hiPSC) technologies offer a unique resource for modeling neurological diseases. However, iPSC models are fraught with technical limitations including abnormal aggregation and inefficient maturation of differentiated neurons. These problems are in part due to the absence of synergistic cues of the native extracellular matrix (ECM). We report on the use of three artificial ECMs based on peptide amphiphile (PA) supramolecular nanofibers. All nanofibers display the laminin-derived IKVAV signal on their surface but differ in the nature of their non-bioactive domains. Here, we find that nanofibers with greater intensity of internal supramolecular motion have enhanced bioactivity toward hiPSC-derived motor and cortical neurons. Proteomic, biochemical, and functional assays reveal that highly mobile PA scaffolds caused enhanced β1-integrin pathway activation, reduced aggregation, increased arborization, and matured electrophysiological activity of neurons. Our work highlights the importance of designing biomimetic ECMs to study the development, function, and dysfunction of human neurons.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Non-perturbative renormalization scheme for the C P -odd three-gluon operator

We define a regularization-independent momentum-subtraction scheme for the C P -odd three-gluon operator at dimension six. This operator appears in effective field theories for heavy physics beyond the Standard Model, describing the indirect effect of new sources of C P -violation at low energies. In a hadronic context, it induces permanent electric dipole moments. The hadronic matrix elements of the three-gluon operator are non-perturbative objects that should ideally be evaluated with lattice QCD. We define a non-perturbative renormalization scheme that can be implemented on the lattice and we compute the scheme transformation to \( \overline{\mathrm{MS}} \) at one loop. Our calculation can be used as an interface to future lattice-QCD calculations of the matrix elements of the three-gluon operator, in order to obtain theoretically robust constraints on physics beyond the Standard Model from measurements of the neutron electric dipole moment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Electronic Structure Theory and Novel Materials

This grant supported research on electronic structure and materials theory, with focus on three main issues: (i) novel techniques to deal with correlation in the electronic ground-state, (ii) topological materials, (iii) the phase diagram of lattice spin models. Regarding (i), we applied to the homogeneous electron liquid an approach that we previously developed in the context of molecular systems. In this scheme the electronic occupation probabilities and the natural spin orbitals are used to construct an approximate two-body density matrix for the electronic ground-state. Regarding (ii) we used standard electronic structure methods based on density functional theory to model topological materials and interpret experimental observations. Finally, regarding (iii) we further developed a numerical approach to compute the renormalized couplings within real space renormalization group theory in the context of lattice spin models. The main findings were the following. (i) We found that with our approximate two-body density matrix, which works well for small molecules, is not sufficiently accurate for condensed phase systems. Missing a systematic way of improving on the adopted approximations, we decided not to pursue this approach. (ii) We performed two studies. In one, we investigated the influence of Te defects on the topological properties of a WTe2 monolayer, finding that while Te vacancies, even in modest concentration, destroy the topological character, Te adatoms do not, consistent with a recent experiment. In another study, we predicted Weyl semimetal character and strong anomalous Hall effect in the Heusler compensated ferrimagnet Ti2MnAl. (iii) We developed a new Monte Carlo method to do real space renormalization group calculations for lattice spin models. We subsequently extended the scheme to deal with lattice spin models in presence of quenched disorder, finding that the approach can distinguish systems with finite and strong disorder. In the finite disorder case, the method allows one to find with good approximation the critical coupling distribution and the critical exponents.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A graph signal processing‐based multiple model Kalman filter ( GSP‐MMKF ) tool for predictive analytics: An air separation unit process application

Abstract The industrial Air Separations Unit (ASU) is a complicated and tightly operated process. The use of dynamic process analytics is also a key element of safe and economic operation of these processes, with increasing focus on predictive analytics to take preemptive actions. With the availability of real‐time data from hundreds of sensors, the data analysis process should also consider the topology of the data, as seen in sensor networks. In this paper, a novel tool is presented that considers the complex connectivity patterns in the sensor network and uses local adaptive disturbance estimations to predict global network‐scale trends. The paper introduces the emerging field of Graph Signal Processing (GSP) and presents a rigorous derivation of the tool starting from the extraction of the sensor‐network (in a graph theoretical sense) from the data. This network, which is in the form of a matrix, is then used to derive a Kalman‐filter type of state‐space model driven by input disturbances. Multiple disturbance models (e.g., step, ramp, periodic) are included to allow the model to have different kinds of disturbance propagation. Each graph node (representing the sensors used) dynamically adapts to the most recent detected disturbance individually. These estimated disturbances are propagated to the global network using the graph. Modifications to ensure stability are also discussed. The fidelity of the tool is tested on certain downtime events and the paper concludes by discussing the advantages of the method and planned future improvements.

Ghosh, Sambit↗
DOE OSTI2022

Fast increased fidelity samplers for approximate Bayesian Gaussian process regression

Gaussian processes (GPs) are common components in Bayesian non-parametric models having a rich methodological literature and strong theoretical grounding. The use of exact GPs in Bayesian models is limited to problems containing several thousand observations due to their prohibitive computational demands. We develop a posterior sampling algorithm using H-matrix approximations that scales at O(n log 2 n). We show that this approximation’s Kullback-Leibler divergence to the true posterior can be made arbitrarily small. Though multidimensional GPs could be used with our algorithm, d-dimensional surfaces are modeled as tensor products of univariate GPs to minimize the cost of matrix construction and maximize computational efficiency. We illustrate the performance of this fast increased fidelity approximate GP, FIFA-GP, using both simulated and non-synthetic data sets

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Quantifying and Controlling Entanglement in the Quantum Magnet Cs 2 CoCl 4

The lack of methods to experimentally detect and quantify entanglement in quantum matter impedes our ability to identify materials hosting highly entangled phases, such as quantum spin liquids. We thus investigate the feasibility of using inelastic neutron scattering (INS) to implement a model-independent measurement protocol for entanglement based on three entanglement witnesses: one-tangle, two-tangle, and quantum Fisher information (QFI). We perform high-resolution INS measurements on Cs 2 CoCl 4 , a close realization of the S = 1/2 transverse-field XXZ spin chain, where we can control entanglement using the magnetic field, and compare with density-matrix renormalization group calculations for validation. The three witnesses allow us to infer entanglement properties and make deductions about the quantum state in the material. Overall, we find QFI to be a particularly robust experimental probe of entanglement, whereas the one and two-tangles require more careful analysis. Our results lay the foundation for a general entanglement detection protocol for quantum spin systems.

1-dimensional spin chains↗
DOE OSTI2025

Towards excitations and dynamical quantities in correlated lattices with density matrix embedding theory

Density matrix embedding theory (DMET) provides a framework to describe ground-state expectation values in strongly correlated systems, but its extension to dynamical quantities is still an open problem. We show one route to obtaining excitations and dynamical spectral functions by using the techniques of DMET to approximate the matrix elements that arise in a single-mode inspired excitation ansatz. We demonstrate this approach in the one-dimensional Hubbard model, comparing the neutral excitations, single-particle density of states, charge, and spin dynamical structure factors to benchmarks from the Bethe ansatz and density matrix renormalization group. Finally, our work highlights the potential of these ideas in building computationally efficient approaches for dynamical quantities.

1-dimensional systems↗
DOE OSTI2021

Empirical Estimation of Shortest Route Length along U.S. Interstate Highways Based on Great Circle Distance

In this study, 98 regression models were specified for easily estimating shortest distances based on great circle distances along the U.S. interstate highways nationwide and for each of the continental 48 states. This allows transportation professionals to quickly generate distance, or even distance matrix, without expending significant efforts on complicated shortest path calculations. For simple usage by all professionals, all models are present in the simple linear regression form. Only one explanatory variable, the great circle distance, is considered to calculate the route distance. For each geographic scope (i.e., the national or one of the states), two different models were considered, with and without the intercept. Based on the adjusted R-squared, it was observed that models without intercepts generally have better fitness. Additionally, all these models generally have good fitness with the linear regression relationship between the great circle distance and route distance. At the state level, significant variations in the slope coefficients between the state-level models were also observed. Furthermore, a preliminary analysis of the effect of highway density on this variation was conducted.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Measuring σ 8 using DESI Legacy Imaging Surveys Emission-Line galaxies and Planck CMB lensing, and the impact of dust on parameter inference

Measuring the growth of structure is a powerful probe for studying the dark sector, especially in light of the σ 8 tension between primary CMB anisotropy and low-redshift surveys. This paper provides a new measurement of the amplitude of the matter power spectrum, σ 8 , using galaxy-galaxy and galaxy-CMB lensing power spectra of Dark Energy Spectroscopic Instrument Legacy Imaging Surveys Emission-Line Galaxies and the Planck 2018 CMB lensing map. We create an ELG catalog composed of 24 million galaxies and with a purity of 85%, covering a redshift range 0 < z < 3, with z mean = 1.09. We implement several novel systematic corrections, such as jointly modeling the contribution of imaging systematics and photometric redshift uncertainties to the covariance matrix. We also study the impacts of various dust maps on cosmological parameter inference. We measure the cross-power spectra over f sky = 0.25 with a signal-to-background ratio of up to 30σ. We find that the choice of dust maps to account for imaging systematics in estimating the ELG overdensity field has a significant impact on the final estimated values of σ 8 and Ω M , with far-infrared emission-based dust maps preferring σ 8 to be as low as 0.702 ± 0.030, and stellar-reddening-based dust maps preferring as high as 0.719 ± 0.030. The highest preferred value is at ∼ 3 σ tension with the Planck primary anisotropy results. These findings indicate a need for tomographic analyses at high redshifts and joint modeling of systematics.

79 ASTRONOMY AND ASTROPHYSICS↗