Diffusion in arrays of obstacles: beyond homogenisation

by   Yahya Farah, et al.

We revisit the classical problem of diffusion of a scalar (or heat) released in a two-dimensional medium with an embedded periodic array of impermeable obstacles such as perforations. Homogenisation theory provides a coarse-grained description of the scalar at large times and predicts that it diffuses with a certain effective diffusivity, so the concentration is approximately Gaussian. We improve on this by developing a large-deviation approximation which also captures the non-Gaussian tails of the concentration through a rate function obtained by solving a family of eigenvalue problems. We focus on cylindrical obstacles and on the dense limit, when the obstacles occupy a large area fraction and non-Gaussianity is most marked. We derive an asymptotic approximation for the rate function in this limit, valid uniformly over a wide range of distances. We use finite-element implementations to solve the eigenvalue problems yielding the rate function for arbitrary obstacle area fractions and an elliptic boundary-value problem arising in the asymptotics calculation. Comparison between numerical results and asymptotic predictions confirm the validity of the latter.



There are no comments yet.


page 6

page 14


Numerical solution of boundary value problems for the eikonal equation in an anisotropic medium

A Dirichlet problem is considered for the eikonal equation in an anisotr...

On the approximation of dispersive electromagnetic eigenvalue problems in 2D

We consider time-harmonic electromagnetic wave equations in composites o...

On the finite element approximation of fourth order singularly perturbed eigenvalue problems

We consider fourth order singularly perturbed eigenvalue problems in one...

Frequency extraction for BEM-matrices arising from the 3D scalar Helmholtz equation

The discretisation of boundary integral equations for the scalar Helmhol...

NEP: a module for the parallel solution of nonlinear eigenvalue problems in SLEPc

SLEPc is a parallel library for the solution of various types of large-s...

On the rate of concentration of maxima in Gaussian arrays

Recently in Gao and Stoev (2018) it was established that the concentrati...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.