• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.20-24
• /
• 2011

In this paper, numerical simulation of cavitation flow for modified NACA66 hydrofoil was made by using the multi-phase RANS equation based on pseudo-compressibility. The Homogeneous mixture model comprised of the mixture continuity, mixture momentum and liquid volume fraction equations was utilized. A vertex-centered finite-volume method was used in conjunction 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing The Spalart-Allmaras one equation model was employed for the closure of turbulence. Reasonable agreements were obtained between the calculation results and the experiment for pressure coefficients on the hydrofoil surface.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.5B
• /
• pp.429-439
• /
• 2009

When Catastrophic extreme flood occurs due to dam break, the response time for flood warning is much shorter than for natural floods. Numerical models can be powerful tools to predict behaviors in flood wave propagation and to provide the information about the flooded area, wave front arrival time and water depth and so on. But flood wave propagation due to dam break can be a process of difficult mathematical characterization since the flood wave includes discontinuous flow and dry bed propagation. Nevertheless, a lot of numerical models using finite volume method have been recently developed to simulate flood inundation due to dam break. As Finite volume methods are based on the integral form of the conservation equations, finite volume model can easily capture discontinuous flows and shock wave. In this study the numerical model using Riemann approximate solvers and finite volume method applied to the conservative form for two-dimensional shallow water equation was developed. The MUSCL scheme with surface gradient method for reconstruction of conservation variables in continuity and momentum equations is used in the predictor-corrector procedure and the scheme is second order accurate both in space and time. The developed finite volume model is applied to 2D partial dam break flows and dam break flows with triangular bump and validated by comparing numerical solution with laboratory measurements data and other researcher's data.

• Journal of Korea Water Resources Association
• /
• v.51 no.10
• /
• pp.853-862
• /
• 2018

An adaptive cut-cell grid based 2D inundation analysis model, K-Flood, is developed in this study. Cut cell grid method divides a grid into a flow area and a non-flow area depending the characteristics of the flows. With adaptive mesh refinement technique cut cell method can represent complex flow area using relatively small number of cells. In recent years, the urban inundation modeling using high resolution and fine quality data is increasing to achieve more accurate flood analysis or flood forecasting. K-Flood has potential to simulate such complex urban inundation using efficient grid generation technique. A finite volume numerical scheme of second order accuracy for space and time was applied. For verification of K-Flood, 1) shockwave reflex simulation by circular cylinder, 2) urban flood experiment simulation, 3) Malpasset dam collapse simulation are performed and the results are compared with observed data and previous simulation results.

• Journal of the Korean Society of Visualization
• /
• v.20 no.3
• /
• pp.49-57
• /
• 2022

The interaction of planar shock wave with rectangular water column is investigated numerically. The flow phenomenon like reflection, transmission, cavitation, recirculation of shock wave, and large negative pressure due to expansion waves was discussed qualitatively and quantitatively. The numerical simulation was performed in a shock tube with a water column, and planar shock was initiated with a pressure ratio of 10. Three cases of the water column with different thicknesses, namely 0.5D, 1D, and 2D, were installed and studied. Water naturally has a higher acoustic impedance than air and mitigates the shock wave considerably. The numerical simulations were modelled using Eulerian and Volume of fluids multiphase models. The Eulerian model assumes the water as a finite structure and can visualize the shockwave propagation inside the water column. Through the volume of fluids model, the stages of breakup of the water column and mitigation effects of water were addressed. The numerical model was validated against the experimental results. The computational results show that the installation of a water column significantly impacts the mitigation of shock wave.

• Structural Engineering and Mechanics
• /
• v.12 no.2
• /
• pp.119-134
• /
• 2001

This paper deals with an experimental and numerical study of a mono-strand wedge anchor head mechanism. First, the experimental program is presented and monitored data such as wedge slippage, anchor deflection and strain distributions along external peripheral surfaces of the anchor head are presented and discussed. In accordance with the experimental set up, these data concern only the global behaviour of the mechanism and cannot provide valuable information such as internal stress-strains distributions, stress concentrations and percentage of yielded volume. Therefore, the second part of this paper deals with the development of an efficient numerical finite element model capable of providing mechanism of the core information. The numerical model which includes all kinematics/material/contact non-linearities is first calibrated using experimental data. Subsequently, a numerical study of the anchorage mechanism is performed and its behaviour is compared to the behaviour of a slightly geometrically modified mechanism where the external diameter has been increased by 5 mm. Finally, different topics influencing the anchorage mechanism behaviour are addressed such as lubrication and wedge shape.

