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

Scalable generalized screening for high-order terms in the many-body expansion: Algorithm, open-source implementation, and demonstration

The many-body expansion lies at the heart of numerous fragment-based methods that are intended to sidestep the nonlinear scaling of ab initio quantum chemistry, making electronic structure calculations feasible in large systems. In principle, inclusion of higher-order n-body terms ought to improve the accuracy in a controllable way, but unfavorable combinatorics often defeats this in practice and applications with n ≥ 4 are rare. Here, we outline an algorithm to overcome this combinatorial bottleneck, based on a bottom-up approach to energy-based screening. This is implemented within a new open-source software application (“Fragme∩t”), which is integrated with a lightweight semi-empirical method that is used to cull subsystems, attenuating the combinatorial growth of higher-order terms in the graph that is used to manage the calculations. This facilitates applications of unprecedented size, and we report four-body calculations in (H2O)64 clusters that afford relative energies within 0.1 kcal/mol/monomer of the supersystem result using less than 10% of the unique subsystems. We also report n-body calculations in (H2O)20 clusters up to n = 8, at which point the expansion terminates naturally due to screening. These are the largest n-body calculations reported to date using ab initio electronic structure theory, and they confirm that high-order n-body terms are mostly artifacts of basis-set superposition error.

Chemistry↗
DOE OSTI2021

McMillan electron lens in a system with space charge

Space charge (SC) forces of a circulating beam in a ring have both linear (defocusing) and nonlinear components, due to a nonuniform beam distribution. The linear component of SC forces produces a betatron tune shift, which is the largest for a zero-amplitude particle, while the nonlinear component produces an amplitude-dependent betatron tune spread. These SC effects are responsible for several undesirable phenomena in accelerators: emittance growth, particle losses, beam halo, etc. In this paper, we investigate the possibility to mitigate the distributed SC forces by a thin McMillan lens, providing an axially-symmetric kick, which is qualitatively opposite to the accumulated effect of beam's own SC. Experimentally, the proposed concept can be tested in Fermilab's IOTA ring. A thin McMillan lens can be implemented by a short (70) insertion of an electron beam with a specifically chosen density distribution in transverse directions. In this article, to test if McMillan lenses can reduce the tune spread induced by SC, we make several simulations with a 6-D particle tracking code, Synergia. We choose such beam and lattice parameters that the SC tune spread is roughly 0.5 and the emittance growth due to the half-integer resonance is clearly observed without the SC compensation. Then, we focus on reducing the emittance growth by adjusting the bare betatron tunes using the ring quadrupoles, and reducing the tune spread by the McMillan lenses. The results of reducing a large tune spread (≈ 0.5), reported here, are not perfect, but substantial. There is still room for further investigation. The simulations performed so far in this work indicate that McMillan lenses can cope with an SC tune spread of ≈ 0.1 per lens.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

NN-OpInf

SAND2026-18878O The NN-OpInf tool is a PyTorch-based approach to operator inference that uses composable, structure-preserving neural networks to represent nonlinear operators. Operator inference is a machine learning method for inferring low-dimensional systems from data and polynomial models for system dynamics. However, many systems do not conform to polynomial structures, which NN-OpInf addresses by parameterizing operators with neural networks. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗
DOE OSTI2021

Closed-Form Solutions for the Equations of Motion of the Heavy Symmetrical Top with One Point Fixed

The equations of motion (EOM) for the heavy symmetrical top with one point fixed are highly nonlinear. The literature describes the numerical methods that are used to resolve this classical system, including modern tools, such as the Runge–Kutta fourth–order method. Finding the derivate of closed-form solutions for the EOM is more difficult and, as mentioned in the literature, discovering the solution is not always possible for all the EOM. Fortunately, a few examples are available that serve as a guide to move further in this topic. The purpose of this paper is to find a methodology that will produce the solutions for a given subset of EOMs that fulfill certain requisites. This paper summarizes the literature available on this topic and then follows with the derivation of the EOM using the Euler–Lagrange method. The Routhian method will be used to reduce the size of the expression, and it continues with the formulation of the classical cubic function, ƒ(u), through a novel process. The roots of ƒ(u) are of the utmost importance in finding the EOM closed-form solution, and once the final roots are selected, the general method that will produce the closed-form solutions is presented. Two sets of examples are included to show the validity of the process, and comparisons of the results from the closed-form solutions vs. the numerical results for these examples are shown.

