• 한국가시화정보학회지
• /
• 제18권3호
• /
• pp.42-51
• /
• 2020

This study performs a numerical simulation of lid driven rectangular cavity flow with different aspect ratios of k = 0.5 to 4 under Reynolds 100, 1,000, 10,000 by using multi-relaxation time (MRT) Lattice Boltzmann Method (LBM). In order to achieve better convergence, well-posed boundary conditions in the domain should be defined such as no-slip condition on side and bottom solid wall surfaces and uniform horizontal velocity on the top of the cavity. This study focuses on the flow inside different shape of rectangular cavity with the aim to observe the effect of the Reynolds number and aspect ratio on the flow characteristics and primary/secondary vortex formation. In order to validate the study, the results have been compared with existing works. The result shows that the Reynolds number and the aspect ratio both has substantial effects on the flow inside the lid-driven rectangular cavity.

• 한국해양공학회지
• /
• 제36권1호
• /
• pp.11-20
• /
• 2022

Generally, tsunamis are generated by the rapid crustal movements of the ocean floor. Other factors of tsunami generation include landslides on coastal and ocean floor slopes, glacier collapses, and meteorite collisions. In this study, two numerical analyses were conducted to examine the formation, propagation, and deformation properties of landslide tsunamis. First, LS-DYNA was adopted to simulate the formation and propagation processes of tsunamis generated by dropping rigid bodies. The generated tsunamis had smaller wave heights and wider waveforms during their propagation, and their waveforms and flow velocities resembled those of theoretical solitary waves after a certain distance. Second, after the formation of the landslide tsunami, a tsunami based on the solitary wave approximation theory was generated in a numerical wave tank (NWT) with a computational domain that considered the stability/steady phase. The comparison of two numerical analysis results over a certain distance indicated that the waveform and flow velocity were approximately equal, and the maximum wave pressures acting on the upright wall also exhibited similar distributions. Therefore, an effective numerical model such as LS-DYNA was necessary to analyze the formation and initial deformations of the landslide tsunami, while an NWT with the wave generation method based on the solitary wave approximation theory was sufficient above a certain distance.

• 대한조선학회논문집
• /
• 제59권3호
• /
• pp.141-148
• /
• 2022

A hybrid turbulence model has developed by combining a sub-grid scale model using dynamic k equation in LES with k-𝜔 SST model of RANS equation. To ascertain potential applicability of the hybrid turbulence model, fully developed turbulent channel flows at Re_𝜏=180 have been simulated of which computational domain has a top wall with coarse cells and a bottom wall with fine cells. The streamwise mean velocity and turbulent intensity profiles showed a good agreement with DNS data when using the hybrid model rather than using a single model in k-𝜔 SST or dynamic k equation models. Computational simulations of turbulent flows around KVLCC2 with a pre-swirl duct have been mainly performed using the hybrid turbulence model. Compared to the results obtained from RANS simulation with k-𝜔 SST model as well as LES with dynamic k equation SGS model, turbulent wakes of the duct in the present simulation using the hybrid turbulence model were very similar to that of LES. Also, the resistances acting on hull, rudder and duct in hybrid turbulence model were similar to those in RANS simulation whereas the viscous forces acting on the hull in LES had a significant error due to coarse cells inappropriate to the sub-grid scale model.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제10권5호
• /
• pp.545-565
• /
• 2018

A numerical code is developed based on potential flow theory to investigate nonlinear sloshing in rectangular Liquefied Natural Gas (LNG) tanks under forced excitation. Using this code, internal free-surface elevation and sloshing loads on liquid tanks can be obtained both in time domain and frequency domain. In the mathematical model, acceleration potential is solved in the calculation of pressure on tanks and the artificial damping model is adopted to account for energy dissipation during sloshing. The Boundary Element Method (BEM) is used to solve boundary value problems of both velocity potential and acceleration potential. Numerical calculation results are compared with published results to determine the efficiency and accuracy of the numerical code. Sloshing properties in partially filled rectangular and membrane tank under translational and rotational excitations are investigated. It is found that sloshing under horizontal and rotational excitations share similar properties. The first resonant mode and excitation frequency are the dominant response frequencies. Resonant sloshing will be excited when vertical excitation lies in the instability region. For liquid tank under rotational excitation, sloshing responses including amplitude and phase are sensitive to the location of the center of rotation. Moreover, experimental tests were conducted to analyze viscous effects on sloshing and to validate the feasibility of artificial damping models. The results show that the artificial damping model with modifying wall boundary conditions has better applicability in simulating sloshing under different fill levels and excitations.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 춘계학술대회논문집D
• /
• pp.97-103
• /
• 2001

