Search NASA⌕ Search

SEARCH · Search NASA

Results for “Non linear”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Integral Field Probe for Mapping of Curved Magnets

The Single Stretched Wire (SSW) method allows highly precise integral field measurements by recording voltage across a tensioned wire mounted to 2-axis linear stages at either end of the magnet aperture. However, traditional SSW probes are not well suited for curved accelerator magnets, which are essential for steering charged particles along arced trajectories in storage rings or beamlines. The tension required to eliminate sag demands a purely straight path, making them incompatible with non-linear magnet geometries. To address this limitation for curved magnets, a modified approach was developed using a segmented, 3D-printed support structure that incorporates a pre-shaped anti-sag curve. Under its own weight and that of the wire bundle, the structure deforms to lie flat while conforming to the curvature of the magnet in the horizontal plane. The optimal geometry of the probe was derived using an iterative process combining FEA simulations in Ansys Mechanical with testing of various carbon fiber-reinforced filaments. The printed and assembled probe was successfully used to measure the SDD-059 magnet at Fermilab, yielding promising results.

Jakopin, Alexander [Northern Illinois U.]↗
DOE OSTI2025

Scrambling Dynamics with Imperfections in a Solvable Model

We study how probes of quantum scrambling dynamics respond to two kinds of imperfections -- unequal forward and backward evolutions and decoherence -- in a solvable Brownian circuit model. We calculate a ``renormalized'' out-of-time-order correlator (ROTOC) in the model with $N$ qubits, and we show that the circuit-averaged ROTOC is controlled by an effective probability distribution in operator weight space which obeys a system of $N$ non-linear equations of motion. These equations can be easily solved numerically for large system sizes which are beyond the reach of exact methods. Moreover, for an operator initially concentrated on weight one $w_0=1$, we provide an exact solution to the equations in the thermodynamic limit of many qubits that is valid for all times, all non-vanishing perturbation strengths $p\gtrsim 1/\sqrt{N}$, and all decoherence strengths. We also show that a generic initial condition $w_0 >1$ leads to a metastable state that eventually collapses to the $w_0=1$ case after a lifetime $\sim \log(N/w_0)$. Our results highlight situations where it is still possible to extract the unperturbed chaos exponent even in the presence of imperfections, and we comment on the applications of our results to existing experiments with nuclear spins and to future scrambling experiments.

FOS: Physical sciences↗
DOE OSTI2026

Evolution of the Antarctic Ice Sheet from 2000–2300 and beyond: model sensitivity and uncertainty analysis using MPAS-Albany Land Ice

We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.

58 GEOSCIENCES↗
DOE OSTI2026

Unpolarized GPDs at small x and non-zero skewness

We study the small-x asymptotics of unpolarized generalized parton distributions (GPDs) and generalized transverse momentum distributions (GTMDs). Unlike the previous works in the literature, we consider the case of non-zero (but small) skewness while allowing for non-linear contributions to the evolution equations. We first show that unpolarized GPDs and GTMDs at small x are related to the eikonal dipole amplitude N, whose small-x evolution is given by the BK/JIMWLK evolution equations, and to the odderon amplitude $\mathscr{O}$, whose evolution is also known in the literature. We then show that the effect of non-zero skewness ξ ≠ 0 is to modify the value of the evolution parameter (rapidity) in the arguments for the dipole amplitudes N and $\mathscr{O}$ from Y = ln(1/x) to Y = ln min{1/|x|, 1/|ξ|}.

Generalized parton distributions (GPDs)↗
DOE OSTI2026

Complex Electrical Conductivity of a Single‐Fractured Rock: Fracture‐and‐Matrix Coupling Mechanism and Aperture Size Predictions

Fractured rocks play a crucial role in myriad natural and engineered systems, including Earth's critical zone, oil/gas/geothermal reservoirs, and geological CO 2 /H 2 /waste storage systems. While complex electrical conductivity is extensively used to estimate the pore and grain sizes of conventional porous rocks and soils, it is rarely used to predict the aperture size of fractured rocks and this remains poorly understood. Here, integrating theory, simulations, and experiments, we show that under external fields, fractured rocks follow the fracture-and-matrix coupling to make the bulk complex conductivity non-linear with respect to water conductivity. We find that the relaxation time and quadrature conductivity for porous media do not apply to fractured rocks, but, instead, reasonably accurate predictions of aperture size can be made based on the true formation factor. This study unravels the fundamental mechanism governing conduction and polarization of fractured rocks and paves the way for the non-invasive investigation of global fractured rocks.

