• Journal of Korean Port Research
• /
• v.6 no.2
• /
• pp.25-31
• /
• 1992

Herein we investigate the mass transport velocity caused by the viscosity near the ocean structure such as circular pile and inclined breakwaters. The mass transport velocity which is represented by the sum of the Eulerian velocity and the stokes drift were derived by Carter, Liu and Mei(1973). The tangential components of the inviscid velocity field at the bottom needed in the calculation of the mass transport velocity is obtained by solving the scattering problem due to breakwaters. The matched asymptotic expansion technique is employed to obtain the inviscid flow fields scattered by inclined breakwaters. The numerical results show that heary sediments tends to be deposited near the center of breakwaters and that the narrowing of the entrance width results in reduction of the magnitude of mass transport.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.559-573
• /
• 2008

In this study, we newly proposed 3-D nonlinear wave driver utilizing the Navier-Stokes Eq. the numerical integration of which is carried out using SPH (Smoothed Particle Hydrodynamics), an internal wave generation with the source function of Gaussian distribution and an energy absorbing layer. For the verification of new 3-D nonlinear wave driver, we numerically simulate the sloshing problem within a parabolic water basin triggered by a Gaussian hump and uniformly inclined water surface by Thacker (1981). It turns out that the qualitative behavior of sloshing caused by relaxing the external force which makes a free surface convex or uniformly inclined is successfully simulated even though phase error is visible and an inundation height shrinks as numerical simulation more proceeds. For the more severe test, we also simulate the nonlinear shoaling and refraction over uniform beach of wedge shape. It is shown that numerically simulated waves are less refracted than the linear counterpart by Hamiltonian ray theory due to nonlinearity, energy dissipation at the bottom and side walls, energy loss induced by breaking, and the hydraulic jump occurring when breaking waves encounter a down-rush by the preceding wave.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.9
• /
• pp.597-604
• /
• 2016

A parallel algorithm of bi-conjugate gradient method was developed based on CUDA for parallel computation of the incompressible Navier-Stokes equations. The governing equations were discretized using splitting P2P1 finite element method. Asymmetric stenotic flow problem was solved to validate the proposed algorithm, and then the parallel performance of the GPU was examined by measuring the elapsed times. Further, the GPU performance for sparse matrix-vector multiplication was also investigated with a matrix of fluid-structure interaction problem. A kernel was generated to simultaneously compute the inner product of each row of sparse matrix and a vector. In addition, the kernel was optimized to improve the performance by using both parallel reduction and memory coalescing. In the kernel construction, the effect of warp on the parallel performance of the present CUDA was also examined. The present GPU computation was more than 7 times faster than the single CPU by double precision.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.29 no.5
• /
• pp.217-227
• /
• 2017

The object of this study is to estimate the net volume transport and the residual flow that changed by space and time at southern part of Yeomha channel, Gyeonggi Bay. The cross-section observation was conducted at the mid-part (Line2) and the southern end (Line1) of Yeomha channel for 13 hours during neap and spring-tides, respectively. The Lagrange flux is calculated as the sum of Eulerian flux and Stokes drift, and the residual flow is calculated by using least square method. It is necessary to unify the spatial area of the observed cross-section and average time during the tidal cycle. In order to unify the cross-sectional area containing such a large vertical tidal variation, it was necessary to convert into sigma coordinate system by horizontally and vertically for every hour. The converted sigma coordinate system is estimated to be 3~5% error when compared with the z-level coordinate system which shows that there is no problem for analyzing the data. As a result, the cross-sectional residual flow shows a southward flow pattern in both spring and neap tides at Line2, and also have characteristic of the spatial residual flow fluctuation: it northwards in the main line direction and southwards at the end of both side of the waterway. It was confirmed that the residual flow characteristics at Line2 were changed by the net pressure due to the sea level difference. The analysis of the net volume transport showed that it tends to southwards at $576m^3s^{-1}$, $67m^3s^{-1}$ in each spring tide and neap tide at Line2. On the other hand, in the control Line1, it has tendency to northwards at $359m^3s^{-1}$ and $248m^3s^{-1}$. Based on the difference between the two observation lines, it is estimated that net volume transport will be out flow about $935m^3s^{-1}$ at spring tide stage and about $315m^3s^{-1}$ at neap tide stage as the intertidal zone between Yeongjong Island and Ganghwa Island. In other words, the difference of pressure gradient and Stokes drift during spring and neap tide is main causes of variation for residual current and net volume transport.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.9
• /
• pp.1209-1217
• /
• 2002

