Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Generalized functions”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A minimum assumption approach to MEG sensor array design

Objective. Our objective is to formulate the problem of the magnetoencephalographic (MEG) sensor array design as a well-posed engineering problem of accurately measuring the neuronal magnetic fields. This is in contrast to the traditional approach that formulates the sensor array design problem in terms of neurobiological interpretability the sensor array measurements. Approach. We use the vector spherical harmonics (VSH) formalism to define a figure-of-merit for an MEG sensor array. We start with an observation that, under certain reasonable assumptions, any array of m perfectly noiseless sensors will attain exactly the same performance, regardless of the sensors' locations and orientations (with the exception of a negligible set of singularly bad sensor configurations). We proceed to the conclusion that under the aforementioned assumptions, the only difference between different array configurations is the effect of (sensor) noise on their performance. We then propose a figure-of-merit that quantifies, with a single number, how much the sensor array in question amplifies the sensor noise. Main results. We derive a formula for intuitively meaningful, yet mathematically rigorous figure-of-merit that summarizes how desirable a particular sensor array design is. We demonstrate that this figure-of-merit is well-behaved enough to be used as a cost function for a general-purpose nonlinear optimization methods such as simulated annealing. We also show that sensor array configurations obtained by such optimizations exhibit properties that are typically expected of 'high-quality' MEG sensor arrays, e.g. high channel information capacity. Significance. Our work paves the way toward designing better MEG sensor arrays by isolating the engineering problem of measuring the neuromagnetic fields out of the bigger problem of studying brain function through neuromagnetic measurements.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Transition Metal Dichalcogenide MoS 2 : Oxygen and Fluorine Functionalization for Selective Plasma Processing

Low-temperature plasma processing is a promising technique for tailoring transition metal dichalcogenides (TMDs). For chalcogen substitution processing, a key challenge is to identify the ion energy window that enables selective chalcogen removal while preserving the metal lattice. Using ab initio molecular dynamics (AIMD), we demonstrate that oxygen and fluorine functionalization widen the processing window by significantly lowering the sulfur sputtering energy threshold (E sputt,S ) of MoS 2 from ∼30 to ∼10 eV via formation of sputtering products such as SO 2 and SF n . Additionally, we show that experimentally relevant cryogenic temperatures strongly affect E sputt,S (T). The dependence is confirmed via AIMD and also predicted by a mechanistic parameter-free theory, suggesting that E sputt (T) generalizes to other TMDs, functionalizations, and surface impact conditions. Our results highlight oxygen/fluorine functionalization, ionic impact angle, and material temperature to be key control parameters for selective, damage-controlled chalcogen removal in TMD processing.

