SEARCH · Search NASA
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.
In this software, we implement an adaptive, multivariate partitioning algorithm for solving mixed-integer nonlinear programs (MINLP) to global optimality. The algorithm combines ideas that exploit the structure of convex relaxations to MINLPs and bound tightening procedures
It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.
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.
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.
Quasilinear turbulent transport models are a successful tool for prediction of core tokamak plasma profiles in many regimes. Their success hinges on the reproduction of local nonlinear gyrokinetic fluxes. We focus on significant progress in the quasilinear gyrokinetic transport model QuaLiKiz (Bourdelle et al 2016 Plasma Phys. Control. Fusion 58 014036), which employs an approximated solution of the mode structures to significantly speed up computation time compared to full linear gyrokinetic solvers. Optimisation of the dispersion relation solution algorithm within integrated modelling applications leads to flux calculations $\times {10}^{6-7}$ faster than local nonlinear simulations. This allows tractable simulation of flux-driven dynamic profile evolution including all transport channels: ion and electron heat, main particles, impurities, and momentum. Furthermore, QuaLiKiz now includes the impact of rotation and temperature anisotropy induced poloidal asymmetry on heavy impurity transport, important for W-transport applications. Application within the JETTO integrated modelling code results in 1 s of JET plasma simulation within 10 h using 10 CPUs. Simultaneous predictions of core density, temperature, and toroidal rotation profiles for both JET hybrid and baseline experiments are presented, covering both ion and electron turbulence scales. The simulations are successfully compared to measured profiles, with agreement mostly in the 5%–25% range according to standard figures of merit. QuaLiKiz is now open source and available at www.qualikiz.com.
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.
A computerized analysis of the nonlinear behavior of fibrous composite laminates including axial loading, thermal loading, temperature dependent properties, and edge effects is presented. Ramberg-Osgood approximations are used to represent lamina stress-strain behavior and percent retention curves are employed to model the variation of properties with temperature. Balanced, symmetric laminates comprised of either boron/epoxy, graphite/epoxy, or borsic-aluminum are analyzed using a quasi-three-dimensional finite element analysis. Results are presented for the interlaminar stress distributions in cross-ply, angle-ply, and more complex laminates. Nonlinear stress-strain curves for a variety of composite laminates in tension and compression are obtained and compared to other existing theories and experimental results.
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.
A parallel trajectory optimization algorithm is being developed. One possible mission is to provide real-time, on-line guidance for the National Aerospace Plane. The algorithm solves a discrete-time problem via the augmented Lagrangian nonlinear programming algorithm. The algorithm exploits the dynamic programming structure of the problem to achieve parallelism in calculating cost functions, gradients, constraints, Jacobians, Hessian approximations, search directions, and merit functions. Special additions to the augmented Lagrangian algorithm achieve robust convergence, achieve (almost) superlinear local convergence, and deal with constraint curvature efficiency. The algorithm can handle control and state inequality constraints such as angle-of-attack and dynamic pressure constraints. Portions of the algorithm have been tested. The nonlinear programming core algorithm performs well on a variety of static test problems and on an orbit transfer problem. The parallel search direction algorithm can reduce wall clock time by a factor of 10 for this part of the computation task.
The evolution of the interface between two fluids confined in a rectangular cavity is investigated numerically to predict its transient behavior. These computations address mixing characteristics of fluids under microgravity conditions, particularly g-jitter conditions resulting from aircraft vibration or crew motion, and have applications in solution crystal growth. The two-dimensional formulation employs the Boussinesq approximation and treats the mixing of two fluids as an initial value problem with a prescribed concentration field. For fluid mixing of practical applications inside a cavity, it is shown that nonlinear convective transport can dominate over both viscous and molecular diffusion. However, viscous diffusion can become important for low Reynolds number or certain cavity sizes. A stacking phenomenon is shown to occur for aspect ratios (Ar) approximately between 5 to 10. For aspect ratios in the neighborhood of 0.1 to 0.2, inner cells evolve at the interface with increasing time. For a square cavity of Ar = 1, chaotic mixing of the fluid appears to occur for a Stokes-Reynolds number greater than about 5.
We study the influence of a nonstationary nonlinear response of a medium on self-compression of soliton-like laser pulses in active fibres with a finite gain bandwidth. Based on the variational approach, we qualitatively analyse the self-action of the wave packet in the system under consideration in order to classify the main evolution regimes and to determine the minimum achievable laser pulse duration during self-compression. The existence of stable soliton-type structures is shown in the framework of the parabolic approximation of the gain profile (in the approximation of the Gnizburg – Landau equation). An analysis of the self-action of laser pulses in the framework of the nonlinear Schrödinger equation with a sign-constant gain profile demonstrate a qualitative change in the dynamics of the wave field in the case of a nonstationary nonlinear response that shifts the laser pulse spectrum from the amplification region and stops the pulse compression. Expressions for a minimum duration of a soliton-like laser pulse are obtained as a function of the problem parameters, which are in good agreement with the results of numerical simulation. (paper)
An application of the extended Kalman filter to delayed systems to estimate the state and time delay is presented. Two nonlinear estimators are discussed and the results compared with those of the Kalman filter. For all the filters considered, the hereditary system was treated with the delay in the pure form and by using Pade approximations of the delay. A summary of the convergence properties of the filters studied is given. The results indicate that the linear filter applied to the delayed system performs inadequately while the nonlinear filters provide reasonable estimates of both the state and the parameters.
Various concepts for advanced syspension systems have been proposed for counteracting gravity loads in ground vibration testing of large space structures. Approximating the flight modes of a low-frequency flexible structure in a ground test requires a very soft suspension system. The dynamic behavior of a passive zero-spring-rate mechanism, sometimes used for such ground testing, is analyzed. This mechanism reduces the stiffness inherent in suspending a test specimen by cables. However, the mechanism is shown to be sensitive to imperfections. Imperfections can initiate nonlinear behavior which becomes more pronounced at lower operating frequencies. Furthermore, large pendular motion of the suspension system couples with the vertical motion, producing additional nonlinearity.
This paper presents a hierarchical adaptive nonlinear iteration method (HANIM) implemented in the NASA computational fluid dynamics code, FUN3D, to improve robustness and computational efficiency. In contrast to the legacy FUN3D iterative solver that relies on an approximate Jacobian, a simple multicolor Gauss-Seidel point-implicit iteration scheme, and linear Courant-Friedrichs-Lewy number (CFL) ramping, HANIM is based upon a hierarchy of modules including preconditioner, generalized conjugate residual, realizability check, nonlinear control, and CFL adaption modules. HANIM performance is systematically compared with the performance of the legacy solver of FUN3D and a baseline solver based on a preconditioner alone. Iterative solutions are compared for three benchmark cases: a subsonic separated flow around a hemisphere cylinder, a supersonic flow through a long duct, and a subsonic flow over the NASA wing-fuselage juncture model. Two Reynolds-averaged Navier-Stokes turbulence models are used in these computations, namely, the negative variant of the linear one-equation Spalart-Allmaras model and its nonlinear extension based on quadratic constitutive relations.
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.
Whistler mode chorus waves are particularly important in outer radiation belt dynamics due to their key role in controlling the acceleration and scattering of electrons over a very wide energy range. The efficiency of wave-particle resonant interactions is defined by whistler wave properties which have been described by the approximation of plane linear waves propagating through the cold plasma of the inner magnetosphere. However, recent observations of extremely high-amplitude whistlers suggest the importance of nonlinear wave-particle interactions for the dynamics of the outer radiation belt. Oblique chorus waves observed in the inner magnetosphere often exhibit drastically nonsinusoidal (with significant power in the higher harmonics) waveforms of the parallel electric field, presumably due to the feedback from hot resonant electrons. We have considered the nature and properties of such nonlinear whistler waves observed by the Van Allen Probes and Time History of Events and Macroscale Interactions define during Substorms in the inner magnetosphere, and we show that the significant enhancement of the wave electrostatic component can result from whistler wave coupling with the beam-driven electrostatic mode through the resonant interaction with hot electron beams. Being modulated by a whistler wave, the electron beam generates a driven electrostatic mode significantly enhancing the parallel electric field of the initial whistler wave. Finally, we confirm this mechanism using a self-consistent particle-in-cell simulation. The nonlinear electrostatic component manifests properties of the beam-driven electron acoustic mode and can be responsible for effective electron acceleration in the inhomogeneous magnetic field.
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.
Abstract We here propose a new statistical method to interpret nonstationary running correlations by decomposing them into a stationary part and a first‐order Taylor expansion approximation for the nonstationary part. Then, this method is applied to investigate the nonstationary behavior of the El Niño–Southern Oscillation (ENSO)‐East Asian winter monsoon (EAWM) relationship, which exhibits prominent multidecadal variations. It is demonstrated that the first‐order approximation of the nonstationary part can be expressed to a large extent by the impact of the nonlinear interaction between the Atlantic Multidecadal Oscillation (AMO) and ENSO (AMO*Niño3.4) on the EAWM. Therefore, the nonstationarity in the ENSO‐EAWM relationship comes predominantly from the impact of an AMO modulation on the ENSO‐EAWM teleconnection via this key nonlinear interaction. This general method can be applied to investigate nonstationary relationships that are often observed between various different climate phenomena.