• 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.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.16 no.1
• /
• pp.118-123
• /
• 2017

It is well known that water jetting is now widely used in the advanced cutting processes of polymers, metals, glass, ceramics, and composite materials because of some advantages, such as heatless and non-contacting cutting different from the laser beam machining. In this paper, we proposed the simulation model of waterjet by lengths and the inner spiral structure of the nozzle. The simulation results show that the outlet velocity of the nozzle is faster than the inlet. Furthermore, we found rapid velocity reduction after passing through the outlet. The nozzle of diameter ${\phi}500$ and length 70mm, shows the optimal fluid width and velocity distribution. Also, the nozzle with inner spiral structure shows a Gaussian distribution of velocity and this model is almost twice as fast as the model without spiral structure, within the effective standoff distance (2.5 mm). In the future, when inserting abrasive material into the waterjet, we plan to analyze the fluid flow and the particle behavior through a simulation model.

• International Journal of Aerospace System Engineering
• /
• v.3 no.1
• /
• pp.30-38
• /
• 2016

The The turbulent mixing characteristics of Kerosene-LOx in a coaxial swirl injector 100 bar have been numerically investigated. Turbulent model is based on large eddy simulation with real-fluid transport and thermodynamics. The effects of equation of state (EOS), chamber pressure are evaluated in a point of the mixing efficiency and pressure fluctuations. The dominant frequency is same as the hairpin vortex shedding frequency generated by film wave at the LOx post.

• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.320-334
• /
• 2020

This study presents an application of hydro-mechanical coupled Particle Flow Code 3D (PFC3D) to simulation of fluid injection induced fault slip experiment conducted in Mont Terri Switzerland as a part of a task in an international research project DECOVALEX-2019. We also aimed as identifying the current limitations of the modelling method and issues for further development. A fluid flow algorithm was developed and implemented in a 3D pore-pipe network model in a 3D bonded particle assembly using PFC3D v5, and was applied to Mont Terri Step 2 minor fault activation experiment. The simulated results showed that the injected fluid migrates through the permeable fault zone and induces fault deformation, demonstrating a full hydro-mechanical coupled behavior. The simulated results were, however, partially matching with the field measurement. The simulated pressure build-up at the monitoring location showed linear and progressive increase, whereas the field measurement showed an abrupt increase associated with the fault slip We conclude that such difference between the modelling and the field test is due to the structure of the fault in the model which was represented as a combination of damage zone and core fractures. The modelled fault is likely larger in size than the real fault in Mont Terri site. Therefore, the modelled fault allows several path ways of fluid flow from the injection location to the pressure monitoring location, leading to smooth pressure build-up at the monitoring location while the injection pressure increases, and an early start of pressure decay even before the injection pressure reaches the maximum. We also conclude that the clay filling in the real fault could have acted as a fluid barrier which may have resulted in formation of fluid over-pressurization locally in the fault. Unlike the pressure result, the simulated fault deformations were matching with the field measurements. A better way of modelling a heterogeneous clay-filled fault structure with a narrow zone should be studied further to improve the applicability of the modelling method to fluid injection induced fault activation.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.5 no.2
• /
• pp.71-99
• /
• 2001

The model and analytical method for solving the problem of coupled fluid flow in the reservoir/well system is presented. The 3-D drainage area is composed of three connected media: the tubing, the annuli as a super conducting collector, and the reservoir itself. To couple these three types of fluid flows a non-overlapping Dirichlet-Neumann domain decomposition method is developed. The method allows us to apply an analytical hybrid simulator for accurate evaluation of the impact of main geometrical and hydrodynamic parameters of the 3-D system on the pressure drop along the horizontal well and its production index.

• 한국HCI학회:학술대회논문집
• /
• 2006.02a
• /
• pp.1120-1125
• /
• 2006