Polyachenko, Yury [Princeton Plasma Physics Labora↗
DOE OSTI2024

Magnetism in the two-dimensional dipolar XY model

Here, motivated by a recent experiment on a square-lattice Rydberg atom array realizing a long-range dipolar XY model [C. Chen et al., Nature (London) 616, 691 (2023)], we numerically study the model's equilibrium properties. We obtain the phase diagram, critical properties, entropies, variance of the magnetization, and site-resolved correlation functions. We consider both ferromagnetic and antiferromagnetic interactions and apply quantum Monte Carlo and pseudo-Majorana functional renormalization group techniques, generalizing the latter to a U⁡(1) symmetric setting. Our simulations perform extensive thermometry in dipolar Rydberg atom arrays and establish conditions for adiabaticity and thermodynamic equilibrium. On the ferromagnetic side of the experiment, we determine the entropy per particle S/N≈0.5, close to the one at the critical temperature, S c /N=0.585⁢(15). The simulations suggest the presence of an out-of-equilibrium plateau at large distances in the correlation function, thus motivating future studies on the nonequilibrium dynamics of the system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Atomic Layer Deposition for Materials-Based H2 Storage: Opportunities and Limitations

The transportation demands in our growing hydrogen (H2) economy requires robust storage systems. The current commercially implemented technology in fuel cell cars relies on the well-established technology of hydrogen gas compressed to 350-700 bar, depending on the application. The compressed gas tanks in use today are bulky and cost intensive. To address this challenge, material-based storage is one of the long-term alternatives considered and constitutes the focus of this talk. Material-based storage is broadly defined as hydrogen bound to solid materials, with its binding strength varying from physisorption to porous materials, such as zeolites and metal organic frameworks, to chemisorption in (complex) metal hydrides. The ultimate targets set by the U.S. Department of Energy for this technology include a system gravimetric capacity of 6.5 wt% and volumetric capacity of 40 g/L at 100 bar, operating temperatures ranging between -40 C and +40 C and adsorption/desorption timescales of < 5 min. Storage in the form of physisorbed or chemisorbed hydrogen has guided the materials research, which metal- organic framework and (complex) metal hydrides being the most promising materials classes. A variety of these materials have met one or more of the targets, but it has remained elusive for a single material system to meet all these stringent requirements. Nano-encapsulation and low-concentration chemical additives have previously been employed separately to overcome such challenges. Functionalization via atomic layer deposition (ALD), however, offers unique characteristics that make it suitable for both, nano-encapsulation and "doping" with low-concentration additives. This deposition technique has sub-monolayer thickness control, is highly conformal in high-surface area materials and is self-limiting, i.e., once the gas-phase precursor reacts with the available surface sites, the surface reactions stop. In this presentation, we will showcase examples where ALD, more generally vapor-phase functionalization, on (complex) metal hydrides and organic frameworks has improved the material properties for H2 storage. In our first study, Al2O3 was deposited on magnesium borohydride, Mg(BH4)2, at room temperature using trimethylaluminum (TMA) and water. From our findings, encouraging initial results were obtained: the H2 desorption temperature was lowered by 60-120 degrees C, the desorbed gravimetric H2 capacity at temperatures < 250 degrees C was doubled, and the desorption kinetics increased by a factor of ~6 compared to uncoated Mg(BH4)2. However, hydrolysis reactions caused by residual surface -OH groups from the ALD water-pulse degraded the sample substantially. Through this study, the use of TMA was observed to be highly reactive with the Mg(BH4)2 surface species, leading to the strategy of exposing the sample only to TMA. The relative mole fraction of the vapor-phase additive can be precisely controlled with the number of TMA pulses and pulse duration. By tuning these parameters, we show in our second study that 10 pulses of TMA at ambient conditions were able to decrease the H2 desorption temperature by ~100 degrees C while retaining <95 % of its H2 capacity. We applied this approach with other additives such as BBr3, TiCl4 and tetrahydrofuran, and demonstrate that this unique approach opens the door to a new class of molecular additives and catalysts which cannot easily be introduced with conventional mechano-chemical or solvent-based techniques for (complex) metal-hydrides. In the case of metal organic frameworks, ALD is a promising technique to functionalize the pores with metal atoms or functional groups able to tune the gas selectivity and binding energy, for which ~15 kJ/mol has been established as the optimal value for H2 storage in sorbent materials. Vapor-phase techniques such as ALD have numerous benefits over other functionalization tools for H2 storage materials opening innumerable opportunities to the field. To exploit these opportunities, the current limitations on room temperature and water-less ALD processes needs to be overcome which will greatly expand the possibilities of encapsulation and incorporation of additives for organic frameworks and (complex) metal hydrides.

atomic layer deposition↗
NASA NTRS1985

Effects of temperature on remotely sensed mineral absorption features

The present paper provides the results from a series of laboratory experiments which were conducted to investigate the variable spectral behavior of diagnostic Fe(++) crystal field absorptions in common rock-forming minerals as a function of temperature. The general goal of the study is to refine existing techniques of mineralogic and chemical characterization, based on specific optical behavior which has been measured and correlated in a systematic way. It is pointed out that the presented work, taken together with a number of recent studies, provides new information on the manner in which specific features of mafic minerals combine in rocks and soil spectra, taking into account also variations of identification and subsequent interpretation as a function of planetary surface temperature. The obtained results are directly applicable to interpretation of remotely sensed data for many solar system objects, particularly Mars, asteroids, the moon, earth, and Mercury.

Singer, R. B.↗
DOE OSTI2025

Back to functional hydrides: Effects of neutron-irradiated microstructure on hydrogen retention in yttrium hydride

Functional hydrides are promising candidates for advanced nuclear reactors, particularly in portable or transportable applications, due to their high hydrogen-retention capabilities, enabling efficient neutron moderation, and compact reactor design. However, hydrogen mobility in hydrides at elevated irradiation temperatures poses significant technological challenges, necessitating a comprehensive understanding of their irradiation behavior. Furthermore, this study investigated the microstructural and chemical stability of neutron-irradiated yttrium hydrides to assess their hydrogen-retention capacity. A targeted literature review was also conducted to contextualize neutron-irradiation effects on functional hydrides with regards to structural stability and hydrogen retention. Experimental characterizations revealed that, at high temperatures, irradiated hydrides retained their phase stability, which was likely enhanced by irradiation-induced microstructure evolution. Notably, an amorphous yttrium and oxygen -rich surface layer was present at the free surface of the hydride. Its thickness decreased while a continuous crystalline Y-O-rich layer was formed with increasing neutron damage. Additionally, the number density of dislocation loops and cavities generally increased as a function of neutron dose. First-principles calculations of hydrogen behavior within yttrium vacancy clusters in yttrium hydrides demonstrated vacancy-size-dependent hydrogen stability and configuration, highlighting the role of vacancy geometry in regulating hydrogen retention. Thus, the presence of irradiation-induced dislocation loops and cavities were hypothesized to improve hydrogen retention. Collectively, these findings advance the understanding of hydride behavior under neutron irradiation as well as their technological readiness for portable or transportable nuclear reactors.

Defect clusters↗
NASA NTRS2018

GMI-IPS: Python Processing Software for Aircraft Campaigns

NASA's Atmospheric Tomography Mission (ATom) seeks to understand the impact of anthropogenic air pollution on gases in the Earth's atmosphere. Four flight campaigns are being deployed on a seasonal basis to establish a continuous global-scale data set intended to improve the representation of chemically reactive gases in global atmospheric chemistry models. The Global Modeling Initiative (GMI), is creating chemical transport simulations on a global scale for each of the ATom flight campaigns. To meet the computational demands required to translate the GMI simulation data to grids associated with the flights from the ATom campaigns, the GMI ICARTT Processing Software (GMI-IPS) has been developed and is providing key functionality for data processing and analysis in this ongoing effort. The GMI-IPS is written in Python and provides computational kernels for data interpolation and visualization tasks on GMI simulation data. A key feature of the GMI-IPS, is its ability to read ICARTT files, a text-based file format for airborne instrument data, and extract the required flight information that defines regional and temporal grid parameters associated with an ATom flight. Perhaps most importantly, the GMI-IPS creates ICARTT files containing GMI simulated data, which are used in collaboration with ATom instrument teams and other modeling groups. The initial main task of the GMI-IPS is to interpolate GMI model data to the finer temporal resolution (1-10 seconds) of a given flight. The model data includes basic fields such as temperature and pressure, but the main focus of this effort is to provide species concentrations of chemical gases for ATom flights. The software, which uses parallel computation techniques for data intensive tasks, linearly interpolates each of the model fields to the time resolution of the flight. The temporally interpolated data is then saved to disk, and is used to create additional derived quantities. In order to translate the GMI model data to the spatial grid of the flight path as defined by the pressure, latitude, and longitude points at each flight time record, a weighted average is then calculated from the nearest neighbors in two dimensions (latitude, longitude). Using SciPya's Regular Grid Interpolator, interpolation functions are generated for the GMI model grid and the calculated weighted averages. The flight path points are then extracted from the ATom ICARTT instrument file, and are sent to the multi-dimensional interpolating functions to generate GMI field quantities along the spatial path of the flight. The interpolated field quantities are then written to a ICARTT data file, which is stored for further manipulation. The GMI-IPS is aware of a generic ATom ICARTT header format, containing basic information for all flight campaigns. The GMI-IPS includes logic to edit metadata for the derived field quantities, as well as modify the generic header data such as processing dates and associated instrument files. The ICARTT interpolated data is then appended to the modified header data, and the ICARTT processing is complete for the given flight and ready for collaboration. The output ICARTT data adheres to the ICARTT file format standards V1.1. The visualization component of the GMI-IPS uses Matplotlib extensively and has several functions ranging in complexity. First, it creates a model background curtain for the flight (time versus model eta levels) with the interpolated flight data superimposed on the curtain. Secondly, it creates a time-series plot of the interpolated flight data. Lastly, the visualization component creates averaged 2D model slices (longitude versus latitude) with overlaid flight track circles at key pressure levels. The GMI-IPS consists of a handful of classes and supporting functionality that have been generalized to be compatible with any ICARTT file that adheres to the base class definition. The base class represents a generic ICARTT entry, only defining a single time entry and 3D spatial positioning parameters. Other classes inherit from this base class; several classes for input ICARTT instrument files, which contain the necessary flight positioning information as a basis for data processing, as well as other classes for output ICARTT files, which contain the interpolated model data. Utility classes provide functionality for routine procedures such as: comparing field names among ICARTT files, reading ICARTT entries from a data file and storing them in data structures, and returning a reduced spatial grid based on a collection of ICARTT entries. Although the GMI-IPS is compatible with GMI model data, it can be adapted with reasonable effort for any simulation that creates Hierarchical Data Format (HDF) files. The same can be said of its adaptability to ICARTT files outside of the context of the ATom mission. The GMI-IPS contains just under 30,000 lines of code, eight classes, and a dozen drivers and utility programs. It is maintained with GIT source code management and has been used to deliver processed GMI model data for the ATom campaigns that have taken place to date.

Damon, M. R.↗
DOE OSTI2021

Whetstone v.0.9.2

Whetstone is a deep learning software library for spiking/binary threshold neuromorphic hardware. Built to integrate with Keras models, Whetstone provides a collection of layers, callbacks, and utility functions designed to allow the training of integrate-and-fire neurons rather than traditional analog activation functions (e.g. rectified linear units). Generally, training a network for neuromorphic hardware is challenging due to the discontinuity of the activation function. Whetstone overcomes this challenge by using ever-closer continuous approximations of a threshold activation, modified during training time. This allows rapid and consistent training for ultra-low powered neuromorphic hardware.

Severa, William↗
DOE OSTI2022

A comprehensive and fair comparison of two neural operators (with practical extensions) based on $\mathrm{FAIR}$ data

Neural operators can learn nonlinear mappings between function spaces and offer a new simulation paradigm for real-time prediction of complex dynamics for realistic diverse applications as well as for system identification in science and engineering. Herein, we investigate the performance of two neural operators, which have shown promising results so far, and we develop new practical extensions that will make them more accurate and robust and importantly more suitable for industrial-complexity applications. The first neural operator, DeepONet, was published in 2019 (Lu et al., 2019), and its original architecture was based on the universal approximation theorem of Chen & Chen (1995). The second one, named Fourier Neural Operator or FNO, was published in 2020, and it is based on parameterizing the integral kernel in the Fourier space. DeepONet is represented by a summation of products of neural networks (NNs), corresponding to the branch NN for the input function and the trunk NN for the output function; both NNs are general architectures, e.g., the branch NN can be replaced with a CNN or a ResNet. According to Kovachki et al. (2021), FNO in its continuous form can be viewed conceptually as a DeepONet with a specific architecture of the branch NN and a trunk NN represented by a trigonometric basis. In order to compare FNO with DeepONet computationally for realistic setups, we develop several extensions of FNO that can deal with complex geometric domains as well as mappings where the input and output function spaces are of different dimensions. We also develop an extended DeepONet with special features that provide inductive bias and accelerate training, and we present a faster implementation of DeepONet with cost comparable to the computational cost of FNO, which is based on the Fast Fourier Transform. Here we consider 16 different benchmarks to demonstrate the relative performance of the two neural operators, including instability wave analysis in hypersonic boundary layers, prediction of the vorticity field of a flapping airfoil, porous media simulations in complex-geometry domains, etc. We follow the guiding principles of FAIR (Findability, Accessibility, Interoperability, and Reusability) for scientific data management and stewardship. The performance of DeepONet and FNO is comparable for relatively simple settings, but for complex geometries the performance of FNO deteriorates greatly. We also compare theoretically the two neural operators and obtain similar error estimates for DeepONet and FNO under the same regularity assumptions.

42 ENGINEERING↗
NASA NTRS1996

Refining Linear Fuzzy Rules by Reinforcement Learning

Linear fuzzy rules are increasingly being used in the development of fuzzy logic systems. Radial basis functions have also been used in the antecedents of the rules for clustering in product space which can automatically generate a set of linear fuzzy rules from an input/output data set. Manual methods are usually used in refining these rules. This paper presents a method for refining the parameters of these rules using reinforcement learning which can be applied in domains where supervised input-output data is not available and reinforcements are received only after a long sequence of actions. This is shown for a generalization of radial basis functions. The formation of fuzzy rules from data and their automatic refinement is an important step in closing the gap between the application of reinforcement learning methods in the domains where only some limited input-output data is available.

Berenji, Hamid R.↗
DOE OSTI2023

Shadow energy functionals and potentials in Born–Oppenheimer molecular dynamics

In Born–Oppenheimer molecular dynamics (BOMD) simulations based on the density functional theory (DFT), the potential energy and the interatomic forces are calculated from an electronic ground state density that is determined by an iterative self-consistent field optimization procedure, which, in practice, never is fully converged. The calculated energies and forces are, therefore, only approximate, which may lead to an unphysical energy drift and instabilities. Here, we discuss an alternative shadow BOMD approach that is based on backward error analysis. Instead of calculating approximate solutions for an underlying exact regular Born–Oppenheimer potential, we do the opposite. Instead, we calculate the exact electron density, energies, and forces, but for an underlying approximate shadow Born–Oppenheimer potential energy surface. In this way, the calculated forces are conservative with respect to the approximate shadow potential and generate accurate molecular trajectories with long-term energy stabilities. We show how such shadow Born–Oppenheimer potentials can be constructed at different levels of accuracy as a function of the integration time step, δt, from the constrained minimization of a sequence of systematically improvable, but approximate, shadow energy density functionals. For each energy functional, there is a corresponding ground state Born–Oppenheimer potential. These pairs of shadow energy functionals and potentials are higher-level generalizations of the original “zeroth-level” shadow energy functionals and potentials used in extended Lagrangian BOMD. The proposed shadow energy functionals and potentials are useful only within this extended dynamical framework, where also the electronic degrees of freedom are propagated as dynamical field variables together with the atomic positions and velocities. The theory is quite general and can be applied to MD simulations using approximate DFT, Hartree–Fock, or semi-empirical methods, as well as to coarse-grained flexible charge models.

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

Structural and electronic properties of rare-earth chromites: A computational and experimental study

Here, in this work, the structural, optical, and electronic properties of rare-earth perovskites of the general formula R⁢ CrO 3 , where R⁢ represents the rare-earth Gd, Tb, Dy, Ho, Er, and Tm, have been studied in detail. These compounds were synthesized through a facile citrate route. X-ray diffraction, Raman spectroscopy, and UV-Visible spectroscopy were utilized to reveal the structural evolutions in R⁢ ⁢CrO 3 . The lattice parameters, Cr 3+ –O 2 –Cr 3+ bond angle, and CrO 6 octahedral distortions were found to strongly depend on the ionic radii of R⁢ . First-principles calculations based on density-functional theory within the generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) and strongly constrained and appropriately normed (SCAN) meta-GGA were also employed to calculate the structural and electronic properties of R⁢ ⁢CrO 3 . The ground-state energy, lattice constants, electronic structures, and density of states of R⁢ ⁢CrO 3 were calculated. These provide some insights into the electronic characteristics of the R⁢ ⁢CrO 3 compounds. The calculated values of lattice parameters and band gaps with Hubbard U correction (SCAN+ U ) agree well with values measured experimentally and show more accuracy in predicting the ground-state crystal structure and band structure compared to PBE+ U approximation. The band gap of R⁢ ⁢CrO 3 is found to be independent of the ionic radii of R⁢ from both experiments and calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

