Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Functional coupling”

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,387 records · Page 77

Incorporation of Damage and Failure into an Orthotropic Elasto-Plastic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

A material model which incorporates several key capabilities which have been identified by the aerospace community as lacking in the composite impact models currently available in LS-DYNA(Registered Trademark) is under development. In particular, the material model, which is being implemented as MAT 213 into a tailored version of LS-DYNA being jointly developed by the FAA and NASA, incorporates both plasticity and damage within the material model, utilizes experimentally based tabulated input to define the evolution of plasticity and damage as opposed to specifying discrete input parameters (such as modulus and strength), and is able to analyze the response of composites composed with a variety of fiber architectures. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. The capability to account for the rate and temperature dependent deformation response of composites has also been incorporated into the material model. For the damage model, a strain equivalent formulation is utilized to allow for the uncoupling of the deformation and damage analyses. In the damage model, a diagonal damage tensor is defined to account for the directionally dependent variation of damage. However, in composites it has been found that loading in one direction can lead to damage in multiple coordinate directions. To account for this phenomena, the terms in the damage matrix are semi-coupled such that the damage in a particular coordinate direction is a function of the stresses and plastic strains in all of the coordinate directions. The onset of material failure, and thus element deletion, is being developed to be a function of the stresses and plastic strains in the various coordinate directions. Systematic procedures are being developed to generate the required input parameters based on the results of experimental tests.

Finite Element Method↗

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Replica symmetry breaking for the integrable two-site Sachdev–Ye–Kitaev model

We analyze a two-body non-Hermitian two-site Sachdev–Ye–Kitaev (SYK) model with the couplings of one site complex conjugated to the other site. This model, with no explicit coupling between the sites, shows an infinite number of second-order phase transitions, which is a consequence of the factorization of the partition function into a product over Matsubara frequencies. We calculate the quenched free energy in two different ways: first in terms of the single-particle energies and second by solving the Schwinger–Dyson equations of the two-site model. The first calculation can be done entirely in terms of a one-site model. The conjugate replica enters due to non-analyticities when Matsubara frequencies enter the spectral support of the coupling matrix. The second calculation is based on the replica trick of the two-site partition function. Both methods give the same result. The free-fermion partition function can also be rephrased as a matrix model for the coupling matrix. Up to minor details, this model is the random matrix model that describes the chiral phase transition of QCD, and the order parameter of the two-body model corresponds to the chiral condensate of QCD. Comparing to the corresponding four-body model, we are able to determine which features of the free energy are due to the chaotic nature of the four-body model. The high-temperature phase of both models is entropy dominated, and in both cases, the free energy is determined by the spectral density. The chaotic four-body SYK model has a low-temperature phase whose free energy is almost temperature-independent, signaling an effective gap of the theory even though the actual spectrum does not exhibit a gap. On the other hand, the low-temperature free energy of the two-body SYK model is not flat; in fact, it oscillates to arbitrarily low temperature. This indicates a less desirable feature that the entropy of the two-body model is not always positive in the low-temperature phase, which most likely is a consequence of the non-hermiticity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Real-Time Testbed for Smart Grid Recloser Controller

The growing need for a low voltage recloser has become apparent due to the rise in requirements for a smart grid. This includes more detailed management of power flow forward (towards load) and backward (towards generation), source synchronization in real time, more indepth fault responses, and the use of green energy. The SEL-651R-2 relay is a device that can manage these needs, especially in fault response and synchronization, and is commonly used in systems called microgrids. Microgrids are distribution level systems that are able to operate separated from the main grid, are typically installed much closer to the load(s), and are fed by distributed energy resources (DERs), such as wind, solar or diesel generators. The SEL-651R-2 is normally used in the field with presets operative settings, but the Western Michigan University (WMU) Center for Interdisciplinary Research on Secure, Efficient and Sustainable Energy Technology (WMU InterEnergy Center) wished to test this device in its range of capabilities for microgrid application. A Hardware-In-the-Loop (HIL) testbed was implemented and used through the Real Time Digital Simulator (RTDS) using the RSCAD software to test the SEL-651R-2's use cases and functions. The testbed includes a microgrid with interconnection to a larger main grid, and the relay is meant to control the recloser at the point of common coupling (PCC) between the main grid and microgrid. The testbed shows how basic protections, reclosing, and synchronization checks function when handling faults that affect both the microgrid and the main grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Navier-Stokes calculations with a coupled strongly implicit method. Part 2: Spline solutions

A coupled strongly implicit method is combined with a deferred-corrector spline solver for the vorticity-stream function form of the Navier-Stokes equation. Solutions for cavity, channel and cylinder flows are obtained with the fourth-order spline 4 procedure. The strongly coupled spline corrector method converges as rapidly as the finite difference calculations and also allows for arbitrary large time increments for the Reynolds numbers considered. In some cases fourth-order smoothing or filtering is required in order to suppress high frequency oscillations.

Rubin, S. G.↗

Magnon-Magnon Interaction Induced by Dynamic Coupling in a Hybrid Magnonic Crystal

We report a combined experimental and numerical investigation of spin-wave dynamics in a hybrid magnonic crystal consisting of a CoFeB artificial spin ice (ASI) of stadium-shaped nanoelements patterned atop a continuous NiFe film separated by a 5 nm Al 2 O 3 spacer. Using Brillouin light scattering spectroscopy, we probe the frequency dependence of thermal spin waves as functions of applied magnetic field and wavevector, revealing the decisive role of interlayer dipolar coupling in the magnetization dynamics. Micromagnetic simulations complement the experiments, showing a strong interplay between ASI edge modes and backward volume modes in the NiFe film. The contrast in saturation magnetization between CoFeB and NiFe enhances this coupling, leading to a pronounced hybridization manifested as a triplet of peaks in the BLS spectra predicted by simulations and observed experimentally. This magnon−magnon coupling persists over a wide magnetic field range, shaping both the spin-wave dispersion at fixed fields and the full frequency-field response throughout the magnetic hysteresis loop. Our findings establish how ASI geometry can selectively enhance specific spin-wave wavelengths in the underlying film, thereby boosting their amplitude and identifying them as preferential channels for spin wave transmission and manipulation.

Brillouin light scattering spectroscopy↗

AMReX and pyAMReX: Looking beyond the exascale computing project

AMReX is a software framework for the development of block-structured mesh applications with adaptive mesh refinement (AMR). AMReX was initially developed and supported by the AMReX Co-Design Center as part of the U.S. DOE Exascale Computing Project (ECP), and is continuing to grow post-ECP. In addition to adding new functionality and performance improvements to the core AMReX framework, we have also developed a Python binding, pyAMReX, that provides a bridge between AMReX-based application codes and the data science ecosystem. pyAMReX provides zero-copy application GPU data access for AI/ML, in situ analysis and application coupling, and enables rapid, massively parallel prototyping. In this paper we review the overall functionality of AMReX and pyAMReX, focusing on new developments, new functionality, and optimizations of key operations. We also summarize capabilities of ECP projects that used AMReX and provide an overview of new, non-ECP applications.

Myers, Andrew↗

Oxygen Vacancy Injection as a Pathway to Enhancing Electromechanical Response in Ferroelectrics

Since their discovery in late 1940s, perovskite ferroelectric materials have become one of the central objects of condensed matter physics and materials science due to the broad spectrum of functional behaviors they exhibit, including electro-optical phenomena and strong electromechanical coupling. In such disordered materials, the static properties of defects such as oxygen vacancies are well explored but the dynamic effects are less understood. In this work, the first observation of enhanced electromechanical response in BaTiO3 thin films is reported driven via dynamic local oxygen vacancy control in piezoresponse force microscopy (PFM). A persistence in peizoelectricity past the bulk Curie temperature and an enhanced electromechanical response due to a created internal electric field that further enhances the intrinsic electrostriction are explicitly demonstrated. The findings are supported by a series of temperature dependent band excitation PFM in ultrahigh vacuum and a combination of modeling techniques including finite element modeling, reactive force field, and density functional theory. Furthermore, this study shows the pivotal role that dynamics of vacancies in complex oxides can play in determining functional properties and thus provides a new route toward– achieving enhanced ferroic response with higher functional temperature windows in ferroelectrics and other ferroic materials.

