Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXACT SOLUTIONS”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Fibonacci fast convergence for neutrino oscillations in matter

Understanding neutrino oscillations in matter requires a non-trivial diagonalization of the Hamiltonian. As the exact solution is very complicated, many approximation schemes have been pursued. Here we show that one scheme, systematically applying rotations to change to a better basis, converges exponentially fast wherein the rate of convergence follows the Fibonacci sequence. We find that the convergence rate of this procedure depends very sensitively on the initial choices of the rotations as well as the mechanism of selecting the pivots. We then apply this scheme for neutrino oscillations in matter and discover that the optimal convergence rate is found using the following simple strategy: first apply the vacuum (2-3) rotation and then use the largest off-diagonal element as the pivot for each of the following rotations. The Fibonacci convergence rate presented here may be extendable to systems beyond neutrino oscillations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Delivery drone route planning over a battery swapping network

Many enterprises invest on drone delivery research and development to drop off packages at consumers’ doorsteps in a matter of minutes. We study delivery drone route planning over a battery swapping network allowing farther reach by penetrating current battery capacity constraints. A mixed-integer nonlinear programming model is created to plan efficient drone routing over the swapping machines by minimizing the delivery lead time. We develop an exact solution method, evaluate its performance, and compare it with a straightforward nonlinear solver application. A case study highlights the applicability of the model. Data and source code to the solver are publicly shared.

drone battery swapping↗
DOE OSTI2024

Richtmyer–Meshkov instability when a shock is reflected for fluids with arbitrary equation of state

First predicted by Richtmyer in 1960 and experimentally confirmed by Meshkov in 1969, the Richtmyer–Meshkov instability (RMI) is crucial in fields such as physics, astrophysics, inertial confinement fusion and high-energy-density physics. These disciplines often deal with strong shocks moving through condensed materials or high-pressure plasmas that exhibit non-ideal equations of state (EoS), thus requiring theoretical models with realistic fluid EoS for accurate RMI simulations. Approximate formulae for asymptotic growth rates, like those proposed by Richtmyer, are helpful but rely on heuristic prescriptions for compressible materials. These prescriptions can sometimes approximate the RMI growth rate well, but their accuracy remains uncertain without exact solutions, as the fully compressible RMI growth rate is influenced by both vorticity deposited during shock refraction and multiple sonic wave refractions. This study advances previous work by presenting an analytic, fully compressible theory of RMI for reflected shocks with arbitrary EoS. It compares theoretical predictions with heuristic prescriptions using ideal gas, van der Waals gas and three-term constitutive equations for simple metals, the latter being analysed with detailed and simplified ideal-gas-like EoS. We additionally offer an alternative explicit approximate formula for the asymptotic growth rate. The comprehensive model also incorporates the effects of constant-amplitude acoustic waves at the interface, associated with the D'yakov–Kontorovich instability in shocks.

Napieralski, Mario (ORCID:0009000692344901)↗
DOE OSTI2021

Multiple fermion scattering in the weakly coupled spin-chain compound YbAlO 3

The Heisenberg antiferromagnetic spin-1/2 chain, originally introduced almost a century ago, is one of the best studied models in quantum mechanics due to its exact solution, but nevertheless it continues to present new discoveries. Its low-energy physics is described by the Tomonaga-Luttinger liquid of spinless fermions, similar to the conduction electrons in one-dimensional metals. In this work we investigate the Heisenberg spin-chain compound YbAlO3 and show that the weak interchain coupling causes Umklapp scattering between the left- and right-moving fermions and stabilizes an incommensurate spin-density wave order at q = 2k F under finite magnetic fields. These Umklapp processes open a route to multiple coherent scattering of fermions, which results in the formation of satellites at integer multiples of the incommensurate fundamental wavevector Q = nq. Our work provides surprising and profound insight into bandstructure control for emergent fermions in quantum materials, and shows how neutron diffraction can be applied to investigate the phenomenon of coherent multiple scattering in metals through the proxy of quantum magnetic systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Dynamic gain and frequency comb formation in exceptional-point lasers

Abstract Exceptional points (EPs)—singularities in the parameter space of non-Hermitian systems where two nearby eigenmodes coalesce—feature unique properties with applications such as sensitivity enhancement and chiral emission. Existing realizations of EP lasers operate with static populations in the gain medium. By analyzing the full-wave Maxwell–Bloch equations, here we show that in a laser operating sufficiently close to an EP, the nonlinear gain will spontaneously induce a multi-spectral multi-modal instability above a pump threshold, which initiates an oscillating population inversion and generates a frequency comb. The efficiency of comb generation is enhanced by both the spectral degeneracy and the spatial coalescence of modes near an EP. Such an “EP comb” has a widely tunable repetition rate, self-starts without external modulators or a continuous-wave pump, and can be realized with an ultra-compact footprint. We develop an exact solution of the Maxwell–Bloch equations with an oscillating inversion, describing all spatiotemporal properties of the EP comb as a limit cycle. We numerically illustrate this phenomenon in a 5-μm-long gain-loss coupled AlGaAs cavity and adjust the EP comb repetition rate from 20 to 27 GHz. This work provides a rigorous spatiotemporal description of the rich laser behaviors that arise from the interplay between the non-Hermiticity, nonlinearity, and dynamics of a gain medium.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Bose–Einstein condensation of a two-magnon bound state in a spin-1 triangular lattice