Laos, Hector↗
DOE OSTI2022

Optimizing molecular light absorption in the strong coupling regime for solar energy harvesting

The strong coupling of optical absorbers (e.g., molecules or semiconductors) to confined photonic modes fundamentally alters the physical properties of the coupled system via the formation of hybrid light-matter states. One potential application of strong light-matter coupling relies on exploiting it to localize light-induced charge excitation processes to small volumes of material. Applications that would benefit from this localization include thin-film photovoltaics, photodetection, photocatalysis, and others, where the overall performance depends on the ability of a material to amplify light absorption (i.e., the formation of electron-hole pairs) within specific locations in space. This contribution investigates how strong light-matter coupling affects light absorption rates in molecular absorbers coupled to photonic nanostructures. Here, our results show that the molecular light absorption efficiencies are highest in configurations where the strongly coupled molecules interact directly with the incoming photon flux. We also identify a nonlinear dependence in the molecular absorption as a function of concentration, unique to the strongly coupled systems. Based on these results, we propose design principles for engineering nanostructured systems that allow for high efficiencies of charge carrier localization into strongly coupled absorbers.

14 SOLAR ENERGY↗
DOE OSTI2024

Direct Observation of Circularly Polarized Nonlinear Optical Activities in Chiral Hybrid Lead Halides

Circularly polarized light emission is a crucial application in imaging, sensing, and photonics. However, utilizing low-energy photons to excite materials, as opposed to high-energy light excitation, can facilitate deep-tissue imaging and sensing applications. The challenge lies in finding materials capable of directly generating circularly polarized nonlinear optical effects. In this study, we introduce a chiral hybrid lead halide (CHLH) material system, R/S-DPEDPb 3 Br 8 ·H 2 O (DPED = 1,2-diphenylethylenediammonium), which can directly produce circularly polarized second harmonic generation (CP-SHG) through linearly polarized infrared light excitation, exhibiting a polarization efficiency as high as 37% at room temperature. To understand the spin relaxation mechanisms behind the high polarization efficiency, we utilized two models, so-called D’yakonov–Perel’ (DP) and Bir–Aronov–Pikus (BAP) mechanisms. The unique zigzag inorganic frameworks within the hybrid structure are believed to reduce the dielectric confinement and exciton binding energy, thus enhancing spin polarization, especially in regions with a high excitation pump fluence based on the DP mechanism. In the case of low excitation pump fluence, the BAP mechanism dominates, as evidenced by the observed decrease in the polarization ratio from CP-SHG measurement. Using density functional theory analysis, we elucidate how the distinctive 8-coordination environment of lead bromide building blocks effectively suppresses spin–orbit coupling at the conduction band minimum. This suppression significantly diminishes spin-splitting, thereby slowing the spin relaxation rate.

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

Equations of Motion for the Vertical Rigid-Body Rotor: Linear and Nonlinear Cases

Centuries ago, the prolific mathematician Leonhard Euler (1707–1783) wrote down the equations of motion (EOM) for the heavy symmetrical top with one point fixed. The resulting set of equations turned out to be nonlinear and had a limited number of closed-form solutions.Today, tools such as transfer matrix and finite elements enable the calculation of the rotordynamic properties for rotor-bearing systems. Some of these tools rely on the “linearized” version of the EOM to calculate the eigenvalues, unbalance response, or transients in these systems.In fact, industry standards mandate that rotors be precisely balanced to have safe operational characteristics. However, in some cases, the nonlinear aspect of the EOM should be considered.The purpose of this chapter is to show examples of how the linear vs. nonlinear formulations differ. This chapter also shows how excessive unbalance is capable of dramatically altering the behavior of the system and can produce chaotic motions associated with the “jump” phenomenon.

Laos, Hector↗
DOE OSTI2024

MHD-induced beta limits in the Large Helical Device

Using the extended-magnetohydrodynamics code, M3D-C1, we perform a systematic numerical study of the effect of externally applied heating on the achievable plasma beta in a ten field-period heliotron. Heat sources of varying intensity are applied to a vacuum magnetic field that is representative of the standard configuration of the Large Helical Device, with R 0 = 3.66 m, where R 0 is the radial position of the magnetic axis in vacuum. As the system is driven to a state that is unstable to low-n magnetohydrodynamic (MHD) modes, nonlinear mode interactions lead to the formation of chaotic magnetic fields. With sufficiently strong heating, a collapse of the electron temperature profile is observed. This demonstrates the necessity of simulating the self-consistent evolution of plasma profiles, without imposing assumptions on the structure of the magnetic field, to accurately determine transport properties in stellarator plasmas. It also highlights the value of these advanced simulation capabilities for accelerating the development of high-performance stellarator operating scenarios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Time-lens photon Doppler velocimetry (TL-PDV)

Here, we describe a time lens (TL) to expand the dynamic range of photon Doppler velocimetry (PDV) systems. The principle and preliminary design of a TL-PDV system are explained and shown to be feasible through simulations. In a PDV system, an interferometer is used for measuring frequency shifts due to the Doppler effect from the target motion. However, the sampling rate of the electronics could limit the velocity range of a PDV system. A four-wave-mixing (FWM) TL applies a quadratic temporal phase to an optical signal within a nonlinear FWM medium (such as an integrated photonic waveguide or a highly nonlinear optical fiber). By spectrally isolating the mixing product, termed the idler, and with appropriate lengths of dispersion prior to and after this FWM TL, a temporally magnified version of the input signal is generated. Therefore, the frequency shifts of PDV can be “slowed down” with the magnification factor M of the TL. M = 1 corresponds to a regular PDV system without a TL. M = 10 has been shown to be feasible for a TL-PDV system. The use of this effect for PDV can expand the velocity measurement range and allow for the use of lower bandwidth electronics. TL-PDV will open up new avenues for various dynamic material experiments.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Coupling magnons to an opto-electronic parametric oscillator

Hybrid magnonic systems have emerged as versatile modular components for quantum signal transduction and sensing applications owing to their capability of connecting distinct quantum platforms. To date, the majority of the magnonic systems have been explored in a local, near-field scheme, due to the close proximity required for realizing a strong coupling between magnons and other excitations. This constraint greatly limits the applicability of magnons in developing remotely coupled, distributed quantum network systems. On the contrary, opto-electronic architectures hosting self-sustained oscillations have been a unique platform for long-haul signal transmission and processing. Here, we integrated an opto-electronic oscillator with a magnonic oscillator consisting of a microwave waveguide and a Y 3 Fe 5 O 12 (YIG) sphere, and demonstrated strong and coherent coupling between YIG’s magnon modes and the opto-electronic oscillator’s characteristic photon modes—revealing the hallmark anti-crossing gap in the measured spectrum. In particular, the photon mode is produced on-demand via a nonlinear, parametric process as stipulated by an external seed pump. Both the internal cavity phase and the external pump phase can be precisely tuned to stabilize either degenerate or nondegenerate auto-oscillations. Our result lays out a new, hybrid platform for investigating the long-distance coupling and nonlinearity in coherent magnonic phenomena, which may be found useful in constructing the future “distributed hybrid magnonic systems.”

36 MATERIALS SCIENCE↗
DOE OSTI2020

Pattern formation in a coupled membrane-bulk reaction-diffusion model for intracellular polarization and oscillations

