 Award ID(s):
 1907667
 NSFPAR ID:
 10326231
 Date Published:
 Journal Name:
 Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences
 Volume:
 380
 Issue:
 2226
 ISSN:
 1364503X
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

null (Ed.)Nonlinear differential equations model diverse phenomena but are notoriously difficult to solve. While there has been extensive previous work on efficient quantum algorithms for linear differential equations, the linearity of quantum mechanics has limited analogous progress for the nonlinear case. Despite this obstacle, we develop a quantum algorithm for dissipative quadratic ndimensional ordinary differential equations. Assuming R < 1 , where R is a parameter characterizing the ratio of the nonlinearity and forcing to the linear dissipation, this algorithm has complexity T 2 q poly ( log T , log n , log 1 / ϵ ) / ϵ , where T is the evolution time, ϵ is the allowed error, and q measures decay of the solution. This is an exponential improvement over the best previous quantum algorithms, whose complexity is exponential in T. While exponential decay precludes efficiency, driven equations can avoid this issue despite the presence of dissipation. Our algorithm uses the method of Carleman linearization, for which we give a convergence theorem. This method maps a system of nonlinear differential equations to an infinitedimensional system of linear differential equations, which we discretize, truncate, and solve using the forward Euler method and the quantum linear system algorithm. We also provide a lower bound on the worstcase complexity of quantum algorithms for general quadratic differential equations, showing that the problem is intractable for R ≥ 2 . Finally, we discuss potential applications, showing that the R < 1 condition can be satisfied in realistic epidemiological models and giving numerical evidence that the method may describe a model of fluid dynamics even for larger values of R.more » « less

We consider the interaction between a free flowing fluid and a porous medium flow, where the free flowing fluid is described using the time dependent Stokes equations, and the porous medium flow is described using Darcy’s law in the primal formulation. To solve this problem numerically, we use a diffuse interface approach, where the weak form of the coupled problem is written on an extended domain which contains both Stokes and Darcy regions. This is achieved using a phasefield function which equals one in the Stokes region and zero in the Darcy region, and smoothly transitions between these two values on a diffuse region of width (more » « less
ϵ ) around the StokesDarcy interface. We prove convergence of the diffuse interface formulation to the standard, sharp interface formulation, and derive rates of convergence. This is performed by deriving a priori error estimates for discretizations of the diffuse interface method, and by analyzing the modeling error of the diffuse interface approach at the continuous level. The convergence rates are also shown computationally in a numerical example. 
Abstract We study Bayesian data assimilation (filtering) for timeevolution Partial differential equations (PDEs), for which the underlying forward problem may be very unstable or illposed. Such PDEs, which include the Navier–Stokes equations of fluid dynamics, are characterized by a high sensitivity of solutions to perturbations of the initial data, a lack of rigorous global wellposedness results as well as possible nonconvergence of numerical approximations. Under very mild and readily verifiable general hypotheses on the forward solution operator of such PDEs, we prove that the posterior measure expressing the solution of the Bayesian filtering problem is stable with respect to perturbations of the noisy measurements, and we provide quantitative estimates on the convergence of approximate Bayesian filtering distributions computed from numerical approximations. For the Navier–Stokes equations, our results imply uniform stability of the filtering problem even at arbitrarily small viscosity, when the underlying forward problem may become illposed, as well as the compactness of numerical approximants in a suitable metric on timeparametrized probability measures.more » « less

This study aims to bridge length scales in immiscible multiphase flow simulation by connecting two published governing equations at the porescale and continuumscale through a novel validation framework. We employ Niessner and Hassnaizadeh's [“A model for twophase flow in porous media including fluidfluid interfacial area,” Water Resour. Res. 44(8), W08439 (2008)] continuumscale model for multiphase flow in porous media, combined with the geometric equation of state of McClure et al. [“Modeling geometric state for fluids in porous media: Evolution of the Euler characteristic,” Transp. Porous Med. 133(2), 229–250 (2020)]. Porescale fluid configurations simulated with the latticeBoltzmann method are used to validate the continuumscale results. We propose a mapping from the continuumscale to porescale utilizing a generalized additive model to predict nonwetting phase Euler characteristics during imbibition, effectively bridging the continuumtopore length scale gap. Continuumscale simulated measures of specific interfacial area, saturation, and capillary pressure are directly compared to upscaled porescale simulation results. This research develops a numerical framework capable of capturing multiscale flow equations establishing a connection between porescale and continuumscale simulations.

In numerical simulations of binary neutron star systems, the equation of state of the dense neutron star matter is an important factor in determining both the physical realism and the numerical accuracy of the simulations. Some equations of state used in simulations are 𝐶2 or smoother in the pressure/density relationship function, such as a polytropic equation of state, but may not have the flexibility to model stars or remnants of different masses while keeping their radii within known astrophysical constraints. Other equations of state, such as tabular or piecewise polytropic, may be flexible enough to model additional physics and multiple stars’ masses and radii within known constraints, but are not as smooth, resulting in additional numerical error. We will study in this paper a recently developed family of equation of state, using a spectral expansion with sufficient free parameters to allow for a larger flexibility than current polytropic equations of state, and with sufficient smoothness to reduce numerical errors compared to tabulated or piecewise polytropic equations of state. We perform simulations at three mass ratios with a common chirp mass, using two distinct spectral equations of state, and at multiple numerical resolutions. We evaluate the gravitational waves produced from these simulations, comparing the phase error between resolutions and equations of state, as well as with respect to analytical models. From our simulations, we estimate that the phase difference at the merger for binaries with a dimensionless weighted tidal deformability difference greater than ΔΛ≈55 can be captured by the spectral Einstein code for these equations of state.more » « less