Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear systems”

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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Shadow Lagrangian dynamics for superfluidity

Motivated by a similar approach for Born-Oppenheimer molecular dynamics, this paper proposes an extended "shadow" Lagrangian density for quantum states of superfluids. The extended Lagrangian contains an additional field variable that is forced to follow the wave function of the quantum state through a rapidly oscillating extended harmonic oscillator. By considering the adiabatic limit for large frequencies of the harmonic oscillator, we can derive the two equations of motions, a Schrödinger-type equation for the quantum state and a wave equation for the extended field variable. The equations are coupled in a nonlinear way, but each equation individually is linear with respect to the variable that it defines. The computational advantage of this new system is that it can be easily discretized using linear time stepping methods, where we propose to use a Crank-Nicolson-type approach for the Schrödinger equation and an extended leapfrog scheme for the wave equation. Furthermore, the difference between the quantum state and the extended field variable defines a consistency error that should go to zero if the frequency tends to infinity. By coupling the time-step size in our discretization to the frequency of the harmonic oscillator we can extract an easily computable consistency error indicator that can be used to estimate the numerical error without additional costs. The findings are illustrated in numerical experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

High-speed fiber-coupled terahertz time-domain spectroscopy system using a bias-free plasmonic source

We demonstrate a high-speed fiber-coupled terahertz time-domain spectroscopy (THz-TDS) system using a bias-free plasmonic terahertz source. The system is powered by a compact, low-power femtosecond fiber laser oscillator that drives a dual-branch fiber amplifier to pump and probe a plasmonic terahertz source and detector, respectively. A polarization-maintaining, fiber-coupled optical delay line is employed to precisely control the temporal offset between the pump and probe beams through a compact, robust platform. Operating in the self-similar amplification regime, the fiber amplifiers are engineered to mitigate nonlinear pulse distortion and dispersion while delivering robust, high-energy pulses at a 1560-nm wavelength, with an output power exceeding 326 mW. In reflection mode, the THz-TDS system achieves a peak dynamic range of 81 dB and a bandwidth greater than 3.5 THz at a scan rate of 3.75 Hz. The utility of the system for rapid sensing is further demonstrated by quantifying the porosity and density of pharmaceutical tablets through measurements of their physical thickness and effective refractive index.

