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 145 records · Page 8

Surrogate models for linear response

Linear response theory is a well-established method in physics and chemistry for exploring excitations of many-body systems. In particular, the quasiparticle random-phase approximation (QRPA) provides a powerful microscopic framework by building excitations on top of the mean-field vacuum; however, its high computational cost limits model calibration and uncertainty quantification studies. Here, we present two complementary QRPA surrogate models and apply them to study response functions of finite nuclei. One is a reduced-order model that exploits the underlying QRPA structure, while the other utilizes the recently developed parametric matrix model algorithm to construct a map between the system’s Hamiltonian and observables. Our benchmark applications, the calculation of the electric dipole polarizability of 180 Yb and the 𝛽-decay half-life of 80 Ni, show that both emulators can achieve 0.1%–1% accuracy while offering a 6–7 orders of magnitude speedup compared to state-of-the-art QRPA solvers. These results demonstrate that the developed QRPA emulators are well positioned to enable Bayesian calibration and large-scale studies of computationally expensive physics models describing the properties of many-body systems.

Beta decay↗

FEOTS v0.0.0: a new offline code for the fast equilibration of tracers in the ocean

Abstract. In this paper we introduce a new software framework for the offline calculation of tracer transport in the ocean. The Fast Equilibration of Ocean Tracers Software (FEOTS) is an end-to-end set of tools to efficiently calculate tracer distributions on a global or regional sub-domain using transport operators diagnosed from a comprehensive ocean model. To the best of our knowledge, this is the first application of a transport matrix model to an eddying ocean state. While a Newton–Krylov-based equilibration capability is still under development and not presented here, we demonstrate in this paper the transient modeling capabilities of FEOTS in an application focused on the Argentine Basin, where intense eddy activity and the Zapiola Anticyclone lead to strong mixing of water masses. The demonstration illustrates progress in developing offline passive tracer simulation capabilities, while highlighting the challenges of the impulse response functions approach in capturing tracer transports by a non-linear advection scheme. Our future work will focus on improving the computational efficiency of the code to reduce time-to-solution, using different basis functions to better represent non-linear advection operators, applying FEOTS to a parent model with unstructured grids (Ocean Model for Prediction Across Scales, MPAS-Ocean), and fully implementing a Newton–Krylov steady-state solver.

54 ENVIRONMENTAL SCIENCES↗

Analysis of a new multispecies tumor growth model coupling 3D phase-fields with a 1D vascular network

In this study, we present and analyze a mathematical model for tumor growth incorporating ECM erosion, interstitial flow, and the effect of vascular flow and nutrient transport. The model is of phase-field or diffused-interface type in which multiple phases of cell species and other constituents are separated by smooth evolving interfaces. The model involves a mesoscale version of Darcy’s law to capture the flow mechanism in the tissue matrix. Modeling flow and transport processes in the vasculature supplying the healthy and cancerous tissue, one-dimensional (1D) equations are considered. Since the models governing the transport and flow processes are defined together with cell species models on a three-dimensional (3D) domain, we obtain a 3D–1D coupled model.

3D–1D coupled blood flow models↗

Dual Phase Synergy Enabled Large Elastic Strains of Nanoinclusions in a Dislocation Slip Matrix Composite

Freestanding nanomaterials (such as nanowires, nanoribbons and nanotubes) are known to exhibit ultra-large elastic strains and ultra-high strengths. However, harnessing their superior intrinsic mechanical properties in bulk composites has proven to be difficult. A recent breakthrough has overcome this difficulty by using a martensitic phase transforming matrix, in which ultra-large elastic strains approaching the theoretical limit is achieved in Nb nanowires embedded in the matrix. This discovery, whereas breaking a long-standing challenge, still limits our ability of harnessing the exceptional properties of nanomaterials and developing ultrahigh strength bulk materials to a narrow selection of phase transforming alloy matrices. In this study, we investigated the possibility to harness the intrinsic mechanical properties of nanoinclusions in conventional dislocation slip matrix based on a principle of synergy between the inclusion and the matrix. The small spacing between the densely populated hard and dislocation-impenetrable nanoinclusions departmentalize the plastic matrix into small domains to effectively impede dislocation motion within the matrix, inducing significant strengthening and large local elastic strains of the matrix, which in turn induced large elastic strains in the nanoinclusions. This dual phase synergy is verified in a Ti 3 Sn inclusions/B2-NiTi(Fe) plastic matrix model materials system. The maximum elastic strain of Ti 3 Sn inclusion obtained in the dislocation slip matrix is comparable to that achieved in a phase transforming matrix. As a result, this finding opens new opportunities for the development of high-strength nanocomposites.

