Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential inversion”

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 469 records · Page 26

Electron Cyclotron Emission–Based Separatrix Identification in ITER with OMFIT Synthetic Modeling

Accurate determination of the separatrix location is essential for understanding edge plasma behavior and optimizing confinement in tokamaks, especially in next-generation devices such as ITER. In this study, a synthetic microwave diagnostics module was developed and implemented in the OMFIT framework to assess the feasibility of an electron cyclotron emission–based separatrix detection method in ITER plasmas. Simulations were carried out using ITER H-mode equilibrium scenarios with different plasma density profiles and different pedestal widths. Here, the results show that the electron emission temperature profiles consistently exhibit an inversion pattern near the edge, with a well-defined minimum point that could serve as a proxy for the separatrix location. However, unlike in the DIII-D, the minimum point in ITER is systematically offset by approximately 2 cm into the scrape-off layer, independent of density or pedestal width, which is within the radial resolution range (2 to 5 cm) determined by the 500-MHz channel spacing. While the method does not provide the exact separatrix location, it offers a reliable indicator of the boundary region and has potential applications for real-time boundary monitoring in ITER and other future fusion devices.

Electron cyclotron emission↗

Topological features in the ferromagnetic Weyl semimetal CeAlSi: Role of domain walls

In the ferromagnetic (FM) Weyl semimetal CeAlSi both space-inversion and time-reversal symmetries are broken. Our quantum oscillation (QO) data indicate that the FM ordering modifies the Fermi-surface topology and also leads to an unusual drop in the QO amplitude. In the FM phase, we find a pressure-induced suppression of the anomalous and the loop Hall effects. This cannot be explained based on the electronic band structure or magnetic structure, both of which are nearly pressure independent. Instead, we show that a simplified model describing the scattering of Weyl fermions off FM domain walls can potentially explain the observed topological features. Our study highlights the importance of domain walls for understanding transport in FM Weyl semimetals.

36 MATERIALS SCIENCE↗

TRINIDI: Time-of-Flight Resonance Imaging With Neutrons for Isotopic Density Inference

Accurate reconstruction of 2D and 3D isotope densities is a desired capability with great potential impact in applications such as evaluation and development of next-generation nuclear fuels. Neutron time-of-flight (TOF) resonance imaging offers a potential approach by exploiting the characteristic neutron absorption spectra of each isotope. However, it is a major challenge to compute quantitatively accurate images due to a variety of confounding effects such as severe Poisson noise, background scatter, beam non-uniformity, absorption non-linearity, and extended source pulse duration. We present the TRINIDI algorithm which is based on a two-step process in which we first estimate the neutron flux and background counts, and then reconstruct the areal densities of each isotope and pixel. Both components are based on the inversion of a forward model that accounts for the highly non-linear absorption, energy-dependent emission profile, and Poisson noise, while also modeling the substantial spatio-temporal variation of the background and flux. Further, to do this, we formulate the non-linear inverse problem as two optimization problems that are solved in sequence. We demonstrate on both synthetic and measured data that TRINIDI can reconstruct quantitatively accurate 2D views of isotopic areal density that can then be reconstructed into quantitatively accurate 3D volumes of isotopic volumetric density.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Landscape Topography and Regional Drought Alters Dust Microbiomes in the Sierra Nevada of California

Dust provides an ecologically significant input of nutrients, especially in slowly eroding ecosystems where chemical weathering intensity limits nutrient inputs from underlying bedrock. In addition to nutrient inputs, incoming dust is a vector for dispersing dust-associated microorganisms. While little is known about dust-microbial dispersal, dust deposits may have transformative effects on ecosystems far from where the dust was emitted. Using molecular analyses, we examined spatiotemporal variation in incoming dust microbiomes along an elevational gradient within the Sierra Nevada of California. We sampled throughout two dry seasons and found that dust microbiomes differed by elevation across two summer dry seasons (2014 and 2015), which corresponded to competing droughts in dust source areas. Dust microbial taxa richness decreased with elevation and was inversely proportional to dust heterogeneity. Likewise, dust phosphorus content increased with elevation. At lower elevations, early season dust microbiomes were more diverse than those found later in the year. The relative abundances of microbial groups shifted during the summer dry season. Furthermore, mutualistic fungal diversity increased with elevation, which may have corresponded with the biogeography of their plant hosts. Although dust fungal pathogen diversity was equivalent across elevations, elevation and sampling month interactions for the relative abundance, diversity, and richness of fungal pathogens suggest that these pathogens differed temporally across elevations, with potential implications for humans and wildlife. This study shows that landscape topography and droughts in source locations may alter the composition and diversity of ecologically relevant dust-associated microorganisms.

