SEARCH · Engineering Papers

Results for “Nonlinear processes”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Probabilistic-learning-based stochastic surrogate model from small incomplete datasets for nonlinear dynamical systems

We consider a high-dimensional nonlinear computational model of a dynamical system, parameterized by a vector-valued control parameter, in the presence of uncertainties represented by an uncontrolled parameter modeled by a vector-valued random variable, and possibly with stochastic excitation. The objective is to construct a statistical surrogate model where the input is any deterministic value of the control parameter, and the output is a vector-valued observation of the computational model, which is a random vector whose probability measure is updated using a target dataset. To construct this statistical surrogate model, the stochastic response of the computational model must be built, which is a vector-valued time-discretized stochastic process in high dimension, depending on the control parameter. It is assumed that the computational cost of a single evaluation of the deterministic model is high. For the probabilistic updating, we consider a subset of the components of the observation of the computational model, defined as the “identification observation” of the computational model, for which a small target dataset is available. Therefore, the target dataset is associated with partial observability, corresponding to an incomplete data case. Given a prior probability model of the random control and uncontrolled parameters, a training dataset is constructed, consisting of realizations of the random triplet composed of the stochastic response, the random identification observation, and the random control parameter. Since the computational cost of a single evaluation of the deterministic model is assumed to be large, the training dataset is also of small size. The main challenges in this problem are the high dimensionality, partial observability leading to incomplete data in the target dataset for the identification observation of the computational model (which is not sufficient to identify the computational stochastic responses), and the availability of a small training dataset. To address these challenges, we propose a methodology based on statistical methods for constructing necessary reduced representations, direct probabilistic learning under constraints using probabilistic learning on manifolds (PLoM) constrained by the target dataset, and the use of a weak formulation of the Fourier transform of probability measures. Statistical conditioning is also employed to explore the learned dataset. The constructed predictive statistical surrogate model can be implemented in the context of online computation. Here, we apply this approach to a problem of nonlinear stochastic dynamics in high dimensions within the framework of deformable solids mechanics.

Engineering↗

DOE OSTI2022

Determination of the nonlinear thermo-optic coefficient of silicon nitride and oxide using an effective index method

There is little literature characterizing the temperature-dependent thermo-optic coefficient (TOC) for low pressure chemical vapor deposition (LPCVD) silicon nitride or plasma enhanced chemical vapor deposition (PECVD) silicon dioxide at temperatures above 300 K. In this study, we characterize these material TOC’s from approximately 300-460 K, yielding values of (2.51 ± 0.08) · 10 −5 K −1 for silicon nitride and (5.67 ± 0.53) · 10 −6 K −1 for silicon oxide at room temperature (300 K). We use a simplified experimental setup and apply an analytical technique to account for thermal expansion during the extraction process. We also show that the waveguide geometry and method used to determine the resonant wavelength have a substantial impact on the precision of our results, a fact which can be used to improve the precision of numerous ring resonator index sensing experiments.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Upper-Hybrid Waves Driven by Meandering Electrons Around Magnetic Reconnection X Line

Magnetic reconnection is a fundamental process in collisionless space plasma environment, and plasma waves relevant to the kinetic interactions can have a significant impact on the multiscale behavior of reconnection. Here, we present Magnetospheric Multiscale (MMS) observations during an encounter of an X line of symmetric magnetic reconnection in the magnetotail. The X line is characterized by reversals of ion and electron jets and electromagnetic fields, agyrotropic electron velocity distribution functions (VDFs), and an electron-scale current sheet. MMS observe large-amplitude nonlinear upper-hybrid (UH) waves on both sides of the neutral line, and the wave amplitudes have highly localized distribution along the normal direction. Here, the inbound meandering electrons drive the UH waves, releasing the free energy stored from the reconnection electric field along the meandering trajectories. The interaction between the meandering electrons and the UH waves may modify the balance of the reconnection electric field around the X line.

58 GEOSCIENCES↗

DOE OSTI2022

