دانلود رایگان مقاله شبیه ساز عنصر طیفی سلول تابع جسم صلب سیال

Abstract

In this work, we propose an immersed-boundary-type simulation method for the two-way coupling problems between a rigid body and a fluid. In this simulation tool, the entire fluid-solid domain is treated as an incompressible fluid with non-uniform density and the no-slip boundary condition at the rigid body surface is enforced by the penalization method. The fluid solver is developed in use of the spectral element method for the spatial discretization and the mixed explicit/implicit scheme for the temporal discretization. An additional Lagrangian mesh is employed and attached to the rigid body in order to trace the rigid body and to perform the area integration over the solid domain. Besides, a so-called sub-cell scheme is developed to smooth the discontinuity at the fluid–solid interface. The validity and accuracy of the proposed simulation method were examined well by applying it to the sedimentation problems of circular, triangular, square, as well as elliptic cylinders in a channel. An accuracy of 2nd order was observed, probably due to the use of triangular Lagrangian elements and consequently a piecewise-linear approximation of the rigid body shape and also due to the 2nd-order interpolation involved in the sub-cell scheme.

5. Conclusion

A 2D simulation method was developed to solve the two-way coupling problems between an incompressible viscous fluid and a rigid body. It is developed based on the spectral element method and the immersed boundary concept. The whole system is modeled as a nonuniform density flow and the rigid body motion is enforced by the penalization method. A Lagrangian mesh is employed to trace the rigid body and a so-called sub-cell approach is constructed to deal with the discontinuity at the solid–fluid interface. The proposed scheme was found to be 2nd-order accurate and its geometric flexibility was confirmed. With this geometric flexibility arising from the use of the spectral element method as well as the use of a Lagrangian mesh, the present scheme is applicable to flow induced by rigid bodies of arbitrary shapes in channels with complicated wall boundaries. Numerical instabilities associated with high density-ratio and high Reynoldsnumber flow need further improving however.