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

A robust and efficient generalized real cascade model

This paper describes a new separation cascade model that allows for arbitrary cascade connectivity and includes the effects of mixing, material losses and generalized separation element performance. The model also accommodates mixtures with an arbitrary number of component species. The mathematical and numerical framework of this model includes both time-dependent and steady state cascade simulation. Robust numerical approximations of the continuous mathematical model are accomplished using the Sundials suite of nonlinear algebraic and differential-algebraic equation solvers. Verification of the model is accomplished using manufactured solutions while validation is demonstrated by comparison to published experimental results and benchmarks.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Estimating the state of a geophysical system with sparse observations: time delay methods to achieve accurate initial states for prediction

The problem of forecasting the behavior of a complex dynamical system through analysis of observational time-series data becomes difficult when the system expresses chaotic behavior and the measurements are sparse, in both space and/or time. Despite the fact that this situation is quite typical across many fields, including numerical weather prediction, the issue of whether the available observations are "sufficient" for generating successful forecasts is still not well understood. An analysis by Whartenby et al. (2013) found that in the context of the nonlinear shallow water equations on a β plane, standard nudging techniques require observing approximately 70 % of the full set of state variables. Here we examine the same system using a method introduced by Rey et al. (2014a), which generalizes standard nudging methods to utilize time delayed measurements. Here, we show that in certain circumstances, it provides a sizable reduction in the number of observations required to construct accurate estimates and high-quality predictions. In particular, we find that this estimate of 70 % can be reduced to about 33 % using time delays, and even further if Lagrangian drifter locations are also used as measurements.

58 GEOSCIENCES↗
DOE OSTI2024

A Data-Driven Approach for High-Impedance Fault Localization in Distribution Systems: Preprint

Accurate and quick identification of high-impedance faults (HIFs) is critical for the reliable operation of distribution systems. Unlike other faults in power grids, HIFs are very difficult to detect by conventional overcurrent relays due to the low fault current. Although HIFs can be affected by various factors, the voltage-current characteristics can substantially imply how the system responds to the disturbance and thus provides opportunities to effectively localize HIFs. In this work, we propose a data-driven approach for the identification of HIF events. To tackle the nonlinearity of the voltage-current trajectory, first, we formulate optimization problems to approximate the trajectory with piecewise functions. Then we collect the function features of all segments as inputs and use the support vector machine approach to efficiently identify HIFs at different locations. Numerical studies on the IEEE 123-node test feeder demonstrate the validity and accuracy of the proposed approach for real-time HIF identification.

explainable artificial intelligence↗
DOE OSTI2016

Local Voltage Control in Distribution Networks: A Game-Theoretic Perspective

Inverter-based voltage regulation is gaining importance to alleviate emerging reliability and power-quality concerns related to distribution systems with high penetration of photovoltaic (PV) systems. This paper seeks contribution in the domain of reactive power compensation by establishing stability of local Volt/VAr controllers. In lieu of the approximate linear surrogate used in the existing work, the paper establishes existence and uniqueness of an equilibrium point using nonlinear AC power flow model. Key to this end is to consider a nonlinear dynamical system with non-incremental local Volt/VAr control, cast the Volt/VAr dynamics as a game, and leverage the fixed-point theorem as well as pertinent contraction mapping argument. Numerical examples are provided to complement the analytical results.

14 SOLAR ENERGY↗
DOE OSTI2016

The anisotropic redistribution of free energy for gyrokinetic plasma turbulence in a Z-pinch

For a Z-pinch geometry, we report on the nonlinear redistribution of free energy across scales perpendicular to the magnetic guide field, for a turbulent plasma described in the framework of gyrokinetics. The analysis is performed using a local flux-surface approximation, in a regime dominated by electrostatic fluctuations driven by the entropy mode, with both ion and electron species being treated kinetically. To explore the anisotropic nature of the free energy redistribution caused by the emergence of zonal flows, we use a polar coordinate representation for the field-perpendicular directions and define an angular density for the scale flux. Positive values for the classically defined (angle integrated) scale flux, which denote a direct energy cascade, are shown to be also composed of negative angular sections, a fact that impacts our understanding of the backscatter of energy and the way in which it enters the modeling of sub-grid scales for turbulence. A definition for the flux of free energy across each perpendicular direction is introduced as well, which shows that the redistribution of energy in the presence of zonal flows is highly anisotropic.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Galerkin v. least-squares Petrov–Galerkin projection in nonlinear model reduction

