Concrete compressive strain behavior and magnitudes under uniaxial fatigue loading
Not provided.
Construction & Building Technology↗
SEARCH · Engineering Papers
Results for “nonlinear analysis”
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.
An initial exploration of Bayesian model calibration for estimating the composition of rocks and soils on Mars
The Mars Curiosity rover carries an instrument, ChemCam, designed to measure the composition of surface rocks and soil using laser-induced breakdown spectroscopy (LIBS). The measured spectra from this instrument must be analyzed to identify the component elements in the target sample, as well as their relative proportions. This process, which we call disaggregation, is complicated by so-called matrix effects, which describe nonlinear changes in the relative heights of emission lines as an unknown function of composition due to atomic interactions within the LIBS plasma. In this work, we explore the use of the plasma physics code ATOMIC, developed at Los Alamos National Laboratory, for the disaggregation task. ATOMIC has recently been used to model LIBS spectra and can robustly reproduce matrix effects from first principles. The ability of ATOMIC to predict LIBS spectra presents an exciting opportunity to perform disaggregation in a manner not yet tried in the LIBS community, namely via Bayesian model calibration. However, using it directly to solve our inverse problem is computationally intractable due to the large parameter space and the computation time required to produce a single output. Therefore, we also explore the use of emulators as a fast solution for this analysis. We discuss a proof of concept Gaussian process emulator for disaggregating two-element compounds of sodium and copper. The training and test datasets were simulated with ATOMIC using a Latin hypercube design. After testing the performance of the emulator, we successfully recover the composition of 25 test spectra with Bayesian model calibration.
97 MATHEMATICS AND COMPUTING↗
A Generative Model for Realistic Galaxy Cluster X-Ray Morphologies
Abstract The X-ray morphologies of clusters of galaxies display significant variations, reflecting their dynamical histories and the nonlinear dependence of X-ray emissivity on the density of the intracluster gas. Qualitative and quantitative assessments of X-ray morphology have long been considered a proxy for determining whether clusters are dynamically active or “relaxed.” Conversely, the use of circularly or elliptically symmetric models for cluster emission can be complicated by the variety of complex features realized in nature, spanning scales from megaparsecs down to the resolution limit of current X-ray observatories. In this work, we use mock X-ray images from simulated clusters from The Three Hundred project to define a basis set of cluster image features. We take advantage of the clusters’ approximate self-similarity to minimize the differences between images before encoding the remaining diversity through a distribution of high-order polynomial coefficients. Principal component analysis then provides an orthogonal basis for this distribution, corresponding to natural perturbations from an average model. This representation allows novel, realistically complex X-ray cluster images to be easily generated, and we provide code to do so. The approach provides a simple way to generate training data for cluster image analysis algorithms and could be straightforwardly adapted to generate clusters displaying specific types of features or selected by physical characteristics available in the original simulations.
79 ASTRONOMY AND ASTROPHYSICS↗
Tidal frequencies and quasiperiodic subsurface water level variations dominate redox dynamics in a salt marsh system
Salt marshes are hotspots of nutrient processing en route to sensitive coastal environments. Whilst our understanding of these systems has improved over the years, we still have limited knowledge of the spatiotemporal variability of critical biogeochemical drivers within salt marshes. Sea-level rise will continue to force change on salt marsh functioning, highlighting the urgency of filling this knowledge gap. Our study was conducted in a central California estuary experiencing extensive marsh drowning and relative sea-level rise, making it a model system for such an investigation. Here we instrumented three marsh positions subjected to different degrees of tidal inundation (6.7%, 8.9%, and 11.2% of the time for the upper, middle, and lower marsh positions, respectively), providing locations with varied biogeochemical characteristics and hydrological interactions at the site. We continuously monitored redox potential (Eh) at depths of 0.1, 0.3, and 0.5 m, subsurface water levels (WL), and temperature at 0.7 m depth at each marsh position. To understand how drivers of subsurface biogeochemical processes fluctuate across tidal cycles, we used wavelet analyses to explain the interactions between Eh and WL. We found that tidal forcing significantly affects key drivers of biogeochemical processes by imparting controls on Eh variability, likely driving subsurface hydro-biogeochemistry of the salt marsh. Wavelet coherence showed that the Eh-WL relationship is nonlinear, and their lead–lag relationship is variable. We found that precipitation events perturb Eh at depth over timescales of hours, even though WL shows relatively minimal change during events. Furthermore, this work highlights the importance of high frequency in situ measurements, such as Eh, to help explain factors that govern subsurface biogeochemistry and hydrological processes in salt marshes.
54 ENVIRONMENTAL SCIENCES↗
Analysis of degradation in residential battery energy storage systems for rate-based use-cases
This article examines the impact of residential battery energy storage (BES) systems' operational modes on the life (i.e. usable energy capacity) of the battery under several climatic conditions and battery chemistries. The sharp increase in residential BES installations has been a result of decreasing costs of batteries, increase in rate structures motivated applications such as solar self-consumption and time-of-use energy management, and customers purchasing these systems for backup power. While these different modes of operations provide a combination of increased bill savings, reliability of supply, and energy sustainability to the customer, their operational characteristics vary significantly between use-cases. Functional life of BES systems is known to be strongly dependent on their operating conditions. In this article, we analyze the operation of residential BES systems under different rate-based use-cases, for different battery chemistries and cell designs, and under different environmental conditions. This is conducted by simulating the control of BES operations using rate-based cycling algorithms and analyzing the prognosis of multiple battery lifetime models that consider complex nonlinear dependencies of operational stress factors such as state-of-charge, depth-of-discharge, and temperature on degradation. Significant variations in battery life are observed owing to the differences in characteristics of the uses-cases coupled with environmental conditions and battery chemistries. Such differences lead to the conclusion that choice of residential BES technologies and chemistries should account for their intended use-cases.
25 ENERGY STORAGE↗
Impact of Advection Schemes on Tracer Interrelationships in Large-Eddy Simulations of Deep Convection
This study investigates the preservation of tracer interrelationships during advection in large-eddy simulations of an idealized deep convective cloud, which is particularly relevant to chemistry, aerosol, and cloud microphysics models. Employing the Cloud Model 1, advection is represented using third-, fifth-, and seventh-order weighted essentially non-oscillatory schemes. As a simplified analogy for cloud hydrometeors and aerosols, several inert passive tracers following linear and nonlinear relationships are initialized after the cloud reaches ~6-km depth. Numerical mixing in the simulated turbulent convective clouds leads to significant deviations from the initial nonlinear relationships between tracers. In these simulations, a considerable fraction of the grid points where the tracers’ nonlinear relationships are altered from advection are classified as unrealistic (e.g., ~13% for the environmental tracers on average), including errors from range-preserving unmixing and overshooting. Errors in the sum of three tracers are also relatively large, ranging between ~1% and 16% for 5% of the grid points in and near the cloud. The magnitude of unrealistic mixing and errors in the sum of three tracers generally increase with the order of accuracy of the advection scheme. These results are consistent across model grid spacings ranging from 50 to 200 m, and across three different flow realizations for each combination of grid spacing and advection scheme tested. Tests employing a previously proposed scalar normalization procedure show substantially reduced errors in the sum of three tracers with a relatively small negative impact on other tracer relationships. This analysis, therefore, suggests efficacy of the normalization procedure when applied to turbulent three-dimensional cloud simulations.
54 ENVIRONMENTAL SCIENCES↗
Suhl instabilities in nanoscopic spheroids
We simulate the magnetization dynamics of a permalloy spheroid of nanoscopic size in zero external field, such that both dipolar and exchange interactions are important. Low excitation power is used to obtain the frequencies and mode patterns of many normal modes. At higher power, non-linear three and four mode couplings between magnons carrying orbital angular momentum are observed to give rise to Suhl instabilities. Suhl’s analysis of the selection rules governing the allowed processes is extended to initial states other than uniform precession. These rules are studied and confirmed by the simulations. Both down- and up-conversion are seen as well as three and four-mode processes. As a result, general trends are inferred for preferred instabilities among those that are allowed, although the thresholds for some instabilities appear to be very high.
3D nanomagnetism↗
Statistical relationships across epigenomes using large-scale hierarchical clustering
Recent advances in genomics and sequencing platforms have revolutionized our ability to create immense data sets, particularly for studying epigenetic regulation of gene expression. However, the avalanche of epigenomic data is difficult to parse for biological interpretation given nonlinear complex patterns and relationships. This attractive challenge in epigenomic data lends itself to machine learning for discerning infectivity and susceptibility. In this study, we explore over 3000 epigenomes of uninfected individuals and provide a framework to characterize the relationships among epigenetic modifiers, their modifiers, genetic loci, and specific immune cell types across all chromosomes using hierarchical clustering. Hierarchical clustering of epigenomic data revealed consistent epigenetic patterns across chromosomes, demonstrating that variation due to epigenetic modifiers is greater than variation between cell types. Gene Ontology and KEGG pathway analyses indicated significant enrichment of genes involved in chromatin remodeling, mRNA splicing, immune responses, and the regulation of microRNAs and snoRNAs. Epigenetic modifiers frequently formed biologically relevant clusters, including the cohesin complex, RNA Polymerase II transcription factors, and PRC2 complex members. These clustering behaviors remained consistent across all chromosomes, supported by entropy analysis and high Adjusted Rand Index scores, indicating robust cross-chromosomal similarity. Co-occurrence analysis further revealed specific sets of modifiers that consistently appeared together within clusters, reflecting shared biological functions and interactions. Validation using another dataset confirmed the reproducibility of these clustering patterns and modifier co-occurrence relationships, underscoring the reliability and generalizability of the methodology.
97 MATHEMATICS AND COMPUTING↗
Single crystal growth, chemical defects, magnetic and transport properties of antiferromagnetic topological insulators ( Ge 1 – δ – x Mn x ) 2 Bi 2 Te 5 ( x ≤ 0.47 , 0.11 ≤ δ ≤ 0.20 )
Magnetic topological insulators provide a platform for emergent phenomena arising from the interplay between magnetism and band topology. Here we report the single crystal growth, crystal structure, magnetic and transport properties, as well as the neutron scattering studies of topological insulator series (Ge 1–δ–x Mn x ) 2 Bi 2 Te 5 (x ≤ 0.47, 0.11 ≤ δ ≤ 0.20). Upon doping up to x = 0.47, the lattice parameter c decreases by 0.8%, while the lattice parameter a remains nearly unchanged. Significant Ge vacancies and Ge/Bi site mixing are revealed via elemental analysis as well as refinements of the neutron and x-ray diffraction data, resulting in holes dominating the charge transport. At x = 0.47, below 10.8 K, a bilayer A-type antiferromagnetic ordered state emerges, featuring an ordered moment of 3.0(3) μ B /Mn at 5 K, with the c axis as the easy axis. Magnetization data unveils a much stronger effective interlayer antiferromagnetic exchange interaction and a much smaller uniaxial anisotropy compared to MnBi 2 Te 4 . We attribute the former to the shorter nearest-neighbor Mn-Mn interlayer superexchange path and the latter to the smaller ligand-field splitting in (Ge 1–δ–x Mn x ) 2 Bi 2 Te 5 . Finally, our study demonstrates that this series of materials holds promise for the investigation of the layer Hall effect and quantum metric nonlinear Hall effect.
36 MATERIALS SCIENCE↗
Detecting Permafrost Active Layer Thickness Change From Nonlinear Baseflow Recession
We report permafrost underlies about one fifth of the global land area and affects ground stability, freshwater runoff, soil chemistry, and surface-atmosphere gas exchange. The depth of thawed ground overlying permafrost (active layer thickness) has broadly increased across the Arctic in recent decades, coincident with a period of increased streamflow, especially the lowest flows (baseflow). Mechanistic links between active layer thickness and baseflow have recently been explored using linear reservoir theory, but most watersheds behave as nonlinear reservoirs. We derive theoretical nonlinear relationships between long-term average saturated soil thickness $\bar{η}$(proxy for active layer thickness) and long-term average baseflow. When applied to 38 years of daily streamflow data for the Kuparuk River basin on the North Slope of Alaska, the theory predicts $\bar{η}$ increased 0.17 ± 0.22[2σ] cm a -1 between 1983 and 2020 (6.4 ± cm total). The rate of increase nearly doubled to 0.29 ± 0.31 cm a -1 between 1990 and 2020, during which time local field measurements from Circumpolar Active Layer Monitoring sites indicate the active layer increased 0.31 ± 0.22 cm a -1 . The predicted rate of increase more than doubled again between 2002 and 2020, outpacing a near doubling of observed active layer thickening, consistent with trends in terrestrial water storage inferred from Gravity Recovery and Climate Experiment satellite gravimetry and Modern-Era Retrospective Analysis for Research and Applications climate reanalysis. Overall, hydrologic change is accelerating in the Kuparuk River basin, and we provide a theoretical framework for estimating basin-scale changes in active layer water storage from streamflow measurements.
54 ENVIRONMENTAL SCIENCES↗
Performance and transport in the ARC tokamak
The ARC TM tokamak, a high-field (𝐵 𝑇 = 11.4 T) fusion power plant, under development by Commonwealth Fusion Systems, is studied using a suite of integrated modelling tools to predict its fusion power generation (𝑃𝑓𝑢𝑠), transport and confinement properties. Analysis is based off an ARC operational point scoped first with zero-dimensional (0-D) plasma operational contour (POPCON) modelling to produce 1.13 GW of fusion power. A suite of integrated modelling tools (TRANSP, ASTRA and TORAX) were applied to predict the performance and kinetic profiles of the ARC design point, yielding a range of predicted performance spanning from ∼900 to 1300 MW in rough quantitative agreement with POPCON predictions. The sensitivity of these results to uncertain modelling inputs was probed using scans of pedestal boundary conditions around EPED-predicted values (total pressure and temperature ratios), tungsten concentration and seperatrix density around their nominal assumptions. Pedestal pressure and pedestal top (𝑇 𝑖 /𝑇 𝑒 ) play a large role in 1.5-dimensional performance predictions, able to modify the predicted 𝑃 𝑓𝑢𝑠 by a factor of 2 within reasonable assumptions. High-fidelity core nonlinear gyrokinetic profile predictions, performed using CGYRO (Candy et al. 2016 J. Comput. Phys., vol. 324, pp. 73–93) coupled with the PORTALS (Rodriguez-Fernandez et al. 2024 Nucl. Fusion, vol. 64, 076034; Phys. Plasmas, vol. 31, 2024, 062501) framework, yield substantially lower performance (𝑃 𝑓𝑢𝑠 =677 MW) compared with 0-D and medium-fidelity modelling for nominal assumptions, showing that there is non-negligible uncertainty between models and that future work on SPARC may help resolve discrepancies. Lower overall performance results from significantly reduced volume-averaged densities and temperatures, along with reduced levels of density and temperature peaking. Turbulence and transport are largely dominated by ion temperature gradient across the profile, confirmed by both linear stability and the response of the nonlinear fluxes to changes in gradients, with some impact of kinetic ballooning modes in the deep core. This work represents one of the most complete scoping of potential fusion power plant conditions performed to date. The extensive integrated modelling provides confidence in ARC performance approaching 1 GW, while nonlinear gyrokinetic modelling results in open questions into the physics of density and temperature peaking in fusion-power-plant-relevant operational space. A discussion of results and the role that the SPARC tokamak (Creely et al. 2020 J. Plasma Phys., vol. 86, 865860502) will play in informing ARC design, performance and operation is presented.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
The impact of alkali‐ion intercalation on redox chemistry and mechanical deformations: Case study on intercalation of Li, Na, and K ions into FePO 4 cathode
Abstract Batteries made of charge carriers from Earth‐crust abundant materials (e.g., Na, K, and Mg) have received extensive attention as an alternative to Li‐ion batteries for grid storage. However, a lack of understanding of the behavior of these larger ions in the electrode materials hinders the development of electrode structures suitable for these large ions. In this study, we investigate the impact of alkali ions (Li, Na, and K) on the redox chemistry and mechanical deformations of iron phosphate composite cathodes by using electrochemical techniques and in situ digital image correlation. Na‐ion and Li‐ion intercalation demonstrate a nearly linear correlation between electrochemical strains and the state of charge and discharge. The strain development shows nonlinear dependance on the state of charge and discharge for K ions. Strain rate calculations show that K ion intercalation results in a progressive increase in the strain rate for all cycles. Li and Na intercalation induce nearly constant strain rates with the exception of the first discharge cycle of Na intercalation. When the same amount of ions are inserted into the electrode, the electrode shows the lowest strain generation upon Li intercalation compared to larger alkali ions. Na and K ions induce similar volumetric changes in the electrode when the state of charge and discharge is around 30%. Although the electrode experiences larger absolute strain generation at the end of the discharge cycles upon Na intercalation, strain rates were found to be greater for K ions. Potential‐dependent behaviors also demonstrate more sluggish redox reactions during K intercalation, compared to Li and Na. Our quantitative analysis suggests that the strain rate, rather than the absolute value of strain, is the critical factor in amorphization of the crystalline electrode.
25 ENERGY STORAGE↗
Performance Improvements of the Griffin Solvers in FY24
The Griffin code is a MOOSE-based reactor physics application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, we have made significant efforts to improve the performance of transport solver options and cross-section generation for the efficient use of Griffin in advanced reactor applications. For the HFEM-PN solver, the residual evaluations of HFEM kernels were optimized by utilizing the pre- computed averaged cross sections for individual elements. Numerical integration involving the evaluation of basis functions at quadrature points was bypassed by facilitating precomputed element mass matrices for response matrices. Red-black iterations were improved by introducing a new generalized minimum residual based solver. The memory usage of response matrix storage was significantly reduced by applying basis function rotations on interfaces and calculating volumetric odd-parity moments on the fly. Additionally, the adjoint flux and transient calculation capabilities of the HFEM-PN solver were successfully implemented and verified using the TWIGL benchmark problem. For the DFEM-SN solver, memory footprint and computation time were significantly reduced by not treating angular flux vectors as the MOOSE nonlinear system vectors. Specifically for IQS, scalar adjoint weighting was introduced to further eliminate angular adjoint flux storage in the MOOSE auxiliary system. It was demonstrated through the three-dimensional Advanced Burner Test Reactor core problem that the memory usage for transient calculations with the IQS method was reduced by over 7.5× compared to before the optimizations. For the self-shielding application programming interface, a new double-heterogeneity treatment method, named the Bell Function-Based Analytic Two-Region Slowing Down Method, was developed to efficiently flux-volume homogenize TRISO particles with the matrix. Additionally, optimizations were made to hyper- fine group (HFG) slowing down calculations by pretabulating collision probability coefficients and grouping isotopes, significantly reducing the computational time for calculating scattering sources per HFG. Lastly, the pin power reconstruction module was extended to account for temporal behavior in a microreactor analysis problem, specifically for a control drum transient. Verification tests for each of these improvements demonstrated significant performance enhancements and memory reduction.
22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
Analysis of DESI × DES using the Lagrangian effective theory of LSS
In this work we use Lagrangian perturbation theory to analyze the harmonic space galaxy clustering signal of the Bright Galaxy Survey (BGS) and luminous red galaxies (LRGs) targeted by the dark energy spectroscopic instrument (DESI), combined with the galaxy-galaxy lensing signal measured around these galaxies using Dark Energy Survey Year 3 source galaxies. The BGS and LRG galaxies are extremely well characterized by DESI spectroscopy and, as a result, lens galaxy redshift uncertainty and photometric systematics contribute negligibly to the error budget of our “2 × 2-point” analysis. On the modeling side, this work represents the first application of the SPINOSAURUS code, implementing an effective field theory model for galaxy intrinsic alignments, and we additionally introduce a new scheme ( MAIAR ) for marginalizing over the large uncertainties in the redshift evolution of the intrinsic alignment signal. Furthermore, this is the first application of a hybrid effective field theory model for galaxy bias based on the A emulus 𝜈 simulations. Our main result is a measurement of the amplitude of the lensing signal, 𝑆 8 = 𝜎 8 (Ω 𝑚 /0.3) 0.5 = 0.850$^{+0.042}_{−0.050}$, consistent with values of this parameter derived from the primary cosmic microwave background. This constraint is artificially improved by a factor of 51% if we assume a more standard, but restrictive parametrization for the redshift evolution and sample dependence of the intrinsic alignment signal, and 63% if we additionally assume the nonlinear alignment model. Furthermore, we show that when fixing the cosmological model to the best-fit values from Planck PR4 there is > 5𝜎 evidence for a deviation of the evolution of the intrinsic alignment signal from the functional form that is usually assumed in cosmic shear and galaxy-galaxy lensing studies.
79 ASTRONOMY AND ASTROPHYSICS↗
A scalable exponential-DG approach for nonlinear conservation laws: With application to Burger and Euler equations
In this work, we propose an Exponential DG framework for partial differential equations. We decompose 7 governing equations into linear and nonlinear parts to which we apply the discontinuous Galerkin 8 (DG) spatial discretization. In particular, we construct the linear part using Jacobian that effectively 9 capture stiff characteristics in the system. The former is integrated analytically, whereas the latter 10 is approximated. This approach i) is stable with a large Courant number (Cr > 1); ii) supports 11 high-order solutions both in time and space; iii) is computationally favorable compared to IMEX 12 DG methods with no preconditioner; iv) becomes comparable to explicit RKDG methods on uniform 13 mesh and beneficial on non-uniform grid for Euler equations; v) is scalable in a modern massively 14 parallel computing architecture due to its explicit nature of exponential time integrators and com15 pact communication stencil of DG method. Numerical results demonstrate the performance of our 16 proposed methods through various examples. We also discuss the stability and convergence analysis 17 for our exponential DG scheme in the context of Burgers equation.
42 ENGINEERING↗
In Situ Observation of Ion Migration in a Ferroelectric Ionic Conductor Rb-KTP during Thermal Annealing
Ion exchange in Rb-doped KTiOPO 4 has facilitated significant advancements in ferroelectric domain engineering, yet understanding the underlying mechanisms remains in its infancy. We perform time-of-flight secondary ion mass spectrometry analysis on multiple periodically ion-exchanged and periodically poled Rb-doped KTiOPO 4 samples under different temperatures and annealing durations. The results are compared between annealing in air, which involved ex situ annealing before periodic poling, and vacuum annealing conducted in situ after periodic poling. The Rb + diffusion profile after periodic ion exchange forms a tooth-shaped pattern. We show that in situ annealing causes a surface pinning effect on the nonpolar face, limiting Rb + migration along the polar axis at the surface. Once the pinned layer is removed through milling, the underlying Rb + diffusion is distinctively different from the surface. Additionally, the rate of Rb + diffusion during in situ annealing is linear, while the periodic domain structures remain stable after annealing. These results contribute to understanding the ionic diffusion process in a ferroelectric ionic conductor and using ion exchange to tailor the linear and nonlinear properties of KTiOPO 4 .
36 MATERIALS SCIENCE↗
Quantum kinetic modeling of KEEN waves in a warm-dense regime
We report the first fully kinetic, quantum study of kinetic electrostatic electron nonlinear (KEEN) waves, showing that quantum diffraction systematically erodes the classical trapping mechanism, narrows harmonic locking to the fundamental, and hastens post-drive decay. Electrons are evolved with a second-order Strang-split 1D1V Wigner–Poisson solver that couples conservative semi-Lagrangian WENO advection to an analytic Fourier space update for the non-local Wigner term, while ions remain classical. We focus on collisionless dynamics in a weakly coupled regime, providing a controlled baseline before collisional extensions. Short, frequency-tuned ponderomotive pulses drive KEEN formation in a uniform Maxwellian plasma; as the dimensionless quantum parameter H rises from the classical limit to values relevant to warm-dense matter, doped semiconductors, and 2D electron systems, the drive threshold increases, higher harmonics are damped, trapped electron vortices diffuse, and the subplasma electrostatic energy relaxes to a lower stationary level, as confirmed by continuous wavelet analysis. These microscopic changes carry macroscopic weight. Ignition-scale capsules now compress matter to regimes where the electron de Broglie wavelength rivals the Debye length, making classical kinetic descriptions insufficient. By extending KEEN physics into this quantum domain, our results offer a potential diagnostic of non-equilibrium electron dynamics for next-generation inertial-confinement designs and high-energy-density platforms, indicating that predictive fusion modeling may benefit from the integration of kinetic fidelity with quantum effects.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