<html><body><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; color: #000000"><div>Hi Neks,<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>I have come up with a new problem regarding the mesh deformation. <br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>I was trying to make a rigid motion near the walls of a 3D cylinder . I took the advantage of ocyl2.usrfile to make that. But I am getting a jacobian error when I am runnning the simulation.</div><div><br data-mce-bogus="1"></div><div>I slightly changed the boundary condition part in the ocyl2.usr<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div> nface = 2*ndim<br> do e=1,nelv<br> do f=1,nface<br><br>c Set Dirichlet for mesh velocity on all non-interior boundaries<br> if ( (cbc(f,e,1).ne.'E ' .<strong>or. (cbc(f,e,1).ne.'P ' ) )</strong> <br> $ call facev(msk,e,f,z0,nx1,ny1,nz1)<br><br>c Set inhomogeneous Dirichlet data on cylinder<br> if (cbc(f,e,1).eq.'mv ') then<br> call fcaver(xavg,xm1,e,f)<br> if (xavg.gt.-1.2) call facev(tmp,e,f,m1,nx1,ny1,nz1)<br> endif<br> enddo<br> enddo<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Since I use a 3D domain the right and left walls are Periodic bcs. So I also added one more condition (in bold letters above) to make the smoothing function.</div><div><br data-mce-bogus="1"></div><div>I feel it is correct but it throwing error once it starts running.</div><div><br data-mce-bogus="1"></div><div>I made the 3D geometry by using n2to3 on ocyl.rea file and gave PeriodicBoundary conditions. <br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Could you please tell me how to solve this in 3D<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Thank you <br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Sijo GEORGE<br data-mce-bogus="1"></div></div></body></html>