Least-squares Petrov–Galerkin (LSPG) model-reduction techniques such as the Gauss–Newton with Approximated Tensors (GNAT) method have shown promise, as they have generated stable, accurate solutions for large-scale turbulent, compressible flow problems where standard Galerkin techniques have failed. Furthermore, there has been limited comparative analysis of the two approaches. This is due in part to difficulties arising from the fact that Galerkin techniques perform optimal projection associated with residual minimization at the time-continuous level, while LSPG techniques do so at the time-discrete level. This work provides a detailed theoretical and computational comparison of the two techniques for two common classes of time integrators: linear multistep schemes and Runge–Kutta schemes. We present a number of new findings, including conditions under which the LSPG ROM has a time-continuous representation, conditions under which the two techniques are equivalent, and time-discrete error bounds for the two approaches. Perhaps most surprisingly, we demonstrate both theoretically and computationally that decreasing the time step does not necessarily decrease the error for the LSPG ROM; instead, the time step should be ‘matched’ to the spectral content of the reduced basis. In numerical experiments carried out on a turbulent compressible-flow problem with over one million unknowns, we show that increasing the time step to an intermediate value decreases both the error and the simulation time of the LSPG reduced-order model by an order of magnitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

A memory-distributed quasi-Newton solver for nonlinear programming problems with a small number of general constraints

We approach the problem of parallelizing state-of-the-art nonlinear programming optimization algorithms.Specifically, we focus on parallelizing quasi-Newton interior-point methods that use limited-memory secant Hessian approximations. Such interior-point methods are known to have better convergence properties and to be more effective on large-scale problems than gradient-based and derivative-free optimization algorithms. We target nonlinear and potentially nonconvex optimization problems with an arbitrary number of bound constraints and a small number of general equality and inequality constraints on the optimization variables. These problems occur for example in the form of optimal control, optimal design, and inverse problems governed by ordinary or partial differential equations, whenever they are expressed in a “reduced-space” optimization approach. We introduce and analyze the time and space complexity of a decomposition method for solving the quasi-Newton linear systems that leverages the fact that the quasi-Newton Hessian matrix has a small number of dense blocks that border a low-rank update of a diagonal matrix. This enables an efficient parallelization on memory-distributed computers of the iterations of the optimization algorithm, a state-of-the-art filter line-search interior-point algorithm by Wächter et. al. We illustrate the efficiency of the proposed method by solving structural topology optimization problems on up to 4608 cores on a parallel machine.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Aeroelastic oscillations of a cantilever with structural nonlinearities: theory and numerical simulation.

Our study details the derivation of the nonlinear equations of motion for the axial, biaxial bending and torsional vibrations of an aeroelastic cantilever undergoing rigid body (pitch) rotation at the base. The primary attenstion is focussed on the geometric nonlinearities of the system, whereby the aeroelastic load is modeled by the theory of linear quasisteady aerodynamics. This modelling effort is intended to mimic the wind-tunnel experimental setup at the Royal Military College of Canada. While the derivation closely follows the work of Hodges and Dowell [1] for rotor blades, this aeroelastic system contains new inertial terms which stem from the fundamentally different kinematics than those exhibited by helicopter or wind turbine blades. Using the Hamilton’s principle, a set of coupled nonlinear partial differential equations (PDEs) and an ordinary differential equation (ODE) are derived which describes the coupled axial-bending-bending-torsion-pitch motion of the aeroelastic cantilever with the pitch rotation. The finite dimensional approximation of the coupled system of PDEs are obtained using the Galerkin projection, leading to a coupled system of ODEs. Subsequently, these nonlinear ODEs are solved numerically using the built-in MATLAB implicit ODE solver and the associated numerical results are compared with those obtained using Houbolt’s method. It is demonstrated that the system undergoes coalescence flutter, leading to a limit cycle oscillation (LCO) due to coupling between the rigid body pitching mode and teh flexible mode arising from the flapwise bending motion.

