• Journal of Ship and Ocean Technology
• /
• v.4 no.2
• /
• pp.24-37
• /
• 2000

The purpose of this paper is to study the upstream waves in front of an advancing arbitrary hull shape in a restricted water channel. Conventionally, in a restricted water channel, shallow water effects are amplified because of the finite water depth and width. When the effects of shallow water and the restricted channel width are severe, upstream waves propagate forward from the fore-body of the advancing hull. In this study, numerical simulations are carried out for the relevant analysis of the flow phenomena by the draft variation of advancing hull in a restricted water channel. Numerical simulations are done with a finite-difference method based on the MAC scheme in a rectangular grid system.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.2
• /
• pp.196-204
• /
• 1994

Studied was impact assessment of sea dike construction in Saemankeum. To represent complexity of the geometry and topography of the region a flexible finite grid system are adopted. Combined fine and coarse meshes based on automatic mesh generator were applied in pre-processing. A nonlinear periodic finite element model. TEANL, was implied in this analysis, which gave good results compared to the observed data. It was predicted that the front region of dike connecting Shinsi-Karyeok-Daehang-Pyunsan will behave as a closed rectangular bay with wide width, which may affect significantly to the circulation and dispersion mechanism in the region.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.36 no.3
• /
• pp.144-154
• /
• 1994

An efficient tidal model, TIDE which is an iterative type, nonlinear finite element model has developed for the analysis of the tidal movement in the coastal area which is characterized by irregular boundaries and bottom topography. Traditional time domain finite element models have been in difficulties with requirement for high eddy viscosity coefficients and small time steps to insure numerical instability. These problems are overcome by operating in the frequency domain with an elaborate grid system by combining the triangular and quadrilateral shape grids. Furthermore, in order to handle non-linearity which will be more significant in the shallow region, an iterative scheme with least square error minimization algorithm has been implemented in the model. The results of TIDE model are agreed with the analytical solutions in a rectangular channel under the condition of tidal waves entering the channel closed at one end.

• Journal of the Korean Society of Safety
• /
• v.15 no.3
• /
• pp.14-22
• /
• 2000

A numerical analysis has been performed on the two-dimensional rectangular gallium melting problem using the enthalpy method. The major advantage of this method is that the physical domain is discretized with fixed grids without transforming variables and the interface conditions of phase change are accounted for the definition of suitable source terms in the governing equations. But in the fixed method, there is some ambiguity in defining the porosity constant which has no physical interpretation. If the velocity correction is included in the momentum equation, for the appropriate range of porosity constant, the realistic predictions are obtained. The object of the present work is to predict the phase change within the molten steel with thin riser slab using the modified enthalpy-porosity method. The computational procedures for predicting velocity and temperature are based on the finite volume method and the non-staggered grid system. The influence of natural convection on the melting process is considered. A comparison with the experimental results shows that the modified method is better than the previous one.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.4 no.1
• /
• pp.9-15
• /
• 2006

For the accurate safety assessment of potential radioactive waste disposal site which is located in the crystalline rock it is important to simulate the mass transportation through engineered and natural barrier system precisely, characterized by porous and fractured media respectively. In this work the methods to construct discrete fracture network for the analysis of flow and mass transport through fractured-porous medium are described. The probability density function is adopted in generating fracture properties for the realistic representation of real fractured rock. In order to investigate the intersection between a porous and a fractured medium described by a 2 dimensional rectangular and a cuboid grid respectively, an additional imaginary fracture is adopted at the face of a porous medium intersected by a fracture. In order to construct large scale flow paths an effective method to find interconnected fractures and algorithms of swift detecting connectivities between fractures or porous medium and fractures are proposed. These methods are expected to contribute to the development of numerical program for the simulation of radioactive nuclide transport through fractured-porous medium from radioactive waste disposal site.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.5
• /
• pp.1435-1442
• /
• 2014

This study presented numerical simulations of smooth-bed flows in the rectangular open-channel using the source code by OpenFOAM. For the analysis of the turbulent flow, Large Eddy Simulations were carried out and the dynamic sub-grid scale model proposed by Germano et al. (1991) is used to model the residual stress term. In order to analyze the coherent structure, the uw quadrant method proposed by Lu and Willmarth (1973) is used and the contribution rate and the fraction time of the instantaneous Reynolds stress are obtained in the Reynolds stress. The results by the present study are analyzed and compared with data from previous laboratory studies and direct numerical simulations. It is found that the contribution rate of the ejection events is larger than that of sweep events over the buffer layer in the open-channel flow over the smooth bed, however, the frequency of the sweep event is higher than that of the ejection events.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.39-46
• /
• 2011

A numerical study of a turbulent natural convection in an enclosure with the lattice Boltzmann method (LBM) is presented. The primary emphasis of the present study is placed on investigation of accuracy and numerical stability of the LBM for the turbulent natural convection flow. A HYBRID method in which the thermal equation is solved by the conventional Reynolds averaged Navier-Stokes equation method while the conservation of mass and momentum equations are resolved by the LBM is employed in the present study. The elliptic-relaxation model is employed for the turbulence model and the turbulent heat fluxes are treated by the algebraic flux model. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with the deferred correction way to ensure accuracy and stability of solutions. The present LBM is applied to the prediction of a turbulent natural convection in a rectangular cavity and the computed results are compared with the experimental data commonly used for the validation of turbulence models and those by the conventional finite-volume method. It is shown that the LBM with the present HYBRID thermal model predicts the mean velocity components and turbulent quantities which are as good as those by the conventional finite-volume method. It is also found that the accuracy and stability of the solution is significantly affected by the treatment of the convection term, especially near the wall.

• 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 Ocean Engineering and Technology
• /
• v.33 no.2
• /
• pp.131-138
• /
• 2019

This study focused on the roll motion characteristics of a two-dimensional (2D) rectangular floating structure under regular beam sea conditions. An experiment was conducted in a 2D wave tank for a roll free decay test in calm water and the roll motion in a range of regular waves with and without heave motion to investigate the motion response and heave influence on the roll motion. A numerical study was carried out using Reynolds-averaged Navier Stokes (RANS)-based CFD simulations. A grid convergence test was conducted to accurately capture the wave condition on the free surface based on the overset mesh and wave forcing method. It was found in the roll free decay test that the numerical results agreed well with the experimental results for the natural roll period and roll damping coefficient. It was also observed that the heave motion had an impact on the roll motion, and the responses of the heave and roll motion from the CFD simulations were in reasonable agreement with those from the experiment.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.59-72
• /
• 2022

In this work the numerical results of the flow around a 5:1 rectangular cylinder at Reynolds numbers 3 000 and 40 000, zero angle of attack and smooth incoming flow condition are presented. Implicit Large Eddy Simulations (ILES) have been performed with a high-order accurate spatial scheme and an implicit high-order accurate time integration method. The spatial approximation is based on a discontinuous Galerkin (dG) method, while the time integration exploits a linearly-implicit Rosenbrock-type Runge-Kutta scheme. The aim of this work is to show the feasibility of high-fidelity flow simulations with a moderate number of DOFs and large time step sizes. Moreover, the effect of different parameters, i.e., dimension of the computational domain, mesh type, grid resolution, boundary conditions, time step size and polynomial approximation, on the results accuracy is investigated. Our best dG result at Re=3 000 perfectly agrees with a reference DNS obtained using Nek5000 and about 40 times more degrees of freedom. The Re=40 000 computations, which are strongly under-resolved, show a reasonable correspondence with the experimental data of Mannini et al. (2017) and the LES of Zhang and Xu (2020).