Investigating the Influence of Ni, ZrO 2 , and Y 2 O 3 from SOFC Anodes on Siloxane Deposition

Siloxanes, as a type of impurity in biogas, can poison the Ni-YSZ anode of SOFCs. However, the influence of individual components of the anode, such as Ni, ZrO 2 , and Y 2 O 3 , on the siloxane deposition process has not been investigated extensively. In this study, Ni, ZrO 2 , and Y 2 O 3 pellets were exposed to H 2 + N 2 + H 2 O + D4 (octamethylcyclotetrasiloxane, 2.5 ppmv) and H 2 + N 2 + D4 (2.5 ppmv siloxane) gas mixtures at 750 °C to investigate their affinity and tolerance for siloxane degradation. Surface morphology analysis and electrochemical analysis including electrochemical impedance spectroscopy (EIS), related distribution of relaxation times (DRT) analysis and equivalent circuit modeling with complex nonlinear least square (CNLS) fitting were conducted. Here, a microstructure parameter—tortuosity factor to porosity ratio $\tau /\varepsilon $ calculated by diffusion polarization resistance was utilized for siloxane deposition evaluation. After comparing pellets surface morphology changes before and after experiments and $\tau /\varepsilon $ change following the contamination, Ni is considered as a major factor in siloxane deposition reactions in Ni-YSZ anode.

25 ENERGY STORAGE↗

DOE OSTI2021

Parallel projection—An improved return mapping algorithm for finite element modeling of shape memory alloys

Here, we present a novel finite element analysis of inelastic structures containing Shape Memory Alloys (SMAs). Phenomenological constitutive models for SMAs lead to material nonlinearities, that require substantial computational effort to resolve. Finite element analysis methods, which rely on Gauss quadrature integration schemes, must solve two sets of coupled differential equations: one at the global level and the other at the local, i.e. Gauss point level. In contrast to the conventional return mapping algorithm, which solves these two sets of coupled differential equations separately using a nested Newton procedure, we propose a scheme to solve the local and global differential equations simultaneously. In the process we also derive closed-form expressions used to update the internal/constitutive state variables, and unify the popular closest-point and cutting plane methods with our formulas. Numerical testing indicates that our method allows for larger thermomechanical loading steps and provides increased computational efficiency, over the standard return mapping algorithm.

42 ENGINEERING↗

DOE OSTI2025

Optimizing the optimizer for physics-informed neural networks and Kolmogorov-Arnold networks

Physics-Informed Neural Networks (PINNs) have revolutionized the computation of PDE solutions by integrating partial differential equations (PDEs) into the neural network’s training process as soft constraints, becoming an important component of the scientific machine learning (SciML) ecosystem. More recently, physics-informed Kolmogorv-Arnold networks (PIKANs) have also shown to be effective and comparable in accuracy with PINNs. In their current implementation, both PINNs and PIKANs are mainly optimized using first-order methods like Adam, as well as quasi-Newton methods such as BFGS and its low-memory variant, L-BFGS. However, these optimizers often struggle with highly nonlinear and non-convex loss landscapes, leading to challenges such as slow convergence, local minima entrapment, and (non)degenerate saddle points. In this study, we investigate the performance of Self- Scaled BFGS (SSBFGS), Self-Scaled Broyden (SSBroyden) methods and other advanced quasi-Newton schemes, including BFGS and L-BFGS with different line search strategies. These methods dynamically rescale updates based on historical gradient information, thus enhancing training efficiency and accuracy. We systematically compare these optimizers – using both PINNs and PIKANs – on key challenging PDEs, including the Burgers, Allen-Cahn, Kuramoto-Sivashinsky, Ginzburg-Landau, and Stokes equations. Additionally, we evaluate the performance of SSBFGS and SSBroyden for Deep Operator Network (DeepONet) architectures, demonstrating their effectiveness for data-driven operator learning. Our findings provide state-of-the-art results with orders-of-magnitude accuracy improvements without the use of adaptive weights or any other enhancements typically employed in PINNs. More broadly, our work reveal insights into the effectiveness of quasi-Newton optimization strategies in significantly improving the convergence and accurate generalization of PINNs and PIKANs.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Numerical study of coronal plasma jet formation

