Related papers: Computing electrostatic potentials using regulariz…
We develop a structure-preserving numerical discretization for the electrostatic Euler-Poisson equations with a constant magnetic field. The scheme preserves positivity of the density, positivity of the internal energy and a minimum…
In this paper we have derived explicitly computable bounds on the error in energy norms for the fully nonlinear Poisson-Boltzmann equation. Together with the computable bounds, we have also obtained efficient error indicators which can…
For low dimension systems admitting a moment equation representation (MER), the development of an effective eigenenergy bounding theory applicable to all discrete states had remained elusive, until now. Whereas Handy et al (1988 Phys. Rev.…
Fast and accurate resolution of electromagnetic problems via the \ac{BEM} is oftentimes challenged by conditioning issues occurring in three distinct regimes: (i) when the frequency decreases and the discretization density remains constant,…
Machine learning has been successfully applied to various fields of scientific computing in recent years. In this work, we propose a sparse radial basis function neural network method to solve elliptic partial differential equations (PDEs)…
Applications in quantitative finance such as optimal trade execution, risk management of options, and optimal asset allocation involve the solution of high dimensional and nonlinear Partial Differential Equations (PDEs). The connection…
Probabilistic encoding introduces Gaussian noise into neural networks, enabling a smooth transition from deterministic to uncertain states and enhancing generalization ability. However, the randomness of Gaussian noise distorts point-based…
We propose and analyze an extended Fourier pseudospectral (eFP) method for the spatial discretization of the Gross-Pitaevskii equation (GPE) with low regularity potential by treating the potential in an extended window for its discrete…
This paper presents a quantum algorithm for solving the fractional Poisson equation \((-\Delta)^s u = f\) with \(s \in (0,1)\) on bounded domains. The proposed approach combines rational approximation techniques with quantum linear system…
An efficient and easy-to-implement method is proposed to regularize integral equations in the 3D boundary element method (BEM). The method takes advantage of an assumed three-noded triangle discretization of the boundary surfaces. The…
A new recovery technique based on discretization corrected particle strength exchange (DC PSE) operators is developed in this paper. DC PSE is a collocation method that can be used to compute derivatives directly at nodal points, instead of…
In this article we propose and investigate a hierarchy of mathematical models based on partial differential equations (PDE) and ordinary differential equations (ODE) for the simulation of the biophysical phenomena occurring in the…
The numerical approximation of partial differential equations (PDEs) poses formidable challenges in high dimensions since classical grid-based methods suffer from the so-called curse of dimensionality. Recent attempts rely on a combination…
In this article, we propose and analyze a fully coupled, nonlinear, and energy-stable virtual element method (VEM) for solving the coupled Poisson-Nernst-Planck (PNP) and Navier--Stokes (NS) equations modeling microfluidic and…
We develop a variational regularisation framework that enables analytical solutions of the stationary de~Broglie--Bohm wave equation. The formulation begins with a Fisher-information-augmented action functional for the probability density…
We provide an integration of the universal, perturbative explicitly correlated [2]$_\text{R12}$-correction in the context of the Variational Quantum Eigensolver (VQE). This approach is able to increase the accuracy of the underlying…
Simulation of plasmas in electromagnetic fields requires numerical solution of a kinetic equation that describes the time evolution of the particle distribution function. In this paper we propose a finite volume scheme based on integral…
Dimension reduction techniques are often used when the high-dimensional tensor has relatively low intrinsic rank compared to the ambient dimension of the tensor. The CANDECOMP/PARAFAC (CP) tensor completion is a widely used approach to find…
The object of this paper is a one-dimensional generalized porous media equation (PDE) with possibly discontinuous coefficient $\beta$, which is well-posed as an evolution problem in $L^1(\mathbb{R})$. In some recent papers of Blanchard et…
We present teraflop-scale calculations of biomolecular electrostatics enabled by the combination of algorithmic and hardware acceleration. The algorithmic acceleration is achieved with the fast multipole method (FMM) in conjunction with a…