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,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Assessing sensitivities of climate model weighting to multiple methods, variables, and domains in the south-central United States

Given the increasing use of climate projections and multi-model ensemble weighting for a diverse array of applications, this project assesses the sensitivities of climate model weighting strategies and their resulting ensemble means to multiple components, such as the weighting schemes, climate variables, or spatial domains of interest. The purpose of this study is to assess the sensitivities associated with multi-model weighting strategies. The analysis makes use of global climate models from the Coupled Model Intercomparison Project Phase 5 (CMIP5) and their statistically downscaled counterparts created with the localized constructed analogs (LOCA) method. This work focuses on historical and projected future mean precipitation and daily high temperatures of the south-central United States. Results suggest that the model weights and the corresponding weighted model means can be sensitive to the weighting strategy that is applied. For instance, when estimating model weights based on Louisiana precipitation, the weighted projections show a wetter and cooler south-central domain in the future compared to other weighting strategies. Alternatively, for example, when estimating model weights based on New Mexico temperature, the weighted projections show a drier and warmer south-central domain in the future. However, when considering the entire south-central domain in estimating the model weights, the weighted future projections show a compromise in the precipitation and temperature estimates. As for uncertainty, our matrix of results provided a more certain picture of future climate compared to the spread in the original model ensemble. If future impact assessments utilize weighting strategies, then our findings suggest that how the specific weighting strategy is used with climate projections may depend on the needs of an impact assessment or adaptation plan.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

VANGL2 interacts with integrin αv to regulate matrix metalloproteinase activity and cell adhesion to the extracellular matrix

Highlights: • VANGL2 regulates cell adhesion to ECM proteins. • VANGL2 function influences cell surface integrin αvβ3 expression. • Inhibition of MMP activity rescues the adhesion defect in VANGL2 knockdown cells. • Integrin αv is required for VANGL2-dependent MMP2 activation. • VANGL2 binds integrin αvβ3. Planar cell polarity (PCP) proteins are implicated in a variety of morphogenetic processes including embryonic cell migration and potentially cancer progression. During zebrafish gastrulation, the transmembrane protein Vang-like 2 (VANGL2) is required for PCP and directed cell migration. These cell behaviors occur in the context of a fibrillar extracellular matrix (ECM). While it is thought that interactions with the ECM regulate cell migration, it is unclear how PCP proteins such as VANGL2 influence these events. Using an in vitro cell culture model system, we previously showed that human VANGL2 negatively regulates membrane type-1 matrix metalloproteinase (MMP14) and activation of secreted matrix metalloproteinase 2 (MMP2). Here, we investigated the functional relationship between VANGL2, integrin αvβ3, and MMP2 activation. We provide evidence that VANGL2 regulates cell surface integrin αvβ3 expression and adhesion to fibronectin, laminin, and vitronectin. Inhibition of MMP14/MMP2 activity suppressed the cell adhesion defect in VANGL2 knockdown cells. Furthermore, our data show that MMP14 and integrin αv are required for increased proteolysis by VANGL2 knockdown cells. Lastly, we have identified integrin αvβ3 as a novel VANGL2 binding partner. Together, these findings begin to dissect the molecular underpinnings of how VANGL2 regulates MMP activity and cell adhesion to the ECM.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

A pathway to unveiling neutrinoless ββ decay nuclear matrix elements via γγ decay

We investigate the experimental feasibility of detecting second-order double-magnetic dipole (γγ-M1M1) decays from double isobaric analog states (DIAS), which have recently been found to be strongly correlated with the nuclear matrix elements of neutrinoless ββ decay. Using the nuclear shell model, we compute theoretical branching ratios for γγ-M1M1 decays and compare them with other competing processes, such as single-γ decay and proton emission, which represent the dominant decay channels. We also estimate the potential competition from internal conversion and internal pair creation, which can influence the decay dynamics. Additionally, we propose an experimental strategy based on using LaBr3 scintillators to identify γγ-M1M1 transitions from the DIAS amidst the background of the competing processes. Our approach emphasizes the challenges of isolating the rare γγ-M1M1 decay and suggests ways to enhance the experimental detection sensitivity. Our simulations suggest that it may be possible to access experimentally γγ-M1M1 decays from DIAS, shedding light on the neutrinoless ββ decay nuclear matrix elements.

Romeo, Beatriz↗
DOE OSTI2024

