Search NASA⌕ Search

SEARCH · Search NASA

Results for “high non-linearity”

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 OSTI2024

Observation of critical scaling in spin glasses below T c using thermoremanent magnetization

Time-dependent thermoremanent magnetization (TRM) studies have been instrumental in probing energy dynamics within the spin glass phase. In this paper, we review the evolution of the TRM experiment over the last half century and discuss some aspects related to how it has been used in the understanding of spin glasses. We also report on recent experiments using high-resolution DC SQUID magnetometry to probe the TRM at temperatures less than but near to the transition temperature T c . These experiments have been performed as a function of waiting time, temperature, and five different magnetic fields. We find that as the transition temperature is approached from below, the characteristic time scale of TRM is suppressed up to several orders of magnitude in time. In the highest-temperature region, we find that the waiting time effect subsides, and a waiting time-independent crossover line is reached. We also find that increasing the magnetic field further suppresses the crossover line. Using a first-principles energy argument across the crossover line, we derive an equation that is an excellent fit to the crossover lines for all magnetic fields probed. The data show strong evidence for critical slowing down and an H = 0 Oe phase transition.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Impact of U-10Mo HALEU fuel element tolerances on the Massachusetts Institute of Technology reactor safety and operational performance – Thermal hydraulics

The U.S. is coordinating efforts for the conversion of six U.S. High Performance Research Reactors (USHPRR) including one critical facility from highly enriched uranium (HEU) to low-enriched uranium (LEU). In order to continue the mission of these reactors, including the Massachusetts of Institute of Technology Reactor (MITR), and achieve similar performance, high-assay low-enriched uranium (HALEU) with a high-density metallic alloy of uranium with 10 wt% molybdenum (U-10Mo) is being evaluated. The impact of the fabrication specification and tolerances was assessed following the preliminary design of the MITR LEU fuel elements using the U-10Mo monolithic alloy. This research focuses on the analysis of fabrication specification impact on thermal hydraulics (TH) characteristic of the MITR LEU core as a function of the variation of the relevant fuel specification parameters (e.g., coolant channel gap thickness, fuel plate thickness, etc.). The analyses are performed based on an all-fresh LEU fuel conversion plan identified in a preliminary safety analysis report submitted to the Nuclear Regulatory Commission. The reactor power margin to the onset of nucleate boiling (ONB) is assessed under the limiting safety system settings (LSSS), where a scram occurs, to ensure there is sufficient margin to the reactor safety limit, which is defined by the onset of flow instability that occurs after the ONB. The best estimate plus uncertainty approach is employed to analyze this TH characteristic, which yields realistic results while maintaining adequate conservatism, utilizing a statistical uncertainty propagation method with the STAT7 code. The TH characteristic is analyzed as a function of the variability of the specification parameters resulting from the fabrication process. The main findings of this study show that the MITR core can meet the TH safety and operational requirements at the all-LEU initial core startup (cycle 1), selected transition cycles (most reactive cycle and most limiting cycle: cycle 3 and 5, respectively) and equilibrium (cycle 14) cores under all limiting fabrication parameter combinations considered. In addition, the analyses show that the dependency of the core power margin to ONB on those specification parameters that have the most direct impact on TH performance is non-linear but monotonically decreasing within the specification tolerances. The third order polynomial fit curves are reported in detail for selected limiting cases and can serve as a powerful tool for future MITR fuel management in cases such as when HALEU supply is established that may allow additional cycle length or other operational benefits.

Conversion↗
DOE OSTI2025

Cosmology From CMB Lensing and Delensed EE Power Spectra Using 2019-2020 SPT-3G Polarization Data

