Effective diffusion on riemannian fiber bundlesThanks: Partially supported by CONACyT grant 135106.

Abstract.

The purpose of this paper is to provide equations to model the evolution of effective diffusion over a Riemannian fiber bundle (under the hypothesis of infinite diffusion rate along compact fibers). These equations are obtained by projecting the diffusion equation onto the base manifold of the fiber bundle. The projection (or dimensional reduction) is achieved by integrating the diffusion equation along the fibers of the bundle. This work generalizes an put into a general framework previous work on effective diffusion over channels and the interfaces between curved surfaces.

1. Introduction

Understanding spatially constrained diffusion is of fundamental importance in various sciences, such as biology, chemistry and nano-technology. However, solving the diffusion equation in arbitrarily constrained geometries is a very difficult task. One way to tackle it consists in reducing the degrees of freedom of the problem by considering only the main direction(s) of transport.

Figure 1.1. Region between to transversal cross section of a channel.

For example, the study of diffusion on thin channels can be carried out by reducing several spatial degrees of freedom to a single one by means of a projection method. More concretely, consider a process in a channel modeled by a density function $P$ that obeys the diffusion equation

$$\frac{\partial P}{\partial t}(x,t)=D_{0}\Delta P(x,t),$$

subject to the restriction that there is no density flow along the channel’s wall(s). We can construct an effective density function by letting

(1.1) $$\rho(u,t)=\lim_{h\mapsto 0}\frac{\hbox{Total\, concentration\, of\,\,}P\hbox{\, in\,\,}R(h)}{h},$$

where the region $R(h)$ is the section of the channel within two transversal cross sections that are an arc distance of $h$ apart over a base curve that “follows” the channel’s geometry (see Figure 1.1), and the variable $u$ is the arc-length parameter on this curve. It turns out that this effective density function $\rho$ obeys in an approximate manner an equation of the form (known as a generalized Fick-Jacobs equation)

(1.2) $$\frac{\partial\rho}{\partial t}(u,t)=\frac{\partial}{\partial u}\left(\sigma(u)\mathcal{D}(u)\frac{\partial}{\partial u}\left(\frac{\rho(u,t)}{\sigma(u)}\right)\right),$$

where $\sigma$ is given by

(1.3) $$\sigma(u)=\lim_{h\mapsto 0}\frac{\hbox{Area}(R(h))}{h}.$$

The function $\mathcal{D}$ is known as the effective diffusion coefficient and it encapsulates the effect of the channel’s geometry on the diffusion process along the base curve. Much work has been done (see [1, 8, 3, 4, 12, 6, 5, 7, 10, 11]) to find explicit formulas for $\mathcal{D}$ in terms of geometrical quantities associated to the channel, so that the Fick-Jacobs equation models the evolution of $\rho$ as closely as possible . We can distinguish two cases.

1. Infinite transverse diffusion rate. In this case it is assumed that the density function $P$ stabilizes instantly in the transversal directions of the channel. In mathematical terms this means that $P$ is constant along these transversal directions. This assumption results in an effective diffusion coefficient that depends on the curvature function of the base curve and $0$-th order geometrical quantities of the cross sections, such as width or area.

2. Finite transverse diffusion rate. In this case the finite time of transversal stabilization of $P$ is taken into account. This is characterized mathematically in that the resulting formulas for $\mathcal{D}$ involve the curvature function of the base curve, and tangential and curvature information of the channel’s wall(s).

The selection of the base curve is very important in the dimension-reduction technique described above. This is demonstrated by the fact (see [11, 10]) that if for channels of constant width the base curve is chosen properly, then the formulas for $\mathcal{D}$ coincide for the finite and infinite transversal diffusion rate cases.

Motivated by the above discussion, we know describe the main purpose of the paper. We develop a very general theory, in the infinite transversal diffusion rate case, for projecting the diffusion equation in a space of dimension $n$ to a base space of dimension $m$ with $m<n$. We do this in the context of the theory of fibre bundles. Such objects have a total space $E$, a base space $M$ and a projection map $\pi:E\rightarrow M$. In the case discussed above, the total space $E$ is the channel, $M$ is the base curve, and $\pi$ sends points on the transversal cross sections (fibers) to their base point in the base curve. The process of passing from the density $P$ to the reduced density $\rho$ is a particular case of a very well known construction: that of integrating a differentiable form along the fibers of the bundle. Using these tools we are able to give global and coordinate-free proofs of all our results. In this general setting, the effective diffusion $\mathcal{D}$ becomes an endomorphism of the tangent bundle of $M$, i.e for every $x$ in $M$ we have that $\mathcal{D}(x)$ is a linear map in the tangent space $T_{x}M$ of $M$. We will compute the effective diffusion using local frames instead of local coordinates, which results in a simpler and more geometric way of doing calculations. The paper is organized as follows

