Dynamic Mesh Update in FEniCS for Lagrangian Electrospinning Simulation

[KIMO_D]

I am solving a system of transient coupled PDEs: x(t, s), y(t, s), and z(t, s), along with an ODE for stress(t), in a Lagrangian framework for electrospinning. The model tracks a particle’s trajectory, and initially, the fiber follows a straight-line path before transitioning into a series of loops due to whipping.

I need to know if it’s possible to update the computational mesh dynamically after each calculation step to accommodate the evolving fiber shape

[dokken]

Yes, the mesh geometry can evolve throughout the simulation.
See for instance:

Moving submesh diverges from its designated trajectory in a two-phase problem

This is not the right thing to do, as this assumes alot of things about the ordering of degrees of freedom. Here is a minimal example that maps a displacement field from a mesh onto a submesh: """ Map degrees of freedom to the mesh geometry nodes and corresponding submesh. Author: Jørgen S. Dokken and Henrik N.T. Finsberg SPDX-License-Identifier: MIT """ from mpi4py import MPI import dolfinx import numpy as np def vertex_to_dof_map_vectorized(V): """Create a map from the vertices of th…