From CMB polarization data alone we reconstruct the CMB lensing power spectrum, comparable in overall constraining power to previous temperature-based reconstructions, and an unlensed E -mode power spectrum, with clear detections of the third through tenth acoustic peaks. The observations, taken in 2019 and 2020 with the South Pole Telescope (SPT) and the SPT-3G camera, cover 1500 deg 2 at 95, 150, and 220 GHz with arcminute resolution and roughly 4.9 µ K-arcmin coadded noise in polarization. The power spectrum estimates, together with systematic parameter estimates and a joint covariance matrix, follow from a Bayesian analysis using the Marginal Unbiased Score Expansion (MUSE) method. The E -mode spectrum at ℓ > 2000 and lensing spectrum at L > 350 are the most precise to date. Assuming the ΛCDM model, and using only these SPT data and priors on τ and absolute calibration from Planck, we find H 0 = 66.81 ± 0.81 km/s/Mpc, comparable in precision to the Planck determination and in 5.4 σ tension with the most precise H 0 inference derived via the distance ladder. We also find S 8 ≡ σ 8 (Ω m /0.3) 0.5 = 0.850 ± 0.017, providing further independent evidence of a slight tension with low-redshift structure probes. The ΛCDM model provides a good simultaneous fit to the combined Planck, ACT, and SPT data, and thus passes a powerful test. Combining these CMB datasets with BAO observations, we explore extensions to the ΛCDM model. We find that the effective number of neutrino species, spatial curvature, and primordial helium fraction are consistent with standard model values, and that the 95% confidence upper limit on the neutrino mass sum is 0.075 eV, close to the minimum sum expected from observations of solar and atmospheric neutrino oscillations. The SPT data are consistent with the somewhat weak (< 3 σ ) preference for excess lensing power seen in Planck and ACT data relative to predictions of the ΛCDM model given the combined Planck, ACT, and BAO data sets. Finally, we also detect at greater than 3 σ the influence of non-linear evolution in the CMB lensing power spectrum and discuss it in the context of the S 8 tension. Forthcoming SPT-3G analyses will feature deeper and wider observations in temperature and polarization, providing even tighter constraints and more powerful tests of the ΛCDM model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Smooth splitting and zeros from on-shell recursion

We describe a new approach to understanding the origins of recently discovered “hidden zeros” and “smooth splitting” of tree-level amplitudes in Tr ϕ 3 , Non-Linear Sigma Model (NLSM), Yang-Mill-Scalar (YMS) and the special Galileon. Introducing a new type of linear shift in kinematic space we demonstrate that the mysterious splitting formulae follow from a simple contour integration argument in the style of on-shell recursion. The argument makes use of only standard notions of tree-level factorization on propagators, but assumes improved UV behavior in the form of the absence of a residue at infinity. In the case of Tr ϕ 3 and NLSM this is proven by identifying our shift as a special case of a more general construction called a g-vector shift; in the case of YMS it remains an unproven conjecture. This recursive perspective leads to numerous new results: we derive generalizations of the splitting formulae on more relaxed near-zero kinematics, including interesting new kinematic limits in which the amplitude splits into a triple-product; we also demonstrate that the uncolored special Galileon model has improved UV scaling and hence also splits. We also investigate the possible realization of hidden zeros in four dimensions. The conditions under which the dimensionality constraints are compatible with zero kinematics is investigated in detail for Tr ϕ 3 and YMS; for the latter we find they can be realized only with certain restrictions on external helicity states. The realizable 4d zeros are proven by a similar recursive argument based on BCFW and is found to generalize to a new class of intrinsically 4d “helicity zeros” present in all sectors of YM and also gravity.

effective field theories↗
DOE OSTI2025

Horn Location Sensors for LBNF

The Horn Location Sensor (HLS) system—named for the magnetic focusing horns it monitors, along with other critical beamline elements – is a high-precision alignment system developed for the Long-Baseline Neutrino Facility (LBNF) to support the Deep Underground Neutrino Experiment (DUNE). With minimal maintenance, the HLS can operate reliably in environments with high radiation, and it ensures that key components such as the protective baffle, focusing horns, and beam position monitors are aligned correctly – each essential for maintaining one of the world’s most intense muon-neutrino beams. A high-precision hydrostatic level sensor, a linear variable differential transformer, and INVAR rods are all used in the system to monitor vertical motion and tilt with sub-millimeter accuracy. A novel Sweep Tracker interferometer enhances calibration fidelity by correcting for non-linearities in laser wavelength and scan rate in real time. A critical part of DUNE’s precision alignment and flux prediction requirements, the HLS system initially supports beam power of up to 1.2 MW and can be upgraded to 2.4 MW.

Frequency scanning interferometry↗
DOE OSTI2025

Pore Resolved Simulations of Joule Heating in Fibrous Media using an Embedded Boundary Method