The effects of the electrohydrodynamic (EHD) flow and turbulent diffusion on the collection efficiency of a model ESP composed of the plates with a cavity were studied through numerical computation. The electric field and ion space charge density were calculated by the Poisson equation of the electrical potential and the current continuity equation. The EHD flow field was solved by the continuity and momentum equations of the gas phase including the electrical body force induced by the movement of ions under the electric field. The RNG $k-{\varepsilon}$ model was used to analyze the turbulent flow. The particle concentration distribution was calculated from the convective diffusion equation of the particle phase. As the ion space charge increased, the particulate collection efficiency increased because the electrical potential increased over the entire domain in the ESP. The collection efficiency decreased and then increased, i.e. had a minimum value, as the EHD circulating flow became stronger when the electrical migration velocity of the charged particle was low. However, the collection efficiency decreased with the stronger EHD flow when the electrical migration of the particle was higher relatively. The collection efficiency of the model ESP increased as the turbulent diffusivity of the particle increased when the electrical migration velocity of the particle was low. However, the collection efficiency decreased for increasing the turbulent diffusivity when the electrical migration of the particle was higher relatively.

• 한국전산구조공학회논문집
• /
• 제34권3호
• /
• pp.167-174
• /
• 2021

최근 국내 원자력발전소의 격납건물 벽체와 Containment Liner Plate(CLP) 사이에서 다양한 크기의 공극이 발견됨에 따라 원전 격납건물의 보수를 위해 내부 공극의 분포와 크기를 정밀하게 평가할 수 있는 진단기법의 개발이 요구되고 있다. 이에 따라 이 연구에서는 격납건물 벽체에서의 탄성파 전파거동을 계산하는 2차원 유한요소해석 기법을 제시한다. 격납건물 벽체를 기반으로 해석영역을 구성하고 경계면에서의 반사파를 제거하기 위해 수치적 파동흡수 경계층인 perfectly matched layer를 도입하였다. Galerkin 기반 혼합유한요소법을 이용해 2차원 유한영역에서 탄성파 파동방정식의 해를 구하여 충격하중에 대한 격납건물 벽체의 변위와 응력을 계산하였다. 제시한 수치적 기법을 이용하여 격납건물 콘크리트 벽체의 CLP 부착 유무와 공동의 위치 및 크기 변화에 따른 탄성파 전파거동을 살펴보았다. 이 연구의 결과는 원전 격납건물 내부의 공동을 진단하는 탄성파 전체파형 역해석 기법 개발에 활용될 수 있다.

• 한국전산구조공학회논문집
• /
• 제35권5호
• /
• pp.299-308
• /
• 2022

본 논문은 슬로싱 상태에 놓인 포화 상태 액체수소탱크에서 열 유속 및 BOG(Boil-off gas)의 경향을 다루고 있다. 특히, 액체-기체간의 침투 및 혼합에 의한 열 교환에 관심을 두었다. 먼저, VOF(Volume of fluid)와 Eulerian 기반의 다상 유동모델로 모형 슬로싱 실험을 모사하여 압력을 예측하고 계측된 값과 비교하였다. 자유 수면 및 충격 압력 실험 결과와 해석 결과를 비교하였으며, 유체의 속도 예측에서 정확할 수 있음을 간접적으로 증명하였다. 그리고 2차원의 Type-C 원통형 수소탱크를 대상으로 다상열유동해석을 수행하였다. 이때 포화상태에 놓인 액체 및 기체수소를 가정하고, 해석을 통해 각 상간의 혼합에 의한 열 교환의 수준을 확인하고자 하였다. 단, 상간의 열 교환만을 관심으로 두고 있었으므로 질량전달 및 기화모델은 해석에서 제외하였다. 최종적으로 상의 혼합으로 인해 액체수소로 유입되는 열 유속의 기여도에 대하여 정리하였다. 또한 액체수소로 유입되는 열 유속과 집중 질량 기반의 간이식을 통해 BOG 발생량 및 경향을 예측하고 분석하였다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 추계학술대회논문집B
• /
• 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.

