Change search
ReferencesLink to record
Permanent link

Direct link
Algebraic multigrid for a mass-consistent wind model, the Nordic Urban Dispersion model
Umeå University, Faculty of Science and Technology, Department of Physics.
2015 (English)Independent thesis Advanced level (degree of Master (Two Years)), 300 HE creditsStudent thesis
Abstract [en]

In preparation for, and for decision support during, CBRN (chemical, biological, radiological and nuclear) emergencies it is essential to know how such an event would turn out, so that one can prepare a possible evacuation. Afterwards it might be good to know how to backtrack and see what caused the emergency, and in the case of e.g. a gas leak, where did it begin? The Swedish Defence Research Agency (FOI) develops models for such scenarios.

In this thesis FOI's model, "The Nordic Urban Dispersion model" (NUD), has been studied. The system of equations set up by this model was originally solved using Intel's PARDISO solver, which is a direct solver. An evaluation on how an iterative multigrid method would work to solve the system has been done in this thesis. The wind model is a mass-consistent model which sets up a diagnostic initial wind field. The final wind field is later minimized under the constraint of the continuity equation. The minimization problem is solved using Lagrange multipliers and the system turns into a Poisson-like problem.

The iterative algebraic multigrid solver (AMG) which has been evaluated had difficulties solving the problem of an asymmetric system matrix generated by NUD. The AMG solver was then tried on a symmetric discrete Poisson problem instead, and the solution turns out to be the same as for the PARDISO solver. A comparison was made between the AMG and PARDISO solver, and for the discrete Poisson case the AMG solver turned out on top for both larger system size and less computational time.

To try out the solvers for the original NUD case a modification of the boundary conditions was made to make the system matrix symmetric. This modification turns the problem into a mathematical problem rather than a physical one, as the wind fields generated are not physically correct. For this modified case both the solvers get the same solution in essentially the same computational time. A method of how to in the future solve the original (asymmetric) problem, by modifying the discretization of the boundary conditions, has been discussed.

Place, publisher, year, edition, pages
Keyword [en]
Mass-consistent model, urban environment, numerical method, algebraic multigrid
National Category
Computational Mathematics Physical Sciences
URN: urn:nbn:se:umu:diva-105499OAI: diva2:825911
External cooperation
Totalförsvarets forskningsinstitut, FOI
Subject / course
Examensarbete i teknisk fysik
Educational program
Master of Science Programme in Engineering Physics
Available from: 2015-06-24 Created: 2015-06-24 Last updated: 2015-06-24Bibliographically approved

Open Access in DiVA

fulltext(1798 kB)106 downloads
File information
File name FULLTEXT01.pdfFile size 1798 kBChecksum SHA-512
Type fulltextMimetype application/pdf

Search in DiVA

By author/editor
Pogulis, Markus
By organisation
Department of Physics
Computational MathematicsPhysical Sciences

Search outside of DiVA

GoogleGoogle Scholar
Total: 106 downloads
The number of downloads is the sum of all downloads of full texts. It may include eg previous versions that are now no longer available

Total: 236 hits
ReferencesLink to record
Permanent link

Direct link