> hi Roy, > thanx for the explaination. But how did you solve ((u > * grad)u, v)_Omega? It's a square term. I heard there > are some other methods, streamline, least square FEM > ... I would like to hear your comments.
By square I guess you mean asymmetric? The convection term is asymmetric and causes problems when convection is stronger than diffusion, e.g. Cell Reynolds numbers (u*h/nu) exceed 2. In that case you need to do something line the streamline/upwind Petrov-Galerkin method, Galerkin-Least squares, etc... We sidestep that problem in ex13 by keeping the cell Reynolds numbers low enough that it is not an issue. -Ben ------------------------------------------------------------------------- This SF.net email is sponsored by: Microsoft Defy all challenges. Microsoft(R) Visual Studio 2008. http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
