Chemical Model for Thermal Treatment of Sewage Sludge

: Sewage sludge is here studied as a valuable source for processing or energy conversation thanks to its high nutrition and energy content. However, various origins of the wastewater, different water cleaning technologies, and seasonal and regional dependencies lead to the high variability of the sewage sludge properties. In thermal treatment units, that is, incineration, gasiﬁcation and pyrolysis, sewage sludge serves as feedstock or fuel, hence a proper characterization and a mathematical description of the sewage sludge are required to estimate product streams and to formulate numerical simulations and optimization methods. The presented work introduces a surrogate concept that allows replication of sewage sludge’s ultimate composition, moisture, and ash content. The surrogate approach aims to model the decomposition of any sewage sludge sample, opposite to the established determination of kinetic rates for individual samples. Based on chemical solid surrogate species and corresponding reaction mechanisms, the thermal decomposition path is described. Sewage sludge is represented by a combination of lignocellulosic species, proteins, sugars, lipids, and representative inorganic species. The devolatilization and heterogeneous reactions are formulated such that they can be used together with a detailed gas-phase model, including tar oxidation and emission models for nitrogen and sulfur oxides, recently proposed by the authors. The developed chemical model is applied using a zero-dimensional gasiﬁcation reactor in order to model weight loss within the thermogravimetric analysis, pyrolysis, gasiﬁcation and combustion conditions. Weight loss, the composition of product gases, and emission release (nitrogen and sulfur oxides) are captured well by the model. The ﬂexible surrogate approach allows us to represent various sewage sludge samples.


