On Some Developments and Evaluation of an Eulerian-Lagrangian Method for the Transport Equation

The modelling of typical engineering problems such as water-jet cooling of hot-rolled steel strip products in industry, directly involves the solution of a transport (advection-diffusion) equation for the cooling characteristics of the strip. The non-linear nature of the heat conduction involved, aggravates the diffculty of the problem. Traditional Finite Difference techniques for the solution of this advection dominated transport equation incur severe Courant number stability restrictions as well as instabilities in the presence of temperature discontinuities. Eulerian-Lagrangian Methods (ELM's) solve the transport equation in Lagrangian form `along' backward characteristics effec- tively decoupling the advection and diffusion terms but retaining the convenience of fixed computational grids. Typical interpolation methods used to obtain the values at the feet of characteristic lines lead to spurious oscillations, numerical diffusion, peak clipping and phase errors. Through the use of `peak tracking', by the forward-tracking of Eulerian nodal points, this paper attempts to alleviate these errors. A comparison of 1-D benchmark tests from the Convection-Diffusion Forum as well as appropriate error measures, are shown to produce appreciable improvements over the standard methods for a range of timesteps, very large Peclet numbers and Courant numbers in excess of one

Numerical methods and tangible interfaces for pollutant dispersion simulation

Dissertação apresentada para obtenção do Grau de Doutor em Engenharia do Ambiente, pela Universidade Nova de Lisboa, Faculdade de Ciências e TecnologiaThe first main objective of this thesis is to reduce numerical errors in advection-diffusion modelling. This is accomplished by presenting DisPar methods, a class of numerical schemes for advection-diffusion or transport problems, based on a particle displacement distribution for Markov processes. The development and analyses of explicit and implicit DisPar formulations applied to one and two dimensional uniform grids are presented. The first explicit method, called DisPar-1, is based on the development of a discrete probability distribution for a particle displacement, whose numerical values are evaluated by analysing average and variance. These two statistical parameters depend on the physical conditions (velocity, dispersion coefficients and flows). The second explicit method,DisPar-k, is an extension of the previous one and it is developed for one and two dimensions. Besides average and variance, this method is also based on a specific number of particle displacement moments. These moments are obtained by the relation between the advection-diffusion and the Fokker-Planck equation, assuming a Gaussian distribution for the particle displacement distribution. The number of particle displacement moments directly affects the spatial accuracy of the method, and it is possible to achieve good results for pure-advection situations. The comparison with other methods showed that the main DisPar disadvantage is the presence of oscillations in the vicinity of step concentration profiles. However, the models that avoid those oscillations generally require complex and expensive computational techniques, and do not perform so well as DisPar in Gaussian plume transport. The application of the 2-D DisPar to the Tagus estuary demonstrates the model capacity of representing mass transport under complex flows. Finally, an implicit version of DisPar is also developed and tested in linear conditions, and similar results were obtained in terms of truncation error and particle transport methods. The second main objective of this thesis, to contribute to modelling cost reduction, is accomplished by presenting TangiTable, a tangible interface for pollutant dispersion simulation composed by a personal computer, a camera, a video projector and a table. In this system, a virtual environment is projected on the table, where the users place objects representing infrastructures that affect the water of an existent river and the air quality. The environment and the pollution dispersion along the river are then projected on the table. TangiTable usability was tested in a public exhibition and the feedback was very positive. Future uses include public participation and collaborative work applications.Fundação para a Ciência e Tecnologia - scholarship contract BD/5064/2001 and the research contract MGS/33998/99-0

Modelling the Eddy Pattern in the harbour of Zeebrugge

Hydrodynamic

Integrated 2D-3D free surface hydro-environmental modelling

An integrated horizontally two- and fully three-dimensional numerical model system has been developed based on a combined unstructured and σ-coordinate grid to simulate the flow and water quality process in large water bodies with a focus on the three dimensional behaviours at specific areas. The model is based on the time dependent Reynolds-Averaged Navier-Stokes equations with a non-hydrostatic pressure distribution and a baroclinic force being incorporated in the three dimensional (3D) model. The two sub models interact dynamically during the solution procedure with no time-step restriction due to integration. The main idea is to use a fractional step algorithm for each model and then integrate the two models fraction by fraction. Hybrid 2D-3D finite volume cells have been introduced for the link nodes which are partly in the 2D domain and partly in the 3D domain. Thus an interpolation/averaging procedure at the interface and domain overlapping is no longer needed. The 3D model uses the projection method for pressure calculation. The advection equation is solved by the semi-Lagrangian method. Other components are solved via the finite element - finite volume (FV) method. The water surface is determined implicitly through a global matrix equation created by assembling the domain's matrices. The cell integrals are calculated analytically to eliminate a common source of numerical diffusion due to the use of approximation techniques for the FV integrals. The horizontal gradients of the density and shear stresses are calculated on true horizontal planes, in order to avoid artificial velocity and diffusion in highly stratified flows. Neumann interpolation elements with virtual nodes have been introduced at Neumann type of boundaries for more accuracy. The integrated model has been verified using analytical solutions and benchmark test cases, including the Ekman velocity distribution, wind driven circulation, lock exchange and integrated 2D-3D flows in basin. The results show the model is capable of the model for accurate simulation and implicit 2D-3D integration. Keywords: integrated modelling, hydrodynamic numerical model, non-hydrostatic, unstructured mesh, hybrid finite element finite volume method.EThOS - Electronic Theses Online ServiceGBUnited Kingdo

Integrated 2D-3D free surface hydro-environmental modelling

An integrated horizontally two- and fully three-dimensional numerical model system has been developed based on a combined unstructured and σ-coordinate grid to simulate the flow and water quality process in large water bodies with a focus on the three dimensional behaviours at specific areas. The model is based on the time dependent Reynolds-Averaged Navier-Stokes equations with a non-hydrostatic pressure distribution and a baroclinic force being incorporated in the three dimensional (3D) model. The two sub models interact dynamically during the solution procedure with no time-step restriction due to integration. The main idea is to use a fractional step algorithm for each model and then integrate the two models fraction by fraction. Hybrid 2D-3D finite volume cells have been introduced for the link nodes which are partly in the 2D domain and partly in the 3D domain. Thus an interpolation/averaging procedure at the interface and domain overlapping is no longer needed. The 3D model uses the projection method for pressure calculation. The advection equation is solved by the semi-Lagrangian method. Other components are solved via the finite element - finite volume (FV) method. The water surface is determined implicitly through a global matrix equation created by assembling the domain's matrices. The cell integrals are calculated analytically to eliminate a common source of numerical diffusion due to the use of approximation techniques for the FV integrals. The horizontal gradients of the density and shear stresses are calculated on true horizontal planes, in order to avoid artificial velocity and diffusion in highly stratified flows. Neumann interpolation elements with virtual nodes have been introduced at Neumann type of boundaries for more accuracy. The integrated model has been verified using analytical solutions and benchmark test cases, including the Ekman velocity distribution, wind driven circulation, lock exchange and integrated 2D-3D flows in basin. The results show the model is capable of the model for accurate simulation and implicit 2D-3D integration. Keywords: integrated modelling, hydrodynamic numerical model, non-hydrostatic, unstructured mesh, hybrid finite element finite volume method