<div dir="ltr">Why does this petsc-maint message have a messed up Reply-To: header?<div><br></div><div>   Matt</div><div><br><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Yann Jobic</b> <span dir="ltr"><<a href="mailto:yann.jobic@univ-amu.fr">yann.jobic@univ-amu.fr</a>></span><br>Date: Wed, Nov 22, 2017 at 12:39 PM<br>Subject: [petsc-users] 2 Dirichlet conditions for one Element in PetscFE<br>To: PETSc <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>><br><br><br>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>
<br>
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>
</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>