Computational study of high order numerical schemes for fluid-structure interaction in gas dynamics.

Abstract

Solving the fluid-structure interaction (FSI) problems is particularly challenging. This is because the coupling of the fluid and structure may require different solvers in different points of the solution domain, and with different mesh requirements. In this thesis, a partitioned approach is considered. Two solvers are employed to deal with each part of the problem (fluid and structure), where the interaction process is realized via exchanging information from the fluid-structure interface in a staggered fashion. One of the advantages of this approach is that we can take advantage of the existing algorithms that have been used for solving fluid or structural problems, which leads a reduction in the code development time, Hou (2012). However, it requires careful implementation so that spurious results in terms of stability and accuracy can be avoided. We found that most fluid-structure interaction computations through a staggered approach are based on at most second order time integration methods. In this thesis we studied the performance of some high order fluid and structure dynamic methods, when applied in a staggered approach to an FSI problem in a structure prediction way by combining predictors with time integration schemes to obtain stable schemes. Nonlinear Euler equations for gas dynamics were investigated and the analysis was realized through the piston problem. An adapted one-dimensional high order finite volume WENO₃ scheme for nonlinear hyperbolic conservation laws-Dumbser (2007a), Dumbser et al. 2007b)-was considered and a numerical flux was proposed. The numerical results of the proposed method show the non-oscillatory property when compared with traditional numerical methods such as the Local Lax-Friedrichs. So far to our knowledge, the WENO₃₋ as proposed in this work- has not been applied to FSI problems. Thus, it was proposed to discretize the fluid domain in space, and in order to adapt it to a moving mesh was reformulated to couple with an Arbitrary Lagrangian Eulerian (ALE) approach. To integrate in time the structure we started by using Newmark schemes as well as the trapezoidal-rule backward differentiation formulae of order 2 (TR − BDF2). Two study cases were carried out by taking into account the transient effects on the fluid behaviour. In the first case, we only consider the structural mass in the dynamic coupled system and in the second case, a quasi-steady fluid was considered. In order to test the performance of the structural solvers, simulations were carried out, firstly, without the contribution of fluid mass, and then a comparative study of the performance of various structure solvers in a staggered approach framework were realized in order to study the temporal accuracy for the partitioned fluid-structure interaction coupling. For a quasi-steady fluid case, the oscillation frequency of the coupled system was successfully estimated using the TR-BDF2 scheme, and the coupled system was solved for various Courant numbers in a structural predictor fashion. The results showed better performance of the TR-BDF2 scheme. Newmark’s schemes as well as the TR-BDF2 are only second order accuracy. However, the Newmark (average acceleration) is traditionally preferred by researchers as a structure solver in a staggered approach for FSI problems, although higher order schemes do exist. van Zuijlen (2004), in his partitioned algorithm proposed the explicit singly diagonal implicit Runge-Kutta (ESDIRK) family of schemes of order 3 to 5 to integrate both fluid and structure. Therefore in this work, these schemes were considered and applied as structural solvers. Their performance was studied through numerical experiments, and comparisons were realized with the performance of the traditional Newmark’s schemes. The results show that although their computational cost is high, they present a high order of accuracy.