Reaction-diffusion systems have been widely used to study spatio-temporal phenomena in cell biology, such as cell polarization. Coupled bulk-surface models naturally include compartmentalization of cytosolic and membrane-bound polarity molecules. Here we study the distribution of the polarity protein Cdc42 in a mass-conserved membrane-bulk model, and explore the effects of diffusion and spatial dimensionality on spatio-temporal pattern formation. We first analyze a one-dimensional (1-D) model for Cdc42 oscillations in fission yeast, consisting of two diffusion equations in the bulk domain coupled to nonlinear ODEs for binding kinetics at each end of the cell. In 1-D, our analysis reveals the existence of symmetric and asymmetric steady states, as well as anti-phase relaxation oscillations typical of slow-fast systems. We then extend our analysis to a two-dimensional (2-D) model with circular bulk geometry, for which species can either diffuse inside the cell or become bound to the membrane and undergo a nonlinear reaction-diffusion process. We also consider a nonlocal system of PDEs approximating the dynamics of the 2-D membrane-bulk model in the limit of fast bulk diffusion. In all three model variants we find that mass conservation selects perturbations of spatial modes that simply redistribute mass. In 1-D, only anti-phase oscillations between the two ends of the cell can occur, and in-phase oscillations are excluded. In higher dimensions, no radially symmetric oscillations are observed. Instead, the only instabilities are symmetry-breaking, either corresponding to stationary Turing instabilities, leading to the formation of stationary patterns, or to oscillatory Turing instabilities, leading to traveling and standing waves. Codimension-two Bogdanov—Takens bifurcations occur when the two distinct instabilities coincide, causing traveling waves to slow down and to eventually become stationary patterns. Our work clarifies the effect of geometry and dimensionality on behaviors observed in ma.ss-conserved cell polarity models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Charge Conservation beyond Uniformity: Spatially Inhomogeneous Electromagnetic Response in Periodic Solids

Nonlinear electromagnetic response functions have reemerged as a crucial tool for studying quantum materials, due to recently appreciated connections between optical response functions, quantum geometry, and band topology. Most attention has been paid to responses to spatially uniform electric fields, relevant to low-energy optical experiments in conventional solid state materials. However, magnetic and magnetoelectric phenomena are naturally connected by responses to spatially varying electric fields due to Maxwell’s equations. Furthermore, in the emerging field of moiré materials, characteristic lattice scales are much longer, allowing spatial variation of optical electric fields to potentially have a measurable effect in experiments. In order to address these issues, we develop a formalism for computing linear and nonlinear responses to spatially inhomogeneous electromagnetic fields. Starting with the continuity equation, we derive an expression for the second-quantized current operator that is manifestly conserved and model independent. Crucially, our formalism makes no assumptions on the form of the microscopic Hamiltonian and so is applicable to model Hamiltonians derived from tight-binding or calculations. We then develop a diagrammatic Kubo formalism for computing the wave vector dependence of linear and nonlinear conductivities, using Ward identities to fix the value of the diamagnetic current order by order in the vector potential. We apply our formula to compute the magnitude of the Kerr effect at oblique incidence for a model of a moiré-Chern insulator and demonstrate the experimental relevance of spatially inhomogeneous fields in these systems. We further show how our formalism allows us to compute the (orbital) magnetic multipole moments and magnetic susceptibilities in insulators. Turning to nonlinear response, we use our formalism to compute the second-order transverse response to spatially varying transverse electric fields in our moiré-Chern insulator model, with an eye toward the next generation of experiments in these systems. Published by the American Physical Society 2024

Physics↗
DOE OSTI2022

HomPINNs: Homotopy physics-informed neural networks for learning multiple solutions of nonlinear elliptic differential equations

Physics-informed neural networks (PINNs) based machine learning is an emerging framework for solving nonlinear differential equations. However, due to the implicit regularity of neural network structure, PINNs can only find the flattest solution in most cases by minimizing the loss functions. In this paper, we combine PINNs with the homotopy continuation method, a classical numerical method to compute isolated roots of polynomial systems, and propose a new deep learning framework, named homotopy physics-informed neural networks (HomPINNs), for solving multiple solutions of nonlinear elliptic differential equations. The implementation of an HomPINN is a homotopy process that is composed of the training of a fully connected neural network, named the starting neural network, and training processes of several PINNs with different tracking parameters. The starting neural network is to approximate a starting function constructed by the trivial solutions, while other PINNs are to minimize the loss functions defined by boundary condition and homotopy functions, varying with different tracking parameters. These training processes are regraded as different steps of a homotopy process, and a PINN is initialized by the well-trained neural network of the previous step, while the first starting neural network is initialized using the default initialization method. Finally, several numerical examples are presented to show the efficiency of our proposed HomPINNs, including reaction-diffusion equations with a heart-shaped domain.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Preparation for quantum simulation of the ( 1 + 1 ) -dimensional O(3) nonlinear σ model using cold atoms