36 MATERIALS SCIENCE↗

Spin pumping from antiferromagnetic insulator spin-orbit-proximitized by adjacent heavy metal: a first-principles Floquet-nonequilibrium Green function study

Motivated by the recent experiment on spin pumping from sub-THz radiation-driven uniaxial antiferromagnetic insulator (AFI) MnF 2 into heavy metal (HM) Pt hosting strong spin-orbit (SO) coupling, we compute and compare pumped spin currents in Cu/MnF 2 /Cu and Pt/MnF 2 /Cu heterostructures. Recent theories of spin pumping by AFI have relied on simplistic Hamiltonians (such as tight-binding) and the scattering approach to quantum transport yielding the so-called interfacial spin mixing conductance (SMC), but the concept of SMC ceases to be applicable when SO coupling is present directly at the interface. In contrast, we use a more general first-principles quantum transport approach which combines noncollinear density functional theory with Floquet-nonequilibrium Green's functions in order to take into account: SO-proximitized AFI as a new type of quantum material, different from isolated AFI and brought about by AFI hybridization with adjacent HM layer; strong SO coupling at the interface; and evanescent wavefunctions penetrating from Pt or Cu into AFI layer to make its interfacial region conducting rather than insulating as in the isolated AFI. The DC component of pumped spin current $I_\mathrm{DC}^{S_z}$ vs. precession cone angle $\theta_{\boldsymbol{l}}$ of the Néel vector l of AFI does not follow putative $I^{S_z}_\mathrm{DC} \propto \sin^2 \theta_{\boldsymbol{l}}$, except for very small angles $\theta_{\boldsymbol{l}} \lesssim 10^\circ$ for which we define an effective SMC from the prefactor and find that it doubles from MnF2/Cu to MnF2/Pt interface. In addition, the angular dependence $I^{S_z}_\mathrm{DC}(\theta_{\boldsymbol{l}})$ differs for opposite directions of precession of the Néel vector, leading to twice as large SMC for the right-handed than for the left-handed chirality of the precession mode.

36 MATERIALS SCIENCE↗

The QCD Adler Function and the Muon g − 2 Anomaly from Renormalons

We describe the Adler function in Quantum Chromodynamics using a transseries representation within a resurgent framework. The approach is based on a Borel-Ecalle resummation of the infrared renormalons combined with an effective running for the strong coupling. The new approach is flexible enough to give values in agreement with the current Adler function determinations. We then apply our finding to the muon’s anomalous magnetic moment studying the possibility of saturating, solely in terms of the vacuum polarization function, the current discrepancy between the best Standard Model value for the muon’s anomalous magnetic moment and the experimental value obtained by the most recent muon g−2 collaboration. The latter shows that the Adler function’s new representation can also be consistent with recent lattice determinations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electronic structure and energetics of a heterodimeric BChl g ′/Chl a ′ special pair generated by exposure of Heliomicrobium modesticaldum to dioxygen

Heliobacteria are anoxygenic phototrophs that have a Type I homodimeric reaction center containing bacteriochlorophyll g (BChl g). Previous experimental studies have shown that in the presence of light and dioxygen, BChl g is converted into 8 1 -OH-chlorophyll a F (hereafter Chl a F ), with an accompanying loss of light-driven charge separation. These studies suggest that the reaction center only loses the ability to transfer electrons once both BChl g′ molecules of the P 800 special pair have been converted to Chl a F ′. Here, the present work confirms that the partially converted BChl g′/Chl a F ′ special pair remains functional in samples exposed to dioxygen by demonstrating its presence using hyperfine couplings obtained from Q-band 1 H ENDOR, 2D 14 N HYSCORE and DFT methods. The DFT calculations of the BChl g′/BChl g′ homodimeric primary donor, which are based on the recently published X-ray crystal structure, predict that the unpaired electron spin is equally delocalized over both BChl g′ molecules and provide an excellent match to the experimental hyperfine couplings of the anaerobic samples. Exposure to dioxygen leads to substantial changes in the hyperfine interactions, indicative of greater localization of the unpaired electron spin. The measured hyperfine couplings are reproduced in the DFT calculations by replacing one of the BChl g′ molecules of the primary donor with a Chl a F ′ molecule. The calculations reveal that the spin density becomes localized on BChl g′ in the heterodimeric primary donor. Time-dependent DFT calculations demonstrate that conversion of either or both of the accessory BChl g molecules and/or one of the BChl g′ molecules of P800 to Chl a F ′ results in minor effects on the energy of the charge-separated states. In contrast, if both of the BChl g′ molecules of P 800 are converted a large increase in the energy of the charge-separated state occurs. This suggests that the reaction center remains functional when only one half of the dimer is converted, however, conversion of both halves of the P 800 dimer leads to loss of function

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of spin-orbit coupling and very large supercells on the description of acceptors in CdTe

