Abstract and Applied Analysis

Volume 2013 (2013), Article ID 474608, 7 pages

http://dx.doi.org/10.1155/2013/474608

Research Article

Schwarz Waveform Relaxation for Heat Equations with Nonlinear Dynamical Boundary Conditions

Shu-Lin Wu

School of Science, Sichuan University of Science and Engineering, Zigong, Sichuan 643000, China

Received 4 July 2013; Revised 17 November 2013; Accepted 1 December 2013

Academic Editor: Changbum Chun

Copyright © 2013 Shu-Lin Wu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We are interested in solving heat equations with nonlinear dynamical boundary conditions by using domain decomposition methods. In the classical framework, one first discretizes the time direction and then solves a sequence of state steady problems by the domain decomposition method. In this paper, we consider the heat equations at spacetime continuous level and study a Schwarz waveform relaxation algorithm for parallel computation purpose. We prove the linear convergence of the algorithm on long time intervals and show how the convergence rate depends on the size of overlap and the nonlinearity of the nonlinear boundary functions. Numerical experiments are presented to verify our theoretical conclusions.

1. Introduction

Schwarz waveform relaxation (SWR) algorithm can be regarded as a combination of the Schwarz domain decomposition (DD) method [1, 2] and the waveform relaxation method [3, 4]. It is characterized by first partitioning the spatial domain into several subdomains and then solving subproblems simultaneously inside each subdomain through iterations (see [5–7] for the original idea of this approach). In the classical DD framework, one first discretizes the interested problem in time direction by using an implicit time integrator and then solves the state steady problems by using the DD iterative process. Therefore, compared to the classical DD approach, the SWR algorithm possesses the advantages that to be able to treat different subdomains numerically differently with an adapted procedure for each subdomain.

The SWR algorithm has been studied widely and deeply by many authors during the past decade. For example, Gander and his colleagues have investigated the convergence of the algorithm for linear heat equations [5, 8], for the linear (advection) reaction diffusion equation [9–11]. The nonlinear problems have also been investigated; see [6] for the reaction diffusion equation with a nonlinear reaction term and [12] for the nonlinear conservation laws. Nowadays, the SWR algorithm has been becoming more and more popular in scientific and engineering computing and has been adopted to solve complex problems arising from physics and engineering. For example, Martin [13] applied the algorithm to the shallow-water problem; Dolean et al. [14] and El Bouajaji et al. [15] investigated the application of the algorithm for the Maxwell’s equations and Zhang and Jiang [16] applied the algorithm to time-periodic parabolic problems.

However, all of the aforementioned problems are initial value problems or initial boundary value problem with Dirichlet boundary conditions. In this paper, we consider the heat equations with nonlinear dynamical boundary conditions and according to our best knowledge there are no results about the SWR algorithm for this type of problems. This type of boundary conditions arises for many known equations of mathematical physics. They are motivated by problems in diffusion phenomena [17], reaction diffusion systems in phase-transition phenomena [18], models in climatology [19], and many others. For physical interpretations of dynamical boundary conditions for problem (1), the interested reader can refer to [20].

Here, we focus on the convergence properties of the algorithm on sufficient long time intervals and the dependence of the convergence rate on the overlap size and the nonlinearity of the boundary nonlinear functions. The remainder of this paper is organized as follows. In Section 2, we introduce the problem studied in this paper and the SWR iterative formula for this problem. Our main results are presented in Section 3, where we prove convergence of the SWR algorithm and investigate how the convergence rate depends on the overlap size and the nonlinearity of the boundary nonlinear functions. Section 4 provides some numerical experiments to validate our theoretical finding. Finally, in Section 5 we finish this paper by some conclusion remarks.

2. Problem Description and the SWR Algorithm

We consider the following problem on space domain :

where , are smooth functions and satisfying growth assumption

With this assumption, we know from Corollary  3.6 in [21] that the solution of (1) is globally defined and bounded. We decompose into two overlapping subdomains , with and and . Following the spirit of [5, 6], the two-subdomain SWR algorithm for (1) can be written as

where is the iteration index and are initial approximations to the solution of (1).