The 1+1D O(3) non-linear σ-model is a model system for future quantum lattice simulations of other asymptotically-free theories, such as non-Abelian gauge theories. We find that utilizing dimensional reduction can make efficient use of two-dimensional layouts presently available on cold atom quantum simulators. A new definition of the renormalized coupling is introduced, which is applicable to systems with open boundary conditions and can be measured using analog quantum simulators. Monte Carlo and tensor network calculations are performed to determine the quantum resources required to reproduce perturbative short-distance observables. In particular, we show that a rectangular array of 48 Rydberg atoms with existing quantum hardware capabilities should be able to adiabatically prepare low-energy states of the perturbatively-matched theory. Furthermore, these states can then be used to simulate non-perturbative observables in the continuum limit that lie beyond the reach of classical computers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Evaluating route to impact convergence of the harmonic balance method for piecewise-smooth systems

Here in this work, we investigate the applicability of the harmonic balance method (HBM) to predict periodic solutions of a single degree-of-freedom forced Duffing oscillator with freeplay nonlinearity. By studying the route to impact, which refers to a parametric study as the contact stiffness increases from soft to hard, the convergence behavior of the HBM can be understood in terms of the strength of the non-smooth forcing term. HBM results are compared to time-integration results to facilitate an evaluation of the accuracy of nonlinear periodic responses. An additional contribution of this study is to perform convergence and stability analysis specifically for isolas generated by the non-smooth nonlinearity. Residual error analysis is used to determine the approximate number of harmonics required to get results accurate to a given error tolerance. Hill’s method and Floquet theory are employed to compute the stability of periodic solutions and identify the types of bifurcations in the system.

42 ENGINEERING↗
DOE OSTI2021

Up‐And‐Coming Advances in Optical and Microwave Nonreciprocity: From Classical to Quantum Realm

Reciprocity is a fundamental physical principle that roots in the time‐reversal symmetry of physical laws. It allows making predictions on any arbitrary complex system's response and operation and hence simplifies the analysis. However, there are many practical situations in which it is advantageous to break reciprocity, e.g., isolators preventing wave scattering back to lasers and generators, full‐duplex systems for multiplexing transmission and receiving in the same channel, nonreciprocal cavity excitation, and protection of fragile states of superconductor quantum computers from thermal noise. The most widespread approach to time‐reversal symmetry breaking and nonreciprocity based on magnetic field biasing suffers from bulkiness, cost ineffectiveness, and loss, motivating researchers and engineers to search for more practical approaches. Herein, the up‐and‐coming advances in optical nonreciprocity, including new materials (Weyl semimetals, topological insulators, metasurfaces), active structures, time‐modulation, parity‐time (PT)‐symmetry breaking, nonlinearity combined with a structural asymmetry, quantum nonlinearity, unidirectional gain and loss, chiral quantum states and valley polarization are overviewed. A general description of nonreciprocal systems is provided and the pros and cons of the mentioned approaches to nonreciprocity are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

When physics-informed data analytics outperforms black-box machine learning: A case study in thickness control for additive manufacturing

Aerosol jet printing (AJP) has emerged as a promising noncontact additive manufacturing method for high-resolution printing for a wide range of material systems. A key challenge limiting the broader adoption of AJP in the material science community is the lack of methods to precisely control thickness. Herein, we develop a model-based design of experiment (MBDoE) framework that integrates physics-informed models, nonlinear regression, and information criteria to postulate, select and calibrate the best model to describe and optimize the AJP manufacturing process. Starting with already available data from system commissioning (e.g., prior single variable sensitivity analysis), four candidate physics-informed models are postulated and trained. MBDoE identifies a single additional optimal experiment to validate these predictive models with quantified uncertainties, which are then used to determine the best experimental conditions to control printed film thickness. As a comparative benchmark, the analysis is repeated using the same dataset with nonparametric Gaussian process regression (GPR) model that does not incorporate physical information. Using MBDoE principles, we find that only five experiments are necessary to calibrate the nonlinear physics-informed parametric model, and with said limited data, this model outperforms the black-box machine learning GPR model. This key result underscores an emerging trend in the data science community: incorporating physical information into predictive models often drastically reduces the data requirements. Leveraging MBDoE further increased the data efficiency. By design, the proposed data science framework is general in nature and can be easily extended to other experimental and additive manufacturing systems beyond AJP.

Aerosol jet printing↗