42 ENGINEERING↗
DOE OSTI2015

On the Existence of Regular Global Attractor for p-Laplacian Evolution Equation

In this study, we consider the nonlinear evolution equation of parabolic type u{sub t}−div(|∇u|{sup p−2}∇u)+f(u)=g. We analyze the long time dynamics (in the sense of global attractors) under very general conditions on the nonlinearity f. Since we do not assume any polynomial growth condition on it, the main difficulty arises at first in the proof of well-posedness. In fact, the very first contribution to this problem is a pioneering paper (Efendiev and Ôtano, Differ Int Equ 20:1201–1209, 2007) where the well-posedness result has been shown by exploiting the technique from the theory of maximal monotone operators. However, from some physical aspects, to obtain the solution in variational sense might be demanding which requires limiting procedure on the approximate solutions. In this work, we are interested in variational (weak) solution. The critical issue in the proof of well-posedness is to deal with the limiting procedure on f which is overcome utilizing the weak convergence tecniques in Orlicz spaces (Geredeli and Khanmamedov, Commun Pure Appl Anal 12:735–754, 2013; Krasnosel’skiĭ and Rutickiĭ, Convex functions and Orlicz spaces, 1961). Then, proving the existence of the global attractors in L{sup 2}(Ω) and in more regular space W{sub 0}{sup 1,p}(Ω), we show that they coincide. In addition, if f is monotone and g=0, we give an explicit estimate of the decay rate to zero of the solution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Three-dimensional modeling and numerical predictions of multimodal nonlinear behavior in damaged concrete blocks

In this work, the multimodal nonlinear elastic behavior of concrete, which is representative of a consolidated granular material, is modeled numerically. Starting from a local three-dimensional softening law, the initial stiffness properties are re-estimated according to the local strain field. The experiments deal with samples of thermally damaged concrete blocks successively excited around their first three modes of vibration. The geometry of these samples cannot be described by a one-dimensional approximation in these experiments where compressional and shear motions are strongly coupled. Despite this added complexity, the nonlinear behavior for the three modes of vibration of the samples is well captured by the simulations using a single scalar nonlinear parameter appropriately integrated into the elasticity equations. Finally, it is shown that without sufficient attention paid to the latter, the conclusions would have brought erroneous statements such as nonlinearity dispersion or strain type dependence.

58 GEOSCIENCES↗
DOE OSTI2015

Hydromagnetic flow of third grade nanofluid with viscous dissipation and flux conditions

This article investigates the magnetohydrodynamic flow of third grade nanofluid with thermophoresis and Brownian motion effects. Energy equation is considered in the presence of thermal radiation and viscous dissipation. Rosseland’s approximation is employed for thermal radiation. The heat and concentration flux conditions are taken into account. The governing nonlinear mathematical expressions of velocity, temperature and concentration are converted into dimensionless expressions via transformations. Series solutions of the dimensionless velocity, temperature and concentration are developed. Convergence of the constructed solutions is checked out both graphically and numerically. Effects of interesting physical parameters on the temperature and concentration are plotted and discussed in detail. Numerical values of skin-friction coefficient are computed for the hydrodynamic and hydromagnetic flow cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Temperature dependence of the Goos–Hänchen shift in the nonlinear metal-dielectric nanocomposites

In this paper the temperature-dependent Goos–Hänchen (GH) shift of the nonlinear nanocomposites made of the ellipsoidal nonlinear metal nanoparticles and ellipsoidal dielectric nanoparticles is investigated theoretically in the range of the visible wavelengths. We used the effective medium approximation in the quasi-static regime to calculate the effective permittivity of a composite (ε{sub eff}) and the stationary phase method is utilized to investigate the GH shifts. Spectral representation theory has been used to obtain the local electric field in the nonlinear metal nanoparticles. The temperature impact on the GH shift is explained by the temperature-dependent permittivity of nonlinear metal nanoparticles in the range of 300–1200 K. Our numerical analysis revealed that in addition to the magnitude of the GH shift, the lower and upper threshold of the GH bistability as a nonlinear response will be affected by a temperature increment. So that, the width of the bistability, i.e., the difference between the low and high threshold of bistability, can be tuned by temperature. It is important to note that the GH shift can be switched from positive to negative values in high temperatures for a given physical parameters. This effect can open the interesting possibility of optical switching in the considered structures. Moreover, the occurrence of optical bistability as a function of the volume fraction of nanoparticles and depolarization factor (indicating the shape of nanoparticles) are discussed. Finally, the influence of the wavelength and intensity of the incident light on the nonlinear GH shift is presented. Optimizing the results of the present work may be practically utilized in the designing of optical switching-based devices including metal-dielectric nanocomposites.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2017

