Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear correlation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Nonlinear viscoelastic response of silicone additively manufactured direct ink write (DIW) foams under repetitive compression

To investigate dynamic fatigue behavior of foam in military protective applications, such as helmets, additively manufactured (AM) foams were compressively strained into the plateau region using a reduced design of experiments. A simple power law was found to govern the decline in dynamic stiffness (complex modulus) as the foams underwent the purchase order requirement of 10,000 cycles of small deformation in the plateau region. This rate of decline was newly found to correlate with the degree of nonlinearity in the material’s deformation, quantified using total harmonic distortion. Materials with low nonlinearity exhibited relatively stable stiffness across cycles, while those with high nonlinearity experienced greater losses. The observed nonlinearity depended on both applied stress and strain rate. A strong linear correlation (R 2 = 0.78) was identified between second-order nonlinearity and the time-dependent stiffness response. Two lattice structures were examined: face-centered tetragonal (FCT) and simple cubic (SC). The SC material exhibited higher total harmonic distortion (5%) and lower stiffness retention than the FCT (2%). In conclusion, these results suggest that for cyclic compression applications in a wide variety of industries such as packaging, personal protective equipment, or aerospace, selecting materials with lower stress and greater structural uniformity can enhance the stability of dynamic performance.

42 ENGINEERING↗
DOE OSTI2021

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗
DOE OSTI2026

Nonlinear electrical transport near the metal–insulator transition in V 4 O 7 thin films

Nonlinear electrical transport associated with correlated electronic states has been widely investigated in transition-metal oxides near metal–insulator transitions. Here, in this study, we investigate nonlinear transport and threshold switching in sputter-deposited V 4 O 7 thin films grown on fused silica substrates. Temperature-dependent transport measurements show a metal–insulator transition near 240 K with negligible thermal hysteresis, defining the temperature scale that governs the nonlinear electrical response. Current–voltage measurements reveal reproducible bipolar threshold switching over a broad temperature range, with the threshold voltage decreasing systematically as the transition temperature is approached. Electrothermal finite-element simulations reproduce the measured switching characteristics and show that switching occurs when localized Joule heating drives a confined region of the device toward the transition temperature. The simulated peak local temperature at threshold lies close to the intrinsic metal–insulator transition temperature of V 4 O 7 , establishing that the nonlinear electrical response originates from electrothermal feedback acting on the strongly temperature-dependent conductivity near the transition. These results identify V 4 O 7 thin films as a model correlated-oxide system in which volatile threshold switching emerges from proximity to an extended metal–insulator transition.

25 ENERGY STORAGE↗
DOE OSTI2020

Spectral fluctuations in the Sachdev-Ye-Kitaev model

We present a detailed quantitative analysis of spectral correlations in the Sachdev-Ye-Kitaev (SYK) model. We find that the deviations from universal Random Matrix Theory (RMT) behavior are due to a small number of long-wavelength fluctuations (of the order of the number of Majorana fermions N ) from one realization of the ensemble to the next one. These modes can be parameterized effectively in terms of Q-Hermite orthogonal polynomials, the main contribution being due to scale fluctuations for which we give a simple analytical estimate. Our numerical results for N = 32 show that only the lowest eight polynomials are needed to eliminate the nonuniversal part of the spectral fluctuations. The covariance matrix of the coefficients of this expansion can be obtained analytically from low-order double-trace moments. We evaluate the covariance matrix of the first six moments and find that it agrees with the numerics. We also analyze the spectral correlations in terms of a nonlinear σ-model, which is derived through a Fierz transformation, and evaluate the one and two-point spectral correlation functions to two-loop order. The wide correlator is given by the sum of the universal RMT result and corrections whose lowest-order term corresponds to scale fluctuations. However, the loop expansion of the σ-model results in an ill-behaved expansion of the resolvent, and it gives universal RMT fluctuations not only for q= 4 or higher even q-body interactions, but also for the q= 2 SYK model albeit with a much smaller Thouless energy while the correct result in this case should have been Poisson statistics. In our numerical studies we analyze the number variance and spectral form factor for N = 32 and q = 4. We show that the quadratic deviation of the number variance for large energies appears as a peak for small times in the spectral form factor. After eliminating the long-wavelength fluctuations, we find quantitative agreement with RMT up to an exponentially large number of level spacings for the number variance or exponentially short times in the case of the spectral form factor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Measurements of the Five-Dimensional Phase Space Distribution of an Intense Ion Beam