58 GEOSCIENCES↗
DOE OSTI2024

Emergency department visits in California associated with wildfire PM 2.5 : differing risk across individuals and communities

The threats to human health from wildfires and wildfire smoke (WFS) in the United States (US) are increasing due to continued climate change. A growing body of literature has documented important adverse health effects of WFS exposure, but there is insufficient evidence regarding how risk related to WFS exposure varies across individual or community level characteristics. To address this evidence gap, we utilized a large nationwide database of healthcare utilization claims for emergency department (ED) visits in California across multiple wildfire seasons (May through November, 2012–2019) and quantified the health impacts of fine particulate matter <2.5 μm (PM 2.5 ) air pollution attributable to WFS, overall and among subgroups of the population. We aggregated daily counts of ED visits to the level of the Zip Code Tabulation Area (ZCTA) and used a time-stratified case-crossover design and distributed lag non-linear models to estimate the association between WFS and relative risk of ED visits. We further assessed how the association with WFS varied across subgroups defined by age, race, social vulnerability, and residential air conditioning (AC) prevalence. Over a 7 day period, PM 2.5 from WFS was associated with elevated risk of ED visits for all causes (1.04% (0.32%, 1.71%)), non-accidental causes (2.93% (2.16%, 3.70%)), and respiratory disease (15.17% (12.86%, 17.52%)), but not with ED visits for cardiovascular diseases (1.06% (–1.88%, 4.08%)). Analysis across subgroups revealed potential differences in susceptibility by age, race, and AC prevalence, but not across subgroups defined by ZCTA-level Social Vulnerability Index scores. These results suggest that PM 2.5 from WFS is associated with higher rates of all cause, non-accidental, and respiratory ED visits with important heterogeneity across certain subgroups. Notably, lower availability of residential AC was associated with higher health risks related to wildfire activity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Aerosol Influences on Cloud Water: Insights From ARM EPCAPE Observations With Explainable Machine Learning

This study employs an explainable machine learning (ML) framework (XGBoost-SHapley Additive exPlanations analysis) to investigate controlling factors on cloud liquid water path (LWP) using EPCAPE observations near the California coast. Aerosols are found to be the dominant factor explaining LWP variability, surpassing meteorological factors (MFs). By isolating aerosol effects from meteorological influences, the ML reveals a negative linear relationship between LWP and cloud droplet number concentration (Nd) in log space, likely driven by entrainment drying via evaporation-entrainment feedback. This aligns with the negative regime of the inverted-V relationship reported in previous studies, while no positive LWP responses are found due to a limited number of precipitating cases in EPCAPE. Furthermore, the sensitivity of LWP to Nd shows a non-linear dependence on MFs like moisture contrast between surface and free troposphere and lower-tropospheric stability. This occurs due to the interplay between the MFs' direct effects on entrainment drying and indirect effects through LWP adjustments.

EPCAPE observations↗
DOE OSTI2024

A High-Performance Discrete-Element Framework for Simulating Flow and Jamming of Moisture Bearing Biomass Feedstocks

We developed and verified a high-performance open-source discrete element method (DEM) solver with simultaneously-supported feedstock-specific interaction models, including bonded-sphere, liquid bridge, cohesion, and non-linear contact models. Our solver uses parallel data structures on hybrid central and graphics processing unit (CPU/GPU) architectures, with favorable strong scaling performance observed for large problem sizes comprised of (100 M particles), and 4X single-node GPU speedup. The particles for corn stover feedstock were conceptualized and calibrated based on experimental measurements and results. Sensitivity analyses demonstrate that the mass flow rate from a wedge hopper is governed primarily by moisture content, friction coefficient, and cohesion energy density. The model is used to reproduce experimentally observed hopper jamming results, highlighting that the experimental no-flow trends can only be achieved by using non-spherical particles, liquid bridge and cohesion models, highlighting the importance of using concurrent feedstock specialized models for the effective representation of biomass material handling problems.

bioenergy↗
DOE OSTI2026

Adaptive Control for Load-Following of Boiling Water Reactors Part I: Linear Systems and Fully-Observable Dynamics