In ordered magnets, the elementary excitations are spin waves (magnons), which obey Bose–Einstein statistics. Similarly to Cooper pairs in superconductors, magnons can be paired into bound states under attractive interactions. The Zeeman coupling to a magnetic field is able to tune the particle density through a quantum critical point, beyond which a ‘hidden order’ is predicted to exist. Here, in this work, we report direct observation of the Bose–Einstein condensation of the two-magnon bound state in Na 2 BaNi(PO 4 ) 2 . Comprehensive thermodynamic measurements confirmed the two-dimensional Bose–Einstein condensation quantum critical point at the saturation field. Inelastic neutron scattering experiments were performed to establish the microscopic model. An exact solution revealed stable two-magnon bound states that were further confirmed by electron spin resonance and nuclear magnetic resonance experiments, demonstrating that the quantum critical point is due to the pair condensation, and the phase below the saturation field is likely the long-sought-after spin nematic phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

The transition from field emission to collisional space-charge limited current with nonzero initial velocity

Abstract Multiple electron emission mechanisms often contribute in electron devices, motivating theoretical studies characterizing the transitions between them. Previous studies unified thermionic and field emission, defined by the Richardson-Laue-Dushman (RLD) and Fowler–Nordheim (FN) equations, respectively, with the Child-Langmuir (CL) law for vacuum space-charge limited current (SCLC); another study unified FN and CL with the Mott-Gurney (MG) law for collisional SCLC. However, thermionic emission, which introduces a nonzero injection velocity, may also occur in gas, motivating this analysis to unify RLD, FN, CL, and MG. We exactly calculate the current density as a function of applied voltage over a range of injection velocity (i.e., temperature), mobility, and gap distance. This exact solution approaches RLD, FN, and generalized CL (GCL) and MG (GMG) for nonzero injection velocity under appropriate limits. For nonzero initial velocity, GMG approaches zero for sufficiently small applied voltage and mobility, making these gaps always space-charge limited by either GMG at low voltage or GCL at high voltage. The third-order nexus between FN, GMG, and GCL changes negligibly from the zero initial velocity calculation over ten orders of magnitude of applied voltage. These results provide a closed form solution for GMG and guidance on thermionic emission in a collisional gap.

42 ENGINEERING↗
DOE OSTI2020

Modeling of surface mechanical behaviors of soft elastic solids: theory and examples

Surfaces of soft solids can have significant surface stress, extensional modulus and bending stiffness. Previous theoretical studies have usually examined cases in which both the surface stress and bending stiffness are constant, assuming small deformation. In this work we consider a general formulation in which the surface can support large deformation and carry both surface stresses and surface bending moments. We show that the large deformation theory can be reduced to the classical linear theory (Shuttleworth equation). We obtain exact solutions for problems of an inflated cylindrical shell and bending of a plate with a finite thickness. Our analysis illustrates the different manners in which surface stiffening and surface bending stabilize these structures. We discuss how the complex surface constitutive behaviors affect the stress field of the bulk. Our calculation provides insights into effects of strain-dependent surface stress and surface bending in the large deformation regime, and can be used as a model to implement surface finite elements to study large deformation of complex structures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Stability analysis of the Eulerian–Lagrangian finite volume methods for nonlinear hyperbolic equations in one space dimension

In this paper, we construct a novel Eulerian–Lagrangian finite volume (ELFV) method for nonlinear scalar hyperbolic equations in one space dimension. It is well known that the exact solutions to such problems may contain shocks though the initial conditions are smooth, and direct numerical methods may suffer from restricted time step sizes. To relieve the restriction, we propose an ELFV method, where the space-time domain was separated by the partition lines originated from the cell interfaces whose slopes are obtained following the Rakine–Hugoniot junmp condition. Unfortunately, to avoid the intersection of the partition lines, the time step sizes are still limited. To fix this gap, we detect effective troubled cells (ETCs) and carefully design the influence region of each ETC, within which the partitioned space-time regions are merged together to form a new one. Then with the new partition of the space-time domain, we theoretically prove that the proposed first-order scheme with Euler forward time discretization is total-variation-diminishing and maximum-principle-preserving with at least twice larger time step constraints than the classical first order Eulerian method for Burgers’ equation. Numerical experiments verify the optimality of the designed time step sizes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Variational control forces for enhanced sampling of nonequilibrium molecular dynamics simulations

We introduce a variational algorithm to estimate the likelihood of a rare event within a nonequilibrium molecular dynamics simulation through the evaluation of an optimal control force. Optimization of a control force within a chosen basis is made possible by explicit forms for the gradients of a cost function in terms of the susceptibility of driven trajectories to changes in variational parameters. We consider probabilities of time-integrated dynamical observables as characterized by their large deviation functions and find that in many cases, the variational estimate is quantitatively accurate. Additionally, we provide expressions to exactly correct the variational estimate that can be evaluated directly. We benchmark this algorithm against the numerically exact solution of a model of a driven particle in a periodic potential, where the control force can be represented with a complete basis. Additionally, we then demonstrate the utility of the algorithm in a model of repulsive particles on a line, which undergo a dynamical phase transition, resulting in singular changes to the form of the optimal control force. In both systems, we find fast convergence and are able to evaluate large deviation functions with significant increases in statistical efficiency over alternative Monte Carlo approaches.

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

Driving a magnetized domain wall in an antiferromagnet by magnons

We theoretically study the interaction of magnons, quanta of spin waves, and a domain wall in a one-dimensional easy-axis antiferromagnet in the presence of an external magnetic field applied along the easy axis. To this end, we begin by obtaining the exact solution for spin waves in the background of a domain wall magnetized by an external field. The finite magnetization inside the domain wall is shown to give rise to reflection of magnons scattering off the domain wall, deviating from the well-known result of reflection-free magnons in the absence of a magnetic field. For practical applications of the predicted reflection of magnons, we show that the magnon reflection contributes to the thermally driven domain-wall motion. Our work leads us to envision that inducing a finite magnetization in antiferromagnetic solitons such as vortices and skyrmions can be used to engender phenomena that do not occur in the absence of magnetization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Quantum simulations employing connected moments expansions

Further advancement of quantum computing (QC) is contingent on enabling many-body models that avoid deep circuits and excessive use of CNOT gates. To this end, we develop a QC approach employing finite-order connected moment expansions (CMX) and affordable procedures for initial state preparation. We demonstrate the performance of our approach employing several quantum variants of CMX through the classical emulations on the H2 molecule potential energy surface and the Anderson model with a broad range of correlation strength. The results show that our approach is robust and flexible. Good agreements with exact solutions can be maintained even at the dissociation and strong correlation limits.

quantum computing, quantum chemistry, Electron cor↗
DOE OSTI2024

Collisional space-charge-limited current with monoenergetic velocity: From Child–Langmuir to Mott–Gurney

All theories coupling electron emission theories ultimately approach the space-charge-limited current (SCLC) in vacuum, given by the Child–Langmuir (CL) law, for sufficiently high voltage, or the Mott–Gurney (MG) law for finite electron mobility and high (but not infinite) voltage. These analyses demonstrate the presence of an SCLC regime that cannot be described by either CL or MG. Here, we derive an exact solution for SCLC for general electron mobility and nonzero velocity. We recover the traditional CL with nonzero initial velocity at high voltage. For low mobility (or infinite collision frequency), we derive corrections to the MG law that depend on the ratio of initial velocity to the product of collision frequency and gap distance or initial velocity to drift velocity for low and high voltage, respectively. Increasing collisionality decreases the correction to SCLC for nonzero velocity, indicating that these corrections are less important for low-mobility materials (e.g., solids) than high-mobility materials (e.g., air or vacuum). For a given gap distance (collision frequency), increasing the collision frequency (gap distance) increases the voltage necessary to make the gap appear more like vacuum. These results provide a generalized SCLC for both collisionality and initial velocity when assessing the transitions between electron emission mechanisms.

Physics↗
DOE OSTI2024

Quantum Drude oscillators coupled with Coulomb potential as an efficient model for bonded and non-covalent interactions in atomic dimers

The quantum Drude oscillator (QDO) model has been widely used as an efficient surrogate to describe the electric response properties of matter as well as long-range interactions in molecules and materials. Most commonly, QDOs are coupled within the dipole approximation so that the Hamiltonian can be exactly diagonalized, which forms the basis for the many-body dispersion method [Phys. Rev. Lett. 108, 236402 (2012)]. The dipole coupling is efficient and allows us to study non-covalent many-body effects in systems with thousands of atoms. However, there are two limitations: (i) the need to regularize the interaction at short distances with empirical damping functions and (ii) the lack of multipolar effects in the coupling potential. In this work, we convincingly address both limitations of the dipole-coupled QDO model by presenting a numerically exact solution of the Coulomb-coupled QDO model by means of quantum Monte Carlo methods. We calculate the potential-energy surfaces of homogeneous QDO dimers, analyzing their properties as a function of the three tunable parameters: frequency, reduced mass, and charge. We study the coupled-QDO model behavior at short distances and show how to parameterize this model to enable an effective description of chemical bonds, such as the covalent bond in the H2 molecule.

