76 research outputs found
PDEs with Compressed Solutions
Sparsity plays a central role in recent developments in signal processing,
linear algebra, statistics, optimization, and other fields. In these
developments, sparsity is promoted through the addition of an norm (or
related quantity) as a constraint or penalty in a variational principle. We
apply this approach to partial differential equations that come from a
variational quantity, either by minimization (to obtain an elliptic PDE) or by
gradient flow (to obtain a parabolic PDE). Also, we show that some PDEs can be
rewritten in an form, such as the divisible sandpile problem and
signum-Gordon. Addition of an term in the variational principle leads to
a modified PDE where a subgradient term appears. It is known that modified PDEs
of this form will often have solutions with compact support, which corresponds
to the discrete solution being sparse. We show that this is advantageous
numerically through the use of efficient algorithms for solving based
problems.Comment: 21 pages, 15 figure
The elastic field of a surface step: The Marchenko–Parshin formula in the linear case
AbstractStrain has significance for both the growth characteristics and material properties of thin epitaxial films. In this work, the method of lattice statics is applied to an epitaxial system with cubic symmetry, using harmonic potentials. The energy density and force balance equations are written using a finite difference formalism that clearly shows their consistency with continuum elasticity. For simplicity, the atomic interactions are assumed to be maximally localized. For a layered material system with a material/vacuum interface and with surface steps, force balance equations are derived, and intrinsic surface stress at the material/vacuum interface is included by treating the atoms at the surface as having different elastic properties. By defining the strain relative to an appropriately chosen nonequilibrium lattice, as in the method of eigenstrains, analytic formulas in terms of microscopic parameters are found for the local force field near a step and for the macroscopic monopole and dipole moment forces due to a step. These results provide an atomistic validation of the Marchenko–Parshin formula for the dipole moment in terms of the elastic surface stress
Sequential quasi-Monte Carlo: Introduction for Non-Experts, Dimension Reduction, Application to Partly Observed Diffusion Processes
SMC (Sequential Monte Carlo) is a class of Monte Carlo algorithms for
filtering and related sequential problems. Gerber and Chopin (2015) introduced
SQMC (Sequential quasi-Monte Carlo), a QMC version of SMC. This paper has two
objectives: (a) to introduce Sequential Monte Carlo to the QMC community, whose
members are usually less familiar with state-space models and particle
filtering; (b) to extend SQMC to the filtering of continuous-time state-space
models, where the latent process is a diffusion. A recurring point in the paper
will be the notion of dimension reduction, that is how to implement SQMC in
such a way that it provides good performance despite the high dimension of the
problem.Comment: To be published in the proceedings of MCMQMC 201
Optimal time decay of the non cut-off Boltzmann equation in the whole space
In this paper we study the large-time behavior of perturbative classical
solutions to the hard and soft potential Boltzmann equation without the angular
cut-off assumption in the whole space \threed_x with \DgE. We use the
existence theory of global in time nearby Maxwellian solutions from
\cite{gsNonCutA,gsNonCut0}. It has been a longstanding open problem to
determine the large time decay rates for the soft potential Boltzmann equation
in the whole space, with or without the angular cut-off assumption
\cite{MR677262,MR2847536}. For perturbative initial data, we prove that
solutions converge to the global Maxwellian with the optimal large-time decay
rate of O(t^{-\frac{\Ndim}{2}+\frac{\Ndim}{2r}}) in the
L^2_\vel(L^r_x)-norm for any .Comment: 31 pages, final version to appear in KR
Lattice Boltzmann simulations of soft matter systems
This article concerns numerical simulations of the dynamics of particles
immersed in a continuum solvent. As prototypical systems, we consider colloidal
dispersions of spherical particles and solutions of uncharged polymers. After a
brief explanation of the concept of hydrodynamic interactions, we give a
general overview over the various simulation methods that have been developed
to cope with the resulting computational problems. We then focus on the
approach we have developed, which couples a system of particles to a lattice
Boltzmann model representing the solvent degrees of freedom. The standard D3Q19
lattice Boltzmann model is derived and explained in depth, followed by a
detailed discussion of complementary methods for the coupling of solvent and
solute. Colloidal dispersions are best described in terms of extended particles
with appropriate boundary conditions at the surfaces, while particles with
internal degrees of freedom are easier to simulate as an arrangement of mass
points with frictional coupling to the solvent. In both cases, particular care
has been taken to simulate thermal fluctuations in a consistent way. The
usefulness of this methodology is illustrated by studies from our own research,
where the dynamics of colloidal and polymeric systems has been investigated in
both equilibrium and nonequilibrium situations.Comment: Review article, submitted to Advances in Polymer Science. 16 figures,
76 page
- …