• In section 2, we show how to reduce the continuity equation in the fiber bundle $E$ to a reduced continuity equation in its base space $M$. We prove that if Fick’s law holds on $E$ (for a constant diffusion coefficient $D_{0}$) and if we have infinite diffusion rate in fiber direction, then the reduced continuity equation becomes a diffusion equation in $M$ (see Proposition 3). This last equation involves an effective diffusion coefficient $\mathcal{D}$, which is a bundle endomorphism of the tangent space of $M$.

• In section 3 we compute the effective diffusion for channels of constant width over arbitrary curves on the plane. We have obtained this result previously in [11], and our re-derivation of the formula serves as a test case of our general theory.

• In section 4 we compute the effective diffusion endomorphism $\mathcal{D}$ corresponding to the interface of two equidistant surfaces in 3-dimensional space. We do this by showing that the principal directions of the base surface are eigenvectors of $\mathcal{D}$, and then computing the eigenvalues of $\mathcal{D}$ in these directions. Using this result, we show that the calculation of $\mathcal{D}$ made by Ogawa in [9] for an elliptical cylinder, is just an approximation to ours obtained by only considering the two first terms in a series expansion of the function arctanh (the inverse of the hyperbolic tangent).

• In section 5 we compute the effective diffusion on the surface (not the interior) of a circular channel over an arbitrary curve in $\mathbb{R}^{3}$. The main point of the calculation is to illustrate how our techniques still apply to fibre bundles whose fibers are manifolds without boundary (in this case, circles).

• In the Appendix we make a brief review of the geometrical concepts needed for the construction of the effective diffusion endomorphism.

2. The effective diffusion equation

Let $E$ be fiber bundle over an $m$-dimensional manifold $M$, having compact fibers (with or without boundary) of dimension $k$, and projection map $\pi:E\rightarrow M$. We will assume that $M$ and $E$ are orientable, and with Riemannian metrics $<,>_{E}$ and $<,>_{M}.$

The continuity equation

The continuity equation on the fibre bundle $E$ is given by

(2.1) $$\frac{\partial P}{\partial t}+\hbox{div}(J)=0,$$

where the density function $P$ is a time dependent function on $E$, and the density flow $J$ is a time dependent vector field in $E$. The divergence of $J$ is given by

(2.2) $$\hbox{div}(J)=(*(d(*J^{\flat})))^{\sharp},$$

where $d$ is the differential operator acting on differentiable forms in $E$, $*$ is the Hodge star operator, and the $\flat$-operator converts vector fields to 1-forms (see the Appendix).

Integrating the continuity equation along the fibers

Since we are assuming that the fibers of $E$ are compact, by integrating along them we can define the operator $\pi_{*}$ sending $l$-forms in $E$ to $(l-k$)-forms in $M$. For any $l-$form $\omega$ on $E$ we define

$$\pi_{*}(\omega)_{x}(X_{1},\ldots,X_{l-k})=\int_{\pi^{-1}(x)}\beta_{\omega}$$

where

$$\beta_{\omega}(Y_{1},\ldots,Y_{k})=\omega(\tilde{X}_{1},\ldots,\tilde{X}_{l-k},Y_{1},\ldots,Y_{k})$$

and $\tilde{X}_{i}$ is a lift of $X_{i}$, i.e $D\pi(\tilde{X}_{i})=X_{i}$. Observe that $\tilde{X}_{i}$ is defined on the whole fiber $\pi^{-1}(x).$

By applying $*$ to 2.1 and using formula 2.2 we get

(2.3) $$\frac{\partial}{\partial t}(*P)+d(*J^{\flat})=0.$$

If we apply $\pi_{*}$ and then $*$ to the above equation we obtain

(2.4) $$\frac{\partial\rho}{\partial t}+*(\pi_{*}(d(*J^{\flat})))=0,$$

where

(2.5) $$\rho=*(\pi_{*}(*P))=\frac{\pi_{*}(P\mu_{E})}{\mu_{M}},$$

and $\mu_{E}$ and $\mu_{M}$ be the metric volume forms in $E$ and $M$.

Remark.

Formula 2.5 is a generalization of formula 1.1.

Neumann Boundary conditions

We will need adequate boundary condition on $J$ in order to write equation 2.4 as a continuity equation in $M$. If the fibers of $E$ are manifolds without boundary then $\pi_{*}$ commutes with $d$ (see [2, pg. 62]), i.e for any differential form $\omega$ we have that

(2.6) $$d(\pi_{*}\omega)=\pi_{*}(d\omega)$$

When the fibers of $E$ are manifolds with boundary we will assume that $\omega$ vanishes on all the vectors perpendicular to the boundary $\partial E$ of $E$. This last condition is equivalent to the assuming that there is no density flow across $\partial E$, i.e $J$ is parallel to the boundary of $E$. This last condition ensures that formula 2.6 still holds.

Integrating the density flow along the fibers

We will now use formula 2.6 to show that we can write 2.4 as a continuity equation in $M$. We need to find a time dependent vector field $j$ in $M$ such that

