• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.4
• /
• pp.865-875
• /
• 1988

A two-dimensional stepped wall jet was numerically investigated by applying three different models : One is the standard k-.epsilon. and the other is the modified k-.epsilon. model which takes account of the streamline curvature effect by modifying the Reynolds shear stress and a source term in the dissipation equation, and a third is curvature dependent third-order correlation model. In order to test the influences of the numerical result, both the upwind scheme and the skew-upwind scheme were sued for the computations. By comparing the numerical results with available experiments, it was found that the modified k-.epsilon. model gives best overall prediction accuracy only when the numerical diffusion is eliminated by using the skew-upwind scheme. The numerical scheme was found to have more pronounced effect on the accuracy of the turbulence computation than the turbulence models.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.409-413
• /
• 2005

Even though the relative importance of length scale of flow system allow us to simplify three dimensional flow problem to one or two dimensional representation, many systems still require three dimensional analysis. The objective of this study is to develop an efficient and accurate finite element model for analyzing and predicting three dimensional flow features in natural rivers and to offend to model spreading of pollutants and transport of sediments in the future. Firstly, three dimensional Reynolds averaged Navier-Stokes equations with the hydrostatic pressure assumption in generalized curvilinear coordinates were combined with the kinematic free-surface condition. Secondly. to simulate realistic high Reynolds number flow, the model employed the Streamline Upwind/Petrov-Galerkin(SU/PG) scheme as a weighting function for the finite element method in conjunction with an appropriate turbulence model(Smagorinsky scheme for the horizontal plain and Mellor-Yamada scheme for the vertical direction). Several tests is performed for the purpose of validation and verification of the developed model. A simple rectangular channel, 5-shaped and U-shaped channel are used for tests and comparisons are made with RMA-10 model. Runs for each case is converged stably without a oscillation and calculated water-surface deformation, longitudinal and transversal velocities, and velocity vector fields are in good agreement with the results of RMA-10 model.

• Journal of Korea Water Resources Association
• /
• v.38 no.5 s.154
• /
• pp.333-343
• /
• 2005

The upwind model is well known to simulate shockwaves, but it is rarely applied to natural rivers because of problems caused by the source terms. Although several methods have been developed to deal with the source terms, none of them has been applied to natural rivers. This paper deals with application of the upwind model to the natural river. An implicit upwind model is applied to a hypothetical irregular channel and a natural river with highly irregular bed, width, and hydraulic structures. Different types of the flows including steady-state flow, flood wave, dam-break wave, and bore are simulated to test accuracy and applicability of the implicit upwind model. It is proved that the model can simulate various types of flows in natural rivers with high accuracy and robustness.

• Journal of Environmental Health Sciences
• /
• v.31 no.6
• /
• pp.489-497
• /
• 2005

A numerical method for the mathematical water modeling in 2-dimensional flow has been developed. The model based on a split operator technique, in which, the advection term is calculated using the upwind scheme. The diffusion term is one- dimensionalized and calculated using Crank-Nicholson's implicit finite difference scheme to reduce the numerical errors from large time steps and variable spacings. It also provides a relatively simple and economic method for more accurate simulation of pollutant dispersion. Water depths and flow velocities in the Boreyong reservoir during the normal water periods were predicted by numerical experiments with a 2-dimensional flow model so as to provide current field data for the study of advection and diffusion of pollutants. Developed 2-dimensional water quality model is applied to Boreyong reservoir to simulate a spatial and periodical changes of water quality.

• 한국방재학회:학술대회논문집
• /
• 2008.02a
• /
• pp.233-236
• /
• 2008

To analyze the saline intrusion in the place, such as an estuary, the three-dimensional numerical model is developed. In this study, the advection terms of the governing equations are discretized by upwind scheme. By using an explicit scheme for the longitudinal direction and an implicit scheme for the vertical direction, the numerical model is free from the restriction of temporal step size caused by a relatively small grid ratio. The equation of state is used to consider the density, and the scalar transport equation for salinity is employed the third order TVD to scheme to prevent unphysical oscillation near discontinuity. In order to verify saline intrusion, the numerical model is conducted to compare the previous model in the lock exchange. The present model generally show a good agreement with the previous one.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2006.05a
• /
• pp.425-429
• /
• 2006

