Adaptivity in CEA’s fluid elements in EUROPLEXUS

The present work completes the implementation of adaptivity routines by extending them to CEA's fluid finite elements both in 2D (TRIA and CAR1) and in 3D (TETR and CUBE). The CAR1 is treated like other 2D quadrilaterals (Q41L, FL24) as far as geometrical issues are concerned. In addition to the solid case, the activation of adaptivity for fluids requires the suitable treatment of transport terms which arise in the Eulerian or ALE forms of the governing equations. For the CEA's fluid finite elements mentioned above (TRIA, CAR1, TETR and CUBE) this is done in routines tr2me.ff (for the 2D case) and tr3me.ff (for the 3D case), respectively. Therefore, most modifications for the current implementation are concentrated in those two routines. Actually, a special version of the routines is written, valid for the mesh adaptive case.