Augmented Lagrangian method for FEMCFD
Hi, i am looking for someone who uses the augmented lagrangian method with a finite element method to simulate incompressible fluid flow. I am actualy using a penalty method and i am wondering if this method is better. Thank you best regards

Re: Augmented Lagrangian method for FEMCFD
Hi,
I have used Augmented Lagrangian Method (ALM) for solving incompressible fluid flows (especially nonNewtonian) in the past. The discretization method was a Finite Volume method, but that is likely not to make any difference with your FEM since at the end you get symmetric matrix system that looks like [ A B ](v)=(f) [ Bt 0 ](p) (g) (the two lines are together and represent a single matrix system, e.g. it corresponds to : A.v+B.p=f and Bt.v=g) where A is the velocity matrix, B pressure matrix, v velocity vector, p pressure vector, f and g the second members of your system With the ALM, you replace this problem by the following one [ A+rBBt B ](v)=(f+rBg) [ Bt 0 ](p) (g) where r stands for the Lagrangian parameter. The larger r, the faster the convergence is achieved. But since in the algorithm, you are to solve matrix subsystems like : (A+rBBt)(x)=b Those systems become badconditioned as r increases. Therefore, if u use a direct method to solve this subsystem (as for instance Cholevski factorization), then you don't have to take care of the magnitude of r. On the other hand, if u choose to solve this subsystem using iterative solver (conjugate gradient), then the larger r, the poorer the convergence rate. Finally, the ALM is a very powerful method provided you can factortize the augmented matrix A+rBBt once at the beginning of the procedure in its Cholevski form. Otherwise, you can still use it but it requires to pay attention to the magnitude of r. Do not hesitate to contact me if you need additional informations I hope it might help u Best regards, Anthony 
All times are GMT 4. The time now is 12:36. 