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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Exactly Solvable Model for Strongly Interacting Electrons in a Magnetic Field

States of strongly interacting particles are of fundamental interest in physics and can produce exotic emergent phenomena and topological structures. We consider here two-dimensional electrons in a magnetic field, and, departing from the standard practice of restricting to the lowest LL, introduce a model short-range interaction that is infinitely strong compared to the cyclotron energy. Here, we demonstrate that this model lends itself to an exact solution for the ground as well as excited states at arbitrary filling factors v<1/2p and produces a fractional quantum Hall effect at fractions of the form v=n/(2pn+1), where n and p are integers. The fractional quantum Hall states of our model share many topological properties with the corresponding Coulomb ground states in the lowest Landau level, such as the edge physics and the fractional charge of the excitations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1997

Least-Squares Spectral Element Solutions to the CAA Workshop Benchmark Problems

This paper presents computed results for some of the CAA benchmark problems via the acoustic solver developed at Rocketdyne CFD Technology Center under the corporate agreement between Boeing North American, Inc. and NASA for the Aerospace Industry Technology Program. The calculations are considered as benchmark testing of the functionality, accuracy, and performance of the solver. Results of these computations demonstrate that the solver is capable of solving the propagation of aeroacoustic signals. Testing of sound generation and on more realistic problems is now pursued for the industrial applications of this solver. Numerical calculations were performed for the second problem of Category 1 of the current workshop problems for an acoustic pulse scattered from a rigid circular cylinder, and for two of the first CAA workshop problems, i. e., the first problem of Category 1 for the propagation of a linear wave and the first problem of Category 4 for an acoustic pulse reflected from a rigid wall in a uniform flow of Mach 0.5. The aim for including the last two problems in this workshop is to test the effectiveness of some boundary conditions set up in the solver. Numerical results of the last two benchmark problems have been compared with their corresponding exact solutions and the comparisons are excellent. This demonstrates the high fidelity of the solver in handling wave propagation problems. This feature lends the method quite attractive in developing a computational acoustic solver for calculating the aero/hydrodynamic noise in a violent flow environment.

Lin, Wen H.↗

NASA NTRS2004

Three-Dimensional High-Order Spectral Volume Method for Solving Maxwell's Equations on Unstructured Grids

A three-dimensional, high-order, conservative, and efficient discontinuous spectral volume (SV) method for the solutions of Maxwell's equations on unstructured grids is presented. The concept of discontinuous 2nd high-order loca1 representations to achieve conservation and high accuracy is utilized in a manner similar to the Discontinuous Galerkin (DG) method, but instead of using a Galerkin finite-element formulation, the SV method is based on a finite-volume approach to attain a simpler formulation. Conventional unstructured finite-volume methods require data reconstruction based on the least-squares formulation using neighboring cell data. Since each unknown employs a different stencil, one must repeat the least-squares inversion for every cell at each time step, or to store the inversion coefficients. In a high-order, three-dimensional computation, the former would involve impractically large CPU time, while for the latter the memory requirement becomes prohibitive. In the SV method, one starts with a relatively coarse grid of triangles or tetrahedra, called spectral volumes (SVs), and partition each SV into a number of structured subcells, called control volumes (CVs), that support a polynomial expansion of a desired degree of precision. The unknowns are cell averages over CVs. If all the SVs are partitioned in a geometrically similar manner, the reconstruction becomes universal as a weighted sum of unknowns, and only a few universal coefficients need to be stored for the surface integrals over CV faces. Since the solution is discontinuous across the SV boundaries, a Riemann solver is thus necessary to maintain conservation. In the paper, multi-parameter and symmetric SV partitions, up to quartic for triangle and cubic for tetrahedron, are first presented. The corresponding weight coefficients for CV face integrals in terms of CV cell averages for each partition are analytically determined. These discretization formulas are then applied to the integral form of the Maxwell equations. All numerical procedures for outer boundary, material interface, zonal interface, and interior SV face are unified with a single characteristic formulation. The load balancing in a massive parallel computing environment is therefore easier to achieve. A parameter is introduced in the Riemann solver to control the strength of the smoothing term. Important aspects of the data structure and its effects to communication and the optimum use of cache memory are discussed. Results will be presented for plane TE and TM waves incident on a perfectly conducting cylinder for up to fifth order of accuracy, and a plane wave incident on a perfectly conducting sphere for up to fourth order of accuracy. Comparisons are made with exact solutions for these cases.