• Journal of Korea Water Resources Association
• /
• v.51 no.10
• /
• pp.895-903
• /
• 2018

There are a large number of weirs installed in rivers of Korea, and these characteristics are not common in other countries. When the flow passes through a structure such as a weir, discontinuous flow occurs. In terms of numerical simulation, it affects the numerical instability due to the balance between the flow term and the source term. In order to solve these problems, many researchers used empirical formulas or numerical scheme simplification. Recently, researches have been conducted to use more accurate numerical scheme. K-River was developed to reflect the characteristics of domestic rivers and calculate the discontinuous flow more accurately. For the verification of K-River, 1) numerical experiment simulations with a bump in the bed, 2) laboratory experiment of hydraulic jump simulation, 3) real river were performed. K-River verified its applicability by simulating results similar to the exact solution and observed value in all simulations.

• Journal of Korea Water Resources Association
• /
• v.47 no.7
• /
• pp.599-614
• /
• 2014

In this study, two dimensional finite volume model was parallelized to improve computing time, which has been developed to be able to apply for the mixed meshes of triangle and quadrilateral. MPI scheme which is free from limitation of the number of cores was applied, and non-blocking point-to-point communication was used for fluxes and time steps calculation domain. The developed model is applied to analyze dam break in a L-shaped experimental channel with $90^{\circ}$ bend and Malpasset dam breach event to calibrate the consistency between parallelized model and existing model and examine the speed-up and efficiency of computing time. Computational speed-up about the size of the input data was considered by simulating 4 cases classified by the number of meshes, Consequently, the simulation results reached a satisfactory accuracy compared to measured data and the results from existing model, and achieved more than 3 times benefit of computational speed-up against computing time of existing model. Simulation results of 3 cases classified by the size of input data lead us to the conclusion that it is important to use proper size of input data and the number of process in order to minimize the communication overhead.

• Computers and Concrete
• /
• v.5 no.3
• /
• pp.175-194
• /
• 2008

Increasingly numerical (finite element) modeling of concrete hinges on our ability to develop a representative volume element with all its heterogeneity properly discretized. Yet, despite all the sophistication of the ensuing numerical models, the initial discretization has been for the most part simplistic. Whenever the heterogeneity of the concrete is to be accounted for, a mesh is often manually crafted through the arbitrary inclusion of the particles (aggregates and/or voids) in an ad-hoc manner. This paper develops a mathematical strategy to precisely address this limitation. Algorithms for the random generation and placement of elliptical (2D) or ellipsoid (3D) inclusions, with possibly radiating cracks, in a virtual concrete model are presented. Collision detection algorithms are extensively used.

• Steel and Composite Structures
• /
• v.27 no.5
• /
• pp.567-576
• /
• 2018

The objective of this work is to analyze large deflections of a fiber reinforced composite cantilever beam under point loads. In the solution of the problem, finite element method is used in conjunction with two dimensional (2-D) continuum model. It is known that large deflection problems are geometrically nonlinear problems. The considered non-linear problem is solved considering the total Lagrangian approach with Newton-Raphson iteration method. In the numerical results, the effects of the volume fraction and orientation angles of the fibre on the large deflections of the composite beam are examined and discussed. Also, the difference between the geometrically linear and nonlinear analysis of fiber reinforced composite beam is investigated in detail.

• Wind and Structures
• /
• v.5 no.2_3_4
• /
• pp.151-164
• /
• 2002

Predictions of the pedestrian level wind speeds for the downtown area of Auckland that have been obtained by wind tunnel and computational fluid dynamic (CFD) modelling are presented. The wind tunnel method involves the observation of erosion patterns as the wind speed is progressively increased. The computational solutions are mean flow calculations, which were obtained by using the finite volume code PHOENICS and the $k-{\varepsilon}$ turbulence model. The results for a variety of wind directions are compared, and it is observed that while the patterns are similar there are noticeable differences. A possible explanation for these differences arises because the tunnel prediction technique is sensitivity to gust wind speeds while the CFD method predicts mean wind speeds. It is shown that in many cases the computational model indicates high mean wind speeds near the corner of a building while the erosion patterns are consistent with eddies being shed from the edge of the building and swept downstream.