Predicting accurate shallow donor and acceptor levels in semiconductors has been quite challenging using periodic boundary conditions as implemented in current density functional theory codes. The reason is that the wave functions associated with the shallow centers are quite extended and are not fully contained in the supercell (typically with a few hundred atoms) so impurities in the periodically repeated image cells interact with each other. Errors of ~0.1 eV are expected, and these are of the same order of magnitude as the ionization energies themselves. In the case of acceptors in CdTe, this problem is exacerbated by the strong spin-orbit coupling that split the Te-related states at the top of the valence band, making the calculations at least 8 times more expensive. CdTe is an important solar-cell material with record high efficiency of 22%. One of the main limiting factors to increasing the efficiency towards the theoretical limit of ~30% is the often reported very low hole concentration. Pushing to the limit of computational capability by using very large supercells with spin-orbit coupling we report the results of hybrid functional calculations of group-V acceptors in CdTe. We show that extrapolation to the dilute limit leads to an interpretation of the experimental data that is qualitatively different from previous DFT and hybrid functional calculation reports. We find that the group-V impurities indeed behave as shallow acceptors and that the corresponding compensating AX-centers are unstable and do not limit p-type doping. We address the differences between our results and previous theoretical predictions and show that our calculated ionization energies predict hole concentrations that are in excellent agreement with recent temperature-dependent Hall measurements on high-quality single-crystal samples.

14 SOLAR ENERGY↗

Impact of symmetry breaking and spin-orbit coupling on the band gap of halide perovskites

Halide perovskite (HP) materials have recently emerged as a class of semiconductors with immense promise for various optoelectronic applications, ranging from solar cells to light-emitting diodes. One of the unique attributes of HPs is their tunable band gaps with different factors governing their value. Furthermore, the first factor is related to relativistic corrections [“mass-Darwin,” connected to the ns 2 lone pairs, and spin-orbit coupling (SOC)] that induce an orbital shift or degeneracy splitting, resulting in a band-gap reduction. The second factor involves the structural configuration: in HPs the local symmetry of each Wyckoff position tends to be broken, inducing an opening of the band gap. Based on high-throughput density functional theory calculations, this paper systematically studies a possible self-cancelation on the band-gap correction for HPs when the polymorphous configuration—structural effects—and the SOC—electronic effects—are included. Our results indicate that the nature of interplay between SOC and symmetry breaking (SB) is that they are independent decoupling effects to describe the band-gap magnitude in halide perovskites. As a result of that, we observe a transitivity of the band-gap description; i.e., if we know the band gap of halide perovskites without SB and SOC, we can independently add the effects of band-gap reduction due to SOC and band-gap opening due to SB, regardless of the order in which these effects are considered.

14 SOLAR ENERGY↗

A continuum framework for coupled solid deformation–fluid flow through anisotropic elastoplastic porous media

We present a continuum framework for coupled solid deformation–fluid flow in anisotropic elastoplastic porous media. A thermodynamic formulation of the coupled processes gives rise to an anisotropic Biot tensor that is a function of drained elastic tangent moduli tensor of the solid skeleton and the intrinsic bulk modulus of the solid constituent. Two effective stress measures emerge from the formulation, namely, σ', which is energy-conjugate to the elastic strain, and σ'', which is energy-conjugate to the plastic strain. For the special case of transverse isotropy that is commonly encountered in natural rocks, the Biot tensor can be expressed in terms of its normal and tangential components to the bedding plane, along with a microstructure tensor. Apart from its thermodynamic consistency, an advantage of this new formulation is that standard mixed finite element formulation can be employed to discretize the domain and solve initial boundary-value problems. We conduct plane strain simulations of coupled solid deformation–fluid flow in a transversely isotropic porous medium to demonstrate the impacts of material anisotropy, stress history, and the Biot tensor on the system response.

