Lagrangian movement algorithm


Move an object inside a mesh keeping the mesh conformity.

Input data

Nodal values of the displacement (velocity) over boundaries of reference 10.


  1. Extension of the velocity field:
    1. Creation of a submesh containing 20 layers of elements around the boundary to move;
    2. Propagation calling a linear elasticity solver with Dirichlet boundary conditions on the submesh;
  2. Dichotomy loop:
    1. Computation of the largest valid motion along the velocity field (with point relocation only);
    2. Mesh motion;
    3. Remeshing depending on the lagrangian mode (point relocation, point relocation + edge swapping, point relocation, edge swapping, point insertion and collapse);
  3. go back to step 1.

Comments are closed.