Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear reduction”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Multilevel Techniques for Compression and Reduction of Scientific Data---The Unstructured Case

Previous work on multilevel techniques for compression and reduction of scientific data is extended to the case of data given on unstructured meshes in two and three dimensions. The centerpiece of the work is a decomposition algorithm which is shown to be optimal, in terms of both storage and operational complexity, applicable to unstructured grids in both two and three dimensions, and which implicitly gives a Riesz basis that can be exploited to reduce the data while maintaining rigorous bounds on the loss incurred. The flexibility of the approach is illustrated by applications to potential flow around an airfoil and the effect of compression on quantities of interest relevant to airfoil design; compression of computational simulation of a nonlinear reaction-diffusion system with special attention given to the problem of time series reduction; and, data from a simulation of magnetically confined plasma in a fusion reactor reduced so as to preserve the electric field computed from the data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Predicting nonequilibrium Green’s function dynamics and photoemission spectra via nonlinear integral operator learning

Understanding the dynamics of nonequilibrium quantum many-body systems is an important research topic in a wide range of fields across condensed matter physics, quantum optics, and high-energy physics. However, numerical studies of large-scale nonequilibrium phenomena in realistic materials face serious challenges due to intrinsic high-dimensionality of quantum many-body problems and the absence of time-invariance. The nonequilibrium properties of many-body systems can be described by the dynamics of the correlator, or the Green's function of the system, whose time evolution is given by a high-dimensional system of integro-differential equations, known as the Kadanoff–Baym equations (KBEs). The time-convolution term in KBEs, which needs to be recalculated at each time step, makes it difficult to perform long-time numerical simulation. In this paper, we develop an operator-learning framework based on recurrent neural networks (RNNs) to address this challenge. We utilize RNNs to learn the nonlinear mapping between Green's functions and convolution integrals in KBEs. By using the learned operators as a surrogate model in the KBE solver, we obtain a general machine-learning scheme for predicting the dynamics of nonequilibrium Green's functions. Besides significant savings per each time step, the new methodology reduces the temporal computational complexity from $O(N_t^3)$ to $O(N_t)$ where N t is the number of steps taken in a simulation, thereby making it possible to study large many-body problems which are currently infeasible with conventional KBE solvers. Through various numerical examples, we demonstrate the effectiveness of the operator-learning based approach in providing accurate predictions of physical observables such as the reduced density matrix and time-resolved photoemission spectra. Moreover, our framework exhibits clear numerical convergence and can be easily parallelized, thereby facilitating many possible further developments and applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Multifrequency stimulated Raman scattering of light in a calcite single crystal

Using pulses of multifrequency stimulated Raman scattering (SRS) in calcite, nonlinear photoluminescence in a stilbene molecular crystal is excited. When multifrequency SRS is excited by the radiation of a Nd{sup 3+} : YAG laser with a wavelength of 1064 nm, eleven anti-Stokes components are observed in the visible spectrum with an average frequency shift of 1086 cm{sup −1} between them. Using the second harmonic of the Nd {sup 3+} : YAG laser radiation allowed three anti-Stokes and four Stokes components to be recorded. The large spectral range of the frequency comb facilitated effective reduction in the duration and increase in the intensity of the radiation pulse of multifrequency SRS. (nonlinear optical phenomena)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Robust Output Feedback Control Design for Inertia Emulation by Wind Turbine Generators

Wind generation has gained widespread use as a renewable energy source. Most wind turbines and other renewables connected to the grid through converters result in a reduction in the natural inertial response to grid frequency changes. The doubly-fed induction generator (DFIG) can be controlled to compensate for this reduction and, in fact, provide faster response than traditional synchronous machines. This paper proposes to design observer based output feedback linear quadratic regulator (LQR) and H control laws to realize the inertia emulation function and deliver fast frequency support. Furthermore, the aim is to track the reference speed by a diesel synchronous generator (DSG) in order to reach the desired inertia. The control signal is computed based on a reduced order model using the balanced truncation technique. A comparison with selective modal analysis (SMA) and balanced truncation model reduction techniques is presented. Comprehensive results show the effective emulation of synthetic inertia by implementing the control laws on a nonlinear three- phase diesel-wind system. The proposed technique is analyzed for different short circuit ratio (SCR) scenarios.

17 WIND ENERGY↗
DOE OSTI2022

Preconditioning for Hyper-reduction in Reduced Order Models

