• Geomechanics and Engineering
• /
• 제29권3호
• /
• pp.269-279
• /
• 2022

In this research, a series of dynamic numerical analysis were carried out for deep underground building structures under the various earthquake conditions. Dynamic numerical analysis model was developed based on the PLAXIS2D and calibrated with centrifuge test data from Kim et al. (2016). The hardening soil model with small strain stiffness (HSSMALL) was adopted for soil constitutive model, and interface elements was employed at the interface between plate and soil elements to simulate dynamic interaction effect. In addition, parametric study was performed for fixed condition and embedded depth. Finally, the dynamic behavior of underground building structure was thoroughly analyzed and evaluated.

• 한국환경과학회지
• /
• 제23권7호
• /
• pp.1223-1232
• /
• 2014

Density currents have been easily observed in environmental flows, for instance turbidity currents and pollutant plumes in the oceans and rivers. In this study, we explored the propagation dynamics of density currents using the FLOW-3D computational fluid dynamics code. The renormalization group (RNG) $k-{\varepsilon}$ scheme, a turbulence numerical technique, is employed in a Reynold-averaged Navier-Stokes framework (RANS). The numerical simulations focused on two different types of intrusive density flows: (1) propagating into a two-layer ambient fluid; (2) propagating into a linearly stratified fluid. In the study of intrusive density flows into a two-layer ambient fluid, intrusive speeds were compared with laboratory experiments and analytical solutions. The numerical model shows good quantitative agreement for predicting propagation speed of the density currents. We also numerically reproduced the effect of the ratio of current depth to the overall depth of fluid. The numerical model provided excellent agreement with the analytical values. It was also clearly demonstrated that RNG $k-{\varepsilon}$ scheme within RANS framework is able to accurately simulate the dynamics of density currents. Simulations intruding into a continuously stratified fluid with the various buoyancy frequencies are carried out. These simulations demonstrate that three different propagation patterns can be developed according to the value of $h_n/H$ : (1) underflows developed with $h_n/H=0$ ; (2) overflows developed when $h_n/H=1$ ; (3) intrusive interflow occurred with the condition of 0 < $h_n/H$ < 1.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제6권2호
• /
• pp.219-235
• /
• 2014

With the development of the technology of underwater moving bodies, the need for developing the knowledge of surface effect interaction of free surface and underwater moving bodies is increased. Hence, the two-phase flow is a subject which is interesting for many researchers all around the world. In this paper, the non-linear free surface deformations which occur during the water-exit of a circular cylinder due to its buoyancy are solved using finite volume discretization based code, and using Volume of Fluid (VOF) scheme for solving two phase flow. Dynamic mesh model is used to simulate dynamic motion of the cylinder. In addition, the effect of cylinder mass in presence of an external force is studied. Moreover, the oblique exit and entry of a circular cylinder with two exit angles is simulated. At last, water-exit of a circular cylinder in six degrees of freedom is simulated in 3D using parallel processing. The simulation errors of present work (using VOF method) for maximum velocity and height of a circular cylinder are less than the corresponding errors of level set method reported by previous researchers. Oblique exit shows interesting results; formation of waves caused by exit of the cylinder, wave motion in horizontal direction and the air trapped between the waves are observable. In 3D simulation the visualization of water motion on the top surface of the cylinder and the free surface breaking on the front and back faces of the 3D cylinder at the exit phase are observed which cannot be seen in 2D simulation. Comparing the results, 3D simulation shows better agreement with experimental data, specially in the maximum height position of the cylinder.

• 한국산학기술학회논문지
• /
• 제15권1호
• /
• pp.495-505
• /
• 2014

강변저류지는 홍수량을 일시적으로 분담하여 홍수위험도를 저감시키는 시설로 국내에서도 대하천 정비사업을 통하여 여주, 나주, 영월 지역의 홍수 저감을 위하여 강변저류지가 새로이 설치되었다. 본 연구에서는 신규 도입된 강변저류지 중 여주 강변저류지를 대상으로 1차원 및 2차원 수치모의를 수행하여 대하천에서의 강변저류지에 의한 홍수위 저감효과를 분석하였다. 1차원 모의는 HEC-RAS 및 FLDWAV 모형의 저류지 기능을 적용하여 부정류 모의를 통한 홍수위 저감효과를 분석하였으며, 2차원 모의는 CCHE2D 모형을 이용하여 1차원 모형과 동일한 홍수 사상을 적용한 후 홍수위 저감효과를 분석하였다. 1차원 모의에서는 월류제 부근에서 첨두홍수위가 가장 크게 감소하였으며, 2차원 모의에서는 모의 구간의 가장 상류에서 홍수위 저감 효과가 가장 크게 나타났다. 또한 첨두홍수위는 1차원 모형에서 최대 0.13 m, 2차원 모형에서는 최대 0.20 m가 감소하는 것으로 나타났다. 검토 결과 2차원 모형이 하도 흐름에 대한 지형의 물리적 특성이 반영되어 보다 합리적인 모의 결과를 얻는 것이 가능한 것으로 보이며, 횡월류위어의 부근에서 하도 내 흐름 거동뿐만 아니라 월류 유량이 유입된 후 강변저류지 내 거동을 검토하기에는 1차원 모형보다 유리한 것으로 나타났다.

• Wind and Structures
• /
• 제22권1호
• /
• pp.17-41
• /
• 2016

To address the uncertainty of the flight trajectories caused by the turbulence and gustiness of the wind field over the roof and in the wake of a building, a 3-D probabilistic trajectory model of flat-type wind-borne debris is developed in this study. The core of this methodology is a 6 degree-of-freedom deterministic model, derived from the governing equations of motion of the debris, and a Monte Carlo simulation engine used to account for the uncertainty resulting from vertical and lateral gust wind velocity components. The influence of several parameters, including initial wind speed, time step, gust sampling frequency, number of Monte Carlo simulations, and the extreme gust factor, on the accuracy of the proposed model is examined. For the purpose of validation and calibration, the simulated results from the 3-D probabilistic trajectory model are compared against the available wind tunnel test data. Results show that the maximum relative error between the simulated and wind tunnel test results of the average longitudinal position is about 20%, implying that the probabilistic model provides a reliable and effective means to predict the 3-D flight of the plate-type wind-borne debris.

• 한국전산유체공학회지
• /
• 제20권2호
• /
• pp.46-51
• /
• 2015

In this paper, we have developed numerical model for particulated flow through narrow slit using Eulerian-Granular method. Commercial software (FLUENT) was utilized as simulation tool and main focus was to identify the effect from various numerical options for modeling of solid particles as continuos phase in granular flow. Gidaspow model was chosen as basic model for solid viscosity and drag model. And lun-et-al model was used as solid pressure and radial distribution model, respectively. Several other model options in FLUENT were tested considering the cross related effect. Mass flow rate of the particulate through the slit was measured to compare. Due to the high volume density of the stacked particulates above the slit, effect from various numerical options were not significant. The numerical results from basic model were also compared with experimental results and showed very good agreement. The effects from the characteristics of particles such as diameter, angle of internal friction, and collision coefficient were also analyzed for future design of velocity resistance layer in solar thermal absorber. Angle of internal friction was found to be the dominat variable for the particle mass flow rate considerably. More defined 3D model along with energy equation for complete solar thermal absorber design is currently underway.

• 한국산학기술학회논문지
• /
• 제18권3호
• /
• pp.14-21
• /
• 2017

본 최근 지구의 온난화로 인하여 극한 홍수가 자주 발생하고 있으며, 기존 도시 유역의 우수배제시설 용량부족 등으로 도시홍수가 빈번하게 발생하고 있으므로, 주요 범람지역의 홍수량을 우회하거나 저류하여 홍수를 방지하기 위한 수로터널의 개발이 요구된다. 본 연구에서는 교통기능과 수로기능을 동시에 갖춘 다기능 수로 터널의 설계 기준을 개발하기 위한 수리 실험 및 Flow3D를 이용한 수치모의을 수행하였다. 수치모의결과 동일한 수로 터널 구간 내 발생하는 마찰손실의 크기는 수치모의로 도출된 마찰손실이 이론적으로 계산한 마찰손실보다 더 크게 발생함을 관측하였으며, 이는 수로의 형상이 비원형인 경우에는 관의 기하학적 형상에 의한 흐름구조의 변화로 추가적인 마찰손실이 발생하는 것이 원인으로 판단된다. 마찰손실의 증가는 난류보다 층류에서 두드러졌다. 따라서 터널의 홍수량 흐름 시 마찰손실계수가 터널의 형상에 좌우되며, 실무에서 정확한 설계를 위해 방수로 터널의 형상을 주의 깊게 고려해야 한다는 결론을 내렸다. 이는 실제 방수로 터널 설계에 활용될 수 있는 기본 정보를 제공할 수 있을 것으로 보인다.

• 대한기계학회논문집B
• /
• 제28권12호
• /
• pp.1608-1615
• /
• 2004

CFD data compression methods based on Full-3D and Quasi-1D supercompact multiwavelets are presented. Supercompact wavelets method provide advantageous benefit that it allows higher order accurate representation with compact support. Therefore it avoids unnecessary interaction with remotely located data across singularities such as shock. Full-3D wavelets entails appropriate cross-derivative scaling function & wavelets, hence it can allow highly accurate multi-spatial data representation. Quasi-1D method adopt 1D multiresolution by alternating the directions rather than solving huge transformation matrix in Full-3D method. Hence efficient and relatively handy data processing can be conducted. Several numerical tests show swift data processing as well as high data compression ratio for CFD simulation data.

• 정보과학회 논문지
• /
• 제42권3호
• /
• pp.281-288
• /
• 2015

화산재해 피해 예측 시스템은 기상과 화산분화 시뮬레이션 결과를 기반으로 화산재해대응을 위한 판단을 도와주는 시스템이다. 이 시스템에서 Fall3D라는 프로그램은 기상정보를 바탕으로 화산분화 이후 화산재의 확산예측결과를 생성하고 기상정보를 생성하기 위해 WRF라는 기상수치예보모델을 사용한다. 두 시뮬레이션의 프로그램을 수정하지 않고, 전체 실행시간을 줄이기 위해서는 WRF의 기상예측모델의 시간별 부분결과가 발생할 때마다 Fall3D를 부분수행 할 수 있도록 하는 파이프라이닝 방식을 생각할 수 있다. 이를 위해서 Fall3D와 같은 후속계산은 선행계산의 부분결과가 생성될 때까지 일시정지하고, 계산에 필요한 정보가 발생하면 재개할 수 있어야한다. 비록 Fall3D가 이런 일시정지와 재개기능을 가지고 있지는 않지만 그 이전 계산을 이어서 진행할 수 있는 재시작기능을 활용하여 파이프라이닝 효과를 낼 수 있다. 본 논문에서는 이러한 실행 형태를 제어할 수 있는 워크플로우를 제안한다.

• 한국지진공학회논문집
• /
• 제21권5호
• /
• pp.255-264
• /
• 2017

Dynamic numerical simulation of pile-supported slab track system embedded in a soft soil and embankment was performed. 3D model was formulated in a time domain to consider the non-linearity of soil by utilizing FLAC 3D, which is a finite difference method program. Soil non-linearity was simulated by adopting the hysteric damping model and liner elements, which could consider soil-pile interface. The long period seismic loads, Hachinohe type strong motions, were applied for estimating seismic respose of the system, Parametric study was carried out by changing subsoil layer profile, embankment height and seismic loading conditions. The most of horizontal permanent displacement was initiated by slope failure. Increase of the embedded height and thickness of the soft soil layer leads increase of member forces of PHC piles; bending moment, and axial force. Finally, basic guidelines for designing pile-supported slab track system under seismic loading are recommended based on the analysis results.