/[escript]/branches/doubleplusgood/doc/user/darcyflux.tex
ViewVC logotype

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

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

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}
73  adding the the gradient of the divergence of the flux:  adding the gradient of the divergence of the flux:
74  \begin{equation}\label{DARCY PROBLEM POST FLUX}  \begin{equation}\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

  ViewVC Help
Powered by ViewVC 1.1.26