Development of a Shed-Ice Trajectory Simulation in FENSAP-ICE

A numerical approach to the trajectory simulation of break-up ice in the in-flight icing code FENSAP-ICE is investigated. At each time step, the displacement and rotation of the moving domains is computed from a 6-DOFs analysis of the forces and moments. The moving domains are amalgamated into the fixed background mesh by a hole-cutting and stitching algorithm, producing a continuous unstructured hybrid mesh, eliminating interpolation between domains and ensuring that the fluxes are fully conserved across the entire mesh. To maintain good load balancing on parallel computers, the finite element CFD flow solver uses an efficient parallel iterative matrix solver and domain decomposition to partition the computational domain into equal-sized subdomains. Test cases used to validate and demonstrate the features of the computational algorithm are shown.