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

Investigating Backward-Facing Step(BFS) flow is important in that it is a representative case for separation flows in various engineering flow systems. There have been a wide range of experimental, theoretical, and numerical studies to investigate the flow characteristics over BFS, such as flow separation, reattachment length and recirculation zone. However, most of such previous studies were concentrated only on the perpendicular step angle. In this study, several numerical investigations on the flow pattern over BFS with various step angles (10° ~ 90°) and expansion ratios (1.48, 2 and 3.27) under different Reynolds numbers (5000 ~ 64000) were carried out, mainly focused on the reattachment length. The numerical simulations were performed using an open source 3D CFD software, OpenFOAM, in which the velocity profiles and turbulence intensities are calculated by RANS (Reynolds Averaged Navier-Stokes equation) and 3D LES (Large Eddy Simulation) turbulence models. Overall, it shows a good agreement between simulations and the experimental data by Ruck and Makiola (1993). In comparison with the results obtained from RANS and 3D LES, it was shown that 3D LES model can capture much better and more details on the velocity profiles, turbulence intensities, and reattachment length behind the step for relatively low Reynolds number(Re < 11000) cases. However, the simulation results by both of RANS and 3D LES showed very good agreement with the experimental data for the high Reynolds number cases(Re > 11000). For Re > 11000, the reattachment length is no longer dependent on the Reynolds number, and it tends to be nearly constant for the step angles larger than 30°.) Based on the calibrated and validated numerical simulations, several additional numerical simulations were also conducted with higher Reynolds number and another expansion ratio which were not considered in the experiments by Ruck and Makiola (1993).

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2001년도 춘계 학술대회논문집
• /
• pp.9-14
• /
• 2001

Astrophysics has been one of the disciplines which utilize actively supercomputers for their researches. In this paper, studies of astrophysical flows by large scale numerical simulations in supercomputers are discussed. The followings are described: 1) characteristics of astrophysical flows, 2) numerical codes to study astrophysical flows, 3) performance of the codes in parallel supercomputers. A couple of examples of such studies are briefly presented: 1) three-dimensional evolution of the nonlinear Kelvin-Helmholtz instability in a magnetized medium, 2) three-dimensional simulations of astrophysical jets in a magnetized medium.

• 대한수학회논문집
• /
• 제22권3호
• /
• pp.453-464
• /
• 2007

A computationally efficient numerical scheme is presented for the phase-field model of two-phase systems for anisotropic interfacial energy. The scheme is solved by using a nonlinear multigrid method. When the coefficient for the anisotropic interfacial energy is sufficiently high, the interface of the system shows corners or missing crystallographic orientations. Numerical simulations with high and low anisotropic coefficients show excellent agreement with exact equilibrium shapes. We also present spinodal decomposition, which shows the robustness of the pro-posed scheme.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1997년도 추계 학술대회논문집
• /
• pp.10-15
• /
• 1997

Unsteady natural convection of an enclosed fluid has been one of the fundamental thermo-fluid problems, of which dynamic relevance is found in many engineering applications. Together with the inherent coupling between the boundary layers and the interior core, and strong interaction between flow and temperature fields, the unsteadiness poses serious hurdles for analytical and experimental approaches. With the recent development of computers and solution algorithms, computational fluid dynamics has become the prevailing tool to tackle the underlying problems. In this presentation, a few examples of numerical studies are introduced. The usefulness and potential of numerical simulations in investigating unsteady natural convection are elaborated.

• Structural Engineering and Mechanics
• /
• 제40권5호
• /
• pp.617-636
• /
• 2011

We present the numerical implementation, simulation, and validation of the high-velocity impact experiments that have been described in the companion article. In this part, numerical investigations and simulations performed to mimic the tests are presented. The experiments were analyzed by the explicit integration-based software ABAQUS for improved simulations. Targets were modeled with a damaged plasticity model for concrete. Computational results of residual velocity and crater dimensions yielded acceptable results.

