• Journal of Ship and Ocean Technology
• /
• 제4권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.

• 한국해안해양공학회지
• /
• 제6권2호
• /
• pp.196-204
• /
• 1994

도서가 산재해 있고 수심의 변화가 불규칙한 해역을 공간적으로 이산화하기 위하여 전처리 과정에서 가변 격자를 자동으로 발생하여 새만금 해역에 적용함으로써 방조제 축조의 영향을 해석하였다. 비선형의 주기적 유한요소모형(TEANL)이 해석의 도구로 이용되었는데 실측유속과 비교적 잘 일치되는 양호한 결과를 얻을 수 있었다. 모형 적용결과 건설중민 새만금 방조제의 신시도와 가력도 그리고 대항리, 변산반도를 잇는 방조제 전면 해역은 폭이 넓은 일종의 구형 폐쇄만으로 변화되어 이 일대 해역의 순환 및 확산 환경에 심각한 영향을 미칠 수 있을 것으로 예측되었다.

• 한국농공학회지
• /
• 제36권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.

• 한국안전학회지
• /
• 제15권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.

• 방사성폐기물학회지
• /
• 제4권1호
• /
• pp.9-15
• /
• 2006

결정질 암반에 위치한 가상의 방사성폐기물처분장에 대한 정확한 안전성 평가를 수행하기 위해서는 다공암반으로 대표되는 공학방벽 및 결정질 암반으로 대표되는 자연방벽을 모두 고려한 매질에서의 물질 거동 특성을 정확하게 모사하는 것이 중요하다. 본 논문에서는 단열-다공암반내 유동 및 물질전달의 정확한 해석을 위한 연결망 구축방법에 대하여 서술하였다. 실제 단열암반을 사실적으로 모사하기 위하여 단열 물성자료는 확률밀도함수를 이용하여 생성하였다. 2차원의 사각격자로 모사된 단열암반과 육면체로 모사된 다공암반간 교차선의 원활한 검색을 위하여 단열암반과 연결전 다공암반의 교차면에 가상의 단열암반을 추가적으로 도입하였다. 전체유동경로를 구성하기 저하여 생성된 단열들 간의 교차선을 효율적으로 검색할 수 있는 방법 및 단열암반 및 단열-다공암반간의 연결도를 신속히 확인하는 알고리즘을 제안하였다. 이러한 방법들은 추후 방사성폐기물처분장에서의 단열-다공 암반을 통한 핵종 이동 특성을 모사할 수 있는 수치코드 개발에 많은 도움을 줄 것으로 사료된다.

• 대한토목학회논문집
• /
• 제34권5호
• /
• pp.1435-1442
• /
• 2014

본 연구에서는 OpenFOAM에서 제공하는 소스코드를 이용하여 매끄러운 하상의 직사각형 개수로 흐름에 대해 수치모의를 수행하였다. 난류 해석을 위해 LES를 수행하였는데, 비등방성 잔여 응력항을 모델링하기 위해서 Germano et al. (1991)이 제시한 Dynamic Subgrid-scale 모형을 이용하였다. 조직구조를 분석하기 위하여 Lu and Willmarth (1973)가 제시한 uw 사분면기법을 이용하여 순간레이놀즈 응력이 레이놀즈 응력에 미치는 영향을 기여율과 시간비로 나누어 분석하였다. LES 모의 결과를 토대로 기존 실험 및 DNS 모의 결과와 비교하고 분석하였다. 매끈한 하상을 가진 개수로 흐름에서 완충층 이후의 구간에서 분출현상이 쓸기현상에 비해 레이놀즈 응력의 양의 생성에 기여하는 바가 크지만, 분출현상에 비해 쓸기현상의 발생빈도가 큰 것으로 확인되었다.

• 한국전산유체공학회지
• /
• 제16권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.

• 한국항공우주학회지
• /
• 제44권4호
• /
• pp.290-297
• /
• 2016

사각덕트에서 난류 유동장으로 분사되는 액체 제트의 액주 분열과 미립화 현상에 관한 LES를 수행하였다. 기체상태의 공기 유동 해석에 오일러리안 해법을 사용하고, 액적 추적을 위하여 라그랑지안 해법을 사용하여 기체-액체간 이상유동(two phase flow) 해석을 수행하였다. 액적 분열 모델, 아격자 스케일 모델 및 공간 차분법에 따른 액적 분열을 조사하였다. 액체 제트의 침투깊이를 경험식과 비교하였으며 경험식보다 약간 높음을 알 수 있었다. 제트 후류에서 사우터 평균직경에 대한 분석을 수행하였다.

• 한국해양공학회지
• /
• 제33권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
• /
• 제34권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).