59 BASIC BIOLOGICAL SCIENCES↗

Evolution of artificial intelligence for application in contemporary materials science

Abstract Contemporary materials science has seen an increasing application of various artificial intelligence techniques in an attempt to accelerate the materials discovery process using forward modeling for predictive analysis and inverse modeling for optimization and design. Over the last decade or so, the increasing availability of computational power and large materials datasets has led to a continuous evolution in the complexity of the techniques used to advance the frontier. In this Review, we provide a high-level overview of the evolution of artificial intelligence in contemporary materials science for the task of materials property prediction in forward modeling. Each stage of evolution is accompanied by an outline of some of the commonly used methodologies and applications. We conclude the work by providing potential future ideas for further development of artificial intelligence in materials science to facilitate the discovery, design, and deployment workflow. Graphical abstract

Materials Science↗

Direct cross-section measurement of the weak r -process Sr 88 ( α , n ) Zr 91 reaction in ν -driven winds of core-collapse supernovae

About half of the heavy elements beyond iron are known to be produced by the rapid neutron capture process, known as the r process. However, the astrophysical site producing the r process is still uncertain. Chemical abundances observed in several cosmic sites indicate that different mechanisms should be at play. For instance, the abundances around silver measured in a subset of metal-poor stars indicate the presence of a weak r process. This process may be active in neutrino-driven winds of core collapse supernovae where (α,n) reactions dominate the synthesis of Ζ ≈ 40 elements in the expelled materials. Scarcely measured, the rates of (α,n) reactions are determined from statistical Hauser-Feshbach calculations with α-optical-model potentials, which are still poorly constrained. Further, the uncertainties of the (α,n) reaction rates therefore make a significant contribution to the uncertainties of the abundances determined from stellar modeling. In this work, the 88 Sr (α,n)⁢ 91 Zr reaction which impacts the weak r-process abundances has been probed at astrophysics energy for the first time; directly measuring the total cross sections at astrophysical energies of 8.37–13.09 MeV in the center of mass (3.8–7.5 GK). Two measurements were performed at ATLAS with the electrically segmented ionization chamber MUSIC, in inverse kinematics, while following the active target technique. The cross sections of this α-induced reaction on 88 Sr, located at the shell closure N = 50, have been found to be lower than expected, by a factor of 3, despite recent statistical calculations validated by measurements on neighboring nuclei. This result encourages more experimental investigations of (α,n) reactions, at Ν = 50 and towards the neutron-rich side, to further test the predictive power and reliability of such calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Towards Construction of a Novel Nanometer-Resolution MeV-STEM for Imaging Thick Frozen Biological Samples

Driven by life-science applications, a mega-electron-volt Scanning Transmission Electron Microscope (MeV-STEM) has been proposed here to image thick frozen biological samples as a conventional Transmission Electron Microscope (TEM) may not be suitable to image samples thicker than 300–500 nm and various volume electron microscopy (EM) techniques either suffering from low resolution, or low speed. The high penetration of inelastic scattering signals of MeV electrons could make the MeV-STEM an appropriate microscope for biological samples as thick as 10 μm or more with a nanoscale resolution, considering the effect of electron energy, beam broadening, and low-dose limit on resolution. The best resolution is inversely related to the sample thickness and changes from 6 nm to 24 nm when the sample thickness increases from 1 μm to 10 μm. To achieve such a resolution in STEM, the imaging electrons must be focused on the specimen with a nm size and an mrad semi-convergence angle. This requires an electron beam emittance of a few picometers, which is ~1000 times smaller than the presently achieved nm emittance, in conjunction with less than 10 -4 energy spread and 1 nA current. We numerically simulated two different approaches that are potentially applicable to build a compact MeV-STEM instrument: (1) DC (Direct Current) gun, aperture, superconducting radio-frequency (SRF) cavities, and STEM column; (2) SRF gun, aperture, SRF cavities, and STEM column. Beam dynamic simulations show promising results, which meet the needs of an MeV-STEM, a few-picometer emittance, less than 10 -4 energy spread, and 0.1–1 nA current from both options. Also, we designed a compact STEM column based on permanent quadrupole quintuplet, not only to demagnify the beam size from 1 μm at the source point to 2 nm at the specimen but also to provide the freedom of changing the magnifications at the specimen and a scanning system to raster the electron beam across the sample with a step size of 2 nm and the repetition rate of 1 MHz. This makes it possible to build a compact MeV-STEM and use it to study thick, large-volume samples in cell biology.

