[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [ff3d-users] Variational form of potential flow with free surface on
Re: [ff3d-users] Variational form of potential flow with free surface on mesh
Thu, 08 Jan 2009 20:11:17 +0000
Thunderbird 18.104.22.168 (X11/20090105)
last May you gave me a pointer to how to work with periodic boundary
conditions in ff3d. 8 months later I have found the time to try it.....
I'm not able to get it to work, my file is attached at the end of this
message can you see what I am doing wrong?
I would like to understand how it is implemented in the code have been
trying to understand code but can't manage to follow it. Can you point
me in the right directions, so far I have been looking through
MeshPeriodizer.cpp and VerticesCorrespondance.cpp , Mesh.cpp and
I can see that the nodes of the end boundary of the structured mesh are
mapped to the nodes of the start boundary but I am struggling to see how
this is carried forwards to the solver.
Are there any papers which describe the method you have followed?
thanks for your help
vertex a = (0,0,0);
vertex b = (10,10,10);
vector n = (25,50,10);
mesh M = structured(n,a,b);
//this next line is the only change I made to the static solution
M = periodic(M,0:1);
solve(phi) in M
Thomas Ward wrote:
Many thanks Stephane,
I hope to try the periodic conditions next week, its not totally clear
to me how the periodic mesh becomes a periodic boundary condition but I
will have a try first and then ask if I can't work it out.
Stephane Del Pino wrote:
Yes this is not a big deal in ff3d. Assume that M is your structured
I'm about to try using periodic BC, I note from the mailing list that
you were playing around with this a couple of years ago, do you have any
sample files or other documentation or pointers to the source code?
M = periodic(M,0:1,2:3,4:5);
will create a triperiodic mesh where
0 represents xmin,
1 represents xmax,
2 represents ymin,
3 represents ymax,
4 represents zmin,
5 represents zmax,
M = periodic(M,2:3);
would create a periodic mesh in the direction y. Note that up to now M
must be a structured mesh.
- Re: [ff3d-users] Variational form of potential flow with free surface on mesh,
Thomas Ward <=