This thesis describes the development of a two-dimensional, high-order, fluid-structure interaction (FSI) solver. The well-established spectral difference (SD) method is used for spatial discretization of the Euler equations over deforming, unstructured quadrilateral grids. The Geometric Conservation Law (GCL) is incorporated into the conservative Euler equations, before discretization. After simplification, the equations reduce to a form, in the computational domain, identical to the equations in the physical domain. In this form, the equations can be integrated implicitly in time, without the requirement of any additional source term, to guarantee free-stream preservation. The fluid and structure sub-systems are individually integrated in time using the explicit first stage, single diagonal, diagonally implicit Runge-Kutta (ESDIRK) method. As the first step to solving the coupled, non-linear Euler equations, implicit in time, we linearize the governing equations. The resulting linearized simultaneous equations are then solved sequentially using lower-upper symmetric Gauss-Seidel (LU-SGS) relaxation sweeps. The fluid and structure sub-systems are loosely coupled and the coupling term is integrated in time using an explicit RK method, resulting in an implicit-explicit (IMEX) RK coupling. The spatial accuracy and the free-stream preserving ability of the solver are demonstrated by testing a supersonic, isentropic vortex in a curved channel. Next, the temporal accuracy of the solver is established using an Euler vortex propagation test case. It is also demonstrated that the four-stage ESDIRK is capable of handling time-steps 50 times larger than the four-stage explicit RK. In each of these cases, third- and fourth-order SD for spatial discretization and second-order backward difference (BDF2) and third-order, four stage ESDIRK for time integration were tested. Since the loose (explicit) FSI coupling restricts permissible structural deformation, we limit ourselves to small harmonic oscillations resulting from linearized perturbed Euler equations. The interaction between a linear piston and an inviscid, compressible fluid is simulated to demonstrate that the IMEX coupling does not contaminate the spatial or temporal accuracy of the implemented high-order methods. Through rigorous testing, this development is expected to lay a foundation for a powerful computational framework for various fluid-structure interaction problems.
|Advisor:||Liang, Chunlei, Baginski, Frank|
|Commitee:||Lee, James D., Wickenheiser, Adam|
|School:||The George Washington University|
|Department:||Mechanical and Aerospace Engineering|
|School Location:||United States -- District of Columbia|
|Source:||MAI 55/05M(E), Masters Abstracts International|
|Keywords:||Computational fluid dynamics, Fluid structure interaction, High-order methods, Implicit-explicit time-stepping|
Copyright in each Dissertation and Thesis is retained by the author. All Rights Reserved
The supplemental file or files you are about to download were provided to ProQuest by the author as part of a
dissertation or thesis. The supplemental files are provided "AS IS" without warranty. ProQuest is not responsible for the
content, format or impact on the supplemental file(s) on our system. in some cases, the file type may be unknown or
may be a .exe file. We recommend caution as you open such files.
Copyright of the original materials contained in the supplemental file is retained by the author and your access to the
supplemental files is subject to the ProQuest Terms and Conditions of use.
Depending on the size of the file(s) you are downloading, the system may take some time to download them. Please be