• Journal of Korea Water Resources Association
• /
• v.50 no.5
• /
• pp.289-302
• /
• 2017

DFB (Divergence Form for Bed slope source term) was rigorously derived and the error of mDFB using mean water depth at the cell face in DFB was clearly demonstrated. In addition, DFB technique turned out to be an exact method to the bed slope source term. The existing volume/free-surface relationship to the PSC (Partially Submerged Cell) has been corrected. It was discussed that treatment for the partially submerged edge is required to satisfy the C-property in PSC. It is expected that this study will provides a more accurate means in analyzing the shallow water equations with the approximate Riemann solver.

• Journal of the Korean Society of Visualization
• /
• v.20 no.2
• /
• pp.55-62
• /
• 2022

We experimentally investigate a rising toroidal bubble impacting a free surface. The toroidal bubble is created by releasing pulsed air. By adjusting the volume and circulation of the toroidal bubble, the characteristics of interactions between the toroidal bubble and the free surface are identified. Because of the impact by the toroidal bubble, the free surface is convexly deformed upwards above the center point of the toroidal bubble, while the edge of the deformed free surface is pulled down. When the circulation of the bubble becomes stronger, the surface which was pulled down breaks eventually, and air above the free surface is entrained into water, forming an unstable toroidal bubble. The deformations at the center and edge of the free surface are in a linear relationship with the Froude number and the Weber number, respectively.

• Journal of Korea Water Resources Association
• /
• v.46 no.2
• /
• pp.139-153
• /
• 2013

For analyzing shallow-water flows over the uneven bottom, the HLLL scheme and the divergence form for bed slope source term (DFB) technique, respectively were applied to the flux gradient and the bottom gradient source terms in a finite-volume model for the shallow water equations. And also the model incorporated the volume/free-surface relationship (VFR) to consider the partially submerged cells (PSC). It was identified that a simpler version of the weighted surface-depth gradient method in the MUSCL was equivalent to the original one in the accuracy for 1D steady flows. It was verified that the flux gradient term and the bottom gradient source term were well-balanced exactly by the VFR for the 1D PSC. The VFR for the triangular PSC settled the problem which the governing equations were not well-balanced by the DFB technique for the 2D PSC. There were good agreements in simulations and experiments for 2D dam-break flows over a triangular sill and a round bump. In addition, the partial dam-break flow was successfully simulated for flooding of roughnesses in an irregular bottom as well as a sloping one. Therefore, this model is expected to be applied to the real river with uneven topography.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.2
• /
• pp.149-159
• /
• 2017

In this study, the numerical investigation of ship motions and added resistance at constant forward velocity of KVLCC2 model is presented. Finite volume CFD code is used to calculate three dimensional, incompressible, unsteady RANS equations. Numerical computations show that reliable numerical results can be obtained in head waves. In the numerical analyses, body attached mesh method is used to simulate the ship motions. Free surface is simulated by using VOF method. The relationship between the turbulence viscosity and the velocities are obtained through the standard ${\kappa}-{\varepsilon}$ turbulence model. The numerical results are examined in terms of ship resistance, ship motions and added resistance. The validation studies are carried out by comparing the present results obtained for the KVLCC2 hull from the literature. It is shown that, ship resistance, pitch and heave motions in regular head waves can be estimated accurately, although, added resistance can be predicted with some error.

• Bulletin of the Korean Chemical Society
• /
• v.24 no.12
• /
• pp.1742-1750
• /
• 2003

In order to characterize the hydrophobic parameters of N-acetyl amino acid amides in 1-octanol/water, a theoretical calculation was carried out using a solvation free energy density model. The hydrophobicity parameters of the molecules are obtained with the consideration of the solvation free energy over the solvent volume surrounding the solute, using a grid model. Our method can account for the solvent accessible surface area of the molecules according to conformational variations. Through a comparison of the hydrophobicity of our calculation and that of other experimental/theoretical works, the solvation free energy density model is proven to be a useful tool for the evaluation of the hydrophobicity of amino acids and peptides. In order to evaluate the solvation free energy density model as a method of calculating the activity of drugs using the hydrophobicity of its building blocks, the contracture of Bradykinin potentiating pentapeptide was also predicted from the hydrophobicity of each residue. The solvation free energy density model can be used to employ descriptors for the prediction of peptide activities in drug discovery, as well as to calculate the hydrophobicity of amino acids.

