Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear equations”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Uniform-density Bose-Einstein condensates of the Gross-Pitaevskii equation found by solving the inverse problem for the confining potential

Here, in this work, we consider a “reverse-engineering” approach to construct confining potentials that support exact, constant density kovaton solutions to the classical Gross-Pitaevskii equation (GPE) also known as the nonlinear Schr¨odinger equation (NLSE). In the one-dimensional case, the exact solution is the sum of stationary kink and anti-kink solutions, i.e. a kovaton, and in the overlapping region, the density is constant. In higher dimensions, the exact solutions are generalizations of this wave function. In the absence of self-interactions, the confining potential is similar to a smoothed out finite square well with minima also at the edges. When self-interactions are added, a term proportional to ±gψ*ψ gets added to the confining potential and ±gM, where M is the norm, gets added to the total energy. In the realm of stability analysis, we find (linearly) stable solutions in the case with repulsive self-interactions which also are stable to self-similar deformations. For attractive interactions, however, the minima at the edges of the potential get deeper and a barrier in the center forms as we increase the norm. This leads to instabilities at a critical value of M (related to the number of particles in the BEC). Comparing the stability criteria from Derrick’s theorem and Bogoliubov-de Gennes analysis stability results, we find that both predict stability for repulsive self-interactions and instability at a critical mass M for attractive interactions. However, the numerical analysis gives a much lower critical mass. The numerical analysis shows further that the initial instabilities violate the symmetry x → -x assumed by Derrick’s theorem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Back-to-Back Inclusive Dijets in Deep Inelastic Scattering at Small x : Complete NLO Results and Predictions

We compute the back-to-back dijet cross section in deep inelastic scattering at small x to next-to-leading order (NLO) in the color glass condensate effective field theory. Our result can be factorized into a convolution of the Weizsäcker-Williams gluon transverse-momentum-dependent distribution function (WW gluon TMD) with a universal soft factor and an NLO coefficient function. The soft factor includes both double and single logarithms in the ratio of the relative transverse momentum P ⊥ of the dijet pair to the dijet momentum imbalance q ⊥ ; its renormalization group (RG) evolution is resummed into the Sudakov factor. Likewise, the WW TMD obeys a nonlinear RG equation in x that is kinematically constrained to satisfy both the lifetime and rapidity ordering of the projectile. Exact analytical expressions are obtained for the NLO coefficient function of transversely and longitudinally polarized photons. Our results allow for the first quantitative separation of the dynamics of Sudakov suppression from that of gluon saturation. They can be extended to other final states and provide a framework for precision tests of novel QCD many-body dynamics at the Electron-Ion Collider. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Efficient GPU Implementation of Automatic Differentiation for Computational Fluid Dynamics

Many scientific and engineering applications require repeated calculations of derivatives of output functions with respect to input parameters. Automatic Differentiation (AD) is a method that automates derivative calculations and can significantly speed up code development. In Computational Fluid Dynamics (CFD), derivatives of flux functions with respect to state variables (Jacobian) are needed for efficient solutions of the nonlinear governing equations. AD of flux functions on graphics processing units (GPUs) is challenging as flux computations involve many intermediate variables that create high register pressure and require significant memory traffic because of the need to store the derivatives. This paper presents a forward-mode AD method based on multivariate dual numbers that addresses these challenges and simultaneously reduces the floating-point operation count. The dimension of the multivariate dual numbers is optimized for performance. The flux computations are restructured to minimize the number of temporary variables and reduce register pressure. For effective utilization of memory bandwidth, shared memory is used to store the local flux Jacobian. This AD implementation is compared with several other Jacobian implementations on an NVIDIA V100 GPU (V100). For three-dimensional perfect-gas compressible-flow equations implemented in a practical CFD code, the AD implementation of a flux Jacobian based on multivariate dual numbers of dimension 5 outperforms all other GPU AD implementations on V100. Its performance is comparable with the optimized hand-differentiated version. Finally, the implementation achieves 75% of the peak floating-point throughput and 61 % of the peak global device memory bandwidth usage.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Towards Smart Grids Enhanced Situation Awareness: A Bi-Level Quasi-Static State Estimation Model

