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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Spatiotemporal spread of perturbations in power-law models at low temperatures: Exact results for classical out-of-time-order correlators

Here, we present exact results for the classical version of the out-of-time-order commutator (OTOC) for a family of power-law models consisting of N particles in one dimension and confined by an external harmonic potential. These particles are interacting via power-law interaction of the form ∝ $Σ_{i,j=(i≠j)}^{N}$ |x i -x j | -k for every k>1 where x i is the position of the i th particle. We present numerical results for the OTOC for finite N at low temperatures and short enough times so that the system is well approximated by the linearized dynamics around the many-body ground state. In the large-N limit, we compute the ground-state dispersion relation in the absence of external harmonic potential exactly and use it to arrive at analytical results for OTOC. We find excellent agreement between our analytical results and the numerics. We further obtain analytical results in the limit where only linear and leading nonlinear (in momentum) terms in the dispersion relation are included. The resulting OTOC is in agreement with numerics in the vicinity of the edge of the “light cone.” We find remarkably distinct features in OTOC below and above k=3 in terms of going from non-Airy behavior (1 < k < 3) to an Airy universality class (k>3). We present certain additional rich features for the case k=2 that stem from the underlying integrability of the Calogero-Moser model. We present a field theory approach that also assists in understanding certain aspects of OTOC such as the sound speed. Our findings are a step forward towards a more general understanding of the spatiotemporal spread of perturbations in long-range interacting systems.

1-dimensional systems↗
DOE OSTI2025

Low‐dimensional manifold learning for uncertainty quantification in complex multi‐scale stochastic systems

Broadly speaking, the goals of the project are to develop techniques to use manifold learning to develop reduced‐order and surrogate models for "hyper‐reduction" of very high‐dimensional complex multi‐scale systems. This is being achieved by employing a newly proposed form of manifold projection and learning that leverages recent advancements in computational geometry and data‐driven modeling. In particular, we are applying a manifold projection technique to project the solutions of very high‐dimensional systems onto the so‐called Grassmannmanifold, a Reimannian manifold comprised of orthonormal matrices. We then apply data‐driven machine learning techniques to classify the solutions on the manifold (e.g. clustering techniques) according to their proximity on the manifold and leverage a further nonlinear dimension reduction to organize the structured data on the manifold. Finally, we are developing novel techniques that enable us to directly interpolate the hyper‐reduced data such that we can predict the solution of the complex, high‐ dimensional system without need to call the full expensive computational model. Given their adherence to the underlying structure of the solution of the physical system, it is expected that these approximate solutions will be sufficiently constrained so as to (approximately) adhere to physical principles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Synthesis, crystal structure refinement, and nonlinear-optical properties of CaB{sub 3}O{sub 5}(OH): Comparative crystal chemistry of calcium triborates

Calcium triborate CaB{sub 3}O5(OH) obtained by hydrothermal synthesis in the Ca(OH){sub 2}–H{sub 3}BO{sub 3}–Na{sub 2}CO{sub 3}–KCl system is studied by single-crystal X-ray diffraction. The parameters of the orthorhombic unit cell are as follows: a = 13.490(1), b = 6.9576(3), and c = 4.3930(2) Å; V = 412.32(3) Å{sup 3} and space group Pna2{sub 1}. The structure is refined in the anisotropic approximation of the atomic displacement parameters to R = 4.28% using 972 vertical bar F vertical bar > 4σ(F). It is confirmed that the crystal structure of Ca triborate CaB{sub 3}O{sub 5}(OH) is identical to that described earlier. The hydrogen atom is localized. An SHG signal stronger than that of the quartz standard is registered. The phase transition of calcium triborate into calciborite is found on heating. The comparative crystal-chemical analysis of a series of borates with the general chemical formula 2CaO · 3B{sub 2}O{sub 3} · nH{sub 2}O (n = 0–13) with the constant CaO: B{sub 2}O{sub 3}= 2: 3 ratio and variable content of water is performed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Competing mechanisms at vibrated interfaces of density-contrast fluids

