Hi all
Sorry for the re-post. As Ruhollah suggested, I would like to repeat my
question in detail.
The equation I faced is a simple 3-D transport equation with smooth initial
condition on a relatively large domain \Omega:
u_t + div(a u) =0
In my analysis, u becomes singular as t goes to infinity. So I would like to
apply the adaptive finite element method (Crank-Nicolson scheme for time
discretization). My question is how to preserve the mass (\int_\Omega u dx
=const) during this approach.
Basically, people apply the L^2- projection to transform numerical solution
u_{n-1} on an "old" mesh to a new refined mesh and make
the result the solution at time n-1. For this transport equation, I can not
use the L^2-projection as it violate the conservation law.
Does anyone have some good ideas or suggestions for this case?
Best,
Liang
_______________________________________________
