Interface Tracking and Solid-Fluid Coupling Techniques with Coastal Engineering Applications

Abstract

Multi-material physics arise in an innumerable amount of engineering problems. A broadly scoped numerical model is developed and described in this thesis to simulate the dynamic interaction of multi-fluid and solid systems. It is particularly aimed at modelling the interaction of two immiscible fluids with solid structures in a coastal engineering context; however it can be extended to other similar areas of research. The Navier Stokes equations governing the fluids are solved using a combination of finite element (FEM) and control volume finite element (CVFE) discretisations. The sharp interface between the fluids is obtained through the compressive transport of material properties (e.g. material concentration). This behaviour is achieved through the CVFE method and a conveniently limited flux calculation scheme based on the Hyper-C method by Leonard (1991). Analytical and validation test cases are provided, consisting of steady and unsteady flows. To further enhance the method, improve accuracy, and exploit Lagrangian benefits, a novel moving mesh method is also introduced and tested. It is essentially an Arbitrary Lagrangian Eulerian method in which the grid velocity is defined by semi-explicitly solving an iterative functional minimisation problem. A multi-phase approach is used to introduce solid structure modelling. In this approach, solution of the velocity field for the fluid phase is obtained using Model B as explained by Gidaspow (1994, page 151). Interaction between the fluid phase and the solids is achieved through the means of a source term included in the fluid momentum equations. The interacting force is calculated through integration of this source term and adding a buoyancy contribution. The resulting force is passed to an external solid-dynamics model such as the Discrete Element Method (DEM), or the combined Finite Discrete Element Method (FEMDEM). The versatility and novelty of this combined modelling approach stems from its ability to capture the fluid interaction with particles of random size and shape. Each of the three main components of this thesis: the advection scheme, the moving mesh method, and the solid interaction are individually validated, and examples of randomly shaped and sized particles are shown. To conclude the work, the methods are combined together in the context of coastal engineering applications, where the complex coupled problem of waves impacting on breakwater amour units is chosen to demonstrate the simulation possibilities. The three components developed in this thesis significantly extend the application range of already powerful tools, such as Fluidity, for fluids-modelling and finite discrete element solids-modelling tools by bringing them together for the first time