Smart Grid situational awareness is provided by Energy Management Systems. A core process of these systems is State Estimation. The great majority of state estimators model the Smart Grid through a set of nonlinear algebraic equations, named the measurement model. Problem formulation considers the Gauss solution. Several model improvements have been presented regarding the Gauss solution, aiming between others to provide measurement noise robustness to the state estimation process. While considerable effort has been focused on such developments, state estimation is still constrained by the implicit modelling error, and thus inevitably vulnerable to cyber-threats. In this work, a state estimation bi-level formal model is presented towards Smart Grids enhanced situational awareness considering the concepts of synthetic measurements and innovation. Comparative test results with the state-of-the-art on the IEEE 14-bus system are presented highlighting improved situational awareness to bad data. Easy-to-implement model, without hard-to-derive parameters, built-on the classic weighted least squares solution, highlight potential aspects for real-life implementation.

cyber security↗
DOE OSTI2020

Efficient and Robust Dynamic Simulation of Power Systems With Holomorphic Embedding

Dynamic simulation is vitally important in power system analysis, but traditional approaches based on numerical integration over small time steps are time-consuming. Also, the Newton-Raphson method suffers from difficulty in convergence when solving nonlinear algebraic equations. In this paper, we propose a novel dynamic simulation approach based on holomorphic embedding. By obtaining a high-order approximation of system dynamics, it achieves a much larger time step and thus enhances the computational efficiency significantly. In addition, the new approach avoids non-convergence issues in solving algebraic equations, which improves robustness. The approach includes flexible modeling of synchronous generators and controllers, and we propose a method for modeling generator coordinate transformations. The approach is tested on the IEEE 39-bus, 10-generator system and a Polish 2383-bus, 327-generator system. The results demonstrate promising computational efficiency and satisfactory numerical robustness for the analysis of large-scale power systems.

42 ENGINEERING↗
DOE OSTI2023

Microscopic activated dynamics theory of the shear rheology and stress overshoot in ultradense glass-forming fluids and colloidal suspensions

We formulate a particle and force level, activated dynamics-based statistical mechanical theory for the continuous startup nonlinear shear rheology of ultradense glass-forming hard sphere fluids and colloidal suspensions in the context of the elastically collective nonlinear Langevin equation approach and a generalized Maxwell model constitutive equation. Activated structural relaxation is described as a coupled local-nonlocal event involving caging and longer range collective elasticity which controls the characteristic stress relaxation time. Theoretical predictions for the deformation-induced enhancement of mobility, the onset of relaxation acceleration at remarkably low values of stress, strain, or shear rate, apparent power law thinning of the steady-state structural relaxation time and viscosity, a nonvanishing activation barrier in the shear thinning regime, an apparent Herschel–Buckley form of the shear rate dependence of the steady-state shear stress, exponential growth of different measures of a yield or flow stress with packing fraction, and reduced fragility and dynamic heterogeneity under deformation were previously shown to be in good agreement with experiments. The central new question we address here is the defining feature of the transient response—the stress overshoot. In contrast to the steady-state flow regime, understanding the transient response requires an explicit treatment of the coupled nonequilibrium evolution of structure, elastic modulus, and stress relaxation time. We formulate a new quantitative model for this aspect in a physically motivated and computationally tractable manner. Theoretical predictions for the stress overshoot are shown to be in good agreement with experimental observations in the metastable ultradense regime of hard sphere colloidal suspensions as a function of shear rate and packing fraction, and accounting for deformation-assisted activated motion appears to be crucial for both the transient and steady-state responses.

Mechanics↗
DOE OSTI2021

Multilevel Spectral Coarsening for Graph Laplacian Problems with Application to Reservoir Simulation

We extend previously developed two-level coarsening procedures for graph Laplacian problems written in a mixed saddle point form to the fully recursive multilevel case. The resulting hierarchy of discretizations gives rise to a hierarchy of upscaled models, in the sense that they provide approximation in the natural norms (in the mixed setting). This property enables us to utilize them in three applications: (i) as an accurate reduced model, (ii) as a tool in multilevel Monte Carlo simulations (in application to finite volume discretizations), and (iii) for providing a sequence of nonlinear operators in a full approximation scheme for solving nonlinear pressure equations discretized by the conservative two-point flux approximation. Finally, we illustrate the potential of the proposed multilevel technique in all three applications on a number of popular benchmark problems used in reservoir simulation.

multilevel Monte Carlo↗
DOE OSTI2022

PETSc TSAdjoint: A Discrete Adjoint ODE Solver for First-Order and Second-Order Sensitivity Analysis

