• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.193-200
• /
• 2009

Recently, there have been efforts to construct hybrids among the existing methodologies for multiphase flow such as VOF, Level Set, and Front Tracking with the intention of facilitating simulations of general three-dimensional problems. As one of the hybrid method, we have developed the Level Contour Reconstruction Method (LCRM) for general three-dimensional multiphase flows including phase change. The main idea was focused on simplicity and a robust algorithm especially for the three-dimensional case. It combines characteristics of both Front Tracking and Level Set methods. While retaining an explicitly tracked interface using interfacial elements, the calculation of a vector distance function plays a crucial role in the periodic reconstruction of the interface elements in the LCRM method to maintain excellent mass conservation and interface fidelity. In addition, compact curvature formulation is incorporated for the calculation of the surface tension force thereby reducing parasitic currents to a negligible level.

• Computers and Concrete
• /
• v.8 no.3
• /
• pp.327-341
• /
• 2011

Three-dimensional graphic objects created by MATLAB are exported to the AUTOCAD program through the MATLAB handle functions. The imported SAT format files are used to produce the finite element mesh for MSC.PATRAN. Based on the Monte-Carlo random sample principle, the material heterogeneity of cement composites with randomly distributed fibers is described by the WEIBULL distribution function. In this paper, a concept called "soft region" including micro-defects, micro-voids, etc. is put forward for the simulation of crack propagation in fiber-reinforced cement composites. The performance of the numerical model is demonstrated by several examples involving crack initiation and growth in the composites under three-dimensional stress conditions: tensile loading; compressive loading and crack growth along a bimaterial interface.

• Korea-Australia Rheology Journal
• /
• v.12 no.3_4
• /
• pp.165-173
• /
• 2000

Three-dimensional numerical simulation of isothermal/nonisothermal coextrusion process of two immiscible polymers through a rectangular channel has been done using the finite element method. The encapsulation phenomenon with the less viscous layer encapsulating the more viscous layer was investigated with the generalized Newtonian fluids. The interface position around the symmetric plane obtained by numerical simulation nearly coincided with the one observed in experiments, but the degree of encapsulation was less than the one observed experimentally. Open boundary condition method was found to be applied to the simulation of nonisothermal coextrusion process, however, the results are not far from those using the fully developed boundary condition, because the temperature development along the downstream direction is very slow in the case of convection dominated flow. When the inlet velocity is increased, the interface profile does not change in isothermal flow, while it moves upward in nonisothermal situation. The degree of encapsulation decreases along the downstream direction in nonisothermal flow. When the inlet temperature increases compared to the wall temperature, the outlet interface moves downward and the degree of encapsulation increases. The difference of degree of encapsulation between the simulation and the experiments seems to arise from the viscoelastic effect of the materials. It was concluded that the nonisothermal effect alone does not explain the complex coextrusion process and the viscoelastic effect needs to be considered.

• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2005.06a
• /
• pp.210-215
• /
• 2005

A three-dimensional numerical simulation is performed to investigate on a low Reynolds number mixed convection in a horizontal rectangular channel with the upper part cooled and the lower part heated uniformly. The three-dimensional governing equations are solved using a finite volume method. For convective term, the central differencing scheme is used and for the pressure correction, the PISO algorithm is used. Solutions are obtained for A=4, Pr=0.72, 10, 909, the Reynolds number ranging from $2.1{\times}10^{-2}$ to $1.2{\times}10^{-1}$, the Rayleigh number is $3.5{\times}10^4$. It is found that vortex roll structures of mixed convection in horizontal rectangular channel can be classified into three roll structures which affected by Prandtl number and Reynolds number.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.10 no.2
• /
• pp.217-226
• /
• 1998