Jiang, Xinghe [University of California, Los Angel↗
DOE OSTI2020

Morphology and Transport of Multivalent Cation-Exchanged Ionomer Membranes Using Perfluorosulfonic Acid–Ce Z+ as a Model System

Perfluorosulfonic acids (PFSAs) are commonly used as solid polymer electrolyte membranes (PEMs) in electrochemical energy devices, where they are vulnerable to attack by radical species during operation. A popular strategy to combat this problem is to introduce radical scavengers like cerium (Ce) ions that neutralize these species before they attack the PFSA. Such cation doping creates a multi-ion system, in which understanding the mechanisms of cation solvation and transport becomes important for the effective design and utilization of PFSA–cation systems. Ce ions also provide a representative model system for multication-exchanged ionomers in electrochemical systems. In this study, hydration and conductivity measurements, along with X-ray fluorescence and scattering, are employed to elucidate how Ce ion exchange alters PFSA’s ionic solvation, as well as nano- and mesoscale morphologies, which ultimately control its ion transport properties. A molecular transport model is used to deconvolute the impact of Ce ions on the local solvation structure of water in the membrane from mesoscale changes of the transport pathways. The combined experimental and theoretical analysis reveals a nonlinear decrease in conductivity driven by cation solvation at the molecular level and morphological changes at longer length scales. Migration–diffusion coupling, its nonlinear dependence on ion exchange and hydration, and its overall implications for ionomer performance are also discussed. Finally, these findings have the potential to be translated into other mixed cation–ionomer systems for a wide range of energy and environmental devices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Lyapunov-Based Generalized Dc-Side Controller Design for PV-Connected Systems

The objective of this paper is to realize a universal dc-side controller for photovoltaic (PV) systems where the control is agnostic to the downstream converter configuration. To achieve this, the downstream power converter and its controls are manipulated into an effective power control loop that is then cast into a generalized multi-loop design framework. On the dc side, a nonlinear small-signal model of the PV input is realized exclusively in terms of PV datasheet parameters (i.e., open-circuit voltage, short-circuit current, and maximum power point). Finally, a linear controller is used to modulate the dc-side PV system with the generic downstream power controller. A Lyapunov candidate is proposed to analyze the stability of the interconnected system and provide a streamlined approach for the controller design. The proposed design is validated on a 1k VA experimental setup that interfaces a PV module to the grid.

generalized PV controller↗
DOE OSTI2022

A Lyapunov-Based Generalized Dc-Side Controller Design for PV-Connected Systems: Preprint

The objective of this paper is to realize a universal dc-side controller for photovoltaic (PV) systems where the control is agnostic to the downstream converter configuration. To achieve this, the downstream power converter and its controls are manipulated into an effective power control loop that is then cast into a generalized multi-loop design framework. On the dc side, a nonlinear small-signal model of the PV input is realized exclusively in terms of PV datasheet parameters (i.e., open-circuit voltage, short-circuit current, and maximum power point). Finally, a linear controller is used to modulate the dc-side PV system with the generic downstream power controller. A Lyapunov candidate is proposed to analyze the stability of the interconnected system and provide a streamlined approach for the controller design. The proposed design is validated on a 1-kVA experimental setup that interfaces a PV module to the grid.

dc-dc control↗
DOE OSTI2021

Tunable plasmon-enhanced second-order optical nonlinearity in transition metal dichalcogenide nanotriangles

The development of nanomaterials with a large nonlinear susceptibility is essential for nonlinear nanophotonics. Here we show that transition metal dichalcogenide (TMD) nanotriangles have a large effective second-order susceptibility [χ (2) ] at midinfrared to near-infrared frequencies owing to their broken centrosymmetry. χ (2) is calculated within the density-matrix formalism that accounts for dissipation and screening. χ (2) peaks in the vicinity of both two-photon resonances (specified by the geometry) and plasmon resonances (tunable via the carrier density). Aligning the resonances yields the values of χ (2) as high as 10 –6 m/V. These findings underscore the potential of TMD nanotriangles for nonlinear nanophotonics, particularly second-harmonic generation.

0-dimensional systems↗
DOE OSTI2023

Multigrid Reduction in Time for Chaotic Dynamical Systems

As CPU clock speeds have stagnated and high performance computers continue to have ever higher core counts, increased parallelism is needed to take advantage of these new architectures. Traditional serial time-marching schemes can be a significant bottleneck, as many types of simulations require large numbers of time-steps which must be computed sequentially. Parallel-in-time schemes, such as the Multigrid Reduction in Time (MGRIT) method, remedy this by parallelizing across time-steps and have shown promising results for parabolic problems. However, chaotic problems have proved more difficult, since chaotic initial value problems (IVPs) are inherently ill-conditioned. MGRIT relies on a hierarchy of successively coarser time-grids to iteratively correct the solution on the finest time-grid, but due to the nature of chaotic systems, small inaccuracies on the coarser levels can be greatly magnified and lead to poor coarse-grid corrections. Here we introduce a modified MGRIT algorithm based on an existing quadratically converging nonlinear extension to the multigrid Full Approximation Scheme (FAS), as well as a novel time-coarsening scheme. Together, these approaches better capture long-term chaotic behavior on coarse-grids and greatly improve convergence of MGRIT for chaotic IVPs. Further, we introduce a novel low-memory variant of the algorithm for solving chaotic PDEs with MGRIT which not only solves the IVP, but also provides estimates for the unstable Lyapunov vectors of the system. Finally, we provide supporting numerical results for the Lorenz system and demonstrate parallel speedup for the chaotic Kuramoto–Sivashinsky PDE over a significantly longer time-domain than in previous works.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Unitary Qubit Lattice Algorithms for Plasma Physics

This final technical report summarizes research conducted under DOE Award DE-SC0021653 to develop unitary Quantum Lattice Algorithms for modeling electromagnetic wave propagation and scattering in complex media, including plasmas. The project developed and validated quantum-inspired formulations of Maxwell's equations that preserve unitary evolution and can be evaluated on classical high-performance computing systems while providing a foundation for future quantum-computing implementations. Major accomplishments include the development of two- and three-dimensional algorithms for electromagnetic scattering; scalable, distributed-memory implementations demonstrated on the Perlmutter supercomputer; formulations for nonlinear lossless fluid dynamics and cold, lossless, inhomogeneous magnetized plasmas; and an explicit quantum algorithm for a time-discretized Lorenz model. Simulations reproduced a range of characteristic wave phenomena, including transient effects that are not readily apparent in conventional frequency-domain studies, demonstrating the effectiveness of the proposed approach for modeling complex electromagnetic and plasma systems. The work establishes a unified theoretical and computational framework for quantum and quantum-inspired simulation and provides a foundation for future implementation on fault-tolerant quantum systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Hybrid large-eddy simulation with adaptive downstream anisotropic eddy viscosity model

Turbulence pervades most flows of engineering interest, and its prediction remains a challenge on both accuracy and cost. One promising predictive approach that reduces cost combines large-eddy simulation (LES) with simulation based on Reynolds averaged Navier–Stokes equations (RANS). This study presents a method to overcome stability and accuracy issues associated with these hybrid LES-RANS methods. The method developed involves extracting the Reynolds stresses from the upstream LES solution and then using that information to convert the downstream RANS model from a scalar eddy-viscosity model to an anisotropic nonlinear eddy-viscosity model. The method developed differs from the downstream tensorial eddy-viscosity model by being independent of the coordinate system. The method developed was evaluated by computing film cooling of a flat plate with the coolant injected through one row of circular holes. Results obtained show instabilities at the LES-to-RANS interface to be eliminated. Also, the method developed yielded solutions that compare reasonably well with those from LES, even though a significant portion of the flow is computed by the adapted anisotropic RANS model instead of LES, which significantly reduced the number of grid points and computational time needed. Since modification of the downstream RANS model is based on information extracted from the upstream LES solution, the method developed is adaptive to the problem being studied.

42 ENGINEERING↗
DOE OSTI2023

Ultra‐Low‐Loss Silicon Nitride Photonics Based on Deposited Films Compatible with Foundries

Abstract The fabrication processes of silicon nitride (Si 3 N 4 ) photonic devices used in foundries require low temperature deposition, which typically leads to high propagation losses. Here, it is shown that propagation loss as low as 0.42 dB cm −1 can be achieved using foundry compatible processes by solely reducing waveguide surface roughness. By postprocessing the fabricated devices using rapid thermal anneal (RTA) and furnace anneal, propagation losses down to 0.28 dB cm −1 and 0.06 dB cm −1 , respectively, are achieved. These low losses are comparable to the conventional devices using high temperature, high‐stress LPCVD films. The dispersion of the devices is also tuned, and it is proved that these devices can be used for linear and nonlinear applications. Low threshold parametric oscillation, broadband frequency combs, and narrow‐linewidth laser are demonstrated. This work demonstrates the feasibility of scalable photonic systems based on foundries.

Ji, Xingchen↗
DOE OSTI2025

Density Functional Tight-Binding Enables Tractable Studies of Quantum Plasmonics

Routine investigations of plasmonic phenomena at the quantum level present a formidable computational challenge due to the large system sizes and ultrafast time scales involved. This Feature Article highlights the use of density functional tight-binding (DFTB), particularly its real-time time-dependent formulation (RT-TDDFTB), as a tractable approach to study plasmonic nanostructures from a quantum mechanical purview. We begin by outlining the theoretical framework and limitations of DFTB, emphasizing its efficiency in modeling systems with thousands of atoms over picosecond time scales. Applications of RT-TDDFTB are then explored in the context of optical absorption, nonlinear harmonic generation, and plasmon-mediated photocatalysis. We demonstrate how DFTB can reconcile classical and quantum descriptions of plasmonic behavior, capturing key phenomena such as size-dependent plasmon shifts and plasmon coupling in nanoparticle assemblies. Lastly, we showcase DFTB’s ability to model hot carrier generation and reaction dynamics in plasmon-driven H2 dissociation, underscoring its potential to model photocatalytic processes. Collectively, these studies establish DFTB as a powerful, yet computationally efficient tool to probe the emergent physics of materials at the limits of space and time.

Electrical energy↗
DOE OSTI2021

Quantum Machine-Learning for Eigenstate Filtration in Two-Dimensional Materials

Quantum machine-learning algorithms have emerged to be a promising alternative to their classical counterparts as they leverage the power of quantum computers. Such algorithms have been developed to solve problems like electronic structure calculations of molecular systems and spin models in magnetic systems. However, the discussion in all these recipes focuses specifically on targeting the ground state. Herein we demonstrate a quantum algorithm that can filter any energy eigenstate of the system based on either symmetry properties or a predefined choice of the user. The workhorse of our technique is a shallow neural network encoding the desired state of the system with the amplitude computed by sampling the Gibbs–Boltzmann distribution using a quantum circuit and the phase information obtained classically from the nonlinear activation of a separate set of neurons. We show that the resource requirements of our algorithm are strictly quadratic. To demonstrate its efficacy, we use state filtration in monolayer transition metal dichalcogenides which are hitherto unexplored in any flavor of quantum simulations. We implement our algorithm not only on quantum simulators but also on actual IBM-Q quantum devices and show good agreement with the results procured from conventional electronic structure calculations. We thus expect our protocol to provide a new alternative in exploring the band structures of exquisite materials to usual electronic structure methods or machine-learning techniques that are implementable solely on a classical computer.

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

Nonlinear multi-magnon scattering in artificial spin ice

Abstract Magnons, the quantum-mechanical fundamental excitations of magnetic solids, are bosons whose number does not need to be conserved in scattering processes. Microwave-induced parametric magnon processes, often called Suhl instabilities, have been believed to occur in magnetic thin films only, where quasi-continuous magnon bands exist. Here, we reveal the existence of such nonlinear magnon-magnon scattering processes and their coherence in ensembles of magnetic nanostructures known as artificial spin ice. We find that these systems exhibit effective scattering processes akin to those observed in continuous magnetic thin films. We utilize a combined microwave and microfocused Brillouin light scattering measurement approach to investigate the evolution of their modes. Scattering events occur between resonance frequencies that are determined by each nanomagnet’s mode volume and profile. Comparison with numerical simulations reveals that frequency doubling is enabled by exciting a subset of nanomagnets that, in turn, act as nanosized antennas, an effect that is akin to scattering in continuous films. Moreover, our results suggest that tunable directional scattering is possible in these structures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Real-Time Movement-Based Traffic Volume Prediction at Signalized Intersections

The traffic volume of each movement at signalized intersections can provide valuable information on real-time traffic conditions that enable traffic control systems to dynamically respond to the fluctuated traffic demands. Real-time movement-based traffic volume prediction is challenging due to various nonlinear spatial relationships at different locations/approaches and the complicated underlying temporal dependencies. In this study, a novel deep intersection spatial-temporal network (DISTN) is developed for real-time movement-based traffic volume prediction at signalized intersections, which considers both spatial and temporal features by the convolutional neural network (CNN) and long short-term memory (LSTM), respectively. In addition, the within-day, daily, and weekly periodic trends of traffic volume are also considered in the proposed model. This is the first time that a deep-learning method has been applied for movement-based traffic volume prediction at signalized intersections. In the numerical experiment, the proposed model is evaluated using real-world data and simulation data to demonstrate its effectiveness. The impacts of various structures of traffic networks on the proposed model are also discussed. Results show that the proposed model outperforms some of the state-of-the-art volume prediction methods currently in the literature.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Shock tube experiments on the three-layer Richtmyer–Meshkov instability

Here, a vertical shock tube is used for experiments on the three-layer Richtmyer–Meshkov instability. Two closely spaced membrane-less interfaces are formed by the flow of two different sects of three gases: one with air above CO 2 above SF 6 and the other with helium above air above SF 6 . The lightest of the three gases enters the shock tube at the top of the driven section and flows downward. Conversely, the heaviest gas enters at the bottom of the shock tube and flows upward while the intermediate density gas enters at the middle through porous plates. All three gases are allowed to escape through holes at the layer location, leaving an approximately 30-mm layer of intermediate-density gas suspended between the lightest gas from above and the heaviest gas from below. A single-mode, two-dimensional initial perturbation is then imposed on the lower interface by oscillating the shock tube in the horizontal direction. The flow is visualized by seeding the intermediate gas with particles and illuminating it with a pulsed laser. Image sequences are then captured using high-speed video cameras. Perturbation amplitude measurements are made from the three-layer system and compared with measurements from 2, two-layer systems. It is observed that the presence of the upper, initially flat interface produces a decrease in growth of instability amplitude in the nonlinear phase over an equivalent single-interface configuration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Static actuator-sharing algorithm for concurrent control of multiple plasma properties

Simultaneous regulation of multiple properties in next-generation tokamaks like ITER and fusion pilot plant may require the integration of different plasma control algorithms. Such integration requires the conversion of individual controller commands into physical actuator requests while accounting for the coupling between different plasma properties. This work proposes a tokamak and scenario-agnostic actuator-sharing algorithm (ASA) to perform the above-mentioned command-request conversion and, hence, integrate multiple plasma controllers. The proposed algorithm implicitly solves a quadratic programming (QP) problem formulated to account for the saturation limits and the relation between the controller commands and physical actuator requests. Since the constraints arising in the QP program are linear, the proposed ASA is highly computationally efficient and can be implemented in the tokamak plasma control system in real time. Furthermore, the proposed algorithm is designed to handle real-time changes in the control objectives and actuators’ availability. Nonlinear simulations carried out using the Control Oriented Transport SIMulator illustrate the effectiveness of the proposed algorithm in achieving multiple control objectives simultaneously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Transport, correlations, and chaos in a classical disordered anharmonic chain

Here, we explore transport properties in a disordered nonlinear chain of classical harmonic oscillators, and thereby identify a regime exhibiting behavior analogous to that seen in quantum many-body-localized systems. Through extensive numerical simulations of this system connected at its ends to heat baths at different temperatures, we computed the heat current and the temperature profile in the nonequilibrium steady state as a function of system size $\textit{N}$, disorder strength Δ, and temperature $\textit{T}$. The conductivity $κ_N$, obtained for finite length ($\textit{N}$), saturates to a value $κ_∞$ > 0 in the large $\textit{N}$ limit, for all values of disorder strength Δ and temperature $\textit{T}$ > 0. We show evidence that for any Δ > 0 the conductivity goes to zero faster than any power of $\textit{T}$ in the ($\textit{T}$/Δ) → 0 limit, and find that the form $κ_∞$ ~ $e^{{–B|\text{ln}(CΔ/T)|}^{3}}$ fits our data. This form has earlier been suggested by a theory based on the dynamics of multioscillator chaotic islands. The finite-size effect can be $κ_N < κ_∞$ due to boundary resistance when the bulk conductivity is high (the weak disorder case), or $κ_N > κ_∞$ due to direct bath-to-bath coupling through bulk localized modes when the bulk is weakly conducting (the strong disorder case). We also present results on equilibrium dynamical correlation functions and on the role of chaos on transport properties. Finally, we explore the differences in the growth and propagation of chaos in the weak and strong chaos regimes by studying the classical version of the out-of-time-ordered commutator.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Interactions Between Climate Mean and Variability Drive Future Agroecosystem Vulnerability

ABSTRACT Agriculture is crucial for global food supply and dominates the Earth's land surface. It is unknown, however, how slow but relentless changes in climate mean state, versus random extreme conditions arising from changing variability , will affect agroecosystems' carbon fluxes, energy fluxes, and crop production. We used an advanced weather generator to partition changes in mean climate state versus variability for both temperature and precipitation, producing forcing data to drive factorial‐design simulations of US Midwest agricultural regions in the Energy Exascale Earth System Model. We found that an increase in temperature mean lowers stored carbon, plant productivity, and crop yield, and tends to convert agroecosystems from a carbon sink to a source, as expected; it also can cause local to regional cooling in the earth system model through its effects on the Bowen Ratio. The combined effect of mean and variability changes on carbon fluxes and pools was nonlinear, that is, greater than each individual case. For instance, gross primary production reduces by 9%, 1%, and 13% due to change in mean temperature, change in temperature variability, and change in both temperature mean and variability, respectively. Overall, the scenario with change in both temperature and precipitation means leads to the largest reduction in carbon fluxes (−16% gross primary production), carbon pools (−35% vegetation carbon), and crop yields (−33% and −22% median reduction in yield for corn and soybean, respectively). By unambiguously parsing the effects of changing climate mean versus variability and quantifying their nonadditive impacts, this study lays a foundation for more robust understanding and prediction of agroecosystems' vulnerability to 21st‐century climate change.

54 ENVIRONMENTAL SCIENCES↗