You are viewing a javascript disabled version of the site. Please enable Javascript for this site to function properly.
Go to headerGo to navigationGo to searchGo to contentsGo to footer
In content section. Select this link to jump to navigation

Simultaneous diffusion and homogenization asymptotic for the linear Boltzmann equation

Abstract

This article is on the simultaneous diffusion approximation and homogenization of the linear Boltzmann equation when both the mean free path ε and the heterogeneity length scale η vanish. No periodicity assumption is made on the scattering coefficient of the background material. There is an assumption made on the heterogeneity length scale η that it scales as εβ for β(0,). In one space dimension, we prove that the solutions to the kinetic model converge to the solutions of an effective diffusion equation for any β2 in the ε0 limit. In any arbitrary phase space dimension, under a smallness assumption of a certain quotient involving the scattering coefficient in the H12 norm, we again prove that the solutions to the kinetic model converge to the solutions of an effective diffusion equation in the ε0 limit.

1.Introduction

The unknown quantity modeled by equations in kinetic theory is the probability distribution function f(t,x,v) of a population of particles which is a function of time, position and velocity. The linear Boltzmann equation describes the evolution of the distribution function modeling the collision of a population of particles with a background medium:

tf(t,x,v)+v·xf(t,x,v)=σ(x)(Vf(t,x,w)dμ(w)f(t,x,v)),
where μ is a Borel probability measure on the space of velocities V. The coefficient σ(x) is the scattering coefficient of the background material. Our objective is to perform an asymptotic analysis when the background medium is inhomogeneous – say, a composite material with microstructure. We consider two asymptotically small parameters associated with the above kinetic model: (a) mean free path ε of the particles between two interactions; (b) the scale of heterogeneity η of the background medium – it can be the average distance between two neighboring inhomogeneities or the average size of the inhomogeneities.

In the above kinetic model, the scattering coefficient σ(x) represents the background medium. The inhomogeneous nature of the background medium implies that smaller the parameter η is, more rapid the oscillations are in σ(x). As is standard in the theory of homogenization, we consider a family of scattering coefficients indexed by η, i.e. ση(x) and study an associated family of solutions to the kinetic model. We also wish to study the evolution of the local equilibria for the above kinetic model. This corresponds to scaling the above kinetic model using parabolic scaling with the parameter ε. The objective of this article is to study the following scaled linear Boltzmann equation

εtfε,η(t,x,v)+v·xfε,η(t,x,v)=ση(x)ε(Vfε,η(t,x,w)dμ(w)fε,η(t,x,v))
in the simultaneous limit as both the scaling parameters ε, η vanish. The diffusion approximation of the linear Boltzmann equation corresponds to the ε0 limit in the above scaled equation (see [6] and references therein for the state-of-the-art on the techniques used in the diffusion approximation of linear transport equations). Therefore in the regime εη, one can first perform the diffusion asymptotic for a fixed η yielding a parabolic equation with heterogeneous coefficient ση(x). The η0 limit corresponds to deriving an homogenized equation for the thus obtained heterogeneous parabolic equation. This can either be performed using the asymptotic expansions method [11] or the two-scale convergence method [1,21] when ση(x)=σ(x/η) is η-periodic. The method of H-convergence [19] (see also [20]) can be used when the family ση(x) is a general family of L coefficients (see [2] for a pedagogical exposition of the method of H-convergence). In the regime ηε, we need to homogenize the linear Boltzmann equation in the limit η0 for a fixed ε (see [12] on the homogenization of linear transport equations). This shall be followed by a diffusion asymptotic in the ε0 limit. These points have already been observed in an expository article by F. Golse [13].

The present article addresses the issue of simultaneous limit procedure when both the small parameters ε and η vanish. This problem has been addressed in [5,9,10,22] when the heterogeneous scattering coefficient is periodic. A first work in this direction goes back to the work of R. Sentis [22] where the heterogeneity length scale η is related to the mean free path as η=εβ with β<1. Recently, there has been a revival of this problem. The works [5,9] address this problem in the periodic setting and in the regime εη. The approach in [5,9] is to introduce a new parameter εη and study some cell problems involving the new small parameter εη. They extensively use the method of two-scale convergence. We also cite [3] where the spectral problem associated with the scaled linear Boltzmann equation is studied when η=ε. For the simultaneous limit procedures in the case η=ε, we cite [16,17] (see also Chapter 7 in the lecture notes [4]).

Our work, essentially, goes beyond the periodic setting in the spirit of the compensated compactness theory [18,19] developed by F. Murat and L. Tartar in the context of homogenization of elliptic and parabolic problems in the late 1970s. All the results and computations in this article are presented for a special case of the stationary linear Boltzmann, i.e. the probability density function is supposed to be time independent. All the results can be straightaway generalized to the time-dependent setting (see Remark 4).