42 ENGINEERING↗

Modeling creep behavior in ceramic matrix composites

Here, in this work, a three-dimensional viscoplasticity formulation with progressive damage is developed and used to investigate the complex time-dependent constituent load transfer and progressive damage behavior in ceramic matrix composites (CMCs) subjected to creep. The viscoplasticity formulation is based on Hill's orthotropic plastic potential, an associative flow rule, and the Norton-Bailey creep power law with Arrhenius temperature dependence. A fracture mechanics-informed isotropic matrix damage model is used to account for CMC brittle matrix damage initiation and propagation, in which two scalar damage variables capture the effects of matrix porosity as well as matrix property degradation due to matrix crack initiation and propagation. The Curtin progressive fiber damage model is utilized to simulate progressive fiber failure. The creep-damage formulation is subsequently implemented as a constitutive model in the generalized method of cells (GMC) micromechanics formulation to simulate time-dependent deformation and material damage under creep loading conditions. The developed framework is used to simulate creep of single fiber SiC/SiC microcomposites. Simulation results are in excellent agreement with experimental and numerical data available in the literature.

36 MATERIALS SCIENCE↗

The Yang-Mills duals of small AdS black holes

Abstract We study the largeNmatrix model for the index of 4d$$ \mathcal{N} $$ N = 4 Yang-Mills theory and its truncations to understand the dual AdS 5 black holes. Numerical studies of the truncated models provide insights on the black hole physics, some of which we investigate analytically with the full Yang-Mills matrix model. In particular, we find many branches of saddle points which describe the known black hole solutions. We analytically construct the saddle points dual to the small black holes whose sizes are much smaller than the AdS radius. They include the asymptotically flat BMPV black holes embedded in large AdS with novel thermodynamic instabilities.

Physics↗

Tip‐Enhanced Imaging and Control of Infrared Strong Light‐Matter Interaction

Optical antenna resonators enable control of light‐matter interactions on the nano‐scale via electron–photon hybrid states in strong coupling. Specifically, mid‐infrared (MIR) nano‐antennas coupled to saturable intersubband transitions in multi‐quantum‐well (MQW) semiconductor heterostructures allow for the coupling strength to be tuned through antenna resonance and field intensity. Here, in this study, tip‐enhanced nano‐scale variation of antenna‐MQW coupling across the antenna is demonstrated, with a spatially‐dependent coupling strength $g_{\textrm{aq}}$ varying from 73 (strong coupling) to 24 cm -1 (weak coupling). This behavior is modeled based on the spatially dependent local constructive and destructive interference between tip and antenna fields. Using a quantum‐mechanical density‐matrix model of the MQW system with its designed values of transition dipole moment, doping density, and population decay time, the picosecond IR pulse coupling to intersubband transitions and the associated tip induced strong‐field saturation effects are described. These results present a new regime of nonlinear IR light‐matter control based on the dynamic manipulation of quantum hybrid states on the nanoscale and in the infrared, with a perspective regarding extension to molecular vibrations.

Wang, Yueying↗

A supersymmetric SYK model with a curious low energy behavior

We consider N = 2,4 supersymmetric SYK models that have a peculiar low energy behavior, with the entropy going like S = S 0 + (constant)T a , where a ≠ 1. The large N equations for these models are a generalization of equations that have been previously studied as an unjustified truncation of the planar diagrams describing the BFSS matrix quantum mechanics or other related matrix models. Here we reanalyze these equations in order to better understand the low energy physics of these models. We find that the scalar fields develop large expectation values which explore the low energy valleys in the potential. The low energy physics is dominated by quadratic fluctuations around these values. These models were previously conjectured to have a spin glass phase. We did not find any evidence for this phase by using the usual diagnostics, such as searching for replica symmetry breaking solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