Joule heating has been regarded as an energy-efficient and sustainable method for heating materials and gases at large scales. The modeling of local temperature effects at pore-resolved scales for such systems, however, has been difficult to achieve due to challenges in coupling thermo-chemical processes in complex porous media and in large representative volume elements (RVEs). To this end, we developed an electro-thermal model at the pore scale to study Joule heating effects in large heterogeneous systems with different microstructures. This was achieved using the level set method to implicitly delineate distinct regions within the domain, and an embedded boundary method to facilitate heat exchange across the fluid-solid interface. Moreover, we applied this method to investigate unsteady non-linear electro-thermal effects in non-woven fibrous graphite conductors for RVEs with characteristic lengths of 2 mm, with different fiber orientations, porosity (80% – 90%) and fiber diameters (10 – 20µm). The coupled equations were solved numerically and they produced peak temperatures greater than 2000 K resulting in heating rates as high as 80,000 K/s. Moreover, the results depended strongly on the microstructure of the fiber skeleton and current density. Geometries with large fibers (∼ 20µm) had the highest average and peak temperatures with the mean temperature increasing by 3.9 % while the peak temperature increased by 9.9 %. Anisotropic domains on the other hand had the lowest mean and peak temperatures with peak and mean temperatures of 2293 K and 1437.7K respectively representing a corresponding 12.1% and 5.1% drop in the temperatures. An increase in porosity from 80% to 90%, however, led to an increase in the peak temperature by 5.1%.

Joule heating↗
DOE OSTI2024

GPAW: An open Python package for electronic structure calculations

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe–Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn–Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Simulations of the churning mode: Toroidally symmetric plasma convection and turbulence around the X-points in a snowflake divertor

Using a reduced MHD model, extended to include field-aligned thermal conduction, we present numerical simulations of the churning mode (CM): a toroidally symmetric, non-linear plasma vortex in the vicinity of the null points in a snowflake (SF) divertor [D. D. Ryutov et al., Phys. Scr. 89, 088002 (2014)]. Simulations are carried out across a range of inter-null separations, d xx ⁠, and inter-null orientations, θ⁠, primarily in conditions relevant to the MAST-U tokamak. We find that, when d xx is small, the CM induces additional transport across the X-points when β pm ≳8%, where β pm is the ratio of the plasma pressure in the null region to poloidal magnetic pressure at the midplane. This transport also increases approximately linearly as d xx is reduced. A diffusive model of this transport is shown to predict the total transport across the null points, where diffusion coefficients of up to ∼ 10 2 m 2 s −1 centered on a small region around the X-points are used. However, the CM also results in significant changes to the flux surfaces in the null region, which is not captured by this diffusive model. The changes in magnetic geometry mean the fractional exhaust power delivered to each divertor leg is highly sensitive to β pm ⁠, d xx ⁠, and θ⁠. For small values of θ⁠, the CM can induce a change in topology, redirecting exhaust power from a secondary divertor leg on the high field side to one on the low field side. Similar behavior is found in the fraction of exhaust power going to the inner and outer divertor. Such changes in the flux surfaces may not be captured by Grad–Shafranov solvers and so may be a source of error in the magnetic reconstruction of SF experiments. We consistently find that the fractional exhaust power going to a secondary divertor leg on the high field side is small, consistent with SF experiments.

Computer simulation↗
DOE OSTI2025

Automated Resonance Fitting for Nuclear Data Evaluation

Global and national efforts to deliver high-quality nuclear data to users have a wide-ranging impact, affecting applications in national security, reactor operations, basic science, medicine, and more. Cross section evaluation is a major part of this effort, combining theory and experimentation to produce recommended values and uncertainties for reaction probabilities. Resonance region evaluation is a specialized type of nuclear data evaluation that can require significant manual effort and months of time from expert scientists. In this article, non-convex non-linear optimization methods are combined with concepts of inferential statistics to infer a resonance model from experimental data in an automated manner that is not dependent on prior evaluation(s). This methodology aims to enhance the workflow of a resonance evaluator by minimizing time, effort, and the potential for bias from prior assumptions, while enhancing reproducibility and documentation, thereby addressing well-known challenges in the field.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Climate change, malaria and neglected tropical diseases: a scoping review

Abstract To explore the effects of climate change on malaria and 20 neglected tropical diseases (NTDs), and potential effect amelioration through mitigation and adaptation, we searched for papers published from January 2010 to October 2023. We descriptively synthesised extracted data. We analysed numbers of papers meeting our inclusion criteria by country and national disease burden, healthcare access and quality index (HAQI), as well as by climate vulnerability score. From 42 693 retrieved records, 1543 full-text papers were assessed. Of 511 papers meeting the inclusion criteria, 185 studied malaria, 181 dengue and chikungunya and 53 leishmaniasis; other NTDs were relatively understudied. Mitigation was considered in 174 papers (34%) and adaption strategies in 24 (5%). Amplitude and direction of effects of climate change on malaria and NTDs are likely to vary by disease and location, be non-linear and evolve over time. Available analyses do not allow confident prediction of the overall global impact of climate change on these diseases. For dengue and chikungunya and the group of non-vector-borne NTDs, the literature privileged consideration of current low-burden countries with a high HAQI. No leishmaniasis papers considered outcomes in East Africa. Comprehensive, collaborative and standardised modelling efforts are needed to better understand how climate change will directly and indirectly affect malaria and NTDs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Optimal Control of SOEC-Based Hydrogen Production Systems for Demand Response Using Deep Reinforcement Learning in Smart Grids