Fluid-fluid interfacial instability and subsequent fluid mixing are ubiquitous in nature and engineering. The hydrodynamic instability of fluid interfaces has long centered on the pressure gradient-driven long-wavelength Rayleigh-Taylor instability and the resonance-induced short-wavelength Faraday instability. However, neither instability alone can explain the dynamics when both mechanisms are present. We identify a previously unseen multi-modal instability emerging from their coexistence. When the denser fluid is polydimethylsiloxane, the mixed region at a high density contrast (Atwood number = 0.9) spans a vibration amplitude range approximately twice the gravitational acceleration. Using Floquet stability analysis, we show how vibrations govern transitions between the RT and Faraday instabilities, leading to contention between these instabilities rather than resonant enhancement. Here, the initial transient growth is represented by the exponential modal growth of the most unstable Floquet exponent, along with its accompanying periodic behavior. Direct numerical simulations validate these findings and track interface breakup into the multiscale and nonlinear regimes. Specifically, we show that growing RT modes nonlinearly suppresses Faraday responses even when the initial growth rate of the Faraday instability is 3.63 times that of RT, so a bidirectional competition hinders their sustained coexistence.

Direct numerical simulations↗
DOE OSTI2025

A New Hybrid Quantum-Classical Algorithm for Solving the Unit Commitment Problem

Solving problems related to planning and operations of large-scale power systems is challenging on classical computers due to their inherent nature as mixed-integer and nonlinear problems. Quantum computing provides new avenues to approach these problems. We develop a hybrid quantum-classical algorithm for the Unit Commitment (UC) problem in power systems which aims at minimizing the total cost while optimally allocating generating units to meet the hourly demand of the power loads. The hybrid algorithm combines a variational quantum algorithm (VQA) with a classical Benders-type heuristic. The resulting algorithm computes approximate solutions to UC in three stages: i) a collection of UC vectors capable meeting the power demand with lowest possible operating costs is generated based on VQA; ii) a classical sequential least squares programming (SLSQP) routine is leveraged to find the optimal power level corresponding to a predetermined number of candidate vectors; iii) in the last stage, the approximate solution of UC along with generating units power level combination is given. To demonstrate the effectiveness of the presented method, three different systems with 3 generating units, 10 generating units, and 26 generating units were tested for different time periods. In addition, convergence of the hybrid quantum-classical algorithm for select time periods is proven out on IonQ's Forte system.

Aboumrad, Willie [IonQ, Inc]↗
DOE OSTI2014

Water transport through functionalized nanotubes with tunable hydrophobicity

Molecular dynamics simulations are used to study the occupancy and flow of water through nanotubes comprised of hydrophobic and hydrophilic atoms, which are arranged on a honeycomb lattice to mimic functionalized carbon nanotubes (CNTs). We consider single-file motion of TIP3P water through narrow channels of (6,6) CNTs with varying fractions (f) of hydrophilic atoms. Various arrangements of hydrophilic atoms are used to create heterogeneous nanotubes with separate hydrophobic/hydrophilic domains along the tube as well as random mixtures of the two types of atoms. The water occupancy inside the nanotube channel is found to vary nonlinearly as a function of f, and a small fraction of hydrophilic atoms (f ≈ 0.4) are sufficient to induce spontaneous and continuous filling of the nanotube. Interestingly, the average number of water molecules inside the channel and water flux through the nanotube are less sensitive to the specific arrangement of hydrophilic atoms than to the fraction, f. Two different regimes are observed for the water flux dependence on f – an approximately linear increase in flux as a function of f for f < 0.4, and almost no change in flux for higher f values, similar to the change in water occupancy. We are able to define an effective interaction strength between nanotube atoms and water's oxygen, based on a linear combination of interaction strengths between hydrophobic and hydrophilic nanotube atoms and water, that can quantitatively capture the observed behavior.

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

