• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.1B
• /
• pp.73-81
• /
• 2009

This paper presents a numerical investigation of the impact of the secondary currents on solute transport in open-channel flows. The RANS model with Reynolds stress model is used for flow modeling, and the GGDH(generalized gradient diffusion hypothesis) model is used to close the scalar transport equation. Using the developed model, the impact of secondary currents on solute transport in open channel flows over smooth-rough strip is investigated. Through numerical experiments, the secondary currents are found to affect the solute spreading, leading a movement of the position of the peak concentration and a skewed distribution of solute concentration. Due to the lateral flow of secondary currents near the free surface, the concentration at the rough strip is found to be larger than that at the smooth strip bed. The solute at the rough strip is more rapidly transported than smooth bed. A magnitude analysis of the solute transport rate in scalar transport equation is also carried out to investigate the effect of secondary currents and scalar flux on the concentration distribution.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.24 no.4
• /
• pp.277-286
• /
• 2012

In the numerical simulation of wave fields using a multi-phase flow model that considers simultaneous flows of materials with different states such as gas, liquid and solid, there is need of an accurate representation of the interface separating the fluids. We adopted an algebraic interface capturing method called tangent of hyperbola for interface-capturing(THINC) method for the capture of the free-surface in computations of multi-phase flow simulations instead of geometrical-type methods such a volume of fluid(VOF) method. The THINC method uses a hyperbolic tangent functions to represent the surface, and compute the numerical flux for the fluid fraction functions. One of the remarkable advantages of THINC method is its easy applicability to incorporate various numerical codes based on Navier-Stokes solver because it does not require the extra geometric reconstruction needed in most of VOF-type methods. Several tests were carried out in order to investigate the advection of interfaces and to verify the applicability of the THINC method to wave fields based on the one-field model for immiscible two-phase flows (TWOPM). The numerical results revealed that the THINC method is able to track the interface between air and water separating the fluids although its algorithm is fairly simple.

• The Transactions of the Korean Institute of Electrical Engineers C
• /
• v.54 no.7
• /
• pp.319-325
• /
• 2005

In order to analyze ground potential rise of grounding system installed in buildings, the hemispherical grounding simulation system has been designed and fabricated as substantial and economical measures. Ground potential rise(GPR) has been measured and analyzed for shapes of grounding electrode using the system in real time. The system is apparatus to have a free reduced scale for conductor size and laying depth of a full scale grounding system and is constructed so that a shape of equipotential surface is nearly identified a free reduced scale with a real scale when a current flows through grounding electrode. The system was composed of a hemispherical water tank, AC Power supply, a movable potentiometer, and test grounding electrodes. The test grounding electrodes were fabricated through reducing grounding electrode installed in real buildings such as rod type, mesh grid type. When a mesh grid type was associated with a rod type, GPR was the lowest value. The proposed results would be applicable to evaluate GPR in the grounding systems. and the analytical data can be used 0 stabilize the electrical installations and prevent the electrical disasters.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1765-1770
• /
• 2003

This paper describes a study of the spin-up of a free-surface fluid in a rectangular container in which an internal cylindrical obstacle is mounted. Experiments and numerical analysis have been carried out for a variety of obstacle position. Increase in the speed of background rotation and near wall position of cylindrical obstacle results in the complex flow structures. Numerical and experimental results agree well with each other and the Ekman-pumping model is also applied to this flow.

• Journal of Advanced Research in Ocean Engineering
• /
• v.4 no.3
• /
• pp.135-145
• /
• 2018

This paper deals with the added resistance of a ship in waves using computational fluid dynamics (CFD). The ship added resistance is one of the key considerations in the design of energy-efficient ship. In this study, the added resistance of a LNG carrier in head waves is computed using a CFD code to consider the nonlinearity and the viscous effects. The unsteady Reynolds Averaged Navier-Stokes equation (RANS) is numerically solved and the volume of fluid (VOF) approach is used to simulate the free surface flows. The length of incident wave varies from half the ship length to twice the ship length. To investigate the nonlinearity effect, both the linear wave condition and the nonlinear wave condition are considered. The heave and pitch motions are calculated along with the added resistance, and the wave contours are obtained. Grid convergence test is conducted thoroughly to achieve the converged motion and resistance values. The calculated results are compared and validated with experimental data.

• Journal of Ship and Ocean Technology
• /
• v.4 no.3
• /
• pp.21-32
• /
• 2000