The response field and the saddle points of quantum mechanical path integrals

Highlights: • Moyal quantum mechanics and Marinov’s path integral. • Classical and semiclassical limits of Marinov’s path integral. • Oscillating functional integrals. • Instantons of the Marinov’s path integral. In quantum statistical mechanics, Moyal’s equation governs the time evolution of Wigner functions and of more general Weyl symbols that represent the density matrix of arbitrary mixed states. A formal solution to Moyal’s equation is given by Marinov’s path integral. In this paper we demonstrate that this path integral can be regarded as the natural link between several conceptual, geometric, and dynamical issues in quantum mechanics. A unifying perspective is achieved by highlighting the pivotal role which the response field, one of the integration variables in Marinov’s integral, plays for pure states even. The discussion focuses on how the integral’s semiclassical approximation relates to its strictly classical limit; unlike for Feynman type path integrals, the latter is well defined in the Marinov case. The topics covered include a random force representation of Marinov’s integral based upon the concept of “Airy averaging”, a related discussion of positivity-violating Wigner functions describing tunneling processes, and the role of the response field in maintaining quantum coherence and enabling interference phenomena. The double slit experiment for electrons and the Bohm–Aharonov effect are analyzed as illustrative examples. Furthermore, a surprising relationship between the instantons of the Marinov path integral over an analytically continued (“Wick rotated”) response field, and the complex instantons of Feynman-type integrals is found. The latter play a prominent role in recent work towards a Picard–Lefschetz theory applicable to oscillatory path integrals and the resurgence program.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1980

