Hi,
I have just created the issue [1] aiming at implementing the slip/no-slip boundary conditions. Is there anybody here with some experience in this field? Several ways come to my mind:
adapt the existing linear combination BC (LCBC) (see [2]) to nonzero right-hand side - the LCBC allow imposing \sum_{i=1..N} a_i u_i = 0 for the dofs u_i and coefficients a_i. We need \sum_{i=1,2,3} n_i u_i = g.
enforce the condition via the Lagrange multipliers, leading to a saddle point system,
?
Any comments would be appreciated! r.
[1] http://code.google.com/p/sfepy/issues/detail?id=82 [2] http://sfepy.kme.zcu.cz/cgi-bin/moin.cgi/SfePy/Examples#RigidBodyMotion
Robert Cimrman wrote:
Hi,
I have just created the issue [1] aiming at implementing the slip/no-slip boundary conditions. Is there anybody here with some experience in this field? Several ways come to my mind:
adapt the existing linear combination BC (LCBC) (see [2]) to nonzero right-hand side - the LCBC allow imposing \sum_{i=1..N} a_i u_i = 0 for the dofs u_i and coefficients a_i. We need \sum_{i=1,2,3} n_i u_i = g.
enforce the condition via the Lagrange multipliers, leading to a saddle point system,
?
Any comments would be appreciated! r.
[1] http://code.google.com/p/sfepy/issues/detail?id=82 [2] http://sfepy.kme.zcu.cz/cgi-bin/moin.cgi/SfePy/Examples#RigidBodyMotion
I have been faster than you - the case u_i n_i = 0 is implemented :)
cheers, r.
On Wed, Mar 11, 2009 at 8:22 AM, Robert Cimrman <cimr...@ntc.zcu.cz> wrote:
Robert Cimrman wrote:
Hi,
I have just created the issue [1] aiming at implementing the slip/no-slip boundary conditions. Is there anybody here with some experience in this field? Several ways come to my mind:
adapt the existing linear combination BC (LCBC) (see [2]) to nonzero right-hand side - the LCBC allow imposing \sum_{i=1..N} a_i u_i = 0 for the dofs u_i and coefficients a_i. We need \sum_{i=1,2,3} n_i u_i = g.
enforce the condition via the Lagrange multipliers, leading to a saddle point system,
?
Any comments would be appreciated! r.
[1] http://code.google.com/p/sfepy/issues/detail?id=82 [2] http://sfepy.kme.zcu.cz/cgi-bin/moin.cgi/SfePy/Examples#RigidBodyMotion
I have been faster than you - the case u_i n_i = 0 is implemented :)
You are faster than light!
O.
participants (2)
-
Ondrej Certik -
Robert Cimrman