• Journal of Power System Engineering
• /
• v.19 no.6
• /
• pp.19-25
• /
• 2015

The unsteady-state, incompressible and three-dimensional large eddy simulation(LES) was carried out to analyze the aerodynamic performance of three-dimensional small-size axial fan(SSAF) with the different depth of bellmouth. The static pressure coefficients analyzed by LES predict a little bit larger than measurements except stall region regardless of the installation depth between SSAF and bellmouth. Moreover, static pressure efficiencies analyzed by LES show about maximum 30% at the actual operating point ranges, but measurements do not. Therefore, if the blades of conventional SSAF have some more rigidity and complete dynamic balance, the aerodynamic performance of SSAF will be some more improved. In consequence, LES shows the best prediction performance in comparison with any other Reynolds averaged Navier-Stokes(RANS) method.

• Journal of Mechanical Science and Technology
• /
• v.15 no.9
• /
• pp.1302-1310
• /
• 2001

This paper introduces modeling and simulation results for pipeline inspection gauge (PIG) with bypass flow control in natural gas pipeline. The dynamic behaviour of the PIG depends on the different pressure across its body and the bypass flow through it. The system dynamics includes: dynamics of driving gas flow behind the PIG, dynamics of expelled gas in front of the PIG, dynamics of bypass flow, and dynamics of the PIG. The bypass flow across the PIG is treated as incompressible flow with the assumption of its Mach number smaller than 0.45. The governing nonlinear hyperbolic partial differential equations for unsteady gas flows are solved by method of characteristics (MOC) with the regular rectangular grid under appropriate initial and boundary conditions. The Runge-Kuta method is used for solving the steady flow equations to get initial flow values and the dynamic equation of the PIG. The sampling time and distance are chosen under Courant-Friedrich-Lewy (CFL) restriction. The simulation is performed with a pipeline segment in the Korea Gas Corporation (KOGAS) low pressure system, Ueijungboo-Sangye line. Simulation results show us that the derived mathematical model and the proposed computational scheme are effective for estimating the position and velocity of the PIG with bypass flow under given operational conditions of pipeline.

• Journal of Korea Water Resources Association
• /
• v.50 no.1
• /
• pp.55-64
• /
• 2017

This study presents a numerical model for simulating turbidity currents using the ULTIMATE scheme. For this, the layer-averaged model is used. The model is applied to laboratory experiments, where the flume is composed of sloping and flat parts, and the characteristics of propagating turbidity currents are investigated. Due to the universal limiter of the ULTIMATE scheme, the frontal part of the turbidity currents at a sharp gradient without numerical oscillations is computed. Simulated turbidity currents propagate super-critically to the end of the flume, and internal hydraulic jumps occur at the break-in-slope after being affected by the downstream boundary. It is found that the hydraulic jumps are computed without numerical oscillations if Courant number is less than 1. In addition, factors that affect propagation velocity of turbidity currents is studied. The particle size less than $9{\mu}m$ does not affect propagation velocity but the buoyancy flux affects clearly. Finally, it is found that the numerical model computes the bed elevation change due to turbidity currents properly. Specifically, a discontinuity in the bed elevation, arisen from the hydraulic jumps and resulting difference in sediment entrainment, is observed.

• Water for future
• /
• v.27 no.4
• /
• pp.123-133
• /
• 1994

The kinematic and diffusion models using simplified momentum equations of the full dynamic equation have been frequently used for numerical flow simulations, because they have several computational advantages compared to the full dynamic model. In this paper, the more generally acceptable application ranges of the kinematic and diffusion finite difference models were investigated based on three major parameters, which are channel bed slopes So, dimensionless depth increasing numbers Gw at upstream boundary and Froude numbers Fr. The applicable ranges were obtained by comparing the relative magnitudes of the local acceleration, convective acceleration, pressure, gravity and friction terms in the full dynamic equation. In the simulations, a Courant number of 0.5 was used and the channel bed slopes were changed from 0.00001 to 0.05. Also, Froude numbers of 0.1, 0.5 and 0.9 were employed. In this paper, it is indicated that the applicable ranges of kinematic models are increased with increasing of Froude numbers. However, the applicable ranges of diffusion models are decreased with increasing of Froude numbers. Finally, 9 figures were proposed as a guideline in the application of kinematic and diffusion finite difference models based upon the allowable deviation compared to the full dynamic model. With applying the proposed criteria, it is expected that the flow simulations in the channels, streams or rivers are more efficiently achieved.

• Journal of Soil and Groundwater Environment
• /
• v.11 no.1
• /
• pp.37-44
• /
• 2006