Similar to the work of R. Sentis [22], we assume that the two scaling parameters are related as η=εβ where β(0,). Note, however, that the results in [22] hold only when β<1 and under periodicity assumption on the heterogeneous scattering coefficient. Our main result in one dimensional setting (Theorem 1) is that the solutions to the linear Boltzmann equation converge to the solutions of an elliptic problem whenever β2. We also obtain an explicit form of the effective diffusion coefficient in the limit equation – given in terms of some velocity averages and the weak-* limit of the heterogeneous scattering coefficient. The one-dimensional setting is very special as the divergence operator coincides with the gradient operator resulting in uniform H1 estimates for certain family of second moments – refer to Section 5 for further details.

Our main result in any arbitrary dimension is that if the heterogeneous coefficient satisfies

σ¯(x)σε(x)σ¯(x)H12(Ω)=O(ε1+)
for some σ¯L(Ω) and where the family σε(x) is nothing but the family ση(x) with η=εβ, then the solutions to the linear Boltzmann equation converge to the solutions of a diffusion equation in physical space. More precisely, we show that the family
fε,η(x,v)ρ(x)weakly in L2(Ω×V;dxdμ) as ε,η0,
where ρ(x) is the L2-weak limit of the family of local densities and that it is the unique solution to a diffusion equation. Our result in arbitrary dimension (Theorem 2) essentially employs the moments approach inspired by [7] and uses the regularity of velocity averages guaranteed by the now well-known velocity averaging lemma [14,15].

Plan of the paper. In Section 2, we present the linear kinetic model, the scaling parameters and the scaling considered in this article. Section 3 gives some uniform (with respect to the scaling parameters ε and η) estimates on the solutions to the linear transport equation and associated velocity averages. In Section 4, we briefly explain the moments method as given in [7] and apply this method to the linear Boltzmann equation. Section 5 is devoted to deriving the limit equation (in the ε0 limit) in one dimensional setting – Theorem 1. In Section 6, we give a result in any arbitrary dimension under a certain assumption on the scattering coefficient. Finally, in Section 7 we give some concluding remarks and perspectives.

2.Stationary linear Boltzmann equation

Let f(x,v) be the distribution function which depends on xΩRd (space position) and vV (velocity). The distribution function models the probability density of mono-kinetic particles interacting with the background medium. The velocity space can be either of the following:

V=Rd;V=Sd1;V=B(0,l):={vRd s.t. |v|l}.
We denote by μ a Borel probability measure on V. We further suppose that
(1)Vvdμ(v)=0.
In order to define the boundary conditions, taking n(x) to be the unit exterior normal to Ω at the point xΩ, we introduce the following notations:
Σ:={(x,v)Ω×V}Phase space Boundary,Σ+:={(x,v)Ω×V s.t. v·n(x)>0}Outgoing Boundary,Σ:={(x,v)Ω×V s.t. v·n(x)<0}Incoming Boundary,Σ0:={(x,v)Ω×V s.t. v·n(x)=0}Grazing set.
Denote by γ+f (respectively γf), the trace of f on Σ+ (respectively on Σ).

The goal is to perform the simultaneous limiting procedure for the scaled problem:

(2a)fε,η(x,v)+1εv·xfε,η(x,v)+1ε2ση(x)Lfε,η(x,v)=g(x)for (x,v)Ω×V,(2b)γfε,η(x,v)=0for (x,v)Σ,
where the linear Boltzmann operator is the following integral operator:
(3)Lg(v):=g(v)Vg(w)dμ(w)for any gL1(V;dμ).
The family of heterogeneous scattering coefficients indexed by η, i.e. ση(x) is assumed to be a family of differentiable function such that there exist uniform (with respect to η) constants a,b,c such that
0<aση(x)b<+xΩ
and
σηL(Ω)cη1.
The source term in (2a) is assumed to be square-integrable in the space variable, i.e.
gL2(Ω)C.
We suppose that the heterogeneity length scale η and the mean free path ε are related as
η=εβfor β>0.
Hence we drop the superscript η in (2a) and index the family of heterogeneous scattering coefficients by ε, i.e. σε(x) which inherits the following bounds from ση given above.
(4)0<aσε(x)b<+xΩandσεL(Ω)cεβ,
where the constants a,b,c are uniform with respect to ε.

Remark 1.

Consider σW1,(Rd) such that for a.e. xRd,

0<aσ(x)b<andσL(Rd)c<.
Now consider the family
σε(x):=σ(xεβ).
Then, the assumptions given in (4) on the heterogeneous coefficient are satisfied by the above defined family of coefficients.

The main objective of this article is to study the following stationary problem in the ε0 limit:

(5a)fε(x,v)+1εv·xfε(x,v)+1ε2σε(x)Lfε(x,v)=g(x)for (x,v)Ω×V,(5b)γfε(x,v)=0for (x,v)Σ.
We shall prove that the entire family fε(x,v) of solutions to the above kinetic equation exhibit the following compactness property:
fε(x,v)ρ(x)weakly in L2(Ω×V;dxdμ(v)),
where ρ(x) is the L2-weak limit of the associated local densities, i.e.
Vfε(x,v)dμ(v)ρ(x)weakly in L2(Ω)
and the above weak limit uniquely solves a second order elliptic equation:
Aρ(x)=g(x)in Ω.
Further details on the elliptic operator A shall be given in Sections 5 and 6.

3.Uniform a priori bounds

In order to perform the asymptotic analysis in the ε0 limit for (5a)–(5b), we derive uniform (with respect to ε) L2-estimates on the solution family fε(x,v) and some associated velocity averages.

We first show that the Dirichlet form associated with the integral operator L in L2(V;dμ) is positive semi-definite.

Lemma 1.

For any ϕL2(V;dμ), we have

Vϕ(v)Lϕ(v)dμ(v)=12V×V(ϕ(v)ϕ(w))2dμ(w)dμ(v)0.

Proof.

By the definition of the Boltzmann operator (3), we have:

Vϕ(v)Lϕ(v)dμ(v)=V|ϕ(v)|2dμ(v)V×Vϕ(v)ϕ(w)dμ(w)dμ(v).
Split the first integral on the right hand side of the above expression into two, thus yielding
12V|ϕ(v)|2dμ(v)+12V|ϕ(w)|2dμ(w)V×Vϕ(v)ϕ(w)dμ(w)dμ(v).
Thanks to μ being a probability density on V, we can rewrite the above expression as
12V×V(|ϕ(v)|2+|ϕ(w)|22ϕ(v)ϕ(w))dμ(w)dμ(v)0.
Hence the result. □

Next, by using the energy approach (i.e. by choosing appropriate multiplier), we prove an entropy inequality associated with the stationary model (5a)–(5b).

Lemma 2.

The solution fε(x,v) to the linear Boltzmann equation (5a)–(5b) satisfies the following entropy inequality:

ΩV|fε(x,v)|2dμ(v)dx+1ε2ΩV×Vσε(x)(fε(x,v)fε(x,w))2dμ(w)dμ(v)dxΩ|g(x)|2dx.

Proof.

Multiply (5a) by fε(x,v) and integrate over V yielding

V|fε(x,v)|2dμ(v)+12εVv·x|fε(x,v)|2dμ(v)+σε(x)ε2Vfε(x,v)Lfε(x,v)dμ(v)=Vg(x)fε(x,v)dμ(v).
Using Lemma 1 for the Dirichlet form and integrating over Ω yields:
ΩV|fε(x,v)|2dμ(v)dx+12εΩVv·x|fε(x,v)|2dμ(v)dx+Ωσε(x)2ε2V×V(fε(x,v)fε(x,w))2dμ(w)dμ(v)dx=ΩVg(x)fε(x,v)dμ(v)dx.
Consider the transport term in the previous expression and perform an integration by parts in the space variable yielding (with dΓ(x) as the surface measure on Ω):
Σ(v·n(x))|γfε(x,v)|2dμ(v)dΓ(x)=Σ+(v·n(x))|γ+fε(x,v)|2dμ(v)dΓ(x)
because of the zero absorption boundary condition (5b) on the incoming phase-space boundary. In the right hand side of the above expression, make the change of variables: vv:=Rxv for each xΩ, where Rxv=v2(v·n(x))n(x) is the reflection operator, yielding:
Σ+(v·n(x))|γ+fε(x,v)|2dvdΓ(x)=Σ(v·n(x))|γfε(x,v)|2dμ(v)dΓ(x)=0
again by the absorption boundary condition (5b). Hence the transport term does not contribute. We are left with the following expression:
ΩV|fε(x,v)|2dμ(v)dx+Ωσε(x)2ε2V×V(fε(x,v)fε(x,w))2dμ(w)dμ(v)dx=ΩVg(x)fε(x,v)dμ(v)dx.
Apply Cauchy–Schwarz inequality and Young’s inequality to the term involving the source term:
ΩVg(x)fε(x,v)dμ(v)dx(ΩV|g(x)|2dμ(v)dx)1/2(ΩV|fε(x,v)|2dμ(v)dx)1/212Ω|g(x)|2dx+12ΩV|fε(x,v)|2dμ(v)dx.
Using this inequality in the previous equality yields the entropy inequality. □

Our next result is to derive some uniform L2-estimates using the entropy inequality. We use the following notation:

h:=Vh(v)dμ(v)for all hL1(V;dμ).

Lemma 3.

Let fε(x,v) be the solution to (5a)–(5b). We have the following estimates:

(6a)fεL2(Ω×V;dxdμ)gL2(Ω),(6b)(fεfε)L2(Ω×V;dxdμ)εagL2(Ω),(6c)fεL2(Ω)gL2(Ω).

Proof.

The uniform estimate (6a) follows directly from the entropy inequality (Lemma 2).