$$\hbox{div}(j)=*(\pi_{*}(d(*J^{\flat}))).$$

Using formula 2.6, and the definition of the divergence of $j$, we can write the above formula as

$$*(d(*j^{\flat}))=*(d(\pi_{*}(*J^{\flat}))).$$

This last equation is satisfied if we let

(2.7) $$j=(-1)^{m-1}(*(\pi_{*}(*J^{\flat})))^{\sharp},$$

where the $\sharp$-operator converts 1-forms to vector fields (see Appendix).

The effective continuity equation

We will refer to the time dependent function $\rho$ given by formula 2.5 as the effective density function, and to the time dependent vector field $j$ given by 2.7 as the effective density flow. We proved above that these objects satisfy the equation

(2.8) $$\frac{\partial\rho}{\partial t}+\hbox{div}(j)=0,$$

which we will refer to as the effective continuity equation.

The diffusion equation

Fick’s law establishes that

$$J(x,t)=-D(x)\nabla P(x,t),$$

where for $x$ in $E$ we have that $D(x)$ is a linear operator from $T_{x}E$ to $T_{x}E$, i.e an endomorphism of $TE_{x}$. The simplest choice of $D$ is to let it be scalar multiplication by a constant $D_{0}$. By using more general $D$’s we can model the inhomogeneity and anisotropy of $E$. Assuming Fick’s law, the continuity equation in $E$ becomes the diffusion equation

(2.9) $$\frac{\partial P}{\partial t}(x,t)=\hbox{div}(D(x)\nabla P(x,t)),$$

In our work we will always assume that the endomorphism $D$ on $TE$ is multiplication by a positive scalar $D_{0}$, in which case the diffusion equation becomes

$$\frac{\partial P}{\partial t}(x,t)=D_{0}\Delta P(x,t),$$

where the laplacian operator $\Delta$ applied to $P$ is

$$\Delta P=\hbox{div}(\nabla P).$$

Problem.

If Fick’s law holds in $E$ for $D=D_{0}$, does the effective density $\rho$ obeys a diffusion equation in $M$?

There is an important case when we can answer the above question positively (see Proposition 3 in the next paragraph).

An effective diffusion equation for infinite fiber diffusion rate

If the fibers of $E$ are “small enough” compared to the “size” of $M$, then it is to be expected that the density $P$ will stabilize faster along the fibers than along $M$. If we assume that this stabilization occurs infinitely fast, then $P$ must be constant along the fibers of $E$. We borrow the nomenclature from the physics literature, and refer to this situation by saying that there is an infinite fiber diffusion rate. Under this assumption, can write

$$P=\pi^{*}Q=Q\circ\pi$$

for a time dependent function $Q$ in $M$. Using 2.5 and the projection formula (see [2, pg. 63]) , we obtain

$$\rho=\frac{\pi_{*}(\pi^{*}(Q)\mu_{E})}{\mu_{M}}=\frac{Q\pi_{*}(\mu_{E})}{\mu_{M}},$$

which allows us to obtain $Q$ in terms of $\rho$ as

(2.10) $$Q=\frac{\rho}{\sigma},$$

where

(2.11) $$\sigma=\frac{\pi_{*}(\mu_{E})}{\mu_{M}}.$$

Observe that if $R$ is a region in $M$ then

(2.12) $$\int_{R}\sigma\mu_{M}=\int_{\pi^{-1}(R)}\mu_{E}=\hbox{vol}_{E}(\pi^{-1}(R)).$$

Remark.

Formula 2.11 is a generalization of that given by 1.3.

If Fick’s law holds in $E$ for $D=D_{0}$, then we have

$$J^{\flat}=-D_{0}d(\pi^{*}Q)=-D_{0}\pi^{*}(dQ),$$

where $\pi^{*}$ is the pull-back of forms under $\pi$. Hence, we can write the effective flow 2.7 as

(2.13) $$j=-(-1)^{m-1}D_{0}(*(\pi_{*}(*(\pi^{*}(dQ)))))^{\sharp}.$$

If we define $\mathcal{D}$ by

(2.14) $$\mathcal{D}=(-1)^{m-1}\left(\frac{D_{0}}{\sigma}\right)(\sharp\circ*\circ\pi_{*}\circ*\circ\pi^{*}\circ\flat),$$

where $\circ$ stands for composition of operators, then by combining formulas 2.10 and 2.13 we obtain

(2.15) $$j=-\sigma\mathcal{D}(\nabla Q)=-\sigma\mathcal{D}\left(\nabla\left(\frac{\rho}{\sigma}\right)\right).$$

The above formulas now lead to the following result.

Proposition 1.

For the case of infinite fiber diffusion rate, the effective density function $\rho$ satisfies the equation

(2.16) $$\frac{\partial\rho}{\partial t}(x,t)=\hbox{div}\left(\sigma(x)\mathcal{D}(x)\left(\nabla\left(\frac{\rho(x,t)}{\sigma(x)}\right)\right)\right),$$