• 한국전기화학회:학술대회논문집
• /
• 2004.06a
• /
• pp.287-292
• /
• 2004

In order to understand the relationship between molecular structure of Metal-Organic Framework(MOF) and capacity of hydrogen absorption, quantum mechanical calculations and grand canonical Monte Carlo simulations have been carried out on a series of MOF-5 having various organic linkers. The calculation results about specific surface area and electron density for various frameworks indicated that the capacity of the hydrogen storage is largely dependent on effective surface area rather than the free volume. Based on the iso-electrostatic potential surface from density functional calculation and the amount of adsorbed hydrogens from grand canonical Monte Carlo calculation, it was also found that the electron localization ground organic linker plays an important role in hydrogen capacity of MOFs.

• Korean Journal of Materials Research
• /
• v.20 no.11
• /
• pp.581-585
• /
• 2010

High strength sheet steels for automobile are seriously compromised by hydrogen embrittlement. This issue has been continuously studied, but the field of interest, which lies between microstructural characteristics and hydrogen behavior with hydrogen charging, has not yet been thoroughly investigated. This study was done to investigate the behavior of hydrogen according to the hydrogen volume fraction on 590MPa grade DP steels, which are developed under hydrogen charging conditions as high strength sheet steels for automobiles. The penetration depths and the mechanical properties, according to charging conditions, were investigated through the distribution of micro-hardness and the microstructural observation of the subsurface zone. It was found that the amount of hydrogen trapping in 590MPa DP steels was related to the austenite volume fraction. It was confirmed that the distribution of micro-hardnesses according to the depth of the subsurface zone under the free surface showed the relationship of the depth of the hydrogen saturation between the charging conditions.

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

계산 격자에 기반하여 천수 흐름을 모의할 때, 그 격자에 담긴 물의 양을 정확하게 파악할 필요가 있다. 예를 들어, 초기조건으로 수위가 부여된다면 계산격자의 기하 특성에 맞추어 흐름 변수인 수심이나 흐름 단면적으로 바꾸어야하기 때문이다. 필요에 따라서는 모의 결과를 수위로 보이거나 격자 속 수심을 계산에 사용할 수도 있으므로 그 역변환도 고려되어야 한다. 2차원의 삼각형 계산격자에 대해서는 물의 부피와 수위 관계(volume/free-surface relationship)가 이미 정확(exact)하게 구명되어 있다(Hwang, 2017, J. KWRA). 그런데 1차원 문제의 횡단면에서 흐름 단면적과 수위의 관계(area/free-surface relationship)는 수위로부터 면적 환산에 대해서는 정확하나 그 역변환은 그렇지 않다. 매 시간 단계에서 갱신된 흐름 단면적으로부터 수위를 환산하기 위해 미리 작성된 면적-수위 자료를 이용한 선형 보간이 적용된다(Goodell, 2011, The RAS Solution). 이때, 환산 정확도는 자료의 해상도에 의존된다. 다행히 하천 횡단면 대부분을 채워 흐르는 홍수모의에서는 이 문제가 그리 심각하지는 않다. 심지어 수위가 복단면 저수로 턱에 걸쳐있어 흐름단면적이 급변하는 경우에도 환산 수위의 정확도는 크게 훼손되지 않는다. 그러나 미미한 환산 오차일지라도 그로 인해 수위가 저수로 턱을 넘거나 그보다 작을 수 있다. 이 경우, 홍수터의 잠김여부에 따라 수면폭(top width)이 실제 계산 결과에 비해 크게 달라질 수밖에 없다. 수면폭 오차는 그것을 이용하여 결정되는 수리 수심(hydraulic depth)이나 평균 하상고(mean bed level)의 산정에도 전파된다. 이 연구에서는 하천 횡단면에서 수위와 흐름 단면적 사이의 환산 정확도를 크게 높일 수 있는 기법을 제시하였다. 먼저 하천 횡단면에서 주어진 수위에 대해 흐름 단면적을 산정할 수 있는 알고리듬을 보였다. 또한, 횡단면에서 수위와 흐름 단면적의 관계가 단조 증가 함수(monotonically increasing function)임에 착안하여 그 역변환에 대해 해 찾기(root finding) 방법의 하나인 Brent 기법을 적용하였다. 이 기법은 주어진 구간에서 도함수가 알려져 있지 않은 경우에 대해서도 효과적으로 해를 찾을 수 있는 것으로 알려져 있다(Press et al., 2002, Numerical Recipes in C, 2nd Ed.). 내성천 하류 수계의 333개 단면에서 수면폭에 대한 상대 오차를 살펴보면, 선형 보간에 의한 기존 방법으로는 면적-수위 자료의 수가 1,000개가 되어도 그 최대치가 1% 이내에 들지 않은 반면, 이 연구에서 제시한 기법으로 면적-수위 자료 없이도 1% 이내로 줄어드는 것을 확인하였다. 다만, 반복 계산에 의한 계산 시간의 증대를 피할 수 없다. 미리 작성된 면적-수위 자료를 이용하면 계산 비용을 줄일 수 있으며, 약 35개의 구간으로 나누었을 때 비용 대비오차가 적절하였다. 이 연구는 한국건설기술연구원(주요사업 과제번호: 20190116-001)의 지원에 의한 것이다.

