• 한국전산유체공학회:학술대회논문집
• /
• 2005.04a
• /
• pp.75-78
• /
• 2005

A numerical investigation has been performed to discuss the radiation-affected steady-laminar natural convection in an enclosure under a large temperature difference. Due to inherent nature of this study, the Boussinesq approximation is no longer valid. Therefore the radiating fluid in an enclosure is treated as a ideal gas. To examine the effects of thermal radiation on thermo-fluid dynamic behaviors in complex geometries, two incomplete partitions are introduced. Based on the results of this study, the dispositions of incomplete partitions with radiatively participating medium are found to incur a distinct difference in fluid-dynamic as well as thermal behavior.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2019.05a
• /
• pp.98-98
• /
• 2019

중력류 또는 밀도류는 주변 유체에 비해 상대적으로 밀도가 큰 유체가 밀도차에 의한 추진력으로 흐르는 것이다. 중력류의 수치모델링에는 두 가지 어려움이 있다. 즉, 적합한 지배방정식을 구성하여 적용하는 것 그리고 난류의 영향을 합리적으로 반영하는 것이다. 기존 중력류 해석을 위한 지배방정식들은 유체의 연속방정식과 운동량 방정식 그리고 밀도 또는 농도의 이송방정식을 조합하여 구성된다. 이들 지배방정식을 이용한 연구들은 대부분 두 유체 사이의 밀도차가 충분히 작아서 밀도 변동(variations)의 영향은 오로지 부력항에서만 유지된다는 Boussinesq 근사에 근거를 둔다. 그리고 이송방정식에서 밀도 또는 농도의 확산계수을 점성계수의 함수로 표현하기 위해서 Schmidt 수를 이용한다. 수치모델링에서 Schimdt 수는 상수값을 적용하지만, 이 값은 밀도의 연직방향 경사에 근거한 부력빈도(buoyancy frequency)와 난류량의 따라 큰 차이를 보이는 것으로 알려져있다. 한편, 표준 통계학적 난류모델과 벽함수를 적용한 수치모델링은 초기 중력에 의해서 무너지는(slumping) 단계를 넘어 관성력으로 추진되는 단계와 점성 효과가 지배적인 단계에서는 정확도에 현저히 낮아지기 때문에 대부분 큰와모의(large-eddy simulation, LES) 또는 DNS(direct numerical simulation)수준의 고해상도(high-resolution) 해석기법을 적용하여 공학적인 문제에 적용하는 데는 한계가 있다. 이 연구에서는 Boussinesq 근사와 Schmidt 수를 사용하지 않으며, LES 보다 적용이 용이한 DES (detached-eddy simulation)기법을 조합한 다상흐름 수치모델을 적용하여 중력류를 해석을 시도하였다. 수치해석결과를 실험값과 함께 기존 수치모델링 기법으로 구한 수치해와 비교분석하여 이 연구에서 개발 및 적용된 수치모델링 기법의 적용성을 평가한다.

• Water for future
• /
• v.18 no.1
• /
• pp.85-94
• /
• 1985

A method is developed for the separation of the major base flow in a river hydrograph combining the numerical techniques and the empirical methods. The linearized Boussinesq equation and the storage function are used to obtain the base flow recession. The shape of base flow curve made by the recharge of the groundwater table aquifer resulting from rainfall in determined by the Singh and Stall's graphical method, and the continuous from for the curve is approximated by the multiple and polynomial regression. this procedure was successfully tested for the separation of base flow and the establishment of hydrograph in a natural watershed. It was found that the direct numerical method applied to the homogeneous linear second order ordinary differential equation system is not suited to obtain the recession curve, and the case that the loss is generated in the partially penetrating stream can not be solved by the method of this study.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.12 no.4
• /
• pp.181-189
• /
• 2000

A weakly nonlinear mild-slope equation has been derived directly from the continuity equation with the aid of the Galerkin's method. The equation is combined with the momentum equations defined at the mean water level. A single component model has also been obtained in terms of the surface displacement. The linearized form is completely identical with the time-dependent mild-slope equation proposed by Smith and Sprinks(1975). For the verification purposes of the present nonlinear model, the degenerate forms were compared with Airy(1845)'s non-dispersive nonlinear wave equation, classical Boussinesq equation, andsecond¬order permanent Stokes waves. In this study, the present nonlinear wave equations are discretized by the approximate factorization techniques so that a tridiagonal matrix solver is used for each direction. Through the comparison with physical experiments, nonlinear wave model capacity was examined and the overall agreement was obtained.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.6 no.4
• /
• pp.337-343
• /
• 1994

Numerical analysis was applied to a simplified two-dimensional Ondol heating model which consists of heating space on the top of it along with radiant and convective heating floor panel, flue channel in the midway and rectangular underground soil region at the bottom. These three components constitute a system thermally coupled at the top and bottom interfaces of the flue channel. Investigated in the present paper are effects with variations of the Reynolds numbers of 100, 200, and 300, Grashof numbers of $0.1{\times}10^6$ and $0.3{\times}10^6$ and aspect ratios of 15 and 20 on the heat transfer and fluid flow characteristics of two-dimensional Ondol heating model by computer simulation.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.43-43
• /
• 2015

천수방정식을 사용하는 초기 수치모형은 프로드수($F_4$)가 변화하는 흐름 즉, 상류방향과 하류방향으로 전파하는 홍수파를 동시에 해석하기 위해 중앙 차분기법이 필요한 상류(sub-critical flow)와 흐름방향에 따른 상류이송(upwinding)기법이 필요한 사류(super-critical flow)가 나타나는 흐름해석에서 어려움이 있었다. 하지만, 근사 Riemann 해법의 등장으로 흐름방향에 관계없이 특성선을 따라 정확한 상향가중기법의 적용이 가능하게 되어, 천수방정식을 지배방정식으로 하는 수치모형이 더욱 실용적으로 적용될 수 있도록 하였다. 따라서, 현재 근사 Riemann 해법은 Godunov 형 유한체적 기법, 불연속 Galerkin 혹은 Petrov-Galerkin 유한요소기법 그리고 Boussinesq 기법에도 적용되고 있으며, 특히 Godunov 형 유한체적기법과 결합한 근사 Riemann 해법은 댐 붕괴, 하천 범람 그리고 도시 및 해안지역 침수에 이르기까지 여러 가지 문제에 폭넓게 적용되고 있다. 지금까지 홍수 모델링에 적용된 Godunov형 유한체적모형은 정형 사각격자나 비정형 삼각격자 중에서 한가지의 격자 종류만을 적용한 연구가 주로 수행되었으며, 유한요소모형과 같이 이 두 가지 격자를 동시에 적용한 연구는 거의 이루어지지 않고 있다. 일반적으로, 삼각격자는 사각격자와 는 달리 연구유역의 경계나 지형이 복잡한 경우에도 큰 노력없이 격자의 생성이 가능하나, 격자와 노드의 수가 사각격자보다 많아 계산시간이 많이 소요되는 단점이 있다. 반면, 사각격자는 하천과 같이 선형으로 변하는 지형에 대해서는 표현하기가 용이하며 계산시간의 효율성도 뛰어나다. 본 연구에서는 하천, 도시 그리고 해안지역에서의 효율적이고 정확한 홍수 모델링을 위해 삼각 및 사각격자 그리고 이 두 격자를 동시에 고려한 하이브리드 격자의 적용이 가능한 Godunov형 2차원 유한체적 모형을 개발하였다. 그리고 개발모형을 정확해가 있는 댐 붕괴 문제, 실측치가 존재하는 실험하도 및 실제하도에 삼각, 사각 그리고 혼합격자를 생성하여 모의를 수행하고, 각 적용 격자에 따른 정확성과 효율성 및 장점과 단점을 연구하였다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.4
• /
• pp.91-102
• /
• 1987