Demonstration Advanced Avionics System (DAAS) functional description

A comprehensive set of general aviation avionics were defined for integration into an advanced hardware mechanization for demonstration in a Cessna 402B aircraft. Block diagrams are shown and system and computer architecture as well as significant hardware elements are described. The multifunction integrated data control center and electronic horizontal situation indicator are discussed. The functions that the DAAS will perform are examined. This function definition is the basis for the DAAS hardware and software design.

Source record↗
DOE OSTI2025

Range-separated thermal hybrid exchange-correlation density functional for accurate band-gap calculations of warm dense matter

Here, in this study, we present a range-separated thermal hybrid exchange-correlation (XC) functional, termed RS-KDT0, which provides highly accurate band-gap predictions at elevated temperature while reducing the computational cost relative to thermal global hybrid functionals. At zero temperature, RS-KDT0 recovers the widely used ground-state range-separated hybrid functional HSE06. To assess its performance across temperature regimes, we conducted static calculations of electronic pressure, band gaps, and band structures. RS-KDT0 exhibits the same level of accuracy as HSE06 in the low-temperature regime and demonstrates improved agreement with finite-temperature GW calculations at elevated temperatures when compared to other thermal XC functionals such as the generalized gradient approximation (GGA)-level KDT16 and the global hybrid KDT0.

