• Ocean Systems Engineering
• /
• v.8 no.4
• /
• pp.381-407
• /
• 2018

This paper presents a hybrid numerical approach, which combines a two-phase Navier-Stokes model (NS) and the fully nonlinear potential theory (FNPT), for modelling wave-structure interaction. The former governs the computational domain near the structure, where the viscous and turbulent effects are significant, and is solved by OpenFOAM/InterDyMFoam which utilising the finite volume method (FVM) with a Volume of Fluid (VOF) for the phase identification. The latter covers the rest of the domain, where the fluid may be considered as incompressible, inviscid and irrotational, and solved by using the Quasi Arbitrary Lagrangian-Eulerian finite element method (QALE-FEM). These two models are weakly coupled using a zonal (spatially hierarchical) approach. Considering the inconsistence of the solutions at the boundaries between two different sub-domains governed by two fundamentally different models, a relaxation (transitional) zone is introduced, where the velocity, pressure and surface elevations are taken as the weighted summation of the solutions by two models. In order to tackle the challenges associated and maximise the computational efficiency, further developments of the QALE-FEM have been made. These include the derivation of an arbitrary Lagrangian-Eulerian FNPT and application of a robust gradient calculation scheme for estimating the velocity. The present hybrid model is applied to the numerical simulation of a fixed horizontal cylinder subjected to a unidirectional wave with or without following current. The convergence property, the optimisation of the relaxation zone, the accuracy and the computational efficiency are discussed. Although the idea of the weakly coupling using the zonal approach is not new, the present hybrid model is the first one to couple the QALE-FEM with OpenFOAM solver and/or to be applied to numerical simulate the wave-structure interaction with presence of current.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1995.10a
• /
• pp.150-156
• /
• 1995

In the finite element analysis of metal forming processes, the updated Lagrangian approach has been widely and effectively used to simulate the non-steady state problems. However some difficulties have arisen from abrupt flow change as in extrusion through square dies. In the present work, a ALE(arbitrary Lagrangian-Eulerian) finite element formulation for deformation analysis are presented for rigid viscoplastic materials. The developed finite element program is applied to the analysis of square die extrusion of a square section. The computational results are compared with those from the updated Lagrangian finite element analysis.

• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.442-447
• /
• 2008

This paper presents a new approach to simulate fluid-solid interaction problems involving non-matching interfaces. The coupling between fluid and solid domains with dissimilar finite element meshes consisting of 4-node quadrilateral elements is achieved by using the interface element method (IEM). Conditions of compatibility between fluid and solid meshes are satisfied exactly by introducing the interface elements defined on interfacing regions. Importantly, a consistent transfer of loads through matching interface element meshes guarantees the present method to be an efficient approach of the solution strategy to fluid-solid interaction problems. An arbitrary Lagrangian-Eulerian (ALE) description is adopted for the fluid domain, while for the solid domain an updated Lagrangian formulation is considered to accommodate finite deformations of an elastic structure. The stabilized equal order velocity-pressure elements for incompressible flows are used in the motion of fluids. Fully coupled equations are solved simultaneously in a single computational domain. Numerical results are presented for fluid-solid interaction problems involving nonmatching interfaces to demonstrate the effectiveness of the methodology.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2410-2414
• /
• 2008

Human heart valves diseased by congenital heart defects, rheumatic fever, bacterial infection, cancer may cause stenosis or insufficiency in the valves. Treatment may be with medication but often involves valve repair or replacement (insertion of an artificial heart valve). Bileaflet mechanical heart valves (BMHVs) are widely implanted to replace the diseased heart valves, but still suffer from complications such as hemolysis, platelet activation, tissue overgrowth and device failure. These complications are closely related to both flow characteristics through the valves and leaflet dynamics. In this study, the physiological flow interacting with the moving leaflets in a bileaflet mechanical heart valve (BMHV) is simulated with a strongly coupled implicit fluid-structure interaction (FSI) method which is newly organized based on the Arbitrary-Lagrangian-Eulerian (ALE) approach and the dynamic mesh method (remeshing) in FLUENT. The simulated results are in good agreement with previous experimental studies. This study shows the applicability of the present FSI model to the complicated physics interacting between fluid flow and moving boundary.

• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.1487-1491
• /
• 2008

Detailed knowledge on the motion of blood cells flowing in micro-channels under simple shear flow and the influence of blood flow is essential to provide a better understanding on the blood rheological properties and blood cell aggregation. The microscopic behavior of red blood cell (RBCs) is numerically investigated using a fluid-structure interaction (FSI) method based on the Arbitrary-Lagrangian-Eulerian (ALE) approach and the dynamic mesh method (smoothing and remeshing) in FLUENT (ANSYS Inc., USA). The employed FSI method could be applied to the motions and deformations of a single blood cell and multiple blood cells, and the primary thrombogenesis caused by platelet aggregation. It is expected that, combined with a sophisticated large-scale computational technique, the simulation method will be useful for understanding the overall properties of blood flow from blood cellular level (microscopic) to the resulting rheological properties of blood as a mass (macroscopic).

• Wind and Structures
• /
• v.15 no.3
• /
• pp.189-208
• /
• 2012