Automation control is a key strategy to improve the economic competitiveness of nuclear power plants. Not only does it help reduce operational costs, but it also extends the value proposition of these plants to nontraditional markets, including unattended operations in remote villages and space. However, the dynamics of the operating environments of nuclear reactors are subject to changes over time, and there are no widely adopted methods to ensure that the automation strategy will remain effective over the extended durations required for these applications. Adaptive control is a discipline that offers the possibility to accommodate such changes online. However, it relies on mathematical assumptions that must be respected to ensure robustness and reliability. In this work, we derive an adaptive control formulation for linear systems in which all states are observable and apply it to an instance of load-follow operation for Boiling Water Reactors. We assumed uncertainty in two factors: the temperature coefficient and the control rod worth, both of which are affected over time by the evolution of the nuclear reactor core environment. With an arbitrary penalty factor of 5, we found that the mean absolute and integral time absolute errors can be reduced by more than 90%, underscoring the strength of adaptive control. To extend the application to more challenges, different uncertainties and load-follow trajectories, as well as new formulations that include non-linearity and partial observability, are currently being developed.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

A machine learning framework for accurate and robust analysis of radiation detector pulses

The microscopic properties of atomic nuclei are used to study various scientific questions. They are essential for understanding the fundamental forces of nature and the chemical evolution of the universe. Detecting decay radiation from radioactive nuclei makes it possible to probe these fundamental nuclear properties. Detector waveform traces may contain additional information about the radiation. Generally, advanced signal processing techniques are needed to extract this additional information, often involving fitting the waveform with model response functions using non-linear least-squares optimization with second-order gradient methods. While this is a powerful technique, it is also computationally expensive, leading to slow processing time, which scales with the volume of data. To address this problem, we have developed a machine learning (ML) approach that infers the characteristics of traces from a model detector response function. In particular, we are interested in classifying whether a single recorded trace consists of one or two pulse constituents and estimating the pulse parameters. Furthermore, our proposed ML method can precisely extract the pulses’ parameters, such as energy and timing information, and accurately classify the pulse multiplicity of a trace. Unlike non-learning-based approaches, our ML approach uses neural networks that are significantly faster at inference, as they do not require any optimization during this stage.

Curve fitting↗
DOE OSTI2024

Considerations for Defining G-Values for Aluminum-Clad Spent Nuclear Fuel

Sealed-canister dry storage of aluminum-clad spent nuclear fuel (ASNF) generated by research reactors is an alternative to current storage and disposition pathways as directed by the U.S. Department of Energy. The major challenge faced for this storage approach is radiolytic H 2 generation, including from the aluminum (oxy)hydroxide layers on the surface of ASNF. Experimental and modeling activities have been carried out to characterize the radiolytic yield as part of a DOE-sponsored research program to develop the technical basis for ASNF dry storage. The G-value is a commonly way to report results of radiolysis testing and is defined as the radiolytic yield of a species (e.g. molecular hydrogen) per unit radiation energy deposited into the material system. An independent technical review of the ASNF dry storage technical basis performed by Pacific Northwest National Laboratory raised questions about differences in G-value definitions used for experiments on ASNF surrogates consisting of aluminum samples with adherent (oxy)hydroxides compared to G-values reported in prior literature and how the magnitudes compared between different studies. Material systems resembling ASNF pose complications for measuring/defining G-values to predict the evolution of H 2 in a sealed canister, including i) accounting for radiolytic yields potentially arising from multiple sources, i.e., residual free (vapor), physisorbed, and chemisorbed/chemically bound waters; ii) deciding what portions of the multi-material system to include in the absorbed energy (radiation dose) calculation, considering possible energy exchange between materials as well as measurement limitations, and iii) capturing variations in G-value associated with non-linear yield vs. dose curves and/or dependence on the cover gas. This report summarizes previous literature information on radiolytic H 2 generation and associated G-values from mixed-material systems (generally oxides in contact with water or organic compounds) and from (oxy)hydroxides/hydrates to compare with the definitions and values for ASNF surrogate samples containing adherent aluminum (oxy)hydroxides.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

The Artificial Scientist: in-Transit Machine Learning of Plasma Simulations

