• 한국해안해양공학회지
• /
• 제11권2호
• /
• pp.95-106
• /
• 1999

마찰의 효과가 고려된 유한 진폭 천수방정식을 이용하여 잠제 주변의 파랑장, 잠제에 의하여 발생되는 반사파와 투과파 뿐만 아니라 시간과 공간에 따른 자유수면 및 흐름의 동수역학적인 거동 특성을 해석할 수 있는 Lax-Wendroff 유한차분 수치모형이 수립되었다. 비선형 규칙파와 고립파가 입사하는 경우에 대한 해석 결과를 기존의 수치해 및 실험자료와 비교하여 수치모형이 만족스럽게 검증되었다. 규칙파에 대한 해석에서 잠제 전면에서 발생되는 반사파와 관련된 시간과 공간에 따른 자유수면의 변화나 흐름특성이 투과파의 그것보다 강하게 발생되는 것으로 나타나, 본 연구에서 해석된 제한적인 조건하에서는 잠제가 외해에서 내습하는 파랑에너지를 효율적으로 차단하는 것으로 밝혀졌다. 그러나 고립파의경우는 대부분의 파랑에너지가 잠제에 의하여 차단되지 않고 투과되는 것으로 나타났다. 마지막으로 파랑장의 동수역학적인 거동특성으로부터 잠제를 피복하고 있는 피복재의 안전성과 관련된 해석이 간접적으로 수행되었다.

• 대한조선학회논문집
• /
• 제35권1호
• /
• pp.1-14
• /
• 1998

2차원 천수방정식에 기초한 조류계산 수치모델을 개발하였다. 복잡한 해안선의 입력을 단순화할 수 있는 직각격자를 사용하였으며, 방정식은 효율이 뛰어난 것으로 알려져 있는 인수분해음해법(Implicit Factorization Method)에 의하여 수치적으로 풀었다. 대류항은 상류차분법으로 처리하여 수치소산이 자동적으로 들어가도록 하였으며, 따라서 수치불안정에 의한 발산을 억제할 수 있었다. 모델의 검증을 위하여 선형화된 천수방정식을 풀고 계산결과를 해석해와 비교하였으며 잘 일치하고 있음을 확인하였다. 수치모델의 응용 예로서 인천항부근의 경인운하 개발예정인 수로에서의 조류 계산결과를 보였다.

• 한국수자원학회논문집
• /
• 제43권9호
• /
• pp.823-833
• /
• 2010

침수식생 개수로 흐름은 식생영역과 상부영역에 서로 다른 흐름구조를 보인다. 즉, 식생영역에서 전단으로 인해 생성되는 난류는 억제되며 비교적 균일한 유속 분포를 보이며 상부영역에서는 일반 개수로 흐름과 유사한 흐름구조를 보인다. 이와 같이 두 상이한흐름구조가 결합된 복잡한 흐름특성으로 인해 침수식생 개수로흐름은 공학적인 관심의 대상이 되어왔다. 본 연구에서는 침수식생 개수로 흐름의 층적분 모형의 비교 분석을 수행하였다. 일반적으로 식생흐름의 층적분 모형은 층의 수에 따라 2층 및 3층모형으로 구분한다. 즉, 전체 수심을 식생영역과 상부영역으로 구분하는 2층모형과 식생영역을 바닥 조도의 영향 유무에 따라 내부 및 외부 식생영역으로 구분하는 3층모형으로 분류된다. 본 연구에서는 2층모형과 3층모형을 비교하였다. 다양한 실험조건에 적용한 결과, 3층모형이 식생영역에서 유속의 변화를 고려할 수 있으나 결과는 레이놀즈응력 분포에 민감하며, 적분된 유속은 2층모형에 의한 예측 결과가 더욱 정확한 것으로 나타났다. 3층모형에서 내부 식생영역의 결과가 전체 흐름구조에 미치는 영향이 무시할 수 있으므로 이 점을 착안하여 식생영역에서 유속 변화가 고려되는 수정 2층모형을 제시하였다. 수정 2층모형에서 가정하는 레이놀즈응력 분포는 상부영역에서는 선형, 식생영역에서는 멱함수 형으로 변화한다. 다양한 조건에 적용한 결과, 수정 2층모형이 대체로 기존의 모형과 비슷한 정도의 예측을 수행하나 식생밀도가 매우 작은 흐름의 경우 예측 결과가 불량한 것으로 나타났다.

• 대기
• /
• 제23권1호
• /
• pp.73-83
• /
• 2013

In recent several years, East Asia, Europe and North America have suffered successive cold winters and a number of historical records on the extreme weathers are replaced with new record-breaking cold events. As a possible explanation, several studies suggested that cryospheric conditions of Northern Hemisphere (NH), i.e. Arctic sea-ice and snow cover over northern part of major continents, are changing significantly and now play an active role for modulating midlatitude atmospheric circulation patterns that could bring cold winters for some regions in midlatitude. In this study, a dynamical seasonal prediction system for NH winter is newly developed using the snow depth initialization technique and statistically predicted sea-ice boundary condition. Since the snow depth shows largest variability in October, entire period of October has been utilized as a training period for the land surface initialization and model land surface during the period is continuously forced by the observed daily atmospheric conditions and snow depths. A simple persistent anomaly decaying toward an averaged sea-ice condition has been used for the statistical prediction of sea-ice boundary conditions. The constructed dynamical prediction system has been tested for winter 2012/13 starting at November 1 using 16 different initial conditions and the results are discussed. Implications and a future direction for further development are also described.

• 한국수자원학회논문집
• /
• 제43권9호
• /
• pp.801-811
• /
• 2010