where $\mathcal{D}:TX\rightarrow TX$ is the vector bundle morphism defined in formula 2.14.

Proof.

This follows from the effective continuity equation 2.8 and formula 2.15. As defined by formula 2.14, $\mathcal{D}$ sends vector fields in $M$ to vector fields in $M$. To prove that for any $x$ in $M$ it is in fact linear map $\mathcal{D}(x):T_{x}M\rightarrow T_{x}M$, we show that $\mathcal{D}(fX)=f\mathcal{D}(X)$ for any smooth function $f:M\rightarrow\mathbb{R}$ and vector field $X$ in $M$. We have that

$$(*\circ\pi^{*}\circ\flat)(fX)=(\pi^{*}f)(\pi*(X^{\flat})).$$

Using this and the projection formula (see [2, pg. 63]), we obtain that

$$\pi_{*}(*\circ\pi^{*}\circ\flat)(fX)=f\pi_{*}(\pi*(X^{\flat})),$$

and hence

$$\mathcal{D}(fX)=f\mathcal{D}(X).$$

∎

Definition 2.

We will refer to the endomorphism $\mathcal{D}$ of $TM$, given by formula 2.14, as the effective diffusion endomorphism.

Remark.

Equation 2.16 is a generalization of 1.2.

Proposition 3.

For the case of infinite fiber diffusion rate, there exists a metric in $M$ such that the effective density function $\rho$ satisfies the diffusion equation

(2.17) $$\frac{\partial\rho}{\partial t}(x,t)=\hbox{div}\left(\mathcal{D}(x)\nabla\rho(x,t)\right),$$

where $\mathcal{D}$ is the effective diffusion endomorphism. Furthermore, this choice of metric is such that for any region $R$ in $M$ we have that

$$\hbox{vol}_{M}(R)=\hbox{vol}_{E}(\pi^{-1}(R)).$$

Proof.

Let $<,>_{E}$ and $<,>_{M}$ be any metrics in $E$ and $M$. Define a new metric $<,>^{\prime}_{M}$ in $M$ by the formula

$$<X,Y>^{\prime}_{M}=\sigma^{2/m}<X,Y>_{M},$$

for $\sigma$ defined in 2.11. Then, since

$$\sigma^{2/m}=\left(\frac{\pi_{*}(\mu_{E})}{\mu_{M}}\right)^{2/m}$$

we have that

$$\sigma^{\prime}_{M}=\frac{\pi_{*}(\mu_{E})}{\mu^{\prime}_{M}}=\frac{\pi_{*}(\mu_{E})}{(\sigma^{2/m})^{m/2}\mu_{M}}=1.$$

Equation 2.17 then follows directly from Proposition 1, and the last part of the Proposition follows from equation 2.12. ∎

3. Constant-width channels on the plane

Figure 3.1. Lifting tangent vector $T$ to vector $(1-\kappa v)T$ along a fiber $\pi^{-1}(x)$

In this section we compute the effective diffusion function of a channel of constant width $w$ over a curve $C$ on the plane. Such a channel can be represented as the set

$$E=\{x+vN(x)|x\in C\hbox{\, and\,}-w/2\leq v\leq w/2\},$$

where $N$ is a unit normal field to $C$. For a reasonable curves (e.g compact) and small $w$, the space $E$ is a fibre bundle with projection map $\pi:E\rightarrow M$ given by $\pi(p)=x$ where

$$p=x+vN(x).$$

Let $T$ be a unit tangent field to $C$ that makes the frame $T,N$ positively oriented. To perform integration over the fibers of $E$ we need to compute the lift of $T$ to $E$. If we define

$$\mathcal{T}(p) = (1-\kappa(x)v)T(x)$$

$$\mathcal{N}(p) = N(x)$$

then $\mathcal{T}$ is such a lift of $T$ (see Figure 3.1), and $D\pi(\mathcal{N})=0$. Let $T^{*}$ be the dual field to $T$, and $\mathcal{T}^{*},\mathcal{N}^{*}$ be the dual frame to $\mathcal{T},\mathcal{N}$. The matrix of $<,>_{E}$ in the frame $\mathcal{T},\mathcal{N}$ is given by

$$g=\left(\begin{array}[]{cc}(1-\kappa v)^{2}&0\\ 0&1\end{array}\right).$$

Computing $\sigma$

The volume element in $E$ is

$$\mu_{E}=(1-\kappa v)\mathcal{T}^{*}\wedge\mathcal{N}^{*},$$

and hence

$$\pi_{*}(\mu_{E})=\left(\int_{-w/2}^{w/2}(1-\kappa v)dv\right)T^{*}=wT^{*}.$$

Since the volume form in $C$ is $T^{*},$ we conclude that

$$\sigma=\frac{wT^{*}}{T^{*}}=w.$$

Computing $\mathcal{D}$

Observe that

$$\pi^{*}(T^{*})=\mathcal{T}^{*},$$

and by formula 7.1 in the Appendix we have that

$$*(\mathcal{T}^{*})=g^{11}\det(g)^{1/2}\mathcal{N}^{*}=(1-\kappa v)^{-1}\mathcal{N}^{*}.$$

Hence

$$\pi_{*}(*(\pi^{*}(T^{*}))) = \int_{-w/2}^{w/2}(1-\kappa v)^{-1}dv$$

$$= \frac{1}{\kappa}\log\left(\frac{1+\kappa w/2}{1-\kappa w/2}\right)$$

Using formula 2.14 we obtain

(3.1) $$\mathcal{D} = \frac{D_{0}}{\kappa w}\log\left(\frac{1+\kappa w/2}{1-\kappa w/2}\right)$$

$$= \left(\frac{2D_{0}}{\kappa w}\right)\hbox{arctanh}(\kappa w/2).$$

We obtained this formula in [11] by different methods.

4. The interface between two equidistant surfaces in 3-d space

Let $S$ be an orientable surface in $\mathbb{R}^{3}$ and $N$ a unit normal field to this surface. For small $w>0$ and a “reasonable surface” the space

(4.1) $$E=\{x+vN(x)|x\in S\hbox{\,\, and\,\,\,}-w/2\leq v\leq w/2\}.$$

is a fibre bundle over $S$ with projection map $\pi:E\rightarrow M$ given by $\pi(p)=x,$ where

$$p=x+vN(x).$$

The surface $S$ has principal directions fields $T_{1}$ and $T_{2}$, with corresponding principal curvatures $\kappa_{1}$ and $\kappa_{2}$. If we define

$$\mathcal{T}_{1}(p) = (1-\kappa_{1}(x)v)T_{1},$$

$$\mathcal{T}_{2}(p) = (1-\kappa_{2}(x)v)T_{2},$$

$$\mathcal{N}(p) = N(x),$$

then $\mathcal{T}_{1}$ and $\mathcal{T}_{2}$ are lifts of $T_{1}$ and $T_{2}$, and $D\pi(\mathcal{N})=0$. The metric in $<,>_{E}$ is represented in the $\mathcal{T}_{1},\mathcal{T}_{2},\mathcal{N}$ frame by

$$g=\left(\begin{array}[]{ccc}(1-\kappa_{1}v)^{2}&0&0\\ 0&(1-\kappa_{2}v)^{2}&0\\ 0&0&1\end{array}\right).$$

We will let $\mathcal{T}_{1}^{*},\mathcal{T}_{2}^{*},\mathcal{N}^{*}$ be the dual frame to $\mathcal{T}_{1},\mathcal{T}_{2},\mathcal{N}$.

Computing $\sigma$

The volume form in $E$ is

$$\mu_{E}=(1-\kappa_{1}v)(1-\kappa_{2}v)\mathcal{T}_{1}^{*}\wedge\mathcal{T}_{2}^{*}\wedge\mathcal{N}^{*},$$

and hence

$$\pi_{*}(\mu_{E}) = \left(\int_{-w/2}^{w/2}(1-\kappa_{1}v)(1-\kappa_{2}v)dv\right)T_{1}^{*}\wedge T_{2}^{*},$$

$$= w(1+\kappa_{1}\kappa_{2}w^{2}/12)T_{1}^{*}\wedge T_{2}^{*}.$$

Since the volume form in $M$ is $T_{1}^{*}\wedge T_{2}^{*}$, we conclude from formula 2.11 that

$$\sigma=w(1+\kappa_{1}\kappa_{2}w^{2}/12).$$

Computing $\mathcal{D}$

Using formula 7.1 in the Appendix we obtain that

$$*(\mathcal{T}_{1}^{*}) = g^{11}\det(g)^{1/2}\mathcal{T}_{2}^{*}\wedge\mathcal{N}^{*}=\left(\frac{1-\kappa_{2}v}{1-\kappa_{1}v}\right)\mathcal{T}_{2}^{*}\wedge\mathcal{N}^{*},$$

$$(\mathcal{T}_{2}^{*}) = -g^{22}\det(g)^{1/2}\mathcal{T}_{1}^{*}\wedge\mathcal{N}^{*}=-\left(\frac{1-\kappa_{1}v}{1-\kappa_{2}v}\right)\mathcal{T}_{1}^{*}\wedge\mathcal{N}^{*}.$$

From these formulas and the identities $*T_{1}^{*}=T_{2}^{*},*T_{2}^{*}=-T_{1}^{*}$, we obtain

$$*(\pi_{*}(*(\mathcal{T}_{1}^{*}))) =- \left(\int_{-w/2}^{w/2}\left(\frac{1-\kappa_{2}v}{1-\kappa_{1}v}\right)dv\right)T_{1}^{*},$$