Predicting the tensile and compressive modulus of electrospun fiber mat‐reinforced hydrogels using the Halpin–Tsai equations

Abstract The reinforcement of mechanically‐weak hydrogels to yield composites with increased stiffness, strength, or toughness is a well‐established approach. In particular, introducing electrospun nanofibers into hydrogels is a common strategy for biomedical applications, as the resulting hierarchical structure mimics biology and allows for control over fiber diameter and alignment and tuning of mechanical properties. However, further study of the link between the constituent materials and the mechanical properties of the composite is uncommon. One potential model to understand the mechanical properties of fiber‐reinforced hydrogels involves the Halpin–Tsai equations, which relate the modulus values of the fibers and hydrogel matrix and the fiber volume fraction, to the modulus of the composite. To assess the application of this model to fiber‐reinforced hydrogels, predicted values were compared with experimental values from mechanical testing of a poly(ethylene glycol) (PEG) matrix reinforced with an electrospun polycaprolactone (PCL) fiber mat. Although the equations described these systems well in tension, providing a facile approach to identify a fiber volume fraction that will achieve a desired modulus, the Halpin–Tsai approach was less successful under compression. This study motivates additional investigation of the role of structural features of hydrogel composites in determining mechanical properties to enable design of materials for specific applications.

Polymer Science↗
DOE OSTI2021

Model-based and Model-free Designs for an Extended Continuous-time LQR with Exogenous Inputs

We present an extended linear quadratic regulator (LQR) design for continuous-time linear time-invariant (LTI) systems in the presence of exogenous inputs. We first propose a model-based solution with cost minimization guarantees for states and inputs using dynamic programming (DP). The control law consists of a combination of the optimal state feedback and an additional optimal term dependent on the exogenous inputs. The control gains for the two components are obtained by solving a set of matrix differential equations. We provide these solutions for both finite horizons and steady-state cases. In the second part of the paper, we formulate a reinforcement learning (RL) based algorithm which does not need any model information except the input matrix, and can compute an approximate steady-state LQR gain using measurements of the states, the control inputs, and the exogenous inputs. Both model-based and data-driven optimal control algorithms are tested with a numerical example under different exogenous inputs showcasing the effectiveness of the designs.

Mukherjee, Sayak↗
DOE OSTI2013

A Multiscale Modeling Approach to Analyze Filament-Wound Composite Pressure Vessels

A multiscale modeling approach to analyze filament-wound composite pressure vessels is developed in this article. The approach, which extends the Nguyen et al. model [J. Comp. Mater. 43 (2009) 217] developed for discontinuous fiber composites to continuous fiber ones, spans three modeling scales. The microscale considers the unidirectional elastic fibers embedded in an elastic-plastic matrix obeying the Ramberg-Osgood relation and J2 deformation theory of plasticity. The mesoscale behavior representing the composite lamina is obtained through an incremental Mori-Tanaka type model and the Eshelby equivalent inclusion method [Proc. Roy. Soc. Lond. A241 (1957) 376]. The implementation of the micro-meso constitutive relations in the ABAQUS® finite element package (via user subroutines) allows the analysis of a filament-wound composite pressure vessel (macroscale) to be performed. Failure of the composite lamina is predicted by a criterion that accounts for the strengths of the fibers and of the matrix as well as of their interface. The developed approach is demonstrated in the analysis of a filament-wound pressure vessel to study the effect of the lamina thickness on the burst pressure. The predictions are favorably compared to the numerical and experimental results by Lifshitz and Dayan [Comp. Struct. 32 (1995) 313].

Filament-wound composite, pressure vessel, multisc↗
DOE OSTI2016

VR-SCOSMO: A smooth conductor-like screening model with charge-dependent radii for modeling chemical reactions