Large-scale simulations or scientific experiments produce petabytes of data per run. This poses massive challenges for I/O and storage when scientific analysis workflows are run manually offline. Unsupervised deep learning-based techniques to extract patterns and non-linear relations from these large amounts of data provide a way to build scientific understanding from raw data, reducing the need for manual pre-selection of analysis steps, but require exascale compute and memory to process the full dataset available. In this paper, we demonstrate a heterogeneous streaming workflow in which plasma simulation data is streamed directly to a Machine Learning (ML) application training a model on the simulation data in-transit, completely circumventing the capacity-constrained filesystem bottleneck. This workflow employs openPMD to provide a high level interface to describe scientific data and also uses ADIOS2, to transfer volumes of data that exceed the capabilities of the filesystem. We employ experience replay to avoid catastrophic forgetting in learning from this non-steady state process in a continual manner and adapt it to improve model convergence while learning in-transit. As a proof-of-concept, we approach the ill-posed inverse problem of predicting particle dynamics from radiation in a particle-incell (PIConGPU) simulation of the Kelvin-Helmholtz instability (KHI). We detail hardware-software co-design challenges as we scale PIConGPU to full Frontier, the Top-1 system as of June 2024 Top500 list.

Kelling, Jeffrey [Helmholtz-Zentrum Dresden Rossen↗
DOE OSTI2025

Stochastic frequency fluctuation super-resolution imaging

The inherent non-linearity of intensity correlation functions can be used to spatially distinguish identical emitters beyond the diffraction limit, as achieved, for example, in super-resolution optical fluctuation imaging (SOFI). Here, we propose a complementary concept based on spectral correlation functions, termed spectral fluctuation super-resolution (SFSR) imaging. Through theoretical and computational analysis, we show that spatially resolving time-frequency correlation functions in the image plane can improve the imaging resolution by a factor of $\sqrt2$ in most cases and up to twofold for strictly two emitters. This improvement is achieved by quantifying the degree of correlation in spectral fluctuations across the spatial domain. Experimentally, SFSR can be implemented using a combination of interferometry and photon-correlation measurements. The method works for non-blinking emitters and stochastic spectral fluctuations with arbitrary temporal statistics. This suggests its utility in super-resolution microscopy of quantum emitters at low temperatures, where spectral diffusion is often more pronounced than emitter blinking.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

A Lake Biogeochemistry Model for Global Methane Emissions: Model Development, Site‐Level Validation, and Global Applicability

Abstract Lakes are important sentinels of climate change and may contribute over 30% of natural methane (CH 4 ) emissions; however, no earth system model (ESM) has represented lake CH 4 dynamics. To fill this gap, we refined a process‐based lake biogeochemical model to simulate global lake CH 4 emissions, including representation of lake bathymetry, oxic methane production (OMP), the effect of water level on ebullition, new non‐linear CH 4 oxidation kinetics, and the coupling of sediment carbon pools with in‐lake primary production and terrigenous carbon loadings. We compiled a lake CH 4 data set for model validation. The model shows promising performance in capturing the seasonal and inter‐annual variabilities of CH 4 emissions at 10 representative lakes for different lake types and the variations in mean annual CH 4 emissions among 106 lakes across the globe. The model reproduces the variations of the observed surface CH 4 diffusion and ebullition along the gradients of lake latitude, depth, and surface area. The results suggest that OMP could play an important role in surface CH 4 diffusion, and its relative importance is higher in less productive and/or deeper lakes. The model performance is improved for capturing CH 4 outgassing events in non‐floodplain lakes and the seasonal variability of CH 4 ebullition in floodplain lakes by representing the effect of water level on ebullition. The model can be integrated into ESMs to constrain global lake CH 4 emissions and climate‐CH 4 feedback.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Feature learning and generalization in deep networks with orthogonal weights

Fully-connected deep neural networks with weights initialized from independent Gaussian distributions can be tuned to criticality, which prevents the exponential growth or decay of signals propagating through the network. However, such networks still exhibit fluctuations that grow linearly with the depth of the network, which may impair the training of networks with width comparable to depth. We show analytically that rectangular networks with tanh activations and weights initialized from the ensemble of orthogonal matrices have corresponding preactivation fluctuations which are independent of depth, to leading order in inverse width. Moreover, we demonstrate numerically that, at initialization, all correlators involving the neural tangent kernel (NTK) and its descendants at leading order in inverse width—which govern the evolution of observables during training—saturate at a depth of ~20, rather than growing without bound as in the case of Gaussian initializations. We speculate that this structure preserves finite-width feature learning while reducing overall noise, thus improving both generalization and training speed in deep networks with depth comparable to width. We provide some experimental justification by relating empirical measurements of the NTK to the superior performance of deep non-linear orthogonal networks trained under full-batch gradient descent on the MNIST and CIFAR-10 classification tasks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

First constraints from marked angular power spectra with Subaru Hyper Suprime-Cam Survey First-Year Data

We present the first application of marked power spectra to weak lensing data, using maps from the Subaru Hyper Suprime-Cam Year 1 (HSC-Y1) survey. Marked convergence fields, constructed by weighting the convergence field with non-linear functions of its smoothed version, are designed to encode higher-order information while remaining computationally tractable. Using simulations tailored to the HSC-Y1 data, we test three mark functions that up- or down-weight different density environments. Our results show that combining multiple types of marked auto and cross-spectra improves constraints on the clustering amplitude parameter S8≡σ8Ωm/0.3 by ≈43 per cent compared to standard two-point power spectra. When applied to the HSC-Y1 data, this translates into a constraint on S8=0.807±0.024⁠. We assess the sensitivity of the marked power spectra to systematics, including baryonic effects, intrinsic alignment, photometric redshifts, and multiplicative shear bias. We note that some of the additional information introduced by the marked field originates from scales smaller than the scale cut, and is partly Gaussian in nature. This does not invalidate our systematic tests. These results demonstrate the promise of marked statistics as a practical and powerful tool for extracting non-Gaussian information from weak lensing surveys.

Cowell, Jessica A. [Oxford U.; Tokyo U., IPMU] (OR↗
DOE OSTI2024

Microfabricated Ion Traps on Sapphire for Larger Trap Areas and Higher Qubit Count

Surface ion traps are a promising platform for quantum computing due to their potential to store large numbers of ions that can be addressed by electrical and optical control signals in order to implement quantum algorithms. Increasing the power of the quantum computer requires increasing the number of ions, but this poses a significant challenge in that it leads to a non-linear increase in on-chip power dissipation. The primary contributor to this power scaling in current devices is the capacitance between the radio frequency (RF) electrode and the metal plane that shields the silicon substrate from the RF signals applied to it. Silicon has traditionally been chosen for the substrate material for compatibility with the processing required for multi-metal-level traps. In this work, we address these capacitance and fabrication challenges by replacing the commonly used silicon substrate with an insulating sapphire substrate to fabricate a multi-metal-level ion trap, while still employing common semiconductor manufacturing techniques. This change in substrate allows the design to remove the metal shielding from the device design, reducing the capacitance of the RF electrode. The electrical characteristics of these traps were measured, specifically trap impedance, capacitance, and voltage breakdown, and compared to nearly identical silicon trap devices. Finally, we used laser cutting techniques to shape a sapphire wafer into bowtie shapes matching silicon traps previously fabricated at Sandia National Labs to explore solutions for integrating sapphire substrates into non-rectangular ion trap designs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Tensile Modeling PVC Gels for Electrohydraulic Actuators

Polyvinyl chloride (PVC)-dibutyl adipate (DBA) gels are a fascinating dielectric elastomer actuator showing promise in soft robotics. When actuated with high voltages, the gel deforms towards the anode. A recent application of PVC gels in electrohydraulic actuators motivates elastic and hyperelastic constitutive relationships for tensile loading modes. PVC gels with plasticizer-to-polymer weight ratios of 2:1, 4:1, 6:1, and 8:1 w/w were evaluated. PVC gels exhibit a linear elastic region up to 25% strain. The elastic modulus decreased with increasing plasticizer content from 288.8 kPa, 56.1 kPa, 24.7 kPa, to 11 kPa. Poisson’s ratio also decreased with increasing plasticizer content from 0.42, 0.43, 0.39, to 0.35. We suggest that the decrease in polymer concentration facilitates a weakly interconnected polymer network susceptible to chain slippage that hinders the network response, thus lowering Poisson’s ratio. Our work suggests that PVC gels can be treated as isotropic and incompressible for large strains and hyperelastic modeling; however, highly plasticized gels tend to act less incompressible at small strains. The power scaling law between the elastic modulus and plasticizer weight ratio showed high agreement, making the elastic modulus deterministic for any plasticizer content. The Neo–Hookean, Mooney–Rivlin, Yeoh, Gent, Ogden, and extended tube hyperelastic constitutive models are investigated. The Yeoh model shows the highest feasibility when evaluated up to 3.5 stretch, showing a maximum normalized root-mean-square-error of 6.85%. Together, these findings establish a constitutive basis for PVC-DBA gels, incorporating small strain elasticity, large strain non-linear behavior, and network analysis while providing suggestive insight into the network structure required for accurately modeling the EPIC.

Polymer Science↗