• 한국환경과학회지
• /
• 제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.

• Geomechanics and Engineering
• /
• 제5권3호
• /
• pp.263-281
• /
• 2013

Geotextiles and geogrids have been in use for several decades in variety of geo-structure applications including foundation of embankments, retaining walls, pavements. Geocells is one such variant in geosynthetic reinforcement of recent years, which provides a three dimensional confinement to the infill material. Although extensive research has been carried on geocell reinforced sand, clay and layered soil subgrades, limited research has been reported on the aggregates/ballast reinforced with geocells. This paper presents the behavior of a railway sleeper subjected to monotonic loading on geocell reinforced aggregates, of size ranging from 20 to 75 mm, overlying soft clay subgrades. Series of tests were conducted in a steel test tank of dimensions $700mm{\times}300mm{\times}700mm$. In addition to the laboratory model tests, numerical simulations were performed using a finite difference code to predict the behavior of geocell reinforced ballast. The results from numerical simulations were compared with the experimental data. The numerical and experimental results manifested the importance that the geocell reinforcement has a significant effect on the ballast behaviour. The results depicted that the stiffness of underlying soft clay subgrade has a significant influence on the behavior of the geocell-aggregate composite material in redistributing the loading system.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 춘계학술대회논문집
• /
• pp.675-678
• /
• 2008

The effect of wave and current interactions on jet-like current flowing against waves was investigated based on numerical simulations. The numerical simulations are conducted by a combined model system of REF/DIF(a wave model) plus SHORECIRC(a current model) and a Boussinesq equation model, FUNWAVE. In the simulations, regular and irregular waves refracted due to the jet-like opposing current were focused along the core region of current, and the jet-like current was earlier spreaded when the waves had larger wave heights. The numerical results show that the rapid change of wave height distribution in transverse direction near current inlet plays a significant role to spread the jet-like current. In other words, the gradients of radiation stress forcing in transverse direction have a more significant effect on the jet-like current than its relatively small gradients forcing in flowing direction, which tend to accelerate the current, do. In conclusion, it is indispensible to take into account the interaction effect of wave transformation and current characteristics when waves meet the opposing jet-like current such as river mouth.

• Journal of Ship and Ocean Technology
• /
• 제5권4호
• /
• pp.29-43
• /
• 2001

In this paper, free surface flows around an advancing two-dimensional rectangular cylinder piercing the free surface are studied using numerical and experimental methods. Especially, wave breaking phenomenon around the cylinder is treated in detail. A series of numerical simulations and experiments were performed for the purpose of comparison. For the numerical simulations, a finite difference method was adopted with a rectangular grid system, and the variation of the free surface was computed by the marker density method. The computational results are compared with the experiments. It is confirmed that the present numerical method is useful for the numerical simulation of nonlinear free surface waves around a piercing body.

• 천문학회지
• /
• 제34권4호
• /
• pp.321-323
• /
• 2001

Compressible, magnetohydrodynamic (MHD) turbulence in two dimension is studied through high-resolution, numerical simulations with the isothermal equation of state. First, hydrodynamic turbulence with Mach number $(M)_{rms}\;\~$1 is generated by enforcing a random force. Next, initial, uniform magnetic field of various strengths with Alfvenic Mach number Ma $\gg$ 1 is added. Then, the simulations are followed until MHD turbulence is fully developed. Such turbulence is expected to exist in a variety of astrophysical environments including clusters of galaxies. Although no dissipation is included explicitly in our simulations, truncation errors produce dissipation which induces numerical resistivity. It mimics a hyper-resistivity in our second-order accurate code. After saturation, the resulting flows are categorized as SF (strong field), WF (weak field), and VWF (very weak field) classes respectively, depending on the average magnetic field strength described with Alfvenic Mach number, $(Ma)_{rms}{\ge}1$, $(Ma)_{rms}{\~}1$, and $(Ma)_{rms}{\gg}1$. The characteristics of each class are discussed.