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

• Journal of Ocean Engineering and Technology
• /
• v.33 no.4
• /
• pp.313-321
• /
• 2019

The prediction of the hydrodynamic performance of a planing hull vessel is an important and challenging topic for computational fluid dynamic (CFD) applications to naval hydrodynamics. In this paper, the resistance and planing attitude analysis for a Fridsma hull, which is a prismatic planing hull, in still water are numerically studied using OpenFOAM. OpenFOAM is an open source code package based on C++ libraries and the finite volume method (FVM) for the discretization of the RANS equation. The volume of fluid method (VOF) is used to capture the water-air interface and the SST ${\kappa}-{\omega}$ model is used for the turbulence simulation. The overset mesh method is used to capture the large motion of the hull at higher speeds. Before the extensive analysis, uncertainty analyses using various time steps and grid sizes were performed for one ship speed case of Fn = 1.19. The results of the present study are compared with those of a model test, other CFD research, and Savitsky's empirical formula. The results of the present study, following the trend of other CFD results, slightly over predict the resistance and under predict the sinkage and, more significantly, the trim.

• Journal of Korea Water Resources Association
• /
• v.56 no.spc1
• /
• pp.1037-1048
• /
• 2023

Fishways are constructed to ensure the fish migration because river-crossing structures such as dams and weirs cut off the stream longitudinal connectivity and influence on aquatic ecosystems. However, the passage efficiency of fishes varies depending on flow characteristics in the fishway and fish species. In this study, three-dimensional numerical simulations are carried out using a RANS model and the volume of fluid method for resolving free surface fluctuations to calculate the turbulent flow in the pool and weir type fishway. The Flow velocity and turbulent kinetic energy in the pool of fishway are analyzed according to variation of the upstream water level and the length of pool. The present numerical simulations reasonably well reproduce the stream flow and plunging flow characteristics in the pool. The simulation results show that the stream flow changes to the plunging flow as the length of the pool increases. When the upstream level increases, the stream flow becomes more evident. Key parameters related to the fish migration within the fishway such as the flow velocity and the turbulent kinetic energy are examined to assess the ascending possibility of fishes.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.51 no.5
• /
• pp.1-8
• /
• 2009

Photobioreactor (PBR) that houses and cultivates microalgae providing a suitable environment for its growth, such as light, nutrients, CO2, heat, etc. is now getting more popular in the last decade. Among the many types of PBRs, the bubble column type is very attractive because of its simple construction and easy operation. However, despite the availability of these PBRs, only a few of them can be practically used for mass production. Many limitations still holdback their use especially during their scale-up. To enlarge the culture volume and productivity while supplying optimum environmental conditions, various PBR structures and process control are needed to be investigated. In this study, computational fluid dynamics (CFD) was economically used to design a bubble-column type PBR taking the place of field experiments. CFD is a promising technique which can simulate the growth and production of microalgae in the PBR. To study bubble column PBR with CFD, the most important factor is the possibility of realizing bubble. In this study, multi-phase models which are generally used to realize bubbles were compared by theoretical approaches and comparing in a 2D simulation. As a result, the VOF (volume of fluid) model was found to be the most effective model to realize the bubbles shape as well as the flow inside PBR which may be induced by bubble injection. Considering the accuracy and economical efficiency, 0.005 second time step size was chosen for 2.5 mm mesh size. These results will be used as criteria for scale-up in the PBR simulation.

• Tribology and Lubricants
• /
• v.33 no.1
• /
• pp.1-8
• /
• 2017

Drag torque reduction in a wet clutch pack is a key aspect of the design process of the dual clutch transmission (DCT) system. In order to reduce the drag torque caused by lubricant shear resistance, recently developed wet clutch pack systems of DCT, as well as automatic transmission and other four-wheel drive (4WD) couplings, frequently utilize wavy wet clutch pads. Therefore, wavy shape of friction pad are made on the groove patterns like waffle pattern for the reduction of drag torque. Additionally, the groove patterns are designed with larger channels at several locations on the friction pad to facilitate faster outflow of lubricant. However, channel performance is a function of the waviness of the friction pad at the location of the particular channel. This is because the discharge sectional area varies according to the waviness amplitude at the location of the particular channel. The higher location of the additional channel on the friction pad results in a larger cross-sectional area, which allows for a larger flow discharge rate. This results in reduction of the drag torque caused by the shear resistance of DCTF, because of marginal volume fraction of fluid (VOF) in the space between the friction pad and separator. This study computes the VOF in the space between the friction pad and separator, the hydrodynamic pressure developed, and the shear resistance of friction torque, by using CFD software (FLUENT). In addition, the study investigates the dependence of these parameters on the location and waviness amplitude of the channel pattern on the friction pad. The paper presents design guidelines on the proper location of high waviness amplitude on wavy friction pads.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.4
• /
• pp.425-429
• /
• 2011

