\part{0}
\pagerange{\pageref{firstpage}--\pageref{lastpage}}
\doi{S0963548301004989}

\begin{document}

\label{firstpage}
\maketitle

\begin{abstract}
This guide is for authors who are preparing papers for the
{\em Journal of Plasma Physics\/}
using the \LaTeXe\ document preparation system and the Cambridge
University Press JPP class file.
\end{abstract}

\tableofcontents

\section{Introduction}

The layout design for the {\em Journal of Plasma Physics\/} has been implemented
as a \LaTeXe\ class file. The JPP class is based on the ARTICLE class as
discussed in the \LaTeX\ manual  (2nd edition).
Commands which differ from the standard \LaTeX\ interface, or which are
provided in addition to the standard interface, are explained in this
guide. This guide is not a substitute for the \LaTeX\ manual itself.

Note that the final printed version of papers will use the Monotype Modern
typeface rather than the Computer Modern available to authors. Also,
the measure in JPP class is different from \LaTeXe\ article class. For
these reasons line lengths and page breaks will change and authors
should not insert hard breaks in their text.

\subsection{Introduction to \LaTeX}

The \LaTeX\ document preparation system is a special version of the
\TeX\ typesetting program.
 \LaTeX\ adds to \TeX\ a collection of commands which allow the author
to concentrate on the logical structure of the document rather than its
visual layout.

\LaTeX\ provides a consistent and comprehensive document preparation
interface.
\LaTeX\ can automatically number equations, figures,
tables, and list entries, as well as sections and
subsections.
Using this numbering system, bibliographic citations, page references
and cross-references to any other numbered entity (\eg
section, equation, figure, list entry) are quite straightforward.

\subsection{General style issues}

Use of \LaTeX\ defaults will result in a pleasing uniformity of layout
and font selection. Authors should resist the temptation to make
\emph{ad hoc} changes to these. Also avoid use of direct formatting unless
really necessary. Papers will be edited as usual, and this process may be
obstructed by the use of inserted line breaks, etc.

For general style issues, authors are referred to the `Preparation of
manuscripts' in the back cover of the journal. The language used in the
journal is British English, and spelling should conform to this.

Use should be made of symbolic references (\verb"\ref") in order to
protect against late changes of order, etc.
 
\subsection{Submission of \LaTeX\ articles}
 
Authors who intend to submit a \LaTeX\ article should obtain a copy of the
JPP class file. This is available by anonymous FTP from
\begin{verbatim}
  ftp.cup.cam.ac.uk
\end{verbatim}
You will find the class file and instructions contained in a
single zip archive \verb"jpp-cls.zip" in the directory
\begin{verbatim}
  pub/texarchive/journals/latex/jpp-cls
\end{verbatim}
The \verb"readme.txt" (which is the same directory) tells you how to
unzip the \verb"jpp-cls.zip" file.
There may also be an `\verb"unpacked"' directory containing all of the files
separately, in case of difficulty. If you cannot obtain the JPP files,
use the standard \verb"article" class, with the default `\verb"10pt"' option.
 
When submitting the final article, ensure that  the following are included and
are clearly labelled.
\begin{enumerate}
  \item A hardcopy printout of the article.
  \item The input file (exactly matching the hardcopy).
  \item A copy of any user-defined macros.
  \item If you have used \textsc{Bib}\TeX, the \verb".bib", \verb".bbl"
        and \verb".bst" files that were used.
  \item Any other files necessary to prepare the article for typesetting.
\end{enumerate}
The files for the \emph{final} article should be supplied on a PC or Macintosh
disk. Please write on the disk what format it is, as this saves time.
Submit the hardcopy and disk by post in the normal way.

\subsection{The JPP document class}

The use of document classes allows a simple change of style (or style option) to
transform the appearance of your document.
The CUP JPP class file preserves the standard \LaTeX\ interface such that any
document which can be  produced using the standard \LaTeX\ ARTICLE class
can also be produced with the JPP class.
However, the measure (or width of text) is slightly different from that
for ARTICLE; therefore line breaks will change and it is possible that
equations may need re-setting.


\section{Using the JPP class file}

First, copy the file \verb"jpp.cls" and \verb"upmath.sty" (see
\S\ref{sec:adpkg}) into the correct subdirectory on your system.
The JPP document style is implemented as a complete document class,
and {\em not\/} a document class option.
In order to use the JPP class, replace \verb"article" by \verb"jpp" in the
\verb"\documentclass" command at the beginning of your document:
\begin{verbatim}
  \documentclass{article}
\end{verbatim}
is replaced by
\begin{verbatim}
  \documentclass{jpp}
\end{verbatim}
Some of the standard document class options should be used.
Author-defined mac\-ros should be inserted before \verb"\begin{document}",
or in a separate file and should be included with the submission,
see \S\ref{secUser}. Authors must not change any of the macro definitions
or parameters in \verb"jpp.cls".

\subsection{Document class options}

In general, the following standard document class options should \emph{not} be
used with the JPP class file:
\begin{itemize}
  \item \texttt{10pt}, \texttt{11pt} and \texttt{12pt} -- unavailable;
  \item \texttt{twoside} is the default (\texttt{oneside} is disabled);
  \item \texttt{onecolumn} is the default (\texttt{twocolumn} is disabled);
  \item \texttt{titlepage} is not required and is disabled;
  \item \texttt{fleqn} and \texttt{leqno} should not be used, and are disabled.
\end{itemize}

The standard document class options should {\em not\/} be used with
the JPP class:
\begin{itemize}
  \item {\tt 10pt}, {\tt 11pt}, {\tt 12pt} -- unavailable.
  \item {\tt draft}, {\tt twoside} (no associated class file) -- {\tt twoside}
        is the default.
  \item {\tt fleqn}, {\tt leqno}, {\tt titlepage}, {\tt twocolumn} --
        unavailable.
\end{itemize}
\ifprodtf
The following new class options are provided:
\begin{itemize}
  \item \texttt{prodtf} -- tells the class file that we want to use the
    production typeface. This automatically resets the odd, even and top
    margins.
\end{itemize}
\fi

\subsection{Additional `packages' supplied with \texttt{jpp.cls}}\label{sec:adpkg}
The following additional package (\verb".sty") files are supplied in the
JPP distribution:
\begin{itemize}
  \item \verb"upmath" -- provides `upright' Greek math characters.
        Requires the AMS \verb"amsbsy" and \verb"amsgen" packages.
\end{itemize}

If your site does not have the AMS Fonts and \AMSLaTeX\ packages
installed, we strongly recommend that your site installs them. With them, you
can produce output which is much closer to the final result. The latest
AMS Fonts/\AMSLaTeX\ distributions can be found
on your nearest CTAN (Comprehensive \TeX\ Archive Network) site.

\subsubsection{The \texttt{UPMATH} package}

The \verb"upmath" package defines the macros \verb"\mathup" and
\verb"\mathbup", which allow access to the symbols in the AMS Euler fonts.

The \verb"upmath" package provides macros for upright lower-case Greek
(\verb"\ualpha"--\verb"\uxi"), and for bold lower-case Greek
(\verb"\ubalpha"--\verb"\ubxi"). The bold upright symbol \verb"\eta" has to
be treated differently, in this case use \verb"\uboldeta".
The \verb"upmath" package also provides \verb"\upartial" and \verb"\ubpartial".

To use the \verb"upmath" package, you need to have the AMS \verb"eurm/b"
fonts installed.

\section{Additional facilities}

In addition to all the standard \LaTeX\ design elements, the JPP class file
includes the following features:
\begin{itemize}
  \item Extended commands for specifying a short version
        of the title and author(s) for the running
        headlines.
  \item Abstract environment
  \item Control of enumerated lists.
  \item An extended \verb"theorem" environment, enabling you to also typeset
        unnumbered theorems, etc.
  \item Ability to create new theorem-like environments for claims,
        conjectures, examples, problems, remarks, etc.
  \item The \verb"subequations" and \verb"subeqnarray" environments for
        sub-numbering equations.
  \item A \verb"proof" environment.
  \item A \verb"acknowledgements" environment.
\end{itemize}
Once you have used these additional facilities in your document,
it can be processed only with \verb"jpp.cls".

\subsection{Titles, authors' names and running headlines}

In the JPP class, the title of the article and the author's name (or authors'
names) are used both at the beginning of the article for the main title and
throughout the article as running headlines at the top of every page.
The title is used on odd-numbered pages (rectos) and the author's name
(with initials only for first names) appears on even-numbered pages
(versos). The \verb"\pagestyle" and \verb"\thispagestyle" commands should
{\em not\/} be used.  Similarly, the commands \verb"\markright" and
\verb"\markboth" should not be necessary.

Although the article title can run to several lines of text,
the running headline must be a single line.
Moreover, the title can also incorporate new line commands
(\eg \verb"\\") but these are not acceptable in a running headline.
To enable you to specify an alternative short title and author's name, the
standard \verb"\title" and \verb"\author" commands have been extended to take
an optional argument to be used as the running headline:
\begin{verbatim}
  \title[A short title]{The full title which can be as long
                         as necessary}
  \author[Initials and last names of all authors]{The full names of
     all the authors, using letterspacing and \and before
     the last name in the list}
  \affiliation{Author's affiliation}
\end{verbatim}

Unlike most other document classes, \verb"jpp.cls" does not follow the
convention of using the second line of the \verb"\author" command to
write the author's affiliation.  Rather, this must be entered using
the \verb"\affiliation" command as shown above.  The \verb"\and"
command does not produce separate author/affiliation pairs, but
only generates the word `and' in small caps as required
by the JPP class.  Use the predefined macro \verb"\ls" to get the
letterspacing between the letters of the author's name.  If there is
more than one affiliation, use \verb"\\[\affilskip]" between lines of the
\verb"\affiliation", inserting footnote numbers manually.

An example will make all this clearer.  To produce

\medskip
\begin{center}
{\Large\bf Interesting title} \par\medskip
{\large A\ls L\ls A\ls N\ns N.\ns O\ls T\ls H\ls E\ls R$^1$\ns
  \and\ns J.\,Q.\ns P\ls U\ls B\ls L\ls I\ls C$^2$}\par
\medskip
{\small $^1$University of Moscow, Moscow, Russia\\[3pt]
$^2$Ngonga University, Nairobi, Kenya}\par
\medskip
{\small (Received 12 June 1992)}
\end{center}
\medskip