Many projection-based reduced order models (pROM) that utilize the governing equation and data to accelerate physical simulations can be applied to nonlinear dynamical systems. To avoid full order model (FOM) scale update for each time step, hyperreduction techniques are developed to sample high dimensional nonlinear terms. Our study aims to investigate if preconditioning the least-squares problem used for the nonlinear approximation can improve the robustness of the condition number of the problem while achieving high accuracy. In our study, we use the row-normalization matrix motivated by the Christoffel function as the preconditioner and solve the corresponding weighted least-squares problem. Numerical results for Lagrangian hydrodynamics examples are analyzed to explore how the preconditioner works compared to existing hyper-reduction techniques.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Molecular Structure and Thermodynamics of CO 2 and Water Adsorption on Mica

The adsorption of CO 2 and water on clay surfaces plays a key role in applications, such as gas storage in saline aquifers and depleted hydrocarbon reservoirs, but is not yet fully understood. Here, we study the adsorption of CO 2 and water vapor using Grand Canonical Monte Carlo and molecular dynamics simulations. At a bulk pressure of 100 bar, pure CO 2 adsorbs strongly on mica and forms extensive layers next to it. CO 2 adsorption is lowered substantially if introducing water vapor above mica and is largely eliminated when the relative humidity (RH) approaches about 60%. When pure water vapor is introduced above a mica surface, a subnanometer thick liquid water film develops on it to form apparent liquid–solid and liquid–vapor interfaces simultaneously. Using the identification of truly interfacial molecules (ITIM) analysis, we delineate how individual water layers develop in this film as RH increases. We highlight that the water film is spatially heterogeneous and the true liquid–vapor interface emerges only at an RH of 60–80%. Introducing 100 bar of CO 2 into the water vapor above the mica surface modulates water adsorption nonlinearly: at RH = 0.01%, the water adsorption is reduced by ∼30%; as RH increases, the reduction is weakened, and eventually, enhancement of water adsorption by about 7% occurs at RH = 90%. These variations are attributed to the interplay of film thinning by high-pressure CO 2 , competition of mica surface sites by CO 2 molecules, and energetic and entropic stabilization of interfacial water by CO 2 molecules.

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

Microturbulence Suppression by Alfvén Eigenmodes in the DIII-D Tokamak

Mitigation and suppression of low-𝑘 turbulence are observed during the nonlinear evolution of toroidicity-induced Alfvén eigenmodes (TAEs) in DIII-D experiments. Turbulence mitigation begins when the dominant TAE starts to depart from the typical shear Alfvén wave polarization. TAE activity then evolves into a new state, characterized by more discrete coherent modes, increased total amplitude, and more localized radial structures. During this transition, a narrow shear flow layer forms, driven by an enhanced Reynolds stress force, with a shearing rate that exceeds the local turbulence decorrelation rate, leading to a full turbulence suppression. Furthermore, these observations indicate the imbalance between Reynolds and Maxwell stress forces during the nonlinear evolution of the TAE and their important roles in shear flow generation and turbulence reduction.

Alfvén waves↗
DOE OSTI2020

Truncated Nonlinear Interferometry for Quantum-Enhanced Atomic Force Microscopy

Nonlinear interferometers that replace beam splitters in Mach-Zehnder interferometers with nonlinear amplifiers for quantum-enhanced phase measurements have drawn increasing interest in recent years, but practical quantum sensors based on nonlinear interferometry remain an outstanding challenge. Here, we demonstrate the first practical application of nonlinear interferometry by measuring the displacement of an atomic force microscope microcantilever with quantum noise reduction of up to 3 dB below the standard quantum limit, corresponding to a quantum-enhanced measurement of beam displacement of 1.7 fm/√Hz. Further, we minimize photon backaction noise while taking advantage of quantum noise reduction by transducing the cantilever displacement signal with a weak squeezed state while using dual homodyne detection with a higher power local oscillator. This approach may enable quantum-enhanced broadband, high-speed scanning probe microscopy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A Tailored Convolutional Neural Network for Nonlinear Manifold Learning of Computational Physics Data Using Unstructured Spatial Discretizations

In this work, we propose a nonlinear manifold learning technique based on deep convolutional autoencoders that is appropriate for model order reduction of physical systems in complex geometries. Convolutional neural networks have proven to be highly advantageous for compressing data arising from systems demonstrating a slow-decaying Kolmogorov n-width. However, these networks are restricted to data on structured meshes. Unstructured meshes are often required for performing analyses of real systems with complex geometry. Our custom graph convolution operators based on the available differential operators for a given spatial discretization effectively extend the application space of deep convolutional autoencoders to systems with arbitrarily complex geometry that are typically discretized using unstructured meshes. We propose sets of convolution operators based on the spatial derivative operators for the underlying spatial discretization, making the method particularly well suited to data arising from the solution of partial differential equations. We demonstrate the method using examples from heat transfer and fluid mechanics and show better than an order of magnitude improvement in accuracy over linear methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Reduction or Enhancement of Stellarator Turbulence by Impurities

A systematic study of the impact of impurities on the turbulent heat fluxes is presented for the stellarator Wendelstein 7-X (W7-X) and, for comparison, the Large Helical Device and ITER. By means of nonlinear multispecies gyrokinetic simulations, it is shown that impurities, depending on the sign of their density gradient, can significantly enhance or reduce turbulent ion heat losses. For the relevant scenario of turbulence reduction, an optimal impurity concentration that minimizes the ion heat diffusivity emerges as a universal feature. This result demonstrates the potential of impurities for controlling turbulence and accessing enhanced confinement regimes in fusion plasmas and, in particular, in W7-X. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electrochemical CO2 Reduction over Metal-/Nitrogen-Doped Graphene Single-Atom Catalysts Modeled Using the Grand-Canonical Density Functional Theory

Renewably driven, electrochemical conversion of carbon dioxide into value-added products is expected to be a critical tool in global decarbonization. However, theoretical studies based on the computational hydrogen electrode largely ignore the nonlinear effects of the applied potential on the calculated results, leading to inaccurate predictions of catalytic behavior or mechanistic pathways. Here, we use grand canonical density functional theory (GC-DFT) to model electrochemical CO2 reduction (CO2R) over metal- and nitrogen-doped graphene catalysts (MNCs) and explicitly include the effects of the applied potential. We used GC-DFT to compute the CO2 to CO reaction intermediate energies at -0.3, -0.7, and -1.2 VSHE catalyzed by MNCs each doped with 1 of the 10 3d block metals coordinated by four pyridinic nitrogen atoms. Our results predict that Sc-, Ti-, Co-, Cu-, and Zn-N4Cs effectively catalyze CO2R at moderate to large reducing potentials (-0.7 to -1.2 VSHE). ZnN4C is a particularly promising electrocatalyst for CO2R to CO both at low and moderate applied potentials based on our thermodynamic analysis. Our findings also explain the observed pH independence of CO production over FeN4C and predict that the rate-determining step of CO2R over FeN4C is not *CO2- formation but rather *CO desorption. Additionally, the GC-DFT-computed density of states analysis illustrates how the electronic states of MNCs and adsorbates change non-uniformly with applied potential, resulting in a significantly increased *CO2- stability relative to other intermediates and demonstrating that the formation of the adsorbed *CO2- anion is critical to CO2R activation. This work demonstrates how GC-DFT paves the way for physically realistic and accurate theoretical simulations of reacting electrochemical systems.

CO2 reduction↗
DOE OSTI2023

Constraint energy minimizing generalized multiscale finite element method for multi-continuum Richards equations

In fluid flow simulation, the multi-continuum model is a useful strategy. When the heterogeneity and contrast of coefficients are high, the system becomes multiscale, and some kinds of reduced order methods are demanded. Combining these techniques with nonlinearity, we will consider in this paper a dual-continuum model which is generalized as a multi-continuum model for a coupled system of nonlinear Richards equations as unsaturated flows, in complex heterogeneous fractured porous media; and we will solve it by a novel multiscale approach utilizing the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM). In particular, such a nonlinear system will be discretized in time and then linearized by Picard iteration (whose global convergence is proved theoretically). Subsequently, we tackle the resulting linearized equations by the CEM-GMsFEM and obtain proper offline multiscale basis functions to span the multiscale space (which contains the pressure solution). More specifically, we first introduce two new sources of samples, and the GMsFEM is used over each coarse block to build local auxiliary multiscale basis functions via solving local spectral problems, that are crucial for detecting high-contrast channels. Second, per oversampled coarse region, local multiscale basis functions are created through the CEM as constrainedly minimizing an energy functional. Various numerical tests for our approach reveal that the error converges with the coarse-grid size and that only few oversampling layers as well as basis functions are needed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Cross phases of temperature-gradient-driven turbulence as a model basis for I -mode particle transport