Sparse maps—A systematic infrastructure for reduced-scaling electronic structure methods. II. Linear scaling domain based pair natural orbital coupled cluster theory

Domain based local pair natural orbital coupled cluster theory with single-, double-, and perturbative triple excitations (DLPNO-CCSD(T)) is a highly efficient local correlation method. It is known to be accurate and robust and can be used in a black box fashion in order to obtain coupled cluster quality total energies for large molecules with several hundred atoms. While previous implementations showed near linear scaling up to a few hundred atoms, several nonlinear scaling steps limited the applicability of the method for very large systems. In this work, these limitations are overcome and a linear scaling DLPNO-CCSD(T) method for closed shell systems is reported. The new implementation is based on the concept of sparse maps that was introduced in Part I of this series [P. Pinski, C. Riplinger, E. F. Valeev, and F. Neese, J. Chem. Phys. 143, 034108 (2015)]. Using the sparse map infrastructure, all essential computational steps (integral transformation and storage, initial guess, pair natural orbital construction, amplitude iterations, triples correction) are achieved in a linear scaling fashion. In addition, a number of additional algorithmic improvements are reported that lead to significant speedups of the method. The new, linear-scaling DLPNO-CCSD(T) implementation typically is 7 times faster than the previous implementation and consumes 4 times less disk space for large three-dimensional systems. For linear systems, the performance gains and memory savings are substantially larger. Calculations with more than 20 000 basis functions and 1000 atoms are reported in this work. In all cases, the time required for the coupled cluster step is comparable to or lower than for the preceding Hartree-Fock calculation, even if this is carried out with the efficient resolution-of-the-identity and chain-of-spheres approximations. The new implementation even reduces the error in absolute correlation energies by about a factor of two, compared to the already accurate previous implementation.

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

Weighted Composition Operators for Learning Nonlinear Dynamics

Operator theoretic methods in dynamical system have been dominated by the use of Koopman operators and their continuous time counterparts, such as Koopman Generators and Liouville Operators. The advantage gained from their use primarily stems from the ability to extract subspaces and eigenfunctions within a space of observables that are invariant with respect to the Koopman operator over that space. When this occurs, a dynamic mode decomposition of the systems state provides a linear model for the dynamical system. Not all Koopman operators have eigenfunctions that may be exploited in this manner. However, the framework can still be leveraged for approximations using other operators. In this setting, we present a different operator for the study of dynamical systems, the weighted composition operator. These operators are compact for a wide range of dynamics and spaces, and through their interactions with occupation kernels and vector valued kernels, they admit an estimation of the underlying dynamics. Here, this manuscript presents a new algorithm for the data driven study of dynamical systems from data, and also provides two numerical experiments where convergence is achieved as a proof of concept.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Solitary waves in a peridynamic elastic solid

The propagation of large amplitude nonlinear waves in a peridynamic solid is ana- lyzed. With an elastic material model that hardens in compression, sufficiently large wave pulses propagate as solitary waves whose velocity can far exceed the linear wave speed. In spite of their large velocity and amplitude, these waves leave the material they pass through with no net change in velocity and stress. They are nondissipative and nondispersive, and they travel unchanged over large distances. An approximate solution for solitary waves is derived that reproduces the main features of these waves observed in computational simulations. We demonstrate, by numerical studies, that waves interact only weakly with each other when they collide. Finally, we found that wavetrains composed of many non-interacting solitary waves form and propagate under certain boundary and initial conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Matrix-Free High-Performance Saddle-Point Solvers for High-Order Problems in \(\boldsymbol{H}(\operatorname{\textbf{div}})\)

