Investigation on meshfree particle methods for fluid structure interaction problems

The research aims to investigate the application of meshfree particle methods for computational modelling of the Fluid Structure Interaction problems with particular emphasis on flow around cylindrical objects and aerofoils. For this purpose, a solution scheme has been developed for solving incompressible, viscous Navier-Stokes (N-S) equations over meshfree particles. Spatial derivatives appearing in N-S equations are dealt with using radial basis functions infinite difference mode (RBF-FD). A comparative study has also been conducted between implicit and explicit in time solution schemes for N-S equations over meshfree nodes. Subsequently, a coupled meshfree and mesh-based solution scheme is proposed, over hybrid fluid grid, for incompressible, viscous flow around stationary as well as moving objects. The aim of this coupled solver is to provide efficiency and flexibility by combining the advantages of both meshfree and mesh-based methods. The coupled solution scheme suggests generating a body conformal meshfree nodal cloud around the solid body in the near field. A static Cartesian grid surrounds the meshfree cloud in the far field. The Meshfree nodes offer flexibility in dealing with solid motion by moving along the solid boundary without necessitating re-meshing. The Cartesian grid, on the other hand, provides improved performance by allowing faster computation owing to the use of efficient mesh based method. Flow equations, in Arbitrary Lagrangian-Eulerian (ALE) formulation, are solved using RBF-FD based scheme over moving meshfree nodes. Conventional infinite differencing is used over static Cartesian grid for flow equations in Eulerian formulation. The coupled solution scheme, on hybrid grid, is employed for closely coupled Fluid Structure Interaction problems. The equations for solid motion are solved using classical Runge-Kutta method. Close coupling between fluid and structural solvers is realized by a sub-iterative prediction-correction algorithm. In order to reduce computational overhead due to sub-iterations, only near field flow (in meshfree zone) is solved during inner iterations. Solution over full fluid domain is sought during outer (time step) iterations only, when the convergence at fluidsolid interface has already been reached. The solution scheme is also applied for high Reynolds number problems. For this purpose, a stabilization term is included in the flow equations to suppress the spurious oscillations. The stabilization term is derived using vimomentum balance equation over control volume and applying higher order Taylor series expansion of momentum flux and fluid forces. In order to avoid ill-conditioning and accuracy problems related to RBF matrices in domains having varying nodal density, use of shape adaptive RBFs are proposed. In that, the shape parameter of the radial basis function is varied according to local nodal density. Moreover, adaptive sizing of influence domain has also been introduced to maintain suitable number of neighbouring particles. These adaptive techniques are found to be useful as they allow much finer nodal distribution at regions of interest enabling accurate capturing of flow gradients and leading to better results. The use of hybrid grid offers flexibility in dealing with moving boundaries. Moreover, in addition to allowing faster computing over Cartesian grid, it also enables using the reduced fluid domain during inner FSI iterations and therefore helps reduce the number of computations in the fluid domain during fluid-solid coupling. The solution scheme was tested for problems relating to flows around static as well as moving cylinders and aerofoils. Flow induced vibrations have been studied with one and two degrees of freedom. The results are found to be in good agreement with previous numerical work and experimental results.