$$(\pi_{*}(*(\mathcal{T}_{2}^{*}))) =- \left(\int_{-w/2}^{w/2}\left(\frac{1-\kappa_{1}v}{1-\kappa_{2}v}\right)dv\right)T_{2}^{*}.$$

By evaluating the above integrals, using formulas

$$\pi*(T_{1}^{*})=\mathcal{T}_{1}^{*},\pi*(T_{2}^{*})=\mathcal{T}_{2}^{*},$$

and the fact that the matrix of the metric $<,>_{M}$ is the identity, formula 2.14 yields

$$\mathcal{D}(T_{1})=\mathcal{D}_{1}T_{1}\hbox{\,\, and\,\,\,}\mathcal{D}(T_{2})=\mathcal{D}_{2}T_{2}.$$

where

(4.2) $$\mathcal{D}_{1} = \frac{D_{0}\left(w\kappa_{1}\kappa_{2}+2(\kappa_{2}-\kappa_{1})\hbox{arctanh}(\kappa_{1}w/2)\right)}{\kappa_{1}^{2}w(1+\kappa_{1}\kappa_{2}w^{2}/12)},$$

(4.3) $$\mathcal{D}_{2} = \frac{D_{0}\left(w\kappa_{1}\kappa_{2}-2(\kappa_{2}-\kappa_{1})\hbox{arctanh}(\kappa_{2}w/2)\right)}{\kappa_{2}^{2}w(1+\kappa_{1}\kappa_{2}w^{2}/12)}.$$

We have just proved the following result.

Proposition 4.

Let $S$ be a surface in $\mathbb{R}^{3}$ with principal direction fields $T_{1}$ and $T_{2}$, and corresponding principal curvatures $\kappa_{1}$ and $\kappa_{2}$. For the bundle $E$ given by 4.1, we have that the fields $T_{1}$ and $T_{2}$ are eigenvectors of $\mathcal{D}$ with corresponding eigenvalues $\mathcal{D}_{1}$ and $\mathcal{D}_{2}$ given 4.2 and 4.3.

Remark.

At an umbilical point $x$ of $S$ (i.e where $\kappa_{1}(x)=\kappa_{2}(x))$ the above proposition is still valid for any pair of orthonormal vectors $T_{1},T_{2}$ in $T_{x}S$.

It is important to observe that if we want to write the effective diffusion equation 2.17 in coordinates, we need to express the principal direction fields $T_{1}$ and $T_{2}$ in terms of the corresponding coordinate fields. The reason for this is that the divergence operator that enters into the effective diffusion equation needs coordinates for its computation.

Remark.

We can express eigenvalues of $\mathcal{D}$ in terms of the gaussian and mean curvatures

$$K=\kappa_{1}\kappa_{2}\hbox{\,\, and\,\,\,}H=\frac{1}{2}(\kappa_{1}+\kappa_{2})$$

by using the identities

$$\kappa_{1}=H+\sqrt{H^{2}-K}\hbox{\,\, and\,\,\,}\kappa_{2}=H-\sqrt{H^{2}-K}.$$

We now discuss some applications of Proposition 4 to specific families of surfaces.

Spheres

In this case we have that the principal curvatures $\kappa_{1}$ and $\kappa_{2}$ satisfy

$$\kappa_{1}=\kappa_{2}=1/r,$$

where $r$ is the radius of the sphere. From formulas 4.2 and 4.3, we obtain that the eigenvalues of $\mathcal{D}$ are

(4.4) $$\mathcal{D}_{1}=\mathcal{D}_{2}=\frac{D_{0}12r^{2}}{12r^{2}+w^{2}}=\frac{D_{0}}{1+\frac{w^{2}}{12r^{2}}}.$$

Recovering the one dimensional case

Let $C$ be any curve on the plane with curvature function $\kappa$, and let $S=C\times\mathbb{R}\subset\mathbb{R}^{3}$. The principal fields of $S$ are the unit tangent $T$ to the curve and $(0,0,1)$, with corresponding eigenvalues $\kappa$ and $0$. Using formulas 4.2 and 4.3 we obtain

$$\mathcal{D}_{1} = \left(\frac{2D_{0}}{\kappa w}\right)\hbox{arctanh}(\kappa w/2),$$

$$\mathcal{D}_{2} = D_{0}.$$

Hence, the eigenvalue $\mathcal{D}_{1}$ of $\mathcal{D}$ coincides with the case of curves in the plane discussed in section 3 (see Formula 3.1).

We will now show how to write equation 2.16 in local coordinates. Consider the coordinates $(s,z)$ in $S$ where $s$ is the arc-length parameter of the curve $C$ and $z$ is standard $z$-coordinate in $\mathbb{R}^{3}$. The coordinate fields

$$\frac{\partial}{\partial s}\hbox{\,\, and\,\,\,}\frac{\partial}{\partial z}$$

are the principal direction fields of $S$, and

$$\mathcal{D}\left(\frac{\partial}{\partial s}\right) = \left(\frac{2D_{0}\hbox{arctanh}(\kappa w/2)}{\kappa w}\right)\frac{\partial}{\partial s},$$

