• 한국해양공학회지
• /
• 제27권2호
• /
• pp.53-58
• /
• 2013

Recently, some fluid-structure interaction (FSI) problems involving the fluid impact loads interacting with structures, such as sloshing, slamming, green-water, etc., have been considered, especially in the ocean engineering field. The governing equations for both an elastic solid model and flow model were originally derived from similar continuum mechanics principles. In this study, an elastic model based on a particle method, the MPS method, was developed for simulating the FSI problems. The developed model was first applied to a simple cantilever deflection problem for verification. Then, the model was coupled with the fluid flow model, the PNU (Pusan National University modified)-MPS method, and applied to the numerical investigation of the coupling effects between a cantilever and a mass of water, which has variable density, free-falling to the end of the cantilever.

• Structural Engineering and Mechanics
• /
• 제46권2호
• /
• pp.197-212
• /
• 2013

With more and more high-rise building being constructed in recent decades, bluff body flow with high Reynolds number and large scale dimensions has become an important topic in theoretical researches and engineering applications. In view of mechanics, the key problems in such flow are high Reynolds number turbulence and fluid-solid interaction. Aiming at such problems, a parallel fluid-structure interaction method based on socket parallel architecture was established and combined with the methods and models of large eddy simulation developed by authors recently. The new method is validated by the full two-way FSI simulations of 1:375 CAARC building model with Re = 70000 and a full scale Taipei101 high-rise building with Re = 1e8, The results obtained show that the proposed method and models is potential to perform high-Reynolds number LES and high-efficiency two-way coupling between detailed fluid dynamics computing and solid structure dynamics computing so that the detailed wind induced responses for high-rise buildings can be resolved practically.

• Wind and Structures
• /
• 제20권3호
• /
• pp.423-448
• /
• 2015

In this paper the unsteady fluid-structure interaction (FSI) problems with large structural displacement are solved by partitioned solution approaches in the arbitrary Lagrangian-Eulerian finite element framework. The incompressible Navier-Stokes equations are solved by the characteristic-based split (CBS) scheme. Both a rigid body and a geometrically nonlinear solid are considered as the structural models. The latter is solved by Newton-Raphson procedure. The equation governing the structural motion is advanced by Newmark-${\beta}$ method in time. The dynamic mesh is updated by using moving submesh approach that cooperates with the ortho-semi-torsional spring analogy method. A mass source term (MST) is introduced into the CBS scheme to satisfy geometric conservation law. Three partitioned coupling strategies are developed to take FSI into account, involving the explicit, implicit and semi-implicit schemes. The semi-implicit scheme is a mixture of the explicit and implicit coupling schemes due to the fluid projection splitting. In this scheme MST is renewed for interfacial elements. Fixed-point algorithm with Aitken's ${\Delta}^2$ method is carried out to couple different solvers within the implicit and semi-implicit schemes. Flow-induced vibrations of a bridge deck and a flexible cantilever behind an obstacle are analyzed to test the performance of the proposed methods. The overall numerical results agree well with the existing data, demonstrating the validity and applicability of the present approaches.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제12권1호
• /
• pp.85-101
• /
• 2020

The optimal parameters for the fluid-structure interaction analysis using the Smoothed Particle Hydrodynamics (SPH) for fluids and finite elements for structures, respectively, are explored, and the effectiveness of the simulations with those parameters is validated by solving several open surface fluid problems. For the optimization of the Equation of State (EOS) and the simulation parameters such as the time step, initial particle spacing, and smoothing length factor, a dam-break problem and deflection of an elastic plate is selected, and the least squares analysis is performed on the simulation results. With the optimal values of the pivotal parameters, the accuracy of the simulation is validated by calculating the exerted force on a moving solid column in the open surface fluid. Overall, the SPH-FEM coupled simulation is very effective to calculate the fluid-structure interaction. However, the relevant parameters should be carefully selected to obtain accurate results.

• International Journal of Fluid Machinery and Systems
• /
• 제6권2호
• /
• pp.87-93
• /
• 2013

In order to study the effect of the fluid-structure interaction (FSI) on the simulation results, the external characteristics and internal flow features of a diffuser pump were analyzed with a two-way flow solid coupling method. And the static and dynamic structure analysis of the blade was also caculated with the FEA method. The steady flow field is based on Reynolds Averaged N-S equations with standard $k-{\varepsilon}$ turbulent model, the unsteady flow field is based on the large eddy simulation, and the structure response is based on elastic transient structural dynamic equation. The results showed that the effect of FSI on the head prediction based on CFD really exists. At the same radius, the van mises stress on the nodes closed shroud and hub was larger than other nodes. A large deformation region existed near inlet side at the middle of blades. The strength of impeller satisfied the strength requirement with static stress analysis based on the fourth strength theory. The dynamic stress varied periodically with the impeller rotating. It was also found that the fundamental frequency of the dynamic stress is the rotating frequency and its harmonic frequency. The frequency of maximum stress amplitude at node 1626 was 7 times of the rotating frequency. The frequency of maximum stress amplitude at node 2328 was 14 times of the rotating frequency. No matter strength failure or fatigue failure, the root of blades near shroud is the key region to analyse.

• 한국해양공학회지
• /
• 제25권4호
• /
• pp.36-41
• /
• 2011