Next, we focus on the uniform estimate (6b). We have

fε(x,·)fε(x)L2(V;dμ)2=V(V(fε(x,v)fε(x,w))dμ(w))2dμ(v).
Applying Cauchy–Schwarz we get:
V(fε(x,v)fε(x,w))dμ(w)(V(fε(x,v)fε(x,w))2dμ(w))1/2
because μ is a probability measure on V. Hence we have
Ωfε(x,·)fε(x)L2(V;dμ)2dxΩV×V(fε(x,v)fε(x,w))2dμ(w)dμ(v)dxε2agL2(Ω)2,
where the last inequality is because of the entropy inequality (Lemma 2).

Our next goal is to prove the uniform estimate (6c). Consider

fεL2(Ω)2=Ω(Vfε(x,v)dμ(v))2dx.
By Cauchy–Schwarz inequality, we have:
(Vfε(x,v)dμ(v))2(V|fε(x,v)|2dμ(v))(V12dμ(v))=V|fε(x,v)|2dμ(v)
because μ is a probability measure on V. Thus we have:
fεL2(Ω)fεL2(Ω×V;dxdμ)gL2(Ω).
Hence the result. □

Next, we prove a crucial estimate on the velocity average vfε. To begin with, we observe that the integral operator L is self-adjoint in L2(V;dμ), i.e.

Vψ(v)Lϕ(v)dμ(v)=Vϕ(v)Lψ(v)dμ(v)for all ϕ,ψL2(V;dμ).
This observation follows from the following successive equalities for the inner product in L2(V;dμ):
Vψ(v)Lϕ(v)dμ(v)=Vψ(v)ϕ(v)dμ(v)Vψ(v)Vϕ(w)dμ(w)dμ(v)=Vϕ(v)ψ(v)dμ(v)Vϕ(v)Vψ(w)dμ(w)dμ(v)=Vϕ(v)Lψ(v)dμ(v).
Next, we give a very important representation for the velocity variable in terms of the integral operator L.

Lemma 4.

For each ith component of the velocity variable, we have:

(7)vi=Lvi.

Proof.

Take ϕ(v)=vi and apply the integral operator L on to ϕ, i.e.

Lvi=viVwidμ(w)=vi
thanks to our assumption (1). □

Lemma 5.

Let fε(x,v) be the solution to (5a)–(5b). We have the following estimate:

(8)vfε[L2(Ω)]dε|v|2gL2(Ω).

Proof.

Consider the velocity average:

1εvfε=1εVvfε(x,v)dμ(v).
Crucial argument is to substitute for the velocity variable in the above expression using Lemma 4 which yields:
1εvfε=1εVfε(x,v)Lvdμ(v)=1εVvLfε(x,v)dμ(v)
because L is self-adjoint. Substituting for the Boltzmann operator in the above expression, we get
1εvfε=V×Vv1ε(fε(x,v)fε(x,w))dμ(w)dμ(v)(V×V|v|2dμ(w)dμ(v))1/2(V×V1ε2(fε(x,v)fε(x,w))2dμ(w)dμ(v))1/2,
where Cauchy–Schwarz inequality is used. Squaring the above inequality, we get
|1εvfε|2|v|2V×V1ε2(fε(x,v)fε(x,w))2dμ(w)dμ(v).
Integrate the above inequality on Ω yielding:
Ω|1εvfε|2dx|v|2ΩV×V1ε2(fε(x,v)fε(x,w))2dμ(w)dμ(v)dx|v|2gL2(Ω)2,
where we have used the entropy inequality (Lemma 2), thus proving the crucial estimate. □

Next, we prove uniform estimates on (vv)fε and its divergence.

Lemma 6.

Let fε(x,v) be the solution to (5a)–(5b). We have the following estimates:

(9a)(vv)fε[L2(Ω)]d×d|vv|2gL2(Ω),(9b)x·(vv)fε[L2(Ω)]dbgL2(Ω).

Proof.

Consider

(vv)fε[L2(Ω)]d×d2=Ω(V(vv)fε(x,v)dμ(v))2dx.
Cauchy–Schwarz inequality yields:
(V(vv)fε(x,v)dμ(v))2|vv|2(V|fε(x,v)|2dμ(v)).
Integrating the above inequality over Ω yields:
(vv)fε[L2(Ω)]d×d2|vv|2ΩV|fε(x,v)|2dμ(v)dx.
Using the uniform estimate (6a) from Lemma 3 yields the estimate (9a).

Next, we focus on the estimate (9b). To that end, multiply the stationary problem (5a) by the ith component of the velocity variable and integrate over V yielding:

εvifε+j=1dvivjfεxj+σε(x)εvifε=0.
The scattering coefficient is bounded in L(Ω). Hence the crucial estimate (8) of Lemma 5 would straightaway imply the following:
for each i{1,,d}j=1dvivjfεxjL2(Ω)bgL2(Ω),
thus proving the estimate (9b). □