In the present paper, the sloshing flow in the liquid holds of a large tanker is simulated using a numerical method. In the fluid domain, the three-dimensional Navier-Stokes equation with free surface is solved using a finite difference method, and the realistic shapes of multi holds are modeled including the internal members. The time-history of the tank motion is obtained using a time-domain program for ship motion. In order to computer the impulsive pressures on internal structures, a concept of buffer zone is adopted near the tank ceiling during impact occurrence. This study demonstrates that the global fluid motion in the multi liquid holds of ships and FPSO's can be simulated using the numerical method and the corresponding local pressure can be predicted with reasonable accuracy.

• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.1
• /
• pp.23-30
• /
• 2004

In the Lagrangian vortex particle method based on the vorticity-velocity formulation for solving the incompressible Navier-Stokes equations, a numerical scheme for calculating pressure fields is presented. Implementation of the numerical method is directly connected with the well-established surface panel methods, just by dealing with the dynamic coupling among vorticity field. Assuming the vorticity and the velocity fields are to be calculated in time domain analysis, the pressure calculation for a complete set of solution at present time step is performed in a similar way to the one used in the Eulerian description. For a validation of the present method, we illustrate the early development of the viscous flow about an impulsive started circular cylinder for Reynolds number 550. The comparative study with the Eulerian finite Volume method provides an extensive understanding and application of the mesh-free Lagrangian vortex methods for numerical simulation of viscous flows around arbitrary bodies of general shape.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.6
• /
• pp.675-682
• /
• 2010

The flood impact pressure acting on a vertical wall resulting from a dam-breaking problem is simulated using a navier-Stokes(N-S) solver. The N-S solver uses Eulerian Finite Volume Method(FVM) along with Volume Of Fluid(VOF) method for 2-D incompressible free surface flows. A Split Lagrangian Advection(SLA) scheme for VOF method is implemented in this paper. The SLA scheme is developed based on an algorithm of Piecewise Linear Interface Calculation(PLIC). The coupling between the continuity and momentum equations is affected by using a well-known Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) algorithm. Several two-dimensional numerical simulations of the dam-breaking problem are presented to validate the accuracy and demonstrate the capability of the present algorithm. The significance of the time step and grid resolution are also discussed. The computational results are compared with experimental data and with computations by other numerical methods. The results showed a favorable agreement of water impact pressure as well as the global fluid motion.

• The Bulletin of The Korean Astronomical Society
• /
• v.42 no.1
• /
• pp.38.1-38.1
• /
• 2017

We report the results of three-dimensional radiation hydrodynamic simulations of star cluster formation in turbulent molecular clouds, with primary attention to how stellar radiation feedback controls the lifetime and net star formation efficiency (SFE) of their natal clouds. We examine the combined effects of photoionization and radiation pressure for a wide range of cloud masses (10^4 - 10^6 Msun) and radii (2 - 80 pc). In all simulations, stars form in densest regions of filaments until feedback becomes strong enough to clear the remaining gas out of the system. We find that the SFE is primarily a function of the initial cloud surface density, Sigma, (SFE increasing from ~7% to ~50% as Sigma increases from ~30 Msun/pc^2 to ~10^3 Msun/pc^2), with weak dependence on the initial cloud mass. Control runs with the same initial conditions but without either radiation pressure or photoionization show that photoionization is the dominant feedback mechanism for clouds typical in normal disk galaxies, while they are equally important for more dense, compact clouds. For low-Sigma clouds, more than 80% of the initial cloud mass is lost by photoevaporation flows off the surface of dense clumps. The cloud becomes unbound within ~0.5-2.5 initial free-fall times after the first star-formation event, implying that cloud dispersal is rapid once massive star formation takes place. We briefly discuss implications and limitations of our work in relation to observations.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.10
• /
• pp.1365-1372
• /
• 2001

In the planar flow casting(PFC) process, the conditions of the melt puddle between nozzle and rotating wheel affect significantly the quality and dimensional uniformity of the downstream ribbon. For stable puddle formation, the nozzle is placed very close to the quenching wheel, so the surface-tension and wall-adhesion forces have an important effect upon the fluid flow.\`In this study the planar flow casting process has been mode]ed using the VOF method for free surface tracking. The transient puddle formation from the present analysis shows good agreements with the previous experimental results. Furthermore, the variation of melt temperature and the corresponding cooling rate of the melt have been examined. The present results also show how the melt puddle can be farmed on the rotating substrate, how the melt flows within the puddle, and how the changes of the process variables affect the puddle formation and its corresponding fluid flow and heat transfer behavior.