Loading [MathJax]/jax/element/mml/optable/GeneralPunctuation.js
Research article

Production of gasoline from municipal solid waste via steam gasification, methanol synthesis, and Methanol-to-Gasoline technologies: A techno-economic assessment

  • Received: 09 September 2020 Accepted: 29 November 2020 Published: 04 December 2020
  • A techno-economic analysis of the process of producing gasoline via steam gasification of municipal solid waste (MSW) using dolomite catalysts, followed by a methanol synthesis and a methanol-to-gasoline process (SGMG) was conducted using an Aspen Plus model and an economic assessment. The process includes six steps: MSW pretreatment, steam gasification of MSW to produce syngas, gas purification through CO2 and acid gas removal, methanol synthesis using syngas, conversion of methanol to gasoline, and gasoline separation. The Aspen Plus model used to estimate the energy and mass balance for sizing the equipment assumed that 2000 dry metric tonnes of MSW per day would be processed. Capital investment required and operating costs of gasoline production were estimated based on the Aspen Plus model's mass and energy balance and an nth plant. The minimum selling price (MSP) of gasoline was calculated using the capital investment and operating costs of the process. The total capital investment of the process was estimated to be $148 million, with an annual operational cost of $56 million. The minimum selling price was determined at $2.40/gal for a 20-year project life. The sensitivity analysis showed that the cost of feedstock has a direct impact on the MSP of gasoline, and MSP could decrease to $1.55/gal if the owners of the plant received disposal fees. Increasing the gasoline yield to 0.125 kg of gasoline/kg of MSW decreased the MSP to $1.86/gal. The ash content in MSW played a vital role in estimating the MSP. If the ash content in the feedstock increased, the MSP increased owing to a decrease in the gasoline yield and an increase in the cost of operations.

    Citation: Hamad AlMohamadi. Production of gasoline from municipal solid waste via steam gasification, methanol synthesis, and Methanol-to-Gasoline technologies: A techno-economic assessment[J]. AIMS Energy, 2021, 9(1): 50-67. doi: 10.3934/energy.2021004

    Related Papers:

    [1] Bin He . Developing a leap-frog meshless methods with radial basis functions for modeling of electromagnetic concentrator. AIMS Mathematics, 2022, 7(9): 17133-17149. doi: 10.3934/math.2022943
    [2] Cheng Chen . Hyperbolic function solutions of time-fractional Kadomtsev-Petviashvili equation with variable-coefficients. AIMS Mathematics, 2022, 7(6): 10378-10386. doi: 10.3934/math.2022578
    [3] Zihan Yue, Wei Jiang, Boying Wu, Biao Zhang . A meshless method based on the Laplace transform for multi-term time-space fractional diffusion equation. AIMS Mathematics, 2024, 9(3): 7040-7062. doi: 10.3934/math.2024343
    [4] Jing Li, Linlin Dai, Kamran, Waqas Nazeer . Numerical solution of multi-term time fractional wave diffusion equation using transform based local meshless method and quadrature. AIMS Mathematics, 2020, 5(6): 5813-5838. doi: 10.3934/math.2020373
    [5] Muhammad Imran Liaqat, Sina Etemad, Shahram Rezapour, Choonkil Park . A novel analytical Aboodh residual power series method for solving linear and nonlinear time-fractional partial differential equations with variable coefficients. AIMS Mathematics, 2022, 7(9): 16917-16948. doi: 10.3934/math.2022929
    [6] A. E. Matouk, Ismail Gad Ameen, Yasmeen Ahmed Gaber . Analyzing the dynamics of fractional spatio-temporal SEIR epidemic model. AIMS Mathematics, 2024, 9(11): 30838-30863. doi: 10.3934/math.20241489
    [7] Ali Khalouta, Abdelouahab Kadem . A new computational for approximate analytical solutions of nonlinear time-fractional wave-like equations with variable coefficients. AIMS Mathematics, 2020, 5(1): 1-14. doi: 10.3934/math.2020001
    [8] Hijaz Ahmad, Muhammad Nawaz Khan, Imtiaz Ahmad, Mohamed Omri, Maged F. Alotaibi . A meshless method for numerical solutions of linear and nonlinear time-fractional Black-Scholes models. AIMS Mathematics, 2023, 8(8): 19677-19698. doi: 10.3934/math.20231003
    [9] Abdul Samad, Imran Siddique, Fahd Jarad . Meshfree numerical integration for some challenging multi-term fractional order PDEs. AIMS Mathematics, 2022, 7(8): 14249-14269. doi: 10.3934/math.2022785
    [10] Abdul Samad, Imran Siddique, Zareen A. Khan . Meshfree numerical approach for some time-space dependent order partial differential equations in porous media. AIMS Mathematics, 2023, 8(6): 13162-13180. doi: 10.3934/math.2023665
  • A techno-economic analysis of the process of producing gasoline via steam gasification of municipal solid waste (MSW) using dolomite catalysts, followed by a methanol synthesis and a methanol-to-gasoline process (SGMG) was conducted using an Aspen Plus model and an economic assessment. The process includes six steps: MSW pretreatment, steam gasification of MSW to produce syngas, gas purification through CO2 and acid gas removal, methanol synthesis using syngas, conversion of methanol to gasoline, and gasoline separation. The Aspen Plus model used to estimate the energy and mass balance for sizing the equipment assumed that 2000 dry metric tonnes of MSW per day would be processed. Capital investment required and operating costs of gasoline production were estimated based on the Aspen Plus model's mass and energy balance and an nth plant. The minimum selling price (MSP) of gasoline was calculated using the capital investment and operating costs of the process. The total capital investment of the process was estimated to be $148 million, with an annual operational cost of $56 million. The minimum selling price was determined at $2.40/gal for a 20-year project life. The sensitivity analysis showed that the cost of feedstock has a direct impact on the MSP of gasoline, and MSP could decrease to $1.55/gal if the owners of the plant received disposal fees. Increasing the gasoline yield to 0.125 kg of gasoline/kg of MSW decreased the MSP to $1.86/gal. The ash content in MSW played a vital role in estimating the MSP. If the ash content in the feedstock increased, the MSP increased owing to a decrease in the gasoline yield and an increase in the cost of operations.


    The time-fractional partial differential equation (TFPDE) is one of the most widely studied fractional partial differential equations in mathematics and physics and has captured continued interests from researchers, due to its widespread applications [1]. For example, it is frequently used to model anomalous diffusion in porous media [2,3,4,5,6], where heterogeneities cause the diffusion rate to deviate from classical models. Additionally, the TFPDEs are employed to describe the propagation dynamics of light beams [7,8] and heat transfer in Maxwell fluids [9]. Therefore, it is very significant to study the properties of this kind of equation not only in practice but also in theory. In this paper, we consider the following TFPDEs with variable coefficients [10]:

    Dαtu(x,t)=di=1(ai(x,t)2ux2i+bi(x,t)uxi)+h(x,t),xΩ,t[t0,tf], (1.1)

    subjected to the initial and boundary conditions:

    u(x,t0)=f(x),xΩ, (1.2)
    u(x,t)=g(t),xΩ,t0<ttf, (1.3)

    where 0<α1, x is a spatial point, t is the time, t0 and tf represent the initial and final time, ΩRd, d=1,2 is a bounded domain with boundary Γ=Ω, ai(x,t), bi(x,t), h(x,t), f(x), and g(t) are known functions. Dαtu(x,t) is the time-fractional derivative, which has several different definitions [11,12]. Commonly used fractional derivatives include the Riemann-Liouville, Caputo, and Grunwald-Letnikov derivatives. In this study, we consider the following Caputo definition:

    Dαtu(x,t)={1Γ(1α)t0(ts)αu(x,s)sds,0<α<1,u(x,t)t,α=1. (1.4)
    Dαtu(x,t)={1Γ(2α)t0(ts)1α2u(x,s)s2ds,1<α<2,u2(x,t)t2,α=2. (1.5)

    The exact analytical solutions to the above TFPDEs are possible only for a few simple cases, but are virtually impossible for most problems, which greatly limits their applications. On the contrary, this enormously promotes the development of various numerical algorithms [13,14,15,16]. It is well known that the finite element, finite difference, boundary element, and meshless/meshfree methods are commonly used approaches for solving various integer/fractional order partial differential equations [17,18,19,20]. Among them, meshless methods have become quite popular in recent years thanks to their ease of implementation and interactive applications. Representative methods in this category include the Galerkin meshless method [21,22], the meshless local Petrov-Galerkin method [23,24], corrected smoothed particle hydrodynamics [25,26], the radial basis function collocation method [27,28], the generalized finite difference method [29,30,31], the backward substitution method [32,33], and localized semi-analytical collocation methods [34,35,36,37]. In most published works, however, these methods require the discretization of the temporal derivatives by using implicit, explicit, Runge-Kutta or any other approaches [38,39], which may be computationally expensive and sometimes mathematically complicated. To address these issues, numerous scholars in recent years have proposed space-time combined meshless methods [40,41,42,43,44] for solving various partial differential equations numerically.

    Inspired by previous work [41], this paper makes a first attempt to extend the recently proposed spatio-temporal meshless method (STMM) to the solutions of TFPDEs with variable coefficients. The STMM method first constructs a multiquadric radial basis function (ST-MQ-RBF) in terms of the space-time distance and then directly approximates the solution of Eqs (1.1)–(1.3) by a linear combination of the presented ST-MQ-RBFs at scattered nodes. Compared to the traditional RBF method, the main advantage of the developed approach lies in its complete avoidance of the discretization process for time-fractional-order derivatives, employing the same RBF approximation in both time and space. As a result, the method is simpler, more straightforward, easier to implement, and relatively more accurate.

    The outline of the paper is as follows: In Section 2, the STMM for spatial and temporal discretization simultaneously of TFPDEs is proposed, where the numerical approximation for time-fractional derivative is described. Section 3 briefly introduces a multiple-scale technique for solving the ill-conditioned system. In Section 4, numerical examples are studied and discussed. Finally, conclusions and remarks are presented in Section 5.

    The basic idea of STMM is to introduce a space-time distance into the MQ-RBF and then approximate the solutions of the partial differential equations by using the new ST-MQ-RBFs in the whole space-time domain without the differential approximation of time derivatives. For the above TFPDEs with variable coefficients, N=Ni+N0+Nb nodes (xi,ti)Ni=1 should be distributed on the space-time domain Ωst=Ω×[t0,tf] as shown in Figure 1, where Ni,N0,Nb represent the numbers of nodes in the domain, on the initial boundary and on the time boundary, respectively. For the ith node, the following MQ-RBF formulation should hold:

    Figure 1.  The distribution of nodes in the space-time domain Ωst=Ω×[t0,tf]: (a) one-dimension; (b) two-dimension.
    u(xi,ti)=Nj=1λjϕ(xi,ti;xj,tj), (2.1)

    where λj are the unknown coefficients, ϕ(xi,ti;xj,tj) are the MQ-RBFs given as follows:

    ϕ(xi,ti;xj,tj)=rij2+c2, (2.2)

    in which c is the shape parameter, rij denotes the space-time distance between the points (xi,ti) and (xj,tj), which can be written as

    rij=dk=1(xikxjk)2+(titj)2. (2.3)

    Like the traditional MQ-RBF method, the variables at interior points should satisfy the governing equation, and the variables at boundary points should satisfy initial and boundary conditions. With this regard, we have

    Nj=1λj(Dαtiϕ(xi,ti;xj,tj)dk=1(ak(x,t)2ux2k+bk(x,t)uxk))=h(xi,ti), i=1,...,Ni, (2.4)
    Nj=1λjϕ(xi,t0;xj,tj)=f(xi),xiΩ,i=Ni+1,...,Ni+N0, (2.5)
    Nj=1λjϕ(xi,ti;xj,tj)=g1(ti),xiΩ,i=Ni+N0+1,...,Ni+N0+Nb. (2.6)

    The spatial derivatives of the STRBFs in Eq (2.4) can be easily calculated. Here, we concern the calculation of the time-fractional derivative Dαtiϕ(xi,ti;xj,tj) in Eq (2.4). For α=1, it can be computed by the following expression:

    Dαtiϕ(xi,ti;xj,tj)=ϕ(xi,ti;xj,tj)ti=titjrij2+c2. (2.7)

    For 0<α<1, from the definition (1.4) of time-fractional derivative, using Gauss numerical integration, we have

    Dαtiϕ(xi,ti;xj,tj)=1Γ(1α)ti0(tis)αϕ(xi,s;xj,tj)sds=1Γ(1α)ti0(stj)(tis)αdk=1(xikxjk)2+(stj)2+c2ds
    =1Γ(1α)11((tiξ+ti)/2tj)(ti(tiξ+ti)/2)αti/2dk=1(xikxjk)2+((tiξ+ti)/2tj)2+c2dξ=1Γ(1α)Ngk=1ωkψ(ξk), (2.8)

    where ωk and ξk are the nodes and weights of the Gaussian quadrature, Ng is the number of Gaussian nodes (in this study, Ng is set as 16), and

    ψ(ξk)=((tiξk+ti)/2tj)(ti(tiξk+ti)/2)αti/2dk=1(xikxjk)2+((tiξk+ti)/2tj)2+c2. (2.9)

    For 1<α<2,

    Dαtiϕ(xi,ti;xj,tj)=1Γ(2α)ti0(tis)1α2ϕ(xi,s;xj,tj)s2ds=1Γ(2α)ti0(dk=1(xikxjk)2+2c2)(tis)1α(dk=1(xikxjk)2+(stj)2+c2)3/2ds=1Γ(2α)11(dk=1(xikxjk)2+2c2)(ti(tiξ+ti)/2)1αti/2(dk=1(xikxjk)2+((tiξ+ti)/2tj)2+c2)3/2dξ=1Γ(2α)Ngk=1ωkψ(ξk), (2.10)
    ψ(ξk)=(dk=1(xikxjk)2+2c2)(ti(tiξ+ti)/2)1αti/2(dk=1(xikxjk)2+((tiξ+ti)/2tj)2+c2)3/2. (2.11)

    By combining Eqs (2.4)–(2.6), a resultant linear system in matrix form can be generalized:

    Aλ=b, (2.12)

    where AN×N is the coefficient matrix, λ=[λ1,λ2,...,λN]T is a vector consisting of the undetermined coefficients; bN×1 is a vector consisting of the right-hand sides in Eqs (2.4)–(2.6). As long as Eq (2.12) is solved efficiently (that is to say, λ can be obtained accurately), the unknown variable at an arbitrary point (x,t) in the considered space-time box can be acquired by

    u(x,t)=Nj=1λjϕ(x,t;xj,tj). (2.13)

    For the RBF method, in general, the algebraic equation (2.12) is ill-conditioned. A pre-conditioner is required to reduce the condition number and to yield a regular solution. The multiple scale technique [45,46,47,48,49] is an effective pre-conditioner tool and thus is applied in this paper. Here, we will give a brief review of this method.

    Let A=[A1A2...AN], where {Ai}Ni=1 are column vectors of the matrix A. By the multiple scale technique, Eq (2.12) can be reformulated as follows:

    Aλ=b, (3.1)

    where

    A=[A1 (3.2)

    and

    \overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\smile}$}}{\boldsymbol{\lambda}} = {\left[ {{{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\smile}$}}{\lambda } }_1}\;\;{{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\smile}$}}{\lambda } }_2}\;\; \cdot \cdot \cdot \;\;{{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\smile}$}}{\lambda } }_N}} \right]^T} = {\left[ {{\lambda _1}{{\left\| {{A_1}} \right\|}_2}\;\;{\lambda _2}{{\left\| {{A_2}} \right\|}_2}\;\; \cdot \cdot \cdot \;\;{\lambda _n}{{\left\| {{A_N}} \right\|}_2}} \right]^T} . (3.3)

    The condition number of the matrix \overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\smile}$}}{\boldsymbol{A}} will be greatly reduced due to the reduction of round-off errors. After the system of Eq (3.1) is solved, the solution \boldsymbol{\lambda} = {({\lambda _1}, {\lambda _2}, ..., {\lambda _n})^T} of Eq (2.12) can be recovered from Eq (3.3) as follows:

    \boldsymbol{\lambda} = {\left[ {{{{{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\smile}$}}{\lambda } }_1}} \mathord{\left/ {\vphantom {{{{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\smile}$}}{\lambda } }_1}} {{{\left\| {{A_1}} \right\|}_2}}}} \right. } {{{\left\| {{A_1}} \right\|}_2}}}\;\;{{{{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\smile}$}}{\lambda } }_2}} \mathord{\left/ {\vphantom {{{{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\smile}$}}{\lambda } }_2}} {{{\left\| {{A_2}} \right\|}_2}}}} \right. } {{{\left\| {{A_2}} \right\|}_2}}}\;\; \cdot \cdot \cdot \;\;{{{{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\smile}$}}{\lambda } }_N}} \mathord{\left/ {\vphantom {{{{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\smile}$}}{\lambda } }_N}} {{{\left\| {{A_N}} \right\|}_2}}}} \right. } {{{\left\| {{A_N}} \right\|}_2}}}} \right]^T} . (3.4)

    It should be pointed out that Eq (3.1) is solved using the MATLAB command 'A\B'.

    In this section, two numerical examples including one- and two-dimensional TFPDEs with both regular and irregular geometries are studied by using regular and irregular nodes to illustrate the efficiency and accuracy of the proposed space-time meshless method. To evaluate the accuracy of the proposed STMM, the following L-error, L2-error and the root-mean-square error (RMSE) are adopted:

    \begin{gathered} {L_\infty }{\text{ - error}} = \mathop {\max }\limits_{1 \leqslant j \leqslant M} \left| {{u_{num}}({x_j}, {t_j}) - {u_{exa}}({x_j}, {t_j})} \right|, \;\;\;\; \hfill \\ {L_2}{\text{ - error}} = \sqrt {\sum\limits_{j = 1}^M {{{\left( {{u_{num}}({x_j}, {t_j}) - {u_{exa}}({x_j}, {t_j})} \right)}^2}} } , \hfill \\ RMSE = \sqrt {\frac{1}{M}\sum\limits_{j = 1}^M {{{\left( {{u_{num}}({x_j}, {t_j}) - {u_{exa}}({x_j}, {t_j})} \right)}^2}} } , \hfill \\ \end{gathered} (4.1)

    where M is the number of test points, {u_{exa}} and {u_{num}} represent the exact and numerical solutions at the jth test point, respectively.

    Example 1. Consider the following two one-dimensional TFPDEs with constant and variable coefficients on the domain {\Omega ^{st}} = [0, 1] \times [0, 1]:

    \mathit{Case I:}\; D_t^\alpha u(x, t) = \frac{{{\partial ^2}u(x, t)}}{{\partial {x^2}}} - \frac{{\partial u(x, t)}}{{\partial x}} + \frac{{2{t^{2 - \alpha }}}}{{\Gamma (3 - \alpha )}} + 2x - 2, (4.2)
    \mathit{Case II:}\; D_t^\alpha u(x, t) = ({x^2} + t + 1)\frac{{{\partial ^2}u(x, t)}}{{\partial {x^2}}} - {t^2}{e^x}\frac{{\partial u(x, t)}}{{\partial x}} + \frac{{2{t^{2 - \alpha }}}}{{\Gamma (3 - \alpha )}} - 2({x^2} + t) + 2{t^2}x{e^x} - 2. (4.3)

    The exact solution of the above two cases is u{\text{ }}(x, t) = {x^2} + {t^2} . The initial and boundary conditions of this problem are

    u(x, 0) = {x^2}, \;\;x \in [0, 1] , (4.4)
    u(0, t) = {t^2}, \;\;u(1, t) = 1 + {t^2}, \;t \in [0, 1] . (4.5)

    First, we consider the TFPDE with constant coefficients in Case I. Taking {N_i} interior nodes as test points, the numerical accuracy of the proposed STMM for Case I is tested by using different selection strategies of shape parameter c. Table 1 lists the {L_\infty } -errors and RMSEs under different values of \alpha , and N = 1840. For details of the three strategies mentioned in Table 1, we refer to [49]. As can be seen, the three techniques of shape parameter optimization are all efficient and accurate, and thus the first approach is adopted in the following calculations. It can also be observed that the accuracy of the algorithm decreases with the increase of \alpha . Moreover, the comparison between the numerical errors obtained by the present method with those reported in [34] is shown in Table 2, where we set \delta t = 0.01, N = 5151, {N_0} = 51, {N_t} = 2 \times 100, {N_i} = 49 \times 100 and \alpha = 0.5 for the comparison. The results indicate that the proposed scheme is more accurate than that in [50]. The number of nodes was kept constant, and Table 3 compares the condition numbers of the coefficient matrix with and without the use of multiscale techniques. It can be seen that the introduction of multiscale techniques significantly reduces the condition number. It is worth noting that the computational results for both cases are essentially the same, so they are not presented here.

    Table 1.  The errors for Case I under different shape parameter optimizations.
    c α=0.2 α=0.6
    L-error RMSE L-error RMSE
    Hardy [51] 4.629e-04 2.629e-04 5.101e-03 3.685e-03
    Fasshauer [52] 2.177e-04 1.129e-04 5.014e-03 3.634e-03
    Frank [53] 1.809e-04 1.018e-04 5.022e-03 3.640e-03

     | Show Table
    DownLoad: CSV
    Table 2.  The comparison of errors for the Case I under different times.
    t Method in Ref. [50] Present method
    L-error L2-error L-error L2-error
    0.1 6.086e-02 2.613e-01 1.709e-04 1.090e-04
    0.5 2.958e-02 1.277e-01 2.155e-03 1.559e-03
    1.0 2.114e-02 9.134e-02 6.688e-03 4.904e-03

     | Show Table
    DownLoad: CSV
    Table 3.  Condition numbers of the STMM with and without the multiple scale technique.
    Fractional orders Condition numbers
    With multiple scale technique Without multiple scale technique
    0.2 1.029e+09 1.264e+09
    0.5 1.036e+09 1.265e+09
    0.8 1.004e+09 1.234e+09

     | Show Table
    DownLoad: CSV

    Next, we consider the TFPDE with variable coefficients in Case II. The convergence of the proposed algorithm is investigated when \alpha = 0.2 and 0.6, and the absolute error curves are illustrated in Figure 2. The results demonstrate that the STMM has high precision and good convergence. It is also noted that the smaller the \alpha , the faster the convergence. As depicted in Figure 3, the graphs of the numerical solution and absolute error with \alpha = 0.5 further confirm the validity and accuracy of the present method for solving the TFPDE with variable coefficients.

    Figure 2.  Convergence curves of the numerical results for Case II with (a) \alpha = 0.2, and (b) \alpha = 0.6.
    Figure 3.  Surfaces of (a) numerical solution and (b) absolute error for Case II with \alpha = 0.5.

    Then, Figure 4 presents the errors with respect to the nodal spacing to show the convergence in space. It can be found that as the nodal spacing decreases, the algorithm demonstrates a clear convergence trend. When the nodal spacing becomes sufficiently small, the calculation accuracy stabilizes.

    Figure 4.  Errors of the STMM with respect to the nodal spacing (h) for Case II with \alpha = 0.5.

    Finally, we examined the impact of the number of Gaussian points on the accuracy of the algorithm. It can be seen from Table 4 that the algorithm shows a convergence trend as the number of Gauss points increases. Considering that the change in error is minimal, using 16 Gauss points is sufficient to ensure accuracy and save computation time.

    Table 4.  The errors for Case II with α = 0.5 under different numbers of Gaussian points.
    {N_g} 4 8 12 16 20
    L-error 1.192e-02 6.532e-03 4.664e-03 3.688e-03 3.097e-03
    RMSE 4.764e-03 2.548e-03 1.758e-03 1.359e-03 1.108e-03

     | Show Table
    DownLoad: CSV

    Example 2. Consider the following two-dimensional time fractional advection-diffusion equation on \Omega = [0, \; 1] \times [0, \; 1] [54]:

    D_t^\alpha u = \frac{{{\partial ^2}u}}{{\partial {x^2}}} + \frac{{{\partial ^2}u}}{{\partial {y^2}}} - \frac{{\partial u}}{{\partial x}} - \frac{{\partial u}}{{\partial y}} + 0.5\Gamma (3 + \alpha ){t^2}{e^{x + y}}, \;\;(x, y) \in \Omega , \;\;t > 0 , (4.6)

    with the following initial and boundary conditions:

    u(x, y, 0) = 0, \;\;(x, y) \in \Omega , (4.7)
    u(x, y, t) = {t^{2 + \alpha }}{e^{x + y}}, \;\;(x, y) \in \partial \Omega , \;\;t > 0 . (4.8)

    The exact solution of this problem can be given by

    u(x, y, t) = {t^{2 + \alpha }}{e^{x + y}} . (4.9)

    To investigate the influence of the node distribution on the numerical solution, we consider two types of node distribution (regular nodes and irregular nodes) as shown in Figure 5. When \alpha = 0.5 and a total of 4560 nodes are used, the analytical solution and absolute errors on the physical domain at t = 1.0 with regular and irregular nodal distributions are illustrated in Figure 6. These figures have proven that the newly proposed meshless approach has very good accuracy even using irregular nodal distributions.

    Figure 5.  Distributions of nodes for the rectangular domain: (a) regular nodal distribution; (b) irregular nodal distribution.
    Figure 6.  Numerical results of the STMM: (a) analytical solution and absolute errors for (b) regular nodal distribution and (c) irregular nodal distribution, where N = 4560, and \alpha = 0.5.

    Figure 7 presents the absolute errors at the point (0.5, 0.5, t) for the various fractional orders. Noted that the STMM has high computational accuracy, and the smaller the \alpha , the higher the accuracy, which is the same as one-dimensional problems. Furthermore, the convergence of the proposed meshless method is investigated in Table 5. As can be observed, the algorithm converges as the number of nodes increases.

    Figure 7.  Absolute errors at the point (0.5, 0.5, t) for different fractional orders obtained by the STMM with regular nodal distribution different and N = 4560.
    Table 5.  The errors for Example 2 under different number of nodes, where \alpha = 0.5.
    N 270 1540 4560 10080 18850
    L-error 4.908 e-02 3.708e-02 3.159e-02 2.758e-02 2.442e-02
    RMSE 3.703e-03 2.190e-03 1.675 e-03 1.422e-03 1.278e-03

     | Show Table
    DownLoad: CSV

    Finally, the proposed methodology is tested on an irregular domain \Omega shown in Figure 8. Using the same exact solution, initial and boundary conditions, we consider the following two-dimensional time fractional PDE with variable coefficients:

    Figure 8.  Distributions of nodes for the irregular domain in Example 2: (a) regular nodal distribution; (b) irregular nodal distribution.
    D_t^\alpha u = ({x^2} + t + 1)\frac{{{\partial ^2}u}}{{\partial {x^2}}} + ({y^2} - t - 1)\frac{{{\partial ^2}u}}{{\partial {y^2}}} + {t^2}{e^x}\frac{{\partial u}}{{\partial x}} + {t^2}{e^y}\frac{{\partial u}}{{\partial y}} + h(x, y, t), \;\;(x, y) \in \Omega , \;\;t > 0 , (4.10)

    where h(x, y, t) = 0.5\Gamma (3 + \alpha){t^2}{e^{x + y}} - ({x^2} + {y^2} + ({e^x} + {e^y}){t^2}){t^{(2 + \alpha)}}{e^{x + y}} . The boundary \partial \Omega of the computational domain is defined by the following parametric equation:

    \partial \Omega = \{ (x, y)|x = \rho (\theta )\cos (\theta ), \;\;y = \rho (\theta )\sin (\theta ), \;\;0 \leqslant \theta \leqslant 2\pi \} , (4.11)

    where \rho (\theta) = \frac{{3 + \cos (\theta - \pi /5)\sin (4\theta)}}{{5 + \sin (2\theta)}}.

    Consider the regular and irregular nodal distributions on the irregular geometry shown in Figure 8. The profiles of absolute errors t = 1.0 for these two cases of nodal distribution are depicted in Figure 9, where N = 2646 and \alpha = 0.5. It can be seen that the higher accuracy is achieved for different types of node distribution, although the regular nodes are slightly better than irregular nodes.

    Figure 9.  Profiles of absolute errors for Example 2: (a) regular nodal distribution; (b) irregular nodal distribution, where \alpha = 0.5 and N = 2646.

    A spatio-temporal meshless method is developed for the numerical solution of the TFPDEs with variable coefficients. The developed scheme can directly approximate the solutions of fractional partial differential equations within a space–time scale framework and is very simple mathematically, truly meshless, free of difference approximation, and computationally cost-effective. The numerical experiments further demonstrate the efficiency, precision, and convergence of the proposed STMM for the TFPDEs with variable coefficients.

    The non-local nature of fractional derivatives causes the solution's behavior to be influenced by the entire temporal or spatial domain, making numerical solutions inherently challenging. The introduction of variable coefficients further increases the problem's nonlinearity, rendering traditional numerical methods cumbersome and accuracy-limited. The method presented in this paper can efficiently approximate the numerical solution for such problems in a straightforward manner without the need for specialized time discretization techniques, making it applicable to various types of variable coefficient fractional differential equations. It provides a robust tool for solving time fractional PDEs with variable coefficients, thus advancing the field and offering new possibilities for tackling complex real-world problems in engineering and science.

    The present scheme can be extended to other fractional partial differential equations, such as the time-space fractional PDEs and complicated high-order and high-dimensional problems. Moreover, the study on local STMM is conducive to the solution of long-time and large-scale dynamic problems. It should be pointed out that the shape parameter in the STMM has a certain influence on the accuracy of the numerical solution, and thus the determination technique of the optimal shape parameter still needs to be discussed in depth. Otherwise, the new ST-RBFs should be studied in order to avert or greatly reduce the influence of shape parameters.

    Xiangyun Qiu: Conceptualization, Methodology, Simulation analysis, Writing–original draft; Xingxing Yue: Methodology, Project administration, Supervision. Both authors have read and approved the final version of the manuscript for publication.

    This work was supported by National Key R & D Program of China (Grant No. 2023YFB2503900).

    The authors declare no competing financial interests.



    [1] Verma M, Godbout S, Brar S, et al. (2012) Biofuels production from biomass by thermochemical conversion technologies. Int J Chem Eng 2012.
    [2] Ji L, Zhang C, Fang J (2017) Economic analysis of converting of waste agricultural biomass into liquid fuel: A case study on a biofuel plant in China. Renewable Sustainable Energy Rev 70: 224-229.
    [3] Sorda G, Banse M, Kemfert C (2010) An overview of biofuel policies across the world. Energy Policy 38: 6977-6988.
    [4] Onel O, Niziolek AM, Hasan MF, et al. (2014) Municipal solid waste to liquid transportation fuels—Part I: Mathematical modeling of a municipal solid waste gasifier. Comput Chem Eng 71: 636-647.
    [5] Easterly JL, Burnham M (1996) Overview of biomass and waste fuel resources for power production. Biomass Bioenergy 10: 79-92.
    [6] Sipra AT, Gao N, Sarwar H (2018) Municipal solid waste (MSW) pyrolysis for bio-fuel production: A review of effects of MSW components and catalysts. Fuel Process Technol 175: 131-147.
    [7] Matsakas L, Gao Q, Jansson S, et al. (2017) Green conversion of municipal solid wastes into fuels and chemicals. EJB 26: 69-83.
    [8] Le Courtois A (2012) Municipal Solid Waste: turning a problem into resource. Private Sector Dev 15: 1-28.
    [9] Magazzino C, Mele M, Schneider N (2020) The relationship between municipal solid waste and greenhouse gas emissions: Evidence from Switzerland. Waste Manage.
    [10] Magazzino C, Mele M, Schneider N, et al. (2020) Waste generation, wealth and GHG emissions from the waste sector: Is Denmark on the path towards circular economy? Sci Total Environ 755: 142510.
    [11] Park YJ, Heo J (2002) Vitrification of fly ash from municipal solid waste incinerator. J Hazard Mater 91: 83-93.
    [12] McKendry P (2002) Energy production from biomass (part 1): overview of biomass. Bioresour Technol 83: 37-46.
    [13] Huang H, Buekens A (2001) Chemical kinetic modeling of de novo synthesis of PCDD/F in municipal waste incinerators. Chemosphere 44: 1505-1510.
    [14] Malkow T (2004) Novel and innovative pyrolysis and gasification technologies for energy efficient and environmentally sound MSW disposal. Waste Manage 24: 53-79.
    [15] Abnisa F, Wan Daud WMA (2014) A review on co-pyrolysis of biomass: an optional technique to obtain a high-grade pyrolysis oil. Energy Convers Manage 87: 71-85.
    [16] Weil S, Hamel S, Krumm W (2006) Hydrogen energy from coupled waste gasification and cement production—a thermochemical concept study. Int J Hydrogen Energy 31: 1674-1689.
    [17] Zhang Y, Nagamori S, Hinchiranan S, et al. (2006) Promotional effects of Al2O3 addition to Co/SiO2 catalysts for Fischer-Tropsch synthesis. Energy Fuels 20: 417-421.
    [18] Chaudhari S, Dalai A, Bakhshi N (2003) Production of hydrogen and/or syngas (H2 + CO) via steam gasification of biomass-derived chars. Energy Fuels 17: 1062-1067.
    [19] Thamavithya M, Dutta A (2008) An investigation of MSW gasification in a spout-fluid bed reactor. Fuel Process Technol 89: 949-957.
    [20] Ponzio A, Kalisz S, Blasiak W (2006) Effect of operating conditions on tar and gas composition in high temperature air/steam gasification (HTAG) of plastic containing waste. Fuel Process Technol 87: 223-233.
    [21] Luo S, Zhou Y, Yi C (2012) Syngas production by catalytic steam gasification of municipal solid waste in fixed-bed reactor. Energy 44: 391-395.
    [22] Phillips SD, Tarud JK, Biddy MJ, et al. (2011) Gasoline from woody biomass via thermochemical gasification, methanol synthesis, and methanol-to-gasoline technologies: a technoeconomic analysis. Ind Eng Chem Res 50: 11734-11745.
    [23] AlMohamadi H, Gunukula S, DeSisto WJ, et al. (2018) Formate‐assisted pyrolysis of biomass: an economic and modeling analysis. Biofuels, Bioprod Biorefin 12: 45-55.
    [24] Hannula I (2016) Hydrogen enhancement potential of synthetic biofuels manufacture in the European context: A techno-economic assessment. Energy 104: 199-212.
    [25] Carrasco JL, Gunukula S, Boateng AA, et al. (2017) Pyrolysis of forest residues: An approach to techno-economics for bio-fuel production. Fuel 193: 477-484.
    [26] Phillips S (2007) Technoeconomic analysis of a lignocellulosic biomass indirect gasification process to make ethanol via mixed alcohols synthesis. Ind Eng Chem Res 46: 8887-8897.
    [27] Alamolhoda S, Vitale G, Hassan A, et al. (2019) Synergetic effects of cerium and nickel in Ce-Ni-MFI catalysts on low-temperature water-gas shift reaction. Fuel 237: 361-372.
    [28] Tijm P, Waller F, Brown D (2001) Methanol technology developments for the new millennium. Appl Catal A: General 221: 275-282.
    [29] Toyir J, Miloua R, Elkadri N, et al. (2009) Sustainable process for the production of methanol from CO2 and H2 using Cu/ZnO-based multicomponent catalyst. Physics Procedia 2: 1075-1079.
    [30] Yurchak S (1988) Development of Mobil's fixed-bed methanul-to-gasoline (MTG) process, In: Anonymous Studies in Surface Science and Catalysis: Elsevier, 251-272.
    [31] Kwak T, Lee S, Park J, et al. (2006) Gasification of municipal solid waste in a pilot plant and its impact on environment. Korean J Chem Eng 23: 954-960.
    [32] Swanson RM, Platon A, Satrio JA, et al. (2010) Techno-economic analysis of biomass-to-liquids production based on gasification. Fuel 89: S11-S19.
    [33] Posso F, Siguencia J, Sánchez J (2019) Use of municipal solid waste (MSW)-derived hydrogen in Ecuador: Potential applications for urban transportation. Waste Biomass Valorization 10: 1529-1537.
    [34] Larson ED, Jin H, Celik FE (2009) Large‐scale gasification‐based coproduction of fuels and electricity from switchgrass. Biofuels, Bioprod Biorefin 3: 174-194.
    [35] Ramirez JA, Rainey TJ (2019) Comparative techno-economic analysis of biofuel production through gasification, thermal liquefaction and pyrolysis of sugarcane bagasse. J Clean Prod 229: 513-527.
    [36] Enerkem, 2020. Available from: https://enerkem.com/products/methanol/.
    [37] Atalla TN, Gasim AA, Hunt LC (2018) Gasoline demand, pricing policy, and social welfare in Saudi Arabia: A quantitative analysis. Energy Policy 114: 123-133.
    [38] Radwan N, Mangi SA (2019) Municipal solid waste management practices and opportunities in Saudi Arabia. Eng, Technol Appl Sci Res 9: 4516-4519.
    [39] Available from: https://www.thefuelprice.com/, 2020.
  • This article has been cited by:

    1. Kolade M. Owolabi, Sonal Jain, Edson Pindza, Eben Mare, Comprehensive Numerical Analysis of Time-Fractional Reaction–Diffusion Models with Applications to Chemical and Biological Phenomena, 2024, 12, 2227-7390, 3251, 10.3390/math12203251
    2. Lin Qiu, Yanjie Wang, Yan Gu, Qing-Hua Qin, Fajie Wang, Adaptive physics-informed neural networks for dynamic coupled thermo-mechanical problems in large-size-ratio functionally graded materials, 2024, 0307904X, 115906, 10.1016/j.apm.2024.115906
    3. Lin Qiu, Fajie Wang, Wenzhen Qu, Ji Lin, Yan Gu, Qing‐Hua Qin, A Hybrid Collocation Method for Long‐Time Simulation of Heat Conduction in Anisotropic Functionally Graded Materials, 2025, 126, 0029-5981, 10.1002/nme.70002
    4. Lin Qiu, Fajie Wang, Yingjie Liang, Qing-Hua Qin, Physics-informed radial basis function network based on Hausdorff fractal distance for solving Hausdorff derivative elliptic problems, 2025, 183, 08981221, 271, 10.1016/j.camwa.2025.02.012
    5. Chih-Yu Liu, Cheng-Yu Ku, Wei-Da Chen, Ying-Fan Lin, Jun-Hong Lin, Solving Inverse Wave Problems Using Spacetime Radial Basis Functions in Neural Networks, 2025, 13, 2227-7390, 725, 10.3390/math13050725
    6. Hao Chang, Fajie Wang, Xingxing Yue, Lin Qiu, Linlin Sun, A 2.5D Generalized Finite Difference Method for Elastic Wave Propagation Problems, 2025, 13, 2227-7390, 1249, 10.3390/math13081249
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(6558) PDF downloads(624) Cited by(7)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog