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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Streamlined Loads Analysis of Floating Wind Turbines With Fiber Rope Mooring Lines

This paper presents an approach for more accurate yet relatively streamlined accounting for the nonlinear characteristics of synthetic fiber rope mooring lines for floating wind turbines. First, we select a minimal set of parameters that can efficiently approximate a mooring line material's quasi-static and dynamic stiffness characteristics. We also develop a set of baseline coefficients for different rope materials based on published product information and research papers. We then expand a quasi-static mooring model to include dynamic mooring line stiffness terms in a way that allows the nonlinear stiffness behavior of fiber ropes to be considered even in quasi-static analyses. For dynamic analysis, we have updated the model MoorDyn, coupled with OpenFAST, to work with the new dynamic mooring line stiffness terms. This includes the implementation of a new viscoelastic approach that allows the tension-strain relation of each mooring line segment to vary between two stiffness values, depending on the strain rate. After presenting the formulation of the approaches for modeling synthetic ropes, they are demonstrated on a floating wind turbine mooring system with taut polyester rope mooring lines. The results are compared with those of other approaches of similar fidelity, including the static-dynamic method and separate simulations with static and dynamic stiffness values. Comparing the results shows the ability of the new method to match the results of the previous methods in a more streamlined manner.

elasticity↗
DOE OSTI2024

Streamlined Loads Analysis of Floating Wind Turbines With Fiber Rope Mooring Lines: Preprint

This paper presents an approach for more accurate yet relatively streamlined accounting for the nonlinear characteristics of synthetic fiber rope mooring lines for floating wind turbines. First, we select a minimal set of parameters that can efficiently approximate a mooring line material's quasi-static and dynamic stiffness characteristics. We also develop a set of baseline coefficients for different rope materials based on published product information and research papers. We then expand a quasi-static mooring model to include dynamic mooring line stiffness terms in a way that allows the nonlinear stiffness behavior of fiber ropes to be considered even in quasi-static analyses. For dynamic analysis, we have updated the model MoorDyn, coupled with OpenFAST, to work with the new dynamic mooring line stiffness terms. This includes the implementation of a new viscoelastic approach that allows the tension-strain relation of each mooring line segment to vary between two stiffness values depending on the strain rate. After presenting the formulation of the approaches for modeling synthetic ropes, they are demonstrated on a floating wind turbine mooring system with taut polyester rope mooring lines. The results are compared with those of other approaches of similar fidelity including the static-dynamic method and separate simulations with static and dynamic stiffness values. Comparing the results shows the ability of the new method to match the results of the previous methods in a more streamlined manner.

elasticity↗
DOE OSTI2023

Optofluidic memory and self-induced nonlinear optical phase change for reservoir computing in silicon photonics

Nanophotonics allows to employ light-matter interaction to induce nonlinear optical effects and realize non-conventional memory and computation capabilities, however to date, light-liquid interaction was not considered as a potential mechanism to achieve computation on a nanoscale. Here, we experimentally demonstrate self-induced phase change effect which relies on the coupling between geometric changes of thin liquid film to optical properties of photonic waveguide modes, and then employ it for neuromorphic computing. In our optofluidic silicon photonics system we utilize thermocapillary-based deformation of thin liquid film capable to induce nonlinear effect which is more than one order of magnitude higher compared to the more traditional heat-based thermo-optical effect, and allowing operation as a nonlinear actuator and memory element, both residing at the same compact spatial region. The resulting dynamics allows to implement Reservoir Computing at spatial region which is approximately five orders of magnitude smaller compared to state-of-the-art experimental liquid-based systems.

42 ENGINEERING↗
DOE OSTI2023

Recent Developments in Security-Constrained AC Optimal Power Flow: Overview of Challenge 1 in the ARPA-E Grid Optimization Competition

In “Recent Developments in Security-Constrained AC Optimal Power Flow: Overview of Challenge 1 in the ARPA-E Grid Optimization Competition,” we review the state of the art in practical algorithms for scheduling power-systems operations in the short term and the results of the recent competition organized by the U.S. Advanced Research Projects Agency–Energy. We explain the mixed-integer nonlinear formulation used in the competition for nonspecialists in electrical engineering, the context and organization of the competition, and the performance of competitors. We find that the collective approaches and results of competitors provide support for efforts to move nonlinear optimization techniques into industrial applications, as they have proven to be a robust and efficient alternative to current linear approximation techniques.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

A composite neural network that learns from multi-fidelity data: Application to function approximation and inverse PDE problems

Presently the training of neural networks relies on data of comparable accuracy but in real applications only a very small set of high-fidelity data is available while inexpensive lower fidelity data may be plentiful. We propose a new composite neural network (NN) that can be trained based on multi-fidelity data. It is comprised of three NNs, with the first NN trained using the low-fidelity data and coupled to two high-fidelity NNs, one with activation functions and another one without, in order to discover and exploit nonlinear and linear correlations, respectively, between the low-fidelity and the high-fidelity data. We first demonstrate the accuracy of the new multi-fidelity NN for approximating some standard benchmark functions but also a 20-dimensional function that is not easy to approximate with other methods, e.g. Gaussian process regression. Subsequently, we extend the recently developed physics-informed neural networks (PINNs) to be trained with multi-fidelity data sets (MPINNs). MPINNs contain four fully-connected neural networks, where the first one approximates the low-fidelity data, while the second and third construct the correlation between the low- and high-fidelity data and produce the multi-fidelity approximation, which is then used in the last NN that encodes the partial differential equations (PDEs). Specifically, by decomposing the correlation into a linear and nonlinear part, the present model is capable of learning both the linear and complex nonlinear correlations between the low- and high-fidelity data adaptively. By training the MPINNs, we can: (1) obtain the correlation between the low- and high-fidelity data, (2) infer the quantities of interest based on a few scattered data, and (3) identify the unknown parameters in the PDEs. In particular, we employ the MPINNs to learn the hydraulic conductivity field for unsaturated flows as well as the reactive models for reactive transport. The results demonstrate that MPINNs can achieve relatively high accuracy based on a very small set of high-fidelity data. Despite the relatively low dimension and limited number of fidelities (two-fidelity levels) for the benchmark problems in the present study, the proposed model can be readily extended to very high-dimensional regression and classification problems involving multi-fidelity data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Accelerating the solution of the S N equations with highly anisotropic scattering using the Fokker-Planck approximation

The discrete ordinates method can model forward-peaked transport problems accurately. However, convergence of discrete ordinates solution can become arbitrarily slow upon use of standard iterative procedures like source iteration. Standard zeroth and first moment-based acceleration methods like nonlinear diffusion acceleration and diffusion synthetic acceleration are ineffective in accelerating such problems because these methods do not correct higher order Legendre-moments of angular flux. In this paper, we explore the idea of using the Fokker-Planck approximation as a preconditioner to accelerate the convergence of solution of forward-peaked transport problems. Finally, we observe that such preconditioning can successfully reduce the iteration count by up to four orders of magnitude, and solver runtime by up to three orders of magnitude when compared to diffusion synthetic acceleration.

42 ENGINEERING↗
DOE OSTI2025

Slender-body approach for computing second-order wave loads in the frequency domain

This work presents a slender-body approach to evaluate the second-order wave loads acting on a floating structure in the frequency domain. The approach is in the same spirit as the common use of Morison’s equation to approximate the wave loads without solving the radiation/diffraction problem. To do so, we employ Rainey’s equation, which can be seen as an extension of the inertial part of Morison’s equation to include nonlinear effects. We introduce modifications to Rainey’s formulation in order to evaluate wave kinematics at the mean body position instead of the original approach of considering instantaneous displacements. We also propose a simple approximation to partially account for wave scattering effects on the second-order loads based on the analytical solution of a surface-piercing bottom-mounted vertical circular cylinder. Though limited to structures composed of cylinders, this slender-body approach is orders of magnitude faster than computing second-order wave coefficients with a radiation/diffraction code. We implemented this approach for difference-frequency (slow drift) loads in an open-source frequency-domain floating wind turbine model. We present comparisons against results obtained with radiation/diffraction theory for three reference floating wind turbine designs: the OC3-Hywind spar, the OC4-DeepCwind semisubmersible, and the VolturnUS-S semisubmersible. In general, the results show that the proposed slender-body approach with the correction to approximate wave scattering effects provides useful estimations of the difference-frequency wave loads and the resulting motions of the floater.

17 WIND ENERGY↗
DOE OSTI2025

Causal machine learning uncovers conditions for convective intensification driven by organic and sulfate aerosols

Aerosols are often hypothesized to invigorate deep convective clouds (DCCs), but observational evidence remains limited and inconclusive. Clarifying this hypothesis is critical for regions vulnerable to thunderstorms and flooding, particularly highly polluted coastal cities. Leveraging a novel causal discovery–inference pipeline and high-resolution observations near Houston, TX, we identify multiple causal pathways among aerosols (mostly organic and sulfate), DCCs, and meteorological factors. However, a direct causal link from aerosols to DCCs is found to be uncommon, occurring in less than 35% of analyzed scenarios, and is characterized by strong conditionality and nonlinearity. When aerosol impacts on DCCs do occur, they can be substantial, enhancing DCC core heights by approximately 1.7 km, with 92% of this effect concentrated in warmer-phase cloud regions. Notably, the presence of sea breezes and the inclusion of all measured aerosol particles each enhance DCCs in over 95% of aerosol-sensitive cases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Instability of an electron-plasma shear layer in an externally imposed strain flow

The E x B shear instability of a two-dimensional (2D) filament (i.e., a thin, rectangular strip perpendicular to the magnetic field) of magnetized pure electron plasma is investigated experimentally in the presence of an externally imposed strain flow. Data are acquired using a specialized Penning–Malmberg trap in which strain flows can be applied in 2D by biasing segmented electrodes surrounding the plasma. The E x B drift dynamics are well-described by the Drift-Poisson equations, which are isomorphic to the 2D Euler equations describing ideal fluids. Thus, the experimental results correspond to the Rayleigh instability of a shear layer in a 2D ideal fluid, where the electron density is analogous to the fluid vorticity. Shear layers are prepared by stretching initially axisymmetric electron vortices using a strong, applied strain flow. The data at early times are in quantitative agreement with a linear model which extends Rayleigh’s work to account for the influence of an external strain flow. In the presence of weak strain, the system approximately maintains a phase relationship that corresponds to an instantaneous Rayleigh eigenmode. The instability develops into the nonlinear regime later in time and at smaller spatial scales as the strain rate is increased. A secondary vortex pairing instability is observed, but it is suppressed when the strain-to-vorticity ratio exceeds roughly 0.025. In this way, vorticity transport perpendicular to the filament is diminished due to the applied strain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Iterative Linearization for Phasor-Defined Optimal Power Dispatch

Optimal power flow (OPF) problems, which dispatch power targets to controllable generating units across a network, must generally account for non-convex constraints on power flow. Furthermore, adapting those problems so as to make them solvable with convex optimization techniques is an area of much academic and operational interest. In this paper, we present a method for solving OPF as a quadratic program by iteratively refining and re-initializing a linearized model of power flow based on the outputs of an associated nonlinear solver. The linear model on which we demonstrate this method is an adapted version of an approximation designed for use with unbalanced distribution networks. As an important benefit, the model allows for the explicit inclusion of nodal voltage phasor values in both the OPF problem's objective and its constraints, which opens the door to the idea of phasor-based control (PBC) design. We show in simulations on the IEEE 13-node test feeder that our method quickly converges to a set of phasor targets that are sufficiently precise for use in operations at the distribution level.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Theoretical study of various nonlinear phenomena in plasma systems and scaling of magneto-inertial-fusion targets

Plasma physics is an exciting field of study with a wide variety of nonlinear processes that come into play. Examples of such processes include the interaction of small-scale turbulence with large-scale plasma structures and the nonlinear saturation of plasma instabilities, for example those of magneto-hydrodynamical nature. During this Truman LDRD project, I studied a collection of nonlinear problems that are of interest to the field of plasma physics. This LDRD report summarizes four main research accomplishments. First, a new statistical model for describing inhomogeneous drift-wave turbulence inter- acting with zonal flows was developed. This new model includes the effects of nonlinear wave-wave collisions, which are expected to change the spectrum of the underlying DW turbulence and therefore the generation of zonal flows. Second, a new mathematical formalism was proposed to systematically apply the non- linear WKB approximation to general field theories, including those often used in fluid dynamics. This formalism represents an interesting tool for studying physical systems that show an explicit scale separation. Third, a weakly nonlinear model was developed to describe the magneto-Rayleigh-Taylor instability. This instability is of paramount importance to understand as it can reduce the performance of magnetic-inertial-fusion (MIF) platforms. The developed models captures the effects of harmonic generation and saturation of the linear growth of the instability. Finally, a framework was proposed for scaling magneto-inertial fusion (MIF) targets to larger pulsed-power drivers. From this framework, a set of scaling rules were derived that conserve the physical regimes of MIF systems when scaling up in peak current. By doing so, deleterious nonlinear processes that affect MIF performance may be kept at bay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Implicit and coupled fluid plasma solver with adaptive Cartesian mesh and its applications to non-equilibrium gas discharges