A new scenario for solar flare eruption in the coronal holes is analyzed by using MHD stability concepts for a spheromak configuration. The stability properties of a spheromak partially embedded into a conducting surface are studied using three dimensional MHD simulations. In agreement with the analytical theory, a large degree of line-tying stabilizes the spheromak's tilt instability, while the elongation has a destabilizing effect. High-resolution nonlinear simulations also demonstrate current sheet formation at the upper surface of the spheromak, where the tilted magnetic field of the spheromak reconnects with the background magnetic field. The calculated stability threshold and the observed magnetic reconnection support a model of coronal jet eruptions where a dome-like magnetic structure grows through flux emergence on the solar surface, tilts, reconnects, and erupts. Countering the effect from elongation, line-tying strongly stabilizes a spheromak growing from a flux-emergence process, suggesting that to accelerate the onset of eruptive coronal jets, there must be magnetic reconnection at the bottom of the spheromak to detach the structure from the solar surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

A Two-Stage Quantum Reinforcement Learning Method for Multi-Objective Transmission Switching

Multi-objective transmission switching (MO-TS) problems involve the strategic reconfiguration of network topology to simultaneously optimize multiple objectives. As the system scale increases, finding feasible solutions becomes increasingly challenging due to the problem's nonlinearity and high computational complexity. To address these challenges, this paper proposes a two-stage quantum reinforcement learning method that leverages potential quantum advantages for MO-TS. In the first stage, candidate switching lines are identified using a graph-theoretical approach to reduce the problem's dimensionality. The second stage introduces a quantum-classical reinforcement learning framework, where a learnable measurement-based CNN-ResVQC architecture is developed to effectively reduce the input dimension for quantum processing, mitigate vanishing gradients, and enhance trainability while improving the quantum circuit's flexibility in modeling complex decision policies for MO-TS. Numerical studies on IEEE 14-bus, 57-bus, and 118-bus systems demonstrate that the proposed algorithm achieves superior training stability and faster convergence with approximately 1% of the network parameters required by classical algorithms, highlighting its effectiveness, efficiency, and scalability. Furthermore, the practicality is validated through its stable convergence under three common quantum noise channels.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2024

Enhancing the Range and Reliability of the Spacer Layer Imaging Method

The spacer layer imaging method (SLIM) is widely used to measure the thickness of additive and lubricant films, in lubricant development and evaluation, and for fundamental research into elastohydrodynamic lubrication and tribofilm formation mechanisms. The film thickness measurement, as implemented on several popular tribometers, provides powerful, non-destructive in-situ mapping of film topography with nanometre-scale height sensitivity. However, the results can be highly sensitive to experimental procedure, machine condition, and image analysis, in some cases reporting unphysical film thickness trends. The prevailing image analysis techniques make it challenging to interrogate these errors, often hiding their multivariate nonlinear behaviour from the user by spatial averaging. Herein, several common ‘silent errors’ in the SLIM measurement, including colour matching to incorrect fringe orders, and colour drift due to the optical properties of the system or film itself, are discussed, with examples. A robust suite of novel a priori and a posteriori methods to address these issues, and to improve the accuracy and reliability of the measurement, are also presented, including a novel, computationally inexpensive circle-finding algorithm for automated image processing. In combination, these methods allow reliable mapping of films up to at least 800 nm in thickness, representing a significant milestone for the utility of SLIM applied to elastohydrodynamic contact.

EHL film geometry↗

DOE OSTI2021

Viscous effects on plasmoid formation from nonlinear resistive tearing growth in a Harris sheet

