
The frequency and severity of flooding in urban areas have escalated in recent years, and the worst affected urban areas are those in Africa. Despite the escalating flood risks accompanied by the growing vulnerability, cities and urban areas in Africa are struggling to build resilience. One of the actors in urban resilience building are the households in the urban settlements. Indeed, their contribution forms an important component of many flood risk response strategies. Nonetheless, the knowledge about this remains limited and is often confined to specific regions or case studies. In this study, we sought to identify the measures put in place to respond to flood risk by the households in Nyamasaria and Manyatta and explored the potential of mainstreaming resilience to flood risk in these settlements. We adopted a survey design. The settlements were purposively sampled, while the sampling procedure for the households involved transect lines established parallel to the major roads within the settlements. A transect walk was made to mark the households, and then simple random sampling was used. Data collection through personal interviews using questionnaires followed. Descriptive statistics were used to analyze the primary data. The findings showed that both settlements suffered frequent flood events, with 79 % of households having suffered inundation by floods in the past. About 46 % of the households have put in place flood risk reduction measures, including using sandbags, digging trenches around the houses, and raising floor levels. There were also cases of households temporarily relocating from the risk areas to safer places. The findings showed that the households' response measures were short-term and were aimed at addressing immediate risks. We concluded that the existence of capacities and actions to reduce flood risk among households provides a viable starting point for mainstreaming the resilience to flood risk in these settlements. We therefore recommend that households need to transition from short-term measures to risk-sensitive resilience measures, including flood-resilient designing of buildings, installation of risk-reducing infrastructure, and risk-sensitive urban planning.
Citation: Adoyo Laji, Jeremiah N. Ayonga. Mainstreaming resilience to flood risk among households in informal settlements in Kisumu City, Kenya[J]. Urban Resilience and Sustainability, 2024, 2(4): 326-347. doi: 10.3934/urs.2024017
[1] | A. Tomar, H. Kumar, M. Ali, H. Gandhi, D. Singh, G. Pathak . Application of symmetry analysis and conservation laws to a fractional-order nonlinear conduction-diffusion model. AIMS Mathematics, 2024, 9(7): 17154-17170. doi: 10.3934/math.2024833 |
[2] | Miao Yang, Lizhen Wang . Lie symmetry group, exact solutions and conservation laws for multi-term time fractional differential equations. AIMS Mathematics, 2023, 8(12): 30038-30058. doi: 10.3934/math.20231536 |
[3] | Khalid K. Ali, K. R. Raslan, Amira Abd-Elall Ibrahim, Mohamed S. Mohamed . On study the fractional Caputo-Fabrizio integro differential equation including the fractional q-integral of the Riemann-Liouville type. AIMS Mathematics, 2023, 8(8): 18206-18222. doi: 10.3934/math.2023925 |
[4] | Ishtiaq Ali, Saeed Islam . Stability analysis of fractional two-dimensional reaction-diffusion model with applications in biological processes. AIMS Mathematics, 2025, 10(5): 11732-11756. doi: 10.3934/math.2025531 |
[5] | Yumei Chen, Jiajie Zhang, Chao Pan . Numerical approximation of a variable-order time fractional advection-reaction-diffusion model via shifted Gegenbauer polynomials. AIMS Mathematics, 2022, 7(8): 15612-15632. doi: 10.3934/math.2022855 |
[6] | Xiaofeng Guo, Jianyu Pan . Approximate inverse preconditioners for linear systems arising from spatial balanced fractional diffusion equations. AIMS Mathematics, 2023, 8(7): 17284-17306. doi: 10.3934/math.2023884 |
[7] | Xingyang Ye, Chuanju Xu . A posteriori error estimates of spectral method for the fractional optimal control problems with non-homogeneous initial conditions. AIMS Mathematics, 2021, 6(11): 12028-12050. doi: 10.3934/math.2021697 |
[8] | Ndolane Sene . Fractional input stability for electrical circuits described by the Riemann-Liouville and the Caputo fractional derivatives. AIMS Mathematics, 2019, 4(1): 147-165. doi: 10.3934/Math.2019.1.147 |
[9] | Zhichao Fang, Ruixia Du, Hong Li, Yang Liu . A two-grid mixed finite volume element method for nonlinear time fractional reaction-diffusion equations. AIMS Mathematics, 2022, 7(2): 1941-1970. doi: 10.3934/math.2022112 |
[10] | M. Manigandan, Subramanian Muthaiah, T. Nandhagopal, R. Vadivel, B. Unyong, N. Gunasekaran . Existence results for coupled system of nonlinear differential equations and inclusions involving sequential derivatives of fractional order. AIMS Mathematics, 2022, 7(1): 723-755. doi: 10.3934/math.2022045 |
The frequency and severity of flooding in urban areas have escalated in recent years, and the worst affected urban areas are those in Africa. Despite the escalating flood risks accompanied by the growing vulnerability, cities and urban areas in Africa are struggling to build resilience. One of the actors in urban resilience building are the households in the urban settlements. Indeed, their contribution forms an important component of many flood risk response strategies. Nonetheless, the knowledge about this remains limited and is often confined to specific regions or case studies. In this study, we sought to identify the measures put in place to respond to flood risk by the households in Nyamasaria and Manyatta and explored the potential of mainstreaming resilience to flood risk in these settlements. We adopted a survey design. The settlements were purposively sampled, while the sampling procedure for the households involved transect lines established parallel to the major roads within the settlements. A transect walk was made to mark the households, and then simple random sampling was used. Data collection through personal interviews using questionnaires followed. Descriptive statistics were used to analyze the primary data. The findings showed that both settlements suffered frequent flood events, with 79 % of households having suffered inundation by floods in the past. About 46 % of the households have put in place flood risk reduction measures, including using sandbags, digging trenches around the houses, and raising floor levels. There were also cases of households temporarily relocating from the risk areas to safer places. The findings showed that the households' response measures were short-term and were aimed at addressing immediate risks. We concluded that the existence of capacities and actions to reduce flood risk among households provides a viable starting point for mainstreaming the resilience to flood risk in these settlements. We therefore recommend that households need to transition from short-term measures to risk-sensitive resilience measures, including flood-resilient designing of buildings, installation of risk-reducing infrastructure, and risk-sensitive urban planning.
Reaction-diffusion systems have attracted a considerable amount of attention in recent years. They arise naturally in various biological, chemical and physical models to describe the spatio-temporal concentration change of one or more species which involve both local reaction and diffusion simultaneously. The reaction-diffusion system consists of a set of partial differential equations (PDEs) to represent the behaviour of each species individually. Recently, the interest in the study of reaction-diffusion systems has been extended to fractional reaction-diffusion systems [1,2,3,4,5] which from one side exhibit self-organization phenomena and from the other side introduce a new parameter to these systems, which is a fractional derivative index. It gives a greater degree of freedom for diversity of self-organization phenomena and it can exhibit new types of solutions that are not possible to find in classical reaction-diffusion equations with integer derivatives [6,7,8]. From numerical point of view, several methods have been developed, including, for example, finite element schemes [9], finite volume [] and finite difference methods [11,12] and spectral ones [13,14].
In this paper, we consider a system of two coupled nonlinear time-fractional reaction-diffusion equations (TF-RDEs)
{∂αtu(t,x)−k1∂xxu(t,x)+f(t,x,u,v)=0,∂αtv(t,x)−k2∂xxv(t,x)+g(t,x,u,v)=0,0<α<1, | (1.1) |
where ∂αt is the Riemann-Liouville fractional derivative operator [15,16,17,18,19,20]
∂αtw(t,x)=1Γ(1−α)∂∂t∫t0w(s,x)(t−s)αds | (1.2) |
u(t,x) and v(x,t) being the field variables with t and x independent variables, k1>0 and k2>0 the diffusion coefficients, and f=f(t,x,u,v) and g=g(t,x,u,v) the nonlinear interaction terms. This kind of fractional problem has been numerically studied in a lot of papers, see, for example [21,22,23,24].
The main aim of this paper is to solve the system of reaction-diffusion equations (1.1) by applying the procedure recently proposed that combines the Lie symmetry analysis with the numerical methods to get exact and numerical solutions. This procedure was applied to a wide class of FDES: time-fractional advection–diffusion–reaction equations [25,26], space-fractional advection–diffusion–reaction equations with linear [27] and nonlinear sources terms [28], time-fractional and space-fractional advection-diffusion-reaction equations with linear and nonlinear sources terms [29], two dimensional time-fractional reaction-diffusion equations [30].
As well known, classical Lie symmetries allow to reduce the partial differential equations to ordinary differential equations, for this reason, recently the Lie symmetries theory has been extended to the FDES [31,32,33,34]. By using this extension of Lie symmetries approach to FPDEs, we find the Lie symmetries admitted by the model (1.1), then Lie transformations that map the coupled nonlinear FT-PDEs to a system of nonlinear fractional ordinary differential equations (FODEs). We find solutions of a system of two nonlinear time-fractional reaction-diffusion equations by solving the reduced FODEs such that the solutions of the reduced equations lead to obtain solutions of the original equations. The original model is written in terms of the Riemann-Liouville fractional derivative in order to determine the Lie symmetries by applying the MAPLE package, that automates the method of finding symmetries for FDEs as proposed in [31,32,33,34]. After the reduction of the model to a fractional ordinary differential equation, we introduce the Caputo derivative, and rewrite the FODE in terms of this derivative so that we are able to compute the numerical solutions. Classical finite difference methods are proposed with the aim of showing the simplicity of application of the procedure and to obtain, at the same time, highly accurate solutions and a low computational cost.
The plan of the paper is the follows: In Section 2, the Lie symmetries theory for FDEs is presented; in Section 3, we describe the reduction of the original model to FODEs and report some examples of mathematical models obtained by suitable choices of involved parameters and functions; the Section 4 is devoted to describe the numerical method used to obtain solutions of the reduced FODEs and, then, of the assigned original model. Finally, we end with the concluding remarks.
In general, for a 2×2 system of (integer order) partial differential equations
Δi(t,x,u,v,…,u(r1),v(r2))=0,i=1,2 | (2.1) |
where t,x are the independent variables, u and v are the dependent variables, and u(r1), v(r2) are all partial derivatives of the u and v with respect to t and x up to the maximum order r1 and r2 (integer order) respectively, we recall that the invertible transformations of the variables t,x,u,v
T=T(t,x,u,v,a),X=X(t,x,u,v,a), | (2.2) |
U=U(t,x,u,v,a),V=V(t,x,u,v,a) | (2.3) |
depending on a continuous parameter a, are said to be one-parameter (a) Lie point symmetry transformations of the Eq (2.1) if the Eq (2.1) has the same form in the new variables T,X,U,V.
The set G of all such transformations forms a continuous group, also known as the group admitted by the Eq (2.1).
According to the Lie theory, by expanding (2.3) in Taylor's series around a=0, we get the infinitesimal transformations
T=t+aξ1(t,x,u,v)+o(a),X=x+aξ2(t,x,u,v)+o(a), | (2.4) |
U=u+aη1(t,x,u,v)+o(a),V=v+aη2(t,x,u,v)+o(a), | (2.5) |
where their infinitesimals ξ1, ξ2, η1 and η2 are given by
ξ1(t,x,u,v)=∂T∂a|a=0,ξ2(t,x,u,v)=∂X∂a|a=0 |
η1(t,x,u,v)=∂U∂a|a=0,η2(t,x,u,v)=∂V∂a|a=0. |
The corresponding operator
Ξ=ξ1(t,x,u,v)∂t+ξ2(t,x,u,v)∂x+η1(t,x,u,v)∂u+η2(t,x,u,v)∂v | (2.6) |
is known in the literature as the infinitesimal operator or generator of the group G.
The point transformations leaving a differential equation (2.1) invariant are found by means of the straightforward Lie's algorithm, requiring that the k-order prolongation of the operator (2.6) acting on (2.1) is zero along the solutions, i.e.:
ΞkΔi=0|Δi=0, | (2.7) |
where k is the maximum order of derivations.
The invariance condition (2.7) provides an overdetermined set of linear differential equations (determining equations) for the infinitesimals whose integration gives the generators of Lie point symmetries admitted by the Eq (2.1).
Recently, in a series of papers of Buckwar, Luchko, Gazizov et al. [31,32,33,34,35], Lie symmetry methods have been extended to to FDEs. By extending transformations (2.5) to the operator of Riemann-Liouville fractional differentiation on derivatives of u with respect to t, Dαtuj, and v with respect to t, Dαtv, we have
Dαˉtu(ˉt,ˉx)=Dαtu(t,x)+aζ1α+o(a), | (2.8) |
Dαˉtv(ˉt,ˉx)=Dαtv(t,x)+aζ2α+o(a) | (2.9) |
where ζ1α and ζ2α the infinitesimals of fractional derivatives, are given by prolongation formula i.e. [32,36],
ζ1α=Dαt(η1)+ξ2Dαt(ux)−Dαt(ξ2ux)+Dαt(Dt(ξ1)u)−Dα+1t(ξ1u)+ξ1Dα+1t(u),ζ2α=Dαt(η2)+ξ2Dαt(vx)−Dαt(ξ2vx)+Dαt(Dt(ξ1)v)−Dα+1t(ξ1v)+ξ1Dα+1t(v) | (2.10) |
where Dt denotes the total derivative. Moreover, in order to conserve the structure of the fractional derivative operator, the following invariance condition is also required
ξ1(t,x,u)|t=0=0. |
Finally, the coefficients of the determining equation depend on all derivatives of variable u, v, Dαtu and Dαtv.
For the determination of the Lie point symmetries of (1.1), the authors apply an algorithm that has been implemented in the MAPLE package FracSym [37,38]; this algorithm uses some routines of the MAPLE symmetry packages DESOLVII [37] and ASP [38]; the routines in FracSym automate the method of finding symmetries for FDEs as proposed in [31,32,33,34]; these are the first routines to automate the search of symmetries for FDEs in MAPLE. A limit of the algorithm FracSym is that only the fractional derivatives of Riemann-Liouville type are considered. Moreover, it is important to note that the generator of the dependent variables u and v are assumed linear in u and v respectively, as usually done in the literature [33,34,35,36,37,38,39,40].
The extend Lie symmetries method leads to get the following infinitesimals for Eq (1.1)
ξ1=4a3t,ξ2=a0+2a3αx,η1=χ1(t,x)+(a1+2a3(α−1))u,η2=χ2(t,x)+a2v, | (3.1) |
where the function χ1=χ1(t,x) and χ2=χ2(t,x) satisfy the constraints
2a3(2t∂tf+αx∂xf+(α−1)u∂uf+(1+α)f)+∂αtχ1−k1∂xxχ1+a0∂xf+(χ1+a1u)∂uf+(χ2+a2v)∂vf−a1f=0,2a3(2t∂tg+αx∂xg+(α−1)u∂ug+2g)+∂αtχ2−k2∂xxχ2+a0∂xg+(χ1+a1u)∂ug+(χ2+a2v)∂vg−a2g=0, | (3.2) |
that link the arbitrary functions χ1=χ1(t,x) and χ2=χ2(t,x) to the nonlinear reaction terms f=f(t,x,u,v) and g=g(t,x,u,v), where ai, i=0,…,3, are the arbitrary parameters.
In the following, the stretching symmetry (i.e. a3) is neglected; then the constraints (3.2) reduce to
∂αtχ1−k1∂xxχ1+a0∂xf+(χ1+a1u)∂uf+(χ2+a2v)∂vf−a1f=0∂αtχ2−k2∂xxχ2+a0∂xg+(χ1+a1u)∂ug+(χ2+a2v)∂vg−a2g=0 | (3.3) |
and the infinitesimals read
ξ1=0,ξ2=a0,η1=χ1+a1u,η2=χ2+a2v. | (3.4) |
We obtain the transformation
T=t,U=u(t,x)e−a1x−∫e−a1xχ1dx,V=v(t,x)e−a2x−∫e−a2xχ2dx , | (3.5) |
where we choose a0=1. By integrating the constraints (3.3) with respect to x, and using the properties of the Riemann-Liouville fractional derivatives, the following forms of f and g are determined
f=ea1x(ϕ1−∫e−a1x(∂αtχ1−k1∂xxχ1)dx) | (3.6) |
g=ea2x(ϕ2−∫e−a2x(∂αtχ2−k2∂xxχ2)dx) | (3.7) |
where ϕ1=ϕ1(T,U,V) and ϕ2=ϕ2(T,U,V) are arbitrary functions of their arguments.
When the transformation (3.5) and the previous forms of f and g are inserted into the system (1.1), it is reduced into the following system of fractional nonlinear ordinary differential equations
DαTU(T)−k1a21U(T)+ϕ1(T,U,V)=0 | (3.8) |
DαTV(T)−k2a22V(T)+ϕ2(T,U,V)=0. | (3.9) |
The choice of the arbitrary functions ϕ1 and ϕ2 characterize the solutions of the Eqs (3.8) and (3.9) that, with a suitable choice of χ1 and χ2, specialize the source terms f and g and, then, the classes of solutions given by (3.5).
In order to perform some example of solutions of the system under study, we decide to choose ϕ1 and ϕ2 as nonlinear functions of two field variables U and V
ϕ1=c1U+c0UV+h1(T) | (3.10) |
ϕ2=c2V+c0UV+h2(T), | (3.11) |
that lead to get a reaction-diffusion model describes the interaction between two filed variables. The functions ϕ1 and ϕ2 characterize the source terms (3.6) and (3.7) of mathematical models that can be used in order to describe several natural phenomena where the fractional order temporal variation of the field variables involves both the reaction and diffusion processes, simultaneously. They have been applied in several contexts such as chemical reactions [41], genes behaviour [42], populations evolution [43], biological pattern formation [44], epidemics [45] and computer virus spreading [46]. Moreover, they are strictly valid to describe spatially distributed chemical and biochemical system, and can also be applied with considerable success to model non-molecular ensembles of interacting and diffusing objects [47]. The functions h1=h1(T) and h2=h2(T) are arbitrary ones of the time variable.
By these assumptions, the FODEs (3.8) and (3.9) read
DαTU(T)+(c1−k1a21)U(T)+c0U(T)V(T)+h1(T)=0 | (3.12) |
DαTV(T)+(c2−k2a22)V(T)+c0U(T)V(T)+h2(T)=0. | (3.13) |
The system (3.12)–(3.13) reduces to the following FODE
DαT(V(T)−U(T))+(c2−k2a22)(V(T)−U(T))+h2(T)−h1(T)=0 | (3.14) |
if we impose that
c1−k1a21=c2−k2a22. |
Under non-vanishing initial condition
limT→0Dα−1T(V(T)−U(T))=b0, | (3.15) |
and by the Laplace transform, it demonstrates that the FODE (3.14) admits the following exact solution (see [18] for more details)
V(T)=U(T)+b0Tα−1Eα,α((k1a21−c1)Tα)−∫T0(T−S)α−1Eα,α((k1a21−c1)(T−S)α)(h2(S)−h1(S))dS | (3.16) |
where Eα,α(t)=∞∑k=0tkΓ(α(k+1)) is the Mittag-Leffler function. Substituting in (3.12), we get the following FODE
DαTU(T)+c0U(T)2+((c1−k1a21)+c0b0Tα−1Eα,α((k1a21−c1)Tα)−c0∫T0(T−S)α−1Eα,α((k1a21−c1)(T−S)α)(h2(S)−h1(S))dS)U(T)+h1(T)=0. | (3.17) |
By solving the (3.17), we find the solution U(T) that leads to compute the solution V(T) by the (3.16) and, then, by the inverse transformations (3.5), the solutions of the proposed system of FPDEs (1.1) with source terms (3.6) and (3.7) characterized by (3.10) and (3.11). In general, Eq (3.17) could be not analytically resolvable, except for suitable choices of parameters and arbitrary functions, then numerical methods are proposed in order to compute the approximation of the exact solutions, as shown in the following Sections.
Starting from Eq (3.17) and setting c1=k1a21 and b0=0, we get the following FODE
DαTU(T)+c0U2(T)−c0Γ(α)(∫T0(T−S)α−1(h2(S)−h1(S))dS)U(T)+h1(T)=0. | (3.18) |
Suitable choices of the functions h1(T) and h2(T), involved in (3.18), are proposed with the aim to perform examples of mathematical models of interest in many fields of the applied sciences.
Example 1: By following choice of the functions h1(T) and h2(T)
h1(T)=−h0c0(1Γ(α)+1)eλT,h2(T)=h0c0(1Γ(α)−1)eλT, |
with h0 and λ arbitrary constants, the Eq (3.18) reads
DαTU(T)+c0U2(T)−2h0Γ(α)−Γ(α,λT)λαΓ(α)2eλTU(T)−h0c0(1Γ(α)+1)eλT=0. | (3.19) |
Computed U(T), the solution V(T), given in (3.16), reads
V(T)=U(T)−2h0c0Γ(α)−Γ(α,λT)λαΓ(α)2eλT. | (3.20) |
Finally, (3.10) and (3.11) assume the following form
ϕ1(T,U,V)=k1a21U+c0UV−h0c0(1Γ(α)+1)eλTϕ2(T,U,V)=k2a22V+c0UV+h0c0(1Γ(α)−1)eλT. |
Found U(T) and V(T), by inverse transformations (3.5), we get the solutions u(t,x) and v(t,x) of target model (1.1) with source terms obtained by inserting ϕ1 and ϕ2 in (3.6) and (3.7).
Remark 1: For α=1, we are able to solve the FODE (3.19) that assumes the following form
U′(T)+c0U2(T)+2h0λ(1−eλT)U(T)−2h0c0eλT=0, |
whose exact solution, when U(0)=0, is
U(T)=2h0eλT−1c0λ | (3.21) |
and substituting in (3.16), we obtain
V(T)=0. |
By inverse transformation (2.3), assumed χ1(x,t)=χ2(x,t)=0, we get the solutions
u(t,x)=2h0ea1xeλt−1c0λ,v(t,x)=0, |
of the system of PDEs with the source terms
f(t,x,u,v)=c1u+c0e−a2xuv−2h0c0eλt+a1x,g(t,x,u,v)=c2v+c0e−a1xuv. |
The exact solution of the model with α=1 is considered in the following numerical examples as a reference for testing the qualitative behavior of the numerical solution of the model (3.19) for values of the α parameter increasing towards 1.
Example 2: By following choice of the functions h1(T) and h2(T)
h1(T)=h0Γ(α−1)Tα+λ24c0eλT(eλT+21−Γ(α)Γ(α)), |
h2(T)=h0Γ(α−1)Tα+λ24c0eλT(eλT−21+Γ(α)Γ(α)), |
with h0 and λ arbitrary constants, the Eq (3.18) reads
DαTU(T)+c0U2(T)+λ2−αΓ(α)−Γ(α,λt)Γ(α)2eλtU(T)+h0Γ(α−1)Tα+λ24c0eλt(eλt+21−Γ(α)Γ(α))=0. | (3.22) |
Computed the solution U(T), the solution V(T), given in (3.16), reads
V(T)=U(T)+λ2−αΓ(α)−Γ(α,λt)c0Γ(α)2eλt | (3.23) |
and (3.10) and (3.11) assume the following form
ϕ1(T,U,V)=k1a21U+c0UV+h0Γ(α−1)Tα+λ24c0eλT(eλT+21−Γ(α)Γ(α))ϕ2(T,U,V)=k2a22V+c0UV+h0Γ(α−1)Tα+λ24c0eλT(eλT−21+Γ(α)Γ(α)). |
Found U(T) and V(T), by inverse transformations (3.5), we get the solutions u(t,x) and v(t,x) of target model (1.1) with source terms obtained by inserting ϕ1 and ϕ2 in (3.6) and (3.7).
Remark 2: For α=1, the FODE (3.22) reduces to
U′(T)+c0U2(T)+λ(eλt−1)U(T)+λ24c0e2λt=0, |
whose the exact solution, when U(0)=b1, is
U(T)=λeλT2c0(2(2b1c0+λ)(2b1c0+λ)eλT−2b1c0+λ−1). | (3.24) |
Substituting in (3.16), it has
V(T)=λ2c0(2(2b1c0−λ)(2b1c0+λ)eλT−2b1c0+λ+eλT), |
and, by the inverse transformation, we obtain the solutions
u(t,x)=λeλt+a1x2c0(2(2b1c0+λ)(2b1c0+λ)eλt−2b1c0+λ−1), |
v(t,x)=ea2xλ2c0(2(2b1c0−λ)(2b1c0+λ)eλt−2b1c0+λ+eλt), |
of the system of PDEs with the following source terms
f(t,x,u,v)=k1a21u+c0e−a2xuv+λ24c0e2λt+a1x,g(t,x,u,v)=k2a22v+c0e−a1xuv+λ24c0e2λt+a2x. |
The exact solution of the model with α=1 is considered in the following numerical examples as a reference for testing the qualitative behavior of the numerical solution of the model (3.22) for values of the α parameter increasing towards 1.
In this Section, we find the numerical solution of the system of FPDEs (1.1) computed by solving the following FODE
DαTU(T)+c0U2(T)−c0Γ(α)(∫T0(T−S)α−1(h2(S)−h1(S))dS)U(T)+h1(T)=0, |
obtained in the previous Section by the suitable choice of the parameters. Computed the numerical solution U(T), we obtain the solution V(T) by (3.16) and, then the solutions u(t,x) and v(t,x) of the target model (1.1) by the inverse transformations (3.5).
We introduce the α-order Caputo fractional derivative of the solution U(T)
CDαTU(T)=1Γ(1−α)∫T0(T−S)−αU′(S)dS |
and its connection with the α-order Riemann-Liouville fractional derivative
CDαTU(T)=DαT(U(T)−U(0)), |
with U(0) initial condition. We remark that the Caputo definition of the fractional derivative allows to define a FIVP whose the initial conditions are given in terms of the field variable and its integer order derivatives in agreement with the clear physical meaning of most of the processes arising in the real world. In terms of the Caputo derivative, the following fractional initial value problem (FIVP) is obtained
CDαTU(T)=F(T,U),U(0)=U0, | (4.1) |
where
F(T,U)=−c0U2(T)+(c0Γ(α)∫T0(T−S)α−1(h2(S)−h1(S))dS)U(T)−h1(T)−U0Γ(α−1)Tα. |
In order to numerically solve the FIVP (4.1), the classical implicit trapezoidal method (PI2 Im) is used. We built a computational uniform mesh of grid points denoted by Tn, with Tn=T0+nΔT and integration step sizes ΔT and N positive integer. We denote by Un the numerical approximation provided by the numerical method of the exact solution U(Tn) at the mesh points Tn, for n=0,⋯,N. The numerical method reads
Un+1=U0+1Γ(α)(β0F(T0,U0)+n+1∑k=1βkF(Tk,Uk)), | (4.2) |
where the coefficient values βk, for k=0,⋯,n+1, are computed as follows
β0=1α(α+1)(Tn+1)α((T1−T0)(α+1)−Tn+1)+(Tn+1−T1)α+1T1−T0,βk=1α(α+1)×k=1,⋯,n,×((Tn+1−Tk−1)α+1−(Tn+1−Tk)α+1Tk−Tk−1−(Tn+1−Tk)α+1−(Tn+1−Tk+1)α+1Tk+1−Tk),βn+1=1α(α+1)(Tn+1−Tn)α . |
The convergence order of the scheme is O((ΔT)min(1+α,2)). Note that, the convergence order of the trapezoidal method usually is 1+α when 0<α<1, and only when α>1 or when the solution is sufficiently smooth we obtain the reached order two [48]. In general, the numerical method (4.2) leads to obtain a nonlinear equation for whose resolution a root-finding solver is needed. The classical Newton method is proposed.
In the follows, we report some test problems with the aim to validate the proposed approach. Note that the proposed mathematical models represent a wide class of mathematical ones that can be used in order to describe several natural phenomena arising in many field of the applied sciences. Different simulation parameters, initial conditions and fractional α parameter values are considered as the input features of the proposed models. The exact solutions, reported in Remarks 1 and 2, are considered as references for testing qualitative behavior of the numerical solutions of the models (3.19) and (3.22) for values of the α parameter increasing towards 1. All numerical simulations are performed on Intel Core i7 by using Matlab 2020 software.
Example 1: In this example, we recall the FIVP (3.19), obtained in the previous Section, and we rewrite it in term of the Caputo fractional derivative choosing U0=0
CDαTU(T)+c0U2(T)−2h0Γ(α)−Γ(α,λT)λαΓ(α)2eλTU(T)−h0c0(1Γ(α)+1)eλT=0,U(0)=0. | (4.3) |
We set the parameters values as follows: c0=5, h0=1 and λ=−1 and choose a computational domain [0,Tmax] with Tmax=10 and N=100 grid points. In Figure 1, we report the numerical results obtained for different values of the α parameter: in the left frame, the numerical solution Un obtained by solving the FIVP (4.3) by using the PI2 Im numerical method; in the right frame, the numerical solution Vn (3.20) obtained by applying the (3.16). The black lines represent the exact solution U(T) given by (3.21) of the model with α=1 and the exact solution V(T) computed by the (3.16). They are reported in order to test the qualitative behavior of the numerical solutions as the α parameter increases towards 1. The solutions reveal a similar behaviour: Both starts from the value 0, both increase at the beginning of the integration process and, after, decreases as the time evolves. Note that, as the value of α increases, the solution U(T) increases unlike the solution V(T) decreases.
In Figure 2, we report the numerical solutions unj and vnj, approximations of the exact solutions obtained by the inverse transformations (3.5),
u(t,x)=U(t)ea1x,v(t,x)=V(t)ea2x, |
solutions of the model (1.1) with source terms
f(t,x,u,v)=k1a21u+c0e−a2xuv−h0c0(1Γ(α)+1)eλt+a1x,g(t,x,u,v)=k2a22v+c0e−a1xuv+h0c0(1Γ(α)−1)eλt+a2x, |
computed by setting k1=k2=0.5, for α=0.5, with c1=k1a21 and c2=k2a22 with a1=a2=−1, on a computational domain [0,Tmax]×[0,Xmax] with Tmax=10, Xmax=1 and N=J=100 grid points.
Example 2: In this example, we recall the FIVP (3.22), obtained in the previous Section, and we rewrite it in term of the Caputo fractional derivative choosing U0=b1=−h0
CDαTU(T)+c0U2(T)+λ2−αΓ(α)−Γ(α,λt)Γ(α)2eλtU(T)+λ24c0eλt(eλt+21−Γ(α)Γ(α))=0,U(0)=U0. | (4.4) |
We choose the initial condition U0=1 and set the parameters values as follows: c0=0.5 and λ=−1 and consider a computational domain [0,Tmax] with Tmax=10 and N=100 grid points. In Figure 3, we report the numerical results obtained for different values of α parameter: in the left frame, the numerical solution Un obtained by solving the FIVP (4.4) by using the PI2 Im numerical method; in the right frame, the numerical solution Vn (3.23) obtained by the (3.16). The black lines represent the exact solution U(T) given by (3.24) of the model with α=1 and the exact solution V(T) computed by the (3.16). They are reported in order to test the qualitative behavior of the numerical solutions as the α parameter increases towards 1. The solutions reveal a different behaviour: both starts from the value 1 but, the solution U(T) immediately decreases as the time evolves, unlike the solution V(T) increases at the beginning of the integration process and, after, decreases. As the value of α increases, the solution U(T) decreases unlike the solution V(T) increases.
In Figure 4, we report the numerical solutions unj and vnj, approximations of the exact solutions obtained by the inverse transformations (3.5), u(t,x)=U(t)ea1x and v(t,x)=V(t)ea2x, solutions of the model (1.1) with source terms
f(t,x,u,v)=k1a21u+c0e−a2xuv+h0Γ(α−1)Tα+λ24c0eλt+a1x(eλt+21−Γ(α)Γ(α)),g(t,x,u,v)=k2a22v+c0e−a1xuv+h0Γ(α−1)Tα+λ24c0eλt+a2x(eλt−21+Γ(α)Γ(α)), | (4.5) |
computed by setting k1=k2=0.5, for α=0.5, with c1=k1a21 and c2=k2a22 with a1=a2=−1, on a computational domain [0,Tmax]×[0,Xmax] with Tmax=10, Xmax=1 and N=J=100 grid points.
For the second test, we set the parameters values as before and c0=5 in order to provide an example of how this parameter, that represent the rate of interaction between the two solutions, affects their behaviour as the time evolves. In Figure 5, we report the numerical results obtained for different values of α parameter: in the left frame, the numerical solution Un obtained by solving the FIVP (4.4) by using the PI2 Im numerical method, in the right frame, the numerical solution Vn (3.23) obtained by the (3.16). The black lines represent the exact solution U(T) given by (3.24) of the model with α=1 and the exact solution V(T) computed by the (3.16). The solutions reveal a very similar behaviour: both starts from the value 1 and both immediately decrease as the time evolves. As the values of α increases, the solution U(T) decreases. A slightly different behaviour appears for the solution V(T) in the last part of the domain.
In Figure 6, we report the numerical solutions unj and vnj, approximations of the exact solutions obtained by the inverse transformations (3.5), solutions of the model (1.1) with source terms given by (4.5), computed by setting k1=k2=0.5, for α=0.5, with c1=k1a21 and c2=k2a22 with a1=a2=−1, on a computational domain [0,Tmax]×[0,Xmax] with Tmax=10, Xmax=1 and N=J=100 grid points.
This paper deals with the solutions of coupled nonlinear time-fractional reaction–diffusion equations, where the fractional derivatives are defined in terms of the Riemann-Liouville ones. We apply the strategy that combines an analytical and numerical approach in order to solve the problem under study, used recently in [25,26,27,28,29,30]. By Lie symmetries, the assigned system of two FPDEs is reduced into a system of two nonlinear FODEs and, under suitable assumptions of the involved parameters and functions, it is reduced into a nonlinear FODE. By introducing the Caputo derivative, we are able to properly handle the obtained FIVP and, then, by the trapezoidal numerical method, to find the solutions. Finally, the inverse transformations deal to compute the solutions of the original model. We remark that the Caputo definition of the fractional derivative allows to define a FIVP whose the initial condition is given in terms of the field variable in agreement with most of the processes with a clear physical meaning. Numerical examples are performed in order to show how the proposed approach good works with nonlinear systems of FPDEs confirming the clear advantage, from numerical point of view, to use a combined strategy: we numerically solve only a FODE whose solution allows to compute the solution of the original model of FPDEs. We remark that the presented test problems, performed in order to confirm the efficiency of the proposed approach, represent a wide class of mathematical models that can be used in order to describe several natural phenomena arising in many field of the applied sciences.
A. J. acknowledges the financial support by G.N.C.S. of I.N.d.A.M.
M.P.S. acknowledges the financial support by G.N.F.M. of I.N.d.A.M.
The authors declare no conflict of interest.
[1] |
Prashar N, Lakra HS, Shaw R, et al. (2023) Urban flood resilience: A comprehensive review of assessment methods, tools, and techniques to manage disaster. Prog Disaster Sci 20: 100299. https://doi.org/10.1016/j.pdisas.2023.100299 doi: 10.1016/j.pdisas.2023.100299
![]() |
[2] |
Neves JL (2024) Urban planning for flood resilience under technical and financial constraints: The role of planners and competence development in building a flood-resilient city in Matola, Mozambique. City Environ Interact 22: 100147. https://doi.org/10.1016/j.cacint.2024.100147 doi: 10.1016/j.cacint.2024.100147
![]() |
[3] | Porio E (2011) Vulnerability, adaptation, and resilience to floods and climate change-related risks among marginal, riverine communities in Metro Manila. Asian J Soc Sci 39: 425–445. Available from: https://www.jstor.org/stable/43498807. |
[4] | Laji A, Ayonga JN, Daudi F (2017) The interplay between urban development patterns and vulnerability to flood risk in Kisumu city, Kenya. J Environ Earth Sci 7: 92–99. |
[5] |
Cobbinah PB (2021) Urban resilience in climate change hotspot. Land Use Policy 100: 104948. https://doi.org/10.1016/j.landusepol.2020.104948 doi: 10.1016/j.landusepol.2020.104948
![]() |
[6] | Bahadur AV, Tanner T (2022) Resilient Reset: Creating Resilient Cities in the Global South, London: Routledge. https://doi.org/10.4324/9780429355066 |
[7] |
Brody SD, Lee Y, Highfield WE (2017) Household adjustment to flood risk: A survey of coastal residents in Texas and Florida, United States. Disasters 41: 566–586. https://doi.org/10.1111/disa.12216 doi: 10.1111/disa.12216
![]() |
[8] | Cao Y, Wilkinson E, Pettinotti L, et al. (2021) An Analytical Review: A Decade of Urban Resilience, New York: UNDP. Available from: https://www.undp.org/publications/analytical-review-decade-urban-resilience. |
[9] |
Kareem B, Lwasa S, Tugume D, et al. (2020) Pathways for resilience to climate change in African cities. Environ Res Lett 15: 073002. https://doi.org/10.1088/1748-9326/ab7951 doi: 10.1088/1748-9326/ab7951
![]() |
[10] |
Kreibich H, Bubeck P, Van Vliet M, et al. (2015) A review of damage-reducing measures to manage fluvial flood risks in a changing climate. Mitig Adapt Strateg Glob Change 20: 967–989. http://doi.org/10.1007/s11027-014-9629-5 doi: 10.1007/s11027-014-9629-5
![]() |
[11] | van Valkengoed AM, Steg L (2019) Climate change adaptation by individuals and households: A psychological perspective. University of Groningen. Available from: https://gca.org. |
[12] | Lavell A, Oppenheimer M, Diop C, et al. (2012) Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation, Cambridge: Cambridge University Press, 25–64. |
[13] | Melissa P, Ebalu O (2012) City Resilience in Africa: A Ten Essentials Pilot, New York: UNISDR. |
[14] | Mulligan J, Harper J, Kipkemboi P, et al. (2016) Community-responsive adaptation to flooding in Kibera, Kenya', in Proceedings of the institution of civil engineers-engineering sustainability, 170: 268–280. https://doi.org/10.1680/jensu.15.00060 |
[15] |
Ncube S, Wilson A, Petersen L, et al. (2023) Understanding resilience capitals, agency and habitus in household experiences of water scarcity, floods and fire in marginalized settlements in the Cape Flats, South Africa. Soc Sci Humanit Open 8: 100710. https://doi.org/10.1016/j.ssaho.2023.100710 doi: 10.1016/j.ssaho.2023.100710
![]() |
[16] | Coaffee J, Lee P (2016) Urban Resilience: Planning for Risk, Crisis and Uncertainty, London: Palgrave. |
[17] | Cutter SL (2021) Urban risks and resilience, In: Urban Informatics, Singapore: Springer. https://doi.org/10.1007/978-981-15-8983-6-13 |
[18] |
Sharifi A, Yamagata Y (2014) Resilient urban planning: Major principles and criteria. Energy Procedia 61: 1491–1495. https://doi.org/10.1016/j.egypro.2014.12.154 doi: 10.1016/j.egypro.2014.12.154
![]() |
[19] |
Twum KO, Abubakari M (2019) Cities and floods: A pragmatic insight into the determinants of households' coping strategies to floods in informal Accra, Ghana. Jàmbá: J Disaster Risk Stud 11: 608. https://doi.org/10.4102/jamba.v11i1.608 doi: 10.4102/jamba.v11i1.608
![]() |
[20] |
Birkmann J (2011) First-and second-order adaptation to natural hazards and extreme events in the context of climate change. Nat Hazard 58: 811–840. https://doi.org/10.1007/s11069-011-9806-8 doi: 10.1007/s11069-011-9806-8
![]() |
[21] | Pelling M (2011) Adaptation to Climate Change: From Resilience to Transformation, London: Routledge. https://doi.org/10.4324/9780203889046 |
[22] |
Forsyth T, Evans N (2013) What is autonomous adaption? Resource scarcity and smallholder agency in Thailand. World Dev 43: 56–66. https://doi.org/10.1016/j.worlddev.2012.11.010 doi: 10.1016/j.worlddev.2012.11.010
![]() |
[23] |
Ahsan MN (2017) Can strategies to cope with hazard shocks be explained by at-risk households' socioeconomic asset profile? Evidence from tropical cyclone-prone coastal Bangladesh. Int J Disaster Risk Sci 8: 46–63. https://doi.org/10.1007/s13753-017-0119-8 doi: 10.1007/s13753-017-0119-8
![]() |
[24] |
Schaer C (2015) Condemned to live with one's feet in water? A case study of community-based strategies and urban maladaptation in flood prone Pikine/Dakar, Senegal. Int J Clim Change Strategies Manage 7: 534–551. https://doi.org/10.1108/IJCCSM-03-2014-0038 doi: 10.1108/IJCCSM-03-2014-0038
![]() |
[25] |
Chhetri N, Stuhlmacher M, Ishtiaque A (2019) Nested pathways to adaptation. Environ Res Commun 1: 015001. https://doi.org/10.1088/2515-7620/aaf9f9 doi: 10.1088/2515-7620/aaf9f9
![]() |
[26] | Birkman J, Tetzlaff G, Zentel KO (2009) Addressing the Challenge: Recommendations and Quality Criteria for Linking Disaster Risk Reduction and Adaptation to Climate Change, Bonn: DKKV Publication Series. |
[27] |
Bulkeley H, Tuts R (2013) Understanding urban vulnerability, adaptation and resilience in the context of climate change. Local Environ 18: 646–662. https://doi.org/10.1080/13549839.2013.788479 doi: 10.1080/13549839.2013.788479
![]() |
[28] | Glavovic BC, Smith GP (2014) Adapting to Climate Change: Lessons from Natural Hazards Planning, Dordrecht: Springer. |
[29] |
Grothman T, Patt A (2005) Adaptive capacity and human cognition: The process of individual adaptation to climate change. Global Environ Change 15: 199–213. https://doi.org/10.1016/j.gloenvcha.2005.01.002 doi: 10.1016/j.gloenvcha.2005.01.002
![]() |
[30] | Ashwill M, Heltberg R (2014) Is there a Community Level Adaptation Deficit?, Washington, DC: Word Bank. |
[31] |
Haque AN, Dodman D, Hossain MM (2014) Individual, communal and institutional responses to climate change by low-income households in Khulna, Bangladesh. Environ Urban 26: 112–129. https://doi.org/10.1177/0956247813518681 doi: 10.1177/0956247813518681
![]() |
[32] |
Limthongsakul S, Nitivattananon V, Arifwidodo SD (2017) Localized flooding and autonomous adaptation in peri-urban Bangkok. Environ Urban 29: 51–68. https://doi.org/10.1177/0956247816683854 doi: 10.1177/0956247816683854
![]() |
[33] |
Kiunsi R (2013) The constraints on climate change adaptation in a city with a large development deficit: The case of Dar es Salaam. Environ Urban 25: 321–337. https://doi.org/10.1177/0956247813489617 doi: 10.1177/0956247813489617
![]() |
[34] |
John R (2020) Flooding in informal Settlements: Potentials and limits for household adaptation in Dar es Salaam City, Tanzania. Am J Clim Change 9: 68–86. https://doi.org/10.4236/ajcc.2020.92006 doi: 10.4236/ajcc.2020.92006
![]() |
[35] | Owusu K, Obour PB (2021) Urban flooding, adaptation strategies, and resilience: Case study of Accra, Ghana, In: African Handbook of Climate Change Adaptation, Cham: Springer International Publishing, 2387–2403. https://doi.org/10.1007/978-3-030-45106-6_249 |
[36] |
Amoako C (2017) Emerging grassroots resilience and flood responses in informal settlements in Accra, Ghana. GeoJournal 83: 949–965. https://doi.org/10.1007/s10708-017-9807-6 doi: 10.1007/s10708-017-9807-6
![]() |
[37] |
Hunter NB, North MA, Roberts DC, et al. (2020) A systematic map of responses to climate impacts in urban Africa. Environ Res Lett 15: 103005. https://doi.org/10.1088/1748-9326/ab9d00 doi: 10.1088/1748-9326/ab9d00
![]() |
[38] |
Okunola OH, Simatele MD, Olowoporoku O (2022) The influence of socioeconomic factors on individual and household adaptation strategies to climate change risks in Port Harcourt, Nigeria. J Integr Environ Sci 19: 273–288. https://doi.org/10.1080/1943815X.2022.2143821 doi: 10.1080/1943815X.2022.2143821
![]() |
[39] |
Twerefou DK, Adu-Danso E, Abbey E, et al. (2019) Choice of household adaptation strategies to flood risk management in Accra, Ghana. City Environ Interact 3: 100023. https://doi.org/10.1016/j.cacint.2020.100023 doi: 10.1016/j.cacint.2020.100023
![]() |
[40] |
Dilling L, Daly ME, Travis WR, et al. (2015) The dynamics of vulnerability: Why adapting to climate variability will not always prepare us for climate change. Wiley Interdiscip Rev Clim Change 6: 413–425. https://doi.org/10.1002/wcc.341 doi: 10.1002/wcc.341
![]() |
[41] |
Zevenbergen C, Gersonius B, Radhakrishan M (2020) Flood resilience. Phil Trans R Soc A 378: 1–7. https://doi.org/10.1098/rsta.2019.0212 doi: 10.1098/rsta.2019.0212
![]() |
[42] |
Satterthwaite D, Archer D, Colenbrander S, et al. (2020) Building resilience to climate change in informal settlements. One Earth 2: 143–156. https://doi.org/10.1016/j.oneear.2020.02.002 doi: 10.1016/j.oneear.2020.02.002
![]() |
[43] | Anyumba G. Kisumu town: History of the built form, planning and environment: 1890–1990, PhD Dissertation, Technische Universiteit van Delf, Netherlands, 1995. |
[44] |
Mireri C, Atekyereza P, Kyessi A, et al. (2007) Environmental risks of urban agriculture in the Lake Victoria drainage basin: A case of Kisumu municipality, Kenya. Habitat Int 31: 375–386. https://doi.org/10.1016/j.habitatint.2007.06.006 doi: 10.1016/j.habitatint.2007.06.006
![]() |
[45] | Associated Programme on Flood Management (APFM) (2008) APFM Technical Document–11: Urban flood risk management. |
[46] | Kenya National Bureau of Statistics (2019) 2019 Kenya population and housing census volume Ⅱ: Distribution of population by administrative units. Available from: https://www.knbs.or.ke/2019-kenya-population-and-housing-census-reports. |
[47] |
Noll B, Filatova T, Need A (2022) One and done? Exploring linkages between households' intended adaptations to climate‐induced floods. Risk Anal 42: 2781–2799. https://doi.org/10.1111/risa.13897 doi: 10.1111/risa.13897
![]() |
[48] |
Liao KH (2012) A theory on urban resilience to floods—a basis for alternative planning practices. Ecol Soc 17: 48. http://doi.org/10.5751/ES-05231-170448 doi: 10.5751/ES-05231-170448
![]() |
[49] | UN Habitat (2005) Situation analysis of informal settlements in Kisumu. Available from: https://unhabitat.org/situation-analysis-of-informal-settlements-in-kisumu. |
1. | Li Tian, Ziqiang Wang, Junying Cao, A high-order numerical scheme for right Caputo fractional differential equations with uniform accuracy, 2022, 30, 2688-1594, 3825, 10.3934/era.2022195 | |
2. | Xumei Zhang, Junying Cao, A high order numerical method for solving Caputo nonlinear fractional ordinary differential equations, 2021, 6, 2473-6988, 13187, 10.3934/math.2021762 | |
3. | P. Prakash, K.S. Priyendhu, M. Lakshmanan, Nonlinear two-component system of time-fractional PDEs in (2+1)-dimensions: Invariant subspace method combined with variable transformation, 2024, 137, 10075704, 108123, 10.1016/j.cnsns.2024.108123 | |
4. | Alessandra Jannelli, Maria Paola Speciale, 2023, Chapter 6, 978-981-19-7715-2, 91, 10.1007/978-981-19-7716-9_6 | |
5. | P. Prakash, K.S. Priyendhu, M. Lakshmanan, Generalized separable solutions for (2+1) and (3+1)-dimensional m-component coupled nonlinear systems of PDEs under three different time-fractional derivatives, 2025, 191, 09600779, 115852, 10.1016/j.chaos.2024.115852 |