42 ENGINEERING↗

SU(3) Topology of Magnon-Phonon Hybridization in 2D Antiferromagnets

Magnon-phonon hybrid excitations are studied theoretically in a two-dimensional antiferromagnet with an easy axis normal to the plane. In this work, we show that two magnon bands and one phonon band are intertwined by the magnetoelastic coupling through a nontrivial SU(3) topology, which can be intuitively perceived by identifying a skyrmion structure in the momentum space. Our results are insensitive to lattice details and generally applicable to two-dimensional antiferromagnets. We show this by developing a continuum theory as the long-wavelength approximation to the tight-binding model. The theoretical results can be probed by measuring the thermal Hall conductance as a function of the temperature and the magnetic field. We envision that the magnetoelastic coupling in antiferromagnets can be a promising venue in search of various topological excitations, which cannot be found in magnetic or elastic models alone.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Data pattern sensitivity in tracking performance of an AC coupled Costas loop with hard-limited in-phase channel

This paper is concerned with data pattern sensitivity in carrier tracking performance of an AC coupled Costas loop with a suppressed BPSK signal. The signal amplitude suppression factor is derived as a function of data 'asymmetry ratio' - the ratio of '1's to the total number of bits in a period of a periodic signal. For an asymmetric pattern, the effect of AC coupling is noticeable whereas there is almost no effect for symmetric square wave. The tracking performance with an asymmetric pattern is worse than that with a symmetric pattern. However, it is also shown that as expected, the tracking performance of a DC coupled loop with an asymmetric pattern is better than that with a symmetric pattern.

Park, Y. H.↗

Reduced‐Order Modeling of Energetic Materials Using Physics‐Aware Recurrent Convolutional Neural Networks in a Latent Space (LatentPARC)

Physics-aware deep learning (PADL) has gained popularity for use in spatiotemporal dynamics simulations, such as those in computational modeling of energetic materials (EM). We show that the challenge PADL methods face while learning complex field evolution problems can be simplified and accelerated by decoupling it into two tasks: learning complex geometric features in evolving fields and modeling dynamics over these features in a lower-dimensional feature space. We build upon our previous work on physics-aware recurrent convolutional neural networks (PARC). PARC embeds knowledge of underlying physics into its neural network architecture for more robust and accurate prediction of evolving physical fields. PARC was shown to effectively learn complex nonlinear features such as the formation of hotspots and coupled shock fronts in various initiation scenarios of EMs, as a function of microstructures, serving effectively as a microstructure-aware burn model. Here, we further accelerate PARC and reduce its computational cost by projecting the original dynamics onto a lower-dimensional invariant manifold, or “latent space.” The projected latent representation encodes the complex geometry of evolving fields (e.g., temperature and pressure) in a set of data-driven features. The reduced dimension of this latent space allows us to learn the dynamics during the initiation of EM with a lighter and more efficient model. We observe a significant decrease in training and inference time while maintaining results comparable to PARC at inference. This work takes steps towards enabling rapid prediction of EM thermomechanics at larger scales and characterization of EM structure–property–performance linkages at a full application scale.

Mathematics and Computing↗

Loops and trees in generic EFTs

We consider aspects of tree and one-loop behavior in a generic 4d EFT of massless scalars, fermions, and vectors, with a particular eye to the high-energy limit of the Standard Model EFT at operator dimensions 6 and 8. First, we classify the possible Lorentz structures of operators and the subset of these that can arise at tree-level in a weakly coupled UV completion, extending the tree/loop classification through dimension 8 using functional methods. Second, we investigate how operators contribute to tree and one-loop helicity amplitudes, exploring the impact of non-renormalization theorems through dimension 8. We further observe that many dimension 6 contributions to helicity amplitudes, including rational parts, vanish exactly at one-loop level. This suggests the impact of helicity selection rules extends beyond one loop in non-supersymmetric EFTs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