%                      Code_Saturne version 1.3
%                      ------------------------
%
%     This file is part of the Code_Saturne Kernel, element of the
%     Code_Saturne CFD tool.
% 
%     Copyright (C) 1998-2007 EDF S.A., France
%
%     contact: saturne-support@edf.fr
% 
%     The Code_Saturne Kernel is free software; you can redistribute it
%     and/or modify it under the terms of the GNU General Public License
%     as published by the Free Software Foundation; either version 2 of
%     the License, or (at your option) any later version.
% 
%     The Code_Saturne Kernel is distributed in the hope that it will be
%     useful, but WITHOUT ANY WARRANTY; without even the implied warranty
%     of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
%     GNU General Public License for more details.
% 
%     You should have received a copy of the GNU General Public License
%     along with the Code_Saturne Kernel; if not, write to the
%     Free Software Foundation, Inc.,
%     51 Franklin St, Fifth Floor,
%     Boston, MA  02110-1301  USA
%
%-----------------------------------------------------------------------
%
\programme{bilsc2}

\vspace{1cm}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Fonction}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

  Dans ce sous-programme, appel\'e par \fort{codits} et \fort{turbke}, sont calcul\'ees les contributions au bilan explicite des termes
convectifs et diffusifs reconstruits
(sur maillages non orthogonaux et si l'utilisateur le d\'esire) du second membre
d'une \'equation de convection/diffusion pour un scalaire $a$. Ces termes s'\'ecrivent de
facon g\'en\'erique \footnote{Ils interviennent notamment dans le second
membre du syst\`eme en incr\'ements pour la cellule $I$ de l'\'etape de pr\'ediction des vitesses : $\mathcal{EM}(\delta\vect{u}^{k+1},I) =
\mathcal{E}(\vect{u}^{n+1/2,k},I)$  ({\it cf}. \fort{navsto} pour plus de d\'etails).}:

\begin{equation}\label{Base_Bilsc2_eq_continue}
\begin{array}{ll}
\mathcal{B_{\mathcal{\beta}}}((\rho\,\vect{u})^{n},a) 
&=\underbrace{-\dive(\,(\rho \vect{u})^{n}a)}_{\text {partie convective}}
+\underbrace{\dive(\,\beta\,\grad a)}_{\text{partie diffusive}}\\
\end{array}
\end{equation}

avec $\rho$, $\vect{u}$, $\beta$ et $a$ variables connues \`a l'instant ${t^n}$. 