Remark 2.

The result of Lemma 6 says that the matrix valued function (vv)fε is in the Hilbert space [H(div;Ω)]d, i.e. each row vector of (vv)fε belongs to

(10)H(div;Ω):={u[L2(Ω)]d such that x·uL2(Ω)}.

4.Moments method

We present a moments based approach to derive the limit behavior for (5a)–(5b) as ε0. This method is essentially borrowed from [7]. For readers’ convenience, we shall present this approach step-by-step.

Step I.

Integrate (5a) over V:

(11)fε+1ε·vfε=g(x).

Step II.

Multiply (5a) by the velocity variable v and integrate over V:

(12)εσε(x)vfε+1σε(x)·(vv)fε+1εvfε=0.

Step III.

Multiply (11) by a test function φ(x)H01(Ω) and integrate over Ω:

(13)1εΩvfε·φ(x)dx=Ωfε(x)φ(x)dxΩg(x)φ(x)dx.

Step IV.

Take dot product of the vector equation (12) with φ(x) and integrate over Ω:

(14)Ω(1σε(x)·(vv)fε)·φ(x)dx+1εΩvfε·φ(x)dx+εΩ1σε(x)vfε(x)·φ(x)dx=0.

Using (13) in (14) yields the following expression in which we need to pass to the limit.

(15)Ω(1σε(x)·(vv)fε)·φ(x)dx=Ωg(x)φ(x)dxΩfε(x)φ(x)dxεΩ1σε(x)vfε(x)·φ(x)dx.
The moments method culminates in passing to the limit as ε0 in (15) using some compactness properties of the family fε(x,v). In this article, this final step of the moments method is achieved in Section 5 for the one-dimensional case and in Section 6 for any arbitrary dimension. Observe that the expression (15) can be treated as a weak formulation for the following second-order differential equation in the x variable:
(16)x·(1σε(x)·(vv)fε)=Gε(x),
where Gε:ΩR is defined as
(17)Gε(x):=g(x)fε(x)+εx·(1σε(x)vfε(x)).
Next, we give a result showing that the family Gε(x) is uniformly bounded in L2(Ω) provided the exponent β takes values in certain interval of [0,).
Lemma 7.

Let fε(x,v) be the solution to (5a)–(5b) and suppose the exponent β2. Let Gε(x) be the family of scalar functions defined by (17). There exists a constant C, independent of ε, such that

GεL2(Ω)C.

Proof.

By chain rule, we have

(18)Gε(x)=g(x)fε(x)+εσε(x)x·vfε(x)εxσε(x)[σε(x)]2·vfε(x).
Substitute for x·vfε using the continuity equation (11) in the above expression yielding:
Gε(x)=g(x)fε(x)+ε2σε(x)g(x)ε2σε(x)fε(x)εxσε(x)[σε(x)]2·vfε(x)=(1+ε2σε(x))(g(x)fε(x))εxσε(x)[σε(x)]2·vfε(x).
Computing the L2-norm for Gε we have:
GεL2(Ω)(1+ε2σε(x))L(gL2(Ω)+fεL2(Ω))+εσε[σε]2Lvfε[L2(Ω)]dC(1+1a)(gL2(Ω)+fεL2(Ω))+C(ca2)ε1βvfε[L2(Ω)]d,
where we have used the growth assumption (4) on the heterogeneous coefficient σε. The assumption on the exponent β (i.e. β2) and the uniform estimates (6c) on fε help us arrive at the uniform L2-bound for Gε(x). □

5.One dimensional setting

In this section, we treat a special setting: both the spatial and velocity domains are one-dimensional, i.e. x(,+) and vV where V is either R or (1,+1). We consider the transport equation for the one particle distribution function fε(x,v):

(19a)fε+1εvdfεdx+σε(x)ε2(fεfε)=g(x)for (x,v)(,+)×V,(19b)fε(x,v)=0 for x=± and vV,
where the heterogeneous coefficient σε(x) are assumed to satisfy the same regularity assumptions as before, i.e. (4).
Theorem 1.

The solution family fε(x,v) to the one-dimensional stationary linear Boltzmann equation (19a)–(19b) exhibits the following compactness property:

fε(x,v)ρ(x)weakly in L2((,+)×V;dxdμ),
where ρ(x) is the unique solution to the stationary diffusion equation:
(20a)ρ(x)ddx(Dσdρdx)=g(x)for x(,+),(20b)ρ(x)=0for x=±,
where D is a constant equal to v2 and σ is the L weak- limit of the sequence σε.

Proof.

The a priori estimates of Lemma 6 in the one-dimensional setting imply that the sequence v2fε is uniformly bounded in H1(,+). Hence we can extract a sub-sequence such that

