• Tribology and Lubricants
• /
• v.38 no.3
• /
• pp.109-114
• /
• 2022

This paper reports the development of an oil flushing system combined with a microbubble generator. Oil flushing plays a crucial role in regulating the lubricant's performance during the lubricant replacement process. Moreover, harmful contaminants, such as sludge, wear particles, and rust, from piping systems or lubrication system can be removed by oil flushing. Oil flushing aims to increase the system's efficiency using a dedicated flushing oil, increasing of the supply pressure and generating a vortex. In addition, it helps the mechanical system or equipment achieve peak performance and reduces the potential for premature failure. However, the contaminant-removal applications of existing oil flushing system are limited. In this research, we aim to improve the performance of oil flushing system by incorporating a microbubble generator, which uses the venture effect to generate microbubbles and mixes them with lubricant. The microbubbles in the blended lubricant remove contaminants from the lubrication system more effectively. Structural mechanics and fluid dynamics are analyzed through fluid-structure interaction (FSI) analysis, and the numerical analysis results are used for the designing the system. The magnitude of the maximum stress is investigated based on the pressure results obtained by the CFD analysis; through the CFD analysis, the mixing ratio of air (bubble) and lubricant is evaluated using the volume of fluid (VOF) model according to the working conditions.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.3
• /
• pp.258-263
• /
• 2008

Performance of injector equiped with screw type swirler which is suitable for supersonic cruise engine combustor was investigated using theoretical, numerical, and experimental methods. Based on discharge coefficient and spray angle which represent the performance of injectors, the geometrical parameters which affect these performance parameters were defined, control group injectors were designed, and variation of performance parameters according to the geometrical parameters were examined. Within the defined range, measured value of performance of injectors was smaller than result of theoretical prediction, and prediction result from numerical simulation using VOF method agreed with the result of experiments very well. The viscous barrier was not observed, and minimum discharge coefficient and maximum spray angle, 0.05 and 104 respectively, was obtained for this type of injector.

• Nuclear Engineering and Technology
• /
• v.56 no.6
• /
• pp.2029-2038
• /
• 2024

A deep understanding of the characteristics and mechanism of geyser boiling and capillary pumping is necessary to optimize a high-temperature sodium heat pipe. In this work, the Volume of Fluid (VOF) two-phase model and the capillary force model in the mesh wick were used to model the complex phase change and fluid flow in the heat pipe. Computational Fluid Dynamics (CFD) simulations successfully predicted the process of bubble nucleation, growth, aggregation, and detachment from the wall in the liquid pool of the evaporation section of the heat pipe in horizontal and tilted states, as well as the reflux phenomenon of capillary suction within the wick. The accuracy and stability of the capillary force model within the wick were verified. In addition, the causes of geyser boiling in heat pipes were analyzed by extracting the oscillation distribution of heat pipe wall temperature. The results show that adding the capillary force model within the wick structure can reasonably simulate the liquid backflow phenomenon at the condensation; Under the horizontal and inclined operating conditions of the heat pipe, the phenomenon of local dry-out will occur, resulting in a sharp increase in local temperature. The speed of bubble detachment and the timely reflux of liquid sodium (condensate) replenishment in the wick play a vital role in the geyser temperature oscillation of the tube wall. The numerical simulation method and the results of this study are anticipated to provide a good reference for the investigation of geyser boiling in high-temperature heat pipes.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.1
• /
• pp.58-66
• /
• 2014

Numerical simulations for the flows containing free surface remain difficult problems because the drastic differences of physical properties of water and air, The difference of densities makes the solution instable in particular. For the stabilities of the solutions, the most typical methods to simulate free surface flows, such as Volume Of Fluid(VOF) and Level-Set(LS) methods, impose transient zones where the physical prosperities are continuously distributed. The thickness of the transient zone is the source of the numerical errors. The other side, marker-density method does not use such a transient zone. In the traditional marker-density method, however, the air velocities of free surface cells are extrapolated from the water velocity, and the pressures on the free surface are extrapolated from the air pressures for the stability of the solution. In this study, the marker-density method is modified for the decrease of such numerical errors. That is, the pressure on the free surface is determined to coincide with the pressure gradient terms of the governing equations, and the velocity of free surface cells are calculated with the governing equations. Two-dimensional steady spilling breakers behind of a submersed hydrofoil and three-dimensional spilling breaker near a wedge shaped ship model are simulated using INHAWAVE-II including the modified marker-density(MMD) method. The results are compared with the results of Fluent V6.3 including VOF method and several published research results.

• Journal of the Society of Naval Architects of Korea
• /
• v.56 no.4
• /
• pp.361-372
• /
• 2019

The flow pattern of air layers and skin-friction drag reduction by air injection are investigated to find the suitable multiphase flow model using unstructured finite-volume CFD solver for the Reynolds-averaged Navier-Stokes equations. In the present computations, two different multiphase flow modeling approaches, such as the Volume of Fluid (VOF) and the Eulerian Multi-Phase (EMP), are adopted to investigate their performances in resolving the two-phase flow pattern and in estimating the frictional drag reduction. First of all, the formation pattern of air layers generated by air injection through a circular opening on the bottom of a flat plate are investigated. These results are then compared with those of MMkiharju's experimental results. Subsequently, the quantitative ratios of skin-friction drag reduction including the behavior of air layers, within turbulent boundary layers in large scale and at high Reynolds number conditions, are investigated under the same conditions as the model test that has been conducted in the US Navy's William B. Morgan Large Cavitation Channel (LCC). From these results, it is found that both VOF and EMP models have similar capability and accuracy in capturing the topology of ventilated air cavities so called'air pockets and branches'. However, EMP model is more favorable in predicting quantitatively the percentage of frictional drag reduction by air injection.

