• Journal of the Korean Society of Hazard Mitigation
• /
• v.8 no.4
• /
• pp.91-100
• /
• 2008

Numerical implementation with a Cartesian cut-cell method is conducted in this study. A Cartesian cut-cell method is an easy and efficient mesh generation methodology for complex geometries. In this method, a background Cartesian grid is employed for most of computational domain and a cut-cell grid is applied for the peculiar grids where the flow characteristics are changed such as solid boundary to enhance the accuracy, applicability and efficiency. Accurate representation of complex geometries can be obtained by using the cut-cell method. The cut-cell grids are constructed with irregular meshes which have various shape and size. Therefore, the finite volume method is applied to numerical discretization on a irregular domain. The HLLC approximate Riemann solver, a Godunov-type finite volume method, is employed to discretize the advection terms in the governing equations. The weighted average flux method applied on the Cartesian cut cell grid for stabilization of the numerical results. To validate the numerical model using the Cartesian cut-cell grids, the model is applied to the rectangular tank problem of which the exact solutions exist. As a comparison of numerical results with the analytical solutions, the numerical scheme well represents flow characteristics such as free surface elevation and velocities in x-and y-directions in a rectangular tank with the Cartesian and cut-cell grids.

• Ocean Systems Engineering
• /
• v.10 no.2
• /
• pp.201-226
• /
• 2020

A systematic numerical comparative study of the performance of semicircular and rectangular submerged breakwaters interacting with solitary waves is the basis of this paper. To accomplish this task, Nwogu's extended Boussinesq model equations are employed to simulate the interaction of the wave with breakwaters. The finite difference technique has been used to discretize the spatial terms while a fourth-order predictor-corrector method is employed for time discretization in our numerical model. The proposed computational scheme uses a staggered-grid system where the first-order spatial derivatives have been discretized with fourth-order accuracy. For validation purposes, five test cases are considered and numerical results have been successfully compared with the existing analytical and experimental results. The performances of the rectangular and semicircular breakwaters have been examined in terms of the wave reflection, transmission, and dissipation coefficients (RTD coefficients) denoted by K_R, K_T, K_D. The latter coefficient K_D emerges due to the non-energy conserving K_R and K_T. Our computational results and graphical illustrations show that the rectangular breakwater has higher reflection coefficients than semicircular breakwater for a fixed crest height, but as the wave height increases, the two reflection coefficients approach each other. un the other hand, the rectangular breakwater has larger dissipation coefficients compared to that of the semicircular breakwater and the difference between them increases as the height of the crest increases. However, the transmission coefficient for the semicircular breakwater is greater than that of the rectangular breakwater and the difference in their transmission coefficients increases with the crest height. Quantitatively, for rectangular breakwaters the reflection coefficients K_R are 5-15% higher while the diffusion coefficients K_D are 3-23% higher than that for the semicircular breakwaters, respectively. The transmission coefficients K_T for rectangular breakwater shows the better performance up to 2.47% than that for the semicircular breakwaters. Based on our computational results, one may conclude that the rectangular breakwater has a better overall performance than the semicircular breakwater. Although the model equations are non-dissipative, the non-energy conserving transmission and reflection coefficients due to wave-breakwater interactions lead to dissipation type contribution.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.10
• /
• pp.961-969
• /
• 2013

In this study, the numerical analysis of a turbulent flow around both a staggered and an inline tube bundle was conducted using ANSYS CFX V.13, a commercial CFD software. The flow was assumed to be steady, incompressible, and isothermal. According to the CFD Best Practice Guideline, the sensitivity study for grid size, accuracy of the discretization scheme for convection term, and turbulence model was conducted, and its result was compared with the experimental data to estimate the applicability of the CFD Best Practice Guideline. It was concluded that the CFD Best Practice Guideline did not always guarantee an improvement in the prediction performance of the commercial CFD software in the field of tube bundle flow.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.4
• /
• pp.290-297
• /
• 2016

A two-phase Large Eddy Simulation(LES) has been conducted to investigate breakup and atomization of a liquid jet in a cross turbulent flow in a rectangular duct. Gas-droplet two-phase flow was solved by a coupled Eulerian-Lagrangian method which tracks every individual particles. Effects of liquid breakup models, sub-grid scale models, and a order of spatial discretization was investigated. The penetration depth in cross flow was comparable with experimental data by varying breakup model and LES scheme. SMD(Sauter Mean Diameter) distribution downstream of jet was analyzed.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.1
• /
• pp.94-101
• /
• 1994