(21)v2fεDρstrongly in L2(,+),(22)ddxv2fεDdρdxweakly in L2(,+),
where D is a constant given by v2 and ρ(x) is the L2 weak limit of the local densities, i.e.
Vfε(x,v)dμ(v)ρ(x)weakly in L2(,+).
The second order differential equation (i.e. the one similar to (16)) that we get via the moments method (see Section 4) in this one-dimensional setting is the following:
(23)ddx(1σε(x)dv2fεdx)=Gε(x),
where Gε(x) is given by
Gε(x):=g(x)fε(x)+εddx(1σε(x)vfε(x)).
Define
(24)ζε(x):=1σε(x)ddxv2fε.
We have the following uniform L2-bound:
(25)ζεL2(,+)1addxv2fεL2(,+)bagL2(,+),
where we have used the assumption (4) that σε is bounded from below and the uniform a priori bound (9b) from Lemma 6. From (23), it follows that ζε(x) solves
(26)ddxζε(x)=Gε(x).
The uniform L2-estimate on Gε (Lemma 7) yields a uniform bound on the derivative of ζε:
(27)dζεdxL2(,+)=GεL2(,+)C.
The estimates (25) and (27) together imply that the sequence ζε is uniformly bounded in H1(Ω). The compact embedding H1(,+)L2(,+) implies that we can extract a sub-sequence and there exists a limit ζ0 such that
(28)ζεζ0strongly in L2(,+).
By the definition (24) of ζε, we have:
σε(x)ζε(x)=ddxv2fε.
Because of the strong convergence in (28), we have the following:
σεζεσζ0weakly in L2(,+),
where σ is the L weak- limit of the sequence σε. Identifying the above weak limit with the weak limit (22), we get:
ζ0(x)=Dσdρdx.
Upon passing to the limit (as ε0) in (26), we have the following equation:
dζ0dx=g(x)ρ(x).
Substituting for ζ0 in the above equation, we get:
ρ(x)ddx(Dσdρdx)=g(x)for x(,+).
The unique solvability of the limit equation (20a)–(20b) follows from the Lax–Milgram theorem. □

6.Arbitrary dimensions

In the previous section, using the moments method, we managed to prove the ε0 limit in the one-dimensional setting under the assumption that the exponent β2. In this section, we prove that the ε0 limit for the linear Boltzmann equation (5a)–(5b) in arbitrary dimensions can be obtained under some smallness criterion on the H12 norm of a quotient involving the heterogeneous coefficient σε(x).

Theorem 2.

Suppose there exists σ¯(x)L(Ω), bounded away from zero, such that

(29)σ¯(x)σε(x)σ¯(x)H1/2(Ω)=O(ε1+).
Then the family of local densities fε(x,·) exhibit the following compactness property:
Vfε(x,v)dμ(v)ρ(x)weakly in L2(Ω),
where the limit local density ρ(x) satisfies the following diffusion equation:
(30)ρ(x)·(Dσ¯(x)ρ(x))=g(x)for xΩ,(31)ρ(x)=0for xΩ
with D a constant matrix equal to vv.

Proof.

Let us rewrite our steady state model problem (5a) as follows:

εfε+v·fε+σ¯(x)ε(fεfε)=(σ¯(x)σε(x))(fεfεε)+εg(x).
Applying the moments approach outlined in Section 4 to the above problem, we arrive at the following weak formulation with a smooth test function Ψ(x) vanishing on Ω:
(32)εΩ1σ¯(x)vfε·Ψ(x)dx+Ω(1σ¯(x)·(vv)fε)·Ψ(x)dx+ΩfεΨ(x)dx=Ω(σ¯(x)σε(x)σ¯(x))vfεεdx+Ωg(x)Ψ(x)dx.
We have for the first term on the left hand side of the above equality:
|εΩ1σ¯(x)vfε·Ψ(x)dx|εσ¯1L(Ω)vfε[L2(Ω)]dΨ[L2(Ω)]dCε2,
thanks to Lemma 5. Next, for the first term on the right hand side of (32) we have
Ω(σ¯(x)σε(x)σ¯(x))vfεεdx=1ε(σ¯σεσ¯)Id,vfε[H1/2(Ω)]d,[H1/2(Ω)]d
which implies:
|Ω(σ¯(x)σε(x)σ¯(x))vfεεdx|1ε(σ¯σεσ¯)Id[H1/2(Ω)]dvfε[H1/2(Ω)]d.
The hypothesis (29) and the velocity-averaging Lemma [14,15] together imply that the above term is of O(ε+). With all the above observation, we are left to pass to the limit in the following expression:
O(ε2)+Ω(1σ¯(x)·(vv)fε)·Ψ(x)dx+ΩfεΨ(x)dx=O(ε+)+Ωg(x)Ψ(x)dx.
Observe that
limε0Ω(1σ¯(x)·(vv)fε)·Ψ(x)dx=limε0i,j=1dΩvivjfεxj(1σ¯(x)xiΨ(x))dx=i,j=1dΩDijρ(x)xj(1σ¯(x)xiΨ(x))dx,
where Dij=vivj. Thus, in the limit, we arrive at the following expression:
ΩDσ¯(x)ρ(x)·Ψ(x)dx+Ωρ(x)Ψ(x)dx=Ωg(x)Ψ(x)dx,
which is nothing but the weak formulation of the limit problem (30)–(31). The unique solvability of (30)–(31) follows again by the application of Lax–Milgram theorem. Hence the entire family converges to the limit local density. □