In the unlikely of nuclear reactor meltdown, the leaked core melt or corium must be contained in a device called core-catcher so that the corium can be cooled and stabilized. The ex-vessel behavior of corium involves complex physical and chemical mechanisms of flow propagation, heat transfer, and reactions with sacrificial substrates. In this study, the detailed characteristics of corium flow and heat transfer were investigated by using a commercial CFD code for VULCANO VE-U7 test reported in the literature. The volume-of-fluid (VOF) model was used to predict the interfacial surface formation of corium and the surrounding air, and the discrete ordinate model was adopted to calculate radiation between corium and the surroundings. It was found that cooling via radiation through the top surface of corium had a dominant effect on the temperature and viscosity profiles at the front of the corium flow.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.12 no.8
• /
• pp.3374-3383
• /
• 2011

In this study, numerical analyses for slamming impact phenomena have been carried out using a 2-dimensional wedge shaped structure having finite deadrise angles. Fluid is assumed incompressible and entry speed of the structure is kept constant. Geo-reconstruct(or PLIC-VOF) scheme is used for the tracking of the deforming free surface. Numerical analyses are carried out for the deadrise angles of $10^{\circ}$, $20^{\circ}$ and $30^{\circ}$. For each deadrise angle, variations are made for the grid size on the wedge bottom and for the entry speed. The magnitude and the location of impact pressure and the total drag force, which is the summation of pressure distributed at the bottom of the structure, are analyzed. Results of the analyses are compared with the results of the Dobrovol'skaya similarity solutions, the asymptotic solution based on the Wagner method and the solution of Boundary Element Method(BEM).

• Journal of Korea Water Resources Association
• /
• v.37 no.6
• /
• pp.461-466
• /
• 2004

This study presents a combined experimental and numerical effort to investigate wave breaking of sinusoidal waves in a surf zone. Numerical predictions are verified by comparing to laboratory measurements. The model solves the Reynolds equations and$textsc{k}$-$\varepsilon$ models for the turbulence analysis. To track the free surface displacement, the volume of fluid method is employed. As the height of incident wave increases, the wave breaking occurs at a closer point of the slope in the numerical model and laboratory experiments with the same depth and period. When a wave breaking occurs, the ratio of wave height becomes larger, with the same wave height and depth, as the period increases.

• The KSFM Journal of Fluid Machinery
• /
• v.15 no.5
• /
• pp.60-66
• /
• 2012

One of commonly physical phenomena encountered in pump sump systems in which its significant influence to the hydraulic performance of pump system plays an important role in the field of fluid engineering, is the appearance of free surface and submerged vortices. In this paper, a study of the vortices behavior and their formative mechanism of asymmetry is considered in this paper by using numerical approach. The Reynolds-Averaged Navier-Stokes (RANS) equations and k-omega Shear Stress Transport turbulence model used to describe the properties of turbulent flows, in company with VOF multiphase model, are implemented by Fluent code with multi-block structured grid system. In the numerical simulation, the calculated elevation of air-water interface and vortex core contours are used to classify visually surface vortices as well as submerged vortices. It is shown that the free surface vortex is identified by the concavity of liquid region from the free surface and swirling flow at that own plane. To investigate the distinctive behavior of these vortices corresponding to each given flow rate at the same water level, some numerical testing of them are considered here in such a manner that the flow pattern of surface vortex are obtained similarly to the obtained results from experiment. Furthermore, the influence due to the change of grid refinement and the variation of depth of the concavity are also considered in this paper. From that, these influential factors will be implemented to design a good pump sump with higher performance in the future.

• Journal of the Semiconductor & Display Technology
• /
• v.14 no.3
• /
• pp.61-66
• /
• 2015

Recently, the major trends of NIL are high throughput and large area patterning. For UV NIL, if it can be proceeded in the non-vacuum environment, which greatly simplifies tool construction and greatly shorten process times. However, one key issue in non-vacuum environment is air bubble formation problem. In this paper, numerical analysis of bubble defect of UV NIL is performed. Fluent, flow analysis focused program was utilized and VOF (Volume of Fluid) skill was applied. For various resist-substrate and resist-mold angles, effects of velocity inlet and residual layer thickness of resist on bubble defect formation were investigated. The numerical analyses show that the increases of velocity inlet and residual layer thickness can cause the bubble defect formation, however the decreases of velocity inlet and residual layer thickness take no difference in the bubble defect formation.