density functional theory↗
NASA NTRS1972

Computer generation of symbolic network functions - A new theory and implementation.

A new method is presented for obtaining network functions in which some, none, or all of the network elements are represented by symbolic parameters (i.e., symbolic network functions). Unlike the topological tree enumeration or signal flow graph methods generally used to derive symbolic network functions, the proposed procedure employs fast, efficient, numerical-type algorithms to determine the contribution of those network branches that are not represented by symbolic parameters. A computer program called NAPPE (for Network Analysis Program using Parameter Extractions) and incorporating all of the concepts discussed has been written. Several examples illustrating the usefulness and efficiency of NAPPE are presented.

Alderson, G. E.↗
DOE OSTI2022

A Decision-Making Machine Learning Approach in Hermite Spectral Approximations of Partial Differential Equations

The accuracy and effectiveness of Hermite spectral methods for the numerical discretization of partial differential equations on unbounded domains are strongly affected by the amplitude of the Gaussian weight function employed to describe the approximation space. This is particularly true if the problem is under-resolved, i.e., there are no enough degrees of freedom. The issue becomes even more crucial when the equation under study is time-dependent, forcing in this way the choice of Hermite functions where the corresponding weight depends on time. In order to adapt dynamically the approximation space, it is here proposed an automatic decision-making process that relies on machine learning techniques, such as deep neural networks and support vector machines. The algorithm is numerically tested with success on a simple 1D problem, but the main goal is its exportability in the context of more serious applications. Here we also show at the end an application in the framework of plasma physics.

97 MATHEMATICS AND COMPUTING↗