An official website of the United States government
Abstract Quantum computers promise a qualitative speedup in solving a broad spectrum of practical optimization problems. The latter can be mapped onto the task of finding low-energy states of spin glasses, which is known to be exceedingly difficult. Using D-Wave’s 5000-qubit quantum processor, we demonstrate that a recently proposed iterative cyclic quantum annealing algorithm can find deep low-energy states in record time. We also find intricate structures in a low-energy landscape of spin glasses, such as a power-law distribution of connected clusters with a small surface energy. These observations offer guidance for further improvement of the optimization algorithms.
more »
« less
This content will become publicly available on March 26, 2026
Extracting the Temperature Analytically in Hydrodynamics Simulations with Gas and Radiation Pressure
Abstract Numerical hydrodynamics simulations of gases dominated by ideal, nondegenerate matter pressure and thermal radiation pressure in equilibrium entail finding the temperature as part of the evolution. Since the temperature is not typically a variable that is evolved independently, it must be extracted from the evolved variables (e.g., the rest-mass density and specific internal energy). This extraction requires solving a quartic equation, which, in many applications, is done numerically using an iterative root-finding method. Here we show instead how the equation can be solved analytically and provide explicit expressions for the solution. We also derive Taylor expansions in limiting regimes and discuss the respective advantages and disadvantages of the iterative versus analytic approaches to solving the quartic.
more »
« less
- PAR ID:
- 10591211
- Publisher / Repository:
- The Astrophysical Journal
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 982
- Issue:
- 2
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 157
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We develop two totally decoupled, linear and second‐order accurate numerical methods that are unconditionally energy stable for solving the Cahn–Hilliard–Darcy equations for two phase flows in porous media or in a Hele‐Shaw cell. The implicit‐explicit Crank–Nicolson leapfrog method is employed for the discretization of the Cahn–Hiliard equation to obtain linear schemes. Furthermore the artificial compression technique and pressure correction methods are utilized, respectively, so that the Cahn–Hiliard equation and the update of the Darcy pressure can be solved independently. We establish unconditionally long time stability of the schemes. Ample numerical experiments are performed to demonstrate the accuracy and robustness of the numerical methods, including simulations of the Rayleigh–Taylor instability, the Saffman–Taylor instability (fingering phenomenon).more » « less
-
Abstract Melt pool dynamics in metal additive manufacturing (AM) is critical to process stability, microstructure formation, and final properties of the printed materials. Physics-based simulation, including computational fluid dynamics (CFD), is the dominant approach to predict melt pool dynamics. However, the physics-based simulation approaches suffer from the inherent issue of very high computational cost. This paper provides a physics-informed machine learning method by integrating the conventional neural networks with the governing physical laws to predict the melt pool dynamics, such as temperature, velocity, and pressure, without using any training data on velocity and pressure. This approach avoids solving the nonlinear Navier–Stokes equation numerically, which significantly reduces the computational cost (if including the cost of velocity data generation). The difficult-to-determine parameters' values of the governing equations can also be inferred through data-driven discovery. In addition, the physics-informed neural network (PINN) architecture has been optimized for efficient model training. The data-efficient PINN model is attributed to the extra penalty by incorporating governing PDEs, initial conditions, and boundary conditions in the PINN model.more » « less
-
Abstract The thermal equation of state (TEOS) for solids is a mathematic model among pressure, temperature and density, and is essential for geophysical, geochemical, and other high pressure–temperature (high P–T) researches. However, in the last few decades, there has been a growing concern about the accuracy of the pressure scales of the calibrants, and efforts have been made to improve it by either introducing a reference standard or building new thermal pressure models. The existing thermal equation of state,P(V,T) = P(V,T0) + Pth(V,T), consists of an isothermal compression and an isochoric heating, while the thermal pressure is the pressure change in the isochoric heating. In this paper, we demonstrate that, for solids in a soft pressure medium in a diamond anvil cell, the thermal pressure can neither be determined from a single heating process, nor from the thermal pressure of its calibrant. To avoid the thermal pressure, we propose to replace the thermal pressure with a well-known thermal expansion model, and integrate it with the isothermal compression model to yields a Birch–Murnaghan-expansion TEOS model, called VPT TEOS. The predicted pressure of MgO and Au at ambient pressure from Birch–Murnaghan-expansion VPT TEOS model matches the experimental pressure of zero (0) GPa very well, while the pressure prediction from the approximated Anderson PVT TEOS exhibit a big deviation and a wrong trend.more » « less
-
ABSTRACT In many rotating fluids, the lowest order force balance is between gravity, pressure, and rotational acceleration (‘GPR’ balance). Terrestrial GPR balance takes the form of geostrophy and hydrostasy, which together yield the terrestrial thermal wind equation. By contrast, stellar GPR balance is an oblateness equation, which determines the departures of the thermal variables from spherical symmetry; its curl yields the ‘stellar thermal wind equation’. In this sense, the stellar thermal wind should be viewed not as a consequence of geostrophy, but of baroclinicity in the oblateness. Here, we treat the full stellar oblateness, including the thermal wind, using pressure coordinates. We derive the generalized stellar thermal wind equation and identify the parameter regime for which it holds. In the case of the Sun, not considering the full oblateness has resulted in conflicting calculations of the theoretical aspherical temperature anomaly. We provide new calculation here and find that the baroclinic anomaly is ∼3–60 times smaller than the barotropic anomaly. Thus, the anomaly from the thermal wind may not be measurable helioseismically, but if measurement were possible, this would potentially yield a new way to bracket the depth of the solar tachocline.more » « less