To better represent the solvation effects observed along reaction pathways, and of ionic species in general, a charge-dependent variable-radii smooth conductor-like screening model (VR-SCOSMO) is developed. This model is implemented and parameterized with a third order density-functional tight binding quantum model, DFTB3/3OB-OPhyd, a quantum method which was developed for organic and biological compounds, utilizing a specific parameterization for phosphate hydrolysis reactions. Unlike most other applications with the DFTB3/3OB model, an auxiliary set of atomic multipoles is constructed from the underlying DFTB3 density matrix which is used to interact the solute with the solvent response surface. The resulting method is variational, produces smooth energies, and has analytic gradients. As a baseline, a conventional SCOSMO model with fixed radii is also parameterized. The SCOSMO and VR-SCOSMO models shown have comparable accuracy in reproducing neutral-molecule absolute solvation free energies; however, the VR-SCOSMO model is shown to reduce the mean unsigned errors (MUEs) of ionic compounds by half (about 2-3 kcal/mol). The VR-SCOSMO model presents similar accuracy as a charge-dependent Poisson-Boltzmann model introduced by Hou et al. [J. Chem. Theory Comput. 6, 2303 (2010)]. VR-SCOSMO is then used to examine the hydrolysis of trimethylphosphate and seven other phosphoryl transesterification reactions with different leaving groups. Two-dimensional energy landscapes are constructed for these reactions and calculated barriers are compared to those obtained from ab initio polarizable continuum calculations and experiment. Results of the VR-SCOSMO model are in good agreement in both cases, capturing the rate-limiting reaction barrier and the nature of the transition state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Two-scale model for the effect of physical aging in elastomers filled with hard nanoparticles

A two-scale model is developed, and solved numerically, to describe the mechanical behavior of elastomers filled with hard nanoparticles. Of particular interest is the slow recovery of the elastic modulus after large-amplitude oscillatory deformation. To account for this effect, the physical aging of the glassy bridges between the filler particles is captured with two thermal degrees of freedom for the matrix material, namely a kinetic and a configurational one. Formulating the two-scale model enriched with aging in a nonequilibrium thermodynamics context, first results in a constitutive relation for the Cauchy stress tensor. Second, the dynamics of physical aging is described, which eventually results in the slow recovery of the elastic modulus with waiting time. The proposed model is investigated numerically under large amplitude oscillatory shear deformation. Of particular interest in this respect is the coupling of the micro-scale dynamics with the physical aging on the macroscopic scale. This coupling is examined in detail, both in an approximate way using a Gaussian approximation, as well as numerically, under specific conditions. It turns out that the CONNFFESSIT approach (Laso and Öttinger 1993 ) can not be employed for the numerical solution of the model under arbitrary loading conditions because of the novel structure of the two-level coupling term. While a procedure for solving the model numerically for the case of strong applied deformation is presented in this paper, other solution methodologies need to be sought for the cases of weak and no applied deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Enhancing scalability of a matrix-free eigensolver for studying many-body localization

We propose several techniques to enhance the parallel scalability of a matrix-free eigensolver designed for studying many-body localization (MBL) of quantum spin chain models with nearest-neighbor interactions and on-site disorder. This type of problem is computationally challenging because the dimension of the associated Hamiltonian matrix grows exponentially with respect to the number of spins L, and we need to average over different realizations of the random disorder to obtain relevant statistical behavior. For each disorder realization, we need to compute eigenvalues from different regions of the spectrum and their corresponding eigenvectors. In previous work, the interior eigenstates for a single eigenvalue problem are computed via the shift-and-invert Lanczos algorithm. Due to the extremely high memory footprint of the LU factorizations, this technique is not well suited for large L’s. For example, we need thousands of compute nodes on modern high performance computing infrastructures to go beyond L = 24. The matrix-free approach does not suffer from this memory bottleneck, however, its scalability is limited by a computation and communication load imbalance. To reduce this imbalance and to significantly enhance the scalability of the matrix-free eigensolver, we reorder the matrix and leverage the consistent space runtime, CSPACER. We also show its efficiency in managing irregular communication patterns at scale compared to optimized MPI non-blocking two-sided and one-sided RMA implementation variants. This effort enables us to study MBL for spin chains with a larger number of spins. The efficiency and effectiveness of the proposed algorithm is demonstrated by computing eigenstates on a massively parallel many-core high performance computer.

METIS↗
DOE OSTI2025

CALPHAD-based Bayesian optimization to accelerate alloy discovery for high-temperature applications

Two crucial properties influencing the performance of high-temperature alloys are coefficient of thermal expansion (CTE) and phase constitution. It is desirable to have alloys with low CTE, which reduces CTE mismatch with the surface oxide and the likelihood of oxide spallation. Reducing the amount of brittle intermetallic phases such as Sigma (σ) enhances alloy ductility and processability. Here, we propose a multi-objective Bayesian Optimization (BO) model to simultaneously minimize CTE (at an operational temperature of 1150 °C) and T σ (temperature when the Sigma phase completely dissolves in the metal matrix), properties which are obtained from high-throughput CALculation of PHAse Diagrams (CALPHAD). The model successfully identifies several alloys with CTE ≤ 2 × 10 –5 /K and T σ ≤ 500 °C by exploring just 7% of the nickel–chromium–cobalt–aluminum–iron (Ni–Cr–Co–Al–Fe) composition space. Such multi-objective alloy design frameworks can be used to inform additive manufacturing experiments and accelerate alloy discovery for high-temperature energy applications.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Utilizing Metalized Fabrics for Liquid and Rip Detection and Localization