R-matrix description of particle energy spectra produced by low-energy 3 H + 3 H reactions

An R-matrix model for three-body final states is presented and applied to a recent measurement of the neutron energy spectrum from the 3 H + 3 H→ 2n + α reaction. The calculation includes the n alpha and n n interactions in the final state, angular momentum conservation, antisymmetrization, and the interference between different channels. A good fit to the measured spectrum is obtained, where clear evidence for the 5 He ground state is observed. The model is also used to predict the alpha-particle spectrum from 3 H + 3 H as well as particle spectra from 3 He + 3 He. The R-matrix approach presented here is very general, and can be adapted to a wide variety of problems with three-body final states.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Entanglement and confinement in coupled quantum systems

We study some general properties of coupled quantum systems. We consider simple interactions between two copies of identical Hamiltonians such as the SYK model, Pauli spin chains with random magnetic field and harmonic oscillators. Such couplings make the ground states close to the thermofield double states of the uncoupled Hamiltonians. For the coupled SYK model, we push the numerical computation further towards the thermodynamic limit so that an extrapolation in the size of the system is possible. We find good agreement between the extrapolated numerical result and the analytic result in the large-q limit. We also consider the coupled gauged matrix model and vector model, and argue that the deconfinement is associated with the loss of the entanglement, similarly to the previous observation for the coupled SYK model. The understanding of the microscopic mechanism of the confinement/deconfinement transition enables us to estimate the quantum entanglement precisely, and backs up the dual gravity interpretation which relates the deconfinement to the disappearance of the wormhole. Our results demonstrate the importance of the entanglement between the color degrees of freedom in the emergence of the bulk geometry from quantum field theory via holography.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Lithium-Ion Battery Life Model with Electrode Cracking and Early-Life Break-in Processes

This paper develops a physically justified reduced-order capacity fade model from accelerated calendar- and cycle-aging data for 32 lithium-ion (Li-ion) graphite/nickel-manganese-cobalt (NMC) cells. The large data set reveals temperature-, charge C-rate-, depth-of-discharge-, and state of charge (SOC)-dependent degradation patterns that would be unobserved in a smaller test matrix. Model structure is informed by incremental capacity analysis that shows loss of lithium inventory and cathode-material loss as the dominant capacity fade mechanisms. The model includes terms attributable to solid-electrolyte interface (SEI) growth, electrode cracking, cycling-driven acceleration of SEI growth, and "break-in" mechanisms that slightly decrease or increase available Li inventory early in life. The study explores what mathematical couplings of these mechanisms best describe calendar aging, cycle aging, and mixed calendar/cycle aging. Various approaches are discussed for extracting relevant stress factors from complex cycling profiles to predict lifetime during real-world battery loads using models trained on constant-current laboratory test results. The complexity of the present human-driven model identification process motivates future work in machine learning to more widely search and statistically discern the optimal model that correctly extrapolates capacity fade based on physical knowledge.

25 ENERGY STORAGE↗

Sparsity-based photoacoustic image reconstruction with a linear array transducer and direct measurement of the forward model

Photoacoustic (PA) imaging is an emerging imaging technique for many clinical applications. One of the challenges posed by clinical translation is that imaging systems often rely on a finite-aperture transducer rather than a full tomography system. This results in imaging artifacts arising from an underdetermined reconstruction of the initial pressure distribution (IPD). Furthermore, clinical applications often require deep imaging, resulting in a low-signal-to-noise ratio for the acquired signal because of strong light attenuation in tissue. Conventional approaches to reconstruct the IPD, such as back projection and time-reversal, do not adequately suppress the artifacts and noise. We propose a sparsity-based optimization approach that improves the reconstruction of IPD in PA imaging with a linear array ultrasound transducer. In simulation studies, the forward model matrix was measured from k-Wave simulations, and the approach was applied to reconstruct simulated point objects and the Shepp–Logan phantom. The results were compared with the conventional back projection, time-reversal approach, frequency-domain reconstruction, and the iterative least-squares approaches. In experimental studies, the forward model of our imaging system is directly measured by scanning a graphite point source through the imaging field of view. Experimental images of graphite inclusions in tissue-mimicking phantoms are reconstructed and compared with the back projection and iterative least-squares approaches. Altogether these results show that our proposed optimization approach can leverage the sparsity of the PA images to improve the reconstruction of the IPD and outperform the existing popular reconstruction approaches.