$$\mathcal{D}\left(\frac{\partial}{\partial z}\right) = D_{0}\frac{\partial}{\partial z}.$$

Since the metric matrix is the identity, we have that

$$\nabla\rho=\frac{\partial\rho}{\partial s}\frac{\partial}{\partial s}+\frac{\partial\rho}{\partial z}\frac{\partial}{\partial z}.$$

We conclude that in this case formula 2.16 becomes

(4.5) $$\frac{\partial\rho}{\partial t}=\left(\frac{2D_{0}\hbox{arctanh}(\kappa w/2)}{\kappa w}\right)\frac{\partial^{2}\rho}{\partial s^{2}}+\frac{\partial^{2}\rho}{\partial z^{2}}.$$

If we use the series expansion

$$\hbox{arctanh}(x)=x+\frac{x^{3}}{3}+\frac{x^{5}}{5}+\ldots,$$

we obtain

$$\frac{2D_{0}\hbox{arctanh}(\kappa w/2)}{\kappa w}=1+\frac{w^{2}\kappa^{2}}{12}+\frac{w^{4}\kappa^{4}}{180}+\cdots$$

If we use only the first two terms of the above series in equation 4.5 we obtain

$$\frac{\partial\rho}{\partial t}=\left(1+\frac{w^{2}\kappa^{2}}{12}\right)\frac{\partial^{2}\rho}{\partial s^{2}}+\frac{\partial^{2}\rho}{\partial z^{2}}.$$

This last formula is that obtained by Ogawa in [9].

The torus

Figure 4.1. Eigenvalues $\mathcal{D}_{1}$ and $\mathcal{D}_{2}$ of the effective diffusion $\mathcal{D}$ for $D_{0}=1$, on a torus with inner radius $r=1$ , outer radius $R=2$ and width $w=1/4$.

For a torus with inner radius $r$ and outer radius $R$ we have that

$$\kappa_{1}=-1/r\hbox{\,\, and\,\,\,}\kappa_{2}=-\frac{\cos(\theta)}{R+r\cos(\theta)},$$

where $\theta$ is the variable parametrizing the parallels of the torus. In Figure 4.1 we show the graphs of $\mathcal{D}_{1}$ and $\mathcal{D}_{2}$ obtained by using the above values of $\kappa_{1}$ and $\kappa_{2}$ for specific values of $r,R$ and $w$.

5. Effective diffusion in the surface of a tube

Let $C$ be a curve in three dimensional space, and let $T,N,B$ be the corresponding Serret-Frenet frame. We will let $E$ be the set of points of the form (for a constant radius $r$)

$$p=x+r\cos(\theta)N+r\sin(\theta)B\hbox{\,\, where\,\,\,}{x\in C}.$$

We want to construct a lift of $T$ to $E$. To do this, consider $x=x(s),T=T(s),N=N(s)$ and $B=B(s)$ as functions of the arc length parameter $s$ of $C$. From the formula

$$\pi(p(s))=x(s)$$

we obtain

$$D\pi\left(\frac{dp}{ds}\right)=T\hbox{\,\, and\,\,\,}D\pi\left(\frac{dp}{d\theta}\right)=0,$$

where (using the Serret-Frenet formulas)

$$\frac{dp}{ds} = (1-\kappa r\cos(\theta))T+\tau r(-\sin(\theta)N+\cos(\theta)B),$$

$$\frac{dp}{d\theta} = -r\sin(\theta)N+r\cos(\theta)B,$$

for $\kappa$ and $\tau$ the curvature and torsion of $C$. Hence, if we define

$$\mathcal{T} = (1-\kappa r\cos(\theta))T,$$

$$\mathcal{N} = -r\sin(\theta)N+r\cos(\theta)B,$$

then $\mathcal{T}$ is lift of $T$ and $D\pi(\mathcal{N})=0$. The metric matrix in the $\mathcal{T},\mathcal{N}$ frame is

$$g=\left(\begin{array}[]{cc}(1-\kappa r\cos(\theta))^{2}&0\\ 0&r^{2}\end{array}\right).$$

Computing $\sigma$

The volume form in $E$ is

$$\mu_{E}=r(1-\kappa r\cos(\theta))\mathcal{T}^{*}\wedge\mathcal{H}^{*},$$

so that

$$\pi_{*}(\mu_{E})=\left(\int_{0}^{2\pi}r(1-r\kappa\cos(\theta))d\theta\right)T^{*}=2\pi rT^{*}$$

and hence

$$\sigma=\frac{2\pi rT^{*}}{T^{*}}=2\pi r.$$

Computing $\mathcal{D}$

We have that

$$\pi^{*}(T^{*})=\mathcal{T}^{*},$$

and

$$*\mathcal{T}^{*}=g^{11}\det(g)^{1/2}\mathcal{H}^{*}=r(1-\kappa r\cos(\theta))^{-1}\mathcal{H}^{*},$$

