• Journal of Korea Water Resources Association
• /
• v.32 no.6
• /
• pp.607-616
• /
• 1999

A fractional step finite difference model for the longitudinal dispersion of nonconservative contaminants is developed. It is based on splitting the longitudinal dispersion equation into a set of three equations each to be solved over a one-third time step. The fourth-order Holly-Preissmann scheme, an analytic solution, and the Crank-Nicholson scheme are used to solve the equations for the pure advection, the first-order decay, and the diffusion, respectively. To test the model, it is applied to simulate the longitudinal dispersion of continuous source released into a nonuniform flow field as well as the dispersion of an instantaneous source in a uniform flow field. The results are compared with the exact solution and those computed by an existing model. Compared to the existing model which uses Euler method for the first-order decay equation, the present model yield more accurate results as the decay coefficient increases.

• Water for future
• /
• v.29 no.4
• /
• pp.187-198
• /
• 1996

A numerical model for unsteady dispersion of horizontal line source in turbulent shear flow is developed. A fractional step finite difference method is used which splits the unsteady two-dimensional advective diffusion equation into the longitudinal advection and the vertical diffusion equations, and solves them alternately for half time intervals by the Holly-Preissmann scheme and the Crank-Nicholson scheme, respectively. The developed numerical model is verified using a semi-analytic solution for steady dispersion in turbulent shear flow. Dispersion of an instantaneous plane source in turbulent shear flow is analyzed using the model. The degree of mixing at the same dimensionless time is almost the same regardless of the friction factor, and the travel distance required to reach a certain degree of mixing is inversely proportional to the square root of the friction factor.

• Journal of Korea Water Resources Association
• /
• v.45 no.4
• /
• pp.349-360
• /
• 2012

A horizontally curvilinear non-hydrostatic free surface model that was applicable to three-dimensional viscous flows was developed. The proposed model employed a top-layer equation to close kinematic free-surface boundary condition, and an isotropic k-${\varepsilon}$ model to close turbulence viscosity in the Reynolds averaged Navier-Stokes equation. The model solved the governing equations with a fractional step method, which solved intermediate velocities in the advection-diffusion step, and corrects these provisional velocities by accounting for source terms including pressure gradient and gravity acceleration. Numerical applications were implemented to the wind-driven currents in a two-dimensional closed basin, the flow in a steep-sided trench, and the flow in a strongly-curved channel accounting for secondary current by the centrifugal force. Through the numerical simulations, the model showed its capability that were in good agreement with experimental data with respect to free surface elevation, velocity, and turbulence characteristics.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.2
• /
• pp.181-188
• /
• 2012

The present study on nitrogen-diluted non-premixed counterflow flames with finite burner diameters experimentally investigates the important role of the outer edge flame in flame extinction. Flame stability diagrams mapping the flame extinction response of nitrogen-diluted non-premixed counterflow flames to varying global strain rates in terms of the burner diameter, burner gap, and velocity ratio are explored. There exists a critical nitrogen mole fraction beyond which the flame cannot be sustained, and also the curves of the critical nitrogen mole fraction versus the global strain rate have C-shapes in terms of burner diameter, burner gap, and velocity ratio. In flames with sufficiently high strain rates, the curves of the critical nitrogen mole fractions versus global strain rate collapse into one curve, and the flames can have the 1-D flame response of typical diffusion flames. Three flame extinction modes are identified: flame extinctions through the shrinkage of the outer edge flame with and without an oscillation of the outer edge flame prior to the extinction and flame extinction through a flame hole at the flame center. The measured flame surface temperature and a numerical evaluation of the fractional contribution of each term in the energy equation show that the radial conductive heat loss at the flame edge destabilizes the outer edge flame, and the conductive and convection heat addition to the outer edge from the trailing diffusion flame stabilizes the outer edge flame. The radial conductive heat loss at the flame edge is the dominant extinction mechanism acting through the shrinkage of the outer edge flame.