In this work, we present a new fluid plasma solver with adaptive Cartesian mesh (ACM) based on a full-Newton (nonlinear, implicit) scheme for non-equilibrium gas discharge plasma. The electrons and ions are described using drift-diffusion approximation coupled to Poisson equation for the electric field. The electron-energy transport equation is solved to account for electron thermal conductivity, Joule heating, and energy loss of electrons in collisions with neutral species. The rate of electron-induced ionization is a function of electron temperature and could also depend on electron density (important for plasma stratification). The ion and gas temperature are kept constant. The transport equations are discretized using a non-isothermal Scharfetter-Gummel scheme to resolve possible large temperature gradients in the sheaths. We demonstrate the new solver for simulations of direct current (DC) and radiofrequency (RF) discharges. The implicit treatment of the coupled equations allows using large time steps. The full-Newton method (FNM) enables fast nonlinear convergence at each time step, offering significantly improved simulation efficiency. We discuss the selection of time steps for solving different plasma problems. The new solver enables solving several problems we could not solve before with existing software: two- and three-dimensional structures of the entire DC discharges including cathode and anode regions, electric field reversals and double-layer formation, the normal cathode spot and an anode ring, moving striations in diffuse and constricted DC discharges, and standing striations in RF discharges. The developed FNM-ACM technique offers many benefits for tackling the disparity of gas discharge plasma systems' time scales and nonlinearity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A Surrogate-Based Asynchronous Decomposition Technique for Realistic Security-Constrained Optimal Power Flow Problems

Here we present a decomposition approach for obtaining good feasible solutions for the security-constrained, alternating-current, optimal power flow (SC-AC-OPF) problem at an industrial scale and under real-world time and computational limits. The approach was designed while preparing and participating in ARPA-E’s Grid Optimization Competition (GOC) Challenge 1. The challenge focused on a near-real-time version of the SC-AC-OPF problem, where a base operating point is optimized, taking into account possible single-element contingencies, after which the system adapts its operating point following the response of automatic frequency droop controllers and voltage regulators. Our solution approach for this problem relies on state-of-the-art nonlinear programming algorithms, and it employs nonconvex relaxations for complementarity constraints, a specialized two-stage decomposition technique with sparse approximations of recourse terms and contingency ranking and prescreening. The paper describes and justifies our approach and outlines the features of its implementation, including functions and derivatives evaluation, warm-starting strategies, and asynchronous parallelism. We discuss the results of the independent benchmark of our approach by ARPA-E’s GOC team in Challenge 1, where it was found to consistently produce high-quality solutions across a wide range of network sizes and difficulty, and conclude by outlining future extensions of the approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Optomechanical interaction in fibre lasers with micro-optomechanical resonance structures

We report the results of studies of fibre lasers with passive modulators based on light-excited micro-optomechanical resonance structures (micro-oscillators). It is shown that in fibre lasers based on active fibres doped with rare-earth elements (Er, Er/Yb, Yb, Nd), the optomechanical interaction of laser radiation with micro-oscillators of various types (fibre-optic, microvolume) leads to self-oscillations of the characteristics of laser radiation at frequencies of relaxation oscillations and intermode beats synchronised with the frequencies of elastic eigenoscillations of micro-oscillators. It is found that in an ultra-long erbium – ytterbium fibre laser with a nonlinear mirror based on a microcantilever, laser photothermal excitation of the second mode of elastic transverse oscillations of the microcantilever makes it possible to perform passive mode locking exclusively due to Q-switching of the laser cavity. Pulsed lasing with a controlled repetition rate (∼76 kHz), a pulse duration of 2 – 5 μs and an output energy of 0.1 μJ per pulse is implemented. Based on a simplified physical model of the indicated fibre lasers with micro-oscillators, we have developed an approximate mathematical model describing the regimes of passive mode locking of fibre lasers with micro-oscillators that play the role of mirrors with a nonlinear reflection coefficient in the laser cavity. The prospects for the development and application of the considered laser systems are discussed. (control of laser radiation parameters)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