Evolving Waves and Turbulence in the Outer Corona and Inner Heliosphere: The Accelerating Expanding Box

Alfvénic fluctuations in the solar wind display many properties reflecting an ongoing nonlinear cascade, e.g., a well-defined spectrum in frequency, together with some characteristics more commonly associated with the linear propagation of waves from the Sun, such as the variation of fluctuation amplitude with distance, dominated by solar wind expansion effects. Therefore, both nonlinearities and expansion must be included simultaneously in any successful model of solar wind turbulence evolution. Because of the disparate spatial scales involved, direct numerical simulations of turbulence in the solar wind represent an arduous task, especially if one wants to go beyond the incompressible approximation. Indeed, most simulations neglect solar wind expansion effects entirely. Here we develop a numerical model to simulate turbulent fluctuations from the outer corona to 1 au and beyond, including the sub-Alfvénic corona. The accelerating expanding box (AEB) extends the validity of previous expanding box models by taking into account both the acceleration of the solar wind and the inhomogeneity of background density and magnetic field. Our method incorporates a background accelerating wind within a magnetic field that naturally follows the Parker spiral evolution using a two-scale analysis in which the macroscopic spatial effect coupling fluctuations with background gradients becomes a time-dependent coupling term in a homogeneous box. In this paper we describe the AEB model in detail and discuss its main properties, illustrating its validity by studying Alfvén wave propagation across the Alfvén critical point.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Fully localised nonlinear energy growth optimals in pipe flow

A new, fully localised, energy growth optimal is found over large times and in long pipe domains at a given mass flow rate. This optimal emerges at a threshold disturbance energy below which a nonlinear version of the known (streamwise-independent) linear optimal [P. J. Schmid and D. S. Henningson, “Optimal energy density growth in Hagen-Poiseuille flow,” J. Fluid Mech. 277, 192–225 (1994)] is selected and appears to remain the optimal up until the critical energy at which transition is triggered. The form of this optimal is similar to that found in short pipes [Pringle et al., “Minimal seeds for shear flow turbulence: Using nonlinear transient growth to touch the edge of chaos,” J. Fluid Mech. 702, 415–443 (2012)], but now with full localisation in the streamwise direction. This fully localised optimal perturbation represents the best approximation yet of the minimal seed (the smallest perturbation which is arbitrarily close to states capable of triggering a turbulent episode) for “real” (laboratory) pipe flows. Dependence of the optimal with respect to several parameters has been computed and establishes that the structure is robust.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Stochastic behavior of nanoscale dielectric wall buckling

The random buckling patterns of nanoscale dielectric walls are analyzed using a nonlinear multi-scale stochastic method that combines experimental measurements with simulations. The dielectric walls, approximately 200 nm tall and 20 nm wide, consist of compliant, low dielectric constant (low-k) fins capped with stiff, compressively stressed TiN lines that provide the driving force for buckling. The deflections of the buckled lines exhibit sinusoidal pseudoperiodicity with amplitude fluctuation and phase decorrelation arising from stochastic variations in wall geometry, properties, and stress state at length scales shorter than the characteristic deflection wavelength of about 1000 nm. The buckling patterns are analyzed and modeled at two length scales: a longer scale (up to 5000 nm) that treats randomness as a longer-scale measurable quantity, and a shorter-scale (down to 20 nm) that treats buckling as a deterministic phenomenon. Statistical simulation is used to join the two length scales. Through this approach, the buckling model is validated and material properties and stress states are inferred. In particular, the stress state of TiN lines in three different systems is determined, along with the elastic moduli of low-k fins and the amplitudes of the small-scale random fluctuations in wall properties—all in the as-processed state. The important case of stochastic effects giving rise to buckling in a deterministically sub-critical buckling state is demonstrated. The nonlinear multiscale stochastic analysis provides guidance for design of low-k structures with acceptable buckling behavior and serves as a template for how randomness that is common to nanoscale phenomena might be measured and analyzed in other contexts.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

