Abstract
Initial-boundary value problems for the one-dimensional linear advection–dispersion equation with decay (LAD) are studied by utilizing a unified method, known as the Fokas method. The method takes advantage of the spectral analysis of both parts of Lax pair and the global algebraic relation coupling all initial and boundary values. We present the explicit analytical solution of the LAD equation posed on the half line and a finite interval with general initial and boundary conditions. In addition, for the case of periodic boundary conditions, we show that the solution of the LAD equation is asymptotically t-periodic for large t if the Dirichlet boundary datum is periodic in t. Furthermore, it can be shown that if the Dirichlet boundary value is asymptotically periodic for large t, then so is the unknown Neumann boundary value, which is uniquely characterized in terms of the given asymptotically periodic Dirichlet boundary datum. The analytical predictions for large t are compared with numerical results showing the excellent agreement.
1 Introduction
In this paper, we study initial-boundary value problems for the one-dimensional linear advection–dispersion equation with decay (LAD) [1]
where D > 0, c ≥ 0 and rc ≥ 0 represents the dispersion, advection coefficients and the first order decay rate. This equation is also known as the convection–diffusion equation with a reaction, according to scientific contexts [2]. The LAD equation is a combination of the linear dispersion and advection equations plus the decay, but which is also an analytical model, describing physical, chemical or biological diffusive transport phenomena in science and engineering [3], [4], [5], [6]. For example, it describes a mathematical model for one-dimensional contaminant transport in porous medium systems with first-order decay. The LAD equation without the decay term was formally solved in [7], [8]. Analytical solutions of the LAD equation, including the decay term, have been developed in several literatures [1], [3], [6] (see also [9] and references therein). It should be noted that most proposed analytical solutions are for the LAD equation with relatively limited initial and boundary conditions, such as homogeneous, constant or exponential boundary conditions in time. Moreover, these solutions can be involved with the complementary error functions. A few explicit analytical solutions are available if more general or complicated initial and boundary conditions are prescribed.
We propose here explicit analytical solutions for initial-boundary value problems of the LAD equation with general boundary conditions by utilizing a unified method, also known as the Fokas method [10], [11] (see the monograph [12] and the pedagogical literature of the method [13]). The Fokas method has been introduced to analyze boundary value problems for nonlinear integrable systems, which can be considered as a significant extension of the inverse scattering transform. Moreover, it has been shown that the method can be extensively applicable to a large class of partial differential equations (PDEs); for example, nonlinear integrable systems and linear evolution equations [11], [14], [15], [16] and linear and nonlinear elliptic PDEs [17], [18], [19], [20], [21] as well as linear and integrable nonlinear discrete equations [22], [23].
The presented method has several advantages. For linear cases, the Fokas method is relatively simple, but very effective, to implement. The main steps of the Fokas method can be summarized as follows: (i) simultaneous analysis of the both parts of Lax pair, which can be considered a novel type of separability [24]; (ii) analysis of the global relation which is an algebraic equation that involves all initial and boundary values. This global relation can be used to determine the unknown boundary values, known as the generalized Dirichlet-to-Neumann map [25] (see also [26], [27] for the recent application of the global relation). Moreover, the Fokas method presents an integral representation of the solution with explicit exponential (x, t)-dependence. Thus, it enables to characterize the asymptotic behaviors of the solution and the unknown boundary values for large t [28], [29], [30], [31], [32], [33].
The purpose of the work is to present new explicit analytical solutions for initial-boundary value problems of the LAD equation with general initial and boundary conditions. More specifically, letting
Note that one can reduce LAD equation (2) to the LAD equation without decay term by changing variable
2 The LAD equation on the half line
In Section 2, we first study the Dirichlet boundary value problem for the LAD equation posed on the half line
with the initial condition
Also, we denote the unknown Neumann boundary value by
The LAD equation can be written as an overdetermined linear system, known as a Lax pair
where
Indeed, the LAD equation (2) is the compatible condition of the Lax pair equation (4) in the sense that
where
and
Then, taking
Let
Applying
where
Employing the inverse Fourier transform in equation (8), the solution q(x, t) can be recovered as
Since the integrand in the second integral of equation (10) is analytic and bounded in the region, where Im k > 0 and Reω(k) > 0, by the Cauchy theorem, we deform the contour (−∞, ∞) to
Note that the representation of the solution equation (11) involves the unknown function
Note that if
Substituting equation (12) into equation (9), we find
The right-hand-side of equation (13) contains the unknown function
Remark 2.1.
We can derive the unknown Neumann boundary value directly from equation (12) as
Indeed, the effective portion of
Note that letting
Then, we deform the contour
The representation of the solution given in equation (14) has the explicit (x, t)-dependence of the exponential form. Thus, it is possible to study the appropriate asymptotics of the solution for large t. Below, we will show that the solution for the LAD equation with a periodic Dirichlet boundary datum is asymptotically periodic for large t. The similar result of the Proposition 2.1 also can be found in [33] for the linearized Korteweg–de Vries equation.
Proposition 2.1.
Let q(x, t) be the solution given by equation (14). Assume that
where
Proof.
Note that
We will estimate each integral in the above equations. For
Then, using the following identity,
we find
for some constant M1 > 0. For
From the identity,
it follows that
for some constant M2 > 0.
Regarding
where
Note that since
Then using equation (19), we find
In order to derive equation (17), let
where
Hence, we find
For
Since
which implies that by equation (22),
for some constant M3 > 0. Thus, the series given in equation (24) converges and so does the sequence
Then, taking n → ∞ in the right-hand-side of the above equation, we find
Next, we study the asymptotically periodic boundary data for large t. We assume that the Dirichlet boundary value is asymptotically periodic for large t in the sense that [29]
where
We first note that by integration by parts, we write the function
Deforming the contour
where we have used the identity: by the residue theorem,
Then, differentiating equation (25) with respect to x and evaluating the resulting equation at x = 0, the function
Proposition 2.2.
Assume that
where
and its Fourier series is given by
with
where
and
Proof.
Let
Note that using integration by parts, we find
which implies that
By the Cauchy theorem, we deform the contour
Thus, we know that
Regarding the integral
Then, substituting the Fourier expansion (27) into (31), we find
For each
Note that as k → ∞
Thus, by the residue theorem, we evaluate the first integral in equation (33) as
For the second integral in equation (33), we use the steepest descent method for asymptotics of integrals. Note that −ω(k) has a critical point at k = −ib/2 and
Hence, we introduce the steepest descent contour CS, where CS is the line passing through the point k = −ib/2 and parallel to the real k-axis (depicted in Figure 1). Deforming the contour
Therefore,
Remark 2.2.
It should be noted that the Fokas method can be applied effectively well to solve boundary problems with more general boundary conditions, such as the Neumann and Robin boundary conditions, namely, (i) Neumann boundary:
Neumann boundary condition. Solving equation (12) for
Thus, substituting equation (34) into equation (9), we find
Robin boundary condition. The Robin boundary value can be expressed in terms of the transforms
where
and then
Using the above equations in equation (9), the effective portion for
3 The LAD equation on a finite interval
In Section 3, we discuss the Fokas method to solve the LAD equation posed on a finite interval, namely,
We denote the initial and boundary values as
We assume that the boundary values
and
Taking
where
Applying the inverse Fourier transform, we find
By the Cauchy theorem, we deform the contour (−∞, ∞) to
Note that it requires to determine the unknown functions
where
Replacing k → −k − ib in equation (41), we find
We solve the linear system of equations (41) and (42) for
where
Remark 3.1.
Since the function
the functions
4 Examples
In Section 4, we consider the following Dirichlet boundary value problem
with the given initial and boundary values
Example 4.1.
It is well-known that the function
We will derive the form of the solution directly from the above conditions. In terms of the transforms, we find
where
where we have deformed the contour (−∞, ∞) to
Example 4.2.
We consider the case
where C is a real constant. The function
which yields
The integral in equation (49) can not be computed explicitly, so we evaluate integral numerically. To this end, we use the efficient analytical-numerical scheme proposed in [34]. By the Cauchy theorem, we deform the contour
Note that along this contour L, the exponential terms in the integral decay rapidly for large k and hence numerical integration converges quickly. The accuracy and order of convergence for this hybrid analytical-numerical scheme has been discussed in [34]. Analytical and numerical solutions are shown in Figure 2 with C = 1. The analytical solution in Figure 2 (as blue solid curve) has been obtained by integrating equation (49) numerically along the contour L via the adaptive Gauss-Kronrod quadrature method with
Example 4.3.
We consider the case of the periodic boundary value
In this case, the solution is given by
where
The integral in equation (51) can be computed by a similar way as discussed in Example 4.2. Figure 3a, b show that the analytical solution (shown in the figure as blue solid curve) is asymptotically periodic for large t and is again in excellent agreement with the numerical results (marked in the figure as red open dots).
Example 4.4.
We consider the case that the Dirichlet boundary value is asymptotically periodic, namely,
From equation (29), it follows that
where
with
Example 4.5.
We consider the LAD equation posed on the finite interval,
with
In this case,
where
Using change of variable k → −k − i for the second integral, the solution q(x, t) can be written as
The solution q(x, t) given in equation (55) is shown in Figure 5(a). We also have displayed the solution curves in Figure 5(b) as blue solid curves for x = 0.1, 0.5 and 0.8, which are in excellent agreement with the numerical solutions (red dashed curve).
5 Concluding remarks
We have demonstrated the Fokas method to solve initial-boundary value problems for the LAD equation posed on the half line and a finite interval with general boundary conditions. In addition to solving the LAD equation, we have discussed the case of periodic boundary conditions, which commonly appears in physics and engineering. In particular, we have characterized the long time asymptotic behaviors of the solution and the unknown boundary value, which can be uniquely determined by the known asymptotically periodic boundary datum. These analytical predictions have been compared with numerical results showing the excellent agreement.
It should be noted that the presented method is relatively simple, but remarkably effective, to implement, finding an explicit integral representation of the solution. The method works equally well for boundary value problems for linear and nonlinear integrable systems. In contrast to classical transform methods such as the Fourier and Laplace transforms, the Fokas method has several advantages. The integral representation of the solution involves explicit (x, t)-dependence of the exponential form, and hence it allows to study the long time asymptotics of the solution. Also, from the integral representation of analytic functions, it makes possible to compute effectively the numerical solution [12], [20], [34]. More importantly, the method can be nonlinearizable in the sense that it is successfully used for analyzing nonlinear integrable PDEs as well as nonlinear lattices [35]. Recently, the method has been extended to study well-posedness, regularity and controllability of PDEs [36], [37], [38].
The LAD equation can be considered as a generalization of the heat equation and hence we expect that it could be possible to study boundary value problems with more complicated boundary data as shown in [39], [40], [41], [42], [43]. Also, note that the explicit representations of the solutions in finite and infinite domains can be used to analyze the effect of the exit and inlet boundary conditions for small Péclet number as discussed in [44]. We will address regarding these issues in the near future.
Funding source: Daegu University
Award Identifier / Grant number: 2017
Acknowledgments
The author thanks anonymous referees for their useful comments and constructive suggestions. The work is supported by the Daegu University Research Grant 2017.
Author contribution: All the authors have accepted responsibility for the entire content of this submitted manuscript and approved submission.
Research funding: This study was funded by Daegu University Research Grant 2017.
Conflict of interest statement: The authors declare no conflicts of interest regarding this article.
The inversion formula for the t-transform. The inversion formula for the spectral function equation (7) can be derived by performing spectral analysis of equation (4b). We consider the following spectral problem
where
Note that
where
where
Letting k → −k − ib in the integral over
Thus, the reconstruction formula for g(t) can be found as
Note that we can replace
vanishes by the Cauchy theorem, since t − s < 0 and the integrand is analytic and bounded in D+. Therefore, we find the following inversion formula for the spectral function
References
[1] M. T. van Genuchten and W. J. Alves, Analytical solutions of the one-dimensional convective-dispersive solute transport equation, Washington US Department of Agriculture, Technical Bulletins 157268, United States Department of Agriculture, Economic Research Service, 1982.Search in Google Scholar
[2] J. Zhong, C. Zeng, Y. Yuan, Y. Zhang, and Y. Zhang, “Numerical solution of the unsteady diffusion-convection-reaction equation based on improved spectral Galerkin method,” AIP Adv., vol. 8, p. 045314, 2018, https://doi.org/10.1063/1.5023332.Search in Google Scholar
[3] J. S. Pérez Guerrero, E. M. Pontedeiro, M. T. van Genuchten, and T. H. Skaggs, “Analytical solution of the one-dimensional advection–dispersion solute transport equation subject to time-dependent boundary conditions,” Chem. Eng. J., vol. 221, pp. 487–491, 2013, https://doi.org/10.1016/j.cej.2013.01.095.Search in Google Scholar
[4] I. Boztosun and A. Charafi, “An analysis of the linear advection–diffusion equation using mesh-free and mesh-dependent methods,” Eng. Anal. Bound. Elem., vol. 26, pp. 889–895, 2002, https://doi.org/10.1016/s0955-7997(02)00053-x.Search in Google Scholar
[5] A. Mojtabi and M. O. Deville, “One-dimensional linear advection–diffusion equation: analytical and finite element solutions,” Comput. Fluids, vol. 107, pp. 189–195, 2015, https://doi.org/10.1016/j.compfluid.2014.11.006.Search in Google Scholar
[6] A. Moranda, R. Cianci, and O. Paladino, “Analytical solutions of one-dimensional contaminant transport in soils with source production-decay,” Soil Syst., vol. 2, no. 40, pp. 1–16, 2018, https://doi.org/10.3390/soilsystems2030040.Search in Google Scholar
[7] D. H. Tang, E. O. Frind, and E. A. Sudicky, “Contaminant transport in fractured porous media: analytical solution for a single fracture,” Water Resour. Res., vol. 17, pp. 555–564, 1981, https://doi.org/10.1029/wr017i003p00555.Search in Google Scholar
[8] G. Dagan, “Theory of solute transport by groundwater,” Annu. Rev. Fluid Mech., vol. 19, pp. 183–215, 1987, https://doi.org/10.1146/annurev.fl.19.010187.001151.Search in Google Scholar
[9] J. S. Pérez Guerrero, L. C. G. Pimentel, T. H. Skaggs, and M. T. van Genuchten, “Analytical solution of the advection–diffusion transport equation using a change-of-variable and integral transform technique,” Int. J. Heat Mass Tran., vol. 52, pp. 3297–3304, 2009, https://doi.org/10.1016/j.ijheatmasstransfer.2009.02.002.Search in Google Scholar
[10] A. S. Fokas, “A unified transform method for solving linear and certain nonlinear PDEs,” Proc. Math. Phys. Eng. Sci., vol. 453, pp. 1411–1443, 1997, https://doi.org/10.1098/rspa.1997.0077.Search in Google Scholar
[11] A. S. Fokas, “Integrable nonlinear evolution equations on the half-line,” Commun. Math. Phys., vol. 230, pp. 1–39, 2002, https://doi.org/10.1007/s00220-002-0681-8.Search in Google Scholar
[12] A. S. Fokas, “A unified approach to boundary value problems,” in CBMS-NSF Regional Conf. Series in Applied Mathematics, Philadelphia, SIAM, 2008, https://doi.org/10.10.1137/1.9780898717068.10.1137/1.9780898717068Search in Google Scholar
[13] B. Deconinck, T. Trogdon, and V. Vasan, “The method of Fokas for solving linear partial differential equations,” SIAM Rev., vol. 56, pp. 159–186, 2014, https://doi.org/10.1137/110821871.Search in Google Scholar
[14] A. S. Fokas, “On the integrability of certain linear and nonlinear partial differential equations,” J. Math. Phys., vol. 41, pp. 4188–4237, 2000, https://doi.org/10.1063/1.533339.Search in Google Scholar
[15] A. S. Fokas, “A new transform method for evolution partial differential equations,” IMA J. Appl. Math., vol. 67, pp. 559–590, 2002, https://doi.org/10.1093/imamat/67.6.559.Search in Google Scholar
[16] A. S. Fokas and B. Pelloni, “A transform method for linear evolution PDEs on a finite interval,” IMA J. Appl. Math., vol. 70, pp. 564–587, 2005, https://doi.org/10.1093/imamat/hxh047.Search in Google Scholar
[17] M. J. Colbrook, N. Flyer, and B. Fornberg, “On the Fokas method for the solution of elliptic problems in both convex and non-convex polygonal domains,” J. Comput. Phys., vol. 374, pp. 996–1016, 2018, https://doi.org/10.1016/j.jcp.2018.08.005.Search in Google Scholar
[18] A. S. Fokas, “Two dimensional linear PDEs in a convex polygon,” Proc. Math. Phys. Eng. Sci., vol. 457, pp. 371–393, 2001, https://doi.org/10.1098/rspa.2000.0671.Search in Google Scholar
[19] B. Pelloni and D. A. Pinotsis, “The elliptic sine-Gordon equation in a half plane,” Nonlinearity, vol. 23, pp. 77–88, 2010, https://doi.org/10.1088/0951-7715/23/1/004.Search in Google Scholar
[20] M. J. Colbrook, A. S. Fokas, and P. Hashemzadeh, “A hybrid analytical-numerical technique for elliptic PDEs,” SIAM J. Sci. Comput., vol. 41, pp. A1066–A1090, 2019, https://doi.org/10.1137/18m1217309.Search in Google Scholar
[21] M. J. Colbrook, “Extending the unified transform: curvilinear polygons and variable coefficient PDEs,” IMA J. Numer. Anal., vol. 40, pp. 976–1004, 2020, https://doi.org/10.1093/imanum/dry085.Search in Google Scholar
[22] G. Biondini and G. Hwang, “Initial-boundary value problems for discrete evolution equations: discrete linear Schrödinger and integrable discrete nonlinear Schrödinger equations,” Inverse Probl., vol. 24, pp. 1–44, 2008,065011, https://doi.org/10.1088/0266-5611/24/6/065011.Search in Google Scholar
[23] B. Moon and G. Hwang, “Discrete linear evolution equations in a finite lattice,” J. Differ. Equ. Appl., vol. 25, pp. 630–646, 2019, https://doi.org/10.1080/10236198.2019.1613386.Search in Google Scholar
[24] A. S. Fokas, “Lax pair: a novel type of separability,” Inverse Probl., vol. 25, p. 123007, 2009, https://doi.org/10.1088/0266-5611/25/12/123007.Search in Google Scholar
[25] A. S. Fokas, “The generalized Dirichlet-to-Neumann map for certain nonlinear evolution PDEs,” Commun. Pure Appl. Math., vol. LVIII, pp. 639–670, 2005, https://doi.org/10.1002/cpa.20076.Search in Google Scholar
[26] D. G. Crowdy and E. Luca, “Solving Wiener-Hopf problems without kernel factorization,” Proc. Math. Phys. Eng. Sci., vol. 470, p. 20140304, 2014, https://doi.org/10.1098/rspa.2014.0304.Search in Google Scholar
[27] M. J. Colbrook, L. J. Ayton, and A. S. Fokas, “The unified transform for mixed boundary condition problems in unbounded domains,” Proc. R. Soc. A., vol. 475, p. 20180605, 2019, https://doi.org/10.1098/rspa.2018.0605.Search in Google Scholar
[28] A. B. de Monvel, A. Kotlyarov, D. Shepelsky, and C. Zheng, “Initial boundary value problems for integrable systems: towards the long-time asymptotics,” Nonlinearity, vol. 23, p. 2483, 2010, https://doi.org/10.1088/0951-7715/23/10/007.Search in Google Scholar
[29] J. Lenells and A. S. Fokas, “The unified method on the half-line: II. NLS on the half-line with t-periodic boundary conditions,” J. Phys. Math. Theor., vol. 45, p. 195202, 2012, https://doi.org/10.1088/1751-8113/45/19/195202.Search in Google Scholar
[30] G. Hwang and A. S. Fokas, “The modified Korteweg–de Vries equation on the half-line with a sine-wave as Dirichlet datum,” J. Nonlinear Math. Phys., vol. 20, pp. 135–157, 2013, https://doi.org/10.1080/14029251.2013.792492.Search in Google Scholar
[31] J. Lenells and A. S. Fokas, “The nonlinear Schrödinger equation with t-periodic data: II. Perturbative results,” Proc. Math. Phys. Eng. Sci., vol. 471, p. 20140926, 2015, https://doi.org/10.1098/rspa.2014.0926.Search in Google Scholar
[32] G. Hwang, “The modified Korteweg–de Vries equation on the quarter plane with t-periodic data,” J. Nonlinear Math. Phys., vol. 24, pp. 620–634, 2017, https://doi.org/10.1080/14029251.2017.1375695.Search in Google Scholar
[33] B. Moon and G. Hwang, “The Korteweg–de Vries equation on the quarter plane with asymptotically t-periodic data via the Fokas method,” Asymptot. Anal., vol. 107, pp. 115–133, 2018, https://doi.org/10.3233/asy-171452.Search in Google Scholar
[34] F. R. J. de Barros, M. J. Colbrook, and A. S. Fokas, “A hybrid analytical-numerical method for solving advection–dispersion problems on a half -line,” Int. J. Heat Mass Tran., vol. 139, pp. 482–491, 2019, https://doi.org/10.1016/j.ijheatmasstransfer.2019.05.018.Search in Google Scholar
[35] M. Toda, “Theory of a nonlinear lattice,” in Springer Series in Solid-State Science, vol. 20, Berlin, Springer-Verlag, 1989.10.1007/978-3-642-83219-2Search in Google Scholar
[36] A. S. Fokas, A. A. Himonas, and D. Mantzavinos, “The nonlinear Schrödinger equation on the half-line,” Trans. Amer. Math. Soc., vol. 369, pp. 681–709, 2017, https://doi.org/10.1090/tran/6734.Search in Google Scholar
[37] A. A. Himonas, D. Mantzavinos, and F. Yan, “The Korteweg–de Vries equation on an interval,” J. Math. Phys., vol. 60, p. 051507, 2019, https://doi.org/10.1063/1.5080366.Search in Google Scholar
[38] K. Kalimeris and T. Özsari, “An elementary proof of the lack of null controllability for the heat equation on the half line,” Appl. Math. Lett., vol. 104, p. 106241, 2020, https://doi.org/10.1016/j.aml.2020.106241.Search in Google Scholar
[39] D. Mantzavinos and A. S. Fokas, “The unified method for the heat equation: I. Non-separable boundary conditions and non-local constraints in one dimension,” Eur. J. Appl. Math., vol. 24, pp. 857–886, 2013, https://doi.org/10.1017/s0956792513000223.Search in Google Scholar
[40] A. S. Fokas and B Pelloni, “Generalized Dirichlet to Neumann map for moving initial-boundary value problems,” J. Math. Phys., vol. 48, p. 013502, 2007, https://doi.org/10.1063/1.2405405.Search in Google Scholar
[41] B. Pelloni and D. A. Smith, “Nonlocal and multipoint boundary value problems for linear evolution equations,” Stud. Appl. Math., vol. 141, pp. 46–88, 2018, https://doi.org/10.1111/sapm.12212.Search in Google Scholar
[42] B. Deconinck, B. Pelloni, and N. E. Sheils, “Non-steady-state heat conduction in composite walls,” Proc. Math. Phys. Eng. Sci., vol. 470, p. 20130605, 2014, https://doi.org/10.1098/rspa.2013.0605.Search in Google Scholar
[43] N. E. Sheils and B. Deconinck, “Heat conduction on the ring: interface problems with periodic boundary conditions,” App. Math. Lett., vol. 37, pp. 107–111, 2014, https://doi.org/10.1016/j.aml.2014.06.006.Search in Google Scholar
[44] M. Massabó, R. Cianci, and O. Paladino, “An analytical solution of the advection dispersion equation in a bounded domain and its application to laboratory experiments,” J. Appl. Math., vol. 2011, p. 493014, 2011, https://doi.org/10.1155/2011/493014.Search in Google Scholar
© 2020 Walter de Gruyter GmbH, Berlin/Boston