Here, we present a new software system PETSc TSAdjoint for first-order and second order adjoint sensitivity analysis of time-dependent nonlinear differential equations. The derivative calculation in PETSc TSAdjoint is essentially a high-level algorithmic differentiation process. The adjoint models are derived by differentiating the timestepping algorithms and implementing them based on the parallel infrastructure in PETSc. Full differentiation of the library code, including MPI routines, is avoided, and users do not need to derive their own adjoint models for their specific applications. PETSc TSAdjoint can compute the first-order derivative, that is, the gradient of a scalar functional, and the Hessian-vector product, which carries second-order derivative information, while requiring minimal input (a few callbacks) from the users. The adjoint model employs optimal checkpointing schemes in a manner that is transparent to users. Finally, usability, efficiency, and scalability are demonstrated through examples from a variety of applications.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

From the Great Wave of Translation to the Force between Quarks

Here, the chance observation of a novel traveling wave in a canal led over time to the formulation of a nonlinear wave equation—the Korteweg–de Vries equation—that describes strikingly robust disturbances now called solitons. The figure of an isolated soliton corresponds to a reflectionless potential that supports a single bound state in the one-dimensional Schrödinger equation. An appropriate combination of individual solitons yields a symmetric reflectionless potential that supports multiple bound states. Thus, the KdV equation opens the path to solving the inverse scattering problem for a collection of bound states. Applied to the quarkonium spectra, this formalism allows the construction of reflectionless approximations to the confining potentials that account for the force between quarks, and to tests of the flavor-independence of the interquark interaction.

Inverse Scattering↗
DOE OSTI2020

Dynamic Response of a Semiactive Suspension System with Hysteretic Nonlinear Energy Sink Based on Random Excitation by means of Computer Simulation

This paper aims to investigate the property and behavior of the hysteretic nonlinear energy sink (HNES) coupled to a half vehicle system which is a nine-degree-of-freedom, nonlinear, and semiactive suspension system in order to improve the ride comfort and increase the stability in shock mitigation by using the computer simulation method. The HNES model is a semiactive suspension device, which comprises the famous Bouc–Wen (B-W) model employed to describe the force produced by both the purely hysteretic spring and linear elastic spring of potentially negative stiffness connected in parallel, for the half vehicle system. Nine nonlinear motion equations of the half vehicle system are derived in terms of the seven displacements and the two dimensionless hysteretic variables, which are integrated numerically by employing the direct time integration method for studying both the variables of vertical displacements, velocities, accelerations, chassis pitch angle, and the ride comfort and driver safety, respectively, based on the bump and random road inputs of the pseudoexcitation method as excitation signal. Simulation results show that, compared with the HNES model and the magnetorheological (MR) model coupled to the half vehicle system, the ride comfort and stability have been evidently improved. A successful validation process has been performed, which indicated that both the ride comfort and driver safety properties of the HNES model coupled to half vehicle significantly improved.

Chen, Hui↗
DOE OSTI2022

Temporal reflection of an optical pulse from a short soliton: impact of Raman scattering

We study temporal reflection of an optical pulse from the refractive-index barrier created by a short pump soliton inside a nonlinear dispersive medium such as an optical fiber. One feature is that the soliton’s speed changes continuously as its spectrum redshifts because of intrapulse Raman scattering. We use the generalized nonlinear Schrödinger equation to find the shape and spectrum of the reflected pulse. Both are affected considerably by the soliton’s trajectory. The reflected pulse can become considerably narrower compared to the incident pulse under conditions that involve a type of temporal focusing. This phenomenon is explained through space–time duality by showing that the temporal situation is analogous to an optical beam incident obliquely on a parabolic mirror. We obtain an approximate analytic expression for the reflected pulse’s spectrum and use it to derive the temporal version of the transformation law for the

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Parallel Implicit Hydrodynamics for High Explosive Burn Calculations

High explosives in hostile environments will require calculational capabilities that model processes, which evolve on timescales from minutes to nanoseconds. Eventually the HE will begin to move metal. To handle this temporal evolution an implicit hydrodynamics coupled to the chemical release of the HE energy is required. In addition, the use of chemical kinetics to model the transition from, the initially, slow heating of a confined high explosive through to deflagration and on to detonation requires many computational zones to model high explosive engineering systems. This requirement means that a fully parallel implicit hydrodynamics is essential. In this paper we present the calculation of a nonlinear matrix equation for the advanced particle pressure that has been made parallel and implemented in our AMR code, BABBO. This new parallel implicit hydrodynamics has been applied to a cookoff problem, as well as, one and two dimensional shock problems. Results are presented and discussed.

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

Parameterized Neural Ordinary Differential Equations: Applications to Computational Physics Problems

