Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear function”

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 37 records · Page 2

Deep Learning Explicit Differentiable Predictive Control Laws for Buildings

We present a differentiable predictive control (DPC) methodology for learning constrained control laws for unknown nonlinear systems. DPC poses an approximate solution to multiparametric programming problems emerging from explicit nonlinear model predictive control (MPC). Contrary to approximate MPC, DPC does not require supervision by an expert controller. Instead, a system dynamics model is learned from a small dataset of recorded observations of the perturbed system's dynamics and the control law is optimized offline by interaction with the learned system model. The DPC method is based on two sequential steps, i) system identification using a constrained neural state-space model, and ii) optimization of an explicit control law parametrized by another neural network in closed-loop simulation with the identified neural state-space model. The combination of a differentiable closed-loop system and penalty methods for constraint handling of system outputs and inputs allows us to optimize the control law's parameters directly by backpropagating economic MPC loss through the learned system model. By incorporating domain knowledge and leveraging established techniques from optimal control, our method leverages deep neural networks as nonlinear function approximators for system identification and control while avoiding concomitant costs of intractably large datasets, and computationally expensive over-parametrized models. The scalability, data efficiency, and constrained optimal control capability of the proposed DPC method are demonstrated in simulation using a multi-zone building emulator.

Drgona, Jan↗

Pore-resolved investigation of turbulent open channel flow over a randomly packed permeable sediment bed

Pore-resolved direct numerical simulations are performed to investigate the interactions between streamflow turbulence and groundwater flow through a randomly packed porous sediment bed for three permeability Reynolds numbers, Re K = 2.56 , 5.17 and 8.94, representative of natural stream or river systems. Time–space averaging is used to quantify the Reynolds stress, form-induced stress, mean flow and shear penetration depths, and mixing length at the sediment–water interface (SWI). Here, the mean flow and shear penetration depths increase with Re K and are found to be nonlinear functions of non-dimensional permeability. The peaks and significant values of the Reynolds stresses, form-induced stresses, and pressure variations are shown to occur in the top layer of the bed, which is also confirmed by conducting simulations of just the top layer as roughness elements over an impermeable wall. The probability distribution functions (p.d.f.s) of normalized local bed stress are found to collapse for all Reynolds numbers, and their root-mean-square fluctuations are assumed to follow logarithmic correlations. The fluctuations in local bed stress and resultant drag and lift forces on sediment grains are mainly a result of the top layer; their p.d.f.s are symmetric with heavy tails, and can be well represented by a non-Gaussian model fit. The bed stress statistics and the pressure data at the SWI potentially can be used in providing better boundary conditions in modelling of incipient motion and reach-scale transport in the hyporheic zone.

turbulence simulation↗

An analytic theory for the degree of Arctic Amplification

Arctic Amplification (AA), the amplified surface warming in the Arctic relative to the global mean, is a robust and impactful feature of climate change. While the basic physical picture of AA has been depicted, a clear understanding of how the degree of AA is determined has not been established. Here, by deciphering the intricate role of atmospheric heat transport (AHT), we build a two-box energy-balance model of AA and derive that the degree of AA is a simple nonlinear function of the Arctic and global feedbacks, the meridional heterogeneity in radiative forcing, and the partial sensitivities of AHT to global mean warming and meridional warming gradient. The formula captures the varying degree of AA in individual climate models and attributes the variation to specific physical factors. It further conveys a concise picture of how essential physics mutually determine the degree of AA and limit the range within 1.5~3.5. Our results articulate AHT as both forcing and feedback to AA, highlighting its partial sensitivities instead of total change as the key parameters for understanding AA. We also find that the effect of lapse rate feedback, a widely-recognized major contributor to AA, is fully offset by the effect of water vapor feedback.

54 ENVIRONMENTAL SCIENCES↗

Data-based filtered dissipation rate modelling for multi-modal turbulent combustion: evaluating a priori model generalizability

Manifold-based models offer a computationally efficient alternative to directly transporting the thermochemical state in computational simulations of turbulent reacting flows, projecting the high-dimensional thermochemical state-space onto a low-dimensional manifold. Recent efforts have yielded a manifold-based model applicable to multi-modal combustion, enabling reconstruction of the thermochemical state from solutions to two-dimensional manifold equations in mixture fraction and generalized progress variable that are parameterised by three scalar dissipation rates. In coarse-grained simulations such as Large Eddy Simulation (LES), closure of the multi-modal manifold equations and subfilter variances/covariance requires closure of three filtered scalar dissipation rates. Here, the present work adopts a data-based approach, providing closure for the three filtered scalar dissipation rates via deep neural networks (DNNs). High-fidelity datasets corresponding to an autoigniting n-dodecane jet flame and a bluff body swirl-stabilized confined lifted spray flame of two aviation fuels (Jet-A and C1) with different ignition propensities are leveraged to generate training data that spans a diverse range of thermodynamic conditions and combustion modes, including low- and high-temperature ignition regimes in addition to premixed and nonpremixed behaviour. A final DNN model is trained to enforce inherent physical constraints by learning nonlinear functional transformations of the three filtered scalar dissipation rates. The generalizability of this constrained DNN model is demonstrated a priori via conditional statistics evaluated on the lifted spray flame with C1–a configuration that had not been included in the training data. Excellent DNN agreement with conditional DNS statistics is observed, and integrated gradients are computed to identify the most sensitive input variables. The similarity of the marginal PDFs of the most informative input variables and outputs across configurations are quantified via the Wasserstein metric, demonstrating that data-based models may successfully generalize to unseen parametric conditions so long as the most informative input variables share similar distributions across training and testing datasets.

Data-based modelling↗

Adaptive machine learning for time-varying systems: low dimensional latent space tuning

Machine learning (ML) tools such as encoder-decoder convolutional neural networks (CNN) can represent incredibly complex nonlinear functions which map between combinations of images and scalars. For example, CNNs can be used to map combinations of accelerator parameters and images which are 2D projections of the 6D phase space distributions of charged particle beams as they are transported between various particle accelerator locations. Despite their strengths, applying ML to time-varying systems, or systems with shifting distributions, is an open problem, especially for large systems for which collecting new data for re-training is impractical or interrupts operations. Particle accelerators are one example of large time-varying systems for which collecting detailed training data requires lengthy dedicated beam measurements which may no longer be available during regular operations. We present a novel method of adaptive ML for time-varying systems. Our approach is to map very high (N ≈ 100k) dimensional inputs (a combination of scalar parameters and images) into the low dimensional (N ≈ 2) latent space at the output of the encoder section of an encoder-decoder CNN. We then actively tune the low dimensional latent space-based representation of complex system dynamics by the addition of an adaptively tuned feedback vector directly before the decoder sections builds back up to our image-based high-dimensional phase space density representations. This method allows us to learn correlations within and to quickly tune the characteristics of incredibly large parameter space systems and to track their evolution in real time based on feedback without massive new data sets for re-training. We demonstrate that our method can accurately predict and track the phase space of charged particle beams at various locations in a particle accelerator by adaptively adjusting in real-time while the unknown input beam distribution of the accelerator is changing in shape, charge, and offset and while the RF system of the accelerator itself is also changing in an unpredictable way. For FACET-II we demonstrate that such an approach has the potential to use transverse deflecting cavity and energy spread spectrum beam measurements to accurately predict 2D projections of the 6D phase space of the electron beam at the plasma wakefield acceleration interaction point where such diagnostics are unavailable.

47 OTHER INSTRUMENTATION↗

Photon Frequency Conversion in High-Q Superconducting Resonators: Axion Electrodynamics, QED, and Nonlinear Meissner Radiation

High-Q superconducting resonators have been proposed and developed as detectors of light-by-light scattering mediated by the hypothesized axion or virtual electron–positron pairs in quantum electrodynamics: the Euler–Heisenberg (EH) interaction. Photon frequency and mode conversion is central to the scheme for detecting such rare events. Superconducting resonators are nonlinear devices. The Meissner screening currents that confine the electromagnetic fields to the vacuum region of a superconducting RF cavity are nonlinear functions of the EM field at the vacuum–superconducting interface, and as a result can generate source currents and frequency conversion of microwave photons in the cavity. In this report we consider photon frequency and mode conversion in superconducting resonators with high quality factors from Meissner currents in single- and dual-cavity setups proposed for axion and QED searches based on light-by-light scattering. In a single cavity with two pump modes, photon frequency conversion by the Meissner screening current dominates photon generation by the EH interaction for cavities with Q≲1012⁠. The Meissner currents also generate background photons that limit the operation of the resonator for axion detection in three-mode, single-cavity setups. We also consider the leakage of photons from pump modes into the signal mode for both axion- and EH-mediated light-by-light scattering. Photon frequency conversion by the EH interaction can compete with Meissner and leakage radiation in ultra-high-Q cavities that are beyond the current state of the art. Meissner radiation and leakage backgrounds can be suppressed in dual-cavity setups with appropriate choices for pump and spectator modes, as well as the single-cavity setup proposed for heterodyne detection of galactic axion dark matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Integrability technique for fluid flow induced deformation of a boundary hair

The deformation of a dense carpet of hair due to Stokes flow in a channel can be described by a nonlinear integrodifferential equation for the shape of a single hair, which possesses several solutions for a given choice of parameters. Although it was posed in a previous study and it bears a resemblance to the pendulum problem from mechanics, this equation has not been analytically solved until now. Despite the presence on an integral with a nonlinear functional dependence on the dependent variable, the system is integrable. We compare the analytically obtained solution to a finite-difference numerical approach, identify the physically realizable solution branch, and briefly study the solution structure through a conserved energylike quantity. Time-dependent fluid-structure interactions are a rich and complex subject to investigate, and we argue that the solution discussed herein can be used as a basis for understanding these systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electron Dynamics Within a MITL Containing a Load

In this paper, we derive the vacuum electric fields within specific cylindrically symmetric magnetically insulated transmission lines (MITLs) in the limit of an infinite speed of light for an arbitrary time-dependent current. We focus our attention on two types of MITLs: the radial MITL and a spherically curved MITL. We then simulate the motion of charged particles, such as electrons, present in these MITLs due to the vacuum fields. In general, the motion of charged particles due to the vacuum fields is highly nonlinear since the fields are nonlinear functions of spatial coordinates and depend on an arbitrary time-dependent current drive. Using guiding center theory, however, one can describe the gross particle kinetics using a combination of $\textbf {E} \times \textbf {B}$ and $\nabla B$ drifts. In addition, we compare our approximate inner MITL field models and particle kinetics with those from a fully electromagnetic simulation code. We find that the agreement between the approximate model and the electromagnetic simulations is excellent.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improving Trustworthiness of Data-Driven Power Grid Contingency Analysis With Bayesian Residual Graph Neural Networks

The evolving energy landscape requires novel tools to efficiently perform contingency analysis and reliability assessment of power grids, potentially in real-time. The high computational cost of traditional power flow solvers limits their applicability in practice. Machine learning (ML) surrogates such as deep neural networks (NNs) accelerate power flow solvers computations, enabling high-order contingency analysis and real-time decision-making by learning highly nonlinear functions and integrating grid topology via graph architectures. However, (graph) NNs lack predictive power away from training data and do not provide predictive confidence estimates. Here, we present a Bayesian residual graph NN that integrates knowledge from low-fidelity data via residual training and embeds granular quantification of uncertainties, improving trustworthiness critical for high-consequence decision-making. Applying Bayesian concepts to NNs is challenging due to the high-dimensionality of both the parameter space, complicating derivation of a meaningful prior, and the output space in large grid systems, requiring enhanced techniques to assess the predicted high-dimensional uncertainties. Our contributions include: (1) Deriving a prior for fully connected and graph NNs that leverages low-fidelity data to guide mean predictions and appropriately control prior predictive uncertainty. (2) Integrating this prior within an ensembling with anchoring scheme for efficient approximate posterior inference. (3) Deriving enhanced metrics to assess accuracy of both the mean and uncertainty predictions in high dimensions, appropriately accounting for correlations propagated through graph layers. The resulting Bayesian residual graph NN is tested on a contingency analysis task for 14-bus and 118-bus grids.

24 - POWER TRANSMISSION AND DISTRIBUTION↗

sf-24-103

pyzag is a library for efficiently training generic models defined with a recursive nonlinear function. We developed it primarily to provide an efficient means to train NEML2 models it can be used to train other, more general types of models.

Messner, Mark [Argonne National Laboratory (ANL), ↗

Lagrangian and Eulerian Supersaturation Statistics in Turbulent Cloudy Rayleigh–Bénard Convection: Applications for LES Subgrid Modeling

Abstract Turbulent fluctuations of scalar and velocity fields are critical for cloud microphysical processes, e.g., droplet activation and size distribution evolution, and can therefore influence cloud radiative forcing and precipitation formation. Lagrangian and Eulerian water vapor, temperature, and supersaturation statistics are investigated in direct numerical simulations (DNS) of turbulent Rayleigh–Bénard convection in the Pi Convection Cloud Chamber to provide a foundation for parameterizing subgrid-scale fluctuations in atmospheric models. A subgrid model for water vapor and temperature variances and covariance and supersaturation variance is proposed, valid for both clear and cloudy conditions. Evaluation of phase change contributions through an a priori test using DNS data shows good performance of the model. Supersaturation is a nonlinear function of temperature and water vapor, and relative external fluxes of water vapor and heat (e.g., during entrainment-mixing and phase change) influence turbulent supersaturation fluctuations. Although supersaturation has autocorrelation and structure functions similar to the independent scalars (temperature and water vapor), the autocorrelation time scale of supersaturation differs. Relative scalar fluxes in DNS without cloud make supersaturation PDFs less skewed than the adiabatic case, where they are highly negatively skewed. However, droplet condensation changes the PDF shape response: it becomes positively skewed for the adiabatic case and negatively skewed when the sidewall relative fluxes are large. Condensation also increases correlations between water vapor and temperature in the presence of relative scalar fluxes but decreases correlations for the adiabatic case. These changes in correlation suppress supersaturation variability for the nonadiabatic cases and increase it for the adiabatic case. Implications of this work for subgrid microphysics modeling using a Lagrangian stochastic scheme are also discussed.

54 ENVIRONMENTAL SCIENCES↗

Probing the Wind Component of Radio Emission in Luminous High-redshift Quasars

We discuss a probe of the contribution of wind-related shocks to the radio emission in otherwise radio-quiet quasars. Given (1) the nonlinear correlation between UV and X-ray luminosity in quasars, (2) that such a correlation leads to higher likelihood of radiation-line-driven winds in more luminous quasars, and (3) that luminous quasars are more abundant at high redshift, deep radio observations of high-redshift quasars are needed to probe potential contributions from accretion disk winds. We target a sample of 50 z ≃ 1.65 color-selected quasars that span the range of expected accretion disk wind properties as traced by broad C iv emission. 3 GHz observations with the Very Large Array to an rms of ≈10 μJy beam -1 probe to star formation rates of ~400 M ⊙ yr -1 , leading to 22 detections. Supplementing these pointed observations are survey data of 388 sources from the LOFAR Two-meter Sky Survey Data Release 1 that reach comparable depth (for a typical radio spectral index), where 123 sources are detected. These combined observations reveal a radio detection fraction that is a nonlinear function of C iv emission-line properties and suggest that the data may require multiple origins of radio emission in radio-quiet quasars. We find evidence for radio emission from weak jets or coronae in radio-quiet quasars with low Eddington ratios, with either (or both) star formation and accretion disk winds playing an important role in optically luminous quasars and correlated with increasing Eddington ratio. Additional pointed radio observations are needed to fully establish the nature of radio emission in radio-quiet quasars.

79 ASTRONOMY AND ASTROPHYSICS↗

Green’s Function Perspective on the Nonlinear Density Response of Quantum Many-Body Systems

We derive equations of motion for higher order density response functions using the theory of thermodynamic Green’s functions. We also derive expressions for the higher order generalized dielectric functions and polarization functions. Moreover, we relate higher order response functions and higher order collision integrals within the Martin–Schwinger hierarchy. We expect our results to be highly relevant to the study of a variety of quantum many-body systems such as matter under extreme temperatures, densities, and pressures.

Density Functional Theory↗

Gauge invariance and Ward identities in nonlinear response theory

We present a formal analysis of nonlinear response functions in terms of correlation functions in real- and imaginary-time domains. In particular, we show that causal nonlinear response functions, expressed in terms of nested commutators in real time, can be obtained from the analytic continuation of time-ordered response functions, which are more easily amenable to diagrammatic calculation. This generalizes the well-known result of linear response theory. We then use gauge invariance arguments to derive exact relations between second-order response functions in density and current channels. These identities, which are non-perturbative in the strength of inter-particle interactions, allow us to establish exact connections between nonlinear optics calculations done in different electromagnetic gauges.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Topological enhancement of nonlinear transport in unconventional point-node semimetals

The topological singularity of the Bloch states close to the Fermi level significantly enhances nonlinear electric responses in topological semimetals. Here, we systematically characterize this enhancement for a large class of topological nodal-point fermions, including those with linear, linear-quadratic, and quadratic dispersions. Specifically, we determine the leading power-law dependence of the nonlinear response functions on the chemical potential μ defined relative to the nodal point. We identify two characteristics that qualitatively improve nonlinear transports compared to those of conventional Dirac and Weyl fermions. First, the type-II (overtilted) spectrum leads to the log μ enhancement of nonlinear response functions having zero scaling dimension with respect to μ, which is not seen in a type-I (moderately or not tilted) spectrum. Second, the anisotropic linear-quadratic dispersion increases the power of small-μ divergence for the nonlinear response tensors along the linearly dispersing direction. Furthermore, our work reveals new experimental signatures of unconventional nodal points in topological semimetals as well as provides a guiding principle for giant nonlinear electric responses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The impacts of convex piecewise linear cost formulations on AC optimal power flow

Despite strong connections through shared application areas, research efforts on power market optimization (e.g., unit commitment) and power network optimization (e.g., optimal power flow) remain largely independent. A notable illustration of this is the treatment of power generation cost functions, where nonlinear network optimization has largely used polynomial representations and market optimization has adopted piecewise linear encodings. This work combines state-of-the-art results from both lines of research to understand the best mathematical formulations of the nonlinear AC optimal power flow problem with piecewise linear generation cost functions. An extensive numerical analysis of non-convex models, linear approximations, and convex relaxations across fifty-four realistic test cases illustrates that nonlinear optimization methods are surprisingly sensitive to the mathematical formulation of piecewise linear functions. The results indicate that a poor formulation choice can slow down algorithm performance by a factor of ten, increasing the runtime from seconds to minutes. Furthermore, these results provide valuable insights into the best formulations of nonlinear optimal power flow problems with piecewise linear cost functions, an important step towards building a new generation of energy markets that incorporate the nonlinear AC power flow model.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Nonlinear susceptibilities for weakly turbulent magnetized plasma: Electrostatic approximation

The plasma weak turbulence theory is a perturbative nonlinear theory, which has been proven to be quite valid in a number of applications. However, the standard weak turbulence theory found in the literature is fully developed for highly idealized unmagnetized plasmas. As many plasmas found in nature and laboratory are immersed in a background static magnetic field, it is necessary to extend the existing discussions to include the effects of ambient magnetic field. Such a task is quite formidable, however, which has prevented fundamental and significant progresses in the subject matter. The central difficulty lies in the formulation of the complete nonlinear response functions for magnetized plasmas. The present paper derives the nonlinear susceptibilities for weakly turbulent magnetized plasmas up to the third order nonlinearity, but in doing so, a substantial reduction in mathematical complexity is achieved by the use of Bessel function addition theorem (or sum rule). The present paper also constructs the weak turbulence wave kinetic equation in a formal sense. For the sake of simplicity, however, the present paper assumes the electrostatic interaction among plasma particles. Fully electromagnetic generalization is a subject of a subsequent paper.

Physics↗

Statistical Learning for Nonlinear Model Reduction from Local Simulations of Stochastic and Particle- and Agent-Based Systems

Stochastic physical systems across the sciences that have very high-dimensional state spaces, with a large number of fast degrees of freedom that force direct simulators to proceed by integration steps that are orders of magnitude smaller than events of interests (e.g., particle collisions). Examples range from molecular motion to dynamics of large populations of cells. A grand challenge in the simulation and understanding of such systems is the systematic construction of accurate, interpretable, reduced models, enabling faster simulations, revealing fundamental properties of the dynamics, and predicting phenomena of interest that the original simulator could not reached with sufficient accuracy or within a given computational budget. In this projected we developed novel statistical estimation/machine learning techniques for analyzing and building empirical reduced models for important families of high-dimensional stochastic systems, in particular: - we developed techniques for estimating interaction kernels in interacting particle- and agent-based systems, which are ubiquitous in Physics, Biology and many other sciences, given observed trajectories of the system; - we developed techniques for nonlinear model reduction for high-dimensional stochastic systems that have a small number of unknown, nonlinear slow variables, and a large number of fast modes, that are possibly of large magnitude, given observed short trajectories of the system in the form of bursts of trajectories from different initial conditions; - we developed novel techniques for estimating linear dynamical systems on graphs when both the dynamics and the underlying graph are unknown, and we have a sparse set of space-time observations; - we considered the problem of estimating an unknown nonlinear observation function of a standard process (e.g. Brownian motion), so that we can recognized if an observed dynamics is "just" a nonlinear version of a known dynamics; we also developed benchmarks for learning algorithms aimed at learning and classifying diffusion processes.

97 MATHEMATICS AND COMPUTING↗