As a model for understanding the type of transport behavior characteristic of the tokamak I mode, cross-phase physics for particle-transport is studied analytically for turbulence dominated by either ion-temperature-gradient (ITG) or electron-temperature-gradient (ETG) instability. I mode is a transport-barrier regime of reduced thermal transport but essentially unaffected particle transport. It is assumed that ITG turbulence applies to the baseline L mode, ETG to I mode, and that E × B flow shear is stronger in I mode, lowering all fluxes. In ITG turbulence, particle transport is governed by trapped electrons. Sensitivity to collisions produces the well-known temperature-gradient-driven pinch that offsets density-gradient-driven outward diffusion, weakening particle transport in L mode. In ETG turbulence, nonadiabatic ions are collisionless. Nonzero transport requires an ion spectrum feature whose magnetic-drift resonance supplies the necessary cross phase. If frequencies of order the ion diamagnetic drift frequency dominate the ion part of the spectrum, as would occur with weakly unstable ITG turbulence, all components of the particle transport are outward and can offset flow-shear-induced flux reductions to produce a flux that is similar to the ITG L-mode particle flux. Nonlinear frequencies are potentially relevant and discussed in relation to I mode.

Physics↗
DOE OSTI2025

Heterogeneous Mixtures of Dictionary Functions to Approximate Subspace Invariance in Koopman Operators: Why Deep Koopman Operators Work

Abstract Koopman operators model nonlinear dynamics as a linear dynamic system acting on a nonlinear function as the state. This nonstandard state is often called a Koopman observable and is usually approximated numerically by a superposition of functions drawn from a dictionary . In a widely used algorithm, extended dynamic mode decomposition (EDMD), the dictionary functions are drawn from a fixed class of functions. Deep learning combined with EDMD has been used to learn novel dictionary functions in an algorithm called deep dynamic mode decomposition (deepDMD). The learned representation both (1) accurately models and (2) scales well with the dimension of the original nonlinear system. In this paper, we analyze the learned dictionaries from deepDMD and explore the theoretical basis for their strong performance. We explore State-Inclusive Logistic Lifting (SILL) dictionary functions to approximate Koopman observables. Error analysis of these dictionary functions show they satisfy a property of subspace approximation, which we define as uniform finite approximate closure. Typically, a Koopman dictionary’s nonlinear functions are homogeneous. In this paper, we discover that structured mixing of heterogeneous dictionary functions drawn from different classes of nonlinear functions achieve the same accuracy and dimensional scaling as the deep-learning-based deepDMD algorithm Yeung et al. ( In: 2019 American Control Conference (ACC), 2019). We specifically show this by building a heterogeneous dictionary comprised of SILL functions and conjunctive radial basis functions (RBFs). This mixed dictionary achieves similar accuracy and dimensional scaling to deepDMD with an order of magnitude reduction in parameters, while maintaining geometric interpretability. These results strengthen the viability of dictionary-based Koopman models to solving high-dimensional nonlinear learning problems.

Johnson, Charles A.↗
DOE OSTI2020

On differentiable local bounds preserving stabilization for Euler equations

This work presents the design of nonlinear stabilization techniques for the finite element discretization of Euler equations in both steady and transient form. Implicit time integration is used in the case of the transient form. A differentiable local bounds preserving method has been developed, which combines a Rusanov artificial diffusion operator and a differentiable shock detector. Nonlinear stabilization schemes are usually stiff and highly nonlinear. This issue is mitigated by the differentiability properties of the proposed method. Moreover, in order to further improve the nonlinear convergence, we also propose a continuation method for a subset of the stabilization parameters. The resulting method has been successfully applied to steady and transient problems with complex shock patterns. Numerical experiments show that it is able to provide sharp and well resolved shocks. Furthermore, the importance of the differentiability is assessed by comparing the new scheme with its non-differentiable counterpart. Numerical experiments suggest that, for up to moderate nonlinear tolerances, the method exhibits improved robustness and nonlinear convergence behavior for steady problems. Additionally, in the case of transient problem, we also observe a reduction in the computational cost.

42 ENGINEERING↗
DOE OSTI2025

Plasmonic Ag nanocomposite phosphate glasses produced via γ-ray irradiation as reduction route

This paper reports on the impact of γ-ray irradiation (10, 100 kGy) on melt-quenched Ag + -doped phosphate glass and the effects of subsequent thermal processing leading to the production of plasmonic Ag nanocomposites. The γ-irradiated glasses were characterized alongside the pristine by differential scanning calorimetry (DSC), Raman spectroscopy, electron paramagnetic resonance (EPR) spectroscopy, optical absorption, and photoluminescence (PL) spectroscopy. DSC characterization showed consistent glass transition temperatures (T g ) before and after γ-irradiation whereas the crystallization temperatures tended to decrease with increasing γ-ray dose. However, a lack of alteration of the glass network structure was supported by Raman spectroscopy. Room temperature EPR spectra clearly showed the formation of phosphorus oxygen hole center (POHC) defects in the undoped host, in addition to another doublet likely associated with a P 3 defect. The presence of paramagnetic silver species encompassing Ag 2+ and 107/109 Ag 0 atoms was also indicated in the silver-activated glass together with POHC defects. Optical absorption spectra were also consistent with the presence of various radiation-induced centers. Further analyzing the glass absorption edge via Tauc plots suggested the formation of electron center (EC) defects in γ-irradiated samples wherein the silver-doped glass exhibited decreasing band gap energies with increasing γ-ray dose. The PL characterization showed the silver-related radio-PL effect was induced exhibiting broad band emission with two maxima around 500 and 625 nm stemming from various molecular Ag$^{x+}_{n}$ clusters. Emission decay analyses revealed that the longer wavelength emission exhibited slower decay. The highest radiation dose of 100 kGy however resulted in weaker emission and faster decay kinetics attributed to energy transfer between the luminescent silver species and POHC defects. Finally subjecting the γ-irradiated Ag-doped glasses to heat treatment near the T g at 490 °C led to the development of the surface plasmon resonance of Ag nanoparticles (NPs) and the vanishing of the Ag$^{x+}_{n}$ clusters luminescence. In conclusion, the presence of the matrix-related EC defects was deemed accountable for the thermally induced reduction and consequent precipitation of Ag NPs making the plasmonic glasses attractive for photonic applications such as nonlinear optics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Chiral Anomaly in Interacting Condensed Matter Systems

The chiral anomaly is a fundamental quantum mechanical phenomenon which is of great importance to both particle physics and condensed matter physics alike. In the context of QED, it manifests as the breaking of chiral symmetry in the presence of electromagnetic fields. It is also known that anomalous chiral symmetry breaking can occur through interactions alone, as is the case for interacting one-dimensional systems. In this Letter, we investigate the interplay between these two modes of anomalous chiral symmetry breaking in the context of interacting Weyl semimetals. Using Fujikawa’s path integral method, we show that the chiral charge continuity equation is modified by the presence of interactions which can be viewed as including the effect of the electric and magnetic fields generated by the interacting quantum matter. This can be understood further using dimensional reduction and a Luttinger liquid description of the lowest Landau level. These effects manifest themselves in the nonlinear response of the system. In particular, we find an interaction-dependent density response due to a change in the magnetic field as well as a contribution to the nonequilibrium and inhomogeneous anomalous Hall response while preserving its equilibrium value.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Sampling Low-Dimensional Markovian Dynamics for Preasymptotically Recovering Reduced Models from Data with Operator Inference

This work introduces a method for learning low-dimensional models from data of high-dimensional black-box dynamical systems. The novelty is that the learned models are exactly the reduced models that are traditionally constructed with classical projection-based model reduction techniques. Thus, the proposed approach learns models that are guaranteed to have the well-studied properties of reduced models known from model reduction, without requiring full knowledge of the governing equations and without requiring the operators of the high-dimensional systems. The key ingredient is a new data sampling scheme to obtain re-projected trajectories of high-dimensional systems that correspond to Markovian dynamics in low-dimensional subspaces. The exact recovery of reduced models from these re-projected trajectories is guaranteed pre-asymptotically under certain conditions for finite amounts of data and for a large class of systems with polynomial nonlinear terms. Numerical results demonstrate that the low-dimensional models learned with the proposed approach match reduced models from traditional model reduction up to numerical errors in practice. In conclusion, the numerical results further indicate that low-dimensional models fitted to re-projected trajectories are predictive even in situations where models fitted to trajectories without re-projection are inaccurate and unstable.

97 MATHEMATICS AND COMPUTING↗