Introduction
Sewage sludge is a by-or waste product accumulating in the wastewater cleaning process. Municipal wastewater combines domestic, industrial and surface water, and hence, the sewage sludge includes leftovers from households, industry, and agriculture. As with other waste streams, nowadays, sewage sludge is considered a resource, and the recovery of materials and energy content are targeted [1]. The sludge typically has high contents of nutrients such as nitrogen (1.5-5.6 wt% dry) and phosphorus (0.6-2.8 wt% dry) [2][3][4], which makes it an attractive fertilizer for agriculture as well as raw material for composting [5,6]. Nevertheless, due to the diverse origin of the wastewater, sewage sludge contains different harmful contaminants. These contaminants consist of inorganic compounds, for example, heavy metals [3,4,6], and organic substances, for example, pesticides, pharmaceuticals [7,8]. Because of this potential hazard to the environment and human health, legislative restrictions are implemented for direct agricultural use or depositing to landfillings [4,9]. At the same time, with a growing world population, the accumulated volume of sewage sludge is rising [6,9]. Thermal treatment of sewage sludge is a measure for waste management and the integration of the resource into the future circular economy [6]. Different thermal treatment units allow for the re-use of the waste, recovering materials or energy and, thanks to the high temperatures, sterilization or destruction of toxic organic compounds [4,10]. The processes are incineration, gasification, pyrolysis, or a combination of these technologies. During incineration, the sludge's energy content is used for heat generation, and the volume of the solid waste stream is reduced by 90% [8]. In gasification devices, the sewage sludge is partially oxidized to yield a high caloric producer gas or syngas used for heat and power generation [8,11] and a char fraction. Pyrolysis is the endothermic conversion under a reducing atmosphere and favors the formation of tar (oils) and char. Thermodynamic conditions and residence time decide on the quality and amount of the products (oil, gas, and char) [8]. Depending on the origin of the wastewater and its treatment process, the properties of sewage sludge vary widely. For all thermal treatment processes, the varying feedstock is one of the key challenges. Here, the most relevant properties are the moisture content, due to the high energy demand for drying [6,12], the ratio of volatile components and ash-since this is correlated to the lower heating value of the sewage sludge [13]-and the content of elements involved in emission formation processes, that is, nitrogen, sulfur and chlorine [14][15][16][17][18]. To handle this challenging fuel, often sewage sludge is mixed with other solid fuels, such as biomass wastes, in co-gasification [19] and co-incineration in the cement industry [4,20].
In order to support the further technological development of thermal treatment units, it is desirable to employ simulations to understand further and optimize the processes. Goals are improved thermal efficiency, amount and quality of the products, and minimal emission formation. Since the variability of the sewage sludge properties is a determining factor within the conversion process, the model formulation must be sensitive to changes in the feedstock. The classical approach in biomass and waste research is to determine kinetic parameters for specific waste fractions and their mixtures [15,[21][22][23][24][25][26][27][28][29] or to determine pseudo representative species for an individual feedstock [30]. These methods are based on empirical material-specific apparent reaction rates and cannot directly be transferred to other feedstocks [31]. Contrary, the flexible representation of various solid fuels using reoccurring characteristic surrogate species offers the description of the needed dependencies of the solid fuel and the gas phase products (gas, tar, and emission precursors) [31][32][33]. The surrogate concept is the state-of-the-art in liquid (e.g., Diesel and gasoline) and gaseous fuels modeling and has previously been used within the solid fuel conversion of coal and biomass [32][33][34][35][36][37][38]. For example, Debiagi et al. [36] achieved the mathematical description of more than 500 biomass and waste fraction samples, such as woods, grass plants, algae and food industry wastes, with nine reference species, including moisture and ash. Recently, the concept was also applied to describe the heterogeneous feedstock of municipal solid waste (MSW) [31]. Žnidarčič et al. [39] empathized the ability of a surrogate model to account for different sewage sludge properties and introduced a surrogate model for the description of the gas phase. However, to the author's best knowledge, the literature lacks a surrogate description and corresponding model ranging from the solid feedstock via its thermal conversion to the resulting gas-phase products and emission precursors.
The presented work aims to formulate a mathematical description of sewage sludge to model thermal treatment units. The model is formulated by selecting existing surrogate components from the literature [31,32,38] and their devolatilization schemes and combing them with corresponding gas-phase chemistry [31]. The model development is described in Section 2.2, discussing the surrogate species selection (Section 2.2.1), their devolatilization scheme (Section 2.2.2), and their combination to a sewage sludge surrogate (Section 2.2.3). The approach is validated towards ultimate composition measurements (Section 3.1) and thermogravimetric (TG) measurements (Section 3.2), for pyrolysis and gasification in a technical scale application (Section 3.3) and the emission release under an incineration atmosphere (Section 3.3).

Numerical Model
The stochastic gasification model available in LOGEresearch version 1.10 [40] is chosen for the chemical model development and validation against different experiments from the literature. It is selected since it is targeted to solve large detailed chemistry schemes and to resolve the devolatilization, heterogeneous reactions of solid and gas phase, and reactions in the gas phase. Depending on the provided chemistry and operating conditions, it allows us to model pyrolysis, gasification, and combustion applications [31,41,42].
This gasification model is based on three physical phases-the solid particles, the pore gas within the particles, and the bulk gas. The reactor's volume is discretized into nondimensional virtual packages. These virtual packages are, in the following, called stochastic particles and include each of the three phases. The conversion process of the solid fuel is governed by stochastic heating, mixing of the gas phase, and chemistry integration. Preheated gas flows and hot reactor walls can heat up the reactor. In the mixing step, heat is also transferred between the wall and solid particles, between solid particles, and between stochastic particles' gas phases. The stochastic mixing mimics the turbulent gas flow and is controlled via the mixing time. The mixing time is correlated to the turbulent kinetic energy and its dissipation ratio and describes the mixing events within a numerical time step. Infinitely fast mixing would result in a homogeneous mixture, while low mixing conserves and evolves inhomogeneity within the reactor [42]. Highly turbulent flows, as in entrained flow gasifiers, are described by mixing times between 0.001 to 0.05 s [40], while for grate-fired applications, 1 s is assumed to be representative [42], which results in one mixing event per time step. The stochastic heat transfer and the mixing lead to a temperature distribution over the stochastic particles. This temperature field leads to locally different (stochastic particle level) predictions of species and heat release. This changes in the discrete realization of the joint probability density function of species and enthalpy that is further evolved in each time step. While the ensemble of the stochastic particles reproduces a distribution in species concentrations (gas and solid phase) and temperatures within the reactor, their concentration within a stochastic particle is assumed to be homogeneous in order to close the chemistry source terms. In the chemistry step, the equations for devolatilization and the surface reactions are first solved together with that of the mass transfer. Second, gas-phase kinetics and radiation are addressed. The mass source term from solid to gas phase is used to update the representative diameters of the solid particles. A more detailed description of the model is provided elsewhere [40][41][42].
The reactor can be run as a single or series of partially stirred reactors [40,41], partially stirred plug flow reactors [40,41], or connected to a reactor network model for grate-fired applications [31,42].

Selection of the Surrogate Species
Sewage sludge typically has a content of 40 wt% to 60 wt% of volatile solid matter (mass percent, dry, ash-included) considered a valuable bio-originated energy source [4]. This volatile matter has due to their domestic, industrial and agriculture origin a complex composition including carbohydrates, protein, oil, lignocellulosic compounds, humic substances and more [3]. The sludges composition depends on the origin and the applied wastewater treatment technology, and varies largely. The ultimate composition (mass percent, dry, ash-included) can vary; for carbon from 20 wt% to 42 wt%, for hydrogen from 2.5 wt% to 6.8 wt%, oxygen from 8 wt% to 38 wt% and nitrogen from 2 wt% to 7 wt%. These ranges are rounded minimum and maximum values reported in the literature [10,43,44] and accessed via the Phyllis2 database [45][46][47][48][49][50][51]. Their variation is summarized in Figure 1. Proteins make up a significant amount of the volatile matter (20 wt% to 40 wt%), followed by lignin with 23 wt% to 30 wt%, and other lignocellulosic compounds and their decomposition product such as sugars [3,4]. For example, Adar et al. [3] have analyzed three sewage sludge samples regarding their organic component composition and reported their average values with 30.09 wt% lignin, 27.40 wt% crude proteins, 11.11 wt% cellulose (9.60 wt% crude cellulose), 4.12 wt% hemicellulose, and 2.45 wt% sugar. Other investigations show that a low oxygen content indicates a high content of oil and grease [52]. In the present numerical scheme, the lignocellulosic compounds are described using the biomass devolatilization mechanism and its respective surrogate species by Ranzi et al. [32]. It includes cellulose (CELL), hemicellulose (HCE), and three types of lignin (LIG C , LIG H , and LIG O ). The three proteins (PROT C , PROT H and PROT O ), introduced by Debiagi et al. [37] represent the group of proteins. A species representing sugars (SUGAR), and a species representing oils and grease (LIPID), is adopted from the same reaction scheme for the thermal decomposition of algae biomass.
The inorganic components present in the solid matter of sewage sludge are also modeled, describing their product gases of reactions of the inorganic components and the gas phase. The model includes the release of ammonia (NH 3 ), carbon dioxide (CO 2 ) and sulfur species, wich are hydrogen sulfide (H 2 S), sulfur dioxide (SO 2 ) and carbonyl sulfide (COS). The reactions are inherited from Debiagi et al. [37] and Netzer et al. [31]. The ash content of 18 wt% to 55 wt% is assumed to be inert and represented by ash (ASH), consisting of only silicon Si. The moisture content is modeled by the surrogate species H 2 O(S) that describes the drying process via an Arrhenius type equation [32]. All surrogate species and their molecular composition are listed in Table 1 and the model development and origin of the surrogate species is summarized in Figure 2.

Reaction Scheme of the Surrogate Species
The devolatilization schemes for the volatile species are adopted from [32,37] and are illustrated in Figure 3. The typical pattern involves a first decomposition step forming a reduced or activated solid species. In this first decomposition step, depending on the surrogate species, char, metaplastic species (mp, gaseous species that are trapped in the solid phase) and products released to the gas phase (tar and gas) are formed. In successive steps, the intermediate solid species are further decomposed. Heterogeneous reactions of the char with the surrounding gas phase and Ahrrenius type formulations for the release of the metaplastic species complete the description. The inorganic species are freed in a one-step reaction, except for sulfur, modeled by a two-step release. In the following, only reactions that are modified from their original scheme to connect devolatilization and gas-phase chemistry are discussed. All modified reactions to connect the devolatilization model [37] and the gas-phase scheme [31] are listed in Table 2. The complete set of equations is provided in the Supplementary Materials.  Introduced modifications are based on further specification by Debiagi et al. [37], the discussions and modifications by Netzer et al. [31] and group distribution methods, and are discussed in the following. The nitrogen-containing tar (tar N ), released by the protein species, is assumed to consist of pyrrole (C 4 H 5 N), pyridine (C 5 H 5 N) and malondialdehyde (C 3 H 4 O 2 ) [31]. These species have been introduced to represent the cyclic molecule structure and the ketone group in the original proposed equimolar mixture of pyrrole, pyridine, and diketopiperazine [37]. The introduced species are modeled using available gas-phase mechanisms [32,53,54]. In the present mechanism, tar N is given by [31]: Oxygen-containing tar species (tar O,i ) in the formulation by Debiagi et al. [37] are written out based on the description in their original mechanism. The anhydrous monomer of polysaccharides (C 6 H 8 O 6 ) is replaced by xylosan included in the hemicellulose chemistry. Xylosan is the dehydrated product of the xylose sugar monomer and is considered representative here since the main purpose in this step is to keep the elemental balance and comparable speciation.
As specified in Debiagi et al. [37], tar O,2 is described by an equimolar mixture of glyoxal, acetol and furfural, and tar O,3 by a mixture of methanol, ethanol and acetic acid. These are all included in the gas phase chemistry and hence the species can be formulated as: No representative of free-fatty acids (FFA) is present in the selected gas-phase chemistry. FFA is decomposed using the group contribution method, using the species present in the gas phase mechanism not to further extend the size of the gas phase mechanism and hence the computational cost. Following group contribution theories, FFA is decomposed representing the long hydrocarbon tail by C 2 H 2 , C 2 H 4 , C 2 H 6 and acetic acid (CH 3 COOH) to represent the carboxyl group (COOH): Figure 4 shows the elemental ratios of the chosen surrogate species and sewage sludge samples [10,[43][44][45][46][47][48][49][50][51]. Depicted are the H/C ratios, O/C ratios, and N/C ratios. The triangles show the linear combinations of the lignocellulosic compounds and the proteins, respectively. One can see that the great majority of the sewage sludge samples lie within the area enclosed by the surrogate species in the O/C and H/C range. Hence, a linear combination of the surrogate species can represent the sewage sludge on an elemental basis. The species SUGAR and LIPID extend the covered range and are needed to represent outliers. Protein contains nitrogen while non of the lignocellulosic component include nitrogen. Using only protein species would overestimate the N/C ratio in sewage sludge, while using only lignocellulosic species underestimates it. However, linear combinations of both groups allow representation of the solid bonded nitrogen in the right amount.

Formulation of the Sewage Sludge Surrogates
The surrogate representation aims to replicate the ultimate composition of the analyzed sewage sludge, including its ash and moisture content. This is beneficial since the elemental composition is linked to the heating value of the fuel, and the amount of nitrogen and sulfur in the solid matter is linked to the emission of nitrogen o oxides (NO X ) and sulfur oxides (SO X ). A linear least-squares fit is employed to replicate the elemental composition [31]. The main targets for the least-squares fit are the dry mass percentages of carbon (C), hydrogen (H), nitrogen (N), and sulfur (S). Requirements for the species mass fractions (Y i ), 0 ≤ Y i ≤ 1 and ∑ Y i = 1, and fixed ratios for the lignocellulosic compounds [36,55] are additional constraints. Here, for sewage sludge, the ratios Bio 1 = 0.6CELL + 0.4HCE, Bio 2 = 0.95LIG O + 0.05LIG C and Bio 3 = 0.95LIG H + 0.05LIG C are applied. After the dry ash-free surrogate is created, the ash and moisture content of the sludge are linearly added.

Representation of Sewage Sludge
The comparison of sewage sludge samples and their surrogate representation is given in Table 3. Listed are sewage sludge samples that are further discussed in the validation part. ID 1 to 16 denotes samples reported by Gómez-Rico et al. [43] that have been analyzed using thermogravimetric analysis. Sample A is converted in a bench-scale rotary kiln by Freda et al. [10] and sample B is analyzed regarding the NO X , and SO X formation by Shao et al. [44]. All surrogates represent the sewage sludge in close agreement. The maximum relative error is found for the sulfur content of sample 9 with 11.71% and hydrogen in sample B with 22.6%. Both values are outliers and are given their absolute deviation of 0.18 wt% and 1.5 wt%, respectively, which is considered acceptable here. The absolute errors of carbon and hydrogen representation are significantly lower, and hence the energy content of the sewage sludge is captured well. The individual composition of the dry surrogates, including ash, is shown in Figure 5 for sample ID 1 to 16. As in detailed speciation of sewage sludge, proteins and linings are the major groups of organic matter. Covering a range from 13 wt% to 48 wt%, proteins are included in similar ranges as reported in literature (18 wt% to 40 wt% [3]). Lignin species vary largely between 3 wt% and 42 wt%. The average with most appearances of around 25 wt% lies well in between the previously mentioned 23 wt% to 29 wt%. Sugars and lipids are needed to fulfill the ultimate composition requirements. With a small amount up to 9 wt%, this appears plausible considering the measurement by Adar et al. [3]. Species representing inorganic materials contribute to the smallest weight fraction but are essential to balance nitrogen and sulfur.  [43]. Sorted, decreasing, by the carbon content in the dry ash containing sewage sludge.

Thermogravimetric Analysis
Gómez-Rico et al. [43] studied the pyrolysis and combustion of various sewage sludge samples (ID 1 to 16 in Table 3). For this, the authors performed thermogravimetric measurements under a helium and a helium/oxygen atmosphere with 80 vol% and 20 vol%, respectively. The analysis was performed using 5 mg sewage sludge samples, a gas flow rate of 0.06 L/min, and a heating rate of 10 K/min. The gasification reactor was used in a plug flow configuration with a low mixing rate of 1 s to represent these measurements [31]. Each surrogate is initialized separately, assigning a particle size of 150 µm and the two atmospheres. The results are shown in Figure 6.   Table 3 versus thermogravimetric measurements by Gómez-Rico et al. [43]. Shown are pyrolysis conditions, denoted with P, using helium and combustion, denoted with C, using a mixture of 20 vol% oxygen and 80 vol% helium. Overall, the model captures the mass-loss rate and the sewage sludge final solid residue. The early devolatilization is delayed for some samples, for example, ID 1, 2 and 4. In the model, this is caused by the lignocellulosic compounds [31]. However, for the majority of the sewage sludge samples, this stage is soundly represented by the model. The characteristic two-stage devolatilization [43] under pyrolysis conditions is represented in the simulation results. The response to different atmospheres (pyrolysis and combustion) is for the majority of the sewage sludge samples reproduced. Examples of a good prediction of both pyrolysis and combustion conditions are samples 1, 3, and 9. For these samples, the relative error of the weight loss prediction at 1050 K is <5% for the pyrolysis conditions and <7% for the combustion atmosphere. The model under-predicts the residue slightly for the sewage sludge samples, with the highest difference in the final solid residue (samples 10, 12, and 13). This under-prediction appears at temperatures >850 K. With an error of 11% and 21% of the weight loss prediction at 1050 K, the most significant deviation is found for samples 10 and 12, respectively, while the initial under-prediction for sample 13 is not apparent in the final weight loss deviation of 1%. The discussed under-prediction is not present for the combustion condition with an error <6%. The devolatilization is close to complete in this temperature range, and the char model dominates the solid conversion. The employed char conversion model consists of only tree reactions, and its revision is targeted in future work. To summarize, with an average deviation of the predicted weight loss at 1050 K of 6.5% and 4.5% for pyrolysis and incineration conditions, the model's aim to predict the rate of solid mass loss to predict gas and tar release are met for various sewage sludge samples.

Producer Gas Composition
For the validation of the obtained gas composition and share of tar and water, the experimental results reported by Freda et al. [10] are chosen. The authors analyze sewage sludge gasification in a bench-scale rotary kiln. The kiln has a total length of 1550 mm, of which 610 mm are electrically heated, with an inner diameter of 80 mm. Freda et al. [10] studied four operating conditions, given in Table 4, with air as a gasification agent. Condition (Cond) 1 corresponds to pyrolysis conditions that favor a high amount of tar species, and the remaining conditions account for gasification conditions. The distribution of the operating conditions allows validating the model response towards  The gasification reactor is run in a plug flow configuration for the simulation of this experimental setup, accounting for constant fuel and air supply at the reactor inlet. The mixing time is set to 0.033 s to represent the kiln rotational speed of 2 rpm. The surrogate A from Table 3 corresponds to the ultimate analysis reported by Freda et al. [10]. The authors further note that the sewage sludge had to be partially dried to prevent the rotary kiln from clogging. In the presented simulation results, a moisture content of 20 wt% is assumed.
All major product gases are predicted in close agreement with the experiments. The amounts of methane (CH 4 ) and the tar yield are however under and over predicted, respectively (Figures 7 and 8). For the tar yield, all species in the gas phase reaction mechanisms with a molar mass equal or larger than benzene (C 6 H 6 , W = 78.11 g/mol) are combined. The resulting under and over prediction shows that the simulation's tar conversion is slower than in the experiment. The thermal and chemical break down of the tar species, e.g., C 11 H 12 O 4 , C 6 H 10 O 5 and C 5 H 5 NO, would result in an increase of small hydrocarbons and consumption of hydrogen that is slightly over predicted. The model yet replicates the trends over changing gasification temperature and air supply, including CH 4 and tar. Further, besides the major gas species, ethane with only a small concentration is in good agreement with the experiment (Figure 7). The gas composition directly correlates to the lower heating value that is important for the further processing of the gas for energy supply. Given the species' good prediction, the LHV aligns well with the measurements (not shown). The combination of the product gases LHV, the gas yield and the LHV of the solid fuel, yield the cold gas efficiency (CGE = (LHV gas × Yield gas )/LHV f uel ). Also for this measure the model aligns well with the measurement.
The species and tar prediction of the model respond well to increasing airflow rate and temperatures (Cond 1 vs. Cond 2 vs. Cond 4) and the related properties, such as LHV and CGE. The model is further sensitive to individual variations in temperature (Cond 2 vs. Cond 3) and in the airflow (Cond 3 vs. Cond 4). The experimentally observed slight variations of water and the dry gas yield over the gasification conditions (Cond 2, Cond 3, and Cond 4) are not replicated by the model. The reason for this is twofold. First, in the simulations, the fuel feed rate and moisture content are the same, while this is not guaranteed in the experiments [10]. Second, in the presented model, the final char yield is reached around 1000 K ( Figure 6) and the sensitivity for the gas-phase species prediction is obtained by the gas-phase model. Based on local oxidizer and elevated temperatures tar and larger gas species are further decomposed to form the main producer gas species. Overall, the amount of gas, tar, and water yield and the product gas composition are well replicated by the surrogate model.

Emission Release Prediction
Shao et al. [44] performed incineration experiments of a sewage sludge sample from Shanghai. The authors analyzed the NO X and sulfur dioxide (SO 2 ) formation in a tubular furnace. Dried sewage sludge with a mass of 0.2 g was burned using a synthetic air mixture (80 vol% N 2 /20 vol% O 2 ) with a flow rate of 2 L/min. The ultimate analysis and the corresponding surrogate are denoted with B in Table 3. Shao et al. [44] conducted the experiments with a residence time of 600 s. To comply with the experiments, the gasification reactor is run as a partially stirred reactor with a low gas-phase mixing of one mixing event per numerical time step. Figure 9 shows the model predictions over the residence time at a constant temperature of 1023 K and Figure 10 the average emission concentrations over the recorded time (600 s) at various temperatures. In the experiment, the recording is started immediately when inserting the sample [44], while the emitted gases travel through the burner tube and to the analyzer. Hence, a signal delay is assumed for the comparison of experiments and simulations. Here, the delay is set to 80 s. Compared to the experiment, NO X , the sum of NO and NO 2 , is released over a too short time frame. The model resolves the characteristic two peaks; however, the second peak is advanced in time. In contrast, the model well captures the trend in NO X release over different temperatures. The release rate of SO 2 is well-replicated in terms of shape and timing (Figure 9) as is the average trend over temperatures ( Figure 10). With higher temperatures, the absolute amount is underpredicted. However, given the small absolute values, the prediction is accounted acceptable, and the reaction rates are not explicitly optimized for sewage sludge because of two aspects. First, due to the little available data in the literature and the small concentrations of the studied species, and hence typical afflicted measurement errors, only trends, not absolute numbers, are targeted by the model. Second, the model development aims to formulate a mathematical description for various biomass and waste streams. The presented model development is hence targeted to extend the capabilities of the municipal solid waste scheme [31] that was validated for the emission release in municipal solid waste and wood incineration [31,56]. Thus, the model's overall performance is prioritized over individual experiments in this context. However, future work will address the current limitations of the model.

Conclusions
The presented work introduces a surrogate formulation for the representation of sewage sludge in modeling. The solid surrogate species combine proteins, lignocellulosic species, sugar, lipid and species to account for the released gases by inorganic species. The devolatilization reactions are modified to comply with a detailed gas-phase mechanism, including tar and gas species reactions, emission pathways for NO X and SO X . It has been shown that the surrogates represent the ultimate composition of sewage sludge in a close agreement for a variety of different cases. Further, the amount of the individual substance classes, for example, protein and lignin, in the surrogate are plausible. Applied together with a zero-dimensional stochastic gasification reactor approach, the surrogate model enables us to predict the solid mass loss and the yield in the product gas, tar, and water. Further, the model captured the major gas composition and the release of emissions (NO X and SO 2 ) under pyrolysis, gasification, and incineration conditions. Future work will analyze the impact of fuel properties on the gas yield and efficiency of sewage sludge gasification as well as the emission formation in staged combustion of sewage sludge.

Data Availability Statement:
The devolatilization scheme is provided as Supplementary Material, and the data presented in this study are available within the article (tables and figures).