umu.sePublications
Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
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
2015.
Keyword [en]
Mass-consistent model, urban environment, numerical method, algebraic multigrid
National Category
Computational Mathematics Physical Sciences
Identifiers
URN: urn:nbn:se:umu:diva-105499OAI: oai:DiVA.org:umu-105499DiVA: diva2:825911
External cooperation
Totalförsvarets forskningsinstitut, FOI
Subject / course
Examensarbete i teknisk fysik
Educational program
Master of Science Programme in Engineering Physics
Supervisors
Examiners
Available from: 2015-06-24 Created: 2015-06-24 Last updated: 2015-06-24Bibliographically approved

Open Access in DiVA

fulltext(1798 kB)254 downloads
File information
File name FULLTEXT01.pdfFile size 1798 kBChecksum SHA-512
fc9318efd8c2446ae6c5cf5138cde0bff7c3bfba4a3fd0526fb704d0ca25c02fa41593d03187bbdaa57216f193faec84bd561eb08f30d3dc1e29f0a2d81bd426
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: 254 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

urn-nbn

Altmetric score

urn-nbn
Total: 350 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf