102 Chapter 3. The Boltzmann Equation in Bounded Domains
(2) Bounce-back reflection, in which the particles bounce back with reversed
velocity.
(3) Specular reflection, in which the particles bounce back specularly.
(4) Diffuse reflection (stochastic), in which the incoming particles are a proba-
bility average of the outgoing particles.
Due to its importance, there have been many contributions in the mathematical
study of different aspects of the Boltzmann boundary value problems [1], [2], [3],
[4], [6], [9], [10], [11], [15], [30], [31], [34], [37], [39], [46], [49], among others. See
also the references in the books [8], [12] and [44].
According to Grad [28, p. 243], one of the basic problems in the Boltz-
mann study is to prove existence and uniqueness of its solutions, as well as their
time-decay toward an absolute Maxwellian, in the presence of compatible phys-
ical boundary conditions in a general domain. In spite of those contributions to
the study of Boltzmann boundary problems, there are fewer mathematical re-
sults of uniqueness, regularity, and time decay-rate for Boltzmann solutions to-
ward a Maxwellian. In [41], it was announced that Boltzmann solutions near a
Maxwellian would decay exponentially to it in a smooth bounded convex domain
with specular reflection boundary conditions. Unfortunately, we are not aware of
any complete proof for such a result [45]. In [30], global stability of the Maxwellian
was established in a convex domain for diffusive boundary conditions. Recently,
important progress has been made in [16] and [47] to establish an almost ex-
ponential decay rate for Boltzmann solutions with large amplitude for general
collision kernels and general boundary conditions, provided that certain a-priori
strong Sobolev estimates can be verified. Even though these estimates had been
established for spatially periodic domains [22], [23] near Maxwellians, their va-
lidity is completely open for the Boltzmann solutions, even local in time, in a
bounded domain. As a matter of fact, this kind of strong Sobolev estimates may
not be expected for a general non-convex domain [23]. This is because even for
simplest kinetic equations with the differential operator v·∇
x
, the phase boundary
∂Ω×R
3
is always characteristic but not uniformly characteristic at the grazing set
γ
0
= {(x, v):x ∈ ∂Ωandv · n(x)=0} where n(x) is the outward normal at x.
Hence it is very challenging and delicate to obtain regularity from the general
theory of hyperbolic PDE. Moreover, in comparison with the half-space problems
studied, for instance in [34], [49], the geometrical complication makes it difficult
to employ spatial Fourier transforms in x. There are many cycles (bouncing char-
acteristics) interacting with the boundary repeatedly, and analysis of such cycles
is one of the key mathematical difficulties.
We aim to develop a unified L
2
−L
∞
theory in the near Maxwellian regime,
to establish exponential decay toward a normalized Maxwellian μ = e
−
1
2
|v|
2
, for all
four basic types of boundary conditions in rather general domains. Consequently,
uniqueness among these solutions can be obtained. For convex domains, these
solutions are shown to be continuous away from the singular grazing set γ
0
.