On differential singular perturbation problems: a simple variational approach

An alternative approach for the analysis and the numerical approximation of singular perturbation problems, using a variational framework, is presented. It is based on the natural idea of minimizing the residual of the differential equation measured in the L{sup 2} norm. By this strategy, the approximation is based in the solution of a set of linear problems giving the descent step directions of the problem. This is the main advantage of our approach, since we can use stable and convergent methods for linear problems (without assuming the knowledge of good initial guesses used in the approximation of the associated non-linear systems necessary in Newton-type methods). Remember that for this type of problems, we should use implicit methods. We prove that our procedure can never get stuck in local minima, and the error decreases until getting to the original solution independent of the perturbation parameter. Finally, we include some numerical examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Indexing habitat opportunity for juvenile anadromous fishes in tidal-fluvial wetland systems

In intertidal wetlands, annual habitat use by migrating juvenile salmon is limited by both physical and biological constraints. Physical limits are set by inundation patterns and water quality parameters, while biological bounds are determined by species specific migration timing and seasonal residency behaviors. We developed a new metric called the fish habitat opportunity index (FHOI) to quantify the total time per calendar year that salmon could access intertidal sites, based on periods of inundation and juvenile salmon seasonal migration. The index was applied using data collected at a recently reconnected tidal freshwater wetland area during 2006-2009. We found the wetland to be accessible to different salmon species for different lengths of time and periods of the year. On average, chum salmon (Oncorhynchus keta) had remarkably consistent opportunities to use marsh habitat over the 4 years sampled, while opportunities were more variable for Chinook (O. tshawytscha), and coho salmon (O. kisutch). Inundation varied nonlinearly with tidal height, and fish had access to the productive marsh edge for approximately 40-50% of the time during periods of access. The FHOI is widely applicable to wetland systems throughout the tidal fluvial continuum, and provides a metric for assessing fish habitat opportunity that is more realistic than simple wetted area calculations. Finally, the FHOI affords a method of a priori comparison among potential restoration sites to help resource managers predict the likely relative benefit to juvenile salmon.

wetland restoration, estuary, habitat opportunity,↗
DOE OSTI2019

Bootstrap-based confidence interval estimation for thermal security region of bulk power grid

Thermal security region (TSR) is a powerful tool for monitoring and controlling the thermal security of bulk power grid with high penetration of renewable energy. One of the major challenges of TSR in practical application is to obtain the exact analytical expression of the TSR boundary (TSRB) because TSRB is a nonlinear hypersurface described by an implicit function. For this reason, the TSRB is usually approximated by hyperplane instead of getting the exact analytical expression. Traditionally, the coefficients of the hyperplane are estimated with a set of original TSRB points using least square estimation (LSE). However, LSE is a point estimation method, which is incapable of evaluating the quality of the estimation. In addition, the accuracy and computational efficiency of the hyperplane approximation are influenced by the number of TSRB points significantly. Here, a bootstrap based confidence interval estimation is proposed to estimate not only the coefficients of TSRB approximation hyperplane, but also the standard deviations and confidence intervals of the coefficients for evaluating the quality and reliability of the approximation results. First, empirical distribution functions (EDFs) of the coefficients of TSRB approximation hyperplane are approximated from a set of original TSRB points by using residual resampling bootstrap method. Then, the EDFs of the coefficients are employed to estimate the coefficients of the approximation hyperplane. Meanwhile, the standard deviations and confidence intervals of the estimated coefficients of the hyperplanes are also calculated for evaluating the quality and reliability of the approximation. The proposed approach is tested on the China Southern Power Grid (CSG). Results of simulations validate the accuracy and efficiency of the proposed approach in approximating the TSRB.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