No simulation of intense beam transport has accurately reproduced measurements at the level of beam halo. One potential explanation of this discrepancy is a lack of knowledge of the initial distribution of particles in six-dimensional (6D) phase space. A direct 6D measurement of an ion beam was recently performed at the Spallation Neutron Source (SNS) Beam Test Facility (BTF), revealing nonlinear transverse-longitudinal correlations in the beam core that affect downstream evolution. Unfortunately, direct 6D measurements are limited in resolution and dynamic range; here, we discuss the use of three slits and one screen to measure a 5D projection of the 6D phase space distribution, overcoming these limitations at the cost of one dimension. We examine the measured 5D distribution before and after transport through the BTF and compare to particle-in-cell simulations. We also discuss the possibility of reconstructing the 6D distribution from 5D and 4D projections.

Hoover, Austin↗
DOE OSTI2025

Magnon-mediated exciton–exciton interaction in a van der Waals antiferromagnet

Excitons are fundamental excitations that govern the optical properties of semiconductors. Interactions between excitons can lead to various emergent phases of matter and large nonlinear optical responses. In most semiconductors, excitons interact via exchange interactions or phase-space filling. Correlated materials that host excitons coupled to other degrees of freedom could offer pathways for controlling these interactions. Here, in this study, we demonstrate magnon-mediated interactions between excitons in CrSBr, an antiferromagnetic semiconductor. These interactions manifest as the dependence of the exciton energy on the exciton density via a magnonic adjustment of the spin canting angle. Our study demonstrates the emergence of quasiparticle-mediated interactions in correlated quantum materials, leading to large nonlinear optical responses and potential device concepts such as magnon-mediated quantum transducers.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Joint analysis of Dark Energy Survey Year 3 data and CMB lensing from SPT and P l a n c k . II. Cross-correlation measurements and cosmological constraints

Cross-correlations of galaxy positions and galaxy shears with maps of gravitational lensing of the cosmic microwave background (CMB) are sensitive to the distribution of large-scale structure in the Universe. Such cross-correlations are also expected to be immune to some of the systematic effects that complicate correlation measurements internal to galaxy surveys. We present measurements and modeling of the cross-correlations between galaxy positions and galaxy lensing measured in the first three years of data from the Dark Energy Survey with CMB lensing maps derived from a combination of data from the 2500 deg 2 SPT-SZ survey conducted with the South Pole Telescope and full-sky data from the Planck satellite. The CMB lensing maps used in this analysis have been constructed in a way that minimizes biases from the thermal Sunyaev Zel’dovich effect, making them well suited for cross-correlation studies. The total signal-to-noise of the cross-correlation measurements is 23.9 (25.7) when using a choice of angular scales optimized for a linear (nonlinear) galaxy bias model. Here we use the cross-correlation measurements to obtain constraints on cosmological parameters. For our fiducial galaxy sample, which consist of four bins of magnitude-selected galaxies, we find constraints of Ω m =0.272$^{+0.032}_{-0.052}$ and S 8 ≡σ8$\sqrt{Ω_m/0.3}$=0.736$^{+0.032}_{-0.028}$ (Ω m =0.245$^{+0.026}_{-0.044}$ and S 8 =0.734$^{+0.035}_{-0.028}$) when assuming linear (nonlinear) galaxy bias in our modeling. Considering only the cross-correlation of galaxy shear with CMB lensing, we find Ω m =0.270$^{+0.043}_{-0.061}$ and S 8 =0.740$^{+0.034}_{-0.029}$. Our constraints on S8 are consistent with recent cosmic shear measurements, but lower than the values preferred by primary CMB measurements from Planck.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Higher-order symmetry plane correlations in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The correlations between event-by-event fluctuations of symmetry planes are measured in Pb-Pb collisions at a center-of-mass energy per nucleon pair $\sqrt{s_{NN}}$ = 5.02 TeV recorded by the ALICE detector at the Large Hadron Collider. This analysis is conducted using the Gaussian estimator technique, which is insensitive to biases from correlations between different flow amplitudes. The study presents, for the first time, the centrality dependence of correlations involving up to five different symmetry planes. The correlation strength varies depending on the harmonic order of the symmetry plane and the collision centrality. Comparisons with measurements from lower energies indicate no significant differences within uncertainties. Additionally, the results are compared with hydrodynamic model calculations. Although the model predictions provide a qualitative explanation of the experimental results, they overestimate the data for some observables. This is particularly true for correlators that are sensitive to the nonlinear response of the medium to initial-state anisotropies in the collision system. As these new correlators provide unique information—independent of flow amplitudes—their usage in future model developments can further constrain the properties of the strongly interacting matter created in ultrarelativistic heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A modal wave-packet model for the multi-mode Richtmyer–Meshkov instability