The comparison principle is central for the convergence analysis of the overlapping Schwarz waveform relaxation algorithm. For the Dirichlet type boundary conditions, this principle is known for long and can be found in every textbook concerning PDEs (see, e.g., [22]) and has been used widely in the SWR community [5, 6, 8, 12]. The comparison principle for dynamical boundary conditions is however much less known, until recently by von Below and de Coster [23]. Rault [24, 25], von Below and Pincet Mailly [26], and Pincet Mailly [27] used this comparison principle to study the blow-up phenomenon for strong nonlinear reaction diffusion equations.

Lemma 1. Let , , and the function satisfy the following differential inequalities

Then, it holds that , for all and .

Proof. See Theorem 2.2 in [23].

Note that, Lemma 1 can be also used to hybrid boundary conditions. For example, the left boundary condition in (4) can be equivalently rewritten as and thus by letting we can also conclude that from (4) by imposing instead of the dynamical boundary condition.

3. Convergence Analysis

In this section, we analyze the convergence of the SWR iteration (3). To this end, we define with , . Subtracting (3) from (1) gives

By using the Taylors’ expansion we get

where are functions of and which lie between and and are defined by

For any function (where ), we define

Theorem 2. Assume that and that the nonlinear boundary functions satisfy

where are constants. Then the error functions of the SWR iteration (3) decay linearly to zero. Specifically, it holds that

where is the called the convergence factor and is defined by

Proof. Consider the following differential equations:

It is easy to find that the solutions to (13) are the steady state solutions

Since , it is easy to know that for all and for all .

Now, let and then by subtracting (6) from (13) we get

By using and , we have and . Therefore, it follows by applying Lemma 1 that ; that is, (). Similarly, the sum are also nonnegative functions and thus the modulus of can be bounded by

and this implies that

From (17) we get

Then, by a simple recursion it can be shown that

where the argument is defined by (11). Now, substituting (19) into (17) and taking supremum with respect to in each inequality gives (10). For the vanishing of the errors , it just needs to notice because and .

Remark 3. Let be the overlap size. Then, by substituting into (11) we have

It is easy to get

Hence, increasing (i.e., increasing the overlap size) will improve the convergence rate of the SWR algorithm. Moreover, the other two derivatives given in (21) imply that the SWR algorithm (3) is especially suitable for problem (1) that the nonlinear boundary functions have very large negative partial derivatives. The limiting case yields

which is the convergence factor of the two-subdomain SWR algorithm for the heat equation (1) with Dirichlet boundary condition (see [5]).

4. Numerical Results

In this section, we present some numerical results to validate our theoretical predictions obtained in Section 3. For (1), we impose initial condition . The source function is chosen as

where . The nonlinear boundary functions are chosen as

With the above choice, it is easy to know that the exact solution to (1) is . We plot the exact solution in Figure 2(a), where we can see that the exact solution is almost zero except around “” due to the power “.” The reason why we choose this function as the reference solution is that it is easy to intuitively see the evolution of the iterates on subdomain during iterations, as shown in Figure 1. We equally decompose the space domain into two subdomains with overlap size , where is an integer and is the space mesh size. Then, for each subproblem we discretize the governing equation by using a central finite difference in space with mesh parameter and a backward Euler method in time with time step . The space partial derivatives appearing in boundaries are discretized by the upwind finite difference method. All of our experiments are performed in the time interval .

(a)

(b)

Figure 1: Influence of the overlap size on the convergence rate of the SWR algorithm. Top: with the solution generated by the SWR algorithm after 2, 4, and 6 iterations; bottom: the similar information for .

(a)

(b)

Figure 2: (a) The profile of on domain . (b) Convergence rate together with the theoretical linear bound of the SWR algorithm with two different overlap sizes.

In this first set of experiments, we investigate the influence of the overlap size on the convergence rate of the SWR algorithm. We choose in (24) and then we show in Figure 1 the solution generated by the SWR algorithm after 2, 4, and 6 iterations with two different overlap sizes (top) and (bottom). Here, we just show the solution profile on the space-time domain for better viewing purpose. The converged solution is plotted in Figure 2 on the left. We can see clearly that the overlap size has a remarkable influence on the convergence rate and in particular with a larger overlap size the algorithm converges faster. This observation confirms our theoretical predication given in Section 3. The measured error (between the converged solution and the iterates, in the sense of norm) together with the linear bound predicted by Theorem 2 is shown in Figure 2 on the right. We see that the theoretical linear bound predicts the actual convergence rate very well.

Next, we study how the parameter in (24) affects the convergence rate of the SWR algorithm. Note that, with the nonlinear boundary functions (24), the arguments and appearing in the convergence factor are and , respectively. We choose three different and plot the measured error in Figure 3. We see that decreasing has a positive effect on the convergence rate, which confirms our theoretical analysis done in Remark 3.

(a)

(b)

Figure 3: Influence of the problem parameter on the convergence rate of the SWR algorithm. (a) ; (b) .

5. Conclusion

We have briefly analyzed the convergence properties of the overlapping SWR algorithm for heat equation with nonlinear dynamical boundary conditions, based on the recently established comparison principle for PDEs with dynamical boundary conditions [23]. The linear convergence on sufficient long time intervals is proved and it is shown that the convergence rate ameliorates as one increases the overlap size. The dependence of the convergence rate on the nonlinearity of the nonlinear boundary functions is also revealed.

Further work includes generalizing the current work to many subdomains for practical parallel computation purpose. Also, to speed up the convergence rate dynamical transmission conditions, such as the popular Robin transmission conditions, should be utilized. Note that the transmission conditions used in this paper are Dirichlet type. For time dependent PDEs (initial problems on unbounded space domain), it is well known that the dynamical transmission conditions always yield significant convergence rate for the SWR algorithm (see, e.g., [9, 10]). Due to the nonlinearity arising from the boundaries and the boundedness of the space domain, it is impossible to straightforwardly generalize the techniques proposed for Dirichlet boundary conditions to dynamical boundary conditions and new notion is eagerly expected.

Acknowledgments

The authors are very grateful to the anonymous referees for the careful reading of a preliminary version of the paper and their valuable suggestions and comments, which really improve the quality of this paper. This work was supported by the NSF of Sichuan University of Science and Engineering (2011RZY04 and 2012RC24), the project of the Key Laboratory of Cambridge and Non-Destructive Inspection of Sichuan Institutes of Higher Education (2013QZY01), and the NSF of China (11301362, 11226312, 11371157, and 91130003).

References

T. P. A. Mathew, Domain Decomposition Methods for the Numerical Solution of Partial Differential Equations, vol. 61, Springer, Berlin, Germany, 2008.

A. Toselli and O. Widlund, Domain decomposition methods—algorithms and theory, vol. 34, Springer, Berlin, Germany, 2005.

E. Lelarasmee, A. E. Ruehli, and A. L. Sangiovanni-Vincentelli, “Waveform relaxation method for time-domain analysis of large scale integrated circuits,” IEEE Transactions on Computer—Aided Design of Integrated Circuits and Systems, vol. 1, no. 3, pp. 131–145, 1983.

Y.-L. Jiang, “On time-domain simulation of lossless transmission lines with nonlinear terminations,” SIAM Journal on Numerical Analysis, vol. 42, no. 3, pp. 1018–1031, 2004.

M. J. Gander and A. M. Stuart, “Space-time continuous analysis of waveform relaxation for the heat equation,” SIAM Journal on Scientific Computing, vol. 19, no. 6, pp. 2014–2031, 1998.

M. J. Gander, “A waveform relaxation algorithm with overlapping splitting for reaction diffusion equations,” Numerical Linear Algebra with Applications, vol. 6, no. 2, pp. 125–145, 1999.

E. Giladi and H. B. Keller, “Space-time domain decomposition for parabolic problems,” Numerische Mathematik, vol. 93, no. 2, pp. 279–313, 2002.

M. J. Gander and H. Zhao, “Overlapping Schwarz waveform relaxation for the heat equation in dimensions,” BIT Numerical Mathematics, vol. 42, no. 4, pp. 779–795, 2002.