This study summarizes advantages and disadvantages of numerical methods and compares ELLAM and LEZOOMPC to develop an efficient numerical modeling technique on contaminant transport. Eulerian-Lagrangian method and Eulerian method are commonly used numerical techniques. However Eulerian-Lagrangian method does not conserve mass globally and fails to treat boundary in a straightforward manner. Also, Eulerian method has restrictions on the size of Courant number and mesh Peclet number because of time truncation error. ELLAM (Eulerian Lagrangian Localized Adjoint Method) which has been popularly used for past 10 years in numerical modeling, is known for overcoming these numerical problems of Eulerian-Lagrangian method and Eulerian method. However, this study investigates advantages and disadvantages of ELLAM and suggests a change for the better. To figure out the disadvantages of ELLAM, the results of ELLAM, LEZOOMPC (Lagrangian-Eulerian ZOOMing Peak and valley Capturing), and visual MODFLOW are compared for four examples having different mesh Peclet numbers. The result of ELLAM generates numerical oscillation at infinite of mesh Peclet number, but that of LEZOOMPC yields accurate simulations. The simulation results suggest that the numerical error of ELLAM could be alleviated by adopting some schemes in LEZOOMPC. In other words, the numerical model which combines ELLAM with backward particle tracking, forward particle tracking, adaptively local zooming, and peak/valley capturing of LEZOOMPC can be developed for not only overcoming the numerical error of ELLAM, but also keeping the numerical advantage of ELLAM.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.3 no.1
• /
• pp.13-24
• /
• 1983

A mathematical model for the two-dimensional unsteady flow was developed by introducing Abbott's implicit finite difference operator and double sweep algorithm, which could be applied to simulate the respose of a harbor against the intrusion of long waves through the entrance connected to open sea. In order to improve its accuracy corresponding to the field phenomena, bottom resistance, Coriolis force, wind effect terms were included and wave direction and radiating effect was considered. The result of seiche test was always stable and the amplitude was accurate. Some phase shift was occured, but it could be reduced by using small values of Courant number and many points per a wave length as well. A comparision with the Ippen and Goda's theoritical and hydraulic experimental works was fulfilled.

• Journal of ILASS-Korea
• /
• v.26 no.3
• /
• pp.127-134
• /
• 2021

Fire sprinkler initial spray was analyzed by Large eddy simulation (LES) and Volume of Fluid (VoF) integrated method. The IsoAdvector geometric VoF was used to identify the liquid-gas interface clearly even with the large Courant-Friedrichs-Lewy number. To reduce the computational costs, sector meshes and Adaptive Mesh Refinement up to level 3 were used. Base mesh size was 1 mm, which is roughly equivalent to the initial sprinkler droplet. Top surface radius of boss and deflector size were modified to investigate the effects of sprinkler head design on primary breakup process. When top surface radius of boss was increased, vertical liquid sheet was formed. This phenomenon reduced the sheet breakup radius, sheet thickness and velocity. Due to reduced liquid sheet thickness, a large amount of ligaments was created from the liquid sheet. As a result, there was a dramatic decrease in volume per surface area, indicating an increase in breakup process. Spray pattern viewed in radial direction also changed when top surface radius of boss increased. When top surface radius of boss was increased, a T-shaped pattern was observed while a V-shaped pattern was observed in all other cases. When the deflector size increases, the spray pattern remains V-shaped, even if the top surface radius of boss increased. Further studies on promoting atomization of the water supplied to the lower part of the sprinkler head in the T-shape pattern should be conducted.

• Atmosphere
• /
• v.33 no.4
• /
• pp.331-341
• /
• 2023

A numerical forecasting models usually predict future states by performing time integration considering fixed static time-steps. A time-step that is too long can cause model instability and failure of forecast simulation, and a time-step that is too short can cause unnecessary time integration calculations. Thus, in numerical models, the time-step size can be determined by the CFL (Courant-Friedrichs-Lewy)-condition, and this condition acts as a necessary condition for finding a numerical solution. A static time-step is defined as using the same fixed time-step for time integration. On the other hand, applying a different time-step for each integration while guaranteeing the stability of the solution in time advancement is called an adaptive time-step. The adaptive time-step algorithm is a method of presenting the maximum usable time-step suitable for each integration based on the CFL-condition for the adaptive time-step. In this paper, the adaptive time-step algorithm is applied for the Korean Integrated Model (KIM) to determine suitable parameters used for the adaptive time-step algorithm through the monthly verifications of 10-day simulations (during January and July 2017) at about 12 km resolution. By comparing the numerical results obtained by applying the 25 second static time-step to KIM in Supercomputer 5 (Nurion), it shows similar results in terms of forecast quality, presents the maximum available time-step for each integration, and improves the calculation efficiency by reducing the number of total time integrations by 19%.

• Journal of Korea Soil Environment Society
• /
• v.5 no.1
• /
• pp.3-12
• /
• 2000

The multi-region model, to describe preferential flow, is an equation representing solute transport in soils by dividing soil into numerous pore groups and using the hydraulic properties of the soil. As the model partial differential equation (PDE) is solved numerically with finite difference methods. a modified equivalent partial differential equation(MEPDE) of the partial differential equation of the multi-region model is derived to analyze the accuracy and consistency of the solution of the model PDE and the Von Neumann method is used to analyze the stability of the finite difference scheme. The evaluation obtained from the MEPDE indicated that the finite difference scheme was found to be consistent with the model PDE and had the second order accuracy The stability analysis is performed to analyze the model PDE with the amplification ratio and the phase lag using the Von Neumann method. The amplification ratio of the finite difference scheme gave non-dissipative results with various Peclet numbers and yielded the most high values as the Peclet number was one. The phase lag showed that the frequency component of the finite difference scheme lagged the true solution. From the result of the stability analysis for the model PDE, it is analyzed that the model domain should be discretized in the range of Pe < 1.0 and Cr < 2.0 to obtain the more accurate solution.