A Simulator for 2D Rigid Body Motion in Viscous Incompressible Flow
Date Issued
2014
Date
2014
Author(s)
Lin, Chien-Da
Abstract
Complicated fluid structure interaction (FSI) problems commonly appear in industrial applications and in the nature. This work is aimed at developing a 2D numerical tool for simulating the two-way interactions between a rigid body of arbitrary shape and a viscous incompressible fluid flow. We apply the spectral element method for the spatial discretization because of its geometric flexibility and its high accuracy. In order to capture the arbitrary shape of the rigid body, a Lagrangian grid is generated and attached to the rigid body in addition to an Eulerian grid for the whole flow field. The no-slip boundary condition at the rigid body surface is enforced by the penalty method, in which a penalization term is added in the momentum equation inside the solid domain. Besides, we propose a so-called “cut-cell method” to calculate the FSI force. This method guarantees the momentum conservation and enhances the efficiency of the integration over the solid domain.
The freely-falling motion of a circular, elliptic, square, triangular, or rod type cylinder under gravity is simulated. In the circular cylinder case, the simulation results confirm the validity and the accuracy of the presently proposed simulator at low Reynolds numbers and low solid-to-fluid density ratios. The elliptic cylinder case shows the same phenomenon observed by others. The simulation results of the elliptic, square, and triangular cylinder cases at low density ratios are also in accordance with those in the literature. The simulation result of the rod-type cylinder case however does not agree with the result in the literature, mostly probably because it has a density ratio as high as 2.
Subjects
流固耦合
寬頻元素法
沉浸邊界法
分胞法
自由落體
Type
thesis
File(s)![Thumbnail Image]()
Loading...
Name
ntu-103-R01522109-1.pdf
Size
23.54 KB
Format
Adobe PDF
Checksum
(MD5):b39455ae141e28783a33d27ccd29fbdb