In this work, a model for multimode perturbations subject to the Richtmyer–Meshkov (RM) instability is presented and compared with simulations and experiments for conditions relevant to inertial confinement fusion. The model utilizes the single mode response to the RM impulse whereby its amplitude h(k, t) first grows with an initial velocity V 0 ∝ kh(k, 0) that eventually decays in time as 1/kV 0 t. Both the growth and saturation stages are subject to nonlinearities since they depend explicitly on the initial amplitude. However, rather than using the individual mode amplitude h(k, t), nonlinearity is taken to occur when the root-mean-square amplitude h rms (k, t) of a wave-packet within wavenumbers k ± δk becomes comparable to 1/k. This is done because nearby sidebands can act in unison for an auto-correlation distance 1/δk beyond nonlinearity as observed in the beam-plasma instability. Thus, the nonlinear saturation amplitude for each mode is reduced from the usual 1/k by a phase space factor that depends on the physical dimensionality, as in the Haan model for the Rayleigh–Taylor instability. In addition, for RM, the average value of kh rms for the initial spectrum is used to calculate a nonlinear factor F NL that reduces V 0 , as observed for single modes. For broadband perturbations, the model describes self-similar growth ∝t θ as successively longer wavelength modes reach saturation. The growing and saturated modes must be discerned because only the former promote θ and are enhanced by reshock and spherical convergence. All of these flows are described here by the model in good agreement with simulations and experiments.

42 ENGINEERING↗
DOE OSTI2025

Scaling up the transcorrelated density matrix renormalization group

Explicitly correlated methods, such as the transcorrelated method which shifts a Jastrow or Gutzwiller correlator from the wave function to the Hamiltonian, are designed for high-accuracy calculations of electronic structures, but their application to larger systems has been hampered by the computational cost. We develop improved techniques for the transcorrelated density-matrix renormalization group (DMRG), in which the ground state of the transcorrelated Hamiltonian is represented as a matrix product state (MPS), and demonstrate large-scale calculations of the ground-state energy of the two-dimensional Fermi-Hubbard model. Our developments stem from three technical inventions: (i) constructing matrix product operators (MPOs) of transcorrelated Hamiltonians with low bond dimension and high sparsity, (ii) exploiting the entanglement structure of the ground states to increase the accuracy of the MPS representation, and (iii) optimizing the nonlinear parameter of the Gutzwiller correlator to mitigate the nonvariational nature of the transcorrelated method. Here, we examine systems of size up to 12×12 lattice sites, four times larger than previous transcorrelated DMRG studies, and demonstrate that transcorrelated DMRG yields significant improvements over standard nontranscorrelated DMRG for equivalent computational effort. Transcorrelated DMRG reduces the error of the ground-state energy by 2.4×–14×, with the smallest improvement seen for a small system at half filling and the largest improvement in a dilute closed-shell system.

Density matrix renormalization group↗
DOE OSTI2022

Accurate parameterization of the kinetic energy functional for calculations using exact-exchange

In this report, electronic structure calculations based on Kohn-Sham density functional theory (KSDFT) that incorporate exact exchange or hybrid functionals are associated with large computational expense, a consequence of the inherent cubic scaling bottleneck and large associated prefactor, which limits the length and time scales that can be accessed. Though orbital-free density functional theory (OFDFT) calculations scale linearly with system size and are associated with significantly smaller prefactor, they are limited by the absence of accurate density-dependent kinetic energy functionals. Therefore, the development of accurate density-dependent kinetic energy functionals is important for OFDFT calculations of large realistic systems. To this end, we propose a method to train kinetic energy functional models at the exact-exchange level of theory by using a dictionary of physically relevant terms that have been proposed in the literature in conjunction with linear or nonlinear regression methods to obtain the fitting coefficients. For our dictionary, we use gradient expansion of the kinetic energy, nonlocal models proposed in the literature and their nonlinear combinations, such as a model that incorporates spatial correlations between higher order derivatives of electron density at two points. The predictive capabilities of these models are assessed by using a variety of model one-dimensional systems that exhibit diverse bonding characteristics, such as a chain of eight hydrogens, LiF, LiH, C 4 H 2 , C 4 N 2 and C 3 O 2 . We show that by using data from model one-dimensional KSDFT calculations performed using the exact-exchange functional for only a few neutral structures, it is possible to generate models with high accuracy for charged systems as well as electron and kinetic energy densities during self-consistent field iterations. In addition, we show that it possible to learn both the orbital dependent terms, i.e., the kinetic energy and the exact-exchange energy, and models that incorporate additional nonlinearities in spatial correlations, such as a quadratic model, are needed to capture subtle features of the kinetic energy density that are present in exact-exchange-based KSDFT calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Continental-Scale Geographic Trends in Barometric-Pumping Efficiency Potential: A North American Case Study