36 MATERIALS SCIENCE↗

Radiation Source Localization Using Surrogate Models Constructed from 3-D Monte Carlo Transport Physics Simulations

Recent research has focused on the development of surrogate models for radiation source localization in a simulated urban domain. We employ the Monte Carlo N-Particle (MCNP) code to provide high- delity simulations of radiation transport within an urban domain. The model is constructed to employ a source location (x, y, z) as input and return the estimated count rate for a set of speci ed detector locations. Because MCNP simulations are computationally expensive, we develop e cient and accurate surrogate models of the detector responses. We construct surrogate models using Gaussian processes (GP) and neural networks (NN) that we train and verify using the MCNP simulations. The trained surrogate models provide an e cient framework for Bayesian inference and experimental design. We employ Delayed Rejection Adaptive Metropolis (DRAM), a Markov Chain Monte Carlo (MCMC) algorithm, to infer the location and intensity of an unknown source. The DRAM results yield a posterior probability distribution for the source's location conditioned on the observed detector count rates. The posterior distribution exhibits regions of high and low probability within the simulated environment identifying potential source locations. In this manner, we can quantify the source location to within at least one of these regions of high probability in the considered cases. Employing these methods, we are able to reduce the space of potential source locations by at least 60%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Broadband Fourier-Transform-Detected EPR at W-Band

This manuscript reports development of wideband Fourier transform (FT) electron paramagnetic resonance (EPR) detection capability at the uniquely high frequency of 94 GHz. Here, this technique is based on the HiPER spectrometer developed at the University of St. Andrews, into which we have integrated an arbitrary waveform generator (AWG) that is used to modulate the output from a solid-state multiplier chain prior to amplification, generating up to 1 kW microwave power with 1 GHz bandwidth. Benchmark experiments are presented for a standard TEMPOL radical, which comprises a 500 MHz broad EPR spectrum at 94 GHz. Using a single frequency-swept microwave pulse, efficient inversion of this spectrum is demonstrated, enabling frequency-dependent studies of the magnetization relaxation in the time-domain, again via broadband echo detection. From these measurements, an anisotropy in the spin–lattice relaxation time, T 1 , can be determined with ease for TEMPOL. In addition, we implement the FT detection scheme for multi-dimensional (electron–electron double resonance, or ELDOR) experiments, demonstrating the full capabilities of the HiPER spectrometer. As an example, we present FT-detected version of the ELDOR-detected nuclear magnetic resonance technique for the TEMPOL radical. We also describe potential applications of the new instrumentation for studies of gadolinium qubits.

36 MATERIALS SCIENCE↗

CF + excitation in the interstellar medium

The detection of CF + in interstellar clouds potentially allows astronomers to infer the elemental fluorine abundance and the ionization fraction in ultraviolet-illuminated molecular gas. Because local thermodynamic equilibrium (LTE) conditions are hardly fulfilled in the interstellar medium (ISM), the accurate determination of the CF + abundance requires one to model its non-LTE excitation via both radiative and collisional processes. Here, we report quantum calculations of rate coefficients for the rotational excitation of CF + in collisions with para- and ortho-H 2 (for temperatures up to 150 K). As an application, we present non-LTE excitation models that reveal population inversion in physical conditions typical of ISM photodissociation regions (PDRs). We successfully applied these models to fit the CF + emission lines previously observed toward the Orion Bar and Horsehead PDRs. The radiative transfer models achieved with these new rate coefficients allow the use of CF + as a powerful probe to study molecular clouds exposed to strong stellar radiation fields.

79 ASTRONOMY AND ASTROPHYSICS↗

Carbon–Oxygen Classical Novae Are Galactic 7 Li Producers as well as Potential Supernova Ia Progenitors

