Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear approximation”

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 OSTI2016

Scalable Optimization Methods for Distribution Networks With High PV Integration

Here this paper proposes a suite of algorithms to determine the active- and reactive-power setpoints for photovoltaic (PV) inverters in distribution networks. The objective is to optimize the operation of the distribution feeder according to a variety of performance objectives and ensure voltage regulation. In general, these algorithms take a form of the widely studied ac optimal power flow (OPF) problem. For the envisioned application domain, nonlinear power-flow constraints render pertinent OPF problems nonconvex and computationally intensive for large systems. To address these concerns, we formulate a quadratic constrained quadratic program (QCQP) by leveraging a linear approximation of the algebraic power-flow equations. Furthermore, simplification from QCQP to a linearly constrained quadratic program is provided under certain conditions. The merits of the proposed approach are demonstrated with simulation results that utilize realistic PV-generation and load-profile data for illustrative distribution-system test feeders.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2018

Neural Networks as Surrogates of Nonlinear High-Dimensional Parameter-to-Prediction Maps

We present a preliminary investigation of the use of Multi-Layer Perceptrons (MLP) and Recurrent Neural Networks (RNNs) as surrogates of parameter-to-prediction maps of computational expensive dynamical models. In particular, we target the approximation of Quantities of Interest (QoIs) derived from the solution of a Partial Differential Equations (PDEs) at different time instants. In order to limit the scope of our study while targeting a relevant application, we focus on the problem of computing variations in the ice sheets mass (our QoI), which is a proxy for global mean sea-level changes. We present a number of neural network formulations and compare their performance with that of Polynomial Chaos Expansions (PCE) constructed on the same data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Velocity space evolution of a minority energetic electron population undergoing the anomalous Doppler instability

The kinetic evolution in velocity space of a minority suprathermal electron population that is undergoing the anomalous Doppler instability (ADI) is investigated using the results from fully nonlinear numerical simulations that self-consistently evolve particles and fields in a plasma. Electron trajectories in phase space during different stages of the ADI are captured, and are analysed in relation to the characteristics of the excited electric fields and of the overall distribution of particles. For some electrons, trapping and mirroring effects are observed during the saturation phase. A relationship between the second order moments of the perpendicular electron distribution function and time is established, and is used to investigate the range of applicability of analytical approximations drawn from classical theory, that involve a quasilinear wave-driven diffusion operator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Neural Networks as Surrogates of Nonlinear High-Dimensional Parameter-to-Prediction Maps.

We present a preliminary investigation of the use of Multi-Layer Perceptrons (MLP) and Recurrent Neural Networks (RNNs) as surrogates of parameter-to-prediction maps of com- putational expensive dynamical models. In particular, we target the approximation of Quan- tities of Interest (QoIs) derived from the solution of a Partial Differential Equations (PDEs) at different time instants. In order to limit the scope of our study while targeting a rele- vant application, we focus on the problem of computing variations in the ice sheets mass (our QoI), which is a proxy for global mean sea-level changes. We present a number of neural network formulations and compare their performance with that of Polynomial Chaos Expansions (PCE) constructed on the same data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Relativistic Electron Increase During Chorus Wave Activities on the 6–8 March 2016 Geomagnetic Storm

There was a geomagnetic storm on 6–8 March 2016, in which Van Allen Probes A and B separated by ~2.5 h measured increase of relativistic electrons with energies approximately several hundred keV to 1 MeV. Simultaneously, chorus waves were measured by both Van Allen Probes and Magnetospheric Multiscale (MMS) mission. Some of the chorus elements were rising tones, possibly due to nonlinear effects. In this paper, these measurements are compared with a nonlinear theory of chorus waves incorporating the inhomogeneity ratio and the field equation. From this theory, a chorus wave profile in time and one-dimensional space is simulated. Test particle calculations are then performed in order to examine the energization rate of electrons. Some electrons are accelerated, although more electrons are decelerated. Finally, the measured time scale of the electron increase is inferred to be consistent with this nonlinear theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Convergence of sum-up rounding schemes for cloaking problems governed by the Helmholtz equation

In this work, we consider the problem of designing a cloak for waves described by the Helmholtz equation from an integer programming point of view. The problem can be modeled as a PDE-constrained optimization problem with integer-valued control inputs that are distributed in the computational domain. A first-discretize-then-optimize approach results in a large-scale mixed-integer nonlinear program that is in general intractable because of the large number of integer variables that arise from the discretization of the domain. Instead, we propose an efficient algorithm that is able to approximate the local infima of the underlying nonconvex infinite-dimensional problem arbitrarily close without the need to solve the discretized finite-dimensional integer programs to optimality. We optimize only the continuous relaxations of the approximations for local minima and then apply the sum-up rounding methodology to obtain integer-valued controls. If the solutions of the discretized continuous relaxations converge to a local minimizer of the continuous relaxation, then the resulting discrete-valued control sequence converges weakly \(^*\) in \(L^\infty\) to the same local minimizer. These approximation properties follow under suitable refinements of the involved discretization grids. Our results use familiar concepts arising from the analytical properties of the underlying PDE and complement previous results, derived from a topology optimization point of view.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Homogenizing atomic dynamics by fractional differential equations

In this paper, we propose two ways to construct fractional differential equations (FDE) for approximating atomic chain dynamics. Taking harmonic chain as an example, we add a power function of fractional order to Taylor expansion of the dispersion relation, and determine the parameters by matching two selected wave numbers. This approximate function leads to an FDE after considering both directions for wave propagation. As an alternative, we consider the symbol of the force term, and approximate it by a similar function. It also induces an FDE. Both approaches produce excellent agreement with the harmonic chain dynamics. The accuracy may be improved by optimizing the selected wave numbers, or starting with higher order Taylor expansions. When resolved in the lattice constant, the resulting FDE's faithfully reproduce the lattice dynamics. When resolved in a coarse grid instead, they systematically generate homogenized algorithms. Numerical tests are performed to verify the proposed approaches. Moreover, FDE's are also constructed for diatomic chain and anharmonic lattice, to illustrate the generality of the proposed approaches. - Highlights: • Match the dispersion relation of linear chain by power function of fractional order, added on low order Taylor expansion. • Design FDE by approximating dispersion relation, or the symbol for the force term. • Solve the FDE by algorithm of spectral accuracy. • Obtain faithful and systematical homogenization when resolved at coarser grid. • Applies to general lattices: complex, nonlinear and multidimensional.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

The effect of Coulomb interactions on nonlinear thermovoltage and thermocurrent in quantum dots

In the present work, we theoretically study the nonlinear regime of charge transport through a quantum dot coupled to the source and drain reservoirs. The investigation is carried out using a nonequilibrium Green’s function formalism beyond the Hartree-Fock approximation. Employed approximations for the relevant Green’s functions allow to trace a transition from Coulomb blockade regime to Kondo regime in the thermoelectric transport. Effects arising when electrons move in response to thermal gradient applied across the system are discussed, including experimentally observed thermovoltage zeros.

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

Dimits shift, avalanche-like bursts, and solitary propagating structures in the two-field flux-balanced Hasegawa–Wakatani model for plasma edge turbulence

We show that the recently introduced two-field flux-balanced Hasegawa–Wakatani (BHW) model captures the key features of drift-wave turbulent transport mediated by zonal flows observed in more complete and accurate gyrokinetic simulations, such as the existence of a nonlinear upshift of the threshold for drift wave turbulence driven transport, often called the Dimits shift, as well as non-local transport with avalanche bursts and solitary propagating structures. Because of the approximations made in the BHW model, these observations are made for the particle flux instead of the heat flux more commonly studied in ion temperature gradient (ITG) driven turbulence in fluid or gyrokinetic codes. Many of these features are not seen in other Hasegawa–Wakatani models, which confirm the critical role of the electron dynamics parallel to the magnetic field lines. To address questions regarding the role of boundary conditions on the drift-wave zonal flow dynamics, we apply our model to both a channel domain geometry and the more typical doubly periodic geometry. We only observe strong soliton-like solutions in the particle flux for the channel geometry, in the vicinity of the boundaries, where strong velocity shear and density gradients are generated, which are absent in the doubly periodic simulations. Changing the aspect ratio of the simulation domain also has a significant effect. In domains which are elongated in the radial direction, more complex multi-scale dynamics takes place, with multiple zonal jets interacting with each other, and large scale avalanches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

CFL Optimized Forward–Backward Runge–Kutta Schemes for the Shallow-Water Equations

Abstract We present the formulation and optimization of a Runge–Kutta-type time-stepping scheme for solving the shallow-water equations, aimed at substantially increasing the effective allowable time step over that of comparable methods. This scheme, called FB-RK(3,2), uses weighted forward–backward averaging of thickness data to advance the momentum equation. The weights for this averaging are chosen with an optimization process that employs a von Neumann–type analysis, ensuring that the weights maximize the admittable Courant number. Through a simplified local truncation error analysis and numerical experiments, we show that the method is at least second-order in time for any choice of weights and exhibits low dispersion and dissipation errors for well-resolved waves. Further, we show that an optimized FB-RK(3,2) can take time steps up to 2.8 times as large as a popular three-stage, third-order strong stability-preserving Runge–Kutta method in a quasi-linear test case. In fully nonlinear shallow-water test cases relevant to oceanic and atmospheric flows, FB-RK(3,2) outperforms SSPRK3 in admittable time step by factors roughly between 1.6 and 2.2, making the scheme approximately twice as computationally efficient with little to no effect on solution quality. Significance Statement The purpose of this work is to develop and optimize time-stepping schemes for models relevant to oceanic and atmospheric flows. Specifically, for the shallow-water equations we optimize for schemes that can take time steps as large as possible while retaining solution quality. We find that our optimized schemes can take time steps between 1.6 and 2.2 times larger than schemes that cost the same number of floating point operations, translating directly to a corresponding speedup. Our ultimate goal is to use these schemes in climate-scale simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Quantitative Nonlinear Optical Polarimetry with High Spatial Resolution

Nonlinear optical microscopy such as in the optical second-harmonic generation (SHG) modality has become a popular tool today for probing materials in the physical and biological sciences. While imaging and spectroscopy are widely used in the microscopy mode, nonlinear polarimetry, which can shed light on materials’ symmetry and microstructure, is relatively underdeveloped. This is partly because quantitative analytical modeling of the optical SHG response for anisotropic crystals and films largely assumes low-numerical aperture (NA) focusing of light, where the plane-wave approximation is sufficient. Tight focusing provides unique benefits in revealing out-of-plane polarization responses, which cannot be detected by near-plane-wave illumination at normal incidence. Here, we outline a method for quantitatively analyzing SHG polarimetry measurements obtained under high-NA focusing within a microscope geometry. Experiments and simulations of a variety of standard samples, from single crystals to thin films, are in good agreement, including measured and simulated spatial SHG maps of ferroelectric domains. A solution to the inverse problem is demonstrated, where the spatial distribution of an SHG tensor with unknown tensor coefficient magnitudes is determined by experimentally measured polarimetry. The ability to extract the out-of-plane component of the nonlinear polarization in normal incidence is demonstrated, which can be valuable for high-resolution polarimetry of 2D materials, thin films, heterostructures, and uniaxial crystals with a strong out-of-plane response.

36 MATERIALS SCIENCE↗
DOE OSTI2015

An adhesive contact mechanics formulation based on atomistically induced surface traction

In this work, we have developed a novel multiscale computational contact formulation based on the generalized Derjuguin approximation for continua that are characterized by atomistically enriched constitutive relations in order to study macroscopic interaction between arbitrarily shaped deformable continua. The proposed adhesive contact formulation makes use of the microscopic interaction forces between individual particles in the interacting bodies. In particular, the double-layer volume integral describing the contact interaction (energy, force vector, matrix) is converted into a double-layer surface integral through a mathematically consistent approach that employs the divergence theorem and a special partitioning technique. The proposed contact model is formulated in the nonlinear continuum mechanics framework and implemented using the standard finite element method. With no large penalty constant, the stiffness matrix of the system will in general be well-conditioned, which is of great significance for quasi-static analysis. Three numerical examples are presented to illustrate the capability of the proposed method. Results indicate that with the same mesh configuration, the finite element computation based on the surface integral approach is faster and more accurate than the volume integral based approach. In addition, the proposed approach is energy preserving even in a very long dynamic simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Uncertainty propagation in feed-forward neural network models

We develop new uncertainty propagation methods for feed-forward neural network architectures with leaky ReLU activation functions subject to random perturbations in the input vectors. In particular, we derive analytical expressions for the probability density function (PDF) of the neural network output and its statistical moments as a function of the input uncertainty and the parameters of the network, i.e., weights and biases. A key finding is that an appropriate linearization of the leaky ReLU activation function yields accurate statistical results even for large perturbations in the input vectors. This can be attributed to the way information propagates through the network. We also propose new analytically tractable Gaussian copula surrogate models to approximate the full joint PDF of the neural network output. To validate our theoretical results, we conduct Monte Carlo simulations and a thorough error analysis on a multi-layer neural network representing a nonlinear integro-differential operator between two polynomial function spaces. Our findings demonstrate excellent agreement between the theoretical predictions and Monte Carlo simulations.

MLP networks↗
DOE OSTI2015

Three-point bending analysis of doubly clamped silicon nanowire beams; Young's modulus, initial stress, and crystal orientation

A non-linear model is introduced describing the force-deflection relation of doubly clamped beams, including initial stress. Several approximations for the exact model are developed and compared, revealing the importance of considering the initial stress during 3-point bending measurements analysis. A novel approximation is found to be better than others, and both the exact model and this approximation are in perfect agreement with finite element simulations. A brief experimental example of silicon nanowires is presented in which the Young's modulus, the initial stress, and the crystallographic growth orientation are extracted by 3-point bending analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Combined magnetic and kinetic control of advanced tokamak steady state scenarios based on semi-empirical modelling

Our paper shows that semi-empirical data-driven models based on a two-time-scale approximation for the magnetic and kinetic control of advanced tokamak (AT) scenarios can be advantageously identified from simulated rather than real data, and used for control design. The method is applied to the combined control of the safety factor profile, q(x), and normalized pressure parameter, beta(N), using DIII-D parameters and actuators (on-axis co-current neutral beam injection (NBI) power, off-axis co-current NBI power, electron cyclotron current drive power, and ohmic coil). The approximate plasma response model was identified from simulated open-loop data obtained using a rapidly converging plasma transport code, METIS, which includes an MHD equilibrium and current diffusion solver, and combines plasma transport nonlinearity with 0D scaling laws and 1.5D ordinary differential equations. The paper discusses the results of closed-loop METIS simulations, using the near-optimal ARTAEMIS control algorithm (MoreauDet al 2013 Nucl. Fusion 53 063020) for steady stateAT operation. With feedforward plus feedback control, the steady state target q-profile and beta(N) are satisfactorily tracked with a time scale of about 10 s, despite large disturbances applied to the feedforward powers and plasma parameters. We also show the robustness of the control algorithm with respect to disturbances of the H&CD actuators and of plasma parameters such as the H-factor, plasma density and effective charge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Black hole formation from axion stars

The classical equations of motion for an axion with potential V (φ)= m {sub a} {sup 2} f {sub a} {sup 2} [1−cos (φ/ f {sub a} )] possess quasi-stable, localized, oscillating solutions, which we refer to as ''axion stars''. We study, for the first time, collapse of axion stars numerically using the full non-linear Einstein equations of general relativity and the full non-perturbative cosine potential. We map regions on an ''axion star stability diagram', parameterized by the initial ADM mass, M {sub ADM}, and axion decay constant, f {sub a} . We identify three regions of the parameter space: i) long-lived oscillating axion star solutions, with a base frequency, m {sub a} , modulated by self-interactions, ii) collapse to a BH and iii) complete dispersal due to gravitational cooling and interactions. We locate the boundaries of these three regions and an approximate ''triple point' ( M {sub TP}, f {sub TP}) ∼ (2.4 M {sub pl}{sup 2}/ m {sub a} ,0.3 M {sub pl}). For f {sub a} below the triple point BH formation proceeds during winding (in the complex U(1) picture) of the axion field near the dispersal phase. This could prevent astrophysical BH formation from axion stars with f {sub a} || M {sub pl}. For larger f {sub a} ∼> f {sub TP}, BH formation occurs through the stable branch and we estimate the mass ratio of the BH to the stable state at the phase boundary to be O(1) within numerical uncertainty. We discuss the observational relevance of our findings for axion stars as BH seeds, which are supermassive in the case of ultralight axions. For the QCD axion, the typical BH mass formed from axion star collapse is M {sub BH} ∼ 3.4 ( f {sub a} /0.6 M {sub pl}){sup 1.2} M {sub ⊙}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Bulk Photovoltaic Effect Driven by Collective Excitations in a Correlated Insulator

We investigate the bulk photovoltaic effect, which rectifies light into electric current, in a collective quantum state with correlation driven electronic ferroelectricity. We show via explicit real-time dynamical calculations that the effect of the applied electric field on the electronic order parameter leads to a strong enhancement of the bulk photovoltaic effect relative to the values obtained in a conventional insulator. The enhancements include both resonant enhancements at sub-band-gap frequencies, arising from excitation of optically active collective modes, and broadband enhancements arising from nonresonant deformations of the electronic order. The deformable electronic order parameter produces an injection current contribution to the bulk photovoltaic effect that is entirely absent in a rigid-band approximation to a time-reversal symmetric material. Furthermore, our findings establish that correlation effects can lead to the bulk photovoltaic effect and demonstrate that the collective behavior of ordered states can yield large nonlinear optical responses.

14 SOLAR ENERGY↗