Lock-in and drag amplification phenomena are studied for a flexible cantilever using a simplified fluid-structure interaction approach. Instead of solving the 3D domain, a simplified setup is devised, in which 2D flow problems are solved on a number of planes parallel to the wind direction and transversal to the structure. On such planes, the incompressible Navier-Stokes equations are solved to estimate the fluid action at different positions of the line-like structure. The fluid flow on each plane is coupled with the structural deformation at the corresponding position, affecting the dynamic behaviour of the system. An Arbitrary Lagrangian-Eulerian (ALE) approach is used to take in account the deformation of the domain, and a fractional-step scheme is used to solve the fluid field. The stabilization of incompressibility and convection is achieved through orthogonal quasi-static subscales, an approach that is believed to provide a first step towards turbulence modelling. In order to model the structural problem, a special one-dimensional element for thin walled cross-section beam is implemented. The standard second-order Bossak method is used for the time integration of the structural problem.

• Explosives and Blasting
• /
• v.37 no.3
• /
• pp.13-24
• /
• 2019

This paper was prepared to investigate the behavior of fragments in underwater torpedo explosion beneath a frigate or surface ship by using an explicit finite element analysis. In this study, a fluid-structure interaction (FSI) methodology, called the multi-material arbitrary Lagrangian-Eulerian (MM-ALE) approach in LS-DYNA, was employed to obtain the responses of the torpedo fragments and frigate hull to the explosion. The Euler models for the analysis were comprised of air, water, and explosive, while the Lagrange models consisted of the fragment and the hull. The focus of this modeling was to examine whether a worst-case fragment could penetrate the frigate hull located close (4.5 m) to the exploding torpedo. The simulation was performed in two separate steps. At first, with the assumption that the expanding skin of the torpedo had been torn apart by consuming 30% of the explosive energy, the initial velocity of the worst-case fragment was sought based on a well-known experimental result concerning the fragment velocity in underwater bomb explosion. Then, the terminal velocity of the worst-case fragment that is expected to occur before the fragment hit the frigate hull was sought in the second step. Under the given conditions, the possible initial velocities of the worst-case fragment were found to be very fast (400 and 1000 m/s). But, the velocity difference between the fragment and the hull was merely 4 m/s at the instant of collision. This result was likely to be due to both the tremendous drag force exerted by the water and the non-failure condition given to the frigate hull. Anyway, at least under the given conditions, it is thought that the worst-case fragment seldom penetrate the frigate hull because there is no significant velocity difference between them.

• Journal of the Society of Naval Architects of Korea
• /
• v.48 no.6
• /
• pp.552-559
• /
• 2011

An interest in Arbitrary Lagrangian Eulerian (ALE) finite element methods has been increased due to more accurate responses in Fluid-Structure Interaction(FSI) problems. The multi-material ALE approach was applied to the prediction of the acceleration response of free-fall lifeboat, and its responses were compared to those of the single-material ALE one. It could be found that even though there was no big difference in the simulation responses of two methods, the single-material and multi-material ALE ones, the latter multi-material ALE method showed a little bit more close response to those of experimental results compared to the former single-material ALE one, especially in the x- and z-direction acceleration responses. Through this study, it could be found that several parameters in the ALE algorithms have to be examined more carefully for a good structural safety assessment of FSI problems.

• Wind and Structures
• /
• v.17 no.3
• /
• pp.263-274
• /
• 2013

A numerical technique for fluid-structure interaction, which is based on the finite element method (FEM) and computational fluid dynamics (CFD), was developed for application to an industrial chimney equipped with a pendulum tuned mass damper (TMD). In order to solve the structural problem, a one-dimensional beam model (Navier-Bernoulli) was considered and, for the dynamical problem, the standard second-order Newmark method was used. Navier-Stokes equations for incompressible flow are solved in several horizontal planes to determine the pressure in the boundary of the corresponding cross-section of the chimney. Forces per unit length were obtained by integrating the pressure and are introduced in the structure using standard FEM interpolation techniques. For the fluid problem, a fractional step scheme based on a second order pressure splitting has been used. In each fluid plane, the displacements have been taken into account considering an Arbitrary Lagrangian Eulerian approach. The stabilization of convection and diffusion terms is achieved by means of quasi-static orthogonal subscales. For each period of time, the fluid problem was solved and the geometry of the mesh of each fluid plane is updated according to the structure displacements. Using this technique, along-wind and across-wind effects have been properly explained. The method was applied to an industrial chimney in three scenarios (with or without TMD and for different damping values) and for two wind speeds, showing different responses.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.5
• /
• pp.323-331
• /
• 2019

In this study, a computational model approach for the modeling of hydrodynamic pressures acting on radial gates during strong earthquakes is proposed. The use of the dynamic layering method with the Arbitrary Lagrangian Eulerian (ALE) algorithm and the SIMPLE method for simulating free reservoir surface flow in addition to moving boundary interfaces between the fluid domain and a structure due to earthquake excitation are suggested. The verification and validation of the proposed approach are realized by comparisons performed using the renowned formulation derived by the experimental results for vertical and inclined dam surfaces subjected to earthquake excitation. A parameter study for the truncated lengths of the two-dimensional fluid domain demonstrates that twice the water level leads to efficient and converged computational results. Finally, numerical simulations for large radial gates with different curvatures subjected to two strong earthquakes are successfully performed using the suggested computational model.