On numerical broadening of particle-size spectra: a condensational growth study using PyMPDATA 1.0
<p>This work discusses the numerical aspects of representing the condensational growth of particles in models of aerosol systems such as atmospheric clouds. It focuses on the Eulerian modelling approach, in which fixed-bin discretisation is used for the probability density function describing...
| Main Authors: | , , , , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Copernicus Publications
2022-05-01
|
| Series: | Geoscientific Model Development |
| Online Access: | https://gmd.copernicus.org/articles/15/3879/2022/gmd-15-3879-2022.pdf |
| _version_ | 1817990319813165056 |
|---|---|
| author | M. A. Olesik J. Banaśkiewicz P. Bartman M. Baumgartner M. Baumgartner S. Unterstrasser S. Arabas S. Arabas |
| author_facet | M. A. Olesik J. Banaśkiewicz P. Bartman M. Baumgartner M. Baumgartner S. Unterstrasser S. Arabas S. Arabas |
| author_sort | M. A. Olesik |
| collection | DOAJ |
| description | <p>This work discusses the numerical aspects of representing the condensational growth of particles in models of aerosol systems such as atmospheric clouds.
It focuses on the Eulerian modelling approach, in which fixed-bin discretisation is used for the
probability density function describing the particle-size spectrum.
Numerical diffusion is inherent to the employment of the fixed-bin discretisation
for solving the arising transport problem (advection equation describing size spectrum evolution).
The focus of this work is on a technique for reducing the numerical diffusion in solutions
based on the upwind scheme: the multidimensional positive definite advection transport algorithm (MPDATA).
Several MPDATA variants are explored including infinite-gauge, non-oscillatory, third-order terms and recursive antidiffusive correction (double-pass donor cell, DPDC) options.
Methodologies for handling coordinate transformations associated with both particle-size spectrum coordinate choice and with numerical grid layout choice are expounded.
Analysis of the performance of the scheme for different discretisation parameters and different settings of the algorithm is performed using (i) an analytically solvable box-model test case and (ii) the single-column kinematic driver (“KiD”) test case in which the size-spectral advection due to condensation is solved simultaneously with the advection in the vertical spatial coordinate, and in which the supersaturation evolution is coupled with the droplet growth through water mass budget.
The box-model problem covers size-spectral dynamics only; no spatial dimension is considered.
The single-column test case involves a numerical solution of a two-dimensional advection problem
(spectral and spatial dimensions).
The discussion presented in the paper covers size-spectral, spatial and temporal convergence as well as computational cost, conservativeness and quantification of the numerical broadening of the particle-size spectrum.
The box-model simulations demonstrate that, compared with upwind solutions, even a 10-fold decrease in the spurious numerical spectral broadening can be obtained by an apt choice of the MPDATA variant (maintaining the same spatial and temporal resolution), yet at an increased computational cost.
Analyses using the single-column test case reveal that the width of the droplet size spectrum is affected by numerical diffusion pertinent to both spatial and spectral advection.
Application of even a single corrective iteration of MPDATA robustly decreases the relative dispersion
of the droplet spectrum, roughly by a factor of 2 at the levels of maximal liquid water content.
Presented simulations are carried out using PyMPDATA – a new open-source Python implementation of MPDATA based on the
Numba just-in-time compilation infrastructure.</p> |
| first_indexed | 2024-04-14T00:57:53Z |
| format | Article |
| id | doaj.art-ee22e9882ea4438aa673386d1b3f16bb |
| institution | Directory Open Access Journal |
| issn | 1991-959X 1991-9603 |
| language | English |
| last_indexed | 2024-04-14T00:57:53Z |
| publishDate | 2022-05-01 |
| publisher | Copernicus Publications |
| record_format | Article |
| series | Geoscientific Model Development |
| spelling | doaj.art-ee22e9882ea4438aa673386d1b3f16bb2022-12-22T02:21:32ZengCopernicus PublicationsGeoscientific Model Development1991-959X1991-96032022-05-01153879389910.5194/gmd-15-3879-2022On numerical broadening of particle-size spectra: a condensational growth study using PyMPDATA 1.0M. A. Olesik0J. Banaśkiewicz1P. Bartman2M. Baumgartner3M. Baumgartner4S. Unterstrasser5S. Arabas6S. Arabas7Faculty of Physics, Astronomy and Applied Computer Science, Jagiellonian University, Kraków, PolandFaculty of Mathematics and Computer Science, Jagiellonian University, Kraków, PolandFaculty of Mathematics and Computer Science, Jagiellonian University, Kraków, PolandZentrum für Datenverarbeitung, Johannes Gutenberg University Mainz, Mainz, GermanyInstitute for Atmospheric Physics, Johannes Gutenberg University Mainz, Mainz, GermanyInstitute of Atmospheric Physics, German Aerospace Center (DLR), Oberpfaffenhofen, GermanyFaculty of Mathematics and Computer Science, Jagiellonian University, Kraków, PolandDepartment of Atmospheric Sciences, University of Illinois at Urbana-Champaign, Urbana, IL, USA<p>This work discusses the numerical aspects of representing the condensational growth of particles in models of aerosol systems such as atmospheric clouds. It focuses on the Eulerian modelling approach, in which fixed-bin discretisation is used for the probability density function describing the particle-size spectrum. Numerical diffusion is inherent to the employment of the fixed-bin discretisation for solving the arising transport problem (advection equation describing size spectrum evolution). The focus of this work is on a technique for reducing the numerical diffusion in solutions based on the upwind scheme: the multidimensional positive definite advection transport algorithm (MPDATA). Several MPDATA variants are explored including infinite-gauge, non-oscillatory, third-order terms and recursive antidiffusive correction (double-pass donor cell, DPDC) options. Methodologies for handling coordinate transformations associated with both particle-size spectrum coordinate choice and with numerical grid layout choice are expounded. Analysis of the performance of the scheme for different discretisation parameters and different settings of the algorithm is performed using (i) an analytically solvable box-model test case and (ii) the single-column kinematic driver (“KiD”) test case in which the size-spectral advection due to condensation is solved simultaneously with the advection in the vertical spatial coordinate, and in which the supersaturation evolution is coupled with the droplet growth through water mass budget. The box-model problem covers size-spectral dynamics only; no spatial dimension is considered. The single-column test case involves a numerical solution of a two-dimensional advection problem (spectral and spatial dimensions). The discussion presented in the paper covers size-spectral, spatial and temporal convergence as well as computational cost, conservativeness and quantification of the numerical broadening of the particle-size spectrum. The box-model simulations demonstrate that, compared with upwind solutions, even a 10-fold decrease in the spurious numerical spectral broadening can be obtained by an apt choice of the MPDATA variant (maintaining the same spatial and temporal resolution), yet at an increased computational cost. Analyses using the single-column test case reveal that the width of the droplet size spectrum is affected by numerical diffusion pertinent to both spatial and spectral advection. Application of even a single corrective iteration of MPDATA robustly decreases the relative dispersion of the droplet spectrum, roughly by a factor of 2 at the levels of maximal liquid water content. Presented simulations are carried out using PyMPDATA – a new open-source Python implementation of MPDATA based on the Numba just-in-time compilation infrastructure.</p>https://gmd.copernicus.org/articles/15/3879/2022/gmd-15-3879-2022.pdf |
| spellingShingle | M. A. Olesik J. Banaśkiewicz P. Bartman M. Baumgartner M. Baumgartner S. Unterstrasser S. Arabas S. Arabas On numerical broadening of particle-size spectra: a condensational growth study using PyMPDATA 1.0 Geoscientific Model Development |
| title | On numerical broadening of particle-size spectra: a condensational growth study using PyMPDATA 1.0 |
| title_full | On numerical broadening of particle-size spectra: a condensational growth study using PyMPDATA 1.0 |
| title_fullStr | On numerical broadening of particle-size spectra: a condensational growth study using PyMPDATA 1.0 |
| title_full_unstemmed | On numerical broadening of particle-size spectra: a condensational growth study using PyMPDATA 1.0 |
| title_short | On numerical broadening of particle-size spectra: a condensational growth study using PyMPDATA 1.0 |
| title_sort | on numerical broadening of particle size spectra a condensational growth study using pympdata 1 0 |
| url | https://gmd.copernicus.org/articles/15/3879/2022/gmd-15-3879-2022.pdf |
| work_keys_str_mv | AT maolesik onnumericalbroadeningofparticlesizespectraacondensationalgrowthstudyusingpympdata10 AT jbanaskiewicz onnumericalbroadeningofparticlesizespectraacondensationalgrowthstudyusingpympdata10 AT pbartman onnumericalbroadeningofparticlesizespectraacondensationalgrowthstudyusingpympdata10 AT mbaumgartner onnumericalbroadeningofparticlesizespectraacondensationalgrowthstudyusingpympdata10 AT mbaumgartner onnumericalbroadeningofparticlesizespectraacondensationalgrowthstudyusingpympdata10 AT sunterstrasser onnumericalbroadeningofparticlesizespectraacondensationalgrowthstudyusingpympdata10 AT sarabas onnumericalbroadeningofparticlesizespectraacondensationalgrowthstudyusingpympdata10 AT sarabas onnumericalbroadeningofparticlesizespectraacondensationalgrowthstudyusingpympdata10 |