<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Dec 5, 2017 at 9:12 AM, Yann Jobic <span dir="ltr"><<a href="mailto:yann.jobic@univ-amu.fr" target="_blank">yann.jobic@univ-amu.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div text="#000000" bgcolor="#FFFFFF">
<p>Hello,</p>
<p>I tried the master of the git directory, as you have merged the
branch.</p>
<p>It works ! Many thanks for the fix.<br>
</p>
<p></p></div></blockquote><div>Great! Feel free to mail any suggestion for improvements.</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div text="#000000" bgcolor="#FFFFFF"><p>Best regards,</p><div><div class="h5">
<p>Yann<br>
</p>
<br>
<div class="m_135884221591823617moz-cite-prefix">Le 22/11/2017 à 18:51, Matthew Knepley
a écrit :<br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div class="gmail_extra">
<div class="gmail_quote">On Wed, Nov 22, 2017 at 12:39 PM,
Yann Jobic <span dir="ltr"><<a href="mailto:yann.jobic@univ-amu.fr" target="_blank">yann.jobic@univ-amu.fr</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br>
<br>
I've found a strange behavior when looking into a bug for
the pressure convergence of a simple Navier-Stokes problem
using PetscFE.<br>
<br>
I followed many examples for labeling boundary faces. I
first use DMPlexMarkBoundaryFaces, (label=1 to the faces).<br>
I find those faces using DMGetStratumIS and searching 1 as
it is the value of the marked boundary faces.<br>
Finally i use DMPlexLabelComplete over the new label.<br>
I then use :<br>
ierr = PetscDSAddBoundary(prob, DM_BC_ESSENTIAL, "in",
"Faces", 0, Ncomp, components, (void (*)(void)) uIn,
NWest, west, NULL);CHKERRQ(ierr);<br>
in order to impose a dirichlet condition for the faces
labeled by the correct value (west=1, south=3,...).<br>
<br>
However, the function "uIn()" is called in all the
Elements containing the boundary faces, and thus impose
the values at nodes that are not in the labeled faces.<br>
Is it a normal behavior ? I then have to test the position
of the node calling uIn, in order to impose the good
value.<br>
I have this problem for a Poiseuille flow, where at 2
corner Elements i have a zero velocity dirichlet condition
(wall) and a In flow velocity one.<br>
</blockquote>
<div><br>
</div>
<div>I believe I have fixed this in
knepley/fix-plex-bc-multiple which should be merged soon.
Do you know how to merge that branch and try?</div>
<div><br>
</div>
<div> Thanks,</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">
The pressure is then very high at the corner nodes of
those 2 Elements.<br>
Do you think my pressure problem comes from there ? (The
velocity field is correct)<br>
<br>
Many thanks,<br>
<br>
Regards,<br>
<br>
Yann<br>
<br>
PS : i'm using those runtime options :<br>
-vel_petscspace_order 2 -pres_petscspace_order 1 \<br>
-ksp_type fgmres -pc_type fieldsplit -pc_fieldsplit_type
schur -pc_fieldsplit_schur_fact_type full \<br>
-fieldsplit_velocity_pc_type lu
-fieldsplit_pressure_ksp_rtol 1.0e-10
-fieldsplit_pressure_pc_type jacobi<br>
<br>
<br>
---<br>
L'absence de virus dans ce courrier électronique a été
vérifiée par le logiciel antivirus Avast.<br>
<a href="https://www.avast.com/antivirus" rel="noreferrer" target="_blank">https://www.avast.com/antiviru<wbr>s</a><br>
<br>
</blockquote>
</div>
<br>
<br clear="all">
<div><br>
</div>
-- <br>
<div class="m_135884221591823617gmail_signature" data-smartmail="gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div>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><br>
</div>
<div><a href="http://www.caam.rice.edu/%7Emk51/" target="_blank">https://www.cse.buffalo.edu/~<wbr>knepley/</a><br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
<br>
</div></div><pre class="m_135884221591823617moz-signature" cols="72"><span class="HOEnZb"><font color="#888888">--
___________________________
Yann JOBIC
HPC engineer
IUSTI-CNRS UMR 7343 - Polytech Marseille
Technopôle de Château Gombert
</font></span><span class=""><a href="https://maps.google.com/?q=5+rue+Enrico+Fermi&entry=gmail&source=g">5 rue Enrico Fermi</a>
13453 Marseille cedex 13
Tel : (33) 4 91 10 69 43
Fax : (33) 4 91 10 69 69
</span></pre><span class="">
<div id="m_135884221591823617DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2"><br> <table style="border-top:1px solid #d3d4de">
<tbody><tr>
<td style="width:55px;padding-top:18px"><a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=emailclient" target="_blank"><img src="https://ipmcdn.avast.com/images/icons/icon-envelope-tick-round-orange-animated-no-repeat-v1.gif" alt="" width="46" height="29" style="width:46px;height:29px"></a></td>
<td style="width:470px;padding-top:17px;color:#41424e;font-size:13px;font-family:Arial,Helvetica,sans-serif;line-height:18px">Garanti sans virus. <a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=emailclient" style="color:#4453ea" target="_blank">www.avast.com</a> </td>
</tr>
</tbody></table>
<a href="#m_135884221591823617_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2" width="1" height="1"> </a></div></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>