Hi All,
In a typical predictor-corrector scheme for flow problems, we define fields
u (velocity), p(pressure) and a(acceleration). While solving for (u,p)
together,
we can use a predictor corrector algorithm that essentially translates to
(M+gamma*dt*N(u)){delta_a}+G{delta_p} = {f} - M{a} - N(u){u} - G{p}
To do this,
1. I declare a (FE_Q<dim>(1), dim+1) system
2. Vectors: U (holding u and p), A (holding a and 0 (for p) ) and DELTA
(holding
delta_a and delta_p)
3. A few steps in the algorithm go as
a. u += (1-gamma)*dt*a [predictor]
b. u += gamma*dt*delta_a [corrector]
p += delta_p
and so on
My question is that: How can I find out the (a) part of A and add it to
(u) part of U
and so on.
I need this since there does not seem to be a block-solver present (I
might be wrong
here)
Thanks,
Badri
--
Badri Velamur Asokan
Research Assoc. Post-doctoral
Computer Science & Mathematics Division
Oak Ridge National Lab
Ph: 865-576-8784
_______________________________________________
