Optimal low order methods of moments for solving the fragmentation equation

SARS-CoV-2 adsorption on suspended solids along a sewerage network: mathematical model formulation, sensitivity analysis, and parametric study

Accounting for SARS-CoV-2 adsorption on solids suspended in wastewater is a necessary step towards the reliable estimation of virus shedding rate in a sewerage system, based on measurements performed at a terminal collection station, i.e., at the entrance of a wastewater treatment plant. This concept is extended herein to include several measurement stations across a city to enable the estimation of spatial distribution of virus shedding rate. This study presents a pioneer general model describing the most relevant physicochemical phenomena with a special effort to reduce the complicated algebra. This is performed both in the topology regime, introducing a discrete-continuous approach, and in the domain of independent variables, introducing a monodisperse moment method to reduce the dimensionality of the resulting population balance equations. The resulting simplified model consists of a large system of ordinary differential equations. A sensitivity analysis is performed with respect to some key parameters for a single pipe topology. Specific numerical techniques are employed for the integration of the model. Finally, a parametric case study for an indicative-yet realistic-sewerage piping system is performed to show how the model is applied to SARS-CoV-2 adsorption on wastewater solids in the presence of other competing species. This is the first model of this kind appearing in scientific literature and a first step towards setting up an inverse problem to assess the spatial distribution of virus shedding rate based on its concentration in wastewater.

On the approximate solutions of fragmentation equations

A numerical model based on the finite volume scheme is proposed to approximate the binary breakage problems. Initially, it is considered that the particle fragments are characterized by a single property, i.e. particle’s volume. We then investigate the extension of the proposed model for solving breakage problems considering two properties of particles. The efficiency to estimate the different moments with good accuracy and simple extension for multi-variable problems are the key features of the proposed method. Moreover, the mathematical convergence analysis is performed for one-dimensional problems. All mathematical findings and numerical results are validated over several test problems. For numerical validation, we propose the extension of Bourgade & Filbet (2008 Math. Comput. 77 , 851–882. ( doi:10.1090/S0025-5718-07-02054-6 )) model for solving two-dimensional pure breakage problems. In this aspect, numerical treatment of the two-dimensional binary breakage models using finite volume methods can be treated to be the first instance in the literature.

A study of the collisional fragmentation problem using the Gamma distribution approximation

The nonlinear fragmentation population balance formulation has been elevated in recent years from a prototype for studying nonlinear integro-differential equations to a vehicle for analyzing and understanding several physicochemical processes of technological interest. The so-called pure collisional fragmentation, which is the particular mode of nonlinear fragmentation induced by collisions between particles, is studied here. It is shown that the corresponding population balance equation admits large time asymptotic (self-similarity) solutions for homogeneous fragmentation and collision functions (kernels). The self-similar solutions are given in closed form for some simple kernels. Based on the shape of the self-similar solutions the method of moments with Gamma distribution approximation is employed for transient solution (from initial state to establishment of the asymptotic shape) of the collisional fragmentation equation. These solutions are presented for several sets of parameters and their behavior is discussed rather extensively. The present study is similar to the one has already been performed for the case of the much simpler linear fragmentation equation [G. Madras, B.J. McCoy, AIChE J. 44 (1998) 647].

Cluster aggregation and fragmentation kinetics model for gelation

Gelation can occur in polymer, hydrogel, and colloid systems that undergo reversible aggregation-fragmentation (crosslinking accompanied by breakage). Gelation, characterized by rapid divergence of weight-average molecular weight and viscosity due to initial network formation, can be reversed if conditions change. In this paper, reversible aggregation and fragmentation in the pre-gelation time period are modeled with distribution kinetics. Moment equations are obtained from the population balance equation, and solved for eight different rate kernels. We identify the cases for which gelation is possible and obtain the critical values for the rate constants that allow gelation. The model provides a good simulation of published experimental data for aggregation and degradation of plasticized wheat gluten during thermo-mechanical treatments. We also evaluate two closure approximations based on Gamma and log-normal distributions, and conclude that log-normal closure predicts all five possible steady states, in agreement with the Vigil-Ziff criterion, and Gamma closure predicts only three. However, Gamma closure approximates the steady state either closely or exactly, whereas log-normal closure only poorly approximates the steady-state distribution.