D. Bennequin, M. J. Gander, and L. Halpern, “A homographic best approximation problem with application to optimized Schwarz waveform relaxation,” Mathematics of Computation, vol. 78, no. 265, pp. 185–223, 2009.

M. J. Gander and L. Halpern, “Optimized Schwarz waveform relaxation methods for advection reaction diffusion problems,” SIAM Journal on Numerical Analysis, vol. 45, no. 2, pp. 666–697, 2007.

V. Martin, “An optimized Schwarz waveform relaxation method for the unsteady convection diffusion equation in two dimensions,” Applied Numerical Mathematics, vol. 52, no. 4, pp. 401–428, 2005.

M. J. Gander and C. Rohde, “Overlapping Schwarz waveform relaxation for convection-dominated nonlinear conservation laws,” SIAM Journal on Scientific Computing, vol. 27, no. 2, pp. 415–439, 2005.

V. Martin, “Schwarz waveform relaxation algorithms for the linear viscous equatorial shallow water equations,” SIAM Journal on Scientific Computing, vol. 31, no. 5, pp. 3595–3625, 2009.

V. Dolean, M. J. Gander, and L. Gerardo-Giorda, “Optimized Schwarz methods for Maxwell's equations,” SIAM Journal on Scientific Computing, vol. 31, no. 3, pp. 2193–2213, 2009.

M. El Bouajaji, V. Dolean, M. J. Gander, and S. Lanteri, “Optimized Schwarz methods for the time-harmonic Maxwell equations with damping,” SIAM Journal on Scientific Computing, vol. 34, no. 4, pp. A2048–A2071, 2012.

H. Zhang and Y. L. Jiang, “Schwarz waveform relaxation methods of parabolic time-periodic problems,” Scientia Sinica Mathematica, vol. 5, pp. 497–516, 2010.

M. Warma, “Quasilinear parabolic equations with nonlinear Wentzell-Robin type boundary conditions,” Journal of Mathematical Analysis and Applications, vol. 336, no. 2, pp. 1132–1148, 2007.

C. Cavaterra, C. G. Gal, M. Grasselli, and A. Miranville, “Phase-field systems with nonlinear coupling and dynamic boundary conditions,” Nonlinear Analysis A: Theory, Methods & Applications, vol. 72, no. 5, pp. 2375–2399, 2010.

R. G. Watts and M. Morantine, “Rapid climatic change and the deep ocean,” Climatic Change, vol. 16, no. 1, pp. 83–97, 1990.

G. R. Goldstein, “Derivation and physical interpretation of general boundary conditions,” Advances in Differential Equations, vol. 11, no. 4, pp. 457–480, 2006.

A. Rodríguez-Bernal and A. Tajdine, “Nonlinear balance for reaction-diffusion equations under nonlinear boundary conditions: dissipativity and blow-up,” Journal of Differential Equations, vol. 169, no. 2, pp. 332–372, 2001.

C. Pao, Nonlinear Parabolic and Elliptic Equations, Plenum, New York, NY, USA, 1992.

J. von Below and C. de Coster, “A qualitative theory for parabolic problems under dynamical boundary conditions,” Journal of Inequalities and Applications, vol. 5, no. 5, pp. 467–486, 2000.

J.-F. Rault, “The Fujita phenomenon in exterior domains under the Robin boundary conditions,” Comptes Rendus Mathématique, vol. 349, no. 19-20, pp. 1059–1061, 2011.

J.-F. Rault, “The Fujita phenomenon in exterior domains under dynamical boundary conditions,” Asymptotic Analysis, vol. 66, no. 1, pp. 1–8, 2010.

J. von Below and G. Pincet Mailly, “Blow up for reaction diffusion equations under dynamical boundary conditions,” Communications in Partial Differential Equations, vol. 28, no. 1-2, pp. 223–247, 2003.

G. Pincet Mailly, “Blow up for nonlinear parabolic equations with time degeneracy under dynamical boundary conditions,” Nonlinear Analysis A: Theory, Methods & Applications, vol. 67, no. 3, pp. 657–667, 2007.