<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 am bit stuck on this occasion for past 2 days, the code below shows a part of the example called ocyl2.usr (from subroutine my_base_meshv(basev,h1,h2,rhs,msk,tmp) ). This code is to generate a base velocity (where basev=1 at wall of object and 0 at the boundary of the domain). I assume that I have a blending coefficient h1.(using Paul's method or volume of elements method). This subroutine working perfectly for the 2D cases. That is I have the base velocity ranges between 0 and 1. But when I tried it on an aero dynamic foil case(its a 3D case where just 1 element in Y direction), the value of base velocity does not even exceed 0. Almsot all values are like -0.0021564... Could you please tell me is there anything that I have to change in this subroutine to make it working. </div><div><br data-mce-bogus="1"></div><div>Tha value of m1 and z0 are just explicity giving the base velocity to the elements right? for example on the W or M.V conditions basev becomes 1 and on boundaries of the domain it becomes 0.0 <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><br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div> m1 = -1.<br> z0 = 0.0<br><br> nface = 2*ndim<br><br><br> do e=1,nelv<br> do f=1,nface<br><br>c Set Dirichlet for mesh velocity on all boundaries<br><br> if (cbc(f,e,1).ne.'E ') call facev(msk,e,f,z0,nx1,ny1,nz1)<br><br>c Inhomogeneous Dirichlet on cylinder only<br><br> if (cbc(f,e,1).eq.'mv ') then<br> call fcaver(xavg,xm1,e,f)<br><br><br> if (xavg.gt.-3.0) call facev(tmp,e,f,m1,nx1,ny1,nz1)<br> endif<br> enddo<br> enddo<br><br> call axhelm (rhs,tmp,h1,h2,1,1)<br><br> tol = 1.e-4<br> imsh = 1<br> ifield = 1<br> call hmholtz('mshv',basev,rhs,h1,h2,msk,vmult,imsh,tol,200,1)<br> call sub2(basev,tmp,n)<br><br>c call outpost(basev,basev,basev,basev,basev,' ')<br>c stop<br><br><br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div><div>Regards<br data-mce-bogus="1"></div><div>Sijo<br data-mce-bogus="1"></div><div><br data-mce-bogus="1"></div></div></body></html>