本文擬開發一套石油洩漏污染之模擬程式。本軟體可讓工程師，當有石油洩漏時，可迅速評估其在土壤中之油污染擴散之情形，以便提供災變善後處理之依據。本計劃所探討之污染擴散數學模式，則以石油在多孔性介質(porous media)的流場來分析。

本文首先建立石油在多孔性介質中流動之三維暫態流場控制方程式，並配合適型座標轉換(body-fitted coordinate)，以有限差分法來求解。本軟體可讓工程師，當有石油洩漏時，可了解及分析在不同時間及位置時，土壤中石油之壓力場及速度場之擴散情形，以便提出對策以解決油污染擴散之問題。

In this paper, interactive soft wares were developed for the diffusion of the pollutants when the leakage of the petroleum occurs. It can let the engineers quickly evaluate the spread of pollutants and then they can take necessary measures to control the oil leakage. The mathematical models for the diffusion of the pollutants are based on the convection flows of oils in porous media.

First, the paper, three dimensional governing equations that account for the conservation of mass and momentum were derived for the convective flow from in a porous medium. The resulting governing equations were solved by the implicit finite difference method incorporated with the body-fitted coordinate transformation. Finally, interactive computer-aided design softwares were developed. The engineers can easily use this code to predict the pressure and velocity distribution in the soil at different times and locations when the petroleum leaks.

1. H. Darcy, “Les Frontaines Publiques de la Ville de Dijon”, Dalmont Paris, 1856.

2. J. C. Slattery, “Flow of Viscoelastic Fluids Through Porous Medium”, AIChE J., Vol. 13, pp. 1066-1071 , 1967.

3. S. Whitaker, “Diffusion and Dispersion in Porous Medium”, AIChE Journal, Vol. 13, pp. 420-427, 1967.

4. P. Forchheimer, “Wasserbewegung Durch Boden”, Forschtlft Ver. D. Ing. , Vol. 45, pp. 1782-1788, 1901.

5. S. Ergun, “Fluid Flow Through Packed Columns”, Chemical Engineering. Progress, Vol. 48, pp. 89-94, 1952.

6. M. Choudary, M. Propster and J. Szekely, “On The Importance of The Interial Terms in The Modeling of The Flow Misdistribution in Packed Beds” , AIChE Journal, Vol. 22, pp. 600-602 , 1976.

7. K. Vafai and C. L. Tien, “Boundary and Inertia Effects on Flow and Heat Transfer in Porous Media”, International Journal of Heat Mass Transfer, Vol. 24, pp. 195-203, 1981.

8. H. C. Brinkman, “A Calculation of The Viscous Force Exerted by a Flowing Fluid on a Dense Swarm of Particles, Applied Science, Vol. 1, pp. 27-34, 1947.

9. J. G. Georgiadis and I. Catton, “An Effective Equation Governing Convective Transport in Porous Media”, ASME Journal of Heat Transfer, Vol. 110, pp. 635-641, 1988.

10. I. A. Demirzic, A Finite Volume Method for Computation of Fluid Flow in Complex Geometries, PHD thesis, University of London, 1982.

11. J. F. Thompson, Numerical Grid Generation Foundations and Applications, North Holland, New York, 1985.

12. R. E. Smith and L. E. Eriksson, “Algebraic Grid Generation”, Computer Methods in Applied Mechanics and Engineering, Vol. 64, pp. 285-300,1987.

13. R. E. Smith and B. L. Weigel, “Analytic and Approximate Boundary Fitted Coordinate Systems for Fluid Flow Simulation”, AIAA Paper 80-0192, Pasadena, California, January 14-16, 1980.

14. J. S. Leu and J. Y. Jang, "The natural convection from a point heat source embedded in a non-Darcian porous medium", International Journal of Heat and Mass Transfer, Vol.38, No.6, pp.1097-1104, 1995.

15. J. S. Leu and J. Y. Jang, “The wall and free plumes above a horizontal line source in non-Darcian porous medium”, International Journal of Heat and Mass Transfer, Vol.37, No.13, pp.1925-1933, 1994.

16. J. Y. Jang and W. J. Chang, “Buoyancy-Induced Inclined Boundary Layer Flow in a Porous Medium Resulting from Combined Heat and Mass Buoyancy Effects”, International Communications in Heat and Mass Transfer, Vol. 15, No. 1, PP. 17-30.1988.

17. J. Y. Jang and W. J. Chang,“ The Flow and Vortex Instability of Horizontal Natural Convection in a Porous Medium Resulting from Combined Heat and Mass Buoyancy Effect”, International Journal of Heat and Mass Transfer, Vol. 31, No. 4, PP. 769-777,1988.