Here, we report on studies of classical nova (CN) explosions where we follow the evolution of thermonuclear runaways (TNRs) on carbon–oxygen (CO) white dwarfs (WDs). We vary both the mass of the WD (from 0.6 M ⊙ to 1.35 M ⊙ ) and the composition of the accreted material. Our simulations are guided by the results of multidimensional studies of TNRs in WDs, which find that sufficient mixing with WD core material occurs after the TNR is well underway, and levels of enrichment are reached that agree with observations of CN ejecta abundances. We use NOVA (our one-dimensional hydrodynamic code) to accrete solar matter until the TNR is ongoing and then switch to a mixed composition (either 25% WD material and 75% solar or 50% WD material and 50% solar). Because the amount of accreted material is inversely proportional to the initial 12 C abundance, by first accreting solar matter the amount of material taking part in the outburst is larger than in those simulations where we assume a mixed composition from the beginning. Our results show large enrichments of 7 Be in the ejected gases, implying that CO CNe may be responsible for a significant fraction (~100 M ⊙ ) of the 7 Li in the galaxy (~1000 M ⊙ ). Although the ejected gases are enriched in WD material, the WDs in these simulations eject less material than they accrete. We predict that the WD is growing in mass as a consequence of the accretion–outburst–accretion cycle, and CO CNe may be an important channel for SN Ia progenitors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Discovery of an Intrinsic Antiferromagnetic Semiconductor EuSc 2 Te 4 With Magnetism‐Driven Nonlinear Transport

Magnetic topological materials have recently emerged as a promising platform for studying quantum geometry by the nonlinear transport in thin film devices. In this work, an antiferromagnetic (AFM) semiconductor EuSc₂Te₄ as the first bulk crystal that exhibits quantum geometry-driven nonlinear transport is reported. This material crystallizes into an orthorhombic lattice with AFM order below 5.2 K and a bandgap of less than 50 meV. The calculated band structure aligns with the angle-resolved photoemission spectroscopy spectrum. The AFM order preserves combined space-time inversion symmetry but breaks both spatial inversion and time-reversal symmetry, leading to the nonlinear Hall effect (NLHE). Nonlinear Hall voltage measured in bulk crystals appears at zero field, peaks near the spin-flop transition as the field increases, and then diminishes as the spin moments align into a ferromagnetic order. This field dependence, along with the scaling analysis of the nonlinear Hall conductivity, suggests that the NLHE of EuSc₂Te₄ involves contributions from quantum metric, in addition to extrinsic contributions, such as spin scattering and junction effects. Furthermore, this NLHE is found to have the functionality of broadband frequency mixing, indicating its potential applications in electronics. This work reveals a new avenue for studying magnetism-induced nonlinear transport in magnetic materials.

36 MATERIALS SCIENCE↗

Nonadiabatic Molecular Dynamics Study of the Relaxation Pathways of Photoexcited Cyclooctatetraene

In the current study, we present nonadiabatic (NAMD) and adiabatic molecular dynamics simulations of the transition-state dynamics of photoexcited cyclooctatetraene (COT). The equilibrium-state structure and absorption spectra are analyzed using the semiempirical Austin Model 1 potential. The NAMD simulations are obtained by a surface-hopping algorithm. We analyzed in detail an active excited to ground state relaxation pathway accompanied by an S 2 /S 3 (D 2d ) → S 1 (D 8h ) → S o (D 4h ) → S o (D 2d ) double-bond shifting mechanism. The simulated excitation lifetime is in good agreement with experiment. The first excited singlet state S1 plays a crucial role in the photochemistry. Here, the obtained critical molecular conformations, energy barrier, and transition-state lifetime results will provide a basis for further investigations of the bond-order inversion and photoswitching process of COT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Early Inference of Nuclear Technology-Directed Research Activities of Authors from Scientific Publications

Nuclear research articles can provide information about early nuclear proliferation indicators such as influential research entities and technology capability levels of a country, but detection of nuclear activities typically occurs after they have started. We investigate the extent to which nuclear research articles can be used to infer whether a research entity will acquire or develop a nuclear technology before it happens. Early detection of nuclear proliferation or technology development indicators from data is challenging due to partial observability, sparse and unlabeled information, and confounding signals from multiple concurrent activities. This paper presents the early detection problem as a sequential decision-making, goal inference problem, where the objective is to characterize and predict an individual’s, organization’s, or a country’s intent (unobserved goal-directed behavior) towards developing a nuclear capability from partially observed sequences of their research publications, using inverse reinforcement learning and Bayesian goal inference methods. A computational framework is presented, and its application demonstrated using 29,196 Scopus records for a case study related to a civil nuclear capability. The case study results serve as a proof-of-concept demonstration for inference of technology-directed research activity of authors who publish in the nuclear domain. The inference method, combined with advanced computing, may be used to assess and monitor activities pertaining to early developmental stages of a nuclear technology or capability, which in turn can help to identify and prioritize activities with nuclear proliferation potential for further investigation.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