Chemistry↗
DOE OSTI2024

Quantum subspace expansion in the presence of hardware noise

Finding ground state energies on current quantum processing units (QPUs) using algorithms such as the variational quantum eigensolver (VQE) continues to pose challenges. Hardware noise severely affects both the expressivity and trainability of parameterized quantum circuits, limiting them to shallow depths in practice. Here, we demonstrate that both issues can be addressed by synergistically integrating VQE with a quantum subspace expansion, allowing for an optimal balance between quantum and classical computing capabilities and costs. We perform a systematic benchmark analysis of the iterative quantum-assisted eigensolver in the presence of hardware noise. We determine ground state energies of 1D and 2D mixed-field Ising spin models on noisy simulators and the IBM QPUs ibmq_quito (5 qubits) and ibmq_guadalupe (16 qubits). To maximize accuracy, we propose a suitable criterion to select the subspace basis vectors according to the trace of the noisy overlap matrix. Finally, we show how to systematically approach the exact solution by performing controlled quantum error mitigation based on probabilistic error reduction on the noisy backend fake_guadalupe.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Verification and validation of detonation-shock-dynamics relations for explosives described by general equation of state and chemical reaction models

Detonation shock dynamics is a powerful method to model the behaviour of High Explosives (HE). However in order to use this method, the underlying relationship between the local radius of curvature and the detonation speed must be known. Previous work has developed methods to calculate this effect using simple, single-step Arrhenius and polytropic gas, models for the chemical reaction and the equation of state, respectively. In recent years, more complex models for both reaction rates and equations of state have been developed which show better agreement with experimental data than these simple models, especially when considering condensed phase explosives.. This work presents the governing equations for solving these problems in a way that is generalised to use arbitrary equations of state as well as reaction models which may have more than a single step and multiple product species. This implementation is verified against exact solutions, demonstrating that the equations were implemented properly. The verified algorithm is then validated against experimental data and high fidelity simulations, showing that it is able to make accurate predictions in a regime where the underlying assumptions of the governing equations are valid. Importantly, this approach has many applications: from creating equivalent detonation shock dynamics models for existing reactive burn calibrations for HE; to developing new functional forms and calibrations of reactive burn models for condensed phase high explosives.

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

Critical points of the random cluster model with Newman–Ziff sampling

Here, we present a method for computing transition points of the random cluster model using a generalization of the Newman–Ziff algorithm, a celebrated technique in numerical percolation, to the random cluster model. The new method is straightforward to implement and works for real cluster weight q > 0. Furthermore, results for an arbitrary number of values of q can be found at once within a single simulation. Because the algorithm used to sweep through bond configurations is identical to that of Newman and Ziff, which was conceived for percolation, the method loses accuracy for large lattices when q > 1. However, by sampling the critical polynomial, accurate estimates of critical points in two dimensions can be found using relatively small lattice sizes, which we demonstrate here by computing critical points for non-integer values of q on the square lattice, to compare with the exact solution, and on the unsolved non-planar square matching lattice. The latter results would be much more difficult to obtain using other techniques.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Quantum Davidson algorithm for excited states

Abstract Excited state properties play a pivotal role in various chemical and physical phenomena, such as charge separation and light emission. However, the primary focus of most existing quantum algorithms has been the ground state, as seen in quantum phase estimation and the variational quantum eigensolver (VQE). Although VQE-type methods have been extended to explore excited states, these methods grapple with optimization challenges. In contrast, the quantum Krylov subspace (QKS) method has been introduced to address both ground and excited states, positioning itself as a cost-effective alternative to quantum phase estimation. However, conventional QKS methodologies depend on a pre-generated subspace through real or imaginary-time evolutions. This subspace is inherently expansive and can be plagued with issues like slow convergence or numerical instabilities, often leading to relatively deep circuits. Our research presents an economic QKS algorithm, which we term the quantum Davidson (QDavidson) algorithm. This innovation hinges on the iterative expansion of the Krylov subspace and the incorporation of a pre-conditioner within the Davidson framework. By using the residues of eigenstates to expand the Krylov subspace, we manage to formulate a compact subspace that aligns closely with the exact solutions. This iterative subspace expansion paves the way for a more rapid convergence in comparison to other QKS techniques, such as the quantum Lanczos. Using quantum simulators, we employ the novel QDavidson algorithm to delve into the excited state properties of various systems, spanning from the Heisenberg spin model to real molecules. Compared to the existing QKS methods, the QDavidson algorithm not only converges swiftly but also demands a significantly shallower circuit. This efficiency establishes the QDavidson method as a pragmatic tool for elucidating both ground and excited state properties on quantum computing platforms.

97 MATHEMATICS AND COMPUTING↗