Here, this work describes the development of matrix-free GPU-accelerated solvers for high-order finite element problems in H(div). The solvers are applicable to grad-div and Darcy problems in saddle-point formulation, and have applications in radiation diffusion and porous media flow problems, among others. Using the interpolation–histopolation basis, efficient matrix-free preconditioners can be constructed for the (1, 1)-block and Schur complement of the block system. With these approximations, block-preconditioned MINRES converges in a number of iterations that is independent of the mesh size and polynomial degree. The approximate Schur complement takes the form of an M-matrix graph Laplacian and therefore can be well-preconditioned by highly scalable algebraic multigrid methods. High-performance GPU-accelerated algorithms for all components of the solution algorithm are developed, discussed, and benchmarked. Numerical results are presented on a number of challenging test cases, including the “crooked pipe” grad-div problem, the SPE10 reservoir modeling benchmark problem, and a nonlinear radiation diffusion test case.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2018

Difference-Frequency Generation of THz Radiation via Parametric Three-Wave Interaction in CdTe and ZnTe Crystals

Difference-frequency generation of terahertz (THz) radiation in a collinear single-crystal scheme with distributed feedback implemented by means of parametric interaction of counterpropagating optical waves in a periodically inhomogeneous medium exhibiting quadratic nonlinearity is investigated. Using ZnTe and CdTe crystals optically pumped at λ = 0.77 and 1.06 μm, respectively, as an example, the possibility of obtaining parametric generation of THz radiation with an optical-to-THz conversion efficiency exceeding 10{sup –3} (with an extremely narrow line of oscillation of less than 4 GHz for ZnTe and less than 9 GHz for CdTe at 1 THz) is demonstrated in the undepleted-pump-intensity approximation. The obtained results can be used for construction of narrow-band lasers in the THz frequency range with an output power (in the quasicw regime of oscillation) of about 10{sup 5} W/cm{sup 2}.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Vibrational coherences in charge-transfer dyes: A non-adiabatic picture

Essential-state models efficiently describe linear and nonlinear spectral properties of different families of charge-transfer chromophores. Here, the essential-state machinery is applied to the calculation of the early-stage dynamics after ultrafast (coherent) excitation of polar and quadrupolar chromophores. The fully non-adiabatic treatment of coupled electronic and vibrational motion allows for a reliable description of the dynamics of these intriguing systems. In particular, the proposed approach is reliable even when the adiabatic and harmonic approximations do not apply, such as for quadrupolar dyes that show a multistable, broken-symmetry excited state. Our approach quite naturally leads to a clear picture for a dynamical Jahn-Teller effect in these systems. The recovery of symmetry due to dynamical effects is however disrupted in polar solvents where a static symmetry lowering is observed. More generally, thermal disorder in polar solvents is responsible for dephasing phenomena, damping the coherent oscillations with particularly important effects in the case of polar dyes.

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

Resolving ECRH deposition broadening due to edge turbulence in DIII-D

Microwave heat pulse propagation experiments have demonstrated a correlation between millimeter-scale turbulence and deposition profile broadening of electron cyclotron (EC) waves on the DIII-D tokamak. In a set of discharges in DIII-D, a variation in edge density fluctuations on the mm-scale is associated with 40%–150% broader deposition profiles, expressed in terms of normalized minor radius, as compared with equilibrium ray tracing. The 1D power profile is determined from transport analysis of the electron temperature response to EC power modulation using perturbative analysis with a square wave power modulation at 20–70 Hz, producing a series of Fourier harmonics that are fit collectively to resolve transport. Fitting an integrated heat flux expressed in the Fourier basis of the modulation to diffusive, convective, and coupled transport terms in a linear model can resolve the broadened EC deposition width from the power perturbation to resolve a broadening in each case. The best fit degree of beam broadening observed scales approximately linearly with the Doppler backscattering measured fluctuation level in the steep gradient region. Quantifying the effect of edge fluctuation broadening on EC current drive power needs of future devices will require 3D full-wave codes that can be validated on the current generation of machines. These DIII-D experiments provide a quantitative measure of fluctuation effects and a dataset to benchmark full-wave simulations that can model and eventually predict nonlinear effects neglected by 1D equilibrium beam and ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Quasinormal modes for non-minimally coupled scalar fields in regular black hole spacetimes: Grey-body factors, area spectrum and shadow radius

