• Journal of Ocean Engineering and Technology
• /
• v.25 no.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.

• Korea-Australia Rheology Journal
• /
• v.19 no.4
• /
• pp.233-242
• /
• 2007

We performed a numerical investigation to find out the optimal choice of the spatial discretization in the distributed-Lagrangian-multiplier/fictitious-domain (DLM/FD) method for the solid/fluid interaction problem. The elastic solid bar attached on the bottom in a pressure-driven channel flow of a Newtonian fluid was selected as a model problem. Our formulation is based on the scheme of Yu (2005) for the interaction between flexible bodies and fluid. A fixed regular rectangular discretization was applied for the description of solid and fluid domain by using the fictitious domain concept. The hydrodynamic interaction between solid and fluid was treated implicitly by the distributed Lagrangian multiplier method. Considering a simplified problem of the Stokes flow and the linearized elasticity, two numerical factors were investigated to clarify their effects and to find the optimum condition: the distribution of Lagrangian multipliers and the solid/fluid interfacial condition. The robustness of this method was verified through the mesh convergence and a pseudo-time step test. We found that the fluid stress in a fictitious solid domain can be neglected and that the Lagrangian multipliers are better to be applied on the entire solid domain. These results will be used to extend our study to systems of elastic particle in the Stokes flow, and of particles in the viscoelastic fluid.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.6
• /
• pp.347-355
• /
• 2016

This study presents a Lattice Boltzmann Method (LBM) coupled with a momentum-exchange approach/fictitious domain (MEA/FD) method for the simulation of particle suspensions. The method combines the advantages of the LB and the FD methods by using two unrelated meshes, namely, a Eulerian mesh for the flow domain and a Lagrangian mesh for the solid domain. The rigid body conditions are enforced by the momentum-exchange scheme in which the desired value of velocity is imposed directly in the particle inner domain by introducing a pseudo body force to satisfy the constraint of rigid body motion, which is the key idea of a fictitious domain (FD) method. The LB-MEA/FD method has been validated by simulating two different cases, and the results have been compared with those through other methods. The numerical evidence illustrated the capability and robustness of the present method for simulating particle suspensions.

• Korea-Australia Rheology Journal
• /
• v.21 no.3
• /
• pp.193-200
• /
• 2009

We present a direct simulation technique for two-dimensional mold-filling simulations of fluids filled with a large number of circular disk-like rigid particles. It is a direct simulation in that the hydrodynamic interaction between particles and fluid is fully considered. We employ a pseudo-concentration method for the evolution of the flow front and the DLM (distributed Lagrangian multipliers)-like fictitious domain method for the implicit treatment of the hydrodynamic interaction. Both methods allow the use of a fixed regular discretization during the entire computation. The discontinuous Galerkin method has been used to solve the concentration evolution equation and the rigid-ring description has been introduced for freely suspended particles. A buffer zone, the gate region of a finite area subject to the uniform velocity profile, has been introduced to put discrete particles into the computational domain avoiding any artificial discontinuity. From example problems of 450 particles, we investigated the particle motion and effects of particles on the flow for both Newtonian and shear-thinning fluid media. We report the prolonged particle movement toward the wall in case of a shear-thinning fluid, which has been interpreted with the shear rate distribution.

• The Journal of the Acoustical Society of Korea
• /
• v.31 no.1
• /
• pp.19-28
• /
• 2012

The time-domain solution from the Kirchhoff integral equation for an exterior problem is not unique at certain eigen-frequencies associated with the fictitious internal modes as happening in frequency-domain analysis. One of the solution methods is the CHIEF (Combined Helmholtz Integral Equation Formulation) approach, which is based on employing additional zero-pressure constraints at some interior points inside the body. Although this method has been widely used in frequency-domain boundary element method due to its simplicity, it was not used in time-domain analysis. In this work, the CHIEF approach is formulated appropriately for time-domain acoustic boundary element method by constraining the unknown surface pressure distribution at the current time, which was obtained by setting the pressure at the interior point to be zero considering the shortest retarded time between boundary nodes and interior point. Sound radiation of a pulsating sphere was used as a test example. By applying the CHIEF method, the low-order fictitious modes could be damped down satisfactorily, thus solving the non-uniqueness problem. However, it was observed that the instability due to high-order fictitious modes, which were beyond the effective frequency, was increased.

• Geophysics and Geophysical Exploration
• /
• v.20 no.3
• /
• pp.121-128
• /
• 2017