• 한국항공우주학회지
• /
• 제39권3호
• /
• pp.232-241
• /
• 2011

하이브리드 로켓은 난류 산화제 유동과 고체 추진제의 기화로 인한 분사 유동 사이의 상호 작용에 의해 복잡한 형태의 혼합 전단층이 존재한다는 특별한 성질을 가지고 있다. 본 논문에서는 유동 간섭에 의해 표면에서 발생하는 진동 유동의 물리적 특성을 연구하기 위하여 압축성 효과를 고려한 질량분사가 있는 덕트 유동의 LES(Large Eddy Simulation) 해석을 수행하였다. 계산 결과에 따르면, 기화 질량이 분출됨에 따라 주유동방향 와류의 특성이 강해지고 국부적으로 발생하는 역류 현상을 근거로 벽면 근방에서 원주방향 와류가 생성됨을 확인하였다. 그리고 시간 특성을 갖고 나타나는 와류 흘림 현상은 혼합 전단층에 기인한 유동 불안정성에 의해 촉진되었으며, 분출유동에 의해 발달한 고유 진동 유동을 의미하는 압력 섭동의 특정 진동수가 $\omega$=8.8에서 검출됨을 확인하였다.

• 에너지공학
• /
• 제8권2호
• /
• pp.224-232
• /
• 1999

캡슐형 잠열재를 이용한 열저장 시스템은 바닥 난방 및 건물 난방에서 매우 효과적인 시스템이다. 이러한 시스템 개발에 필수적인 요소가 열유동 매체가 순환하는 파이프 주변의 캡슐내 온도 분포와 열유동 매체의 유량 등이다. 그러므로 본 연구에서는 3차원 비정상 상태에서 Navier-Stokes 방정식, 난류모델을 비롯한 스칼라 보존 방정식을 적용하여 캡슐 블록의 온도 분포 및 파이프 내의 유동장 해석을 수행하였다. 또한 본 연구와 같이 계산 영역이 특별한 기하학적 현상을 형상(circle+square)인 문제 해결하는데 적용할 수 있는 새로운 격자 생성 기술(MBFGE/CCM)을 개발하였다. 격자계는 파이프에서 원형 격자를 이용하였고, 캡슐 블록에서 사각 격자를 이용하여 다중격자와 미세격자를 결합하여 사용하였다. 본 연구의 목적은 컴퓨터를 이용한 수치해석적 방법을 미세 캡슐을 이용한 축열보드에 적용하여 2종류의 열경계 상태에 대하여 속도와 온도분포를 계산하여 비교분석을 하는 것이다. 온도는 축열 보드의 한 쪽면은 대류면이고 다른 한쪽면은 단열면인 경우(Case 2)보다 양면 모두 단열인 경우(Case 1)일 때 더 높게 상승하였다. 온수 파이프 중심선인 Y=0 에 가까운 영역에서 Case 1과 Case 2사이에 축열 보드 내에서 온도 차이는 확연하게 나타났다. 향후 수치해석의 정확도를 높이고 축열 보드의 열전달 현상을 보다 정확히 계산하기 위해서는 위치 및 시간에 따른 정밀한 온도 측정값이 필요하고 특히 잠열재인 미세 캡슐이 상변화를 하므로 온도 변화에 따른 물질의 비열(C$_{p}$)과 열전달율(λ)을 고려한 방정식이 요구된다.