Liu, Yen↗

DOE OSTI2020

Noninteracting electrons in a prototypical one-dimensional sinusoidal potential

A prototypical model of a one-dimensional metallic monatomic solid containing noninteracting electrons is studied, where the argument of the cosine potential energy, periodic with the lattice, contains the first reciprocal lattice vector G 1 =2π/a, where a is the lattice constant. The time-independent Schrödinger equation can be written in reduced variables as a Mathieu equation for which numerically exact solutions for the band structure and wave functions are obtained. The band structure has band gaps that increase with increasing amplitude q of the cosine potential. In the extended-zone scheme, the energy gaps decrease with increasing index n of the Brillouin-zone boundary ka=nπ, where k is the crystal momentum of the electron. The wave functions of the band electron are derived for various combinations of k and q as complex combinations of the real Mathieu functions with even and odd parity, and the normalization factor is discussed. The wave functions at the bottoms and tops of the bands are found to be real or imaginary, respectively, corresponding to standing waves at these energies. Irrespective of the wave vector k within the first Brillouin zone, the electron probability density is found to be periodic with the lattice. The Fourier components of the wave functions are derived versus q, which reveal multiple reciprocal-lattice-vector components with variable amplitudes in the wave functions unless q = 0. The magnitudes of the Fourier components are found to decrease exponentially as a power of n for n~3 to 45 for ka = π/2 and q = 2, and a precise fit is obtained to the data. The probability densities and probability currents obtained from the wave functions are also discussed. The probability currents are found to be zero for crystal momenta at the tops and bottoms of the energy bands, because the wave functions for these crystal momenta are standing waves. To conlcude, the band structure is calculated from the central equation and compared to the numerically exact band structure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Arbitrary Order Virtual Element Methods for High‐Order Phase‐Field Modeling of Dynamic Fracture

ABSTRACT Accurate modeling of fracture nucleation and propagation in brittle and ductile materials subjected to dynamic loading is important in predicting material damage and failure under extreme conditions. Phase‐field fracture models have garnered a lot of attention in recent years due to their success in representing damage and fracture processes in a wide class of materials and under a variety of loading conditions. Second‐order phase‐field fracture models are by far the most popular among researchers (and increasingly, among practitioners), but fourth‐order models have started to gain broader acceptance since their more recent introduction. The exact solution corresponding to these high‐order phase‐field fracture models has higher regularity. Thus, numerical solutions of the model equations can achieve improved accuracy and higher spatial convergence rates. In this work, we develop a virtual element framework for the high‐order phase‐field model of dynamic fracture. The virtual element method (VEM) can be regarded as a generalization of the classical finite element method. In addition to many other desirable characteristics, the VEM allows computing on polytopal meshes. Here, we use ‐conforming virtual elements and the generalized‐ time integration method for the momentum balance equation, and adopt ‐conforming virtual elements for the high‐order phase‐field equation. We verify our virtual element framework using classical quasi‐static benchmark problems and demonstrate its capabilities with the aid of numerical simulations of dynamic fracture in brittle materials.

42 ENGINEERING↗

DOE OSTI2025

Efficient Simulation of Open Quantum Systems on NISQ Trapped‐Ion Hardware