CO 2 storage site characterization using ensemble-based approaches with deep generative models

Estimating spatially distributed properties such as permeability from available sparse measurements is a great challenge in efficient subsurface CO 2 storage operations. In this paper, a deep generative model that can accurately capture complex subsurface structure is tested with an ensemble-based inversion method for accurate and accelerated characterization of CO 2 storage sites. We chose Wasserstein Generative Adversarial Network with Gradient Penalty (WGAN-GP) for its realistic reservoir property representation and Ensemble Smoother with Multiple Data Assimilation (ES-MDA) for its robust data fitting and uncertainty quantification capability. WGAN-GP are trained to generate high-dimensional permeability fields from a low-dimensional latent space and ES-MDA then updates the latent variables by assimilating available measurements. Several subsurface site characterization examples including Gaussian, channelized, and fractured reservoirs are used to evaluate the accuracy and computational efficiency of the proposed method and the main features of the unknown permeability fields are characterized accurately with reliable uncertainty quantification. Furthermore, the estimation performance is compared with a widely-used variational, i.e., optimization-based, inversion approach, and the proposed approach outperforms the variational inversion method in several benchmark cases. We explain such superior performance by visualizing the objective function in the latent space: because of nonlinear and aggressive dimension reduction via generative modeling, the objective function surface becomes extremely complex while the ensemble approximation can smooth out the multi-modal surface during the minimization. This suggests that the ensemble-based approach works well over the variational approach when combined with deep generative models at the cost of forward model runs unless convergence-ensuring modifications are implemented in the variational inversion.

42 ENGINEERING↗
DOE OSTI2024

Analysis of a Computational Framework for Bayesian Inverse Problems: Ensemble Kalman Updates and MAP Estimators under Mesh Refinement

This paper analyzes a popular computational framework to solve infinite-dimensional Bayesian inverse problems, discretizing the prior and the forward model in a finite-dimensional weighted inner product space. We demonstrate the benefit of working on a weighted space by establishing operator-norm bounds for finite element and graph-based discretizations of Matérn-type priors and deconvolution forward models. For linear-Gaussian inverse problems, we develop a general theory to characterize the error in the approximation to the posterior. We also embed the computational framework into ensemble Kalman methods and MAP estimators for nonlinear inverse problems. Furthermore, our operator-norm bounds for prior discretizations guarantee the scalability and accuracy of these algorithms under mesh refinement.

Bayesian inverse problem↗
DOE OSTI2025

Universal control in bosonic systems with weak Kerr nonlinearities

Resonators with weak single-photon self-Kerr nonlinearities can theoretically be used to prepare Fock states in the presence of a loss much larger than their nonlinearities. Two necessary ingredients are large displacements and a two-photon (parametric) drive. Here, in this study, we find that these systems can be controlled to achieve any desired gate operation in a finite-dimensional subspace (whose dimensionality can be chosen at will). Moreover, we show that the two-photon driving requirement can be relaxed and that full controllability is achievable with only one-photon (linear) drives. We make use of both Trotter-Suzuki decompositions and gradient-based optimization to find control pulses for a desired gate, which reduces the computational overhead by using a small blockaded subspace. We also discuss the infidelity arising from input power limitations in realistic settings, as well as from corrections to the rotating-wave approximation. Our universal control protocol opens the possibility for quantum information processing using a wide range of lossy systems with weak nonlinearities.

Yuan, Ming [Univ. of Chicago, IL (United States)] ↗
DOE OSTI2021

Approximation of the Frequency-Amplitude Curve Using the Homotopy Analysis Method

The Frequency-Amplitude (F-A) curve on power system oscillation under a large disturbance characterizes how a natural oscillation mode transitions to nonlinear oscillations with growing amplitudes and decaying frequencies. The existing formulation of the F-A curve is derived by solving elliptical integrals on oscillation of a single-machine-infinite-bus equivalent about the targeted oscillation mode. The formula is in a form of infinite series and needs to sum a large number of terms for satisfactory accuracy. This paper introduces an explicit, approximate expression obtained from the Homotopy Analysis Method on the F-A curve. The proposed F-A curve expression is derived from an SMIB system and verified on the IEEE 3-machine 9-bus system to show how the oscillation frequency of a dominant mode varies with oscillation amplitude under large disturbances.

frequency-amplitude (F-A) curve↗