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

댐붕괴파 (dam-break flow)나 지진해일에 의해 발생하는 undular bore와 충격파 (shock) 현상을 동수압 및 분산효과를 고려하여 수치모의를 수행하였다. 완전비선형 Boussinesq-type equations 모형을 이용하여, 동수압 및 분산 효과를 고려하였다. 방정식은 4차 정확도의 유한체적법을 이용하여 해석하였고, 시간적으로도 4차정확도의 기법을 이용하여 고차미분항에 대한 수치분산을 억제하였다. 다양한 경우의 1차원과 2차원 공간에서의 수치모의를 수행하고 검증을 수행하였다. 그 결과, 완전비선형 Boussinesq-type equations 모형은 천수방정식 (shallow water equations) 기반의 모형에서 재현이 불가능한 undular bore 등을 재현 하는 등, 전반적으로 천수방정식 기반의 모형 보다 물리적으로도 타당하고 정량적으로도 실험결과와 잘 일치하는 경향을 보였다. 즉, 댐붕괴파나 지진해일 등에 의한 범람 모의에 있어 동수압과 분산 효과의 중요성이 공학적으로도 매우 중요한 고려사항 임이 나타났다.

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

본 연구에서는 ${\sigma}$-좌표계를 기반으로 하는 3차원의 이송확산 모형과 depth-integrated eddy simulation 모형을 결합한 효율적인 3차원 근역 (near-field) 해석모형을 제시하였다. 흐름 모형은 Boussinesq-type equations과 stochastic backscatter model을 기본으로 하고 있다. 이 흐름 모형은 수면의 변화와 바닥으로부터 발생하는 전단력과 파랑의 유동으로부터 발행하는 수심방향의 유속분포를 예측할 수 있다. 이와 같은 흐름 정보를 3차원 ${\sigma}$-좌표계의 이송확산모형에 제공하고 scalar의 이송과 확산에 대한 거동을 계산한다. 기본적인 이송과 이송-확산에 대한 검증 및 개수로에서 정량적 검증과 정성적 검증을 수행하였다. 전반적으로 타당한 결과가 도출되어 모형의 적합성이 있음을 확인하였다.

• 한국수자원학회논문집
• /
• 제46권5호
• /
• pp.491-504
• /
• 2013

본 연구에서는 shelf 지형을 대상으로 불규칙파를 적용한 2차원 쇄파변형 수리실험과 Boussinesq형 방정식을 이용한 수치해석을 수행하였다. 유의파고, 유의파봉고, 유의파곡고, 평균수위 및 안정파고와 같은 쇄파특성을 수리실험자료로부터 분석하였고, 실험결과와 Boussinesq형 방정식의 해석결과를 비교하였으며, 두 결과는 비교적 잘 일치하였다. 일정수심상에서의 불규칙파의 쇄파후 안정파고는 수심의 약 0.56배로 분석되었다.

• 한국해안해양공학회지
• /
• 제18권2호
• /
• pp.166-177
• /
• 2006

본 연구에서는 Boussinesq방정식을 지배방정식으로 한 기존의 강비선형 모형과 고차 모형(다층모형)의 정확도를 각각 검토하였으며, 사용된 다층모형은 Lynett과 Liu(2004a)에 의해 개발된 것이다. 수치모형으로 shelf 지형에서의 쇄파를 모의하였으며, 실험에 적용된 파랑의 비선형성(${k_0}{A_0}$)은 0.029~0.180 범위이다. 전반적으로 2층 모형의 해석결과가 실험결과와 잘 일치하였다. 1층모형의 경우에는 쇄파점 가까운 곳에서부터 파고가 빨리 증폭이 된 반면, 2층모형의 파고증폭 정도는 수리실험결과와 잘 부합하였다.

• 한국항만학회지
• /
• 제12권1호
• /
• pp.75-86
• /
• 1998

To design a coastal structure in the nearshore region, engineers must have means to estimate wave climate. Waves, approaching the surf zone from offshore, experience changes caused by combined effects of bathymetric variations, interference of man-made structure, and nonlinear interactions among wave trains. This paper has attempted to find out the effects of two of the more subtle phenomena involving nonlinear shallow water waves, amplitude dispersion and secondary wave generation. Boussinesq-type equations can be used to model the nonlinear transformation of surface waves in shallow water due to effect of shoaling, refraction, diffraction, and reflection. In this paper, generalized Boussinesq equations under the complex bottom condition is derived using the depth averaged velocity with the series expansion of the velocity potential as a product of powers of the depth of flow. A time stepping finite difference method is used to solve the derived equation. Numerical results are compared to hydraulic model results. The result with the non-linear dispersive wave equation can describe an interesting transformation a sinusoidal wave to one with a cnoidal aspect of a rapid degradation into modulated high frequency waves and transient secondary waves in an intermediate region. The amplitude dispersion of the primary wave crest results in a convex wave front after passing through the shoal and the secondary waves generated by the shoal diffracted in a radial manner into surrounding waters.