\noindent you would type
{\small\begin{verbatim}
  \title[Short title]{Interesting title}
  \author[A. N. Other and J. Q. Public]
         {A\ls L\ls A\ls N\ns N.\ns O\ls T\ls H\ls E\ls R$^1$\ns
          \and\ns J.\,Q.\ns P\ls U\ls B\ls L\ls I\ls C$^2$}
  \affiliation{$^1$University of Moscow, Moscow, Russia\\[\affilskip]
               $^2$Ngonga University, Nairobi, Kenya}
  \date{12 June 1992}
\end{verbatim}}
The \LaTeX\ \verb"\thanks" command can be used inside the \verb"\author" command
to insert an authors `Current address' information as a footnote, appearing at the
bottom of the title page.

\subsection{Abstract}

The JPP class provides for an abstract, produced by
\begin{verbatim}
  \begin{abstract}
  ...
  \end{abstract}
\end{verbatim}

This should appear just \emph{before} the first \verb"\section" command.

\subsection{Lists}

The JPP style provides the three standard list environments plus an additional
unnumbered list:
\begin{itemize}
  \item Numbered lists, created using the \verb"enumerate" environment.
  \item Bulleted lists, created using the \verb"itemize" environment.
  \item Labelled lists, created using the \verb"description" environment.
\end{itemize}
The \verb"enumerate" environment numbers each list item with an italic
letter in parentheses;
alternative styles can be achieved by inserting a redefinition of the
number labelling command after the \verb"\begin{enumerate}". For example, a
list numbered with roman numerals inside parentheses can be produced by the
following commands:
\begin{verbatim}
  \begin{enumerate}
  \renewcommand{\theenumi}{\roman{enumi}}
  \item first item
         :
  \end{enumerate}
\end{verbatim}
This produces the following list:
\begin{enumerate}
  \renewcommand{\theenumi}{\roman{enumi}}
  \item first item
  \item second item
  \item \etc
\end{enumerate}

In the last example, the item alignment is uneven because
the standard list indention is designed to be sufficient for arabic
numerals rather than the wider roman numerals. In order to enable
different labels to be used more easily, the \verb"enumerate" environment
in the JPP style can be given an optional argument which (like a standard
\verb"thebibliography" environment) specifies the \emph{widest label}. For
example,
\begin{enumerate}[(iii)]
\renewcommand{\theenumi}{\roman{enumi}}
  \item first item
  \item second item
  \item \etc
\end{enumerate}
was produced by the following input:
\begin{verbatim}
  \begin{enumerate}[(iii)]
  \renewcommand{\theenumi}{\roman{enumi}}
  \item first item
          :
  \end{enumerate}
\end{verbatim}

Once you have used the optional argument on the \verb"enumerate"
environment, do not process your document with a standard \LaTeX\ class file.

\subsection{Theorem environments}\label{sectTheor}

The \verb"\newtheorem" command works as described in the \LaTeX\ manual,
but produces spacing and caption typefaces required to the JPP style.  The
preferred numbering scheme is for theorems to be numbered within sections,
as 1.1, 1.2, 1.3, etc., but other numbering schemes are permissible and
may be implemented as described in the \LaTeX\ manual.  In order to allow
authors maximum flexibility in numbering and naming, {\em no\/} theorem-like
environments are defined in \verb"jpp.cls".  Rather, you have to define
each one yourself.  Theorem-like environments include Theorem, Definition,
Lemma, Corollary, and Proposition.  For example,
\begin{verbatim}
  \newtheorem{lemma}{Lemma}[section]
  \renewcommand{\thelemma}{\Roman{lemma}}
\end{verbatim}

The new macro which follows can also produce numbered and unnumbered
environments.

\subsection{Definitions}

The \verb"\newdefinition" command can be used to create environments for
claims, conjectures, examples, problems, remarks, etc. These are typeset
in the way as theorems, except the text is typeset in roman instead of italic.

You can have an environment created by \verb"\newdefinition" number in the
same sequence as a `Theorem' by adding the optional \verb"[theorem]" argument
as below:
\begin{verbatim}
  \newdefinition{definition}[theorem]{Definition}
  :
  \begin{definition}
  This is a definition.
  \end{definition}
  \begin{definition*}
  This is an unnumbered definition.
  \end{definition*}
\end{verbatim}
which produces the following:
\begin{definition}
This is a definition.
\end{definition}
\begin{definition*}
This is an unnumbered definition.
\end{definition*}

\subsection{Proof environment}

The standard \LaTeX\ constructs do not include a proof environment to
follow a theorem, lemma etc. (see also $\S$\ref{sectTheor}), and so one has
been added for the JPP style. 

Note the use in the following examples of an optional argument in square
braces which may contain any information you may wish to add.
For example,
\begin{verbatim}
  \begin{theorem}[2]\label{thenv}
  Let the scalar function $T(x,y,t,\bmath{\omega})$ be a conserved
  density for solutions of {\rm(9)}. Then the two-component function
  \begin{equation}
  \hbox{\boldit P} = \mathsf{J}{\cal E} T
  \end{equation}
  represents the infinitesimal generator of a symmetry group
  for {\rm (9)}.
  \end{theorem}
  \begin{proof}[Proof of Theorem~\ref{thenv}]
  The assumption about $T$ means that
  \[
  0 \sim \frac{\partial T}{\partial t} +
  {\cal E}T \bmath{\omega}_t
  = \frac{\partial T}{\partial t} + \{ T, H\},
  \]
  where $\partial T / \partial t$ refers to explicit
  dependence on $t$. The skew symmetry of $\mathsf{J}$ hence implies
  \begin{equation}
  \frac{\partial T}{\partial t} \sim  \{ H, T\},
  \end{equation}
  whereupon the operation $\mathsf{J}\cal E$, which commutes with
  $\partial_t$ in its present sense, gives
  \[
  \frac{\partial \hbox{\boldit P}}{\partial t}
  = \mathsf{J}{\cal E}\{ H, T\}.
  \]
  This equation reproduces the characterisation of symmetries that
  was expressed by (19), thus showing P to represent a symmetry
  group.
  \end{proof}
\end{verbatim}
produces the following text:
   \begin{theorem}[2]\label{thenv}
   Let the scalar function $T(x,y,t,\bmath{\omega})$ be a conserved
   density for solutions of {\rm(9)}. Then the two-component function
   \begin{equation}
   \hbox{\boldit P} = \mathsf{J}{\cal E} T
   \end{equation}
   represents the infinitesimal generator of a symmetry group for {\rm (9)}.
   \end{theorem}
   \begin{proof}[Proof of Theorem~\ref{thenv}]
   The assumption about $T$ means that
   \[
   0 \sim \frac{\partial T}{\partial t} +
   {\cal E}T \bmath{\omega}_t
   = \frac{\partial T}{\partial t} + \{ T, H\},
   \]
   where $\partial T / \partial t$ refers to explicit
   dependence on $t$. The skew symmetry of $\mathsf{J}$ hence implies
   \begin{equation}
   \frac{\partial T}{\partial t} \sim  \{ H, T\},
   \end{equation}
   whereupon the operation $\mathsf{J}\cal E$, which commutes with
   $\partial_t$ in its present sense, gives
   \[
   \frac{\partial \hbox{\boldit P}}{\partial t} = \mathsf{J}{\cal E}\{ H, T\}.
   \]
   This equation reproduces the characterisation of symmetries that
   was expressed by (19), thus showing P to represent a symmetry group.
   \end{proof}

The final \usebox{\proofbox} will not be included if the \verb"proof*" environment is
used.

\ifprodtf
\smallskip
Note: If a proof environment is ended by and unnumbered equation, it is
customary for Editors to mark the proof box to be moved up, so that it
aligns with the end of the equation.  This can be achieved by placing a
\verb"\raiseproofboxby" command before the \verb"\end{proof}".
A value of \verb"1.5\baselineskip" seems to give the best results. \eg
\begin{verbatim}
  \begin{proof} \raiseproofboxby{1.5\baselineskip}
    The assumption about $T$ means that
    :
    \[
      \frac{\upartial T}{\upartial t} \sim \{ H, T\},
    \]
  \end{proof}
\end{verbatim}
\fi

\section{Mathematics and units}

The JPP class file will insert the correct space above and below
displayed maths if standard \LaTeX\ commands are used; for example use
\verb"\[ ... \]" and {\em not\/} \verb"$$ ... $$". Do not leave blank
lines above and below displayed equations unless a new paragraph is
really intended.

\subsection{Numbering of equations}

The \verb"subequations" and \verb"subeqnarray" environments have been
incorporated into the JPP class file. Using these two environments, you can
number your equations (1.1{\it a}), (1.1{\it b}) \etc automatically. For
example, you can typeset
  \begin{subequations}
  \begin{equation}
    a_1 \equiv (2\Omega M^2/x)^{\frac{1}{4}} y^{\frac{1}{2}}
  \end{equation}
  and
  \begin{equation}
    a_2 \equiv (x/2\Omega)^{\frac{1}{2}}k_y/M.
  \end{equation}
  \end{subequations}
by using the \verb"subequations" environment as follows:
\begin{verbatim}
  \begin{subequations}
  \begin{equation}
    a_1 \equiv (2\Omega M^2/x)^{\frac{1}{4}} y^{\frac{1}{2}}
  \end{equation}
  and
  \begin{equation}
    a_2 \equiv (x/2\Omega)^{\frac{1}{2}}k_y/M.
  \end{equation}
  \end{subequations}
\end{verbatim}
You may also typeset an \verb"array" such as:
  \begin{subeqnarray}
    \dot{X}    & = & \gamma X - \gamma\delta\eta ,\\
    \dot{\eta} & = & {\textstyle\frac{1}{2}} \delta + 2X\eta .
  \end{subeqnarray}
by using the \verb"subeqnarray" environment as follows:
  \begin{verbatim}
  \begin{subeqnarray}
    \dot{X}    & = & \gamma X - \gamma\delta\eta ,\\
    \dot{\eta} & = & {\textstyle\frac{1}{2}} \delta + 2X\eta .
  \end{subeqnarray}
  \end{verbatim}

\subsection{AMS fonts -- especially if you do not have them}

If you need AMS symbols but do not have the AMS fonts you can ensure that
they will be correctly typeset by taking the following steps. Set
up user-defined macros that can be redefined by the typesetter to use
the correct AMS macros. For example, the blackboard bold symbols, sometimes
called shell or outline characters, are obtained with the AMS macro
\verb"\mathbb{..}". Instead, use a macro definition such as:
\begin{verbatim}
  
   \newcommand\BbbE{\ensuremath{\mathsf{E}}}
\end{verbatim}
This substitutes a sans serif character where you want blackboard bold.
You can typeset the input file and the typesetter is alerted to do the
substitution.

 The following example (which uses the \verb"\providecommand" macro) will work
