Improved Finite Volume Droplet Trajectories for Icing Simulation

During a Lagrangian icing simulation a large number of droplet trajectories are calculated to determine the water catch, and as a result it is important that this procedure is as rapid as possible. In order to arrive at a method with minimum complexity, a finite volume representation is developed for streamlines and extended to incorporate the equations of motion for a droplet, with all cells being crossed in a single timestep. However, since cells vary greatly in size, the method must be implicit to avoid an awkward stability restriction which would otherwise degrade performance. An implicit method is therefore implemented by performing iterations to solve for the crossing of each CFD cell, so that the droplet motion is tightly coupled to the underlying flow. By crossing every cell in a single step, and by using the mesh connectivity to track the droplet motion between cells, any need for costly searches or containment checks is eliminated and the resulting method is efficient. The implicit system is solved using functional iteration, which is feasible for the droplet system (which can be stiff) by using a particular factorisation.