<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Apr 5, 2017 at 1:13 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="HOEnZb"><div class="h5">Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> writes:<br>
<br>
> On Wed, Apr 5, 2017 at 12:03 PM, Jed Brown <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>> wrote:<br>
><br>
>> Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> writes:<br>
>> > As a side note, I think using FV to solve an elliptic equation should be<br>
>> a<br>
>> > felony. Continuous FEM is excellent for this, whereas FV needs<br>
>> > a variety of twisted hacks and is always worse in terms of computation<br>
>> and<br>
>> > accuracy.<br>
>><br>
>> Unless you need exact (no discretization error) local conservation,<br>
>> e.g., for a projection in a staggered grid incompressible flow problem,<br>
>> in which case you can use either FV or mixed FEM (algebraically<br>
>> equivalent to FV in some cases).<br>
>><br>
><br>
> Okay, the words are getting in the way of me understanding. I want to see<br>
> if I can pull something I can use out of the above explanation.<br>
><br>
> First, "locally conservative" bothers me. It does not seem to indicate what<br>
> it really does. I start with the Poisson equation<br>
><br>
>   \Delta p = f<br>
><br>
> So the setup is then that I discretize both the quantity and its derivative<br>
> (I will use mixed FEM style since I know it better)<br>
><br>
>   div  v = f<br>
>   grad p = v<br>
><br>
> Now, you might expect that "local conservation" would give me the exact<br>
> result for<br>
><br>
>   \int_T p<br>
><br>
> everywhere, meaning the integral of p over every cell T.<br>
<br>
</div></div>Since when is pressure a conserved quantity?<br>
<br>
In your notation above, local conservation means<br>
<br>
  \int_T (div v - f) = 0<br>
<br>
I.e., if you have a tracer moving in a source-free velocity field v<br>
solving the above equation, its concentration satisfies<br>
<br>
  c_t + div(c v) = 0<br>
<br>
and it will be conserved element-wise.<br>
</blockquote></div><br>But again that seems like a terrible term. What that statement above means is that globally</div><div class="gmail_extra">I will have no loss, but the individual amounts in each cell are not accurate to machine error,</div><div class="gmail_extra">they are accurate to discretization error because the flux is only accurate to discretization error.</div><div class="gmail_extra"><br></div><div class="gmail_extra">   Matt<br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>