자연 현상에서 나타나는 물이나 바다와 같은 유체를 3 차원으로 시뮬레이션하는데 있어서 가장 중요한 요소는 실시간에 사실적으로 실행 가능하도록하는 것이다. 유체 모델은 특정 상황에 따른 다양한 방정식과 많은 파라미터값에 의해 제어되기 때문에 시뮬레이션하는데 많은 어려움이 따른다. 또한 복잡한 물리 수식을 기반으로 하기 때문에 유체 모델을 시뮬레이션하기 위해서는 많은 수행 시간이 소요된다. 본 논문에서는 실시간 유체와 강체(rigid body) 사이의 상호작용을 표현하기 위해 간략화된 유체 표면 모델(Fluid-Surface Model)을 제안하고, 개선된 계산과정을 통해 보다 빠르게 시뮬레이션하도록 한다. 또한 본 논문에서는 유체의 표면과 강체의 상호작용을 표현하는데 있어서 유체의 항력에 의해서 강체와 충돌시 발생하는 유체 표면의 움직임을 강체 모델의 제어를 통해 나타낸다. 본 논문에서 제안하는 자연스러운 유체 표면 모델은 유체역학적 방법을 사용하여 실시간에 사실적으로 표현된다. 그리고 이러한 유체 표면 모델을 PC 환경에서 사용자와 상호작용 가능하도록 재현하여, 게임이나 애니메이션에서의 유체 모델들에도 적용할 수 있다.

• International Journal of Fluid Machinery and Systems
• /
• v.6 no.3
• /
• pp.121-136
• /
• 2013

This paper reports the 1-D fluid transient simulation results of the discharge flow conditions in a 6-cylinder reciprocating slurry pump. Two discharge manifold configurations are studied comparatively; a case with a hexagon shaped discharge manifold where each cylinder discharges at a single vertex, and a case where all the cylinders discharges are lumped together into a tank shaped manifold. In addition, the study examines the effect of two pulsation mitigation measures in the case of hexagonal manifold; a single inline orifice in one of the hexagon sides and a volumetric dampener at the manifold outlet. The study establishes the pressure and flow fluctuation characteristics of each configuration and decouples the pulsation characteristics of the pump and the discharge manifold.

• Journal of the Korean Vacuum Society
• /
• v.18 no.3
• /
• pp.151-163
• /
• 2009

The domestic development status of Inductively Coupled Plasma (ICP) simulator which is based on fluid model is explained. As each part which composes the unified simulator, electron heating module, charged and neutral particle transport module, surface reaction module including a sheath model, and GUI (Graphic User Interface) with pre- and post-processors are described in order. Also, we present data base status of chemical reaction and physical collision, which has been applied to the recently developed simulator until now. Lastly, some future plans of development are suggested.

• Journal of Advanced Marine Engineering and Technology
• /
• v.27 no.3
• /
• pp.365-372
• /
• 2003

The Flow patterns around two cylinders in various arrangements were studied by a discrete vortex method. The flow for the surface of each cylinder was represented by arranging bound vortices at adequate intervals. The viscous diffusion of fluid was represented by the random walk method. The vortex distributions. streaklines. timelines and velocity vectors around two cylinders were calculated for centre-to-centre pitch ratios of P/D=1.5 and 2.5, attack ang1es of $a=0^{circ}, 30^{circ}, 60^{circ} and 90^{\circ}$. and Reynolds number of Re= 1200. The results of simulation correspond to the photographs by flow visualization and the flow intereference between two cylinders in various arrangements were clearly visualized by a numerical simulation.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.05a
• /
• pp.264-267
• /
• 2010

This research focuses on the development of numerical code to deal with compressible two phase flow around three dimensional objects combined with cavitation model suggested by Weishyy et al. with k-e turbulent model. The cryogenic cavitation is carried out by considering the thermodynamic effect on physical properties of cryogenic fluids in physical point of view and implementing the temperature sensitivity in the energy equation of the government equations in numerical point of view, respectively. The formulation has been extensively validated for both liquid nitrogen and liquid hydrogen by simulating the experiments of Hord on hydrofoils. Then, simulations of cavitating turbopump inducers at their design flow rate are presented. Results over a broad range of Nss numbers extending from single-phase flow conditions through the critical head break down point are discussed. In particular, thermal depression effects arising from cavitation in cryogenic fluids are identified and their impact on the suction performance of the inducer quantified.