CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   Main CFD Forum (http://www.cfd-online.com/Forums/main/)
-   -   Problems with SUPG body force term (http://www.cfd-online.com/Forums/main/10669-problems-supg-body-force-term.html)

 FEM question January 21, 2006 18:51

Problems with SUPG body force term

Dear friends,

I have been getting problems with the SUPG term of the body forces (BF) in a program I'm implementing. I have a SUPG/PSPG finite element incompressible Navier-Stokes program where every time I "turn on" the evaluation of the SUPG term of the body force the program blows up. This problem doesn't happen when I use only the Galerkin term and the result seems to be ok. I know that the SUPG term of the forcing terms introduces more non-linearity in the problem, so, my body forces evaluation is before the predictor-multicorrector loop and are linearized with the time lagged velocity field. So, what could be wrong? My SUPG BF vector is defined by:

f_SUPG = int[ tau_SUPG * (u grad w) * f]

int = integral operator grad = gradient operator tau_SUPG = stabilization parameter u = velocity vector (lagged in time for linearization purposes) w = weighting function f = force vector

being f and w defined by:

f = sum(N_a*f_a) and w = sum(N_b*c_b)

would reach: