[P Sharma]

Hi,
I have to verify vof of openFoam with standard case of circular blob of fluid subjected to a periodic roatating flow field . I took a drop on 2D box and I have to give velocity field in x and y direction. Which is function of time and space u(x,y,t) and v(x,y,t).

u(x,y,t)= (1/π)*cos(πt/T)*((sin(πx))^2)*2πsin(πy)cos(πy)


v(x,y,t)= (1/π)*cos(πt/T)*((sin(πy))^2)*2πsin(πx)cos(πx)

Now, I want to know that how can I implement this initial condition in openFoam.

[Tushar@cfd]

Originally Posted by P Sharma

Hi,
I have to verify vof of openFoam with standard case of circular blob of fluid subjected to a periodic roatating flow field . I took a drop on 2D box and I have to give velocity field in x and y direction. Which is function of time and space u(x,y,t) and v(x,y,t).

u(x,y,t)= (1/π)*cos(πt/T)*((sin(πx))^2)*2πsin(πy)cos(πy)


v(x,y,t)= (1/π)*cos(πt/T)*((sin(πy))^2)*2πsin(πx)cos(πx)

Now, I want to know that how can I implement this initial condition in openFoam.

Are you sure that you want to implement it for initial condition? Or, is it a boundary condition?

You can apply custom boundaries using "groovyBC". Anyways try these..

-
Best Luck!

[P Sharma]

thank you for reply, I am sure I have to to implement it for initial condition........can you tell me how to use "groovyBC" ?

Originally Posted by Tushar@cfd

Are you sure that you want to implement it for initial condition? Or, is it a boundary condition?

You can apply custom boundaries using "groovyBC". Anyways try these..

-
Best Luck!

[kmefun]

http://openfoamwiki.net/index.php/Contrib/groovyBC

[ybapat]

Hello,

I think you can implement this though momentum source term. You have velocity variation, from this get acceleration variation. Use this to add momentum source.

-Yogesh