An algorithm for simulating solid objects suspended in stratified flow Original Research Article

An efficient difference algorithm for computing directly deformationless solid objects suspended in stratified flow in 2D has been developed. The objects are represented by colour functions (or density functions) and predicted by a sharpness preserving scheme that is able to prevent the numerical diffusion across the sharp interface. Pressure distribution is then calculated by a unified solver and the solid object is treated as a mass of material of high sound speed. The motion of the solid object is decomposed into translation and rotation, and the force as well as the torque that cause change in the motion of the solid body are evaluated by an averaging calculation over the region occupied by the solid body. Calculations are conducted on a fixed grid system. Operations for reconstructing moving interfaces or dealing with inner boundary conditions are not necessary.