Abstract In this work, the evolution of a highly unstable m = 1 resistive tearing mode, leading to plasmoid formation in a Harris sheet, is studied in the framework of full MHD model using the Non-Ideal Magnetohydrodynamics with Rotation, Open Discussion simulation. Following the initial nonlinear growth of the primary m = 1 island, the X-point develops into a secondary elongated current sheet that eventually breaks into plasmoids. Two distinctive viscous regimes are found for the plasmoid formation and saturation. In the low viscosity regime (i.e. P r ≲ 1), the plasmoid width increases sharply with viscosity, whereas in the viscosity dominant regime (i.e. P r ≳ 1), the plasmoid size gradually decreases with viscosity. Such a finding quantifies the role of viscosity in modulating the plasmoid formation process through its effects on the plasma flow and the reconnection itself.

Physics↗

DOE OSTI2025

Modeling hydraulic fracture opening and closure with proppant transport and settlement

Hydraulic fracturing is a widely used reservoir stimulation technique for improving fluid circulation in rock formations with extremely low permeability, particularly in enhanced geothermal systems (EGS). To better understand the complex processes involved and improve hydraulic stimulation performance, we have developed ELK (ELectrical fracKing), a MOOSE-based 3D finite element application designed to model the behavior of proppant-fluid mixtures in propagating fractures. ELK integrates both the fluid and proppant components, incorporating particle-driven processes such as gravity settling, particle-particle interactions, and strong density and viscosity contracts, in addition to conventional fluid-driven fracture propagation. In this contribution, we extend ELK to model propped fracture closure, which occurs after the injection phase due to a dramatic drop in the effective stress on the fracture plane. During the shut-in, flowback, and production periods, the fracture width decreases, with the closure behavior depending on proppant concentration.. At low concentrations, closure follows a nonlinear joint law linked to the stiffness of asperities in the fracture walls. While at high concentrations, it is controlled by the properties of packed proppant bed. The extended ELK application is validated against several benchmark examples, including the propagation of an inclined frictional crack, fracture opening and sliding in response to fluid injection, and flowback analysis. We believe that ELK’s enhanced capabilities can serve as a valuable tool for the design and optimization of EGS deployment.

15 - GEOTHERMAL ENERGY↗

DOE OSTI2024

Complexity-calibrated benchmarks for machine learning reveal when prediction algorithms succeed and mislead

Abstract Recurrent neural networks are used to forecast time series in finance, climate, language, and from many other domains. Reservoir computers are a particularly easily trainable form of recurrent neural network. Recently, a “next-generation” reservoir computer was introduced in which the memory trace involves only a finite number of previous symbols. We explore the inherent limitations of finite-past memory traces in this intriguing proposal. A lower bound from Fano’s inequality shows that, on highly non-Markovian processes generated by large probabilistic state machines, next-generation reservoir computers with reasonably long memory traces have an error probability that is at least $$\sim 60\%$$ ∼ 60 % higher than the minimal attainable error probability in predicting the next observation. More generally, it appears that popular recurrent neural networks fall far short of optimally predicting such complex processes. These results highlight the need for a new generation of optimized recurrent neural network architectures. Alongside this finding, we present concentration-of-measure results for randomly-generated but complex processes. One conclusion is that large probabilistic state machines—specifically, large $$\epsilon$$ ϵ -machines—are key to generating challenging and structurally-unbiased stimuli for ground-truthing recurrent neural network architectures.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Ultrasound Characterization of Plastic Bonded Explosives With Varied Binder Content and Density

Advanced explosive formulations and methods of production require rapid characterization for both processes and materials. Here, we present results for ultrasound characterization of a plastic bonded explosive (PBX) formulation as a function of binder weight percent and pressing density using through-transmission sound speed measurements. The explosive in this study was composed of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) combined with the fluoropolymer binder Kynar Ultraflex B. Binder content ranged from 2.6 wt.% to 11.5 wt.% with associated explosive formulation densities from 1.69 g/cc to 1.89 g/cc. We find linear and nonlinear changes in longitudinal and shear sound speeds, peak frequencies, and measured elastic properties for the different PBX formulations. These measured changes can differentiate the tested formulations, highlighting the potential for ultrasound as a nondestructive characterization method for newly-formulated PBXs. We also demonstrate through-transmission measurements using different coupling media and present results from an ultrasound technique called resonant ultrasound spectroscopy (RUS). Both transmission and RUS approaches may allow for rapid and even automated characterization of PBX components suitable for high-throughput manufacturing operations.

Engineering↗

DOE OSTI2020

Partitioned exponential methods for coupled multiphysics systems

Multiphysics problems involving two or more coupled physical phenomena are ubiquitous in science and engineering. This work develops a new partitioned exponential approach for the time integration of multiphysics problems. After a possible semi-discretization in space, the class of problems under consideration is modeled by a system of ordinary differential equations where the right-hand side is a summation of two component functions, each corresponding to a given set of physical processes. The partitioned-exponential methods proposed herein evolve each component of the system via an exponential integrator, and information between partitions is exchanged via coupling terms. Here, the traditional approach to constructing exponential methods, based on the variation-of-constants formula, is not directly applicable to partitioned systems. Rather, our approach to developing new partitioned-exponential families is based on a general-structure additive formulation of the schemes. Two method formulations are considered, one based on a linear-nonlinear splitting of the right hand component functions, and another based on approximate Jacobians. The paper develops classical (non-stiff) order conditions theory for partitioned exponential schemes based on particular families of T-trees and B-series theory. Several practical methods of third order are constructed that extend the Rosenbrock-type and EPIRK families of exponential integrators. Several implementation optimizations specific to the application of these methods to reaction-diffusion systems are also discussed. Numerical experiments reveal that the new partitioned-exponential methods can perform better than traditional unpartitioned exponential methods on some problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Theory of stimulated and spontaneous axion scattering.

We present a theory for nonlinear, resonant excitation of dynamical axions by counterpropagating electromagnetic waves in materials that break both 𝒫 and 𝒯 symmetries. We show that dynamical axions can mediate an exponential growth in the amplitude of the lower frequency (Stokes) beam. We also discuss spontaneous generation of a counterpropagating Stokes mode, enabled by resonant amplification of quantum and thermal fluctuations in the presence of a single pump laser. Remarkably, the amplification can be orders of magnitude larger than that obtained via stimulated Brillouin and Raman scattering processes, and can be modulated with the application of external magnetic fields, making stimulated axion scattering promising for optoelectronics applications.

Smith, M.↗

DOE OSTI2022

Effectively using multifidelity optimization for wind turbine design

Abstract. Wind turbines are complex multidisciplinary systems that are challenging to design because of the tightly coupled interactions between different subsystems. Computational modeling attempts to resolve these couplings so we can efficiently explore new wind turbine systems early in the design process. Low-fidelity models are computationally efficient but make assumptions and simplifications that limit the accuracy of design studies, whereas high-fidelity models capture more of the actual physics but with increased computational cost. This paper details the use of multifidelity methods for optimizing wind turbine designs by using information from both low- and high-fidelity models to find an optimal solution at reduced cost. Specifically, a trust-region approach is used with a novel corrective function built from a nonlinear surrogate model. We find that for a diverse set of design problems – with examples given in rotor blade geometry design, wind turbine controller design, and wind power plant layout optimization – the multifidelity method finds the optimal design using 38 %–58 % of the computational cost of the high-fidelity-only optimization. The success of the multifidelity method in disparate applications suggests that it could be more broadly applied to other wind energy or otherwise generic applications.

17 WIND ENERGY↗

DOE OSTI2025

A Parametric, Data-Driven, Non-Intrusive Reduced-Order Model Framework for Crystal Plasticity Simulations of Voids