The understanding and prediction of the behavior of flow in open channels are important to the solution of a wide variety of practical flow problems in water resources engineering. Recently, frequent drought has increased the necessity of an effective water resources control and management of river flows for reserving instream flow. The objective of this study is to develop an efficient and accurate finite element model based on Streamline Upwind/Petrov-Galerkin(SU/PG) scheme for analyzing and predicting two dimensional flow features in complex natural rivers. Several tests were performed in developed all elements(4-Node, 6-Node, 8-Node elements) for the purpose of validation and verification of the developed model. The U-shaped channel of flow and natural river of flow were performed for tests. The results were compared with these of laboratory experiments and RMA-2 model. Such results showed that solutions of high order elements were better accurate and improved than those of linear elements. Also, the suggested model displayed reasonable velocity distribution compare to RMA-2 model in meandering domain for application of natural river flow. Accordingly, the developed finite element model is feasible and produces reliable results for simulation of two dimensional natural river flow. Also, One contribution of this study is to present that results can lead to significant gain in analyzing the accurate flow behavior associated with hydraulic structure such as weir and water intake station and flow of chute and pool.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.528-531
• /
• 2004

For a simple one-dimensional ignition problem a mathematical model is described to investigate the difficulties in numerical simulations. Some computation results are obtained and comparison is made with analytical solution. Discussions are made on topics such as 1) coordinate transformation, 2) gas-phase and solid-phase analysis; (divergence form of the governing system, a finite-volume discretization, implicit time integration, upwind split flux, spatial accuracy improvement are described. Mass, reagent mass, and energy conservations are solved.), and 3) method to determine quantities on the burning surface (matching). Results obtained for small values of the non-dimensional pressure show a steady-combustion and good agreement with the analytical solution. Numerical instability appeared for larger values of the pressure, discussion on the cause of the problem is made. This effort is a part of a study of flame spread phenomena on solid propellant surface.

• Journal of Mechanical Science and Technology
• /
• v.14 no.9
• /
• pp.985-996
• /
• 2000

The issue of quench is related to safety operation of large-scale superconducting magnet system fabricated by cable-in-conduit conductor. A numerical method is presented to simulate the thermal hydraulic quench characteristics in the superconducting Tokamak magnet system, One-dimensional fluid dynamic equations for supercritical helium and the equation of heat conduction for the conduit are used to describe the thermal hydraulic characteristics in the cable-in-conduit conductor. The high heat transfer approximation between supercritical helium and superconducting strands is taken into account due to strong heating induced flow of supercritical helium. The fully implicit time integration of upwind scheme for finite volume method is utilized to discretize the equations on the staggered mesh. The scheme of a new adaptive mesh is proposed for the moving boundary problem and the time term is discretized by the-implicit scheme. It remarkably reduces the CPU time by local linearization of coefficient and the compressible storage of the large sparse matrix of discretized equations. The discretized equations are solved by the IMSL. The numerical implement is discussed in detail. The validation of this method is demonstrated by comparison of the numerical results with those of the SARUMAN and the QUENCHER and experimental measurements.

• Journal of Korea Water Resources Association
• /
• v.36 no.5
• /
• pp.777-785
• /
• 2003

A numerical model for the solution of two-dimensional free surface flow is developed on unstructured grid. By using fractional step method, the two-dimensional shallow water equations (SWE) are treated as two one-dimensional problems. Thus, it is possible to simulate computational hydraulic problems with higher computational efficiency. The one-dimensional problems are solved using upwind TVD version of second-order Weighted Averaged Flux (WAF) scheme with HLLC approximate Riemann solver. The numerical oscillations which are common with second-order numerical scheme are controlled by exploiting WAF flux limiter, Some idealized test problems are solved using this model and very accurate and stable solutions are obtained. It can be concluded as an efficient implement for the computation of SWE including dam break problems that concerning discontinuities, subcritical and supercritical flows and complex domain.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.15 no.6
• /
• pp.38-46
• /
• 2011

In order to simulate the ignition-gas injection in the interior ballistics, a two-dimensional analytic code for two-phase flows has been developed. The Eulerian-Lagrangian approach and the low-dissipation simple high-resolution upwind scheme(LSHUS) have been adopted in the numerical code for the propellant combustion of the gun propelling charges. The ghost-cell extrapolation method has been used for the moving boundary in the chamber with the projectile movement. The calculation results of the developed code have been compared and verified through those of the dimensionless IBHVG2 code and the previous one-dimensional code. In comparison with the two-phase flows according to the drag models, the numerical analysis of the muzzle velocity has been affected by the drag model.