A remote sensing algorithm for vertically resolved cloud condensation nuclei number concentrations from airborne and spaceborne lidar observations

Cloud condensation nuclei (CCN) are mediators of aerosol–cloud interactions (ACIs), contributing to the largest uncertainties in the understandings of global climate change. We present a novel remote-sensing-based algorithm that quantifies the vertically resolved CCN number concentrations (N CCN ) using aerosol optical properties measured by a multiwavelength lidar. The algorithm considers five distinct aerosol subtypes with bimodal size distributions. The inversion used the lookup tables developed in this study, based on the observations from the Aerosol Robotic Network, to efficiently retrieve optimal particle size distributions from lidar measurements. The method derives dry aerosol optical properties by implementing hygroscopic enhancement factors in lidar measurements. The retrieved optically equivalent particle size distributions and aerosol-type-dependent particle composition are utilized to calculate critical diameters using κ-Köhler theory and N CCN at six supersaturations ranging from 0.07 % to 1.0 %. Sensitivity analyses indicate that uncertainties in extinction coefficients and relative humidity greatly influence the retrieval error in N CCN . The potential of this algorithm is further evaluated by retrieving N CCN using airborne lidar from the NASA ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES) campaign and is validated against simultaneous measurements from the CCN counter. The independent validation with robust correlation demonstrates promising results. Furthermore, the N CCN has been retrieved for the first time using a proposed algorithm from spaceborne lidar – Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) – measurements. The application of this new capability demonstrates the potential for constructing a 3D CCN climatology at a global scale, which helps to better quantify ACI effects and thus reduce the uncertainty in aerosol climate forcing.

54 ENVIRONMENTAL SCIENCES↗

Comparison of Ce( iv )/Th( iv )-alkynyl complexes and observation of a trans -influence ligand series for Ce( iv )

Organometallic cerium(iv) complexes have been challenging to isolate and characterize due to the strongly oxidizing nature of the cerium(iv) cation. Herein, we report two cerium(iv) alkynyl complexes, [Ce(TriNOx)(C[triple bond, length as m-dash]C-SiMe3)] (1-CeTMS) and [Ce(TriNOx)(C[triple bond, length as m-dash]C-Ph)] (1-CePh) (TriNOx3- = tris(2-tert-butylhydroxylaminato)benzylamine), that include terminal alkyne moieties. The isostructural thorium analogue [Th(TriNOx)(C[triple bond, length as m-dash]C-SiMe3)] (1-ThTMS) was also synthesized and compared with 1-CeTMS in bond distance, 13C-NMR spectra, vibrational spectra and electronic structure. The Ce-C bond distances were 2.501(3) Å for 1-CePh and 2.513(5) Å for 1-CeTMS on the shorter end of the few reported CeIV-C single bonds (2.478(3)-2.705(2) Å), possibly indicating significant Ce 5d- and 4f-orbital involvement. 13C-NMR spectroscopy was also consistent with Ce-C covalency, with significantly deshielded resonances ranging from 185-213 ppm. Such 13C-NMR shifts demonstrate a strong influence from spin-orbit coupling (SOC) effects, corroborated by computational studies. Raman analysis showed ν C[triple bond, length as m-dash]C stretching frequencies of 2000 cm-1 (1-CeTMS) and 2052 cm-1 (1-CePh), indicating the cerium(iv)-alkynyl interaction, compared to the parent HC[triple bond, length as m-dash]CPh (IR = 2105 cm-1 and Raman = 2104 cm-1). L3-edge X-ray absorption measurements revealed a predominant Ce(iv) electronic configuration, and magnetic measurements revealed temperature-independent paramagnetism. Electrochemical studies similarly revealed the electron donating ability of the alkynyl ligands, stronger than either fluoride or imido ligands for the Ce(iv)(TriNOx)-framework, with a cerium(iv/iii) reduction potential of E pc = -1.58 to -1.66 V vs. Fc/Fc+. Evidence for a trans-influence has been observed by evaluating a series including previously reported [CeIV(TriNOx)X]+/0 complexes with axial ligands X = THF, I-, Br-, Cl-, F-, -C[triple bond, length as m-dash]C-Ph, -C[triple bond, length as m-dash]C-SiMe3, -NH(3,5-(CF3)2-Ar), -OSiPh3, -N(M(L))(3,5-(CF3)2-Ar) [M(L) = Li(TMEDA), K(DME)2 or Cs(2,2,2-crypt)]. These data stand in contrast with previous reports of an inverse trans-influence at cerium(iv) and point to differences in involvement of cerium 4f- versus 5d-orbitals in the electronic structures of the complexes.