The influence of the internal structure at micrometer length scales on the deformation of polycrystalline materials can be effectively captured using crystal plasticity finite element methods (CPFEM). However, the complexity and nonlinearity of the deformation equations CPFEM solves demand significant computational power and resources to achieve accurate predictions, limiting its broader application. To address this challenge, we have identified a reduced-order representation of the complex data in order to establish a computationally efficient reduced-order models (ROM) and drastically reduce the computational expense of CPFEM. Specifically, in this work, we developed a parametric, data-driven, and non-intrusive ROM framework for CPFEM using proper orthogonal decomposition (POD) and sparse variational Gaussian process (SVGP) regression for single-crystal microstructures under tensile loading conditions. The developed protocol enables one to compress field into a latent/low-dimensional space described by principal component analysis (PCA) via the singular value decomposition (SVD) algorithm. As a result, the high-dimensional data are reduced to a significantly smaller amount of dimensions with POD bases and POD coefficients. Furthermore, we deployed an ensemble of SVGPs—extended from the classical Gaussian process (GP) regression for scalability and handling big data—in a massively parallel manner to train and predict latent POD coefficients using known POD bases from a set of previously obtained simulations results. Lastly, using the predicted POD coefficients, we reconstructed the full-field results and showed reasonable agreement compared with the true values obtained from running CPFEM. The developed framework is validated with a set of CPFEM simulations of a single embedded void in single-crystal aluminum alloy. While the framework is broadly applicable, this work specifically focuses on single-crystal microstructures, a single load case (e.g., tensile), and a specific void geometry (spherical).

Anisotropy↗

DOE OSTI2023

Soliton formation and topology manipulation of coupled spins via ultrafast re-magnetization

The major goal of the project was to explore the properties of magnetic order driven far out of equilibrium by optical excitations. These include fundamental questions related to the interplay of magnetic, structural and electronic degrees of freedom in materials that are optically excited. These problems break down in short term quenching of the magnetization as a results of the energy of the optical pulse studying how the energy of the optical pulse and the angular momentum of the magnetization flow between the different degrees of freedom. This is followed by the longer time re-emergence of the magnetism as the system cools. Within this goal we explored the formation of solitons, both with and without non-trivial topology, via rapid re-magnetization processes after optical-driven ultrafast demagnetization. Guided by theory, we predict that turbulence and modulational instabilities will drive the formation of solitons, including dispersive shock waves, magnon droplets, and skyrmions. In the case of topological defect formation, i.e. skyrmion generation, related processes have long been predicted based on general principles of symmetry-breaking phase transitions; the density of topological defects in a long-range-ordered phase can be controlled by varying the quench rate through the second-order phase transition, i.e. the Kibble-Zurek (KZ) mechanism. We, more broadly, had the goal to study a range of other emergent magnetic behaviors after optical excitation. While great attention has been given to ultrafast demagnetization, less is known of the subsequent spin dynamics and coupling to the lattice. What is increasingly appreciated is that ultrafast demagnetization leads to spin currents that can carry angular momentum from the rapidly demagnetized sample. These spin currents appear in many ways but fundamentally controls demagnetization, drives interactions between regions of the material, and control the spin resulting structure. They can be probed many ways including THz emission which, in turn, gives insight into the demagnetization processes. Simultaneously if there is magneto-elastic coupling ultrafast demagnetization can drive structural excitations that are not expected from thermal energy added to the lattice. We explored a range of complex phenomena that arise from photo-excitation of magnetic systems with the coupling of electronic, magnetic and structural degrees of freedom. We have made a number of fundamental discoveries that are reflected in our publications list with additional work still being prepared for publication. Highlights of this work include (i) ultra-efficient, nonlinear THz surface acoustics (ii) spin-current-mediated rapid magnon localization and coalescence, (iii) spin-wave soliton formation in ferromagnetic FePt nanoparticles, (iv) dynamic phonon coupling in elemental antiferromagnetic Cr, (v) THz emission from Co/Pt bilayers and FeRh/Pt bilayers, (vi) theoretical investigation of spin hydrodynamics, solitons and shock waves, and (vii) ultrafast perturbation of magnetic domains by optical pumping.

36 MATERIALS SCIENCE↗