• Journal of Korea Water Resources Association
• /
• v.48 no.4
• /
• pp.281-290
• /
• 2015

Numerical simulations of free surface flow over a broad-crested rectangular weir are conducted by using the volume of fraction (VOF) method and three different turbulence models, the k-${\varepsilon}$, RNG k-${\omega}$ and k-${\omega}$ SST models. The governing equations are solved by a second-order accurate finite volume method and the grid sensitivity study of solutions is carried out. The numerical results are evaluated by comparing the solutions with experimental and numerical results of Kirkgoz et al. (2008) and some non-dimensionalized experimental results obtained by Moss (1972) and Zachoval et al. (2012). The results show that the present numerical model can reasonably reproduce the experimental results, while three turbulent models yield different numerical predictions of two distinct zones of flow separation, the first zone is in front of the upstream edge of the weir and the second is created immediately behind the upstream edge of the weir where the flow is separated to form the separation bubble. The standard k-${\varepsilon}$ model appears to significantly underestimate the size of both separation zones and the k-${\omega}$ SST model slightly over-estimates the first separation zone in front of the weir. The RNG k-${\varepsilon}$ model predicts both separation zones in overall good agreement with the experimental measurement, while the k-${\omega}$ SST model yields the best numerical prediction of separation bubble at the upstream edge of the weir.

• Journal of Ocean Engineering and Technology
• /
• v.26 no.6
• /
• pp.7-13
• /
• 2012

This paper presents a numerical study of the wave loads acting on offshore structures using a Cartesian-grid-based flow simulation method. Finite volume discretization with a volume-of-fluid (VOF) method is adopted to solve two-phase Navier-Stokes equations. Among the many variations of the VOF method, the CICSAM scheme is applied. The body boundary conditions are satisfied using a porosity function, and wave generation is carried out by using transient (wave or damping) zone approaches. In order to validate the present numerical method, three different basic offshore structures, including a sphere, Pinkster barge, and Wigley model, are numerically investigated. First, diffraction and radiation problems are solved using the present numerical method. The wave exciting and drift forces from the diffraction problems are compared with potential-based solutions. The added mass and wave damping forces from the radiation problems are also compared with the potential results. Next, the wave-induced motion responses of the structures are calculated and compared with the existing experimental data. The comparison results are fairly good, showing the validity of the present numerical method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.10
• /
• pp.979-988
• /
• 2012

The 3D structure of GDL for fuel cells was measured using high-resolution X-ray tomography in order to study material transport in the GDL. A computational algorithm has been developed to remove noise in the 3D image and construct 3D elements representing carbon fibers of GDL, which were used for both structural and fluid analyses. Changes in the pore structure of GDL under various compression levels were calculated, and the corresponding volume meshes were generated to evaluate the anisotropic permeability of gas within GDL as a function of compression. Furthermore, the transfer of liquid water and reactant gases was simulated by using the volume of fluid (VOF) and pore-network model (PNM) techniques. In addition, the simulation results of liquid water transport in GDL were validated by analogous experiments to visualize the diffusion of fluid in porous media. Through this research, a procedure for simulating the material transport in deformed GDL has been developed; this will help in optimizing the clamping force of fuel cell stacks as well as in determining the design parameters of GDL, such as thickness and porosity.

• Ocean Systems Engineering
• /
• v.7 no.4
• /
• pp.435-460
• /
• 2017

The main objective of this study is to investigate the turning and zig-zag maneuvering performance of the well-known naval surface combatant DTMB (David Taylor Model Basin) 5415 hull with URANS (Unsteady Reynolds-averaged Navier-Stokes) method. Numerical simulations of static drift tests have been performed by a commercial RANS solver based on a finite volume method (FVM) in an unsteady manner. The fluid flow is considered as 3-D, incompressible and fully turbulent. Hydrodynamic analyses have been carried out for a fixed Froude number 0.28. During the analyses, the free surface effects have been taken into account using VOF (Volume of Fluid) method and the hull is considered as fixed. First, the code has been validated with the available experimental data in literature. After validation, static drift, static rudder and drift and rudder tests have been simulated. The forces and moments acting on the hull have been computed with URANS approach. Numerical results have been applied to determine the hydrodynamic maneuvering coefficients, such as, velocity terms and rudder terms. The acceleration, angular velocity and cross-coupled terms have been taken from the available experimental data. A computer program has been developed to apply a fast maneuvering simulation technique. Abkowitz's non-linear mathematical model has been used to calculate the forces and moment acting on the hull during the maneuvering motion. Euler method on the other hand has been applied to solve the simultaneous differential equations. Turning and zig-zag maneuvering simulations have been carried out and the maneuvering characteristics have been determined and the numerical simulation results have been compared with the available data in literature. In addition, viscous effects have been investigated using Eulerian approach for several static drift cases.

• Journal of the Semiconductor & Display Technology
• /
• v.9 no.3
• /
• pp.29-34
• /
• 2010

Nanoimprint lithography (NIL) is an emerging technology enabling cost effective and high throughput nanofabrication. To successfully imprint a nanometer scale patterns, the understanding of the mechanism in nanoimprint forming is essential. In this paper, a numerical analysis of polymer flow in thermal NIL was performed. First, a finite element model of the periodic mold structure with prescribed boundary conditions was established. Then, the volume of fluid (VOF) and grid deformation method were utilized to calculate the free surfaces of the polymer flow based on an Eulerian grid system. From the simulation, the velocity fields and the imprinting pressure for constant imprinting velocity in thermal NIL were obtained. The velocity field is significant because it can directly describe the mode of the polymer deformation, which is the key role to determine the mechanism of nanoimprint forming. Effects of different mold shapes and various thicknesses of polymer resist were also investigated.