# Diff of /branches/doubleplusgood/doc/user/darcyflux.tex

revision 4344 by jfenwick, Wed Feb 27 03:42:40 2013 UTC revision 4345 by jfenwick, Fri Mar 29 07:09:41 2013 UTC
# Line 70  However this is a very fast way of calcu Line 70  However this is a very fast way of calcu
70
71  \subsubsection{Global Postprocessing \label{SEC DARCY POST}}  \subsubsection{Global Postprocessing \label{SEC DARCY POST}}
72  An improved flux recovery can be achieved by solving a modified version of equation~\ref{DARCY PROBLEM SIMPLE FLUX}  An improved flux recovery can be achieved by solving a modified version of equation~\ref{DARCY PROBLEM SIMPLE FLUX}
74  \label{DARCY PROBLEM POST FLUX}  \label{DARCY PROBLEM POST FLUX}
75  \kappa^{-1}_{ij} u_{j} -  \kappa^{-1}_{ij} u_{j} -
76  (\lambda \cdot u_{k,k} )_{,i}=  (\lambda \cdot u_{k,k} )_{,i}=
# Line 188  which is available in the \ExampleDirect Line 188  which is available in the \ExampleDirect
188     u,p=mypde.solve(u0=x[1]*[0., -1.], p0=0)     u,p=mypde.solve(u0=x[1]*[0., -1.], p0=0)
189     saveVTK("u.vtu",flux=u, pressure=p)     saveVTK("u.vtu",flux=u, pressure=p)
190  \end{python}  \end{python}
191  In the example the pressure is fixed to the initial presure \var{p0} on the right half of the top face.  In the example the pressure is fixed to the initial pressure \var{p0} on the right half of the top face.
192  The normal flux is set on all other faces. The corresponding values for the flux are set by the initial value  The normal flux is set on all other faces. The corresponding values for the flux are set by the initial value
193  \var{u0}.    \var{u0}.

Legend:
 Removed from v.4344 changed lines Added in v.4345