Next, we make an interesting observation on the smallness of the H12-condition (29) with regard to a rapidly oscillating periodic function. Remark that the smallness assumption (29) of the H12-norm in the periodic setting corresponds to having the exponent β>2.

Lemma 8.

Take σε(x)=2+sin(xεβ) with β>2 and take σ¯=2. Then

1εσ¯σεσ¯H12(0,2π)=O(ε+).

Proof.

As the denominator in the quotient is a constant, we shall ignore it for the calculations to follow. We shall compute the H1-norm by testing against a function φH1:

02πsin(xεβ)φ(x)dx=εβ02πx(cos(xεβ))φ(x)dx=εβ02πcos(xεβ)xφ(x)dxεβ(cos(εβ2π)φ(2π)cos(0)φ(0))
which implies that
sin(xεβ)H1(0,2π)εβ.
As we have
sin(xεβ)L2(0,2π)=O(1),
interpolating between H1 and L2 implies that
sin(xεβ)H12(0,2π)εβ2.
Hence we have
1εσ¯σεH12(0,2π)εβ21.
The choice of β>2 indeed implies that the H12-norm is of O(ε+). □

Even though the result of Lemma 8 is given in one dimension, the proof carries over to any arbitrary dimension.

Remark 3.

Note that the smallness condition on the H1/2-norm in Theorem 2 is quite strong as suggested by Lemma 8 which essentially says that any smooth function depending on the argument xεβ would satisfy the smallness assumption (29) provided β>2. Do note that we have treated the case β2 for the one-dimensional case in Theorem 1. If we were to suppose that the heterogeneous coefficient has the following structure

σε(x)=σ¯(x)+α(ε)hε(x),
i.e. the highly heterogeneous oscillations are of small amplitude, of size α(ε), which vanish in the ε0 limit then the limit procedure can be carried out without reverting to the condition (29). However, this is trivial as the heterogeneities are dying out in the ε0 limit. The smallness assumption (29) in Theorem 2 and the velocity averaging lemma (see Proof of Theorem 2) does allow us to have genuine heterogeneity in the coefficients σε(x).

Remark 4.

Even though the results presented so far – Theorem 1 and Theorem 2 – concern the stationary transport model, analogous results for the associated non-stationary model follow straightaway. Consider

tfε=Tεfε,fε(0)=fin(x,v)
with the operator
Tεfε:=1εv·xfε+σε(x)ε2(Vfε(t,x,w)dμ(w)fε(t,x,v)).
The following holds for the semigroup associated with Tε:
(pTε)1=0eptetTεdt.
Inverting the Laplace transform, we get
etTεfin=12πilimγiγ+iept(pTε)1findp.
Taking the ε0 limit in the previous expression, we get
limε0etTεfin=12πilimγiγ+ieptlimε0(pTε)1findp.
The asymptotic limit obtained in Theorem 2 for the resolvent helps us get
limε0etTεfin=12πilimγiγ+iept(pD)1ρindp=etDρin,
where we have used the following notations:
Du:=x·(Dσ¯(x)xu)andρin(x):=Vfin(x,v)dμ(v).

7.Concluding remarks

We have seen in Theorem 1 that we can get an explicit expression for the effective diffusion coefficient. This is analogous to the H-limits in one dimensional setting in the theory of H-convergence [20]. The theory of H-convergence, however, goes beyond the one-dimensional setting in getting explicit expressions while dealing with laminated materials. Our computations show that this is indeed the case in our setting. Results in this flavor will be in a later publication of the authors [8]. The main handicap of our result in the β2 case is that we are unable to handle dimensions higher than one. As noted in Remark 2, the estimates are in H(div;Ω) for the matrix-valued function (vv)fε. This is in stark contrast to the classical estimates used in the H-convergence theory with regard to elliptic problems. In a work which is in progress [8], the authors have made some progress in getting around the available less regularity information via constructing suitable class of test functions which emphasizes the importance of transport behavior at the scale of the microstructure. This approach employs the famous div-curl lemma (see [12] for a kinetic analogue of the div-curl lemma). Finally, it would be interesting to address this simultaneous limit in the case of linear Fokker–Planck equation. The authors shall address this problem in the near future.

Acknowledgements

This work was initiated during C.B.’s visit to the University of Cambridge. It was supported by the ERC grant matkit. The authors warmly thank Rémi Sentis for his constructive suggestions during the preparation of this article. The authors would also like to thank Thomas Holding for his valuable remarks on a preliminary version of this manuscript and Clément Mouhot for helpful suggestions. H.H. acknowledges the support of the ERC grant MATKIT and the EPSRC programme grant “Mathematical fundamentals of Metamaterials for multiscale Physics and Mechanic” (EP/L024926/1).

References

[1] 

G. Allaire, Homogenization and two-scale convergence, SIAM Journal on Mathematical Analysis 23(6) (1992), 1482–1518. doi:10.1137/0523084.

[2] 

G. Allaire, Shape Optimization by the Homogenization Method, Vol. 146, Springer Science & Business Media, 2002.

[3] 

G. Allaire and G. Bal, Homogenization of the criticality spectral equation in neutron transport, ESAIM: Modélisation Mathématique et Analyse Numérique 33(4) (1999), 721–746. doi:10.1051/m2an:1999160.

[4] 

G. Allaire and F. Golse, Transport et Diffusion, Lecture Notes, Ecole Polytechnique, 2013.

[5] 

G. Bal, N. Ben Abdallah and M. Puel, A corrector theory for diffusion-homogenization limits of linear transport equations, SIAM Journal on Mathematical Analysis 44(6) (2012), 3848–3873. doi:10.1137/110857301.

[6] 

C. Bardos, E. Bernard, F. Golse and R. Sentis, The diffusion approximation for the linear Boltzmann equation with vanishing scattering coefficient, Commun. Math. Sci. 13(3) (2015), 641–671. doi:10.4310/CMS.2015.v13.n3.a3.

[7] 

C. Bardos, F. Golse and D. Levermore, Fluid dynamic limits of kinetic equations. I. Formal derivations, Journal of Statistical Physics 63(1–2) (1991), 323–344. doi:10.1007/BF01026608.

[8] 

C. Bardos and H. Hutridurga, Div-curl approach for diffusion-homogenization limits for linear transport equations, 2016, In preparation.

[9] 

N. Ben Abdallah, M. Puel and M.-S. Vogelius, Diffusion and homogenization limits with separate scales, Multiscale Modeling & Simulation 10(4) (2012), 1148–1179. doi:10.1137/110828964.

[10] 

N. Ben Abdallah and M.-L. Tayeb, Diffusion approximation and homogenization of the semiconductor Boltzmann equation, Multiscale Modeling & Simulation 4(3) (2005), 896–914. doi:10.1137/040611227.

[11] 

A. Bensoussan, J.-L. Lions and G. Papanicolaou, Asymptotic Analysis for Periodic Structures, North-Holland, Amsterdam, 1978.

[12] 

L. Dumas and F. Golse, Homogenization of transport equations, SIAM Journal on Applied Mathematics 60(4) (2000), 1447–1470. doi:10.1137/S0036139997332087.

[13] 

F. Golse, Particle transport in nonhomogeneous media, in: Mathematical Aspects of Fluid and Plasma Dynamics, Springer, 1991, pp. 152–170. doi:10.1007/BFb0091366.

[14] 

F. Golse, P.-L. Lions, B. Perthame and R. Sentis, Regularity of the moments of the solution of a transport equation, Journal of Functional Analysis 76(1) (1988), 110–125. doi:10.1016/0022-1236(88)90051-1.

[15] 

F. Golse, B. Perthame and R. Sentis, Un résultat de compacité pour les équations de transport et application au calcul de la limite de la valeur propre principale de l’un opérateur de transport, C. R. Acad. Sci. Paris Sér. I Math 301(7) (1985), 341–344.

[16] 

T. Goudon and A. Mellet, Homogenization and diffusion asymptotics of the linear Boltzmann equation, ESAIM: Control, Optimisation and Calculus of Variations 9 (2003), 371–398. doi:10.1051/cocv:2003018.

[17] 

T. Goudon and F. Poupaud, Approximation by homogenization and diffusion of kinetic equations, 2001.

[18] 

F. Murat, Compacité par compensation, Annali della Scuola Normale Superiore di Pisa-Classe di Scienze 5(3) (1978), 489–507.

[19] 

F. Murat and L. Tartar, H-convergence, in: Séminaire d’Analyse Fonctionnelle et Numérique de L’Université D’Alger, Mimeographed Notes, 1978.

[20] 

F. Murat and L. Tartar, in: H-Convergence, Topics in the Mathematical Modelling of Composite Materials, Progr. Nonlinear Differential Equations Appl., Vol. 31, 1997, pp. 21–43. doi:10.1007/978-1-4612-2032-9_3.

[21] 

G. Nguetseng, A general convergence result for a functional related to the theory of homogenization, SIAM Journal on Mathematical Analysis 20(3) (1989), 608–623. doi:10.1137/0520043.

[22] 

R. Sentis, Approximation and homogenization of a transport process, SIAM Journal on Applied Mathematics 39(1) (1980), 134–141. doi:10.1137/0139011.