We prove the existence and uniqueness of positive analytical solutions with positive initial data to the mean field equation (the Dyson equation) of the Dyson Brownian motion through the complex Burgers equation with a force term on the upper half complex plane. These solutions converge to a steady state given by Wigner’s semicircle law. A unique global weak solution with nonnegative initial data to the Dyson equation is obtained, and some explicit solutions are given by Wigner’s semicircle laws. We also construct a bi-Hamiltonian structure for the system of real and imaginary components of the complex Burgers equation (coupled Burgers system). We establish a kinetic formulation for the coupled Burgers system and prove the existence and uniqueness of entropy solutions. The coupled Burgers system in Lagrangian variable naturally leads to two interacting particle systems, the Fermi–Pasta–Ulam–Tsingou model with nearest-neighbor interactions, and the Calogero–Moser model. These two particle systems yield the same Lagrangian dynamics in the continuum limit.
more »
« less
Analytic regularity of strong solutions for the complexified stochastic nonlinear Poisson-Boltzmann equation
The nonlinear Poisson-Boltzmann equation (nPBE) is a fundamental partial differential equation (PDE) in electrostatics, widely used in computational biology and chemistry to model potential fields in solvents or plasmas. In this paper, we consider the problem of quantifying the statistical uncertainty of the stochastic nPBE solution under random variations in its coefficients. We establish the existence and uniqueness of solutions of the complexified nPBE using a contraction mapping argument, as conventional convex optimization techniques for the real-valued nPBE do not naturally extend to the complex setting. Using the existence and uniqueness result, we demonstrate that the solutions admit analytic extensions over a well-defined region in the complex hyperplane The analyticity makes the computation for statistics of real-valued quantities of interest amenable to numerical techniques such as sparse grids. Sparse grids are applied to uniformly approximate analytic functions with algebraic to sub-exponential error with respect to the number of knots, thus allowing for efficient approximations of high-dimensional integrals. Our numerical experiments confirm the predicted error behavior.
more »
« less
- Award ID(s):
- 1736392
- PAR ID:
- 10667247
- Publisher / Repository:
- Elsevier
- Date Published:
- Journal Name:
- Computers & Mathematics with Applications
- Volume:
- 204
- Issue:
- C
- ISSN:
- 0898-1221
- Page Range / eLocation ID:
- 97 to 122
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract From a unified vision of vector valued solutions in weighted Banach spaces, this paper establishes the existence and uniqueness for space homogeneous Boltzmann bi-linear systems with conservative collisional forms arising in complex gas dynamical structures. This broader vision is directly applied to dilute multi-component gas mixtures composed of both monatomic and polyatomic gases. Such models can be viewed as extensions of scalar Boltzmann binary elastic flows, as much as monatomic gas mixtures with disparate masses and single polyatomic gases, providing a unified approach for vector valued solutions in weighted Banach spaces. Novel aspects of this work include developing the extension of a general ODE theory in vector valued weighted Banach spaces, precise lower bounds for the collision frequency in terms of the weighted Banach norm, energy identities, angular or compact manifold averaging lemmas which provide coerciveness resulting into global in time stability, a new combinatorics estimate forp-binomial forms producing sharper estimates for thek-moments of bi-linear collisional forms. These techniques enable the Cauchy problem improvement that resolves the model with initial data corresponding to strictly positive and bounded initial vector valued mass and total energy, in addition to only a$$2^+$$ moment determined by the hard potential rates discrepancy, a result comparable in generality to the classical Cauchy theory of the scalar homogeneous Boltzmann equation.more » « less
-
We establish both the uniqueness and the existence of the solutions to a hidden-memory variable-order fractional stochastic partial differential equation, which models, e.g., the stochastic motion of a Brownian particle within a viscous liquid medium varied with fractal dimensions. We also investigate the inverse problem concerning the observations of the solutions, which eliminates the analytic assumptions on the variable orders in the literature of this topic and theoretically guarantees the reliability of the determination and experimental inference.more » « less
-
We study a regularized version of the Landau equation, which was recently introduced in [J. A. Carrillo, J. Hu, L. Wang and J. Wu, A particle method for the homogeneous Landau equation, J. Comput. Phys. X 7 (2020) 100066, 24] to numerically approximate the Landau equation with good accuracy at reasonable computational cost. We develop the existence and uniqueness theory for weak solutions, and we reinforce the numerical findings in the above-mentioned paper by rigorously proving the validity of particle approximations to the regularized Landau equation.more » « less
-
null (Ed.)Seismic waves in earth media usually undergo attenuation, causing energy losses and phase distortions. In the regime of high-frequency asymptotics, a complex-valued eikonal is an essential ingredient for describing wave propagation in attenuating media, where the real and imaginary parts of the eikonal function capture dispersion effects and amplitude attenuation of seismic waves, respectively. Conventionally, such a complex-valued eikonal is mainly computed either by tracing rays exactly in complex space or by tracing rays approximately in real space so that the resulting eikonal is distributed irregularly in real space. However, seismic data processing methods, such as prestack depth migration and tomography, usually require uniformly distributed complex-valued eikonals. Therefore, we have developed a unified framework to Eulerianize several popular approximate real-space ray-tracing methods for complex-valued eikonals so that the real and imaginary parts of the eikonal function satisfy the classic real-space eikonal equation and a novel real-space advection equation, respectively, and we dub the resulting method the Eulerian partial-differential-equation method. We further develop highly efficient high-order methods to solve these two equations by using the factorization idea and the Lax-Friedrichs weighted essentially nonoscillatory schemes. Numerical examples demonstrate that our method yields highly accurate complex-valued eikonals, analogous to those from ray-tracing methods. Our methods can be useful for migration and tomography in attenuating media.more » « less
An official website of the United States government