47 OTHER INSTRUMENTATION↗

Effects of ionizing radiation and HPSE1 inhibition on the invasion of oral tongue carcinoma cells on human extracellular matrices in vitro

Highlights: • Myoma-based ECM models are suitable for in vitro chemoradiation assays. • Blocking of HPSE1 activity combined with IR may induce tongue cancer cell invasion. • Low doses of IR may increase MMP expression and invasion in tongue cancer cells. Chemoradiation is an established approach in the treatment of advanced oral tongue squamous cell carcinoma (OTSCC), but therapy may cause severe side-effects due to signal interchanges between carcinoma and the tumour microenvironment (TME). In this study, we examined the potential use of our human 3D myoma disc and Myogel models in in vitro chemoradiation studies by analysing the effects of ionizing radiation (IR) and the combined effect of heparanase I (HPSE1) inhibitors and IR on OTSCC cell proliferation, invasion and MMP-2 and − 9 production. Finally, we analysed the long-term effects of IR by studying clones of previously irradiated and invaded HSC-3 cells. We found that in both human uterine leiomyoma-based extracellular matrix models IR inhibited the invasion of HSC-3 cells, but blocking HPSE1 activity combined with IR induced their invasion. Low doses of IR increased MMP expression and initiated epithelial-mesenchymal transition in cells cultured on myoma discs. We conclude that myoma models offer consistent methods for testing human carcinoma cell invasion and phenotypic changes during chemoradiation treatment. In addition, we showed that IR had long-term effects on MMP-2 and − 9, which might elicit different HSC-3 invasion responses when cells were under the challenge of HPSE1 inhibitors and IR.

60 APPLIED LIFE SCIENCES↗

Woven ceramic matrix composite surrogate model based on physics-informed recurrent neural network

A recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear constitutive behavior of woven ceramic matrix composites (CMCs) driven by matrix damage at multiple length scales. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the multiscale generalized method of cells (MSGMC) approach coupled with a matrix damage model. This coupling permits simulating the nonlinear behavior of woven CMCs based on constituent response at the micro-, meso-, and macroscales. The multiscale repeating unit cell is loaded under non-monotonic conditions including multiple load / unload cycles and tension / compression. The fiber volume fraction as well as the intra- and intertow void volume fractions are also varied in the generation of training data. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input strain sequence and fiber and void volume fractions, the resulting stress versus strain response while satisfying physical constraints such as positive semi-definiteness of the tangent stiffness matrix and linear elastic unloading. Further, the trained surrogate model effectively matches the stress versus strain response and successfully predicts the tangent modulus throughout the loading regime. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex architectures, nonlinear multiaxial material response, and under non-monotonic loading conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

EXAFS Spectroscopy of Fractional Mixed-Valence Charge Density Wave Systems

We present extended X-ray absorption fine structure (EXAFS) spectra and modeling of a series of structurally tunable quasi-one-dimensional mixed-valence platinum–halide linear chain materials, [Pt(en 2 )][Pt(en 2 )X 2 ](ClO 4 ) 4 with X = Cl, Br, I. The materials exhibit a commensurate charge density wave with fractional charge states on alternating platinum ions in the chain, as well as a Peierls distortion with alternating platinum–halide bond lengths. The amplitude of the charge density wave and, correspondingly, the extent of the Peierls distortion are controlled by the identity of the bridging halide ion. We have carried out ab initio multiple scattering calculations using the FEFF9 code to relate the oriented Pt L III EXAFS spectra to the tunable electronic and structural properties. The spectral modeling reveals distinct photoelectron threshold energy values for the two inequivalent platinum ions in each of the mixed-valence chains, with values that vary systematically with fractional valence state. The difference in the photoelectron threshold energies of the two inequivalent platinum ions within each material correlates directly with the amplitude of the charge density wave, reflecting the decrease in charge density wave strength through the halide series X = Cl, Br, and I. We use dynamical matrix modeling to relate the experimentally determined mean-square relative displacement parameters for the metal–halide bond distances to the chain–axis vibrational modes that modulate the charge density wave structure. In addition, we discuss the EXAFS fitting results for the Pt–I bond lengths in the [Pt(en 2 )][Pt(en 2 )I 2 ](ClO 4 ) 4 complex in comparison to previous, mutually inconsistent structural determinations for this material.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Windowed multipole representation of R -matrix cross sections

