Hi Stefan,<br><br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
What do you mean by "remapping a periodic-side domain in userdat2"?<br></blockquote><div><br>I build a mesh using genbox with dimensions [0,1]^3.  The BC are periodic in x.  Then in userdat 2 I remap the (x,y) cross-section from a square to an annulus:<br>



<br>      call rescale_x(xm1,zero,pi2)<br>      call rescale_x(ym1,zero,one)<br><br>      n = nx1*ny1*nz1*nelv   ! Turn box into annulus<br>      do i=1,n<br>         th = xm1(i,1,1,1)<br>         r  = r0 + (r1-r0)*(1.-ym1(i,1,1,1))<br>



         xm1(i,1,1,1) = r*cos(th)<br>         ym1(i,1,1,1) = r*sin(th)<br>      enddo<br><br>This gives me an annulus centered on the z-axis with interior radius r0 and exterior radius r1.  (I tried centering it on the y-axis instead with the periodic boundary lying in the xy-plane, but the result was the same.)  My velocity BC in what are initially the y and z directions are all SYM, so I need to use the stress formulation.<br>



<br>This procedure worked in 2D for no-slip or stress-free BC, but in 3D it is not working for stress-free BC.  I can visualize the 3D geometry, and the annular cylinder appears to be properly constructed, but the fluid solution blows up.<br>

<br>I understand now that my fractional-annulus cases blow up because non-cyclic periodic BC at the theta-ends are incompatible with stress-free walls, but this should not be an issue for the full 2š annulus, even in 3D.  Can you think of a reason that non-cyclic periodic BC would work in 2D but not in 3D?<br>

<br></div><blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204, 204, 204);padding-left:1ex">
Sure, even with IFCYCLIC=F you'll get a result. But is probably not<br>
the one you had in mind, simply because the BCs are different i.e.<br>
periodic vs cyclic.<br></blockquote><div><br>Ah yes, you`re right that the results are actually different.<br><br>Thanks a lot,<br><br>David<br></div>
</div>