In the present paper, a direct forcing/fictitious domain (DF/FD) level set method is proposed to simulate the FSI (fluid-solid interaction) in two-phase flow. The main idea is to combine the direct-forcing/fictitious domain (DF/FD) method with the level set method in the Cartesian coordinates. The DF/FD method is a non-Lagrange-multiplier version of a distributed Lagrange multiplier/fictitious domain (DLM/FD) method. This method does not sacrifice the accuracy and robustness by employing a discrete ${\delta}$ (Dirac delta) function to transfer quantities between the Eulerian nodes and Lagrangian points explicitly as the immersed boundary method. The advantages of this approach are the simple concept, easy implementation, and utilization of the original governing equation without modification. Simulations of various water-entry problems have been conducted to validate the capability and accuracy of the present method in solving the FSI in two-phase flow. Consequently, the present results are found to be in good agreement with those of previous studies.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• pp.151-158
• /
• 2009

To implement the insects' flapping flight for developing flapping MAVs(micro air vehicles), the unsteady flow characteristics of the insects' forward flight is investigated. In this paper, two-dimensional FSI(Fluid-Structure Interaction) simulations are conducted to examine realistic flow features of insects' flapping flight and to examine the flexibility effects of the insect's wing. The unsteady incompressible Navier-Stokes equations with an artificial compressibility method are implemented as the fluid module while the dynamic finite element equations using a direct integration method are employed as the solid module. In order to exchange physical information to each module, the common refinement method is employed as the data transfer method. Also, a simple and efficient dynamic grid deformation technique based on Delaunay graph mapping is used to deform computational grids. Compared to the earlier researches of two-dimensional rigid wing simulations, key physical phenomena and flow patterns such as vortex pairing and vortex staying can still be observed. For example, lift is mainly generated during downstroke motion by high effective angle of attack caused by translation and lagging motion. A large amount of thrust is generated abruptly at the end of upstroke motion. However, the quantitative aspect of flow field is somewhat different. A flexible wing generates more thrust but less lift than a rigid wing. This is because the net force acting on wing surface is split into two directions due to structural flexibility. As a consequence, thrust and propulsive efficiency was enhanced considerably compared to a rigid wing. From these numerical simulations, it is seen that the wing flexibility yields a significant impact on aerodynamic characteristics.

• 정보저장시스템학회논문집
• /
• 제7권2호
• /
• pp.60-64
• /
• 2011

Chemical mechanical planarization is one of the core processes in fabrication of semiconductors, which are increasingly used for information storage devices like solid state drives. For higher data capacity in storage devices, CMP process is required to show ultimate precision and accuracy. In this work, 2-dimensional finite element models were developed to investigate the effects of the slurry particle impact on microscratch generation and the phenomena generated at pad-particle-wafer contact interface. The results revealed that no plastic deformation and corresponding material removal could be generated by simple impact of slurry particles under real CMP conditions. From the results of finite element simulations, it could be concluded that the pad-particle mixture formed in CMP process would be one of major factors leading to microscratch generation.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2009년도 정기 학술대회
• /
• pp.59-62
• /
• 2009

유체-고체 상호작용(FSI)은 산업전반에서 꼭 필요한 분야이면서도 쉽게 접근하기가 어려운 전산역학 분야의 난제 중 하나이다. 유체-고체 상호작용의 전산해석에서 유체와 고체 사이의 불일치 격자망을 어떻게 처리하는가는 매우 어렵고 민감한 부분이 된다. 운동학적 연속성과 계면을 따른 응력의 평형을 추적하기 위해 유체와 고체의 계면에서는 각각의 영역에서 해석된 물리량들을 다른 영역으로 정확히 전달해야 하는데 대부분의 유체-고체 상호작용의 문제들은 불일치 격자를 가지고 있기 때문에 불일치 격자망을 효과적으로 처리하는 수단이 필요하다. 그래서 넓은 분야에 걸쳐 적용 가능한 유체 고체 상호작용 문제에 대한 효과적인 해석방법의 제안이 큰 의미를 갖는다고 생각한다. 따라서 본 연구에서는 유체-고체 계면의 운동을 이동최소제곱 기반의 변절점 요소를 사용하여 모사함으로써 2차원 유체-고체의 상호작용(FSI)을 위한 새로운 접근방법을 제시하였다.

• 한국산학기술학회논문지
• /
• 제17권9호
• /
• pp.126-132
• /
• 2016

본 연구는 전산유체역학 기법을 이용하여 수소 생산 플랜트의 개질 튜브 공정가스와 버너 가스 온도에 따른 화학반응과 열변형 특성을 분석한다. 개질로 내부의 온도는 약 800 K 내지 1000 K 이상으로 고온으로 유지되기 때문에 튜브의 열변형 문제가 심각하게 발생할 수 있다. 따라서 개질로의 구조건전성을 평가하고 안정된 생산력을 가진 장비를 운영하기 위해서 반응과 열변형 특성에 대한 이해는 필수적이다. 본 연구는 상용 전산해석 코드(ANSYS Fluent/Mechanical V.13.0)를 사용하여, 대류, 전도 및 복사 열전달을 포함한 복합 열전달과 난류유동을 3차원적으로 해석하였다. 특히, 열유동 특성에 따른 연성해석(Fluid-Solid Interaction: FSI)를 수행하였으며 고온 버너가스와 공정가스 운전조건에 따른 반응 특성과 열변형 변화를 분석하였다. 수치해석 결과, 개질 공정가스와 버너 가스의 주입온도가 각각 200 K 감소하면, 수소생성량은 최대 약 4 배, 최소 약 2 배 감소한다. 또한, 공정가스와 버너 가스의 주입온도에 따라 열변형은 최대 약 20%, 최소 약 15% 감소한다.