Skip to main content

Lund University Publications

LUND UNIVERSITY LIBRARIES

Numerical schemes for a moving-boundary convection-diffusion-reaction model of sequencing batch reactors

Bürger, Raimund LU ; Careaga, Julio LU ; Diehl, Stefan LU and Pineda, Romel LU (2023) In ESAIM: Mathematical Modelling and Numerical Analysis 57(5). p.2931-2976
Abstract

Sequencing batch reactors (SBRs) are devices widely used in wastewater treatment, chemical engineering, and other areas. They allow for the sedimentation and compression of solid particles of biomass simultaneously with biochemical reactions with nutrients dissolved in the liquid. The kinetics of these reactions may be given by one of the established activated sludge models (ASMx). An SBR is operated in various stages and is equipped with a movable extraction and fill device and a discharge opening. A one-dimensional model of this unit can be formulated as a moving-boundary problem for a degenerating system of convection-diffusion-reaction equations whose unknowns are the concentrations of the components forming the solid and liquid... (More)

Sequencing batch reactors (SBRs) are devices widely used in wastewater treatment, chemical engineering, and other areas. They allow for the sedimentation and compression of solid particles of biomass simultaneously with biochemical reactions with nutrients dissolved in the liquid. The kinetics of these reactions may be given by one of the established activated sludge models (ASMx). An SBR is operated in various stages and is equipped with a movable extraction and fill device and a discharge opening. A one-dimensional model of this unit can be formulated as a moving-boundary problem for a degenerating system of convection-diffusion-reaction equations whose unknowns are the concentrations of the components forming the solid and liquid phases, respectively. This model is transformed to a fixed computational domain and is discretized by an explicit monotone scheme along with an alternative semi-implicit variant. The semi-implicit variant is based on solving, during each time step, a system of nonlinear equations for the total solids concentration followed by the solution of linear systems for the solid component percentages and liquid component concentrations. It is demonstrated that the semi-implicit scheme is well posed and that both variants produce approximations that satisfy an invariant region principle: solids concentrations are nonnegative and less or equal to a set maximal one, percentages are nonnegative and sum up to one, and substrate concentrations are nonnegative. These properties are achieved under a Courant-Friedrichs-Lewy (CFL) condition that is less restrictive for the semi-implicit than for the explicit variant. Numerical examples with realistic parameters illustrate that as a consequence, the semi-implicit variant is more efficient than the explicit one.

(Less)
Please use this url to cite or link to this publication:
author
; ; and
organization
publishing date
type
Contribution to journal
publication status
published
subject
keywords
Convection-diffusion-reaction model, Invariant region principle, Monotone scheme, Moving boundary, Semi-implicit method, Sequencing batch reactor
in
ESAIM: Mathematical Modelling and Numerical Analysis
volume
57
issue
5
pages
46 pages
publisher
EDP Sciences
external identifiers
  • scopus:85173250413
ISSN
2822-7840
DOI
10.1051/m2an/2023068
language
English
LU publication?
yes
id
2f354a22-11fa-44d8-8985-e9b6088bb233
date added to LUP
2023-12-12 08:30:24
date last changed
2024-02-09 09:49:10
@article{2f354a22-11fa-44d8-8985-e9b6088bb233,
  abstract     = {{<p>Sequencing batch reactors (SBRs) are devices widely used in wastewater treatment, chemical engineering, and other areas. They allow for the sedimentation and compression of solid particles of biomass simultaneously with biochemical reactions with nutrients dissolved in the liquid. The kinetics of these reactions may be given by one of the established activated sludge models (ASMx). An SBR is operated in various stages and is equipped with a movable extraction and fill device and a discharge opening. A one-dimensional model of this unit can be formulated as a moving-boundary problem for a degenerating system of convection-diffusion-reaction equations whose unknowns are the concentrations of the components forming the solid and liquid phases, respectively. This model is transformed to a fixed computational domain and is discretized by an explicit monotone scheme along with an alternative semi-implicit variant. The semi-implicit variant is based on solving, during each time step, a system of nonlinear equations for the total solids concentration followed by the solution of linear systems for the solid component percentages and liquid component concentrations. It is demonstrated that the semi-implicit scheme is well posed and that both variants produce approximations that satisfy an invariant region principle: solids concentrations are nonnegative and less or equal to a set maximal one, percentages are nonnegative and sum up to one, and substrate concentrations are nonnegative. These properties are achieved under a Courant-Friedrichs-Lewy (CFL) condition that is less restrictive for the semi-implicit than for the explicit variant. Numerical examples with realistic parameters illustrate that as a consequence, the semi-implicit variant is more efficient than the explicit one.</p>}},
  author       = {{Bürger, Raimund and Careaga, Julio and Diehl, Stefan and Pineda, Romel}},
  issn         = {{2822-7840}},
  keywords     = {{Convection-diffusion-reaction model; Invariant region principle; Monotone scheme; Moving boundary; Semi-implicit method; Sequencing batch reactor}},
  language     = {{eng}},
  month        = {{09}},
  number       = {{5}},
  pages        = {{2931--2976}},
  publisher    = {{EDP Sciences}},
  series       = {{ESAIM: Mathematical Modelling and Numerical Analysis}},
  title        = {{Numerical schemes for a moving-boundary convection-diffusion-reaction model of sequencing batch reactors}},
  url          = {{http://dx.doi.org/10.1051/m2an/2023068}},
  doi          = {{10.1051/m2an/2023068}},
  volume       = {{57}},
  year         = {{2023}},
}