without modification by the typesetter, because the \verb"\providecommand"
macro will not overwrite any existing \verb"\mathbb" definition.
\begin{verbatim}
  \providecommand\mathbb[1]{\ensuremath{\mathsf{#1}}}
  ...
  \newcommand\BbbE{\mathbb{E}}
\end{verbatim}

Plain \TeX\ provides only \verb"\leq" and \verb"\geq" which typeset the
Computer Modern symbols $\realleq$ and $\realgeq$, respectively. These will be
redefined at typesetting to use the AMS symbols \verb"\geqslant" and
\verb"\leqslant", to give the slanted symbols.

If you wish to use AMS fonts with \LaTeXe\ you must be using at least
version 2.0. Earlier versions are not supported.

\subsection{Typeface}

Sometimes, non-italic symbols are required in maths. This section describes
how these may be obtained using \LaTeX\ and \verb"jpp.cls".

\subsubsection{Roman symbols}\label{roman}

The mathematical operators and constants, such as sin, cos, log and
exp (and many others) are covered by \LaTeXe\ macros which ensure that they are typeset in
roman text, even in math mode: \verb"\sin", \verb"\cos", \verb"\log",
\verb"\exp". Where single letters are concerned (\eg d, i, e) just use
\verb"\mathrm" in math mode, i.e\ \verb"$E=m\mathrm{c}^2$" which typesets as
$E=m\mathrm{c}^2$, giving the correct roman character but with math spacing.
When the term involves more than one character (\eg \Real\ or \Imag )
text-character spacing is required:
\begin{verbatim}
  $\mbox{Re}\;x$
\end{verbatim}
which typesets as $\Real\;x$.

Where such expressions are used repeatedly, macro definitions can reduce
typing and editing. The following examples are included in the preamble
of the input files for this document, \verb"jppguide.tex", and of the
sample article \verb"jppsampl.tex". Authors are encouraged to use them
and others like them.
\begin{verbatim}
  \newcommand\Real{\mbox{Re}}
  \newcommand\Imag{\mbox{Im}}
  \newcommand\Ai{\mbox{Ai}}
  \newcommand\Bi{\mbox{Bi}}
\end{verbatim}

\subsubsection{Multiletter italic symbols}

If multiletter symbols are used in maths mode, for example Reynolds,
Prandtl numbers, \etc the standard maths mode spacing between them is
too large and text-character spacing is required. As described
in \S\ref{roman} (but here for italic letters) use for example
\begin{verbatim}
  \newcommand\Rey{\mbox{\textit{Re}}}
  \newcommand\Pran{\mbox{\textit{Pr}}}
\end{verbatim}

\subsubsection{Upright and sloping Greek symbols}


Like the italic/roman fonts described above, Greek letters for math
variables are printed in the Journal in sloping type, whereas constants,
operators and units are upright. However, \LaTeX\ normally produces sloping
lowercase Greek and upright capitals. The upright lowercase required for
example for `pi' and `mu' (micro) are provided by the \verb"upmath" package
(from the AMS Euler fonts). If you don't have the AMS font set, you will have
to use the normal math italic symbols and the typesetter will substitute the
corresponding upright characters.

You will make this much easier if you can use the macros \verb"\upi" and
\verb"\umu" in your text to indicate the need for the upright characters,
together with the following definitions in the preamble, where they can
easily be redefined by the typesetter.
\begin{verbatim}
  \providecommand\upi{\pi}
  \providecommand\umu{\mu}
\end{verbatim}

Notice the use of \verb"\providecommand", this stops your
definitions overwriting the ones used by the typesetter.

Sloping uppercase Greek characters can be obtained by using the
\verb"\varGamma" \etc macros (which are built into the class file).
The few required upright Greek capital letters (\eg capital delta for
`difference') are given by the relevant \LaTeX\ macro, \eg \verb"\Delta".

\newcommand\lra{\ensuremath{\quad\longrightarrow\quad}}

\subsubsection{Sans serif symbols}

The \verb"\textsf" and \verb"\mathsf" commands change the typeface to
sans serif, giving upright characters. Occasionally, bold-sloping sans serif
is needed. You should use the following supplied macros to obtain
these fonts.\\[6pt]
\verb"  \textsf{text}  " \lra \textsf{text}
  \qquad \verb"\mathsf{math}  " \lra $\mathsf{math}$\\
\verb"  \textsfi{text} " \lra \textsfi{text}
  \qquad \verb"\mathsfi{math} " \lra $\mathsfi{math}$\\
\verb"  \textsfb{text} " \lra \textsfb{text}
  \qquad \verb"\mathsfb{math} " \lra $\mathsfb{math}$\\
\verb"  \textsfbi{text}" \lra \textsfbi{text}
  \qquad \verb"\mathsfbi{math}" \lra $\mathsfbi{math}$\\[6pt]
You can use them like this:
\begin{verbatim}
  \newcommand\ssC{\mathsf{C}}    
  \newcommand\sfsP{\mathsfi{P}}  
  \newcommand\sfbsX{\mathsfbi{X}} 
\end{verbatim}
Note that the bold-slanted font \verb"\testsfbi" and \verb"\mathsfbi"
use the slanted sans serif font \verb"cmssi" -- because there is no 
bold-slanted maths sans serif font in Computer Modern! If you use the
supplied sans-serif text and math commands the typesetter will be
able to substitute the fonts automatically.

\subsubsection{Bold math italic symbols}
 
If you require bold math italic symbols/letters, \LaTeX\ provides several ways
of getting them.  Firstly, \LaTeX's \verb"\boldmath" switch which can be
used like this:
\begin{verbatim}
  $ \mbox{\boldmath $P$} = \mathsf{J}\mathcal{E} T $
\end{verbatim}
As you can see it takes quite alot of typing to achieve a bold math italic P.
Another problem is that you can't use \verb"\boldmath" in math mode, thus the
\verb"\mbox" forces the resulting text into text style (which you may not
want).

You can cut down on the typing by defining the following macro (which
is in the preamble of this guide):
\begin{verbatim}
  \providecommand\boldsymbol[1]{\mbox{\boldmath $#1$}}
\end{verbatim}
Then for the above example you can define:
\begin{verbatim}
  \newcommand\bldP{\boldsymbol{P}}
\end{verbatim}
the achieve the same result. This still doesn't remove the default to
text style problem.

If you have the \verb"amsbsy" package on your system, then you can remove
this limitation as well by placing a \verb"\usepackage{amsbsy}" after the
\verb"\documentclass" line in your document, and then use \verb"\boldsymbol"
for bold math italic symbols (don't define \verb"\boldsymbol" yourself when
using \verb"amsbsy").

\subsubsection{Script characters}

Script characters should be typeset using \LaTeXe's \verb"\mathcal"
command. This produces the Computer Modern symbols such as $\mathcal{E}\,$
and $\mathcal{F}\,$ in your hard copy but the typesetter will substitute
the more florid script characters normally seen in the journal.

\subsection{Skewing of accents}

Accents such as hats, overbars and dots are normally centred over letters,
but when these are italic or sloping greek the accent may need to be moved
to the right so that it is centred over the top of the sloped letter. For
example, \verb"\newcommand\hatp{\skew3\hat{p}}" will produce $\hatp$.

\subsection{Units of measure}

Numbers and their units of measure should be typeset with fixed spaces
that will not break over two lines. This is easily done with user-defined
macros. For example,\linebreak \verb"52\dynpercm" typesets as
52\dynpercm , providing the following macro definition has been included
in the preamble.
\begin{verbatim}
  \newcommand\dynpercm{\nobreak\mbox{$\;$dynes\,cm$^{-1}$}}
\end{verbatim}

\section{User-defined macros}\label{secUser}

If you define your own macros you must ensure that their names do not
conflict with any existing macros in plain \TeX\ or \LaTeX\ (or \AMSLaTeX\ %
if you are using this). You should also place them in the preamble to
your input file, between the \verb"\documentclass" and
\verb"\begin{document}" commands.

Apart from scanning the indexes of the relevant manuals, you can check
whether a macro name is already used in plain \TeX\ or \LaTeX\ by using
the \TeX\ command \verb"show". For instance, run \LaTeX\ interactively
and type \verb"\show\<macro_name>" at the \TeX\ prompt. (Alternatively,
insert the \verb"\show" command into the preamble of an input file and
\TeX\ it.)
\begin{verbatim}
  * \show\Re
\end{verbatim}
produces the following response from \TeX :
\begin{verbatim}
  > \Re=\mathchar"23C.
  <*> \show\Re
\end{verbatim}
By contrast \verb"\Real" is not part of plain \TeX\ or \LaTeX\ and
\verb"\show\Real" generates:
\begin{verbatim}
  > \Real=undefined
  <*> \show\Real
\end{verbatim}
confirming that this name can be assigned to a user-defined macro.

Alternatively, you can just use \verb"\newcommand", which will check for
the existence of the macro you are trying to define.  If the macro exists
\LaTeX\ will respond with:
\begin{verbatim}
  ! LaTeX Error: Command ... already defined.
\end{verbatim}
In this case you should choose another name, and try again.

Such macros must be in a place where they can easily be found and
modified by the journal's editors or typesetter. They must be gathered
together in the preamble of your input file, or in a separate
\verb"macros.tex" file with the command \verb"\input{macros}" in the
preamble. Macro definitions must not be scattered about your document
where they are likely to be completely overlooked by the typesetter.

The same applies to font definitions that are based on Computer Modern
fonts. These must be changed by the typesetter to use the journal's
\TeX\ typefaces. In this case, you should draw
attention to these font definitions on the hard copy that you submit for
publication and by placing a comment in your input file just before the
relevant definitions, for example \verb


\section{Some guidelines for using standard facilities}

The following notes may help you achieve the best effects with the
standard \LaTeX\ facilities that remain in the JPP style.

\subsection{Sections}

\LaTeX\ provides five levels of section headings, only four of which
are defined in the JPP class file:
\begin{itemize}
  \item[] Heading A -- \verb"\section{...}"
  \item[] Heading B -- \verb"\subsection{...}"
  \item[] Heading C -- \verb"\subsubsection{...}"
  \item[] Heading D -- \verb"\paragraph{...}"
\end{itemize}
There is no \verb"subparagraph" in the JPP style.

To obtain non-bold in a bold heading use the usual \LaTeXe\ commands for
changing typeface; for example \verb"\section{Fluctuations in Ca\textsc{ii}}".
  \begin{table}
    \caption{An example table} \label{sample-table}
    \begin{center}
      \begin{tabular}{cccc}
       \hline
         {Figure} & {$hA$} & {$hB$} & {$hC$}\\
       \hline
         2 & $\exp\;(\upi \mathrm{i} x)$
           & $\exp\;(\upi \mathrm{i} y)$ & $0$\\
         3 & $-1$    & $\exp\;(\upi \mathrm{i} x)$ & $1$\\
         4 & $-4+3\mathrm{i}$ & $-4+3\mathrm{i}$ & 1.6\\
         5 & $-2$    & $-2$    & $1.2 \mathrm{i}$
      \end{tabular}
    \end{center}
  \end{table}

\subsection{Tables}

The \verb"table" environment is implemented as described in
the \LaTeX\ manual to provide consecutively numbered floating
inserts for tables.

JPP class will cope with most table positioning problems
and you should not normally use the optional positional qualifiers \verb"t",
\verb"b", \verb"h" on the
\verb"table" environment, as this would override these decisions.
Table captions should appear at the bottom of the table; therefore you
should place the \verb"\caption" command before the
\verb"\begin{tabular}".

The JPP \verb"table" environment will insert rules above and below
the table as required by the JPP style, so you should not attempt
to insert these yourself.
The only time you need to intervene with the rules is when two or more tables
fall one above another. When this happens, you will get two rules together
and too much space between the tables; the solution is to add \verb"\norule"
before the end of upper table and \verb"\followon" between tables as shown here:
\begin{verbatim}
    ...
    \end{center}
    \caption{Dimensionless parameters}\label{dimensionless_p}
    \norule
  \end{table}
    \followon
  \begin{table}
    \begin{center}
    ...
\end{verbatim}

The JPP style dictates that vertical rules should never be used within the
body of the table, and horizontal rules should be used only to span
columns with the same headings. Extra space can be inserted to distinguish
groups of rows or columns.

As an example, table~\ref{sample-table} is produced using the following
commands:
\begin{verbatim}
  \begin{table}
    \caption{An example table} \label{sample-table}
    \begin{center}
      \begin{tabular}{cccc}
       \hline
         {Figure} & {$hA$} & {$hB$} & {$hC$}\\
       \hline
         2 & $\exp\;(\upi \mathrm{i} x)$
           & $\exp\;(\upi \mathrm{i} y)$ & $0$\\
         3 & $-1$    & $\exp\;(\upi \mathrm{i} x)$ & $1$\\
         4 & $-4+3\mathrm{i}$ & $-4+3\mathrm{i}$ & 1.6\\
         5 & $-2$    & $-2$    & $1.2 \mathrm{i}$
      \end{tabular}
    \end{center}
  \end{table}
\end{verbatim}

The \verb"tabular" environment
has been modified for the JPP style in the following ways:
\begin{enumerate}
  \item Additional vertical space is inserted above and below a horizontal
        rule produced by \verb"\hline"
  \item Tables are centred, and span the full width of the page; that is,
        they are similar to the tables that would be produced by 
		\verb"\begin{tabular*}{\textwidth}".
  \item As with normal \LaTeX, if a footnote needs to be inserted into the
        body of the table, the \verb"tabular" environment should be enclosed
        in a \verb"minipage" environment of width \verb"\textwidth". \eg
\begin{verbatim}
  \begin{minipage}{\textwidth}
    \begin{tabular}{cccc}
     \hline
     {Figure} & {$hA$} & {$hB$}\footnote{A table must be
      inside a \verb"minipage" environment if it includes
      table footnotes.}
       & {$hC$}\\
      ...
    \end{tabular}
  \end{minipage}
\end{verbatim}
\end{enumerate}
Commands to
redefine quantities such as \verb"\arraystretch" should be omitted. If
the original tabular facilities are needed, there is a new environment
`\verb"oldtabular"', which has none of the reformatting; it should be used
in exactly the same way.

\subsection{Illustrations (or figures)}

Artwork for all figures must be supplied separately, to be processed by
the Printer, and so \LaTeX's picture environment, for example, cannot
be used. An approximate amount of space should be left, using
the \verb"\vspace" command.

The JPP style will cope with most figure positioning problems and you should not
normally use the optional positional qualifiers \verb"t", \verb"b", \verb"h"
on the \verb"figure" environment, as this would override these decisions.
Figure captions should be below the figure itself, therefore the \verb"\caption"
command should appear after the space left for the illustration within the
\verb"figure" environment.
For example, figure~\ref{sample-figure} is produced using the following
commands:
\begin{figure}
  \vspace{3cm}
  \caption{An example figure with space for artwork.}
  \label{sample-figure}
\end{figure}

\begin{verbatim}
  \begin{figure}
    \vspace{3cm}
    \caption{An example figure with space for artwork.}
    \label{sample-figure}
  \end{figure}
\end{verbatim}

If a figure caption is too long to fit on the same page as its illustration,
the caption may be typeset as `{\small\textsc{Figure X.} \textrm{For caption
see facing page.}}', and the longer caption typeset at the bottom of the
facing page. Authors should not concern themselves unduly with such details,
and may leave pages long.

\subsection{Acknowledgments}

Acknowledgments should appear at the close of your paper, just before
the list of references and any appendices.  Use the \verb"acknowledgments"
environment, which simply inserts an `Acknowledgments' subsection heading.

There is also an \verb"acknowledgment" (or \verb"acknowledegment") environment
for typesetting a single acknowledegment.

\subsection{Appendices}

You should use the standard \LaTeX\ \verb"
\part{3}
\pagerange{119--126}
\date{1 July 1994 and in revised form 30 October 1998}
\setcounter{page}{119}

\begin{document}

\maketitle

\begin{abstract}
Using Stokes flow between eccentric counter-rotating cylinders as a prototype
for bounded nearly parallel lubrication flow, we investigate the effect of a
slender recirculation region within the flow field on cross-stream heat or mass
transport in the important limit of high P\'{e}clet number \Pen\ where the
enhancement over pure conduction heat transfer without recirculation is most
pronounced.  The steady enhancement is estimated with a matched asymptotic
expansion to resolve the diffusive boundary layers at the separatrices which
bound the recirculation region.  The enhancement over pure conduction is shown
to vary as $\epsilon^{\frac{1}{2}}$ at infinite \Pen, where
$\epsilon^{\frac{1}{2}}$
is the characteristic width of the recirculation region.
The enhancement decays from this asymptote as $\Pen\smh$.
\end{abstract}

\section{Introduction}

The use of integral equations to solve exterior problems in linear acoustics,
i.e.\ to solve the Helmholtz equation $(\nabla^2+k^2)\phi=0$ outside
a surface $S$ given that $\phi$ satisfies certain boundary conditions
on $S$, is very common. A good description is provided by \cite{Martin80}.
Integral equations have also been used to solve
the two-dimensional Helmholtz equation that arises in water-wave problems
where there is a constant depth variation.  The problem of wave
oscillations in arbitrarily shaped harbours using such techniques has
been examined (see for example Hwang \& Tuck 1970; Lee 1971).

In a recent paper \cite{Linton92} have shown how radiation
and scattering problems for vertical circular cylinders placed on
the centreline of a channel of finite water depth can be solved
efficiently using the multipole method devised originally by Ursell
(1950). This method was also used by \cite{Callan91}
to prove the existence of trapped modes in the vicinity of such a
cylinder at a discrete wavenumber $k<\upi/2d$ where $2d$ is the channel
width.


Many water-wave/body interaction problems in which the body is a
vertical cylinder with constant cross-section can be simplified by
factoring out the depth dependence. Thus if the boundary conditions
are homogeneous we can write the velocity potential
$\phi(x,y,z,t) =
 \hbox{Re}\{\phi(x,y)\cosh k(z+h){\rm e}^{-{\rm i}\omega t}\}$,
where the $(x,y)$-plane corresponds to the undisturbed free surface and
$z$ is measured vertically upwards with $z=-h$ the bottom of the channel.


Subsequently Callan {\it et al.} (1991) proved the existence of, and
computed the wavenumbers for, the circular cross-section case.
It should be noted however that experimental evidence for acoustic
resonances in the case of the circular cylinder is given by Bearman
\& Graham (1980, pp.~231--232).

\mbox{}\cite{Koch83} provided a theory for determining the trapped-mode frequencies
for the thin plate, based on a modification of the Wiener--Hopf
technique.

The use of channel Green's functions, developed in~\S\,\ref{sec:greensfunc}, allows
the far-field behaviour to be computed in an extremely simple manner,
whilst the integral equation constructed in~\S\,\ref{sec:trapmode} enables the trapped
modes to be computed in~\S\,\ref{sec:torus} and the scattering of an incident plane
wave to be solved in~\S\,\ref{sec:conc}. Appendix~\ref{app:boundcon} contains comparisons
with experiments.


\section{Green's functions}\label{sec:greensfunc}

\subsection{Construction of equations}

We are concerned with problems for which the solution, $\phi$, is
either symmetric or antisymmetric about the centreline of the
waveguide, $y=0$. The first step is the construction of a
symmetric and an antisymmetric Green's function, $\GsPQ$ and $\GaPQ$.
Thus we require
\begin{equation}
  (\nabla^2_P+k^2)G_s=(\nabla^2_P+k^2)G_a=0
  \label{Helm}
\end{equation}
in the fluid, where $\bnabla_\rho$ is a gradient operator, except when $P\equiv Q$,
\begin{equation}
  G_{\rm s},G_{\rm a}\sim 1 / (2\upi)\ln r
  \qquad \mbox{as\ } r\equiv|P-Q|\rightarrow 0,
  \label{singular}
\end{equation}
\begin{equation}
  \frac{\p G_{\rm s}}{\p y}=\frac{\p G_{\rm a}}{\p y}=0
  \qquad \mbox{on\ } y=d,
  \label{wallbc}
\end{equation}
\begin{equation}
  \frac{\p G_{\rm s}}{\p y}=0
  \qquad \mbox{on\ } y=0,
  \label{symbc}
\end{equation}
\begin{equation}
  G_{\rm a}=0
  \qquad \mbox{on\ } y=0,
  \label{asymbc}
\end{equation}
and we require $G_s$ and $G_a$ to behave like outgoing waves as
$|x|\tti$.

One way of constructing $G_s$ or $G_a$ is to replace (\ref{Helm})
and (\ref{singular}) by
\begin{equation}
  (\nabla^2_P+k^2)G_{\rm s}=(\nabla^2_P+k^2)G_{\rm a}
  = -\delta(x-\xi)\delta(y-\eta)
\end{equation}
and to assume initially that $k$ has a positive imaginary part.


\subsection{Further developments}

Using results from \cite{Linton92} we see that this has the integral
representation
\begin{equation}
  -\frac{1}{2\upi} \int_0^{\infty} \gamma^{-1}[{\rm e}^{-k\gamma|y-\eta|}
   + {\rm e}^{-k\gamma(2d-y-\eta)}] \cos k(x-\xi)t\:{\rm d}t,
   \qquad 0<y,\quad \eta<d,
\end{equation}
where
\[
  \gamma(t) = \left\{
    \begin{array}{ll}
      -{\rm i}(1-t^2)\sh, & t\le 1 \\[2pt]
      (t^2-1)\sh,         & t>1.
    \end{array} \right.
\]
In order to satisfy (\ref{symbc}) we add to this the function
\[
  -\frac{1}{2\upi}
   \pvi B(t)\frac{\cosh k\gamma(d-y)}{\gamma\sinh k\gamma d}
   \cos k(x-\xi)t\:{\rm d}t
\]
which satisfies (\ref{Helm}), (\ref{wallbc}) and obtain
\begin{equation}
  B(t) = 2{\rm e}^{-\kgd}\cosh k\gamma(d-\eta).
\end{equation}
Thus the function
\begin{equation}
  G = -\squart{\rm i}(H_0(kr)+H_0(kr_1))
    - \frac{1}{\upi} \pvi\frac{{\rm e}^{-\kgd}}%
    {\gamma\sinh\kgd} \cosh k\gamma(d-y) \cosh k\gamma(d-\eta)
\end{equation}
satisfies (\ref{Helm})--(\ref{symbc}). By writing this function
as a single integral which is even in $\gamma$, it follows that
$G$ is real.



\section{The trapped-mode problem}\label{sec:trapmode}

The unit normal from $D$ to $\p D$ is $\bmath{n}_q=(-y^{\prime}(\theta),
x^{\prime}(\theta))/w(\theta)$.
Now $G_a=\squart Y_0(kr)+\Gat$ where
$r=\{[x(\theta)-x(\psi)]^2 + [y(\theta)-y(\psi)]^2\}\sh$ and $\Gat$ is
regular as $kr\ttz$. In order to evaluate
$\p G_a(\theta,\theta)/\p n_q$ we note that
\[
  \ndq\left(\squart Y_0(kr)\right)\sim\ndq\left(\frac{1}{2\upi}
  \ln r\right)=\frac{1}{2\upi r^2 w(\theta)}[x^{\prime}(\theta)
  (y(\theta)-y(\psi))-y^{\prime}(\theta)(x(\theta)-x(\psi))]
\]
as $kr\ttz$.  Expanding $x(\psi)$ and $y(\psi)$ about the point
$\psi=\theta$ then shows that
\begin{eqnarray}
  \ndq\left(\squart Y_0(kr)\right) & \sim &
    \frac{1}{4\upi w^3(\theta)}
    [x^{\prime\prime}(\theta)y^{\prime}(\theta)-
    y^{\prime\prime}(\theta)x^{\prime}(\theta)] \nonumber\\
  & = & \frac{1}{4\upi w^3(\theta)}
    [\rho^{\prime}(\theta)\rho^{\prime\prime}(\theta)
    - \rho^2(\theta)-2\rho^{\prime 2}(\theta)]
    \quad \mbox{as\ } kr\ttz . \label{inteqpt}
\end{eqnarray}


\subsection{Computation}

For computational purposes we discretize (\ref{inteqpt}) by dividing
the interval ($0,\upi$) into $M$ segments. Thus we write
\begin{equation}
  \shalf\phi(\psi) = \frac{\upi}{M} \sumjm\phi(\theta_j)\ndq
  G_a(\psi,\theta_j)\,w(\theta_j),
  \qquad 0<\psi<\upi,
  \label{phipsi}
\end{equation}
where $\theta_j=(j-\shalf)\upi/M$. Collocating at $\psi=\theta_i$ and
writing $\phi_i=\phi(\theta_i)$ etc.\ gives
\begin{equation}
  \shalf\phi_i = \frac{\upi}{M} \sumjm\phi_j K_{ij}^a w_j,
  \qquad i=1,\,\ldots,\,M,
\end{equation}
where
\begin{equation}
  K_{ij}^a = \left\{
    \begin{array}{ll}
      \p G_a(\theta_i,\theta_j)/\p n_q, & i\neq j \\[2pt]
      \p\Gat(\theta_i,\theta_i)/\p n_q
      + [\rho_i^{\prime}\rho_i^{\prime\prime}-\rho_i^2-2\rho_i^{\prime 2}]
      / 4\upi w_i^3, & i=j.
  \end{array} \right.
\end{equation}
For a trapped mode, therefore, we require the determinant of the
$M\times M$ matrix whose elements are
\[
  \delta_{ij}-\frac{2\upi}{M}K_{ij}^a w_j,
\]
to be zero.
\begin{table}
  \caption{Values of $kd$ at which trapped modes occur when $\rho(\theta)=a$}
  \label{tab:kd}
  \begin{center}
  \begin{tabular}{cccc}
    \hline
      $a/d$ & $M=4$  & $M=8$   & Callan {\it et al.}\\
    \hline
       0.1 & 1.56905 & 1.56905 & 1.56904\\
       0.3 & 1.50484 & 1.50484 & 1.50484\\
       0.5 & 1.39128 & 1.39131 & 1.39131\\
       0.7 & 1.32281 & 1.32287 & 1.32288\\
       0.9 & 1.34479 & 1.35160 & 1.35185
  \end{tabular}
  \end{center}
\end{table}
Table~\ref{tab:kd} shows a comparison of results
obtained from this method using two different truncation parameters
with accurate values obtained using the method of Callan {\it et al.} (1991).

\begin{figure}
  \vspace{16.5pc}
  \caption{Trapped-mode wavenumbers, $kd$, plotted against $a/d$ for
    three ellipses:\protect\\
    \hbox{{--}\,{--}\,{--}\,{--}\,{--}}, $b/a=0.5$; ---$\!$---,
    $b/a=1$; ---\,$\cdot$\,---, $b/a=1.5$.}\label{fig:trapmode}
\end{figure}
An example of the results that are obtained from our method is given
in figure~\ref{fig:trapmode}.
Figure~\ref{fig:shadcon}\,({\it a},{\it b}) shows shaded contour plots of
$\phi$ for these modes, normalized so that the maximum value of $\phi$ on
the body is 1.
\begin{figure}
  \vspace{15pc}
  \caption{Shaded contour plots of the potential $\phi$ for the two trapped
    modes that exist for an ellipse with $a/d=1.5$, $b/d=0.75$.
    ({\it a}) Symmetric about $x=0$, $kd=0.96$;
    ({\it b}) antisymmetric about $x=0$, $kd=1.398$.}\label{fig:shadcon}
\end{figure}


\subsection{Basic properties}

Let
\def\thesection.\arabic{equation}a,b{\thesection.\arabic{equation}a,b}
\begin{equation}
\rho_l = \lim_{\zeta \rightarrow Z^-_l(x)} \rho(x,\zeta), \quad
  \rho_{\rm u} = \lim_{\zeta \rightarrow Z^{+}_{\rm u}(x)} \rho(x,\zeta)
\end{equation}%
\def\thesection.\arabic{equation}a,b{\thesection.\arabic{equation}}%
be the fluid densities immediately below and above the cat's-eyes.  Finally
let $\rho_0$ and $N_0$ be the constant values of the density and the vorticity
inside the cat's-eyes, so that
\begin{equation}
  (\rho(x,\zeta),\phi_{\zeta\zeta}(x,\zeta))=(\rho_0,N_0)
  \quad \mbox{for\ } Z_l(x) < \zeta < Z_{\rm u}(x).
\end{equation}

The Reynolds number \Rey{\it \/} is defined by $u_\tau H/\nu$ ($\nu$ is the kinematic
viscosity), the length given in wall units is denoted by $(\;)_+$, and the
Prandtl number \Pran{\it \/} is set equal to 0.7. In (\ref{Helm}) and (\ref{singular}),
$\tau_{ij}$ and $\tau^\theta_j$ are
\begin{subeqnarray}
  \tau_{ij} & = &
    (\overline{\overline{u}_i \overline{u}_j}
    - \overline{u}_i\overline{u}_j)
    + (\overline{\overline{u}_iu^{SGS}_j
    + u^{SGS}_i\overline{u}_j})
    + \overline{u^{SGS}_iu^{SGS}_j},\\[3pt]
  \tau^\theta_j & = &
    (\overline{\overline{u}_j\overline{\theta}}
    - \overline{u}_j \overline{\theta})
    + (\overline{\overline{u}_j\theta^{SGS}
    + u^{SGS}_j \overline{\theta}})
    + \overline{u^{SGS}_j\theta^{SGS}}.
\end{subeqnarray}

\subsubsection{Calculation of the terms}

The first terms in the right-hand side of (3.5a) and
(3.5b) are the Leonard terms explicitly calculated by
applying the Gaussian filter in the $x$- and $z$-directions in the
Fourier space.

The interface boundary conditions given by (\ref{Helm}) and (\ref{singular}),
which relate the displacement and stress state of the wall at the mean
interface to the disturbance quantities of the flow, can also be reformulated
in terms of the transformed quantities.  The transformed boundary conditions
are summarized below in a matrix form that is convenient for the subsequent
development of the theory:
\begin{equation}
\setlength{\arraycolsep}{0pt}
\renewcommand{\arraystretch}{1.3}
\slsQ_C = \left[
\begin{array}{ccccc}
  -\omega^{-2}V'_w  &  -(\alpha^t\omega)^{-1}  &  0  &  0  &  0  \\
  \displaystyle
  \frac{\beta}{\alpha\omega^2}V'_w  &  0  &  0  &  0  &  {\rm i}\omega^{-1} \\
  {\rm i}\omega^{-1}  &  0  &  0  &  0  &  0  \\
  \displaystyle
  {\rm i}R^{-1}_{\delta}(\alpha^t+\omega^{-1}V''_w)  &  0
    & -({\rm i}\alpha^tR_\delta)^{-1}  &  0  &  0  \\
  \displaystyle
  \frac{{\rm i}\beta}{\alpha\omega}R^{-1}_\delta V''_w  &  0  &  0
    &  0  & 0 \\
  ({\rm i}\alpha^t)^{-1}V'_w  &  (3R^{-1}_{\delta}+c^t({\rm i}\alpha^t)^{-1})
    &  0  &  -(\alpha^t)^{-2}R^{-1}_{\delta}  &  0  \\
\end{array}  \right] .
\label{defQc}
\end{equation}
$\biS^t$ is termed the displacement-stress vector and $\slsQ_C$
the flow-wall coupling matrix.  Subscript $w$ in (\ref{defQc}) denotes
evaluation of the terms at the mean interface.  It is noted that $V''_w = 0$
for the Blasius mean flow.


\subsubsection{Wave propagation in anisotropic compliant layers}

From (\ref{Helm}), the fundamental wave solutions to (\ref{inteqpt}) and
(\ref{phipsi}) for a uniformly thick homogeneous layer in the transformed
variables has the form of
\begin{equation}
\etb^t = \hat{\etb}^t \exp [{\rm i}(\alpha^tx^t_1-\omega t)],
\end{equation}
where $\hat{\etb}^t=\bmath{b}\exp ({\rm i}\gamma x^t_3)$.  For a
non-trivial wave, the substitution of (1.7) into (1.3) and (1.4) yields
the following determinantal equation:
\begin{equation}
\mbox{Det}[\rho\omega^2\delta_{ps}-C^t_{pqrs}k^t_qk^t_r]=0,
\end{equation}
where the wavenumbers $\langle k^t_1,k^t_2,k^t_3\rangle = \langle
\alpha^t,0,\gamma\rangle $ and $\delta_{ps}$ is the Kronecker delta.


\section{Torus translating along an axis of symmetry}\label{sec:torus}

Consider a torus with axes $a,b$ (see figure~\ref{fig:shadcon}), moving along
the $z$-axis. Symmetry considerations imply the following form for the stress
function, given in body coordinates:
\begin{equation}
\bmath{f}(\theta,\psi) = (g(\psi) \cos \theta,g(\psi) \sin \theta,f(\psi)).
\label{eq41}
\end{equation}

Because of symmetry, one can integrate analytically in the $\theta$-direction
obtaining a pair of equations for the coefficients $f$, $g$ in (4.1),
\begin{eqnarray}
f(\psi_1) = \frac{3b}{\upi(2(a+b \cos \psi_1))^{\frac{3}{2}}}
  \int^{2\upi}_0 \frac{(\sin \psi_1 - \sin \psi)(a+b \cos \psi)\sh}%
  {(1 - \cos (\psi_1 - \psi))(2+\alpha)\sh},
\label{eq42}
\end{eqnarray}
\begin{eqnarray}
g(\psi_1) & = & \frac{3}{\upi(2(a+b \cos \psi_1))^{\frac{3}{2}}}
  \int^{2\upi}_0 \left(\frac{a+b \cos \psi}{2+\alpha}\right)\sh
  \left\{ \usebox{\astrutbox} f(\psi)[(\cos \psi_1 - \cos \psi - b \beta_1)S + \beta_1P]
  \right. \nonumber\\
&& \mbox{}\times \frac{\sin \psi_1 - \sin \psi}{1-\cos(\psi_1 - \psi)}
  + g(\psi) \left[\left(2+\alpha - \frac{(\sin \psi_1 - \sin \psi)^2}
  {1- \cos (\psi - \psi_1)} - b^2 \gamma \right) S \right.\nonumber\\
&& \left.\left.\mbox{} + \left( b^2 \cos \psi_1\gamma -
  \frac{a}{b}\alpha \right) F(\thalf\upi, \delta) - (2+\alpha)
  \cos\psi_1 E(\thalf\upi, \delta)\right] \usebox{\astrutbox}\right\} {\rm d}\psi,
\label{eq43}
\end{eqnarray}
\begin{equation}
\alpha = \alpha(\psi,\psi_1) = \frac{b^2(1-\cos(\psi-\psi_1))}%
  {(a+b\cos\psi) (a+b\cos\psi_1)};
  \quad
  \beta - \beta(\psi,\psi_1) = \frac{1-\cos(\psi-\psi_1)}{a+b\cos\psi}.
\end{equation}


\section{Conclusions}\label{sec:conc}

We have shown how integral equations can be used to solve a particular class
of problems concerning obstacles in waveguides, namely the Neumann problem
for bodies symmetric about the centreline of a channel, and two such problems
were considered in detail.


\begin{acknowledgments}
We would like to acknowledge the useful comments of a referee concerning
the solution procedure used in \S\,\ref{sec:conc}.  A.\,N.\,O. is supported by SERC
under grant number GR/F/12345.
\end{acknowledgments}



\part{0}
\pagerange{\pageref{firstpage}--\pageref{lastpage}}
\doi{S0963548301004989}

\begin{document}

\label{firstpage}
\maketitle

\begin{abstract}
The usual theory of plasma relaxation, based on the selective
decay of magnetic energy over the (global) magnetic helicity,
predicts a force-free state for a plasma. Such a force-free state
is  inadequate to describe most  realistic plasma systems
occurring in laboratory and  space plasmas as it produces a zero
pressure gradient and  cannot couple magnetic fields with flow. A
different theory of relaxation has been proposed by many authors,
based on a well-known principle of irreversible thermodynamics,
the principle of minimum entropy production rate which is
equivalent to the minimum dissipation rate (MDR) of energy. We
demonstrate the applicability of minimum dissipative relaxed
states to various self-organized systems of magnetically confined
plasma in the laboratory and in the astrophysical context. Such
relaxed states are shown to produce a number of basic
characteristics of laboratory plasma confinement systems and solar
arcade structure.
\end{abstract}





\section{Introduction and Motivation:}

Self-organization is a natural process [Ortolani and Schnack,
1993] in which a continuous system evolves toward some preferred
states showing a form of order on long scales. Examples of such
self-organization processes are ubiquitous in nature and in
systems studied in the laboratory. These ordered states are very
often remarkably robust, and their detailed structures are mostly
independent of the way the system is prepared.  The final state to
which a system evolves generally depends on the boundary
conditions, inherent geometry of the particular device, but is
relatively independent of the initial conditions.
Self-organization in plasma is generally called Plasma Relaxation.
The seminal theory  for plasma relaxation first proposed by Taylor
(1974) yields a force-free state with zero pressure gradient.
Taylor's theory was the first to explain the experimentally
observed field reversal of the Reversed Field Pinch (RFP).

Self-organized states with zero pressure gradient are often far
from experimentally realizable scenarios. Any realistic magnetic
configuration confining plasma must have a non-zero pressure
gradient. Extensive numerical simulations [ Sato et al., 1996, Zhu
et al., 1995] have established the existence of self-organized
states of plasma with finite pressure, i.e., these states  are
observed to be governed by the magnetohydrodynamic  force balance
equation. Numerical simulations [Dasgupta et al, 1995, Watanabe et
al, 1997] and experiments [Ono et al., 1993 ] have established
that counterhelicity merging of two spheromaks can produce a
Field-Reversed Configuration (FRC). The  FRC has a zero toroidal
magnetic field, and the plasma is confined entirely by poloidal
magnetic field. It has a finite pressure,  nonzero perpendicular
component of current  and is more often characterized by a high
value of plasma beta. The stability and long life categorize FRC
as a relaxed state, which is not obtainable from Taylor model of
relaxation. Since spheromaks  are depicted as a Taylor state
[Rosenbluth and Bussac, 1979],  formation of a FRC from the
merging of two counterhelicity spheromaks is a unique process from
the point of view of plasma relaxation,  where  a non-force free
state (hence non-Taylor state)  emerges from the coalescence of
two Taylor states. In  astrophysical context, particularly solar
physics scenario, a 3-D MHD simulation by Amari and Luciano
[2000], among others, showed that, after the initial helicity
drive, the final ``relaxed state is far from the ... linear
force-free model that would be predicted by Taylor's conjecture"
and they suggested to derive an alternative variational approach.
Such  physical processes call for an alternate model for plasma
relaxation to broadly accommodate a large number of  these
observations.

In this context, an important work by Turner deserves special
attention. Turner first showed (1986) that a relaxed magnetic field
configuration which can support/confine a plasma with a finite
pressure gradient. Adopting a two-fluid approach, a model of
magnetofluid relaxation is constructed for Hall MHD, under the
assumptions of minimum energy (magnetic plus fluid) with the
constraints of global magnetic helicity hybrid helicity, axial
magnetic flux and fluid vorticity. Euler-Lagrange equations resulted
from such variational approach show the coupling between the magnetic
field and the fluid vorticity.  The solutions for such coupled
equations are shown to be achieved by the linear superposition of the
eigen-vectors of the 'curl' operator (force-free states). These
solutions yield magnetic field configurations which can confine plasma
with a finite pressure gradient.


The principle of  ``minimum rate of entropy production" formulated
by Prigogine and others [Prigogine, 1947]  is believed to play a
major role in many problems of irreversible thermodynamics. This
principle states that for any steady irreversible processes, the
rate of entropy production is minimum. For most irreversible
processes  in nature, the minimum rate of entropy production is
equivalent to the minimum rate of energy dissipation.  A
magnetized plasma is such a dissipative system and it is
appropriate to expect that  the principle of minimum dissipation
rate of energy has a major role to play in a magnetically confined
plasma.   It is worth mentioning that Rayleigh [1873] first used
the term ``principle of least dissipation of energy" in his works
on the propagation of elastic waves in matter.  Chandasekhar and
Woltjer [1958]  considered a relaxed state of plasma with minimum
Ohmic dissipation with  the constraint of constant magnetic energy
and obtained an equation for the magnetic field involving a
`double curl' for the magnetic field and remarked that the
solution of such equation is ``much wider" than the usual force
free solution.


Montgomery and Phillips  [1988] first used the principle of
minimum dissipation rate  (MDR)  of energy in an MHD problem to
understand the steady the steady profile of a RFP under the
constraint of a constant rate of supply and dissipation of
helicity with the boundary conditions for a conducting wall.
Farengo et al [Farengo et al, 1994,1995, 2002, Bouzat, 2006] in a
series of papers applied  the principle of MDR to a variety of
problems ranging from flux-core spheromak to tokamak.

In this work, we present an alternative scenario for plasma
relaxation which is based on the principle of minimum dissipation
rate (MDR) rate of energy. This model gives us a non-force free
magnetic field, capable of supporting finite pressure gradient. We
demonstrate the this model can reproduce some of the basic
characteristics of most of the plasma confining devices in the
laboratory  as well as the arcade structure of solar magnetic
field. A two-fluid generalization of this model for an open system
couples magnetic field with flow - so this model is applicable to
other astrophysical situations.

The plan of the paper is as follows: In section 2, we present a
single fluid MDR model for closed system and briefly discuss its
applicarions to some of the laboratory plasma confinement devices,
such as RFP, Spheromak, Tokamak and FRC. We present the result of
a recent numerical simulation to justify our choice of the energy
dissipation rate as an effective minimizer in our variation
problem. In section 3 we describe a generalized version of the MDR
model for open two fluid system and its application to solar
arcade problem. We summarize a numerical extrapolation method for
non-force free coronal magnetic field based on our MDR model.
Section 4 concludes our paper.
\section{Single fluid MDR based model}


To obtain both the field reversal and a finite pressure gradient
for RFP, Dasgupta et al.[ Dasgupta et al., 1998] considered the
relaxation of a slightly resistive and turbulent plasma using the
principle of MDR under the constraint of constant global helicity.
With the global helicity $K_M$ and the dissipation rate $R$
defined as,
\begin{equation}
K_M = \int_V {\bf A}\cdot{\bf B}dV; \qquad R = \int_V \eta {\bf
j}^2 dV, \end{equation}

the variational principle $\delta (R-\Lambda K)$ leads to the
following Euler-Lagrange equation,

\begin{equation}\nabla\times\nabla\times\nabla\times{\bf B}=
\Lambda{\bf B}.\end{equation}

 Solutions of the above equation can be constructed as a
 superposition of the force free equation using
 Chandrasekhar-Kendall (CK)  eigenfunctions [Chandraskhar and
 Kendall, 1957]

\begin{equation}B =  \alpha_1 B_1 + \alpha_2B_2  + \alpha_3B_3,
\end{equation}
where $\alpha_1, \alpha_2, \alpha_3 $ are constants, to be
determined from the boundary conditions, and ${\bf B}_i's$ are
obtained from \begin{equation} \nabla\times{\bf B}_i = \lambda_i
B_i;\quad    (i = 1,2,3);\quad \lambda_i^3 =
\Lambda.\end{equation} We have used the analytic continuation of
CK eigenfunctions in the complex  domain.  However, one can easily
ascertain that the resulting magnetic field B is real. Moreover,
since the superposition of  force free fields with different
eigenvalues $\lambda_i$ is not a force free field, the resultant
magnetic field B is not force free, ${\bf J}\times {\bf B}\neq 0
$.  An explicit form of the solution in cylindrical geometry with
the boundary conditions for an RFP (perfectly conducting wall) has
shown two remarkable features: (1) Such a state can support a
pressure gradient; and (2) Field reversal is found in states that
are not force free.

\begin{figure}[ht]
\epsfig{file=RFP-one.eps, width=10.cm, height=8.cm}
\epsfig{file=RFP-two.eps, width=10.cm, height=8.cm}
\caption{\small\label{fig:RFPI} Left panel: Plot of $F =
B_z(a)/\left<B_z\right> $ and $\Theta =
B_\theta(a)/\left<B_z\right> $ ($<..> $ is the volume average and
$ a$ is the radial distance) showing the field reversal of RFP
obtained from MDR model. The dotted curve represents the plot for
the minimum energy state of Taylor. Right Panel: Pressure profile
of RFP obtained from MDR model.(\textit{from Dasgupta et al.,
1998}).}
\end{figure}

 In Figure \ref{fig:RFPI} a  plot of the field reversal
parameter $F$ against the pinch parameter $\Theta$ is shown. It is
found that axial field reverses at a value $\Theta \simeq 2.2$;
and the pressure profile of the RFP obtained from our MDR model.


A success of this theory is its applicability in a tokamak
configuration [Bhattacharya et al., 2000]. In toroidal geometry,
the CK eigenfunctions are hypergeometric functions in the large
aspect ratio approximation. Equilibria can be constructed under
the assumption that the total current J = 0 at the edge.  The
solution reproduces the q-profile, the toroidal magnetic field,
the pressure profile and $\sigma = {\bf J}\cdot {\bf B} /B^2$ for
tokamak, which are very close to the observed profiles (e.g.,
figure \ref{fig:toka-one}). The solutions allow for a tokamak,
low-q discharge as well as RFP like behavior with a change of
eigenvalue which would essentially determine the amount of
volt-sec associated with the discharge characterizing the relaxed
state. The poloidal Volt-sec/toroidal flux associated with the
discharge can be directly interpreted in terms of the parameter
$\lambda r_0$ ($r_0$ is the minor radius of the torus, $\lambda
=\sqrt[3]{\Lambda}$) and can be shown to increase with $\lambda
r_0$. MDR relaxed state can yield several kinds of solutions with
distinct q-profiles. For low values of $\lambda r_0$, the
q-profile is monotonically increasing and hence the solutions
resemble a tokamak type discarge together with a nearly constant
toroidal field, for larger values of $\lambda r_0$ the toroidal
field reverses at the edge, signifying an RFP-type behavior and
for intermediate values, the relaxed states indicate an ultra
low-q type of discharge exhibiting a q-profile with a pitch
minimum and $0<q<1 $. In figure \ref{fig:toka-one}, we plot the
toroidal magnetic field $B_\phi, \sigma = {\bf J}\cdot{\bf J}/B^2,
q$ and pressure profiles for different values of $\lambda r_0$,
for a fixed value of the aspect ratio = 4.0. The continuous line
corresponds corresponds to $\lambda r_0 = 1.5$, and shows a
monotonic increase with $r$, which is the essential feature of
tokamak discharge. For this case the the toroidal magnetic field
($B_\phi$) is essentially constant and the pressure profile has a
peak at the center of the minor cross section of the torus and
falls to zero at the edge. A non-constant behavior is shown by the
$\sigma$ profile. The dotted line corresponds to $\lambda r_0=
3.0$. In this case the toroidal field reverses sign at the edge,
which corresponds to an RFP like discharge. The q-profile shows a
similar behavior. The $\sigma$ profile has a dip at $r=0$ and
exhibits a bump near the outer region. Such profiles have been
experimentally observed in the ETA-BETA-II device [ Antoni et al,
1983]. For the values $\lambda r_0= 2.5$ the q-profile shows a
nonmonotonic behavior as shown in figure \ref{fig:toka-two}. The
toroidal magnetic field again peaks at the center and slowly
decreases towards the edge, showing a paramagnetic behavior that
is typical of ultra low-q discharge [Yoshida et al., 1986].
Pressure profile peaks at the center and goes to zero at the edge.
Some of these results have been observed in the SINP tokamak
[Lahiri et al, 1996] and were also reported earlier in REPUTE-1
[Yoshida et al., 1986] and in TORIUT-6 [Yamada et al, 1987].


\begin{figure}[ht]
\epsfig{file=Tokamak-1.eps, width=10.cm, height=8.cm}
\caption{\small\label{fig:toka-one} q-profile, $B_\phi, \sigma$,
and pressure profile of (i)Tokamak discharge (continuous line);
(ii) RFP discharge (dotted line); (\textit{from Bhattacharya et
al., 2000}).}
\end{figure}


This model has been utilized to obtain the spheromak configuration
from the spherically symmetric solution under appropriate boundary
conditions relevant to an insulating boundary at the edge
[Dasgupta et al., 2002]. These boundary conditions are are
consistent with the classical spheromak solution used by
Rosenbluth and Bussac [1979]. The most interesting part of this
spheromak solution  is the non-constant profile for ${\bf
J}\cdot{\bf B}/B^2$  plotted against normalized distance from the
magnetic axis.  This profile shows peaks outside the magnetic
axis, and this feature is qualitatively closer to the
experimentally observed profile reported by Hart et al. [Hart et
al., 1983].



\begin{figure}[ht!]
\epsfig{file=Tokamak-II.eps, width=10.cm, height=8.cm}
\caption{\small\label{fig:toka-two} q-profile, $B_\phi, \sigma$,
and pressure profile of  low-q Discharge;  (\textit{from
Bhattacharya et al., 2000}).}
\end{figure}

The field-reversed configuration (FRC)  is a compact toroidal
confinement device [Tuszewski, 1988]  which is formed with
primarily poloidal magnetic field and  a zero or negligibly small
toroidal magnetic field. The plasma beta of the FRC is one of the
highest for any magnetically confined plasma.  FRC can be
considered as a relaxed state [Guo et al, 2005, 2006], but it is a
distinctly non-Taylor state. The FRC configuration can be obtained
for an externally driven dissipative system  which relaxes  with
minimum dissipation.  In Bhattacharya et al., [2001] we show that
the Euler-Lagrange equations for relaxation under MDR with
magnetic energy as the constraint represents an FRC configuration
for one choice of eigenvalues. This solution is characterized by
high beta, zero helicity and a completely null toroidal field and
supports  a non-constant pressure profile (non-force free).
Another choice of eigenvalue leads to a spheromak configuration.
A generalization of this work incorporating flow [Bhattacharya et
al., 2003] and taking magnetic and  flow energies and cross
helicities as constraints also produces the FRC configuration.
This state supports field aligned flows with strong shears  that
may lead to stability and better confinement.





\begin{figure}[ht!]
\epsfig{file=Spheromak.eps, width=10.cm, height=8.cm}
\epsfig{file=Spheroeps-I.eps, width=10.cm, height=8.cm}
\caption{\small\label{fig:sphero} Spheromak configuration (left
panel) and the non-constant $J_{||}/B$ derived from our model
based on MDR (right panel); (\textit{from Dasgupta et al.,
2002}).}
\end{figure}

\subsection{Simulation Results:} To justify the  use of the MDR
principle in plasma relaxation and the choice of the total
dissipation rate as a minimizer in the variational principle we
performed a fully 3-D numerical simulations for a turbulent and
slightly dissipative plasma  using a fully compressible MHD code
[Shaikh et al., 2007] with periodic boundary conditions.
Fluctuations are Fourier expanded and nonlinear terms are
evaluated in real space. Spectral methods provide accurate
representation in Fourier space. In our code there is very little
numerical dissipation, and hence the ideal invariants are
preserved. The selective decay processes (in addition to
dissipation) depend critically on the cascade properties
associated with the rugged MHD invariants that eventually govern
the spectral transfer in the inertial range. This can be
elucidated as follows [Biscamp, 2003]. Magnetic vector potential
in 3D MHD dominates, over the magnetic field fluctuations, at the
smaller Fourier modes, which in turn leads to a domination of the
magnetic helicity invariant over the magnetic energy. On the other
hand, dissipation occurs predominantly at the higher Fourier modes
which give rise to a rapid damping of the energy dissipative
quantity $R$. A heuristic argument for this process can be
formulated in the following way. The decay rates of helicity $K$
and dissipation rate $R= \int_V\eta j^2 dV$ in the dimensionless
form, with the magnetic field Fourier decomposed as ${\bf B}({\bf
k},t)= \Sigma_k {\bf b_k}\exp(i{\bf k}\cdot{\bf r}) $, can be
expressed as,


\begin{equation}\label{comp1}
\frac{d K}{dt}= -\frac{2\eta}{S}\Sigma_k k {\bf b_k}^2;\qquad
\frac{d R}{dt}= -\frac{2\eta^2}{S^2}\Sigma_k k^4 {\bf b_k}^2
\end{equation}
where $S= \tau_R/\tau_A$, is the Lundquist number, $\tau_R,
\tau_A$ are the resistive and Alfv\'{e}n time scales,
respectively. The Lundquist number in our simulations varies
between $10^6$ and $10^7$. We find that at scale lengths for which
$k\approx S^{1/2}$, the decay rate of energy dissipation is $\sim
O(1)$ . But at these scale lengths, helicity dissipation is only
$\sim O(S^{-1/2})\ll 1$. This physical scenario is consistent with
our 3D simulations. Figure \ref{fig:sim-three} shows the time
evolution of the global helicity $K_M$, total magnetic energy
$W_M$ and the total dissipation rate $R$. It is seen that global
helicity remains approximately constant, (decaying very slightly
during the simulation time) while the magnetic energy decays
faster than the global helicity, but the energy dissipation rate
decays even faster than the magnetic energy.  The results of the
simulation show that total dissipation rate can serve as a
minimizer during a relaxation process.




\begin{figure}[ht]
\begin{center}
\epsfig{file=Sim-two.eps, width=10.cm, height=8.cm}
\end{center}
\caption{\small\label{fig:sim-three} Time evolution of decay rates
associated with the turbulent relaxation of the ideal MHD
invariants, viz, global magnetic helicity $K_M= \int_V {\bf
A}\cdot{\bf B}dV$, magnetic energy $W= \int_V B^2/2dV$ and  the
total dissipation $R=\int_V \eta J^2 dV$, all shown in the same
scale. }
\end{figure}


\section{MDR based relaxation model for two fluid plasma with
external drive:} Recently, Bhattacharya and Janaki [2004]
developed a theory of dissipative relaxed states in two-fluid
plasma with external drive.  We describe the plasma (consisting of
ions and electrons) by the two-fluid equations. The basic
advantages of a two-fluid formalism over the single-fluid or MHD
counterpart are: (i) in a two-fluid formalism, one can incorporate
the coupling of plasma flow with the magnetic field in a natural
way, and (ii) the two-fluid description being more general than
MHD, is expected to capture certain results that are otherwise
unattainable from a single-fluid description. It should be
mentioned that although we are invoking the principle of minimum
dissipation rate of energy - instead of the principle of minimum
energy - our approach is closely similar to the approach to plasma
relaxation including Hall terms, as developed by Turner (1986).


In the two-fluid formalism the canonical momentum ${\bf{P}}_j$
canonical vorticity ${\bf{\Omega}}_j$, for the $j$-th species are
defined respectively as,

\begin{equation}
{\bf{P}}_j=m_j{\bf{u}}_j+\frac{q_j}{c}{\bf{A}}, \qquad
{\bf{\Omega}}_j=\nabla\times{\bf{P}}_j, \quad (j= i, e).
\end{equation}
A gauge invariant expression for the generalized helicity is
introduced by defining
\begin{equation}
K_j=\int_V{\bf{P}}_j\cdot{\bf{\Omega}}_j d\tau-
\int_V{\bf{P}}_j^\prime\cdot{\bf{\Omega}}_j^\prime dV;
\end{equation}
which is a relative helicity with respect to reference field
having ${\bf{P}}_j^\prime$  as the canonical momentum, and
${\bf{P}}_\alpha$ and ${\bf{P}}_\alpha^\prime$ are different
inside the volume of interest and the same outside. The gauge
invariance condition holds
\begin{equation}\label{omega}
({\bf{\Omega}}_j-{\bf{\Omega}}_j^\prime) \cdot\hat{n} =0
\Rightarrow ({\bf{B}}-{\bf{B}}^\prime)\cdot\hat{n}=0 \qquad
({\bf{\omega}}_j-{\bf{\omega}}_j^\prime) \cdot\hat{n}=0, \quad
({\bf{\omega}}_j=\nabla\times{\bf{u}}_j),
\end{equation}
under the condition that there is no flow-field coupling at the
surface. The above equations  then represent natural boundary
conditions inherent to the problem.



The helicity balance equation for the generalized helicity,
obtained from the balance equations for the generalized momentum
and vorticity, consists of two terms, the helicity injection terms
and helicity dissipation term (containing resistivity and
viscosity). The global helicity can be taken as a bounded function
following Montgomery  et al.,(1988) and it is possible to form a
variational problem  with the helicity dissipation rates as
constraints. Thus,
\begin{eqnarray}
\frac{d K_\alpha}{dt}&=&-\oint[G_\alpha-{\bf{P}}_\alpha\cdot{\bf
{u}}_\alpha]{\bf{\Omega}}_\alpha^\prime\cdot\hat{n}da-\oint({\bf{P}}_\alpha
\cdot{\bf{\Omega}}_\alpha){\bf{u}}_\alpha\cdot\hat{n}da \nonumber\\
&-&\oint\left(\frac{\partial{\bf{P}}_\alpha^\prime}{\partial t}
 \times{\bf{P}}_\alpha^\prime\right)\cdot\hat{n}da
 -\eta_\alpha^\prime\oint[\{\nabla\times{\bf{B}}+L_\alpha
 {\bf{\omega}}_\alpha)\}\times{\bf{P}}_\alpha]\cdot\hat{n}da \nonumber\\
 &-&2\eta_\alpha^\prime\int[\nabla\times({\bf{B}}+L_\alpha
 {\bf{\omega}}_\alpha)]\cdot{\bf{\Omega}}_\alpha d\tau
 \end{eqnarray}with
\begin{eqnarray}
 & &G_\alpha=\left[q_\alpha\phi+\frac{h_\alpha}{n_\alpha}-\frac{m_\alpha
 u_\alpha^2}{2}\right]; \qquad
\eta_\alpha^\prime=\frac{q_\alpha\eta c}{4\pi} \nonumber\\
 & & L_\alpha=\frac{4\pi}{q_\alpha c
 n_\alpha}\frac{\mu_\alpha}{\eta}\equiv~~{\rm{Prandtl~~number.}}
\end{eqnarray}The minimization integral is then obtained as,
\begin{eqnarray}
 {\cal{I}}&=&\int
 \left[(\nabla\times{\bf{B}})\cdot(\nabla\times{\bf{B}})+\frac{4\pi
 n q L_i}{c}{\bf{\omega}}_i\cdot{\bf{\omega}}_i\right]d\tau \nonumber\\
 &+&\lambda_i\int[\nabla\times({\bf{B}}+L_i{\bf{\omega}}_i)]\cdot{\bf{\Omega_i}}
 d\tau-\lambda_e\int(\nabla\times{\bf{B}})\cdot{\bf{\Omega}}_e
 d\tau
 \end{eqnarray}

\noindent where $L_i = 4\pi\mu_i/ecn_i \eta$ is the magnetic
Prandtl number for ions, $\mu_i$ is the ion viscosity and $\eta$
is the resistivity. We assumed the electron mass $m_e\rightarrow
0. $ The first two terms in the above integral represent the ohmic
and viscous dissipation rates while the second two integrals
represent the ion and electron generalized helicity dissipation
rates. $\lambda_i$ and $\lambda_e$ are the corresponding Lagrange
undetermined multipliers.The Euler-Lagrange equations are written
as,\begin{eqnarray} &
&\nabla\times\nabla\times{\bf{B}}+\frac{(\lambda_i+\lambda_e)q}{c}
\nabla\times){\bf{B}}+\frac{\lambda_i}{2}\left( \frac{q
L_i}{c}+m_i\right)(\nabla\times){\bf{\omega}}_i=\nabla\psi;
\label{mageq}\\
& &\nabla\times\nabla\times{\bf{\omega}}_i+\frac{2\pi n
q}{m_i\lambda_i c}\nabla\times{\bf{\omega}}_i+\frac{1}{2 m_i
L_i}\left( \frac{q L_i}{c}+m_i\right)(\nabla\times)^2{\bf{B}}=
\nabla\chi;\label{floweq}
\end{eqnarray}
where $\psi$ and $\chi $ are arbitrary gauge functions, with
$\nabla^2\psi=\nabla^\chi=0 $. Following Turner (1986) we
can write the exact solution of the coupled equations
(\ref{mageq}) and (\ref{floweq}) as a linear superposition of two
CK eigenfunctions, which are the solutions of
$\nabla\times{\bf{Y}}_k=\lambda_k{\bf{Y}}_k $ i.e.,


\begin{figure}[ht]
\begin{center}
\epsfig{file=singlearcflow.eps, width=10.cm, height=8.cm}
\end{center}
\caption{\small\label{fig:flow} Flow lines in the $x$-$z$ plane
for $\gamma=1.5$ (\textit{from Bhattacharya et al., 2007}).}
\end{figure}

\begin{equation}
{\bf{B}}={\bf{Y}}_1+\alpha{\bf{Y}}_2 \qquad
{\bf{\omega}}_i={\bf{Y}}_1+\beta{\bf{Y}}_2,
\end{equation}
and $\alpha$ and $\beta$ are constants and quantify the non
force-free part, to  be obtained from some boundary conditions.



\subsection{Application to solar arcade problem}



One of the important features  of this formalism is the presence
of a nontrivial  flow and the coupling between the  magnetic field
and the  flow. The solutions of the corresponding Euler-Lagrange
equation highlight the role of nonzero plasma flow in securing a
steady MDR relaxed state  for an open and externally driven
system.  This model has been very successful in describing the
arcade structures observed in Solar corona [Bhattacharya et al.,
2007]. To model solar arcade-type magnetic fields, we make some
simplified assumptions and consider a 2-D magnetic field in the
$x$-$z$ plane with the $x$ direction lying on the photosphere
surface and $z$ is the vertical direction. With appropriate
boundary conditions and after some lengthy algebra, the solutions
for the components of the magnetic fields are obtained as,




\begin{figure}[ht]
\begin{center}
\epsfig{file=solartwo.eps, width=10.cm, height=8.cm}
\end{center}
\caption{\small\label{fig:arcade}Magnetic field lines in the
$x$-$z$ plane corresponding to (i) single arcade structure for
$\gamma=1.5$ in left panel and (ii) triple arcade structure for
$\gamma=3.36$ in right panel (\textit{from Bhattacharya et al.,
2007}).}
\end{figure}


\begin{eqnarray}
& &B_x=\sqrt{\gamma^2-1}\left[2\cos\gamma
x-\frac{2\cos\gamma}{\cos\sqrt{\gamma^2-1}}\cos(x\sqrt{\gamma^2-1})\right]
e^{-z\sqrt{\gamma^2-1}};\nonumber\\
& &B_y=0;\\
& &B_z=-\left[2\gamma\sin\gamma
x-\frac{2\sqrt{\gamma^2-1}\cos\gamma}{\sin\sqrt{\gamma^2-1}}\sin(x\sqrt{\gamma^2-1})\right]
e^{-z\sqrt{\gamma^2-1}}\nonumber;
\end{eqnarray}
with $\gamma = \sqrt{(k^2+\lambda^2)/\lambda^2}$ ($\lambda$ is the
eigenvalue and $k$ is the  wave number associated with the
solution). The same procedures can be adopted for the flow
pattern. In figure \ref{fig:arcade}  the magnetic field lines in
the$x$-$z$ plane is plotted for two different values of the
parameter $\gamma$. The arcade structures are obtained by solving
the field line equations $\frac{dz}{dx}=\frac{B_z}{B_x}$.  It is
interesting to note that the arcade structure changes with the
change of the single parameter $\gamma$. Figure \ref{fig:flow}
plots the flow lines for $\gamma = 1.5$. In cylindrical geometry
the solutions of eqn. (\ref{mageq}) can be written in terms of a
superposition of Bessel functions. Such solutions yield 3-D Flux
ropes, which are plotted in the Figure
({\ref{fig:fluxrope}).


\begin{figure}[ht]
\begin{center}
\epsfig{file=mdr_fr.eps, width=10.cm, height=8.cm}
\end{center}
\caption{\small\label{fig:fluxrope} Plot of 3-D Fluxrope obtained
from the solution eqn.(\ref{mageq}).}
\end{figure}



\subsection{Non-force Free Coronal Magnetic Field Extrapolation
Based on MDR:}  An outstanding problem in solar physics is to
derive the coronal magnetic field structure from measured
photospheric and/or chromospheric magnetograms. A new approach to
deriving three-dimensional non-force free coronal magnetic field
configurations from vector magnetograms, based on the MDR
principle has been recently developed [Hu, Dasgupta, \& Choudhary
(2007)]. In contrast to the principle of minimum energy, which
yields a linear force-free magnetic field, the MDR gives a more
general non-force free magnetic field with flow. The full
MDR-based approach requires two layers of vector magnetic field
measurements. Its exact solution can be expressed as the
superposition of two linear force-free fields with distinct
$\lambda$ parameters, and one potential field. The final solution
is thus decomposed into three linear force-free extrapolations,
with bottom boundary conditions derived from the measured vector
magnetograms, at both photospheric and chromospheric levels. The
semi-analytic test case shown in Figure~\ref{FIGHu} illustrates
the feasibility and the high performance of this method. This
extrapolation is easy to implement and much faster than most other
nonlinear force-free methods [Hu \& Dasgupta 2007]. It also takes
full advantage of multiple layer solar magnetic field
measurements, thus gives a more realistic result.



\begin{figure}[ht]
\begin{center}
\epsfig{file=extrapolation.eps, width=13.cm, height=8.cm}
\end{center}
\caption{3D view of magnetic field lines in a box simulating the
typical extrapolation domain. (a) exact solution; (b) numerical
extrapolation result. Bottom image at height 0 shows the normal
field component distribution with gray scale indicating strongly
negative (black) to positive (white) field. The great ressemblance
between the exact and the numerical solutions justifies the
feasibility and high performance of this method (Hu, Dasgupta, \&
Choudhary 2007).}
\label{FIGHu}
\end{figure}


\section{Conclusion}

We have demonstrated that a relaxation model based on the minimum
dissipation of energy can successfully explain some of the salient
observations both in laboratory and space plasma. An MDR based
model model in its simpler form is capable of yielding basic
characteristics of most of the laboratory plasma confinement
devices, and a generalization of this model for two fluid with
external drive can couple magnetic field and flow and predict
solar arcade structures. The above findings definitely points out
that the MDR relaxed states applied to astrophysical plasmas is a
worthy case for further investigations. In future we plan to apply
this model to investigate magnetic field structure in the sun with
an aim to understanding coronal mass ejection.

\section{Acknowledgement}


We gratefully acknowledge partial supports from NASA LWS grant
NNX07A073G, and NASA grants NNG04GF83G, NNG05GH38G, NNG05GM62G, a
Cluster University of Delaware subcontract BART372159/SG, and NSF
grants ATM0317509, and ATM0428880.









\begin{thereferences}{9}







\bibitem{Amari2000}
Amari, T., and Luciani, J. F., Phys. Rev. Letts., {\bf 84}, 1196,
(2000)

\bibitem{antoni1983}
Antoni, V., S. Martini, S. Ortolini and R. Paccagnella, in
\textit{Workshop on Mirror-based and Field-reversed approaches to
Magnetic Fusion}, Int. School of Plasma Physics, Varenna, Italy,
p. 107, (1983)

\bibitem{bhatta2000}
Bhattacharya, R., M.S. Janaki and B. Dasgupta,  Phys. Plasmas,{\bf
7}, 4801, (2000)
\bibitem{bhatta2001}
Bhattacharya, R., M.S. Janaki and B. Dasgupta,
 Phys. Lett. A , {\bf 291A}, 291, (2001)

\bibitem{bhatta2003}
Bhattacharyya, R., M. S. Janaki and B. Dasgupta, Plasma Phys.
Contr. Fusion, {\bf 45}, 63  (2003).


 \bibitem{bhatta2004}
Bhattacharya, R., and M. S. Janaki, Phys. Plasmas, {\bf 11}, 5615
(2004)

 \bibitem{bhatta2007}
Bhattacharya, R., M. S. Janaki, B. Dasgupta, G. P. Zank, Solar
Phys.  {\bf 240},  63, (2007)


\bibitem{bhiskamp}
 Biskamp, D. 2003, Magnetohydrodynamic Turbulence,
 {\it Cambridge University Press}.



\bibitem{bouzat}
Bouzat, S. and R. Farengo, J. Plasma Phys., {\bf 72}, 443, (2006)

\bibitem{chandra1957}
Chandrasekhar, S. \& P. C. Kendall, Astrophys. J., {\bf 126}, 457
(1957).

\bibitem{chandra1958}
Chandrasekhar, S. \& L. Woltjer, Proc. Natl. Acad. Sci. USA, {\bf
44}, 285 (1958).

\bibitem{dasgupta1995}
Dasgupta, B., T. Sato, T Hayashi, K. Watanabe and T-H Watanabe,
Trans. Fusion Tecnol., {\bf 27}, 374, (1995)

\bibitem{dasgupta1998}
Dasgupta, B., P. Dasgupta, M. S. Janaki, T. Watanabe and T. Sato,
    Phys. Rev. Lett., {\bf 81}, 3144, (1998)


\bibitem{dasgupta2002}
Dasgupta, B., M. S. Janaki, R. Bhattacharya, P. Dasgupta, T.
Watanabe, and T. Sato,     Phys. Rev. E , {\bf E 65}, 046405,
(2002)

\bibitem{farengo1994}
Farengo, R. and J. R Sobehart, Plasma Phys. Contr. Fus., {\bf 36},
465 (1994)

\bibitem{farengo1995}
Farengo, R. and J. R. Sobehart. Phys. Rev. E, {\bf 52}, .2102
(1995)
\bibitem{farengo2002}
Farengo, R. and K. I. Caputi, Plasma Phys.  Contr. Fus., {\bf 44},
1707 (2002)


\bibitem{guo2005}
Guo, H. Y.., A. L. Hoffman, L. C. Steinhauer, and K. E. Miller,
Phys. Rev. Letts, {\bf 95}. 175001, (2005)

\bibitem{guo2005}
Guo, H. Y., A. L. Hoffman, L. C. Steinhauer, K. E. Miller, and R.
D. Milroy,  Phys. Rev. Letts, {\bf 97}. 175001, (2006)

\bibitem{hu2007}
Hu, Q., Dasgupta, B., \& Choudhary, D.P. 2007,  AIP CP, {\bf 932},
376, (2007)

\bibitem{hu2007b}
Hu, Q., and Dasgupta, B. , \textit{Sol. Phys.},
{\bf 247}, 87, (2008)


\bibitem{hart}
Hart, G. W., C. Chin-Fatt, A. W. deSilva, G. C. Goldenbaum, R.
Hess and R. S. Shaw, Phys. Rev. Lett., {\bf 51}, 1558, (1983)

\bibitem{lahiri}
Lahiri, S., A. N. S. Iyenger, S. Mukhopadhyay and R. Pal, Nucl.
Fusion, {\bf 36}, 254, (1996)



\bibitem{mont}
Montgomery D. and L. Phillips, Phys. Rev. A, {\bf 38}, 2953 (1988)




\bibitem{ono}
Ono, Y., A. Morita, and M. Katsurai and M Yamada, Phys. Fluids,
{\bf B5}, 3691, (1993)


\bibitem{orto}
Ortolani, S. and D. D. Schnack, Magnetohydrodynamics of plasma
relaxation, \textit{World Scientific,} (1993)

\bibitem{prigo}
Prigogine, I., \textit{Etude Thermodynamique des
Ph\'{e}nom\`{e}nes Irreversibles, Editions Desoer, Li\`{e}ge,
(1946) }

\bibitem{rayl}
Rayleigh, Proc. Math Soc. London,  {\bf 363}, 357 (1873)

\bibitem{rosen}
Rosenbluth, M. N. and M. N. Bussac, Nucl. Fusion, {\bf 19}, 489,
(1979)


\bibitem{sato}
Sato, T.  and Complexity Simulation Group: S. Bazdenkov,
B.Dasgupta, S. Fujiwara, A. Kageyama, S. Kida, T. Hayashi, R.
Horiuchi, H. Muira, H. Takamaru, Y. Todo, K. Watanabe and T.-H.
Watanabe, Phys. Plasmas. {\bf 3}, 2135 (1996)

\bibitem{dastgeer}
Shaikh, D., B. Dasgupta, G. P., Zank, and Q. Hu,
\textit{Phys. Plasmas}, {\bf 15}, 012306, (2008)


\bibitem{taylor}
  Taylor, J. B., Phys. Rev. Lett. {\bf 33}, 139 (1974)

\bibitem{turner}
Turner, L.,  IEEE Trans. Plasma Sc. {\bf 14} , 849, (1986)


\bibitem{tus}
Tuszewski, M., Nuc. Fusion, {\bf 28}, 2033, (1988)

\bibitem{}
 Watanabe, T.-H.,  T. Sato, and T. Hayashi, Phys. Plasmas {\bf 4},
(1997)


\bibitem{yamada}
Yamada, H., K. Kusano, Y. Kamada, H. Utsumi, Z. Yoshida and N.
Inoue, Nucl. Fusion, {\bf 27}, 1169, (1987)


\bibitem{yoshida}
Yoshida, Z., S. Ishida, K. Hattori, Y. Murakami and J. Morikawa,
J. Phys. Soc. Jpn., {\bf 55}, 554, (1986)


\bibitem{zhu}
Zhu, S.,  R. Horiuchi, and T. Sato, Phys. Rev. E, {\bf 51}, 6047
(1995).


\end{thereferences}









\label{lastpage}
\end{document}