Yang, Qiaomu↗

Coherent diffractive imaging with twisted X-rays: Principles, applications, and outlook

Recent technological breakthroughs in synchrotron and x-ray free electron laser facilities have revolutionized nanoscale structural and dynamic analyses in condensed matter systems. This review provides a comprehensive overview of the advancements in coherent scattering and diffractive imaging techniques, which are now at the forefront of exploring materials science complexities. These techniques, notably Bragg coherent diffractive imaging and x-ray photon correlation spectroscopy, x-ray magnetic dichroism, and x-ray correlation analysis leverage beam coherence to achieve volumetric three-dimensional imaging at unprecedented sub-nanometer resolutions and explore dynamic phenomena within sub-millisecond timeframes. Such capabilities are critical in understanding and developing advanced materials and technologies. Simultaneously, the emergence of chiral crystals—characterized by their unique absence of standard inversion, mirror, or other roto-inversion symmetries—presents both challenges and opportunities. These materials exhibit distinctive interactions with light, leading to phenomena such as molecular optical activity, chiral photonic waveguides, and valley-specific light emissions, which are pivotal in the burgeoning fields of photonic and spintronic devices. This review elucidates how novel x-ray probes can be leveraged to unravel these properties and their implications for future technological applications. A significant focus of this review is the exploration of new avenues in research, particularly the shift from conventional methods to more innovative approaches in studying these chiral materials. Inspired by structured optical beams, the potential of coherent scattering techniques utilizing twisted x-ray beams is examined. This promising direction not only offers higher spatial resolution but also opens the door to previously unattainable insights in materials science. By contextualizing these advancements within the broader scientific landscape and highlighting their practical applications, this review aims to chart a course for future research in this rapidly evolving field.

Physics↗

Preconditioned least‐squares Petrov–Galerkin reduced order models

Abstract In this article, we introduce a methodology for improving the accuracy and efficiency of reduced order models (ROMs) constructed using the least‐squares Petrov–Galerkin (LSPG) projection method through the introduction of preconditioning. Unlike prior related work, which focuses on preconditioning the linear systems arising within the ROM numerical solution procedure to improve linear solver performance, our approach leverages a preconditioning matrix directly within the minimization problem underlying the LSPG formulation. Applying preconditioning in this way has the potential to improve ROM accuracy for several reasons. First, preconditioning the LSPG formulation changes the norm defining the residual minimization, which can improve the residual‐based stability constant bounding the ROM solution's error. The incorporation of a preconditioner into the LSPG formulation can have the additional effect of scaling the components of the residual being minimized to make them roughly of the same magnitude, which can be beneficial when applying the LSPG method to problems with disparate scales (e.g., dimensional equations, multi‐physics problems). Importantly, we demonstrate that an “ideal preconditioned” LSPG ROM (a ROM in which the preconditioner is the inverse of the Jacobian of its corresponding full order model) emulates projection of the full order model solution increment onto the reduced basis. This quantity defines a lower bound on the error of a ROM solution for a given reduced basis. By designing preconditioners that approximate the Jacobian inverse—as is common in designing preconditioners for solving linear systems—it is possible to obtain a ROM whose error approaches this lower bound. The proposed approach is evaluated on several mechanical and thermo‐mechanical problems implemented within the Albany HPC code and run in the predictive regime, with prediction across material parameter space. We demonstrate numerically that the introduction of simple Jacobi, Gauss‐Seidel, and ILU preconditioners into the proper orthogonal decomposition/LSPG formulation reduces significantly the ROM solution error, the reduced Jacobian condition number, the number of nonlinear iterations required to reach convergence, and the wall time (thereby improving efficiency). Moreover, our numerical results reveal that the introduction of preconditioning can deliver a robust and accurate solution for test cases in which the unpreconditioned LSPG method fails to converge.

Lindsay, Payton↗