The use of a two dimensional hierarchical elements are investigated for the incompressible flow computation. The construction of hierarchical elements are explained by both a geometric configuration and a determination of degrees of freedom. Also a systematic treatment of essential boundary values has been developed for the degrees of freedom corresponding to higher order terms. The numerical study for the poisson problem shows that the computation with hierarchical higher order elements can increase the convergence rate and accuracy of finite element solutions in more efficient manner than the use of standard first order element. for Stokes and Cavity flow cases, a mixed version of penalty function approach has been introduced in connection with the hierarchical elements. Solutions from hierarchical elements showed better resolutions with consistent trends in both mesh shapes and the order of elements.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.3
• /
• pp.933-944
• /
• 1996

In this paper, two new techniques, the pattern filling and the refined flow field regeneration, based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible viscous flow with free surfaces. The gorerning equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The penalty and Newton-Raphson methods are used effectively for finite element formulation. The flow front surface and the volume inflow rate are calculated using the pattern filling technique to select an adequate pattern among five filling patterns at each quadrilateral control volume. By the refined flow field regeneration technique, the new flow field which renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. Using the new thchniques to be developed, the dam-breaking problem has been analyzed to predict flow phenomenon of fluid and the predicted front positions versus time have been compared with the reported experimental result.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.239-244
• /
• 2009

The gridless (or meshfree) methods, such as MPS, SPH, FPM an so forth, are feasible and robust for the problems with moving boundary and/or complicated boundary shapes, because these methods do not need to generate a grid system. In this study, a gridless solver, which is based on the combination of moving least square interpolations on a cloud of points with point collocation for evaluating the derivatives of governing equations, is presented for two-dimensional unsteady incompressible Navier-Stokes problem in the low Reynolds number. A MAC-type algorithm was adopted and the Poission equation for the pressure was solved by successively in the moving least square sense. Some weighing functions were tested in order to investigate the up-winding effect for the convection term. Some typical problems were solved by the presented solver for the validation and the results obtained were compared with analytic solutions and the numerical results by conventional CFD methods, such as FVM.

• International Journal of Fluid Machinery and Systems
• /
• v.3 no.2
• /
• pp.150-159
• /
• 2010

Laidback fan shaped film-cooling hole is formulated numerically and optimized with the help of three-dimensional numerical analysis, surrogate methods, and the multi-objective evolutionary algorithm. As Pareto optimal front produces a set of optimal solutions, the trends of objective functions with design variables are predicted by hybrid multi-objective evolutionary algorithm. The problem is defined by four geometric design variables, the injection angle of the hole, the lateral expansion angle of the diffuser, the forward expansion angle of the hole, and the ratio of the length to the diameter of the hole, to maximize the film-cooling effectiveness compromising with the aerodynamic loss. The objective function values are numerically evaluated through Reynolds- averaged Navier-Stokes analysis at the designs that are selected through the Latin hypercube sampling method. Using these numerical simulation results, the Response Surface Approximation model are constructed for each objective function and a hybrid multi-objective evolutionary algorithm is applied to obtain the Pareto optimal front. The clustered points from Pareto optimal front were evaluated by flow analysis. These designs give enhanced objective function values in comparison with the experimental designs.

• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.19-26
• /
• 1999

This paper describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows with free-surface. The Navier-Stokes equations governing the flows have been discretized by means of finite-difference approximations, and the resulting equations have been solved via the SIMPLE-C algorithm. The free-surface is defined by the motion of a set of marker particles and the interface behaviour was investigated by means of a "Lagrangian" technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• Journal of computational fluids engineering
• /
• v.19 no.4
• /
• pp.20-28
• /
• 2014

An efficient solution algorithm for simulating free surface problem is presented. Navier-Stokes equations for variable density incompressible flow are employed as the governing equation on Cartesian meshes. In order to describe the free surface motion efficiently, VOF(Volume Of Fluid) method utilizing THINC(Tangent of Hyperbola for Interface Capturing) scheme is employed. The most time-consuming part of the current free surface flow simulations is the solution step of the linear system, derived by the pressure Poisson equation. To solve a pressure Poisson equation efficiently, the PCG(Preconditioned Conjugate Gradient) method is utilized. This study showed that the proper application of the preconditioner is the key for the efficient solution of the free surface flow when its pressure Poisson equation is solved by the CG method. To demonstrate the efficiency of the current approach, we compared the convergence histories of different algorithms for solving the pressure Poisson equation.