급경사에서의 고립파의 처오름을 예측하기 위해 3차원 수치모형에 내부조파기법을 도입하여 수치모형실험을 수행하였다. 수치모형은 Navier-Stokes 방정식을 유한차분법을 이용하여 계산하는 동수압 모형으로서, 난류의 해석을 위해서 상대적으로 큰 에디(eddy)만을 고려하는 SANS(spatially averaged Navier-Stokes) 방정식을 푸는 LES(large-eddy-simulation) 기반의 수치모형을 사용한다. 엇갈림 격자체계에서 유한차분법을 사용하여 지배방정식을 해석하는 모형으로서 수치기법으로 Two-step projection 기법을 사용하여 SANS 방정식을 풀었으며, Poisson 방정식을Bi-CGSTAB 기법을 이용하여 풀고 압력장을 계산하였다. 또한, 자유수면의 추적을 위하여 2차 정확도의 VOF(volume-of-fluid) 기법을 사용하였다. 먼저 고립파를 3차원 공간의 일정 수심상에서 내부조파하여 해석해와 비교한 후 분산오차에 대해 분석하였다. 그리고 고립파를 내부조파하여 급경사에서의 고립파의 처오름 및 처내림 현상을 예측하고 수리모형 실험결과와 비교 및 분석하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2006년도 PARALLEL CFD 2006
• /
• pp.311-314
• /
• 2006

The supersonic flow around tandem cavities was investigated by three- dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes(RANS) equation with the $\kappa-\omega$ thrbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split using van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge-Kutta method. The aspect ratio of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two-dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the fire cavity flow cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• Nuclear Engineering and Technology
• /
• 제47권6호
• /
• pp.700-711
• /
• 2015

Sodium-cooled fast breeder reactors use liquid sodium as a moderator and coolant to transfer heat from the reactor core. The main hazard associated with sodium is its rapid reaction with water. Sodium-water reaction (SWR) takes place when water or vapor leak into the sodium side through a crack on a heat-transfer tube in a steam generator. If the SWR continues for some time, the SWR will damage the surface of the defective area, causing it to enlarge. This self-enlargement of the crack is called "self-wastage phenomena." A stepwise numerical evaluation model of the self-wastage phenomena was devised using a computational code of multicomponent multiphase flow involving a sodium-water chemical reaction: sodiumwater reaction analysis physics of interdisciplinary multiphase flow (SERAPHIM). The temperature of gas mixture and the concentration of NaOH at the surface of the tube wall are obtained by a numerical calculation using SERAPHIM. Averaged thermophysical properties are used to assess the local wastage depth at the tube surface. By reflecting the wastage depth to the computational grid, the self-wastage phenomena are evaluated. A two-dimensional benchmark analysis of an SWAT (Sodium-Water reAction Test rig) experiment is carried out to evaluate the feasibility of the numerical model. Numerical results show that the geometry and scale of enlarged cracks show good agreement with the experimental result. Enlarged cracks appear to taper inward to a significantly smaller opening on the inside of the tube wall. The enlarged outer diameter of the crack is 4.72 mm, which shows good agreement with the experimental data (4.96 mm).

• Journal of Mechanical Science and Technology
• /
• 제20권8호
• /
• pp.1256-1265
• /
• 2006

The supersonic flows around tandem cavities were investigated by two-dimensional and three-dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes (RANS) equation with the k- ω turbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split with van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge- Kutta method. The aspect ratios of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two- dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the first cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• 한국해안·해양공학회논문집
• /
• 제26권2호
• /
• pp.96-102
• /
• 2014

하구역의 지형변화를 예측하기 위해 3차원 해빈 변형 모델을 통한 계산이 수행되었다. 본 모델은 수심적분을 기초로 한 준 3차원 연안흐름 모듈로 구성되며, 해안선의 변화, 부유사의 이송-확산 효과를 고려할 수 있다. 우선 모델의 성능을 확인하기 위해 3차원 해빈 변화 모델이 하구역 사주 형성에 적용되었다. 다음으로 동해에 인접한 Ara 강 하구에 모델이 적용되었다. Ara 강의 사주의 동계 변화가 재현되었으며 계산결과는 현장 관측 결과와 좋은 일치를 나타냈다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2019년도 학술발표회
• /
• pp.267-267
• /
• 2019

Driftwood is one of serious problems in a river environment. In several countries, such as Indonesia, Japan, and Italy, the driftwood frequently appears in a river basin, and it can alter the channel bed, flow configuration by wood deposition and jam formation. Therefore, the studies related to driftwood have been actively conducted by many researchers to understand the mechanism of driftwood dynamics. In particular, wood motion by collision is one of the difficult issues in the numerical simulation because the calculation for wood collision requires significantly expensive calculation time due to small time step. Thus, this study conducted the numerical simulation in consideration of the wood motion by water flow and wood collision to understand the wood dynamics in terms of computation. We used the 2D (two-dimensional) depth-averaged velocity model, Nays2DH, which is a Eulerian model to calculate the water flow on the generalized coordinate. A Lagrangian type driftwood model, which expresses the driftwood as connected sphere shape particles, was employed to Nays2DH. In addition, the present study considered root wad effect by using larger diameter for a particle at a head of driftwood. An anisotropic bed friction was considered for the sliding motion dependent on stemwise, streamwise and motion directions. We particularly considered changeable draft at each particle and projection area by an angle between stemwise and flow directions to precisely reproduce the wood motions. The simulation results were compared with experimental results to verify the model. As a result, the simulation results showed good agreement with experimental results. Through this study, it would be expected that this model is a useful tool to predict the driftwood effect in the river flow.