Nuclear cross sections are basic inputs to any nuclear computation. Campaigns of experiments are fitted with the parametric R-matrix model of quantum nuclear interactions, and the resulting cross sections are documented—both pointwise and as resonance parameters (with uncertainties)—in standard evaluated nuclear data libraries (ENDF, JEFF, BROND, JENDL, CENDL, TENDL): these constitute our common knowledge of fundamental low-energy nuclear cross sections. In the past decade, a collaborative effort has been deployed to establish a new nuclear cross-section library format—the Windowed Multipole Library—with the goal of considerably reducing the computational cost of cross-section calculations in nuclear transport simulations. This work lays the theoretical foundations underpinning these efforts. From general R-matrix scattering theory, we derive the windowed multipole representation of nuclear cross sections. Though physically and mathematically equivalent to R-matrix cross sections, the windowed multipole representation is particularly well suited for subsequent temperature treatment of angle-integrated cross sections, in particular Doppler broadening, which is the averaging of cross sections over the thermal motion of the target atoms. Doppler broadening is of critical importance in neutron transport applications, as it ensures the stability of many nuclear reactors (negative thermal reactivity). Yet, Doppler broadening of nuclear cross sections has been a considerable bottleneck for nuclear transport computations, often requiring memory-costly pretabulations. We show that the windowed multipole representation can perform accurate Doppler broadening analytically (up to the first reaction threshold), from which we derive cross-section temperature derivatives to any order—all computable on the fly (without precalculations stored in memory). Furthermore, we here establish a way of converting the R-matrix resonance parameters uncertainty (covariance matrices) into windowed multipole parameters uncertainty. We show that generating stochastic nuclear cross sections by sampling from the resulting windowed multipole covariance matrix can reproduce the cross-section uncertainty in the original nuclear data file. The windowed multipole representation is therefore a novel nuclear physics formalism able to generate Doppler broadened stochastic nuclear cross sections on the fly, unlocking breakthrough computational gains for nuclear computations. Through this foundational paper, we hope to make the windowed multipole representation accessible, reproducible, and usable for the nuclear physics community, as well as provide the theoretical basis for future research on expanding its capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Strong coupling expansion of free energy and BPS Wilson loop in $$ \mathcal{N} $$ = 2 superconformal models with fundamental hypermultiplets

As a continuation of the study (in arXiv:2102.07696 and arXiv:2104.12625 ) of strong-coupling expansion of non-planar corrections in \( \mathcal{N} \) = 2 4d superconformal models we consider two special theories with gauge groups SU( N ) and Sp(2 N ). They contain N -independent numbers of hypermultiplets in rank 2 antisymmetric and fundamental representations and are planar-equivalent to the corresponding \( \mathcal{N} \) = 4 SYM theories. These \( \mathcal{N} \) = 2 theories can be realised on a system of N D3-branes with a finite number of D7-branes and O7-plane; the dual string theories should be particular orientifolds of AdS 5 × S 5 superstring. Starting with the localization matrix model representation for the \( \mathcal{N} \) = 2 partition function on S 4 we find exact differential relations between the 1 /N terms in the corresponding free energy F and the \( \frac{1}{2} \) -BPS Wilson loop expectation value \( \left\langle \mathcal{W}\right\rangle \) and also compute their large ’t Hooft coupling ( λ » 1) expansions. The structure of these expansions is different from the previously studied models without fundamental hypermultiplets. In the more tractable Sp(2 N ) case we find an exact resummed expression for the leading strong coupling terms at each order in the 1 /N expansion. We also determine the exponentially suppressed at large λ contributions to the non-planar corrections to F and \( \left\langle \mathcal{W}\right\rangle \) and comment on their resurgence properties. We discuss dual string theory interpretation of these strong coupling expansions.

1/N expansion↗