A numerical study has been carried out for two- and three-dimensional buoyant turbulent flow in a stairwell model. The Reynolds-averaged Navier-Stokes and energy equations are solved with the authors'own computer program. Two models by the Boussinesq approximation and the density-gradient form are used for buoyancy terms in the governing equations. Two- and three-dimensional predictions of the velocity and temperature fields are presented and the results are compared with experimental data. Comparisons have also been made in detail with two-dimensional predictions. Two-dimensional and three-dimensional simulations have predicted the overall features of the flow satisfactorily. A better agreement with experiment is achieved with three-dimensional simulations.

• Journal of Korea Water Resources Association
• /
• v.38 no.11
• /
• pp.973-984
• /
• 2005

A three-dimensional thermal hydrodynamic model is developed for carrying out unsteady simulation of the selective withdrawal of the stably stratified flow in a geometrically complex, natural reservoir The governing equations are discretized on a non-staggered grid using a second-order accurate, finite-volume scheme. The numerical model is validated by applying it to simulate three-dimensional, turbulent, stratified, shear-layer flow case. The numerical predictions appear to capture reasonably well the general shape of velocity and temperature profiles observed in the laboratory experiments, while significant overestimation of the magnitude of velocity profiles is observed in the application to the flow in a natural reservoir. The physics of selective withdrawal as emerge from the numerical simulations are also discussed.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.258-259
• /
• 2003

Flow around Savonius rotor is studied by means of the numerical simulation. Three-dimensional incompressible Navier-Stokes equations are solved numerically. Overgrid system is employed in order to enable the flow calculation of complex geometry. The basic equations in each region are solved by using the standard MAC method. The physical quantities such as the velocity and the pressure among each region are transferred through the overlapping region which is common in each region. Some numerical results of static and rotating rotor will be presented.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.42-48
• /
• 2004

Numerical simulation of viscous compressible flow in turbomachinery cascade involves many problems due to the complex geometry of blade but also flow phenomena. In the present study, numerical investigations have been performed to examine the three-dimensional flow characteristics inside the transonic linear turbine cascades using a commercial code, FLUENT. Multi-block H-type grids are applied to the high-turning turbine rotor blades and comparisons with the experimental data and the numerical results have been done. In addition, the effects of turbulence models on the prediction of the endwall flows are analyzed in the sense of the flow compressibility.

• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2005.11a
• /
• pp.29-30
• /
• 2005

A three-dimensional numerical calculation has been performed to investigate mixed convective vortex flow in rectangular channels(width/height=4) with the upper part cooled and the lower part heated uniformly. In this study, the Prandtl number was 909, the Reynolds number was varied from 0 to $9.6{\times}10^{-2}$ and the Rayleigh number from $10^3$ to $5{\times}10^4$. The governing equations were discretized using the finite volume method. From a parametric study, velocity and temperature distributions were obtained and discussed. It is found that vortex flow of mixed convection in rectangular channels can be classified into three flow patterns which depend on Reynolds and Rayleigh numbers, and the regular vortex structure disappears around Rayleigh number $5{\times}10^4$.

• Computers and Concrete
• /
• v.23 no.1
• /
• pp.69-80
• /
• 2019

Modeling approach for mesoscopic model of concrete depicting mass transportation and physicochemical reaction is important since there is growing demand for accuracy and computational efficiency of numerical simulation. Mesoscopic numerical simulation considering binder, aggregate and Interfacial Transition Zone (ITZ) generally produces huge number of DOFs, which is inapplicable for full structure. In this paper, a three-dimensional multiscale approach describing three-phase structure of concrete was discussed numerically. An effective approach generating random aggregate in polygon based on checking centroid distance was introduced. Moreover, ITZ elements were built by parallel expanding the surface of aggregates on inner side. By combining mesoscopic model including full-graded aggregate and macroscopic model, cases related to diffusivity and thickness of ITZ, volume fraction and grade of aggregate were studied regarding the consideration of multiscale compensation. Results clearly showed that larger analysis model in multiscale model expanded the diffusion space of chloride ion and decreased chloride content in front of rebar. Finally, this paper addressed some worth-noting conclusions about the chloride distribution and rebar corrosion regarding the configuration of, rebar diameter, concrete cover and exposure period.