This paper describes a short-term prediction model for the movement of an oil slick in shallow waters. Under the assumption that the initial movement of the oil slick is governed by spreading and advection, the model has been developed and applied to Kyungki-Bay near Incheon Harbor. The initial spreading is estimated by using an empirical formula. The depth-averaged momentum equations are solved numerically for the volume transport velocities, in which the $M_2$ tide is the main driving source. A staggered grid system is adopted fur spatial discretization and the half-time method is implemented for time marching. The numerical result is visualized with the help of animation and thus the contaminated area is displayed on a monitor in time sequence. The input data are the time, the location and the volume of spill accident as well as environmental data such as md and $M_2$ tide.

• 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 the Computational Structural Engineering Institute of Korea
• /
• v.25 no.6
• /
• pp.559-567
• /
• 2012

A level set based topological shape optimization method for nonlinear structure considering hyper-elastic problems is developed. To relieve significant convergence difficulty in topology optimization of nonlinear structure due to inaccurate tangent stiffness which comes from material penalization of whole domain, explicit boundary for exact tangent stiffness is used by taking advantage of level set function for arbitrary boundary shape. For given arbitrary boundary which is represented by level set function, a Delaunay triangulation scheme is used for current structure discretization instead of using implicit fixed grid. The required velocity field in the actual domain to update the level set equation is determined from the descent direction of Lagrangian derived from optimality conditions. The velocity field outside the actual domain is determined through a velocity extension scheme based on the method suggested by Adalsteinsson and Sethian(1999). The topological derivatives are incorporated into the level set based framework to enable to create holes whenever and wherever necessary during the optimization.

• Journal of Korea Water Resources Association
• /
• v.52 no.5
• /
• pp.361-372
• /
• 2019

The purpose of this study is to develop a two-dimensional land surface flood analysis model based on uniform square grid using the governing equations except for the convective acceleration term in the momentum equation. Finite volume method and implicit method were applied to spatial and temporal discretization. In order to reduce the execution time of the model, parallel computation techniques using CPU were applied. To verify the developed model, the model was compared with the analytical solution and the behavior of the model was evaluated through numerical experiments in the virtual domain. In addition, inundation analyzes were performed at different spatial resolutions for the domestic Janghowon area and the Sebou river area in Morocco, and the results were compared with the analysis results using the CAESER-LISFLOOD (CLF) model. In model verification, simulation results were well matched with the analytical solution, and the flow analyses in the virtual domain were also evaluated to be reasonable. The results of inundation simulations in the Janghowon and the Sebou river area by this study and CLF model were similar with each other and for Janghowon area, the simulation result was also similar to the flooding area of flood hazard map. The different parts in the simulation results of this study and the CLF model were compared and evaluated for each case. The results of this study suggest that the model proposed in this study can simulate the flooding well in the floodplain. However, in case of flood analysis using the model presented in this study, the characteristics and limitations of the model by domain composition method, governing equation and numerical method should be fully considered.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.12
• /
• pp.1175-1183
• /
• 2013

Even if some CFD software developers and its users think that a state-of-the-art CFD software can be used to reasonably solve at least single-phase nuclear reactor safety problems, there remain limitations and uncertainties in the calculation result. From a regulatory perspective, the Korea Institute of Nuclear Safety (KINS) is presently conducting the performance assessment of commercial CFD software for nuclear reactor safety problems. In this study, to examine the prediction performance of commercial CFD software with the porous model in the analysis of the scale-down APR (Advanced Power Reactor Plus) internal flow, a simulation was conducted with the on-board numerical models in ANSYS CFX R.14 and FLUENT R.14. It was concluded that depending on the CFD software, the internal flow distribution of the scale-down APR was locally somewhat different. Although there was a limitation in estimating the prediction performance of the commercial CFD software owing to the limited amount of measured data, CFX R.14 showed more reasonable prediction results in comparison with FLUENT R.14. Meanwhile, owing to the difference in discretization methodology, FLUENT R.14 required more computational memory than CFX R.14 for the same grid system. Therefore, the CFD software suitable to the available computational resource should be selected for massively parallel computations.