This paper proposes a novel technique for utilizing conductive textiles as a distributed sensor for detecting and localizing liquids (e.g., blood), rips (e.g., bullet holes), and potentially biosignals. The proposed technique is verified through both simulation and experimental measurements. Circuit theory is utilized to depict conductive fabric as a bounded, near-infinite grid of resistors. Solutions to the well-known infinite resistance grid problem are used to confirm the accuracy and validity of this modeling approach. Simulations allow for discontinuities to be placed within the resistor matrix to illustrate the effects of bullet holes within the fabric. A real-time experimental system was developed that uses a multiplexed Wheatstone bridge approach to reconstruct the resistor grid across the conductive fabric and detect liquids and rips. The resistor grid model is validated through a comparison of simulated and experimental results. Results suggest accuracy proportional to the electrode spacing in determining the presence and location of discontinuities in conductive fabric samples. Future work is focused on refining the experimental system to provide more accuracy in detecting and localizing events as well as developing a complete prototype that can be deployed for field testing. Potential applications include intelligent clothing, flexible, lightweight sensing systems, and combat wound detection.

Holland, Stephen↗
DOE OSTI2022

TEM-based phase characterization of $\mathrm{U–19Pu–10Zr}$ irradiated in ATR

We report the results of the microstructural examination of the irradiated U–19Pu–10Zr fuel. New insights gained using selective area electron diffraction analysis in a transmission electron microscope provide crystallographic information on the phases previously identified based primarily on the chemical compositions. After irradiation, the microstructure of the fuel separates into γ-(U, Pu, Zr) in the central region with a secondary nanocrystalline UO 2 phase, γ-(U, Pu, Zr) with a secondary ZrO 2 phase in the intermediate, and ζ-(U, Pu) in the outer regions. While the surface microstructural features and porosity indicate separation into three concentric zones consistent with the existing constituent redistribution model, crystallographically the fuel consists of two zones: inner (with γ-(U, Pu, Zr) as a matrix phase) and outer (with ζ-(U, Pu) as a matrix phase). Our results highlight the importance of crystal structure in phase identification and can be used to supplement the existing constituent redistribution model.

extreme environment↗
DOE OSTI2024

Suppressed electric quadrupole collectivity in Si 32

Lying between O 16 and Ca 40 , the s d shell is well described by robust phenomenological and nuclear theories. In this work, however, we highlight an unexplained reduction in electric-quadrupole strength in the rare isotope Si 32 , studied through sub-barrier Coulomb excitation. It is found that the oblate nature of the deformation is well reproduced, while the absolute scale of quadrupole deformation, however, is inhibited by approximately a factor of 2 compared to theoretical predictions. Through comparison with shell-model and calculations, we present a number of possible explanations for this inhibited E 2 strength. By comparing the results of these calculations to multiple observables, we conclude that there is a reduced role for out-of-space excitations in Si 32 , resulting in a reduction in the corrections normally applied to both models. Published by the American Physical Society 2024

20 ≤ A ≤ 38↗
DOE OSTI2014

Modeling the electrical resistivity of deformation processed metal-metal composites

Deformation processed metal–metal (matrix–reinforcement) composites (DMMCs) are high-strength, high-conductivity in situ composites produced by severe plastic deformation. The electrical resistivity of DMMCs is rarely investigated mechanistically and tends to be slightly higher than the rule-of-mixtures prediction. In this paper, we analyze several possible physical mechanisms (i.e. phonons, interfaces, mutual solution, grain boundaries, dislocations) responsible for the electrical resistivity of DMMC systems and how these mechanisms could be affected by processing conditions (i.e. temperature, deformation processing). As an innovation, we identified and assembled the major scattering mechanisms for specific DMMC systems and modeled their electrical resistivity in combination. From this analysis, it appears that filament coarsening rather than dislocation annihilation is primarily responsible for the resistivity drop observed in these materials after annealing and that grain boundary scattering contributes to the resistivity at least at the same magnitude as does interface scattering.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Topological phase diagram of generalized Su-Schrieffer-Heeger models with interactions