so that

$$\pi_{*}(*\mathcal{T}^{*})=\int_{0}^{2\pi}\left(\frac{r}{1-\kappa r\cos(\theta)}\right)d\theta=\frac{2\pi r}{1+r\kappa}\left(\sqrt{\frac{2}{1-r\kappa}-1}\right)T^{*}.$$

Hence

$$\mathcal{D} = \frac{D_{0}}{(1+r\kappa)}\left(\sqrt{\frac{2}{1-r\kappa}-1}\right).$$

$$= D_{0}\sqrt{\frac{1}{1-r^{2}\kappa^{2}}}$$

6. Conclusions

We have shown that the diffusion equation on the total space of a fiber bundle can be projected onto a diffusion equation on its base space, under the hypothesis of infinite diffusion rate along the fibers. We provided a general formula for the effective diffusion endomorphism of the reduced diffusion equation, that we later applied to obtain explicit formulas for diverse fiber bundles. Of particular interest was the computation of the effective diffusion endomorphism associated to the interface of two equidistant surfaces in 3-dimensional space, in terms of the principal curvatures of the base surface.

7. Appendix

The sharp and flat operators

We can see a 1-form $\alpha$ as a vector field $\alpha^{\sharp}$ defined by

$$<\alpha^{\sharp},X>=\alpha(X),$$

where $<,>$ is the metric of the space under consideration. Similarly, a vector field $X$ can be seen as a 1-form $X^{\flat}$ defined by

$$X^{\flat}(Y)=<X,Y>.$$

For a local frame $X_{1},\ldots,X_{n}$ the metric can be expressed as a symmetric matrix $g$ with coefficients

$$g_{ij}=<X_{i},X_{j}>,$$

and we will write

$$g^{ij}=(g^{-1})_{ij}.$$

Let $X^{1},\ldots,X^{n}$ be the 1-forms forming the dual frame to $X_{1},\ldots,X_{n}$, so that $X^{i}(X_{j})=\delta_{j}^{i}$. For a vector field

$$X=\sum_{i=1}^{n}a^{i}X_{i}$$

we have that

$$X^{\flat}=\sum_{i=1}^{n}a_{i}X^{i}\hbox{\,\, where\,\,\,}a_{i}=\sum_{j=1}^{n}g_{ij}a^{j}.$$

For a 1-form

$$\alpha=\sum_{i=1}^{n}\alpha_{i}X^{i}$$

we have that

$$\alpha^{\sharp}=\sum_{i=1}^{n}\alpha^{i}X_{i}\hbox{\,\, where\,\,\,}\alpha^{i}=\sum_{j=1}^{n}g^{ij}\alpha_{j}.$$

Hodge star operator

The Hodge star operator $*$ maps $l$-forms to $(n-l)$ forms, where $n$ is the dimension of the space under consideration, and it is defined so that for $l$-forms $\omega$ and $\eta$ we have that

$$\omega\wedge(*\eta)=<\omega,\eta>\mu,$$

where $\mu$ is the volume form of the metric. For monomials

$$\omega=\alpha_{1}\wedge\ldots\wedge\alpha_{l}\hbox{\,\, and\,\,\,}\eta=\beta_{1}\wedge\ldots\beta_{l},$$

where the $\alpha_{i}$’s and the $\beta_{j}$’s are 1-forms, we have that

$$<\omega,\eta>=\det(<\alpha_{i},\beta_{j}>)\hbox{\,\, where\,\,\,}<\alpha_{i},\beta_{j}>=<\alpha_{i}^{\sharp},\beta_{j}^{\sharp}>.$$

The $*$-operator satisfies the duality relation, where for an $l$-form $\omega$ we have that

$$**\omega=(-1)^{l(n-l)}\omega.$$

If for a local frame of vector fields $X_{1},\ldots,X_{n}$ we have the metric matrix $g_{ij}=\left\langle X_{i},X_{j}\right\rangle$, then the coefficients $g^{ij}$ of $g^{-1}$ are

$$g^{ij}=<X^{i},X^{j}>.$$

The metric volume form is given by

$$\mu=\det(g)^{1/2}X^{1}\wedge\ldots\wedge X^{n}.$$

Using the above formulas we obtain

$$*(X^{1}\wedge\ldots\wedge X^{n})=\det(g)^{-1/2},$$

and

(7.1) $$*X^{i}=\det(g)^{1/2}\sum_{j=1}^{n}g^{ij}\iota_{j}(X^{1}\wedge\ldots\wedge X^{n}),$$

where

$$\iota_{j}(X^{1}\wedge\ldots\wedge X^{n})=(-1)^{j+1}X^{1}\wedge\ldots\wedge X^{i-1}\wedge\hat{X^{j}}\wedge X^{i+1}\wedge\ldots\wedge X^{n},$$

and the hat in $\hat{X^{j}}$ indicates that that term has been removed as a factor in the above wedge product.