Black hole perturbation theory is a useful approach to study interactions between black holes and fundamental fields. A particular class of black hole solutions arising out of modification of Einstein’s general theory of relativity are regular black holes (RBHs) which can be constructed using a nonlinear electrodynamic Lagrangian. Because of their importance, we are interested in studying the behavior of three kinds of such RBHs under perturbations generated by an external field. Indeed, we investigate the quasinormal modes (QNMs) of a massive scalar field propagating near the RBHs which is non-minimally coupled to the Ricci scalar-tensor of background geometry. We will attempt to find the low-lying quasinormal frequencies of the perturbations by using WKB approximation. We shall also study the relationship between the QNMs and some characteristic properties of black holes such as grey-body factors, area quantization, and shadow radius.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Understanding the breakdown of classic two-phase theory and spray atomization at engine-relevant conditions

We present a generalized framework for multi-component liquid injections to understand and predict the breakdown of classic two-phase theory and spray atomization at engine-relevant conditions. The analysis focuses on the thermodynamic structure and the immiscibility state of representative gas-liquid interfaces. The most modern form of Helmholtz energy mixture state equation is utilized which exhibits a unique and physically-consistent behavior over the entire two-phase regime of fluid densities. It is combined with generalized models for non-linear Gradient Theory and for liquid injections to quantify multi-component two-phase interface structures in global thermal equilibrium. Then, the Helmholtz free energy is minimized which determines the interfacial species distribution as a consequence. This minimal free energy state is demonstrated to validate the underlying assumptions of classic two-phase theory and spray atomization. However, under certain engine-relevant conditions for which corroborating experimental data is presented, this requirement for interfacial thermal equilibrium becomes unsustainable. A rigorously derived probability density function quantifies the ability of the interface to develop internal spatial temperature gradients in the presence of significant temperature differences between injected liquid and ambient gas. Then, the interface can no longer be viewed as an isolated system at minimal free energy. Instead, the interfacial dynamics become intimately connected to those of the separated homogeneous phases. Hence, the interface transitions toward a state in local equilibrium whereupon it becomes a dense-fluid mixing layer. A new conceptual view of a transitional liquid injection process emerges from a transition time scale analysis. Close to the nozzle exit, the two-phase interface still remains largely intact and more classic two-phase processes prevail as a consequence. Further downstream, however, the transition to dense- fluid mixing generally occurs before the liquid length is reached. The significance of the presented modeling expressions is established by a direct comparison to a reduced model, which utilizes widely-applied approximations but fundamentally fails to capture the physical complexity discussed in this paper.

42 ENGINEERING↗
DOE OSTI2020

Forward sensitivity approach for estimating eddy viscosity closures in nonlinear model reduction

In this paper, we propose a variational approach to estimate eddy viscosity using forward sensitivity method (FSM) for closure modeling in nonlinear reduced order models. FSM is a data assimilation technique that blends model's predictions with noisy observations to correct initial state and/or model parameters. We apply this approach on a projection based reduced order model (ROM) of the one-dimensional viscous Burgers equation with a square wave defining a moving shock, and the two-dimensional vorticity transport equation formulating a decay of Kraichnan turbulence. We investigate the capability of the approach to approximate an optimal value for eddy viscosity with different measurement configurations. Specifically, we show that our approach can sufficiently assimilate information either through full field or sparse noisy measurements to estimate eddy viscosity closure to cure standard Galerkin reduced order model (GROM) predictions. Therefore, our approach provides a modular framework to correct forecasting error from a sparse observational network on a latent space. We highlight that the proposed GROM-FSM framework is promising for emerging digital twin applications, where real-time sensor measurements can be used to update and optimize surrogate model's parameters.

42 ENGINEERING↗