Solid oxide electrolysis cell (SOEC) hydrogen production technology can range in size from small, appliance-size equipment to large-scale, central production facilities that can be tied directly to renewable or non-greenhouse-gas-emitting forms of electricity production, making it an ideal resource for demand response (DR). The SOEC hydrogen production system is a complex integrated system that encompasses fluid dynamics, electrical dynamics, and electrochemical and thermal dynamics, all of which involve non-linearity and non-convexity. Proper control of the SOEC hydrogen production system is crucial to enable its participation in the DR program. Here, to overcome the difficulty of designing an explicit control law for such nonlinear systems with nonconvex optimization features in DR applications, deep reinforcement learning (DRL) is explored to achieve the optimal control of the SOEC system for DR participation. Specifically, a twin delayed deterministic policy gradient (TD3) control framework is applied to achieve optimal response performance during DR events by considering power tracking error and hydrogen production efficiency with a suitable reward function. Two case studies with grid connections for tracking different DR commands were investigated. The first case study involved operating conditions reaching the boundaries, while the second involved operating conditions within the boundaries. The results showed that the proposed DRL-based control for SOEC can track the DR signal in a timely manner while maintaining high energy efficiency.

08 HYDROGEN↗
DOE OSTI2025

Impact of high-temperature annealing on hafnia-silica composite coatings deposited via ion beam sputtering for high-peak power 1064 nm lasers

The maximum power handling fluence of high-peak and average power laser systems is often limited by the laser damage of the coatings on optical components. Furthermore, these multilayer dielectric coatings are limited in their maximum power handling due to laser-damage-prone defects in the lower optical bandgap, higher optical index material. Some of these defects can be mitigated by thermal annealing to high temperatures, which can greatly reduce the linear absorbative precursors. Typically, hafnia and silica are the materials of choice for high-peak and average power laser systems in the ultraviolet through infrared spectral range; however, hafnia crystalizes readily when annealed at high temperatures. In this study, we prepare composite HfO 2 -SiO 2 coatings by co-sputtering hafnia and silica in an ion beam sputtering system and compare them to pure hafnia-based coatings. We demonstrate that crystallinity in hafnia can be completely suppressed when it is mixed with silica, such as the composite coatings in this study. High reflectors were fabricated and annealed, demonstrating that the multilayer dielectric stacks can survive high-temperature annealing and exhibit an excellent linear absorption of 0.2 +/- 1 ppm at 1064 nm. Short- and long-pulse laser damage was explored, demonstrating the complex relationship between linear absorption and the non-linear absorption which drives pulsed laser damage. These results provide an excellent route to the creation of very low linear absorption optical coatings, which also utilize low scattering materials that are best suited for high-peak and average power applications.

Harthcock, Colin [Lawrence Livermore National Labo↗
DOE OSTI2025

MLSPICE: Machine Learning based SPICE Modeling Platform for Power Magnetics

Electrical power converters are critical to a wide range of applications ranging from renewable integration to transportation electrification, and can be a key factor determining the size, weight, and efficiency of energy conversion systems. Magnetic components are typically the largest and least efficient components in power electronics. While there have been major strides in the modeling and analysis of power semiconductor devices and circuit simulations, the necessary advances in the design of power magnetics have lagged. In this project, we have transformed the modeling and design of power magnetics with machine learning enabled methods and catalyze simultaneous disruptive improvements for ML-based power electronics design tools. A fully automated open-source machine learning based magnetics modeling platform – the MagNet project - with innovations in full stack have been developed to greatly accelerate the design process and provide new insights to magnetic material and geometry design. The ARPA-E funded MagNet platform contains three major building blocks: 1) a ML-Integrated Data Acquisition System (MIDAS): a highly automated data acquisition testbed which is capable of measuring a large number of magnetic cores with a wide range of electrical circuit excitations; 2) a ML-integrated Core Loss Model (MICLM): a machine-learning trained modeling method for modeling the core loss and saturation effects of magnetic materials for arbitrary excitation waveforms; 3) ML-guided Magnetics SPICE Simulation Tool (PMSPICE): a fully integrated CAD tool which can simulate the magnetics in SPICE. It can help the designers to quickly model the linear and non-linear characteristics of magnetic components and evaluate their behavior in SPICE simulations. The developed MagNet system has fully demonstrated the proposed performance target and has been open sourced to the entire power electronics community to advance the modeling and design of power magnetics from many different angles.

36 MATERIALS SCIENCE↗
DOE OSTI2025

SAM Theory Manual

The System Analysis Module (SAM) is an advanced and modern system analysis tool under development at Argonne National Laboratory for advanced non-LWR reactor safety analysis. It aims to provide fast-running, modest-fidelity, whole-plant transient analyses capabilities, which are essential for fast turnaround design scoping and engineering analyses of advanced reactor concepts. While SAM is being developed as a system-level modeling and simulation tool, advanced modeling techniques being implemented include a reduced-order three-dimensional module, pseudo 3-D conjugate heat transfer modeling in reactor core, flexible and multi-scale modeling of heat transfer between fluid and structures, in addition to the advances in software environments and design, and numerical methods. SAM aims to be a generic system-level safety analysis tool for advanced non-LWRs, including Liquid-Metal-cooled fast Reactors (LMR), Molten Salt Reactors (MSR), Fluoride-salt-cooled High-temperature Reactors (FHR), and High-Temperature Gas-cooled Reactors (HTGR). SAM takes advantage of advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability. It utilizes an object-oriented computational framework (MOOSE), and its underlying meshing and finite-element library and linear and non-linear solvers, to leverage the modern advanced software environments and numerical methods. This document provides the theoretical and technical basis of the code to help users understand the underlying physical models (such as governing equations, closure models, and component models), system modeling approaches, numerical discretization and solution methods, and the overall capabilities in SAM. As new code capabilities and features are added, the SAM Theory Manual will be updated periodically to keep it consistent with the state of the development.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

SAM User's Guide

The System Analysis Module (SAM) is a modern system analysis tool being developed at Argonne National Laboratory for advanced non-LWR safety analysis. It aims to provide fast-running, whole-plant transient analyses capability with improved-fidelity for Sodium-cooled Fast Reactors (SFR), Lead-cooled Fast Reactors (LFR), and Molten Salt Reactors (MSR) or Fluoride-cooled High-temperature Reactors (FHR). SAM takes advantage of advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability. It utilizes an object-oriented application framework (MOOSE), and its underlying meshing and finite-element library (libMesh) and linear and non-linear solvers (PETSc), to leverage the modern advanced software environments and numerical methods. This document provides a user’s guide, which will help users understand the input description and core capabilities of the SAM code. A brief overview of the code is presented, as well as how to obtain and run it. The input syntax for various parts of the code is provided. Additionally, a number of example problems, starting with simple unit component problems to problems with increasing complexity, are provided. Because the code is still under active development, this SAM User’s Guide will evolve with periodic updates.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
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 OSTI2024

Relations between Shot Noise, Gain Bandwidth, and Saturation of Instabilities

There are numerous instabilities present in charged particle beams that undergo exponential growth and reach saturation. In various applications, such as free-electron lasers or micro-bunching light sources, achieving saturation is desirable. Conversely, there are applications where these instabilities are utilized as linear broad-band amplifiers for signals embedded in the charged beam. In the latter scenario, the saturation of an instability induces non-linear distortions in the imprinted signal, thereby limiting the useful range of such amplifiers. Accurate evaluation of these instabilities necessitates a complete and comprehensive modeling approach that includes shot noise within the beam. Unfortunately, such modeling is not always feasible or practical. In this paper, we introduce a methodology utilizing the frequency and bandwidth of the instability as key parameters. Through this, we derive an estimation for the range of linear instability growth. Our derivation is conducted in a model-independent manner, making it applicable to a broad spectrum of instabilities. To validate our approach, we employ established and thoroughly benchmarked simulations with a free electron laser (FEL) code as well as self-consistent 3-dimensional simulation of plasma-cascade instability using code SPACE.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Computational framework for modeling of physical process

Techniques, systems, and devices are described for providing a computational frame for estimating high-dimensional stochastic behaviors. In one exemplary aspect, a method for performing numerical estimation includes receiving a set of measurements of a stochastic behavior. The set of correlated measurements follows a non-standard probability distribution and is non-linearly correlated. Also, a non-linear relationship exists between a set of system variables that describes the stochastic behavior and a corresponding set of measurements. The method includes determining, based on the set of measurements, a numerical model of the stochastic behavior. The numerical model comprises a feature space comprising non-correlated features corresponding to the stochastic behavior. The non-correlated features have a dimensionality of M and the set of measurements has a dimensionality of N, M being smaller than N. The method includes generating a set of approximated system variables corresponding to the set of measurements based on the numerical model.

Chen, Xiao↗