Related papers: From local equilibrium to numerical PDE: Metropoli…
We study the relaxation of the Metropolis Monte Carlo algorithm corresponding to a single particle trapped in a one-dimensional confining potential, with even jump distributions that ensure that the dynamics verifies detailed balance.…
Some recent analytical papers have explored limiting behaviors of Landau-deGennes models for liquid crystals in certain extreme ranges of the model parameters: limits of "vanishing elasticity" (in the language of some of these papers) and…
We study a model, introduced initially by Gates and Westcott to describe crystal growth evolution, which belongs to the Anisotropic KPZ universality class. It can be thought of as a $(2+1)$-dimensional generalisation of the well known…
We present a nonlinear dynamical approximation method for time-dependent Partial Differential Equations (PDEs). The approach makes use of parametrized decoder functions, and provides a general, and principled way of understanding and…
We consider Markov models of large-scale networks where nodes are characterized by their local behavior and by a mobility model over a two-dimensional lattice. By assuming random walk, we prove convergence to a system of partial…
We present explicit methods for simulating diffusions whose generator is self-adjoint with respect to a known (but possibly not normalizable) density. These methods exploit this property and combine an optimized Runge-Kutta algorithm with a…
We present an initial implementation of a probabilistic PDE-constrained shape optimization algorithm. Our method is based on a novel probabilistic representation of the shape derivative, which is evaluated using Monte Carlo sampling; and…
Motivated by a model proposed by Peng et al. [Advances in Coll. and Interf. Sci. 206 (2014)] for break-up of tear films on human eyes, we study the dynamics of a generalized thin film model. The governing equations form a fourth-order…
We describe numerical simulations of the stochastic diffusion equation with a conserved charge. We focus on the dynamics in the vicinity of a critical point in the Ising universality class. The model we consider is expected to describe the…
In this short note we consider a nonlinear and spatially nonlocal PDE modelling moisture evolution in a porous medium. We then show that it naturally arises as a description of superdiffusive jump phenomenon occurring in the medium. We…
We consider the problem of solving partial differential equations (PDEs) in domains with complex microparticle geometry that is impractical, or intractable, to model explicitly. Drawing inspiration from volume rendering, we propose tackling…
In recent years several local extrema based methodologies have been proposed to investigate either the nonlinear or the nonstationary time series for scaling analysis. In the present work we study systematically the distribution of the…
As opposed to the distributed control of parabolic PDE's, very few contributions currently exist pertaining to the Dirichlet boundary condition control for parabolic PDE's. This motivates our interest in the Dirichlet boundary condition…
We propose novel less diffusive schemes for conservative one- and two-dimensional hyperbolic systems of nonlinear partial differential equations (PDEs). The main challenges in the development of accurate and robust numerical methods for the…
This work aims at making a comprehensive contribution in the general area of parametric inference for discretely observed diffusion processes. Established approaches for likelihood-based estimation invoke a time-discretisation scheme for…
In this paper we shall consider optimal scaling problems for high-dimensional Metropolis--Hastings algorithms where updates can be chosen to be lower dimensional than the target density itself. We find that the optimal scaling rule for the…
In this paper we present a classical Monte Carlo simulation of the orthorhombic phase of crystalline polyethylene, using an explicit atom force field with unconstrained bond lengths and angles and periodic boundary conditions. We used a…
Inverse problem is ubiquitous in science and engineering, and Bayesian methodologies are often used to infer the underlying parameters. For high dimensional temporal-spatial models, classical Markov chain Monte Carlo (MCMC) methods are…
This work unifies pseudo-time and inexact regularization techniques for nonmonotone classes of partial differential equations, into a regularized pseudo-time framework. Convergence of the residual at the predicted rate is investigated…
Among random sampling methods, Markov Chain Monte Carlo algorithms are foremost. Using a combination of analytical and numerical approaches, we study their convergence properties towards the steady state, within a random walk Metropolis…