\documentclass[11pt,a4paper]{article}
\topmargin=-2cm
\oddsidemargin=-1cm
\evensidemargin=-1cm
\textwidth=16.5cm
\textheight=24cm
\pagestyle{plain}
\begin{document}
\thispagestyle{empty}
\begin{center}
\textbf{\Large Finite volume methods}
\medskip
\textbf{\Large for non smooth solutions of diffusion models;}
\medskip
\textbf{\Large Application to imperfect contact problems}
\vskip 0.5cm
\textbf{\large Philippe ANGOT}
\end{center}
\vskip 0.5cm
\begin{small}
\begin{center}
{\sl Universit\'e de la M\'editerran\'ee} \\
{\sl I.R.P.H.E. Ch\^ateau-Gombert, Equipe Math\'ematiques Num\'eriques} \\
{\sl La Jet\'ee, Technop\^ole de Ch\^ateau-Gombert} \\
{\sl 38 rue F. Joliot Curie, 13451 Marseille Cedex 20 - France} \\
\smallskip
{\sl E-mail : angot@marius.univ-mrs.fr} \\
{\sl URL : http://www-irphe.univ-mrs.fr}
\end{center}
\end{small}
\vskip 2cm
\textbf{Abstract --}
\bigskip
The concept of crack or contact resistance is sometimes introduced
empirically for diffusion pro\-blems with imperfect contact, e.g.
thermal or electrical contact resistance, or also hydraulic resistance
of fracture for flows in porous media. The objective is to take
account of fault lines or too thin layers compared to the largest
scale under study. Hence, at this scale the solution is
\textit{discontinuous}.
Here, we generalize this concept and we formulate for such diffusion
problems a mathematical model with discontinuous coefficients which
includes an \textit{``imperfect contact'' transmission condition}
linking the diffusive flux with the jump of the solution. It ensures
the well-posedness of the associated elliptic or parabolic problems
and we prove the global solvability within a variational framework.
For the numerical solution, we propose an original finite volume
method based on the introduction of \textit{``crack resistances''}
at the faces of the control volumes. In that way, we show how to satisfy
both conservativity and consistency of the numerical fluxes.
The convergence of this scheme is established within the
framework of finite volume theory. Some error estimates in ${\cal
O}(h)$ are obtained for 1D or 2D rectangular meshes provided the
coefficients are piecewise of class ${\cal C}^{1}$ and the exact
solution is piecewise of class ${\cal C}^{2}$. Besides, the condition
number of the discrete operator is analysed and the improvement due to
the diagonal preconditioning is estimated.
Various numerical results illustrate the capabilities and the
efficiency of this method.
\end{document}