> It struck me however, please correct me if I am wrong, that since we
> working with global shape function (fe), the integration
> (fe(old),fe(new))_new domain is merely loop over all cells and its
> simple provide to fe(old) the new domains quadrature points!
Yes, but it's not that simple because you need to find which cell of the
old mesh and where within that cell a quadrature point lies.
This is only simple if the new meshes has originated from the old one by
refinement/coarsening. It is very hard if the two meshes are defined on
different domains (for example because the domain changes shape in each
time step).
> But here comes the obstacle: if the quadrature point for new domain is
> outside the domain for the old, does it automatically give zero for
> fe(old)(new quadrature) or i must do a GridTool->is_within check?
You must check.
> Further, for large deformations one should simply drop the overlap
> functions, the part i am interested is the inertia force preventing
> immediate reaction to the fluid, so then i could actually have a if
> check, doing a quasistatic calculation then ask after, is the
> deformation to large? At that point i choose either step-23, or step-28,
> or my version for large deformations, drop the overlap terms.
I don't know what the overlap terms are you are talking about.
Independently of that, I really believe you should in a first step try to
take the old solution to the new grid using the SolutionTransfer class.
Everything else is too complicated to do in the first step.
Best
W.
-------------------------------------------------------------------------
Wolfgang Bangerth email: [EMAIL PROTECTED]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