In this study, the finite element method using 8-node isoparametric element is developed theoretically and simulated to see the deformation of soil and plate, when the circular plate resting on Boussinesq's soil type is loaded axisymmetrically. The results of numerical analysis using the Mohr-Coulomb yield criterion, and experimental analysis are approximative, assuming that soil is elasto-plastic medium. The paper shows that the plastic zone of soil medium is displayed at the near the edge of plate at the first place; when the plastic zone of soil medium is linked around central axis, the external load is termed by critical load, and then the contact pressure changes abruptly, in this case it is approved to be the risk of shear failure.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.499-504
• /
• 2001

A validation of a 3D CFD model for predicting local subcooling of moderator in the vicinity of calandria tubes in a CANDU reactor is performed. The small scale moderator experiments performed at Sheridan Park Experimental Laboratory(SPEL) in Ontario, Canada[1] is used for the validation. Also a comparison is made between previous CFD analyses based on 2DMOTH and PHOENICS, and the current model analysis for the same SPEL experiment. For the current model, a set of grid structures for the same geometry as the experimental test section is generated and the momentum, heat and continuity equations are solved by CFX-4.3, a CFD code developed by AEA technology. The matrix of calandria tubes is simplified by the porous media approach. The standard $k-\varepsilon$ turbulence model associated with logarithmic wall treatment and SIMPLEC algorithm on the body fitted grid are used and buoyancy effects are accounted for by the Boussinesq approximation. For the test conditions simulated in this study, the flow pattern identified is a buoyancy-dominated flow, which is generated by the interaction between the dominant buoyancy force by heating and inertial momentum forces by the inlet jets. As a result, the current CFD moderator analysis model predicts the moderator temperature reasonably, and the maximum error against the experimental data is kept at less than $2.0^{\circ}C$ over the whole domain. The simulated velocity field matches with the visualization of SPEL experiments quite well.

• Computational Structural Engineering
• /
• v.10 no.4
• /
• pp.327-336
• /
• 1997

When a structure is built on the consolidation ground, the instant elastic deflection occures according to the characteristics of the ground and the load on it. And the corresponding contact pressure is established. But, as time passes, the secondary consolidating deflection is added to the instant elastic deflection, the upper structure, due to its flexural rigidity, resist to the additional curvature. So the variation of the contact pressure occurs. And this new contact pressure exerts influence on the consolidation form again. The new consolidation form exerts influence on the contact pressure in return. This kind of interaction continues till all the consolidation of the ground is finished. So the consolidation problem can not be definded as the linear problem. This paper intends to scheme an approximate iteration method to analyse this non-linear interaction problem between the upper structure and the lower consolidation ground which supports the former.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.6
• /
• pp.468-477
• /
• 2019

The numerical model of Boussinesq approximation, which is mainly used for evaluating the port calmness due to the irregular waves, has a limit of applicability of lattice size in ports such as marinas with narrow port openings of around 30m. The SWASH model controls the partial reflection according to the depth, porosity coefficient and structure size when applying the reflected wave incident on the structure and terrain. In this study, the partial reflection evaluation at the front of the structure according to the bottom shape and the shape of the structure are examined. In order to evaluate the reproducibility of the model due to the diffraction waves entering the term, the area of incidence at right angles and inclination of the structure is constructed and compared with the diffraction theory suggested by Goda et al. (1978). The experimental results of the sectional structure reflectances calculated as the depth mean show reflectances similar to the approximate values of the reflectances presented by Stelling and Ahrens (1981). It is considered that the reflected wave is well reproduced according to the control of the reflected wave at the boundary and the shape and topography of the structure. Compared with previous studies to examine the diffraction of the wave incident from the breakwater opening, the wave incidence angle and the shape of the diffraction wave are very similar to the theoretical values, but both oblique and rectangular incidence In the case where the direction concentration is small, the diffraction degree is underestimated in some sections with the crest ratio of 0.5 to 0.6.