Computational & Technology Resources
an online resource for computational,
engineering & technology publications 

CivilComp Proceedings
ISSN 17593433 CCP: 80
PROCEEDINGS OF THE FOURTH INTERNATIONAL CONFERENCE ON ENGINEERING COMPUTATIONAL TECHNOLOGY Edited by: B.H.V. Topping and C.A. Mota Soares
Paper 6
Simulation of Thermal Outflow of Power Plants in the Huelva Estuary A. German+, M. Espino+, M. Maidana+ and J. Blasco*
+Laboratory of Maritime Engineering,
A. German, M. Espino, M. Maidana, J. Blasco, "Simulation of Thermal Outflow of Power Plants in the Huelva Estuary", in B.H.V. Topping, C.A. Mota Soares, (Editors), "Proceedings of the Fourth International Conference on Engineering Computational Technology", CivilComp Press, Stirlingshire, UK, Paper 6, 2004. doi:10.4203/ccp.80.6
Keywords: pollutant dispersion, convectiondiffusion equation, finite element methods, explicit characteristicbased method, thermal outflow, Huelva estuary.
Summary
A simple explicit, characteristicbased finite element method [1] for the
numerical simulation of the dispersion of thermal outlow in the Huelva estuary is
studied in this presentation. The derivation involves a local Taylor expansion of the
convectiondiffusion equation. The numerical model used, called SEASCAPE [2,3]
which stands for Simple Explicit GAlerkin Procedure to Solve
ConvectionDiffusion EquAtion in the OPen SEa,
was originally designed to simulate the
dispersion of the of pollutants in the open sea [4] and has now been reconfigured to
simulate the dispersion of thermal outflow of some power plants which may be
constructed in the Huelva Estuary in the near future. Numerical results obtained are
presented.
In term of accuracy, stability and convergence properties, the simple explicit CharacteristicGalerkin method was chosen to solve the convectiondiffusion equation. This method has the required properties which include the introduction of the stabilizing term or artificial viscosity that stabilises the convective oscillations. The method starts from a particular time discretization of the differential equation before obtaining the weak form. Time discretization is based on Taylor expansion of the convectiondiffusion equation. The dispersion of contaminants in bodies of water is intrinsically three dimensional but due to high computational cost of threedimensional discretization, we opted to use a quasi3D approximation. The basic idea of the quasi3D model is the combination of the conventional 2DH model and a local 1DV model that describes the vertical distribution of the concentration of the contaminants through its decomposition in base function and profile. The standard Galerkin method of weighted residuals is applied on the time discretized equations where vertical variation of the unknown variables is approximated by a linear combination of known basis function attached to each nodal value or degree of freedom. The magnitude of the variable is discretized horizontally to eliminate the x and y dependence. The horizontal variations of the variables of the problem are interpolated using the isoparametric finite element basis function sufficient to interpolate the diffusive terms. Galerkin method of weighted residuals is applied and integration by parts is employed to reduce the second order partial derivatives to first order. The algorithm is implemented in the numerical model named SEASCAPE to simulate the dispersion of thermal outflow of some power plants that might be constructed in the Huelva estuary. The computational grid used has 3200 quadrilateral elements and 3531 nodes. Data for the mean flow and numerical propagation of tides required to simulate the dynamics of the estuary were culled from previously formulated numerical models such as NAUTILUS [5,6] and MAREAS [7]. The value of the coefficient of diffusivity K_{H} was set at 10m^{2}/s [8]. A homogenous Neumann boundary condition is prescribed on the mouth of Tinto river and Dirichlet boundary conditions are imposed on the locations of the wastewater discharges. Simulation results, when all three power plants are discharging water into the estuary, show that substantial increase in water temperature can be observed at the end of the simulation period. Contaminant dispersion does not reach a stationary condition even after the simulation period has elapsed. The simulation period must be prolonged to see how the thermal outflow will evolve. The numerical results presented show that the numerical model provides a comprehensive and efficient tool to simulate contaminant dispersion. References
purchase the fulltext of this paper (price £20)
go to the previous paper 