• Geomechanics and Engineering
• /
• v.13 no.6
• /
• pp.977-995
• /
• 2017

It is well accepted that rock failure mechanism influence the cutting efficiency and determination of optimum cutting parameters. In this paper, an attempt was made to research the factors that affect the failure mechanism based on discrete element method (DEM). The influences of cutting depth, hydrostatic pressure, cutting velocity, back rake angle and joint set on failure mechanism in rock-cutting are researched by PFC2D. The results show that: the ductile failure occurs at shallow cutting depths, the brittle failure occurs as the depth of cut increases beyond a threshold value. The mean cutting forces have a linear related to the cutting depth if the cutting action is dominated by the ductile mode, however, the mean cutting forces are deviate from the linear relationship while the cutting action is dominated by the brittle mode. The failure mechanism changes from brittle mode with larger chips under atmospheric conditions, to ductile mode with crushed chips under hydrostatic conditions. As the cutting velocity increases, a grow number of micro-cracks are initiated around the cutter and the volume of the chipped fragmentation is decreasing correspondingly. The crack initiates and propagates parallel to the free surface with a smaller rake angle, but with the rake angle increases, the direction of crack initiation and propagation is changed to towards the intact rock. The existence of joint set have significant influence on crack initiation and propagation, it makes the crack prone to propagate along the joint.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.2
• /
• pp.883-898
• /
• 2019

This paper employs computational tools to predict power increase (or speed loss) and propulsion performances in waves of KVLCC2. Two-phase unsteady Reynolds averaged Navier-Stokes equations have been solved using finite volume method; and a realizable k-ε model has been applied for the turbulent closure. The free-surface is obtained by solving a VOF equation. Sliding mesh method is applied to simulate the flow around an operating propeller. Towing and self-propulsion computations in calm water are carried out to obtain the towing force, propeller rotating speed, thrust and torque at the self-propulsion point. Towing computations in waves are performed to obtain the added resistance. The regular short head waves of λ/LPP = 0.6 with 4 wave steepness of H/λ = 0.007, 0.017, 0.023 and 0.033 are taken into account. Four methods to predict speed-power relationship in waves are discussed; Taylor expansion, direct powering, load variation, resistance and thrust identity methods. In the load variation method, the revised ITTC-78 method based on the 'thrust identity' is utilized to predict propulsive performances in full scale. The propulsion performances in waves including propeller rotating speed, thrust, torque, thrust deduction and wake fraction, propeller advance coefficient, hull, propeller open water, relative rotative and propulsive efficiencies, and delivered power are investigated.