This work proposes an extension of neural ordinary differential equations (NODEs) by introducing an additional set of ODE input parameters to NODEs. This extension allows NODEs to learn multiple dynamics specified by the input parameter instances. Our extension is inspired by the concept of parameterized ordinary differential equations, which are widely investigated in computational science and engineering contexts, where characteristics of the governing equations vary over the input parameters. We apply the proposed parameterized NODEs (PNODEs) for learning latent dynamics of complex dynamical processes that arise in computational physics, which is an essential component for enabling rapid numerical simulations for time-critical physics applications. For this, we propose an encoder-decoder-type framework, which models latent dynamics as PNODEs. We demonstrate the effectiveness of PNODEs with important benchmark problems from computational physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Evaluation of Reduced Order Model for HT-9 Creep and Modifications to Current HT-9 Creep Model in BISON

In order to leverage existing reduced order models (ROM) for the modeling of LANL-developed HT-9 cladding constitutive behavior, mechanistic-based models have been implemented into BISON. It is posited that the increased fidelity of the ROM will show improvement over engineering-scale models in terms of comparison against experimental measurements. INL will assist LANL in implementing an HT9 ROM in BISON. Once the HT-9 mechanistic constitutive model can be leveraged in nuclear performance simulations, the numerical results of fast reactor models will be analyzed and compared against the Fuels Irradiation & Physics Database (FIPD) and separate effects tests to attempt a validation process. INL and LANL will work together to perform any necessary improvement that is identified during the implementation and use processes. Successful completion of the milestone will enable the general use of validated mechanistic ROMs for HT-9 cladding in BISON. This report provides a brief introduction to the Los Alamos Reduced Order Model Applied to Nonlinear Constitutive Equations (LAROMANCE) code, mechanical testing used to illustrate its creep and plastic deformation predictions for HT-9 cladding, and a comparison with existing models implemented into BISON. Current models have been modified to more appropriately account for primary thermal creep and results for corresponding testing are included in this report.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Nominal and adversarial synthetic PMU data for standard IEEE test systems

GridSTAGE (Spatio-Temporal Adversarial scenario GEneration) is a framework for the simulation of adversarial scenarios and the generation of multivariate spatio-temporal data in cyber-physical systems. GridSTAGE is developed based on Matlab and leverages Power System Toolbox (PST) where the evolution of the power network is governed by nonlinear differential equations. Using GridSTAGE, one can create several event scenarios that correspond to several operating states of the power network by enabling or disabling any of the following: faults, AGC control, PSS control, exciter control, load changes, generation changes, and different types of cyber-attacks. Standard IEEE bus system data is used to define the power system environment. GridSTAGE emulates the data from PMU and SCADA sensors. The rate of frequency and location of the sensors can be adjusted as well. Detailed instructions on generating data scenarios with different system topologies, attack characteristics, load characteristics, sensor configuration, control parameters are available in the Github repository - https://github.com/pnnl/GridSTAGE. There is no existing adversarial data-generation framework that can incorporate several attack characteristics and yield adversarial PMU data. The GridSTAGE framework currently supports simulation of False Data Injection attacks (such as a ramp, step, random, trapezoidal, multiplicative, replay, freezing) and Denial of Service attacks (such as time-delay, packet-loss) on PMU data. Furthermore, it supports generating spatio-temporal time-series data corresponding to several random load changes across the network or corresponding to several generation changes. A Koopman mode decomposition (KMD) based algorithm to detect and identify the false data attacks in real-time is proposed in https://ieeexplore.ieee.org/document/9303022. Machine learning-based predictive models are developed to capture the dynamics of the underlying power system with a high level of accuracy under various operating conditions for IEEE 68 bus system. The corresponding machine learning models are available at https://github.com/pnnl/grid_prediction.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

The Spiral Galaxies Flat Rotational Velocity Curve Explained by the Constant Group Velocity of a Nonlinear Density Wave

The rotation velocity curves of stars in galaxies, the motions of pairs of galaxies, and the behavior of galaxies in clusters and super-clusters all indicate that there is a lack of mass on different scales in the universe. In this paper, we derive the expression for rotational velocity using the nonlinear density wave theory considering only stellar components and we show that such theory can support the observed flat rotational velocity curve due to the main property of the soliton wave, which is a constant group velocity of the wave. The surface mass density (SMD) function, used to derive gravitational potential gradient and rotational velocity, is not assumed but rather derived as a solution of the nonlinear Srödinger equation, on the contrary to the widely used, in the literature, exponential disk approximation. Three parameters relevant to the curve shape are the intensities of equilibrium SMD, the amplitude of the wave, and total angular velocity or differential rotation, equivalently. Since the shape of the rotational velocity is highly sensitive to the mentioned parameters, this theory eventually provides a method for a very accurate estimation of galaxy mass and angular velocity as well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