next up previous contents
Next: The nonexistence of hydrostatic Up: Some properties of the Previous: Some properties of the   Contents


Extrema of advected scalars

That no temperature field, being a solution of the steady-state energy equation (2.55), can possess a strong relative maximum or minimum at an interior point of its domain of existence follows from the essential elliptic nature of the equation and the absence of sources or sinks.

Theorem 1 (Nonexistence of extrema of advected scalars)   If a scalar field, $s$, is a regular solution of the steady-state advection-diffusion equation
\begin{displaymath}
\mbox{\boldmath$\nabla$}\cdot [-\mathsf{A}\cdot\mbox{\boldmath$\nabla$}s+\sum_{i=1}^N h_{(i)}(s)\mathbf{v}_{(i)}]=0
\end{displaymath} (2.76)

in some domain for some scalar functions $h_{(i)}$, solendoidal vector fields, $\mathbf{v}_{(i)}$, and symmetric, positive definite tensor field, $\mathsf{A}$; all continuously differentiable; then there are no interior strong relative extrema of $s$.

Proof: The idea for this proof, suggested by Prof. Bob Street (1999, pers. comm., 4 Feb.), is to recast the equation in quasilinear elliptic form, for which the result is known.

Carrying out the divergence,

\begin{displaymath}
\textsf{A}\raisebox{0.4ex}{\textbf{:}}\mbox{\boldmath$\nabla...
... \mathbf{v}_{(i)}\cdot h'_{(i)}(s)\mbox{\boldmath$\nabla$}s =0
\end{displaymath} (2.77)

or
\begin{displaymath}
\mathsf{A}\raisebox{0.4ex}{\textbf{:}}\mbox{\boldmath$\nabla...
...boldmath$\nabla$}s)+\mathbf{v}\cdot\mbox{\boldmath$\nabla$}s=0
\end{displaymath} (2.78)

where
\begin{displaymath}
\mathbf{v}\equiv\mbox{\boldmath$\nabla$}\cdot \mathsf{A}-\sum_{i=1}^N h'_{(i)}(s)\mathbf{v}_{(i)}.
\end{displaymath} (2.79)

In Cartesian tensor notation with the summation convention in force, this is
\begin{displaymath}
A_{jk}s_{,jk} + v_{j}s_{,j}=0,
\end{displaymath} (2.80)

which is of the form for which Hopf's Maximum Principle is shown to hold in treatises on partial differential equations (Courant & Hilbert 1962, pp. 320-8; Garabedian 1964, pp. 227-38).$\Box$


Alternative Proof: I offer here an original and quite different proof which I hope, by avoiding the artifice of a comparison function and using vectorial concepts rather than a general calculus of several variables, is more conducive to physical intuition. The use of vectors also ensures that the result is not tied to any particular coordinate system. The reader uncomfortable with a tensorial diffusivity may replace $\textsf{A}\cdot\mbox{\boldmath$\nabla$}s$ by $A\mbox{\boldmath$\nabla$}s$ wherever it appears; this is the special case of isotropy.

The proof is by contradiction: assume that there does exist an interior relative extremum. For definiteness, and without loss of generality, take this to be a minimum.

Construct a family of rays originating at the minimum and terminating when they encounter either:

(i)
a boundary point of the domain; or
(ii)
a stationary point, with respect to the ray, of $s$; i.e. $\mbox{\boldmath$\hat{r}$}\cdot\mbox{\boldmath$\nabla$}s=0$, where $\hat{r}$ is the unit radial vector from the minimum.

Except at the origin, and possibly the rays' termini, $s$ is strictly increasing along the rays:

\begin{displaymath}
\mbox{\boldmath$\hat{r}$}\cdot\mbox{\boldmath$\nabla$}s> 0,
\end{displaymath} (2.81)

by the definitions of a minimum and the rays (ii). Choose a value $s_1$ of $s$ between that at the minimum and the least of those at the rays' termini. Let $S$ be the set of points with $s=s_1$ passed through by the rays.

Each ray intersects $S$ exactly once, and, since $s$ possesses at least two continuous spatial derivatives (Gresho 1988), $S$ is closed and smooth enough to have a well-defined unit outward normal, $\hat{n}$. No ray is tangent to $S$, since then the ray should have terminated, by (ii); thus,

\begin{displaymath}
\mbox{\boldmath$\hat{r}$}\cdot\mbox{\boldmath$\hat{n}$}> 0.
\end{displaymath} (2.82)

Now, by definition of the vector triple product,

\begin{displaymath}
\mbox{\boldmath$\hat{r}$}\times(\mbox{\boldmath$\hat{n}$}\ti...
...{r}$}\cdot\mbox{\boldmath$\hat{n}$})\mbox{\boldmath$\nabla$}s,
\end{displaymath} (2.83)

but $\mbox{\boldmath$\hat{n}$}\times\mbox{\boldmath$\nabla$}s=0$, since the normal of a level surface is parallel to the gradient; therefore,
\begin{displaymath}
\mbox{\boldmath$\hat{n}$}\cdot(\textsf{A}\cdot\mbox{\boldmat...
...h$\nabla$}s\cdot(\textsf{A}\cdot\mbox{\boldmath$\nabla$}s) > 0
\end{displaymath} (2.84)

by (2.81), (2.82) and since $\textsf{A}$ is positive definite.

Integrate the steady-state advection-diffusion equation (2.76) over the volume $V$ enclosed by $S$:

\begin{displaymath}
\int\!\!\!\int\!\!\!\int _{V} \mbox{\boldmath$\nabla$}\cdot ...
...$}s
+\sum_{i=1}^N h_{(i)}(s) \mathbf{v}_{(i)}] \,\mathrm{d}V.
\end{displaymath} (2.85)

Applying the divergence theorem gives:
\begin{displaymath}
\sum_{i=1}^N\int\!\!\!\int _S \mbox{\boldmath$\hat{n}$}\cdot...
...}$}\cdot\mathsf{A}\cdot\mbox{\boldmath$\nabla$}s\,\mathrm{d}S,
\end{displaymath} (2.86)

of which the right hand side is positive by (2.84). The left hand side, however, vanishes;
$\displaystyle \int\!\!\!\int _S \mbox{\boldmath$\hat{n}$}\cdot h_{(i)}(s)\mathbf{v}_{(i)} \,\mathrm{d}S$ $\textstyle =$ $\displaystyle h_{(i)}(s_1)\int\!\!\!\int _S \mbox{\boldmath$\hat{n}$}\cdot\mathbf{v}_{(i)} \,\mathrm{d}S$ (2.87)
  $\textstyle =$ $\displaystyle h_{(i)}(s_1)\int\!\!\!\int\!\!\!\int _V \mbox{\boldmath$\nabla$}\cdot\mathbf{v}_{(i)} \,\mathrm{d}V=0;$ (2.88)

by virtue of the hypotheses on the $\mathbf{v}_{(i)}$.

This is a contradiction, so that the theorem is proved.$\Box$


Notes:


next up previous contents
Next: The nonexistence of hydrostatic Up: Some properties of the Previous: Some properties of the   Contents
Geordie McBain 2001-01-27