[Nek5000-users] Mesh Deformation in 3D

[nek5000-users at lists.mcs.anl.gov]

Hi Neks, 

I have come up with a new problem regarding the mesh deformation. 

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. 

I slightly changed the boundary condition part in the ocyl2.usr 



nface = 2*ndim 
do e=1,nelv 
do f=1,nface 

c Set Dirichlet for mesh velocity on all non-interior boundaries 
if ( (cbc(f,e,1).ne.'E ' . or. (cbc(f,e,1).ne.'P ' ) ) 
$ call facev(msk,e,f,z0,nx1,ny1,nz1) 

c Set inhomogeneous Dirichlet data on cylinder 
if (cbc(f,e,1).eq.'mv ') then 
call fcaver(xavg,xm1,e,f) 
if (xavg.gt.-1.2) call facev(tmp,e,f,m1,nx1,ny1,nz1) 
endif 
enddo 
enddo 

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. 

I feel it is correct but it throwing error once it starts running. 

I made the 3D geometry by using n2to3 on ocyl.rea file and gave PeriodicBoundary conditions. 


Could you please tell me how to solve this in 3D 

Thank you 

Sijo GEORGE 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180717/d3289cf5/attachment.html>