• Journal of Environmental Science International
• /
• v.11 no.12
• /
• pp.1321-1326
• /
• 2002

A numerical model is investigated to predict a behavior of the gaseous volatile organic compounds and a subsurface contamination caused by them in the unsaturated zone. Two dimensional advective-dispersion equation caused by a density difference and two dimensional diffusion equation are computed by a finite difference method in the numerical model. A laboratory experiment is also carried out to compare the results of the numerical model. The dimensions of the experimental plume are 1.2m in length, 0.5m in height, and 0.05m in thickness. In comparing the result of 2 methods used in the numerical model with the one of the experiment respectively, the one of the advective-dispersion equation shows better than the one the diffusion equation.

• Water for future
• /
• v.23 no.1
• /
• pp.119-127
• /
• 1990

The concentration of cohesive suspended sediment is determined by the circulation of water and the material dispersion. The equations of the two-dimensional, depth-integrated dispersive transport are the Reynolds equation, continuity equation, and advection-dispersion equation based on the Fick's law. A finite difference method has been applied to two models of circulation and dispersion transport. The circulation model is solved by the explicit scheme and the dispersion transport model is solved by multi-operational scheme. It is investigated wheter advective terms are included when the equation of circulation is applied to the model. For advection-dispersion equation, it was also investigated about variations of suspended sediment concentration with respect to the critical shear stresses.

• 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 the Korean Society of Safety
• /
• v.37 no.2
• /
• pp.1-9
• /
• 2022

Most factories deal with toxic or flammable chemicals in their industrial processes. These hazardous substances pose a risk of leakage due to accidents, such as fire and explosion. In the event of chemical release, massive casualties and property damage can result; hence, quantitative risk prediction and assessment are necessary. Several methods are available for evaluating chemical dispersion in the atmosphere, and most analyses are considered neutral in dispersion models and under far-field wind condition. The foregoing assumption renders a model valid only after a considerable time has elapsed from the moment chemicals are released or dispersed from a source. Hence, an initial dispersion model is required to assess risk quantitatively and predict the extent of damage because the most dangerous locations are those near a leak source. In this study, the dispersion model for initial consequence analysis was developed with three-dimensional unsteady advective diffusion equation. In this expression, instantaneous leakage is assumed as a puff, and wind velocity is considered as a coordinate transform in the solution. To minimize the buoyant force, ethane is used as leaked fuel, and two different diffusion coefficients are introduced. The calculated concentration field with a molecular diffusion coefficient shows a moving circular iso-line in the horizontal plane. The maximum concentration decreases as time progresses and distance increases. In the case of using a coefficient for turbulent diffusion, the dispersion along the wind velocity direction is enhanced, and an elliptic iso-contour line is found. The result yielded by a widely used commercial program, ALOHA, was compared with the end point of the lower explosion limit. In the future, we plan to build a more accurate and general initial risk assessment model by considering the turbulence diffusion and buoyancy effect on dispersion.