We investigate interacting Su-Schrieffer-Heeger (SSH) chains with two- and three-site unit cells using density matrix renormalization group (DMRG) simulations. By selecting appropriate filling fractions and sweeping across interaction strength 𝐽 𝑧 and dimerization 𝛿, we map out their phase diagrams and identify transition lines via entanglement entropy and magnetization measurements. In the two-site model, we observe the emergence of an interaction-induced antiferromagnetic intermediate phase between the topologically trivial and nontrivial regimes, as well as a critical region at negative 𝐽 𝑧 with suppressed magnetization and finite-size scaling of entanglement entropy. In contrast, the three-site model lacks an intermediate phase and exhibits asymmetric edge localization and antiferromagnetic ordering in both positive and negative 𝐽 𝑧 regimes. We further examine the response of edge states to Ising perturbations. In the two-site model, zero-energy edge modes are topologically protected and remain robust up to a finite interaction strength. However, in the three-site model, where the edge states reside at finite energy, this protection breaks down. Despite this, the edge-localized nature of these states survives in the form of polarized modes whose spatial profiles reflect the noninteracting limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

The relaxation matrix for symmetric tops with inversion symmetry. II. Line mixing effects in the ν 1 band of NH 3

Line mixing effects have been calculated in the ν 1 parallel band of self-broadened NH 3 . The theoretical approach is an extension of a semi-classical model to symmetric-top molecules with inversion symmetry developed in the companion paper [Q. Ma and C. Boulet, J. Chem. Phys. 144, 224303 (2016)]. This model takes into account line coupling effects and hence enables the calculation of the entire relaxation matrix. A detailed analysis of the various coupling mechanisms is carried out for Q and R inversion doublets. The model has been applied to the calculation of the shape of the Q branch and of some R manifolds for which an obvious signature of line mixing effects has been experimentally demonstrated. Comparisons with measurements show that the present formalism leads to an accurate prediction of the available experimental line shapes. Discrepancies between the experimental and theoretical sets of first order mixing parameters are discussed as well as some extensions of both theory and experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Model representations of kerogen structures: An insight from density functional theory calculations and spectroscopic measurements

Molecular structures of kerogen control hydrocarbon production in unconventional reservoirs. Significant progress has been made in developing model representations of various kerogen structures. These models have been widely used for the prediction of gas adsorption and migration in shale matrix. However, using density functional perturbation theory (DFPT) calculations and vibrational spectroscopic measurements, we here show that a large gap may still remain between the existing model representations and actual kerogen structures, therefore calling for new model development. Using DFPT, we calculated Fourier transform infrared (FTIR) spectra for six most widely used kerogen structure models. The computed spectra were then systematically compared to the FTIR absorption spectra collected for kerogen samples isolated from Mancos, Woodford and Marcellus formations representing a wide range of kerogen origin and maturation conditions. Limited agreement between the model predictions and the measurements highlights that the existing kerogen models may still miss some key features in structural representation. A combination of DFPT calculations with spectroscopic measurements may provide a useful diagnostic tool for assessing the adequacy of a proposed structural model as well as for future model development. This approach may eventually help develop comprehensive infrared (IR)-fingerprints for tracing kerogen evolution.

03 NATURAL GAS↗
DOE OSTI2025

Global phase diagram of D-wave superconductivity in the square-lattice t-J model

The Hubbard and closely related t-J models are exciting platforms for unconventional superconductivity (SC). Through state-of-the-art density matrix renormalization group calculations using the grand canonical ensemble, we address open issues regarding the ground-state phase diagram of the extended t-J model on a square lattice in the parameter regime relevant to cuprate superconductors. On large 8-leg cylinders, we demonstrate that the pure t-J model with only nearest-neighbor hoppings and superexchange interactions, for a wide range of doping ( δ = 0.1 to 0.2 ), hosts robust d-wave superconductivity possibly coexisting with weak unidirectional pair density wave. Furthermore, a small next nearest neighbor hopping t 2 suppresses pair and charge density waves, resulting in a uniform d-wave SC phase in both electron- and hole-doped cuprate model systems. Our work validates the t-J model as a proper minimum model for the emergence of superconductivity in cuprate superconductors.

Chen, Feng (ORCID:000000018196030X)↗