Abstract Simulating open quantum systems, which interact with external environments, presents significant challenges on noisy intermediate‐scale quantum (NISQ) devices due to limited qubit resources and noise. In this study, an efficient framework is proposed for simulating open quantum systems on NISQ hardware by leveraging a time‐perturbative Kraus operator representation of the system's dynamics. This approach avoids the computationally expensive Trotterization method and exploits the Lindblad master equation to represent time evolution in a compact form, particularly for systems satisfying specific commutation relations. The efficiency of this method is demonstrated by simulating quantum channels, such as the continuous‐time Pauli channel and damped harmonic oscillators, on NISQ trapped‐ion hardware, including IonQ Harmony and Quantinuum H1‐1. Additionally, hardware‐agnostic error mitigation techniques are introduced, including Pauli channel fitting and quantum depolarizing channel inversion, to enhance the fidelity of quantum simulations. These results show strong agreement between the simulations on real quantum hardware and exact solutions, highlighting the potential of Kraus‐based methods for scalable and accurate simulation of open quantum systems on NISQ devices. This framework opens pathways for simulating more complex systems under realistic conditions in the near term.

Burdine, Colin [Department of Electrical and Compu↗

DOE OSTI2024

Orbital angular momentum small- x evolution: exact results in the large- N c limit

We construct an exact solution to the revised small-x orbital angular momentum (OAM) evolution equations derived in [1], based on an earlier work [2]. These equations are derived in the double logarithmic approximation (summing powers of α s ln 2 (1/x) with α s the strong coupling constant and x the Bjorken x variable) and the large-N c limit, with N c the number of quark colors. From our solution, we extract the small-x, large-N c expressions of the quark and gluon OAM distributions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Two-loop radiative jet function for exclusive B-meson and Higgs decays

The rare radiative B -meson decay \( {B}^{-}\to \gamma {\mathrm{\ell}}^{-}\overline{v} \) and the radiative Higgs-boson decay h → γγ mediated by light-quark loops both receive large logarithmic corrections in QCD, which can be resummed using factorization theorems derived in soft-collinear effective theory. In these factorization theorems the same radiative jet function appears, which is a central object in the study of factorization beyond the leading order in scale ratios. We calculate this function at two-loop order both in momentum space and in a dual space, where its renormalization-group evolution takes on a simpler form. We also derive the two-loop anomalous dimension of the jet function and present the exact solution to its evolution equation at two-loop order. Another important outcome of our analysis is the explicit form of the two-loop anomalous dimension of the B -meson light-cone distribution amplitude in momentum space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Solitons in lattice field theories via tight-binding supersymmetry

Reflectionless potentials play an important role in constructing exact solutions to classical dynamical systems (such as the Korteweg-de Vries equation), non-perturbative solutions of various large-N field theories (such as the Gross-Neveu model), and closely related solitonic solutions to the Bogoliubov-de Gennes equations in the theory of superconductivity. These solutions rely on the inverse scattering method, which reduces these seemingly unrelated problems to identifying reflectionless potentials of an auxiliary one-dimensional quantum scattering problem. There are several ways of constructing these potentials, one of which is quantum mechanical supersymmetry (SUSY). In this paper, motivated by recent experimental platforms, we generalize this framework to develop a theory of lattice solitons. We first briefly review the classical inverse scattering method in the continuum limit, focusing on the Korteweg-de Vries (KdV) equation and SU(N) Gross-Neveu model in the large N limit. We then generalize this methodology to lattice versions of interacting field theories. Our analysis hinges on the use of trace identities, which are relations connecting the potential of an equation of motion to the scattering data. For a discrete Schrödinger operator, such trace identities had been known as far back as Toda; however, we derive a new set of identities for the discrete Dirac operator. We then use these identities in a lattice Gross-Neveu and chiral Gross-Neveu (Nambu-Jona-Lasinio) model to show that lattice solitons correspond to reflectionless potentials associated with the discrete scattering problem. These models are of significance as they are equivalent to a mean-field theory of a lattice superconductor. To explicitly construct these solitons, we generalize supersymmetric quantum mechanics to tight-binding models. We show that a matrix transformation exists that maps a tight-binding model to an isospectral one which shares the same structure and scattering properties. The corresponding soliton solutions have both modulated hopping and onsite potential, the former of which has no analogue in the continuum limit. We explicitly compute both topological and non-topological soliton solutions as well as bound state spectra in the aforementioned models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Collision of localized shocks in AdS 5 as a series expansion in transverse gradients

We introduce a computational framework to more efficiently calculate the collision of localized shocks in five dimensional asymptotically Anti-de Sitter space. We expand the Einstein equations in transverse gradients and find that our numerical results agree well with exact solutions already at first order in the expansion. Moreover, the Einstein equations at first order in transverse gradients can be decoupled into two sets of differential equations. Here, the bulk fields of one of these sets has only a negligible contribution to boundary observables, such that the computation on each time slice can be simplified to the solution of several planar shockwave equations plus four further differential equations for each transverse plane ‘pixel’. At the cost of errors of ≲ 10% at the hydrodynamization time and for low to mid rapidities, useful numerical solutions can be sped up by roughly one order of magnitude.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Double Kerr-Schild spacetimes and the Newman-Penrose map

The Newman-Penrose map, which is closely related to the classical double copy, associates certain exact solutions of Einstein’s equations with self-dual solutions of the vacuum Maxwell equations. Here we initiate an extension of the Newman-Penrose map to a broader class of spacetimes. As an example, we apply principles from the Newman-Penrose map to associate a self-dual gauge field to the Kerr-Taub-NUT-(A)dS spacetime and we show that the result agrees with previously studied examples of classical double copies. The corresponding field strength exhibits a discrete electric-magnetic duality that is distinct from its (Hodge star) self-dual property.

scattering amplitudes↗

DOE OSTI2024

2D Dilaton Gravity and the Weil–Petersson Volumes with Conical Defects

We derive the Weil-Petersson measure on the moduli space of hyperbolic surfaces with defects of arbitrary opening angles and use this to compute its volume. We conjecture a matrix integral computing the corresponding volumes and confirm agreement in simple cases. Here, we combine this mathematical result with the equivariant localization approach to Jackiw-Teitelboim gravity to justify a proposed exact solution of pure 2d dilaton gravity for a large class of dilaton potentials.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

hr-Adaptivity for nonconforming high-order meshes with the target matrix optimization paradigm

In this study, we present an \(hr\) -adaptivity framework for optimization of high-order meshes. This work extends the r -adaptivity method by Dobrev et al. (Comput Fluids, 2020), where we utilized the Target-Matrix Optimization Paradigm (TMOP) to minimize a functional that depends on each element’s current and target geometric parameters: element aspect-ratio , size , skew , and rotation . Since fixed mesh topology limits the ability to achieve the target size and aspect-ratio at each position, in this paper, we augment the r -adaptivity framework with nonconforming adaptive mesh refinement to further reduce the error with respect to the target geometric parameters. The proposed formulation, referred to as \(hr\) -adaptivity, introduces TMOP-based quality estimators to satisfy the aspect-ratio target via anisotropic refinements and size target via isotropic refinements in each element of the mesh. The methodology presented is purely algebraic, extends to both simplices and hexahedra/quadrilaterals of any order, and supports nonconforming isotropic and anisotropic refinements in 2D and 3D. Using a problem with a known exact solution, we demonstrate the effectiveness of \(hr\) -adaptivity over both r - and \(h\) -adaptivity in obtaining similar accuracy in the solution with significantly fewer mesh nodes. We also present several examples that show that \(hr\) -adaptivity can help satisfy geometric targets even when \(r\) -adaptivity fails to do so, due to the topology of the initial mesh.

42 ENGINEERING↗

DOE OSTI2026

Dynamics of McMillan mappings III. Symmetric map with mixed nonlinearity

This article extends the study of the dynamical properties of the symmetric McMillan map, emphasizing its utility in understanding and modeling complex nonlinear systems. Although the map features six parameters, we demonstrate that only two are irreducible: the linearized rotation number at the fixed point and a nonlinear parameter representing the ratio of terms in the biquadratic invariant. Through a detailed analysis, we classify regimes of stable motion, provide exact solutions to the mapping equations, and derive a canonical set of action-angle variables, offering analytical expressions for the rotation number and nonlinear tune shift. We further establish connections between general standard-form mappings and the symmetric McMillan map, using the area-preserving Hénon map and accelerator lattices with thin sextupole magnet as representative case studies. Our results show that, despite being a second-order approximation, the symmetric McMillan map provides a highly accurate depiction of dynamics across a wide range of system parameters, demonstrating its practical relevance in both theoretical and applied contexts.

43 PARTICLE ACCELERATORS↗

DOE OSTI2021

Topological order in Mott insulators

Here in this commemorative issue for Phil Anderson, I reproduce parts of a joint unpublished paper that first got me thinking about topological order in Mott insulators. With Phil, we explored aspects of the momentum distribution function and the non-classical rotational inertia in the Bose Hubbard model close to the quantum phase transition. My interest on this question continued, and later in the article I take up the question of topological order in the Kitaev Mott insulator on the honeycomb lattice with bond-dependent interactions in a finite magnetic field. Kitaev obtained the exact solution for zero field that demonstrated an exotic Z 2 quantum spin liquid phase with topological order. In a finite magnetic field, we discovered a new intermediate U(1) gapless quantum spin liquid phase, sandwiched between the known gapped Kitaev spin liquid and a polarized phase, that harbors fractionalized neutral fermionic excitations, dubbed spinons, that remarkably form a Fermi surface in a Mott insulator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

A priori error analysis of high-order LL* (FOSLL*) finite element methods

A number of non-standard finite element methods have been proposed in recent years, each of which derives from a specific class of PDE-constrained norm minimization problems. The most notable examples are LL* methods. In this work, we argue that all high-order methods in this class should be expected to deliver substandard uniform h-refinement convergence rates. In fact, one may not even see rates proportional to the polynomial order p > 1 when the exact solution is a constant function. Here, we show that the convergence rate is limited by the regularity of an extraneous Lagrange multiplier variable which naturally appears via a saddle-point analysis. In turn, limited convergence rates appear because the regularity of this Lagrange multiplier is determined, in part, by the geometry of the domain. Numerical experiments support our conclusions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Optimization of battery swapping infrastructure for e-commerce drone delivery

Drone delivery is widely-researched to alter the current e-commerce delivery convention for providing short delivery lead times. Yet, flight range of drones, constrained by the available battery technology, sets a milestone toward realizing the sole-drone delivery. To tackle with the flight range limitation, locating automated battery swapping machines (ABSM) have been proposed and a few studies modeled the problem. Using the ABSMs, drones can be loaded with fully-charged batteries along the route to a demand location. Here, we introduce a mixed-integer nonlinear program to model the problem. The objective of the program is to optimally select ABSM locations, determine the delivery-mode choices (drone-only, truck-only, and mixed delivery) of demand locations, find drone delivery routes, and approximate the baseline requirements for the number of drones and batteries needed. The program minimizes the overall delivery system costs including: ABSM, delivery, drone ownership, battery inventory, and service congestion. A cutting-plane method is developed to find exact solutions in finite iterations. Computational experiments showed that the method quickly yields the optimal solution to instances with less than 60 ABSM candidates and 20 demand locations. A case study shows that the optimal drone delivery infrastructure can save almost 20% cost compared to the conventional truck-only delivery. Sensitivity analyses were conducted to reveal the impact of key parameters in the decision-making and found that a decrease in the ABSM and drone costs highly affect the system cost.

25 ENERGY STORAGE↗

DOE OSTI2020

Designing a drone delivery network with automated battery swapping machines

Drones are projected to alter last-mile delivery, but their short travel range is a concern. In this study, we propose a drone delivery network design using automated battery swapping machines (ABSMs) to extend ranges. The design minimizes the long-term delivery costs, including ABSM investment, drone ownership, and cost of the delivery time, and locates ABSMs to serve a set of customers. We build a mixed-integer nonlinear program that captures the nonlinear waiting time of drones at ABSMs. To solve the problem, we create an exact solution algorithm that finds the globally optimal solution using a derivative-supported cutting-plane method. To validate the applicability of our program, we conduct a case study on the Chicago Metropolitan area using cost data from leading ABSM manufacturer and geographical data from the planning and operations language for agent-based regional integrated simulation (more commonly known as POLARIS). A sensitivity analysis identifies that ABSM service times and costs are the key parameters impacting the long-term adoption of drone delivery.

25 ENERGY STORAGE↗