Barometric pumping is a gas transport mechanism that has important implications for many applications involving subsurface gas seepage processes. Here, this study provides the first continental-scale analysis of barometric-pumping efficiency potential based on meteorology. We quantified the barometric-pumping efficiency potential at 1,257 locations across the continental US and Canada. The results provide continental-scale geographic dependencies of barometric-pumping efficiency potential, indicating a significant correlation with latitude and a nonlinear dependence on longitude. The analysis also indicates that variability in barometric-pumping efficiency potential decreases with distance from the coast and as elevation increases. Locations far from the coastline are more likely to have upper mid-range potentials, while higher elevation locations are more likely to have low potentials. The highest barometric-pumping efficiency potentials are mostly found around the Gulf of St. Lawrence around 50°N. Locations along the Atlantic coast exhibit large-scale variations in potentials with a clear increasing trend with latitude.

58 GEOSCIENCES↗
DOE OSTI2026

On how avalanches penetrate the SOL and broaden heat loads

Recent experiments reported a correlation between power law core temperature spectra and Dα emission, suggesting that heat avalanches penetrate the SOL. This paper derives a threshold criterion for avalanche penetration using a reduced model. Avalanches with ($\bigtriangledown$$\tilde{T}$) rms > $\bigtriangledown$$\tilde{T}$ crit at the separatrix are predicted to penetrate, and so broaden the SOL and heat load distribution. $\bigtriangledown$$\tilde{T}$ crit is ~ 1/$τ$ ∥ , where $τ$ ∥ is the parallel heat flow time through the SOL. Penetration occurs when avalanches are strong enough to steepen sufficiently to shock at the separatrix. A positive correlation is found between the nonlinear drive for steepening and the penetration depth. In particular, penetration depth exceeds that of the heuristic drift limit when shocks form. Implications for numerical and physical experiments are also discussed.

SOL heat load↗
DOE OSTI2022

Worldwide Physics-Based Lifetime Prediction of c-Si Modules Due to Solder-Bond Failure

Lifetime prediction of the fielded c-Si solar modules due to location-specific weather conditions has been an important topic of photovoltaic research and the economic viability of solar energy. Data analytic techniques such as the performance ratio method, Statistical clear sky model, and Suns-Vmp methods quantify the degradation from measured data of a solar farm, however, the nonlinear time-dependence and correlated degradations make it difficult to use the empirical degradation rates for ultimate lifetime projection. In this article, we propose a complementary physics-based model to predict the solder bond failure caused by mechanical stress associated with the variations of the temperature. Integrating the worldwide weather information from NASA/NSRDB databases, the model predicts the location-specific output-power degradation and the lifetime of a module due to solder bond failure. The model parameters are calibrated against qualification tests involving thermal cycling of specific batches of modules from a specific technology/manufacturer. The results may be summarized as: 1) Modules installed at higher latitudes show a longer lifetime due to reduced damage accumulation. 2) The reduction of temperature fluctuation close to large bodies of water, such as seashores, increases solder bond lifetime significantly. 3) Relatively speaking, modules installed close to the Tropic of Cancer/Capricorn (23.5 degrees North/South) suffer from a higher solder bond damage and have a shorter lifetime, suggesting a conservative design. This model should serve as a building block of a comprehensive reliability framework that can predict the lifetime of a module that experiences simultaneous and correlated degradation mechanisms involving yellowing, corrosion, and potential-induced degradation.

14 SOLAR ENERGY↗
DOE OSTI2024

An Empirical Model For Intrinsic Alignments: Insights From Cosmological Simulations

