<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Apr 16, 2014 at 5:38 PM, Sang pham van <span dir="ltr"><<a href="mailto:pvsang002@gmail.com" target="_blank">pvsang002@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Yes, It does converge toward 0-flow when I refine the mesh !<div><br></div><div>In the code, coordinates are set by </div>
<div>  ierr = DMDASetUniformCoordinates(da_prop,0.0+0.5*dx,1.0-0.5*dx,0.0+0.5*dy,1.0-0.5*dy,0.,0);CHKERRQ(ierr);<br>
</div><div>Does it set coordinate for elements center here?</div><div>When I want to try non-uniform grid, should I just modify coordinates attached with the DM, or I need more implementation?</div></div></blockquote><div>
<br></div><div>Yes, I think that should be enough.</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>
Many thanks.</div>
<div>Sang</div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Apr 16, 2014 at 6:23 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div>On Wed, Apr 16, 2014 at 4:48 PM, Sang pham van <span dir="ltr"><<a href="mailto:pvsang002@gmail.com" target="_blank">pvsang002@gmail.com</a>></span> wrote:<br>


<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Jed,<div><br></div><div>I modified the ex43 code to enforce no-slip BCs on all boundaries. </div><div>


I run the code with volume force (0,-1) and isoviscosity. The expected result is Vx = Vy = 0 everywhere, and linearly decreasing pressure (from to to bottom). </div>
<div>In the attached is plot of velocity field and pressure, so there is still a (light) flow in middle of the domain. Do you know why the solution is that, and what should I do to get the expected result?</div></div></blockquote>


<div><br></div></div><div>It sounds like it is due to discretization error. Your incompressibility constraint is not verified element-wise</div><div>(I think ex43 is penalized Q1-Q1), so you can have some flow here. Refine it and see if it converges toward</div>


<div>0 flow.</div><div><br></div><div>   Matt</div><div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">
<div>Thank you.</div><div><br></div><div>Sang</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Mar 26, 2014 at 2:38 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>Sang pham van <<a href="mailto:pvsang002@gmail.com" target="_blank">pvsang002@gmail.com</a>> writes:<br>
<br>
> Hi Dave,<br>
> I guess you are the one contributed the ex42 in KSP's examples. I want to<br>
> modify the example to solve for stokes flow driven by volume force in 3D<br>
> duct. Please help me to understand the code by answering the following<br>
> questions:<br>
><br>
> 1. Firstly, just for confirmation, the equations you're solving are:<br>
>      \nu * \nabla \cdot \nabla U - \nabla P = 0 and<br>
<br>
</div>For variable viscosity, it must be formulated as in the example:<br>
<br>
  \nabla\cdot (\nu D U) - \nabla P = 0<br>
<br>
  where  D U = (\nabla U + (\nabla U)^T)/2<br>
<div><br>
>      \nabla \cdot U = 0<br>
><br>
>       where U = (Ux,Uy,Uz), \nu is variable viscosity?<br>
><br>
>  2. Are U and P defined at all nodes? (I googled the Q1Q1 element, it looks<br>
> like a box element with U and P defined at 8 corners).<br>
<br>
</div>Yes.<br>
<div><br>
>  3. Are nodes' coordinate defined though the DA coordinates?<br>
<br>
</div>Yes, though they are set to be uniform.<br>
<div><br>
>  4. How can I enforce noslip BC, and where should I plug in volume force?<br>
<br>
</div>Enforce the Dirichlet condition for the entire node.<br>
</blockquote></div><br></div>
</blockquote></div></div><span><font color="#888888"><br><br clear="all"><span class="HOEnZb"><font color="#888888"><div><br></div>-- <br>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
</font></span></font></span></div></div>
</blockquote></div><br></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>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>