• 한국대기환경학회지
• /
• 제16권5호
• /
• pp.511-520
• /
• 2000

The urban atmosphere is characterized by th difference in surface and atmospheric environment between urban and more natural area. To investigate th climatic effect of land use type and anthropogenic heat of urban on wind field, numerical simulations were carried out under typical summer synoptic condition. The wind model PNU_MCM(Pusan National University Mesoscale Circulation Model) is based on the three-dimensional Boussinesq equations, taking into account the hydrostatic assumption . Since lane-use differs over every subdivision on Pusan the surface energy budget model includes sub0grid parameterization scheme which can calculate the total heat flux over a grid surface composed of different surfaces. The simulated surface wind agrees well with the observed value, and average over 6 days which represent typical summer lan-sea breeze days, August 1998, i.e. negligible gradient winds and almost clear skies. Urbanization makes sea-breeze enhance at day and reduce land-breeze at night. The results show that contribution of land-use type is much larger than that of anthropogenic heat in Pusan.

• Ocean Systems Engineering
• /
• 제10권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.

• Ocean Systems Engineering
• /
• 제11권1호
• /
• pp.83-97
• /
• 2021

The submerged U-shape breakwater interaction with the solitary wave is simulated by the Boussinesq equations using the finite-difference scheme. The wave reflection, transmission, and dissipation (RTD) coefficients are used to investigate the U-shape breakwater's performance for different crest width, Lc1, and indent breakwater height, du. The results show that the submerged breakwater performance for a set of U-shape breakwater with the same cross-section area is related to the length of submerged breakwater crest, Lc1, and the distance between the crests, Lc2 (or the height of du). The breakwater has the maximum performance when the crest length is larger, and at the same time, the distance between them increases. Changing the Lc1 and du of the U-shape breakwaters result in a significant change in the RTD coefficients. Comparison of the U-shape breakwater, having the best performance, with the averaged RTD values shows that the transmission coefficients, Kt, has a better performance of up to 4% in comparison to other breakwaters. Also, the reflection coefficients KR and the diffusion coefficients, Kd shows a better performance of about 30% and 55% on average, respectively. However, the model governing equations are non-dissipative. The non-energy conserving of the transmission and reflection coefficients due to wave and breakwater interaction results in dissipation type contribution. The U-shape breakwater with the best performance is compared with the rectangular breakwater with the same cross-section area to investigate the economic advantages of the U-shape breakwater. The transmission coefficients, Kt, of the U-shape breakwater shows a better performance of 5% higher than the rectangular one. The reflection coefficient, KR, is 60% lower for U-shape in comparison to rectangular one; however, the diffusion coefficients, Kd, of U-shape breakwater is 35% higher than the rectangular breakwater. Therefore, we could say that the U-shape breakwater has a better performance than the rectangular one.

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

흐름과 장파에 의하여 발생되는 난류의 subgrid scale mixing effects를 고려할 수 있는 수심적분형 모형(depth-integrated model)을 제시하였다. 완전비선형의 수심적분형 모형은 약분산(weakly dispersive) 환경에서 흐름의 회전성(rotational)을 고려하도록 perturbation approach를 이용하여 유도되었다. 동일한 방법을 이용하여 수심적분형 이송확산방정식(depth-integrated scalar transport equation)을 유도하였다. 방정식은 4차정확도의 유한체적기법을 이용하여 해석하였으며, 다양한 혼합양상을 보이는 흐름에 대한 수치모의를 수행하였다.

• Nuclear Engineering and Technology
• /
• 제36권6호
• /
• pp.559-570
• /
• 2004

A computational fluid dynamics model for predicting moderator circulation inside the Canada deuterium uranium (CANDU) reactor vessel has been developed to estimate the local subcooling of the moderator in the vicinity of the calandria tubes. The buoyancy effect induced by the internal heating is accounted for by the Boussinesq approximation. The standard $k-{\varepsilon}$ turbulence model with logarithmic wall treatment is applied to predict the turbulent jet flows from the inlet nozzles. The matrix of the calandria tubes in the core region is simplified to a porous media in which the anisotropic hydraulic impedance is modeled using an empirical correlation of pressure loss. The governing equations are solved by DFX-4.4, a commercial CFD code developed by AEA technology. The resultant flow patterns of the constant-z slices containing the inlet nozzles and the outlet port are "mined-type", as observed in the former 2-dimensional experimental investigations. With 103% full power for conservatism, the maximum temperature of the moderator is $82.9^{\circ}C$ at the top of the core region. Considering the hydrostatic pressure change, the minimum subcooling is $24.8^{\circ}C$.