We extend current models of the halo occupation distribution (HOD) to include a flexible, empirical framework for the forward modeling of the intrinsic alignment (IA) of galaxies. A primary goal of this work is to produce mock galaxy catalogs for the purpose of validating existing models and methods for the mitigation of IA in weak lensing measurements. This technique can also be used to produce new, simulation-based predictions for IA and galaxy clustering. Our model is probabilistically formulated, and rests upon the assumption that the orientations of galaxies exhibit a correlation with their host dark matter (sub)halo orientation or with their position within the halo. We examine the necessary components and phenomenology of such a model by considering the alignments between (sub)halos in a cosmological dark matter only simulation. We then validate this model for a realistic galaxy population in a set of simulations in the Illustris-TNG suite. We create an HOD mock with Illustris-like correlations using our method, constraining the associated IA model parameters, with the $\mathcal{X}$$^{2}_{dof}$ between our model’s correlations and those of Illustris matching as closely as 1.4 and 1.1 for orientation–position and orientation–orientation correlation functions, respectively. By modeling the misalignment between galaxies and their host halo, we show that the 3-dimensional two-point position and orientation correlation functions of simulated (sub)halos and galaxies can be accurately reproduced from quasi-linear scales down to 0.1 $\mathcal{h}$ –1 Mpc. We also find evidence for environmental influence on IA within a halo. Our publicly-available software provides a key component enabling efficient determination of Bayesian posteriors on IA model parameters using observational measurements of galaxy-orientation correlation functions in the highly nonlinear regime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Advanced Geomechanical Model to Predict the Impact of CO2-Induced Microstructural Alterations on the Cohesive-Frictional Behavior of Mt. Simon Sandstone

We investigated the influence of CO2-induced geochemical reactions on the cohesive-frictional properties of host rock within the context of CO2 storage in a saline aquifer and focused on the Mt. Simon sandstone. The research objective was to model geo-mechanical changes due to host rock exposure to CO2-saturated brine while accounting for heterogeneity, double-scale porosity, and granular structure. We formulated a three-level multi-scale model for host rocks. We conducted scanning electron microscopy analyses to probe the microstructure and grid nanoindentation to measure the mechanical response. We derived new nonlinear strength upscaling solutions to correlate the effective strength characteristics and the macroscopic yield surface to the micro-structure at the nano-, micro-, and meso-scales. Specifically, our theoretical model links CO2-induced microstructural alterations to a reduction in the size of the yield surface, and a drop in the value of the friction coefficient. In turn, regarding the Illinois Basin Decatur Project, the CO2-induced drop in friction coefficient is linked to an increase in the risk of fault slip and a higher probability of induced microseismicity during and after the end of CO2 underground injection operations. The theoretical model presented is essential for the geo-mechanical modeling of CO2 underground injection operations at multiple length-scales.

58 GEOSCIENCES↗
DOE OSTI2022

Visualization and quantum control of light-accelerated condensates by terahertz multi-dimensional coherent spectroscopy

Abstract Characterizing and controlling high-order correlation of quantum systems is key for developing quantum devices and switching technologies. Although conventional static and ultrafast spectroscopy gives access to collective excitations characterizing quantum states, more exotic correlations cannot be easily separated from other contributions. Here we develop density matrix simulations to show that seventh-order-wave-mixing peaks with distinct temperature and field dependences in two-dimensional terahertz nonlinear spectra reveal light-induced correlations in non-equilibrium superconducting states. Above critical terahertz driving, these emerging peaks split from conventional peaks along the second axis introduced by pump-probe relative phase in two-dimensional frequency space. They are photo-generated by correlations between two-photon fluctuations and interacting quasi-particle and quasi-particle/Higgs superconductor excitations. By photo-inducing persistent symmetry breaking via light-wave propagation, we also demonstrate seventh-order-wave-mixing sensing of Higgs collective modes. Our theory suggests to use multi-dimensional spectroscopy for quantum sensing of light-driven superconductivity and paves a path for quantum operations by few-cycle-THz-periodic photocurrent modulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Encoding the complete electric field of an ultraviolet ultrashort laser pulse in a near-infrared nonlinear-optical signal

We introduce a variation on the cross-correlation frequency-resolved optical gating (XFROG) technique that uses a near-infrared (NIR) nonlinear-optical signal to characterize pulses in the ultraviolet (UV). Using a transient-grating XFROG beam geometry, we create a grating using two copies of the unknown UV pulse and diffract a NIR reference pulse from it. We show that, by varying the delay between the UV pulses creating the grating, the UV pulse intensity-and-phase information can be encoded into a NIR signal. We also implemented a modified generalized-projections phase-retrieval algorithm for retrieving the UV pulses from these spectrograms. We performed proof-of-principle measurements of chirped pulses and double pulses, all at 400 nm. This approach should be extendable deeper into the UV and potentially even into the extreme UV or x-ray range.

47 OTHER INSTRUMENTATION↗