Interpretation of time-domain electromagnetic (TEM) data has been made mostly based on one-dimensional (1-D) inversion scheme in Korea. A proper interpretation of TEM data should employ 3-D TEM forward and inverse modeling algorithms. This study developed a 3-D TEM modeling algorithm using a finite difference time-domain (FDTD) method with staggered grid. In numerically solving Maxwell equations, fictitious displacement current is included based on an explicit FDTD method using a central difference approximation scheme. The developed modeling algorithm simulated a small-coil source configuration to be verified against analytic solutions for homogeneous half-space models. Further, TEM responses for a 3-D anomaly are modeled and analyzed. We expect that it will contribute greatly to the precise interpretation of TEM data.

• International Journal of Aeronautical and Space Sciences
• /
• v.4 no.1
• /
• pp.53-62
• /
• 2003

In this study, nonlinear aeroelastic characteristics of an supersonic missile wing with strong shock interferences are investigated. The missile wing model has a freeplay structural nonlinearity at its pitch axis. To practically consider the effects of freeplay structural nonlinearity, the fictitious mass method is applied to structural vibration analysis based on finite element method. Nonlinear aerodynamic flows with unsteady shock waves are also considered in supersonic flow regions. To solve the nonlinear aeroelastic governing equations including the freeplay effect, a modal-based coupled time-marching technique based on the fictitious mass method is used in the time-domain. Various aeroelastic computations have been performed for the nonlinear wing structure model. Linear and nonlinear aeroelastic analyses have been conducted and compared with each other in supersonic flow regions. Typical nonlinear limit cycle oscillations and phase plots are presented to show the complex vibration phenomena with simultaneous fluid-structure nonlinearities.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2002.11b
• /
• pp.914-920
• /
• 2002

Nonlinear flow-induced vibration characteristics of a generic missile wing (or control surface) are investigated in this study. The wing model has freeplay structural nonlinearity at its pitch axis. Nonlinear aerodynamic flows with unsteady shock waves are considered in the transonic flow region. To practically consider the effects of freeplay structural nonlinearity, the fictitious mass method (FMM) is applied to structural vibration analysis based on a finite element method (FEM). A computational fluid dynamics (CFD) technique is used for computing the nonlinear unsteady aerodynamics of all-movable wings. The aerodynamic analysis is based on the efficient transonic small-disturbance aerodynamic equations of motion using the potential-flow theory. To solve the nonlinear aeroelastic governing equations including the freeplay effect, a modal-based computational structural dynamic (CSD) analysis technique based on fictitious mass method (FMM) is used in time-domain. In addition, CSD and unsteady CFD techniques are simultaneously coupled to give accurate computational results. Various aeroelastic computations have been performed for a generic missile wing model. Linear and nonlinear aeroelastic computations have been conducted and the characteristics of flow-induced vibration are introduced.

• Journal of the Society of Naval Architects of Korea
• /
• v.49 no.3
• /
• pp.213-221
• /
• 2012

The present study numerically investigates the interaction between a free-surface and flow around a circular cylinder over a moving wall. In order to simulate the flow past the circular cylinder over a moving wall near a free-surface, this study has adopted the direct-forcing/fictitious domain (DF/FD) method with the level set method in the Cartesian coordinates. Numerical simulation is performed for a Reynolds numbers of 100 in the range of $0.25{\leq}g/D{\leq}2.00$ and $0.5{\leq}h/D{\leq}2.00$, where g/D and h/D are the gaps between the cylinder and a moving wall and the cylinder and a free-surface normalized by cylinder diameter D, respectively. According to g/D and h/D, the vortex structures have been classified into three patterns of the two-row, one-row, steady elongation. In general, both of g/D and h/D have the large values which mean the cylinder is far away from the wall and the free-surface, two-row vortex structure forms in the wake. When g/D decreases, the two-row vortex structure gradually transfers into the one-row vortex structure. When the g/D reveals the critical value below which the flow becomes steady state, resulting in the steady elongation vortex.

• Tunnel and Underground Space
• /
• v.30 no.1
• /
• pp.15-28
• /
• 2020

In this study, a procedure was introduced to estimate strength and deformation modulus of the 3-D discrete fracture network(DFN) systems using the distinct element method(DEM). Fracture entities were treated as non-persistent square planes in the DFN systems. Systematically generated fictitious fractures having similar mechanical characteristics of intact rock were combined with non-persistent real fractures to create polyhedral blocks in the analysis domain. Strength and deformation modulus for 10 m cube domain of various deterministic and stochastic 3-D DFN systems were estimated using the DEM to explore the applicability of suggested method and to examine the effect of fracture geometry on strength and deformability of DFN systems. The suggested procedures were found to effective in estimating anisotropic strength and deformability of the 3-D DFN systems.