.nvuphnuu ‘ ‘ ‘ V ' V‘ r ' i I . A . V . > . . V. . 77 r ’ > V r r $145er GAN TATE IVERSITY l l Itillllllll‘llllllll‘l 13mm 3 1293 01394 This is to certify that the dissertation entitled The Discretized Korteweg-de Vries Equation presented by Michael J. Nixon has been accepted towards fulfillment of the requirements for Ph.D. degree infifllhfilflflflii Major professor Datefiflagtdf /2711?97 MS U i: an Affirmative Action/Equal Opportunity Institution 0- 12771 LIBRARY Michigan State University PLACE fl RETURN BOX to romovo this checkout from your rooord. TO AVOID FINES rotum on or baton dot. duo. DATE DUE DATE DUE DATE DUE MSU In An Affirmation Action/Equal Opportunity inotiMion Wm: The Discretized Korteweg—de Vries Equation By Michael J. Nixon A Dissertation Submitted to Michigan State University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Department of Mathematics 1997 i‘ll' ’ it'll! l ABSTRACT The Discretized Korteweg—de Vries Equation By Michael J. Nixon In this dissertation, we study a discretized version of the (generalized) Korteweg— de Vries equation, Btu + Bin + u4axu = 0. After a number of estimates, we utilize the Contraction Mapping Principle to prove the global well-posedness of this equation in a certain discrete Banach space. Our results are analogous to those of Kenig, Ponce, and Vega in the continuous setting. However, due to the nature of the Fourier multipliers, the proofs of several of these estimates in the discrete setting require new techniques. Our results yield a numerical procedure for computing the solution. We present a numerical algorithm which is based on successive iterations to obtain a fixed point guaranteed by the Contraction Mapping Principle. This fixed point is the desired solution which we demonstrate with several numerical experiments. To my family iii ACKNOWLEDGMENTS Several individuals deserve my gratitude. First and foremost, I wish to express my deep appreciation to my advisor, Professor Michael Frazier, for his excellent guidance, encouragement, and commitment to all aspects of my graduate studies. For his efforts, I am extremely indebted. I would like to thank the other members of my Thesis Committee: Professors Patricia Lamm, Joel Shapiro, William Sledd, and Zhengfang Zhou. Each of them has contributed significantly to my graduate experience. Furthermore, I would like to thank Professor Bruce Sagan for bringing a very helpful combinatorial identity to our attention. I would also like to convey my sincere appreciation to my fellow graduate student, Richard Hensh, for his effort and thoughtful suggestions in regards to the numerical examples in my thesis. iv TABLE OF CONTENTS INTRODUCTION 1 The Main Result 1.1 Preliminaries ............................... 1.2 Discretization of (K dV)4 ......................... 1.3 A Discretized Version of Duhamel’s Principle .............. 1.4 The Operator (Pm ............................. 1.5 The Space Xh ............................... 1.6 Contraction Mapping Theorem ..................... 1.7 Well-Posedness of (K dV): in X h ..................... 2 The Estimates 2.1 A Discrete Version of Fractional Integration .............. 2.2 A Lemma Involving Fractional Differentiation ............. 2.3 Proof of Lemma 1.17 ........................... 2.4 Discrete Complex Interpolation ..................... 3 Numerical Results 3.1 The Cutoff Function 7,,(n, m) ...................... 3.2 The Coefficients Q[n, m] ......................... 3.3 A Few Examples ............................. BIBLIOGRAPHY 14 14 18 2O 21 22 25 34 36 36 38 47 83 85 85 90 92 99 Introduction Let 8' (R) be the set of tempered distributions on the real line. Let [3 denote a Banach space such that B C 8' (R) The following initial value problems which arise in the study of water waves are called the (generalized) K orteweg-de Vries equations, 8tu+82u+ukaxu=0 xER,tER,k€N (KdV)k U|t=0 2 U0 E B. We say that (K dV),c is globally well-posed in B if there exists a unique u = u(:c, t) such that the following holds: 0 u E C((—oo, +00); 8), i.e., u(-,t) E B for all t E R and t —+ u(-,t) is continuous from R to B. o u(-,0) = no. (-, t + h) — u(-, t) o For t E R, Ilia?) it exists as an element of 8’ (R). Define —) amp, t) = ’12:?) u(., t + h}: — u(., t). o The nonlinear term uk(-, t)8$u(o, t) E 8' (R) for all t E R. o u is a solution to (K dV),c in the distributional sense. 0 The mapping no —> u is continuous from B to C((—oo, +00); 8). We say that (K dV),c is locally well-posed if the above holds with C((—oo, +00); 8) replaced by C((Ao, A1); 8) for some A0 < O and A1 > 0. 1 For F E 8' (R), let F be the Fourier Transform of F, defined by - 1 . mg) = Zr' A F(:c)e"x‘€da:. For 3 E R, the Sobolev Space of order s, denoted H ”(R), is the set of tempered distributions F for which F is a function and llFllip = l. |F(€)I2(1 + was < 00. The homogeneous Sobolev Space of order s, denoted H‘(R), is the set of tempered distributions F for which F is a function and I|F||2-. = [R |F(€)|2€2‘d{ < 00. Obviously, H ’ (R) C H 3 (R) for s 2 0. It is well known that H ‘ (R) and H 3 (R) are Hilbert spaces with inner product in H ‘(R) defined by < HG >= AF(£)5(£)(1+€2)Sd£ and in H 3 (R) by < F,G >= [R F(€)5(£)§2’d€- In their 1976 paper [1], using energy methods, J. Bona and R. Scott proved that (K dV),c is locally well-posed in H ’(R) for s > 3/ 2, for all k E N. By the Sobolev Imbedding Theorem for s > 3/ 2, H ‘(R) C Cl(R), the set of all continuously dif- ferentiable functions on R. Hence, this result is restrictive in the sense that it can not guarantee a solution to (K dV),c if the initial data no has, say, a jump discontinu- ity. However, the homogeneity of the nonlinear term suggests that well-posedness in H ‘(R) may hold for smaller values of 3. That is, if u solves (K dV)k, then for A > 0, so does u,\(:c, t) = A2/ku()\:r, A3t) with initial data u,\(a:, 0) = Az/kuO‘m, 0). Since we are using a scaling argument, let us assume that these solutions are in the homogeneous Sobolev space, H"(R). Note that lluA($,0)llip = AQ/kA’A‘l/2IIU($,0)IIH» This suggests that the optimal s for (K dV)k is NIH | Prlto 8);: Note that 3;, Z 0 if and only if k 2 4. In their 1993 paper [5], Kenig, Ponce, and Vega proved that (K dV),c is globally well-posed in H” (R) for k 2 4 provided that the initial data is small, i.e., Iluol I gm 3 6 for some 6 > 0. If the initial data is not small, they proved local well-posedness in H 3" (R) Note that by Plancherel’s Theorem, H°(R) = L2(R), the set of square integrable functions on R. Thus, for k = 4, the smallest value of k for which Kenig, Ponce, and Vega obtained well-posedness in H ”k (R), we have global well-posedness in the familiar function space L2 (R) Because of this, the remainder of this thesis will focus on the Korteweg-de Vries equation with k = 4 and with B = L2(R), i.e., 6tu+82u+u40xu=0 a: E R,t€R “{th . Ult=o = no G L2(R). We now give a brief outline of the proof in [5] for k = 4. Consider the associated homogeneous linear equation wltzO = 100 E L2(R) and the corresponding inhomogeneous equation 6v+0§v= (IH) t .9 'U]t:0 = 0. One can easily verify that (H) is globally well-posed in L2(R) and the solution is given by wait) = W(t)wo($) = Ae‘<$‘+“3>wo<5)d€. (1) Note that W(t) is a Fourier multiplier operator with multiplier ems, i.e., (W(t)wo)“(€) = e“‘3<120(€)- Using Duhamel’s Principle, one can obtain that (I H ) is globally well-posed in L2(R) for 9 such that g(-, t) E L2(R) for all t E R with the solution given by v(:r, t) = [0‘ we — t')g(-, t')(a:)dt'. (2) To see the connection with (K dV)4, if Su(:c, t) = W(t)u0(:z:) — [0‘ we — t')(u46$u)(o, t')(:z:)dt', (3) then (1) and (2) imply that u solves (K dV)4 if and only if u = Su, provided Sn makes sense as a distribution. In other words, a is a solution of (K dV)4 if and only if u is a fixed point of the nonlinear Operator S. More precisely, for T E (0, +00], define 2T = {veC((—T,T);L2(R))nL°°((-T,T);L2(R))= ||v||L2Lpo < 00 and ”mung; < 00}, where T We l/p IIvIILm = (f, (LT Iv(m.t)l"dt) dx] forlgp,q§00. Ifq=00, then 1/10 may“ sup |v(:v,t)l"drr) tel—737"] with a similar definition for p = 00. The norm for these Banach spaces is ”UHZT =max{ SUP llv('at)l]L2;llUlngL10;llaxvlngotzl- t _ T] T T 1 By proving various estimates with these mixed norms involving the terms W(t)uo($) (4) and [we — t’> O. For arbitrary initial data, they were able to find a T > 0 for which S is a contraction on B, C ZT. Hence, the Contraction Mapping Principle guarantees a solution to u = So in this ball in ZT. Contraction mappings and fixed point iterations are useful tools in the numerical study of differential equations. Can the method developed in [5] be adapted to show that a discrete version of (K dV)4 is well-posed? Is the solution a fixed point of some operator in some space? If so, does this lead to a method that can be implemented numerically? In this dissertation, we answer these questions in the affirmative. The general outline of the proof for the discrete setting is analogous to that of the continuous setting. However, there are major differences that make the discrete proofs different and in some cases quite a bit more involved than their continuous counterparts. The first step in this process is to discretize (K dV)4 in a natural way. To his end, fix h > O and let us first consider the associated linear equation (H). Let w(n, m) be a discrete function defined on Z x Z. The obvious replacement for 63u(x,t) is aghwm, m) (see Definition 1.6). However, the replacement of am is not so clear. Assuming m > O, on the one hand, we could replace 6tu(:r, t) with w(n, m + 1) — w(n, m) h3 and (H) becomes w(n, m +1) - w(n, m) h3 + 83,,,w(n, m) = 0. (6) Note, for homogeneity purposes we are letting h be the step size in the LII-direction, while h3 is the step size in the t-direction. On the other hand, we could replace Btu(a', t) with w(n, m) — w(n,m — 1) h3 and (H) becomes w(n, m) — w(n, m — 1) h3 + 8,3,,hw(n, m) = O. (7) The scheme in (6) is referred to as an explicit scheme, because, given the values w(n, m) at height m, one can solve explicitly for the values of w(n, m + 1) at height m+ 1. The other is called an implicit scheme since the values w(n, m) at height m are determined implicitly by the values w(n, m — 1) at height m — 1. It may seem that the explicit scheme is the more convenient of the two. However, considering the Fourier multiplier which corresponds to ems, we will see that the implicit scheme is the better choice. Here, the Fourier Transform takes functions defined on Z to functions defined on [-7r/h,7r/h] (see Definition 1.2). Let wo be the initial data. If we discretize explicitly, taking the Fourier Transform of both sides in the first variable, (6) becomes w(1)(g m+ 1)- (”(0, m) isin3(h9)a,<1) h3 + ’13 (0, m)=0, which implies am, m +1): (1 — ai”(0,m) and hence, afflw, m) = (1 — isin3 (ham, <1 (0). Notice the multiplier (1 — isin3(h9))’" has magnitude (1 + sin6(h0))'"/2, which blows up for 0 < l0] < 1r/h as m —) 00. Hence, by Plancherel’s Theorem (Proposition 1.3), the norms of w(-, m) in the discrete space 1,2,(Z) (see Definition 1.1) go to infinity with m. This is undesirable both mathematically and physically. Applying the same argument to (7 ), it follows that will?) (1 + isin3(h6))m aw, m) = (for the proof of this see §1.2). Note that the multiplier (1 + i sin3(h0))"" has magni- tude which is bounded by one for all 0 6 [—7r / h, 7r / h] and hence, the norm of w(, m) in l2 h(Z) will be bounded by ||wo||,2 for every m. Therefore, we choose the implicit scheme 1n (7) and the resulting multiplier I 1 + isin3(h0)' (8) The fact that the multiplier is not a pure exponential, like eitfa, causes difficulties in the proofs of the discrete case that do not arise in the continuous case. For example, consider the following estimate proved in [5], which is used in establishing that S is a contraction mapping. Lemma 0.1 Let mo 6 L2(R). Then llaxW(t)wol|LgoL§ S Cllwolltz, where 1/2 anrm = sup (/ WW) . 16R R Proof: Recall, W(t)wo(:r) = [R ei($£+‘€3)wo(€)d€. It follows, after the change of vari- ables 53 = 7, that azW(t)’LUo($) = i/ée'($€+t€3)co( £)d€ : fife itr eixrl flT—2/3+1/31D0(T1/3)dT. Thus, using Plancherel’s Theorem in the t-variable, - 1/3 _ .. 2 6131 7. 1/3w0(T1/3)l d7" 1 Hawumonig = 5 [R c (1. words again using 5 = 71/3, which finishes the proof. This proof is simplified by the fact that the multiplier eit’i3 has magnitude one for all 5 6 R, and the fact that after the substitution 7 = 63, the term e“T is just what is needed to regard the integral as a Fourier transform in the t-variable. The multiplier (1 + isin3(h0))"" is not exactly of this form, so when we proceed as in the proof of Lemma 0.1 we encounter an “error” term after the change of variables. The proof will be reduced to the following inequality (see Proposition 2.9): W . 2 / e“'""r cos'mlrg(r)dr mEZ ’1' 2 s c f; Igmrczr, (9) where the cosIml r is the anticipated extra term. If the proof was completely analogous to the proof of Lemma 0.1, we would not have this term and (9) would become a statement of Plancherel’s Theorem involving Fourier series. Instead, the proof requires some additional techniques including a combinatorial identity (see the proof of Claim 2.11 following Proposition 2.9) pointed out to us by Bruce Sagan. Next, consider m < 0. The implicit scheme in this case is w(n, m +123- w(n, m) + fi’hwm, m) = 0. (10) 10 After applying the argument in §1.2, (10) becomes cow) (1 — isin3(h0))|ml «rifle, m) = and the resulting multiplier is 1 1 — isin3(h0)° (11) Notice that because of our need to control the size of the multiplier, we now have two multipliers, (8) for m > 0 and (11) for m < 0. Recall that the Fourier multiplier in the continuous case associated with W(t) is eitia. Since - 3 - 3 - 3 et(t+s)£ : eit£ €136 , the family of operators {W(t)}teR forms a group under composition, i.e., W(t + s)wo = W(t)W(s)w0 for all s,t E R. This plays an important role in the proof of the key estimate ‘ t I I I la: / W(t—t)g(-.t> 0 h3 (KW): n(n, 0) = new m z 0 1 _ 1 770% m + h)3 r)(n, m) + 3131,1100”: m) + gan,,,(n5)(n, m) = O m < 0’ where 770 is the discrete initial data. If Hh(m)no(n) is the analog of (4) and solves the discrete version of (H) and 1 gAhan,h(715)(n, m) is the analog of (5) (see Definitions 1.12 and 1.14), then the analog of the Operator S in (3) is mm. m) = Hh(m>no(n> — §Ah6.,h(n5>(n. m). It will follow (see §1.4) that 17 solves (K dV)2 if and only if n = q’non- Analogous to the space Zco from [5], we will define a discrete space X h (see Definition 1.16). The primary conclusion of this thesis is the following: 12 Main Result: Suppose the initial data 770 is small, i.e., nEZ 1/2 llnollzg = (h 2 WWW) S 60 for some do > 0. Then there exists r > 0 such that the operator (1),,0 is a contraction mapping in Br = {77 E Xh 3 ”77]th S T}, with 60 and r independent of h > 0. Consequently, by the Contraction Mapping Prin- ciple, there exists a unique solution 17 E B, of (K dVfl. Moreover, n E l°°(l,2,(Z); Z) and the map 770—”) is continuous from lflZ) to Xh. Hence, (K dV)f,i is well-posed in Xh. This result lends itself to a numerical procedure for computing the solution to (K dV)Z. To do this, we first numerically implement the operator ,,0. This will be done by using the coefficients Q[n, m] which have been precomputed and stored (see §3.2). Then, by picking the initial guess to be the zero function on Z x Z, we run the fixed point iteration. If the initial data is small enough, then the main result guarantees that the iterates (bglm) converge to the solution of (K dV): as n —> 00. Also, because the contraction mapping constant is less than one, the convergence is exponential. The fact that the operator (Pm, and the iteration scheme are easy to implement is the motivation for this thesis. This approach differs from the standard approach in several ways. The standard approach is to solve the difference equations iteratively going up one level at a time, 13 while we work on the entire space Z x Z for each iteration. A disadvantage of our approach is the amount of computer storage required. Each global iteration has to be stored in order to define the next iteration, whereas when solving the difference equations level by level, one can discard the data at each level. However, the advan- tages of our method include the rapid computation of the entire solution on the entire grid, and the fact that we avoid problems with ill-conditioning and accumulation of round-off error. In addition, our approach guarantees certain size estimates on the iterations, namely the three norms of the space Xh (see Definition 1.16) including a bound on the l2 norms of each level. Finally, note that our result holds for all h > 0 not just, for example, small h. This distinguishes our approach from another possible approach in which one tries to compare ,,o, for small h, with S, which is known to be contraction map in the continuous setting. Thus what we have obtained is a true difference equation result. In Chapter 1, we state several definitions which are used throughout this the- sis. After this, we show that our choice of discretization yields a nonlinear discrete operator whose fixed point solves (K dV)f,‘. Then, we state the crucial estimates and define the discrete Banach space X h in which we find our solution. Finally, with these estimates, we establish that we have a contraction mapping, which quickly leads to a solution of (K dV)j,‘. This part of the proof follows the same general outline as in [5]. Chapter 2 contains the proofs of the main estimates. It is here that techniques dif- ferent from those in [5] are required to overcome the difficulties noted above related to the the multipliers in the discrete setting. Finally, results of various numerical experiments are presented in Chapter 3. An explanation of the numerical algorithm is included there as well. CHAPTER 1 The Main Result 1 . 1 Preliminaries Let h > 0. In the following definitions, let a and 9 denote functions defined on Z and w denote a function defined on Z x Z. Definition 1.1 o E 1,2,(Z) if Hallzg = (h 2 |0(n)|2)1/2 < 00. 7162 If h represents the step size and 0(n) = f (nh) for, say, f E C (R) 0L2 (R), then Hallfg represents a Riemann sum of I f |2. Hence, ”Gilli ~ I] f H L2 for small h. Definition 1.2 For 0 E [—7r/h,7r/h], let 6,,(0) = h z a(k)e"‘ho. kEZ 14 15 Note that by Fourier inversion, _ 1 "/h . -inh0 0(n) — 27r [Tr/h ah(0)e d0. Proposition 1. 3 If a E l2 ,,,(Z) then 1 “M 1/2 Halli; = (5; lo. (more) . Proof: 1 "/h . 2 _ 1 "/h ikho tho 2n/_./1."”‘(9)'d9 _ 27: /_W/h(hZa(k)e) ()(hzao )eJ kEZ jEZ = 72/1r/hh /h h2 Z Z( 009 —j)ei(k-j)h0d6 kEZjEZ 1 vr/h = -— h2 k 2d 2.1,. EM 9 kEZ = hEIOUCV kEZ = Halli;- Definition 1.4 For a and g E 1,2,(Z), let (o*g)(n =(lc)h£a kEZ be the convolution of or and 9. Remark 1.5 One can see as in the standard case with h = 1 that nEZ kEZ (0 =1: 9),, 2h E (h 2 o(k)g( n-—k)) einho 16 = h 2: 0(k)eikh0h 2 9(7), _ k)ei(n—k)h0 lcEZ nEZ = 5h(9)§h(9) for 0 E [—7r/h,7r/h]. Definition 1.6 Let 0(n +1) — 0(n —- 1) 6ho(n) 2 2h For k E N, define 6§o(n) inductively, i.e., 8§o(n) = chat-10(7),). o(n+3) —3o(n+1)+30(n— 1) —o(n—3). For example, 03,,0(n) 2 8h3 Proposition 1.7 Let a E l§(Z). Then A —-i sin h0 . (6h0)h(0) = ——h(—)-oh(0). Proof: By definition, (0h0)7,(6) = h26h0(k)eikh9 keZ = h: 0(k‘l' 1) ”JUV- 1)eikh0 1Icez 2h 2 _ 220( [9+ 1)€ ikh0__ 5W )eikhe. 2k€Z 2kEZ 17 By setting j = k + 1 in the first sum and setting j = k — 1 in the second sum, we have (aloha) = $2; 0(j)e"'"”(e""‘” — em) singhfl) 5;; (0). =—i By induction on k, it follows that Definition 1.8 For a E lfi(Z), fl 6 C, let sin(h0) l9 h &h(0)e"""9d0. 1 1r/h D5001) : 57; ./;1r/h Remark 1.9 Due to the homogeneity of the associated linear equation, if h is the step size in the x-direction, then h3 will be the step size in the t-direction. This motivates the following definition. Definition 1.10 For 1 g p, q s 00, w E lfilfn,h(Z x Z) if nEZ mez p/q llp (bum:(hzhazlwmmn‘?) ] O h3 (Ken/)2 170210) = now) m z 0 770% m +123- "(na 7’” + aghna, m) + gamrxn, m) = 0 m < 0- We discretize the linear equation implicitly as well, w(n, m) -— w(n, m — 1) + 83,hw(n,m) = 0 m > 0 h3 w(n,0) = 01001) m = O (1.1) 1 _ W(n,m 'l" ’33 W(nam) + ag,hw(n’ m) = 0 m < 0. Assume m > 0 and let w(',m) 6 lf,(Z) for all m 6 Z be a solution to (1.1). By taking ‘ of both sides of (1.1) in the first variable, we have affix, m) - affix, m — 1) isin3(h6) h3 + h3 (12,9)(9, m) = 0, 19 with (22,9)(6, m) = h E w(n, m)ei"h9. 7162 This implies wile. m — 1) 1 + isin3(h6) ’ cow, m) = and hence, .. (1) _ (L00)h(0) to), (01m) ‘ (1 + isin3(h0))'"' If m < 0, reasoning as above, it follows that .(1) __ (030)}:(9) wh (0’ m) — (1 —— isin3(h0))lml' This leads to the following definition. Definition 1.12 For m,n 6 Z and wo E 1,2,(Z), let Hh(m)wo(n) 1 [m e“""” l,2,(Z) is bounded with operator norm less than or equal to 1. Also, Hh(0) = I where I is the identity operator on l,2,(Z). Moreover, by the argument above and Fourier inversion, if w(n, m) = Hh(m)wo(n), then (12 is the unique solution to (1.1) for w(-, m) 6 1?,(Z) for all m E Z. 20 1.3 A Discretized Version of Duhamel’s Principle Proceeding along the same path as [5], we now consider a discretized version of the inhomogeneous linear equation. Definition 1.14 For (1) such that w(-,m) E lflZ) for all m E Z, let h3ZHh(m +1 —j)w(',j)(n) m > 0 j=l Ahw(n, m) = 0 m = O -—h3 i: Hh(m —1-j)w(-,j)(n) m < 0. Proposition 1.15 Let n be such that n(-, m) E 12 ,,(Z) for all m E Z. Let w(n, m)— — Ahn(n, m). Then w solves w(n, m) — w(n, m — 1) + aghwm, m) = n(n, m) m > O h3 wmm+1-wmm ( ’33 ( ) + 33,,hw(n, m) = 17(n, m) m < 0. Moreover, w is unique for w(-, m) E lf,(Z) for all m E Z. Proof: The case m = 0 is trivial by definition of Am. Assume m > 0. Then (“’"ml 3:,“ m " 1) + 63,hw(-,m))i(6> = (ZHh(m+1—j)n(,j)- -EH Hh(m— ])77(' J) +h3§233,hHh( m(+1-J')n ( 9)) (9) j=1 ___ 2m: 1 fl(1)(j)0 _mE—f 1 fi(1)(0 ) j=l( (1+isin3(h0))m+1‘3 h ’ i=1 (1+isin3(h0))m'j h "7 +isin3(h9) Z offlw, j) j=1((1+isin3(h0))'"+1 J 21 _ m 1(1) — Z (1+isin3 (h6))m+1—Jnh (0’ j) j=l( 1(1) (1+isin3 (h6))"‘+1 17711 (6 j) —(1+isin3(h(6))))zl1 m 1 1(1) . +ZSin3( ()h6 )1; (1+isin 3(h6))m+1‘jnh (9’3) «(1) _ 77h (6,m) 3 6(1) — 1+ isin3(h6) —isin (h6) 12:1(1 + isin 3(h6))m+1‘Jnh (0 j) m 1 ~(1) +isin3 ()h6 )E=:1( 1+isin3 (h6))m+1‘ _J-77h (91]) H. 1+ isin 3(h6) ,(1) " 1+isin3(h6) " (9’7") = 17%, m). Thus, (whm) 3:0,?” -1) + 63,hw(',m)) (9) — 6i1)(91m)1 with n£1)(6,m) = h 2 n(k,m)ei"h9. The proposition follows after integrating both kEZ sides against em”. The proof for m < 0 is very similar. 1.4 The Operator (1)170 Our goal is the following: given no E lf,(Z), find 17 such that n(-,m) E lf,(Z) for all m E Z and 776». m) = Hh(m)no(n) - éAhan,h(775)(n1m) ”(n1 0) = 770 (n) (1.2) for m, n E Z. Combining Remark 1.13 and Proposition 1.15, such an n solves (K dVfl. To find an n solving (1.2), we will prove that the operator nun, m) = Hh(m)no(n> — §A16.,1(n5>(n,m> 22 is a contraction mapping on a ball centered at the origin in a particular Banach space. 1.5 The Space X), We now define the discrete Banach space upon which we will show (Pm, is a contraction. Definition 1.16 Let X), = {W(mm) E l°°(l?.(Z); Z) 3 llwllzgzgg‘h < 001||3n,hwllzgozghh < 00}, with ilwiixh = max{81ég ||w(-, m)llz§; llwlll5llo illan,h“1llt°<>z2 } m n m.h n m.h Obviously, X h is a Banach space with the above norm. Using the estimates below, we can find the solution to (1.2) in this Banach space X h. This will be done by proving that (Pm, is a contraction map from a small ball around the origin in the space X), to itself. The estimates needed are the following: Lemma 1.17 For h > 0, 2') ||3n,hHh(m)Tlollt:°l3n,h S cllnollt: and ii) ||63,,,,Ahwlllgozgn.h S Cllwllzgtfmi with c independent of h > 0. 23 Lemma 1.18 For h > 0, 2') lth(m)770llzgz;gh s cunt, and 12') ”mum, s cunning), with c independent of h > 0. Before proceeding, we need the following proposition. Pr0position 1.19 If (.1) E Xh, then 6",),(105) E 12/41253, with ilafl,h(w5)iilz/4l;fl/h9 S 5||wllih Proof: By definition, 5 _ 5 _ anh(w5)=w(n+1,m) w(n 1,m) ’ 2h = w(n + 1,m)2—hw(n —1,m) (w4(n +1,m) + w3(n +1,m)w(n — 1,m) +w2(n + 1, m)w2(n — 1,m) + w(n +1,m)w3(n — 1,m) + w4(n — 1,m)) = 8n,hw(n, m) (w4(n +1,m)+ w3(n +1,m)w(n —1,m) +w2(n +1,m)w2(n --1,m)+ w(n +1,m)w3(n - 1,m) + w4(n - 1,m)). Set V(n, m) = w(n + 1, m) and u(n, m) = w(n — 1, m). It follows that 4 ]]6n,h(w5)]]l:/4110/9 S Z“(amhwlu4—jl‘jii15/4110/9 m,h j=0 n 771,): 4 9/8 4/5 2: h: (h3 Z Ian,hw|10/9|V4—j#j|10/9) i=0 nEZ mE Z 24 24: (h 2 (ha 2 lan’hwlz) 5/8 (h3 E lu4‘jw'l5/2)1/2)4/5 S j=0 nEZ mEZ mEZ 1/2 4 1/2 4/5 3 Sup h32|0n,hw|2 2 fix h32:|1/4"1uj|"’/2 "Ez mEZ j=0 nEZ mEZ 4 S iiwllthAj’ :0 with nEZ mEZ 1/2 4/5 Aj : (h 2 (hit 2 lV4—j/J‘ji5/2) ) . Note that the second inequality above follows from an application of Hblder’s Inequal- . . . 9 1ty on the sum 1n m w1th p = -5- and q = 4' Claim 1.20 For V,u E X), and 0 g j S 4, A: s llVllx,’llull’x,,- Proof:(of Claim 1.20) Fix j. Set p = 4 4 4 and q = 3 (p or q = 00 when dividing by zero). Then using Hblder’s Inequality on each sum, it follows that 4—3-1 t- 4/5 Aj S It: h3 Z lVllo ha 2: l/Jllo nEZ mEZ mEZ . (h;(hs;a.ullo)l/2)*(h;(ha;a.1.m)‘”)i 4—“ 2' llVllthlll‘llxh- |/\ /// 25 Since Ill/Hxh = ||u|lxh = lelxh, it follows that A,- g llwll‘fiq. Therefore, “albino/m9 3 sum, which concludes the proof of Proposition 1.19. 1.6 Contraction Mapping Theorem In this section, assuming Lemmas 1.17 and 1.18, we prove (Pm, is a contraction map- ping on a ball centered about the origin in the Banach Space X h. Theorem 1.21 There exists do > 0 and r > 0 such that, if B, = {w E X), : Iiwilxh S 7'} and ”Halli; _<. 50; then i) (Pm, : B, ——> B, continuously and ii) ||,,0(1/) — ,,o(u)llxh 5 HIV — ullxh, for some )1 < 1 and V,u E B,, with A, do, and r independent of h > 0. Proof(i): Recall, the operator (Pm, has two terms, Hh(m)no(n) and Ahan,h(n5)(n, m). In order to establish the existence of an r > 0 such that (1),", : B, —+ B,, we need six estimates, three involving each term. First, we focus on Hh(m)no(n). Since the multiplier associated with the operator Hh(m) is bounded by one, we have SUP lth(m)770llz§ S ”Wilt; S 50- (1-3) mEZ 26 By Lemma 1.17(i), llan,hHh(m)7iollz;-,ozfmh S Cllflollzg S C<50 (1-4) and by Lemma 1.18(i), lth(m)770|lzgz},9,h S Cllnolltg S 050, (1-5) with 0 independent of h > 0. Combining (1.3), (1.4), and (1.5), we can conclude Hh(m)770 E X}; and lth(m)770||X1. S 050, (1-6) with c independent of h > 0. Next, consider the nonlinear term Ahan,h(n5)(n, m). Fix m > 0. For a sequence 770, let SinnoWaJ') = Hh(.l — m — 1)3n,h770(n)Xlgjgm- Again, by Lemma 1.17(i), llshllouzgozffi S Cllnollzgi with c independent of m > 0. By duality, Sm : 1,1,13,,,(Z x Z) -> l§(Z) is bounded with bounds independent of m > 0 where Smw(n) = h3 i Hh(m +1 —j)(9n,hw(-,j)(n) J—l = Ahan,hw(n, m). Therefore, llAhaann m)|lz§ S Cllwllzizihi (1-7) 27 with c independent of m > 0 and h > 0. Letting w = n5 in (1.7) and taking the supremum over m > 0, it follows that sup iiAhan,h(775)('il’fnlillfI S Cllnsllzgzzh m>0 J. = Cllnlliw, The same proof can be used for m < 0 (by definition, Ahan,h(n5)(n, 0) = 0). Hence, s13; ”Ataturxummo s Cllnll‘},- (1.8) Next, by Lemma 1.17(ii), llan,hAhan,h(775)lil;',°lg1,h '3 lla§,hAh(775)llla°lfn,h |/\ curling", = Cllfllli’gtfiffi Cllnllio- (1-9) |/\ Finally, by Lemma 1.18(ii) and PrOposition 1.19, l/\ llAhang (775) l hang, cilan,h(775)||1:/4z’1:/h9 < 5c]|n||§(h. (1.10) Thus, ifn E Xh, combining (1.8), (1.9), and (1.10), it follows Ahan,h(n5) E X), as well. Furthermore, llAth,h(n5)llx,. S Cllnllio, (1-11) 28 with c independent of h > 0. To conclude the proof of (i), if n E X), with ”77]th S r, then by (1.6) and (1.11) there exists c1 and e; such that “(1)1707“th < 0160 + 627‘s. If (So and r are chosen small enough, then H‘I’noflllxh S 6150 + C275 S 7' which implies (13,70:B, —>B,. Proof(ii): Let u and u be in B, C Xh. To prove (1),,0 is a contraction mapping on B, for some 1' > 0, we need three estimates, one for each norm in the definition of H ' Hxh. First, 1 SUP ”(pilot/('1 m) - q’noManII}, = g SUP llAhan.h(V5 " #5)(‘1’m)llz§ m>0 m>0 1 m . . = g SUP “’13 Z: Hh(m + 1 " J)3n,h(V5 — #5)(',J)llz§ m>0 1'21 S 0|le - milling“, where the last inequality follows from (1.7) with w replaced with v5 — ’15. As before, this can be proven for m g 0 as well. Thus, Sig; i|@00V(°im) — @flo”('1m)lllfi S Cii(V — ”)(V4 + V3” + V2l1'2 + VH3 + ”finial?" h m . 1/2 = ch 2: (h3 2 |(u — [1)(1/4 + V311 + uzuz + 1413 + #4)|2) nEZ mEZ 29 1/10 2/5 S— chZ (ha 2 IV-ul‘“) - (123 Z |u4 +V3u+u2p2+ufl3+fl4|5fl> nEZ mEZ mEZ c (h; (113:; (u — ”(10) W) 1/2 4/5 . (h E (h3 2 lu" + V311. + V2112 + 1413 + “4'5/2) ) 1/5 l/\ nEZ mEZ S CHV — pllxh(Ao + A1 + A2 + A3 + A4), with 1/2 4/5 A) = (’1 2: (ha 2 IV4‘jujl5/2) ) nEZ mEZ for 0 g j g 4. By Claim 1.20, A,- 3 r4. Therefore, 311g; ll‘PnoI/(w m) — ¢,.u(-, mm), s 07‘4”!» — .Ullxh- (1.12) Note, 3”,,(as)(n, m) = 8n,ha(n, m)6(n + 1, m) + oz(n — 1, m)6,,,hfl(n, m). Thus, 1 ”@170” " q’ml‘llzszlo = gllAhan.h(V5 — #Ullzgzgfm 7‘ m,h S Cilan,h(V5 "' ”SHIP/411°” n "uh nEZ mEZ 9/3 C (h 2 (’13 2 i071). [(11 - 160/“ + 1231‘ + V213 + 1413 + #4)] Im/g) ) |/\ c (h 2 (h3 2 I8",h(1/ — u)(n, m) nEZ mEZ 4 3 2 2 3 4 10/9 9/8 4/5 -(v +vu+uu +vu +u)(n+1,m)| ) 30 1162 mEZ +c (h 2 (h3 2“]! — u)(n —1,m) 4 3 2 2 3 4 10/9 9/8 4/5 4%MV+V#+VH+WM+M)MmM = 1+11, where the first inequality follows from Lemma 1.18(ii). By applying Hiilder’s inequal- 9 9 ity to the sum on m with p = E and q = 4’ it follows that 5/8 I g c (h 2 (h3 2 Ian,h(V — u)(n.m)|2) nEZ mEZ 1/2 4/5 5 2 . (I13 2 ](1/4 +V3u+u2u2 +1413 +u4)(n+1,m)] / ) ) mEZ 1/2 3 csup (h3 2 lan,h(V — H)(n:m)l2) nEZ mEZ 1/2 4/5 2 . (h E (h3 E ](1/4 + V311 + V2112 + VH3 + u4)(n, m)]5/ ) ) l/\ CHI/ - [.tllxh(Ao +A1 + A2 + A3 + A4) 3 emu—mu, min with the last inequality following from Claim 1.20. Next, we consider II. By applying H61der’s Inequality to both the sum in m with p = 9 and q = g and to the the sum 4 in n withp=4 and q: 3’ we have 1/2 II g h: (h3 Ely—[1]”) nEZ mEZ 4/3 - (h 2 (h3 2 lamb/4 + V3lt + 112/12 + V113 + 114)]5/4) ) 1/5 3/5 7162 mEZ 31 4 f 4/3 3/5 S IIV-MIIth hZ(h3:I3n.h(I/4_’ WW”) ) j=0 \ nEZ mEZ 3/5 4 I 54 4/3 s Ilv-ulli: (h:(ha:lanturjxmmmum1,m)l’) ] j=0 nEZ mEZ - 4/3 3/5 +(h:(is:(um-1,m)o,,(,o(,,m)]5“) ) - (1.14) 7162 mEZ Claim 1.22 Let V and u be in Xh. Then for an integer j such that O S j S 4, 4/3 3/5 (h: (ha 2]” 4 3&1th ’)5] H) ) SJIIVIIx,jII/1IIx, nEZ mEZ Proof: Fix j such that 1 S j g 4 (the case j = 0 is trivial). Using Héilder’s Inequality I 8 C O. I I on the sum in m w1th p1 = , and q, = —+4’ then usmg Holder’s Inequality agam — .7 on the sum in n with p, = 4_ j and q2= -—_—1-(p,-, q,= 00 for i = 1, 2 when dividing by zero), it follows that (z (W >) nEZ mEZ 4—;1 1— “5 s (1263:2110) -,<(h3:|a.,> 1% 1162 mEZ mEZ 1/2 71 10 51%, 5 s (h: (It: >314”) ) (h: (13: law, )7) ) nEZ mEZ nEZ mEZ 4- S IIVIIX: II mEZ k=0 j—l . 3% 21—2 h 2: (ha 2 (Z Ian,h.u'(na m)||u-7_1-k(n +12mIII/‘k(n —1,m)|) ) nEZ 32 4_ S IIVIIX,J - 1 1— it?“ ‘7_ log—14:2 ’— 2 h: ’33: Ian,h#(n,m)IJ—I+%I#(n+1,m)| 1"“ INN—1 m)I‘Lgi“) k=0 nEZ mEZ : IIVIIX—hjZBjJCa with 17—1 5 4 21—2 _151, log—14:2 _1_0_ BM = (h: (’13 Z lan,hu(n, mMIMI/102+ 1,m)l 1+4 lu(n- 1,m)|1+'i) ) nEZ mEZ ' 4 Fix It. Again, apply Holder’s Inequality to the sum on m with p, = J: and . ,_ _ 1 q1 = j _ 1 and to the sum on n with p2 = 75—1—3 and q2 = 316—. This gives 7.53 1/2 “I?" 10 '—1-—k 1_0k_ Bch S (h: (’13 Z lamb/42) (ha 2 W(n +1,m)| 1" W(n — 1 m)| 1) nEZ mEZ mEZ 1/2 1/2 31 10 '—l—k 10’: .<. sup (M 2 law?) h: (,3 z: |u(n +1,m)I4L—ZJ—l lu(n —1 mu;- ) "EZ mez nEZ mEZ -_1 I; ig—l 20— 56715 S llullx, h 2 (’13 Z lu(n+ 1,m)lm) (h3 E |u(n —1m)|1°) neZ mEZ mEZ 1/2 t—ii 1/2 i S llullx, h: (’13 Z lull") h: (’13 Z lull”) nEZ mEZ nEZ \ mEZ S Willi}, Therefore, 54 IhZIhszI van,“ ')| ) ) ,,017. By definition of (Pm, 17 solves (1.2). Also, by definition of the space X h, 77 E l°°(l2(Z) ; Z). Moreover, we have the following theorem. Theorem 1.23 Let 60 and r be as in Theorem 1.21. Let 170 E lfi(Z) such that ”170”,: g 60 and let 17 6 B, be the unique solution of (K dV)flI guaranteed by Theorem 1.21. Then the map ’70—’77 is continuous from B(0; 60) to B, with B(0;60) = {o : ||0||7§ S 60}. Proof: Fix h > 0. Let am, 770 6 B(0; 60). Then the unique solutions to (K (1V): in B, with initial data we and 170 are w and 17, respectively, which satisfy 1 5 w(n, m) = tho(n) — gAhBhQu )(n, m) and m m) = Hh770(n) — éAhamSXn, m). Thus, w(n, m) — n(n, m) = Hmong — noxn) — §Ahah 0. Let p > 0 which is yet to be determined. Set k1 = k . X{OSInL<.u}, k2 = k - X{lnl>u}’ and let b E l”(Z) with ”bllzp = 1. Then b * k = b * k1+ b * k2, both of which are well defined. This leads to #{n: |(b* k)(n)| > 2A} 3 #{n2 |(b* k1)(n)| > A} + #{n: |(b* k2)(n)| > A} 2' IA+IIA- (2.1) Now, “(1* kllli’p llklllil Clip/2 [AS-TS ,, 3 AP. (2.2) Let p, = p 1 be the conjugate exponent of p. Then “12* k2||z°° S Ilkzllw S mil/q, since p < 2 hence, p’ > 2. Pick [1 such that curl/‘1 = A, i.e., u = ch‘q. Thus, #{nz |(b * k2)(n) > A} = 0, (2.3) which implies I I ,\ = 0. Therefore, combining (2.1), (2.2), and (2.3), it follows that . cup/2 c __ c #{n.|(b*k)(n)| >2A} 3 AP +E— E. . . 1 1 1 Hence, T 18 weak-type (p, q) for 1 S p < q < oo w1th 6:1-7— 5. By the Marcinkiewicz Interpolation Theorem, T is strong-type (p, q) for 1 < p < 2 and q 1 1 1 such that - = — — —. q p 2 38 2.2 A Lemma Involving Fractional Differentiation We now proceed with proving the main estimates. To this end, we make use of the operators D3 and complex interpolation. The following lemma provides half of the necessary results. Lemma 2.2 For 7 E R, - —1 4 i .) |th ” 7Hh(m)770||zgz;gh s Cllflollzg and ii) l |D,:1/2+”A,,w| luv» 11 m.h S C7l|wllly3lrlmha with C, = C(I’yl + 1) where c is independent ofcy and h > 0. Remark 2.3 Before proving (i), we check to see if the homogeneity is correct. As- sume (i) is true for h = 1. Then 1/4 Ithl/HWHMmMOHIgg-g, = (h 2(3111) ngl/4+'7Hh(m)7)o(n)l)4) nEZ mEZ 1 /h 4 1/4 = h sup — I/ . . 1/4 1 I 1r/h . . (h: 770(k)ezkh0)e-mh0 4 = h1/2 _/ h6 —1/4+:7 kez d0 ("26, (3:22 21r -./hls‘“( )' (1+Sgn(m)isin3(h0))|m| >4) "“7 (fi0)h(0)e“""" (1 + sgn(m)i sin3(h0))lm| sin(h6) h d0 1 " - _ - (mane-W = hl/2 _ f 6 l/4+ry ("26:2 (78:21; 27f -1r I sm I (1 + sgn(m)i sin3 0)|m| d6 s ch1/2(Z |17o(n)|2)1/2 n62 = 00% Z Ino(n)12)‘/2- nEZ 39 Proof(i): By Remark 2.3, the case h = 1 implies all other cases with h > 0. Hence, we can assume h = 1 and drop the subscript notation. Note that Z Z D’1/4+i7H(m)770(n)w(n, m) mEZnEZ 1r . * (9 _— = Z 217;] lsingl—1/4+n(1+s 770( ..) 3 lmla(1)(g,m)d9 mez —1r gn(m)ism 0) 7r “(1) = 2%] 70(o)(|sin6I-1/4-‘7 ” (0’7") )d9 mez (1 -— sgn(m)i sin3 0)|”‘| __— z Z 170(n)D‘1/4“7H(-m)w('am)(n)- mEZ nEZ To prove (i), by duality, we need to show that ll 2 D‘1/4+‘7H(m)W(nm)llzg S Clleg/sun, (2-4) mEZ with c independent of ’y. The square of the left hand side of (2.4) equals Z Z ZD—1/4+i7H(m)w(nam)D'1/4+i7H(u)w(n,1/) nEZmEZ uEZ __ 1 1r - —1/2 7;,(1)(g, m)_ 0 depend on m and V. By (2.5) and duality, we need to show that H 2 D_1/2Tm,l’w(ni V)||l;‘,l$,? S Cllwllzi/flgn' (2'6) VEZ Now, 2 D-l/ZTm,,,w(n,z/) = 22—1-e-W1sin01 1/2T(m, u)(0)a<1>(0,u)00 V62 V62 —7I = 22—1]: e””9| Sln 6| 1/2T(m,) V) (0) (Z w(k,V)eIk0d0 VEZ 2’” 1:62 5 Z Z 271/; a“: ">91sin01 1/2T(m, u) )(0) 00] 1100. u) )1. kEZVEZ Claim 2.4 C., in9 1/2+i7 l/tre |sin0| T(m, V) (0) 0 such that |¢Z,n(0)l Z c Inlp (2-10) 44 suppC1C{0:c1‘lI—:—I S 1913 ”III—ELI} (2.11) Before proceeding, we need the following lemma [10, pg. 342]. OD Lemma(van der Corput) Suppose 45 is C°° and compactly supported. Let 1/1 be a real-valued function so that, for some I: E N, k 2 2, 10%| 2 C,c > 0 throughout the support of 01. Then 00 , C l/;0061¢($)¢($)d$ S WI(II¢IIL°°+IID¢IIL1) Applying van der Corput’s Lemma to [5,12, with k = 2 and appealing to (2.10) and (2.11), it follows that I“) < ___—9— F .. f F’ 0 d0 , 2.12 "17-1n11/4p1/4 (ll 11. (supp(1)+ suppcll <11 ( 1 with 0(0) z Isin01‘1’2+"’/\1(0)_ cp02 on supng. This and the fact that |w;,n(0)| S cp|0| implies C plgl _ C — 1512) < 1/2dO 0+____ 3/2d9 ’7 " lnl l/ln |<101p02I0I Inl 1/lnl<|9|I0I + c [04 1.1010011n101czw1 0,, [ml/2 _fl/4 (1+sin60)q C. S InII/2 Since |I,(,2,),| g c, this finishes the proof of (2.15) which implies C In <—7—-. ’7 _ |n|1/2 + 1 The estimate for I I Inn can be obtained from the estimate for In” using a change of variable, namely r = 0 — 71. These facts and (2.8) conclude the proof of Claim 2.4 and hence concludes the proof of Lemma 2.2(i). Proof:(ii) We need to show that 11D;1/2+‘7A,,w111.,.. umh S CIVIIwIIlz/i’l'ln'hv with C, = c(I'yl + 1). Suppose m > 0. Then, by Claim 2.4, |D;I/2+i7Ahw1(n, m)| h3 —inh0 h:_ 271' 1./—://he _ hl/2 “”7 0111110 01 (1 + sgn(m +1 - j)isin3(h0))lm+1-J'I ~ 0|—1/2+i’7 h3(k_2}1_r_/:I i(k— n)0 I SlIl d9 2 Z( w j) e (1 + 0sin3(0))m+1-J kEZ j=1 sin(h0) h d0 47 m ' 0|-1/2+i'7 < hl/2h3 (k ei(lc— n)0 ISln d0 - £32,212”); (1+0s1n3(0))m+1—J 1 g C, 111/2113 |w( (k ,j)| , BI: — 01/2 +1 with C, = c(I'yl + 1) and independent of m. The same proof works for m S 0. Therefore, 1/2 +1”) 7/2 1 SUP ID), Ahw(n, m)| < C h Iw( (k III . mEZ £3.26; lk _ ”ll/2 +1 Hence, by Lemma 2.1, IID;I/2+I7Ahwlllgl:h = (h(z:(sup lD-1/2+i7Ahwl)4)1/4 I n62 mez S 0101/42/28: |w( 71021112313“ 1162 mEZ ‘7 C7IIwIIgg/31'1n h and this concludes the proof of Lemma 2.2(ii). 2.3 Proof of Lemma 1.17 The following lemma implies Lemma 1.17. Lemma 2.6 For 7 E R, 1) IIDh IIIIHh(m)770II1°°12, < CIIIIOIII}: and i2) IIDIIIIAIIWII100:2 S CVIIWII11121 "111,11 ’1th with C, = c(IyI + 1) and c independent of”) and h > 0. 48 Remark 2.7 To see how this implies Lemma 1.17, first consider Lemma 1.17(ii) which says IIaimAthIzg-ezfm S CvIIwIIIngn‘h- The multiplier associated with 6,2”, is — sin2(h0) (1.2 I which is a constant multiple of the multiplier | sin(h0) |2 7,2 associated with D,2,. Hence, Lemma 1.17 (ii) reduces to Lemma 2.6(ii) with 'y = 0. As for Lemma 1.17(i) which says IIan,hHh(m)770IIz;-,ozfmh S CIIUOIIzg, the multiplier associated with 0",), is —i sin(h0) h while the multiplier associated with D}, is | sin(h0)| _h . Let M(h19) = —isgn(h0) for 0 6 [—7r/h,7r/h]. Hence, —i sin(h6) h |sin(h0)| -_— M(h0) h 49 which implies for m E Z 3n,hHh(m)no = DIHhThflo, where T), is the Fourier multiplier on [—1r/h,1r/h] associated with M (h0). Since T), : 1,2,(Z)—1 lflZ) is bounded (with norm less than one for all h > 0), Lemma 2.6(i) implies Lemma 1.17(i). Remark 2.8 Before proving Lemma 2.6(i), we again check to see if the homogeneity is correct. Assume (i) is true for h = 1. Then IIDIIIIIIHMmIUOIIu-szghh |/\ . 2 1/2 1 “/1 ... sin 0. Before we prove Lemma 2.6(ii), we define three functions on Z x Z. Definition 2.14 h3 )3 Hh(m+1—j)w(-,j)(n) ifm>0 j¢m+1 Ahw(n, m) = 0 ifm = 0 -h3 2: Hum —1—j)w(-,j)(n) ifm < 0. #m—1 Definition 2.15 h3 2 sgn(m +1 —j)Hh(m + 1 — j)w(-,j)(n) ifm > 0 j¢m+l Bhw(n,m) = 0 ifm = 0 ha 2 sgn(m —1—j>H,.(m —1—j>w(-.j) ifm < o. j¢m-1 Definition 2.16 h3 2 Hh(m+1—j)w(-,j)(n) ifm>0 11 jz-oo Ehw(n, m) = 0 ifm = 0 —h3 :0: Hh(m — 1 — j)w(-,j)(n) ifm < 0. It follows that Bhw = 2Ahw — Ahw + 2Ehw, 55 which implies DfimB w— — zpfimA w— D§+‘7A,.w + 2D§+i7Eh u). (2.21) Proposition 2.17 HDZHVEIIW I [13°12 m, S 0110011511an with c independent of 7 and h > 0. Remark 2.18 For the last time, we check the homogeneity. Assume Proposition 2)1/2 2.17 is true for h = 1 and let m > 0. Then 2+” wil’wJ) (1 + isin3(h0))m+1*j h49) sin d0 ( h h3 i 1 lfl/he ~inh0 27f —1r/he j=-oo sup(h3 Z 1 Z27r jz-oo A (1) 0 . —in9 2+i'y w ( 3 .7) A. ‘3 |s1n0| (1 + isin3 0)m+1-J' d0 m>0 = h3/2h3h—2 (Z 1/2 S chm/2: Z:|(.u(n,m)|2 nEZ mEZ = own...“ The same proof works for m < 0. Combining these and the case h = 1 implies all other cases for h > 0. Proof: Assume h = 1 and assume m > 0. Then D2+i7Ew(n, m) D1+17H(m)[H(1)(Z(D1H(—j)[x{jso}w(,j)])] . jEZ 56 This implies ||X{m20}Dz+i7Ew||z;-.oen S C 11(1)(ZDlm-J'qusowhfll) (2-22) jEZ (’2‘ S CHZD1H(—j)[X{jgo}w(°,jlllllg (2°23) jez S C||X{mgo}w(n,m)||z,1,zz,, (2-24) —<— cllwlllhlgn: where (2.22) follows from Lemma 2.6(i), (2.23) follows from fact that H (1) is bounded on l2(Z), and finally (2.24) follows from Proposition 2.13. The case m < O is similar. Combining these two cases and applying the triangle inequality finishes the proof of Proposition 2.17. Proposition 2.19 “DizzflvAhWngozfmh S Cllwllzngn'h, with 6 independent of ”y E R and h > 0. Remark 2.20 If the set of operators {Hh(m)}m€z did form a group under composi- tion, i.e., Hh(m1 + m2) = Hh(m1)Hh(m2) (2.25) for all m1, m2 6 Z, then the proof of Proposition 2.19 would proceed exactly as the proof of Proposition 2.17, which used the fact that (2.25) is true if m1 and m2 are the same sign. However, because of our choice of discretization, we do not have the group structure in general. Hence, our proof will require different methods. 57 Proof: By a similar homogeneity argument as in Remark 2.18, we can assume h = 1. Suppose m > 0. Then D2+'7Aw( (,n m) A (1) . w (0,2) d0 j¢m+1 OO 1 —in0 2+i 0 7 . :2 [22:6 (Sin I (1 + sgn(m + 1 — j)isin3 0)Im+1-2I = lim—/ If?! e—mol sin 6|2+i7(1 - isin3 0) < Z (1 — isin3 6)m—jcb(1)(0,j)) d0 < <1r—c J c—>0 271’ c -in0 1 Sin 0|2+i7 '=m+2 WW 2') + lim / e e—>0 27r <|9|<1r—€ | sin 0|2+i7 m _fi— 1+is1n 0 (J.__ ”(1+isin3 0)” | sin 0|2+'Ve" ‘2‘P J)do = 1 — . 1r . m 2 6""0/ e‘mup . , 3 . + c—>0 47r e<|9|<1r-e —1r 1 + i s1n l9 — e"? 1— isin36 — e“ I Sill 0|2+i7e—i2

0, m 6 Z, and 7 E R, set K,(m, 0) = 0. Hence, £216 Ke,7(m1 (10) : K7(m’ (p) for all m E Z and (,0 6 [—7r, 7r]. Claim 2.22 There exists C such that ||K€,7(m1')HL2[—1r,1r] S Ca with C independent of e, m, and '7. Proof:(of Claim 2.22) By definition, . 1 7’ _,- K.,.(m, -) (z) = gJ-J [_J K.,. 0 27f /€<|9|<7r—e e (1 + isin3 9)l+1 ‘ — 0 ifl=—1 1 _- 9 |sin0|2+i7 , -— "n d0 fl<—2. 27r £<|0|<1r—ee (1—isin3l9)(‘(“1 l — 59 By Plancherel’s Theorem, 1/2 (lKe.7(m1')(|L2[—1r,1r] = (Z lKe.~y(m, -)‘(l)l2) (62 |/\ 1/2 1/2 (J2; (KW )(lzfl)+ (;IKW )(l) l 2) (2-26) First, consider the second sum in (2.26). Assuming e is small, 1/2 1 Isinl9|2H7 2 1/2 , .. 2 : _ —im9 (21mm, > ml) (; 2.12mi (1 ”2.3121“ 1 130 . 1/2 1 ‘ 024-17 2 S X -/ e ""9 (3191, may 120 277 6<|9|<7r/4 (1 +ism 0) . 1/2 1 - |sin0|2+'7 2 + — f "m" . . d0 (5 271' 1r/4<|9|<31r/4 e (1 + is1n3 0)1+1 . 1/2 1 _' 9 lSln0|2+t7 2 + — / "" . . d6 (g5 27’ e0 (1 + Sln6 0)(+1 31r/4 1 d6 1/2 — C (fr/4 sin20 ) S C, with C independent of m, c,and 7. Next, turning to I, it follows that 2 1/2 120 1r/4 1 f . . 3 ,f.,m,.(0)d6 —7r/4 (1 + ism 0) 60 e—imOI sin 0l2+i7 1+isin36 "/4 If... ((9):? ”2 «x4 0 1 1/2 I < / —’—’7-—d0 / - C ( -1r/4 sin2 6 C 4/4 1 + sin6 6 sin2 00w —n/4 1 +sm 0 C, with f€.m,.,(6) = X{¢<|o|<7r/4}. By Corollary 2.12, |/\ l/\ with C independent of m, e, and 7. Furthermore, III can be reduced to I (but with 1 — z'sin3 0 in place of 1 + z'sin3 0 which can be handled similarly to I). The first sum in (2.26) can be bounded by C independently of e, m, and 7 using similar arguments as above. /// Claim 2.23 lim Km(m, ) = K1(m, ) c—iO in L2[—7r, 7r], uniformly in m. Remark 2.24 Assume Claim 2.23 for the moment. If w E lll2 then nm’ lim 6—H) f... 64'” (2mm - 33 so) — K702 — 3', WW 90)) ds0 jEZ €—)0 5 lim / :3 IK.,.,(n — j, .p) — K.(n — j. (p)||o(2)(j, sandy) '"jez £401.62 S £136: H Ke,7(n “ j: ) " K‘r(n " j, ') ”L2[—1r,1r] ' H 03(2)“) ) ”UH—war] jEZ =0. = lim 2 j... IK.,7(n - j. w) - K702 - j, w)ll<51(2)(j, w)ld

61 > 0 with 62 < 7r /4 and assume I _>_ 0. Then (K62.’7(m’ ) _ K61,'7(m7 ))A(l) 1 _- 9 |sin6|2+i7 ;— — "n d6. 2.28 277 (v/61<|9l<€2 + v/17—62<|9I<1r—61) 6 (1 + 2181113 0)H'l ( ) By Corollary 2.12, 2 1/2 . 2 1/2 2 d0 5 c / .3;ng 120 €1 0, uniformly in 177.. Since Km(m, .) —-) K7(m, ) pointwise, K6,,(m, ) —) K7(m, ) in L2[—7r, 7r] uniformly in m. Furthermore, it follows that ( i f" —im9 ISin0|2+i7 d0 fl > I 21r —1r e (1+z'sin319)’+1 l " 0 Kv(m. -)‘(l) = 0 ifl = —1 (2.29) 1 7' _- 0 |sin0|2+i7 . — ”’1 d0 fl<—2. 27r/_1re (l—z'sin30)|‘|—1 ‘ — /// 62 Claim 2.25 The Fourier coefficients of K0(0, ) are absolutely summable, i.e., Z|K0(0,)1)| < C. (EZ Proof: Suppose *y = 0 and l 2 0. Let 1/11 E 08° with 1/21 E 1 on (—7r/8, 71/8) and suppwl C (—7T/4,7T/4). Let ’l/J2(0) = 1,!)1(0 — 7T) and $3 =1- 1&1 — 1/12. From (2.29) it follows that "/4 sin20 Ko(0,‘)A(l) = 0 121(0) . . 3 d9 —1r/4 (1 + z Sln 6)‘+1 sin2 0 0 d0 ) (1 + isin3 19)1+1 +C/ 1P3 11/8<|0|<71r/8 51r/4sin2(5l 116 +6 3.. /. M0 )(1+z'sin3 11)!+1 = 1(1) + 11(1) + 111(1). First, sin2 0 0 d9 ) (1 + z'sin3 0)‘+1 1 d6 /1r/8<|0|<71r/8: _0 (1 + sin 60)”2 :Imm = i (=0 3 [11/8 0 such that for :c, y E H, lflfl-fWMSCW-m2 with distance measured on the unit circle. Proposition 2.27 Let 0 < a S 1 and suppose f is a continuous function defined on H which is differentiable for :c E H\{O}. If there exists C such that IfWHSCW”5 then f E Lip°(H). 64 Proof: We need only concern ourselves with a: and y near the origin, say, 2,3; 6 [—7r/4, 71/4]. Suppose |:1:| < lyl. If neither :1: nor y equal 0 and they are on the same y 1 side of the origin, then f (y) — f (z) = :1: f f (t)dt, which implies |f(y) -f(-’r)| |/\ y 1 if 111mm :tC [y ltla‘ldt 0.01/1“ — 1x10) Caly — ail". |/\ |/\ The last inequality follows from the fact that the function no, (2:) = lea E Lip"(H). Now, if :1: = 0, then f(y) — 1(0) = 21/07th since f' has an integrable singularity at the origin and f is continuous there as well. Hence, the argument above works when either a: or y are zero. Finally, suppose :1: and y are on Opposite sides of the origin with :1: negative. Then lf(y)—f($)| S |f(y)—f(0)|+|f(0)—f($)| g o [0” (110-111 + c [0 (1101-1111 S Ca(|y|a+|$|a) S éa|y+lrvlla = Caly-xl“, where the last inequality comes from the equivalence of the l1 and 11/“ norms on R”. This concludes the proof of Proposition 2.27. 65 Recall, h 1' secz'r 1 gm = (“f _f l) 2 tan 7’ \/1 _ tan2/3 7. = 91(7) '92”), with h(r) = arcsin(tanl/3 7'). Since 1121 is smooth, compactly supported, and 1111 E 1 in neighborhood of the origin, it follows that 91 E 08°(I'I). Next, 1 92(7) = —§(1—tan2/3 T)—3/2 tan‘l/3 7 sec2 ’7', which implies 193(T)| S CITI‘”3 0n SUPP¢1(h(T))- Corollary 2.28 g E Lip2/3(l'l). Proof: Since 91 is smooth and |g'2(r)| g Clrl‘1/3, we have lg'(T)| S CITI'l/B- Since 9 is continuous, Proposition 2.27 implies that g E Lip2/3(l'l). /// The following is a theorem of Bernstein [4, pg. 32]. 66 Theorem 2.29 If feLme for some a E (1 / 2, 1], then the Fourier coefiicients of f are absolutely summable, i.e., Z If(n)l s G. 1162 Corollary 2.28 implies that Zlakl <00 kEZ where ak is the kth Fourier coefficient of h 7' sec2 1 1 flfl=wfip 2 tan 7' \/1 _ tan2/3 7. = 91(T) '92(T)o By the proof of Proposition 2.9, (2.31) becomes 1 , (2.32) m 2 31.1021 (:0 0 j: with Bid = $6). After interchanging the order of summation and applying Claim 2.11(ii), (2.32) becomes 2 NW _<_ 2 lilsz Z 31.1 i=0 1:1 (:0 w 2 2 Mil i=0 C 2 lajl jez c. (2.33) l/\ |/\ 67 Again, by a change of variables, 2: II I I (l)| S C. With (2.30) and (2.33), this implies (=0 2 moo->101 < ca 120 The proof showing 2 |(K0(0, ))“(l)| < 00 is nearly identical. Therefore, (<0 :22 |Ko(0. -)‘(l)l < oo. (2.34) which concludes the proof of Claim 2.25. /// Finally, to finish the proof of Proposition 2.19, observe that (2.34) implies K 0(0, (p) is a continuous function of p, which implies that there exists C such that ”(0(0: 90)] S C- (235) Claim 2.30 There exists C such that |K7(m. <19)| S C. with C independent of m, (,0, and 7. Remark 2.31 Claim 2.30 will finish the proof of Proposition 2.19. To see why this is the case, assume there exists a C such that |K,(m, 90)] g C independently of m, (p, and '7. Hence, by (2.27), I... 8—... 2 mm — 1,2)a<2>(1.2|singp|} lsin3 0| 2 0 __ 2 g C/ 0—d0+C/ ( 7r>111) {0:10:3chva Icpl {ale—«191211131 lrl /‘ lgo|s1n 0 {0:I sin3 0|>2| sin gal} ] sin [p] 02 g C / d0 + C {9 101 0. Proof: By a similar homogeneity argument as in Remark 2.18, we can assume h = 1. Suppose m > 0. Then _ . A(1) . D2+"’Bw( (,n m) =2 —/: e""9|sin0]2+iv sgn(m+1 le £911) 1 .d6 J';r£rrz+127r (1+ sgn(m +1 — ])ism 0)lm+ “3| 70 1 . ' 62+” m “(1)0 ' = lim—/ e—zn0_l_31_n._l.__3__ w. .(373) . d9 HO 27r £<|0|<1r-—c 1+zsm 0 j=_°o (1+zsm 0 "“1 , 1 - . - . . °° . . . A . —£1_1+16 -2—7r/€ 0. Let (0 E C8°(R) be an even function so that Co E 1 for |a:| _<_ 2 and suppCo C {le S 3}. Also, let (1 = 1 - Co- Then I (m (.0) < / e-‘marcsin1tr1€'<|x|<1..+..( d2: r +92: +— lyl r, — 9x2 2dx 1 C _ + M (r' + 9x2) dx |/\ |/\ cos(xy) + c [mbfig'defia c 1 WI Izl>¢1 1172 C, |/\ |/\ where C and c are independent of e', m, a, and r' (c is the boundary term). /// Claim 2.37 There exists C such that llgwllu S C(I7I + 1), with C independent of 7 and T(small). Remark 2.38 Since llfe’m,a,r' * g7,TllL°° S llfe'm,a,r'llL°° ° llg'YflllLlr Claims 2.36 and 2.37 imply that lfe'm,a,r' * 9717(ma’)l S C(l’Yl +1) With (2.43) this implies Agmmfl) g C(|7| +1). (2.44) 77 Proof:(of Claim 2.37) Assume r is small. Recall, g (x) z Ix + 1|“’[1 — «)(n: +1)1 1). Thus, lg'm (y) Id?! 1 . d2 _ ’33!— " y. [Re d$2(g%.,)(x)dx 1 CmaX(|7|. 1) + 011180107th / —dy lyIZmaX(l7l.1) y2 g ,T y dy = / [RI 7 ( )I IyISmaxmm) dy +/ lyIZmaX(|7|.1) |/\ |/\ CmaXUVI. 1) |/\ C(|7| + 1). and this concludes the claim. /// 78 Finally, assuming m 75 0 and by change of variables, A:.5L(m, 7') becomes / e—imarcsin(rl/3[x+1])3lx + 1|i7x (0(5):) d2? max(e',1#[)<|x|<3 r' + 9:1:2 \/1 _ 7.2/3“, + D2 = / e—im arcsin(rl/3[x+1])3$l$ + ”in/)0: + 1) (0(3) d2? max(e’,m)<|x|<3 r' + 9x2 \/1 _ 72/3“: + 1)2 + / e—imarcsin(rl/3[x+1])3$lx +1li7l1 _ WI” + 1)] (OCT) d$ max(e',I;10—r)<|x|<3 r' + 9x2 \/1 _ 7.2/30; + U2 1 - r2/3 x +1 2 S C+c / \/ 1(3 ) max(e',]—mi;[)<|x|<3 mr / i (e—imarcsin(rl/3[x+1])) 113 III: + ”fill _ ¢($ + 1)]C0($) d2? dx r' + 9x2 \/1 _ 7.2/3($ + 1)2 < C + ___C_ / e—imarcsin(rl/3[x+1]) _ ImT1/3I max(e’.I.—.‘;I) 2}, it follows that dx (1) m T c |7"11=2(II=3 - 1) - r'(II: -1)| B" ( ’ ) S [Imlzz (1‘ - 102(1133 - 1)2 |$2($ -1)2($3 - 1) — (a? - 1X51?3 -1)2| “($.22 (x — 1)2(x3 — 1)2 dz 5 6 x x c/ I—I-dx + c/ —8dx IIIZ2 17 le22 517 C. (2.48) |/\ |/\ / (e—imarcsin(rl/3x) _ eixma\ I$Ii7($ _' 1) C1($ — 1) I2I5391-1/3.c' 0. To prove (2.55), consider the analytic family of operators rznom) = D;‘/4D§,‘—Z>Hh(m)no(n), o g Rez g 1, with 770 E lflZ). If 2 = i7, then IszUOIllgOIEn'h = IIDl—iS/Mflollzooz? < Cllflollzg: (2-56) nm,h_ by Lemma 2.6(i). If 2 = 1 + i7, then IITz770||141°° = IID_l/4-i5/47Hh(m)flolItglsgm S CIIflong, (257) n m.h by Lemma 2.2(i). Combining (2.56) and (2.57) with Stein’s analytic interpolation theorem [7] and letting x = 4/5 implies (2.55) which finishes the proof of Lemma 1.18(i). 84 Next, recall the statement of Lemma 1.18(ii), llAhwllzgzggh S Cllwllgmwf) (2-58) with c independent of h > 0. Similar to the proof of (2.55), consider the analytic family of operators TM = Dgz/Zpfill‘zlAhw, 0 g Rez g 1, with (.0 defined on Z x Z and compactly supported. If 2 = i7, then “1.421143”, = ”Dr‘s/24141141,, s 0.111111%“ (2.59) by Lemma 2.6(ii). If z = 1 + 27, then IITzwIIIgzsgm = IID-lfl—iszthzgzgh S Cvllwlllg/fi'lmh, (2-60) by Lemma 2.2(ii). In both (2.59) and (2.60), C, = c(1 + |7|) with 0 independent of 7 and h > 0. Again, applying Stein’s analytic interpolation theorem finishes the proof of (2.58) which concludes the proof of Lemma 1.18(ii). CHAPTER 3 Numerical Results In this chapter, we discuss the numerical implementation of the fixed point iteration with the operator (Pm. 3.1 The Cutoff Function 7h(n,rn) The iterates of the contraction mapping (Pm are defined on the entire grid, Z x Z. To go from one iteration to the next, the entire iteration is needed. This is obviously not feasible numerically. Hence, we need to introduce a “smooth” cutoff function which is zero for n or m large. One may think we can simply set the iterates equal to zero for large n or m. However, this may adversely affect the norm estimates of the contraction mapping. This is due to the fact that if 02 is defined on Z x Z and 1ifIn| SN,|m| 3M XN,M(n1m) : 0 else, then it is not necessarily the case that llan.h(XN.Mw)II(;-,°zghh S Cllanmwlllgwgha 85 86 with c independent of h > 0. In fact, it is quite the contrary as h —) 0. Thus, we need to introduce a cutoff function which is one for Inl S N and |m| S M, but decays slowly to zero for large n. Let h > 0 be small. Suppose we would like a solution of (K dV): for (n, m) E Z x Z such that |nh| S N and Imh3| S M where h is chosen so that N is an integer multiple of h and M is an integer multiple of h3. Let 7;,(n, m) be a piecewise linear function such that 7,, E 1 for InhI S N and |mh3I S M and supp7h = {(n, m) E Z x Z: Inhl S 2Nh"1/5, Imh3| S M}. Obviously, 7), can be chosen so that chl/S sup lamh'th S N n,meZ (3.1) for small h. Note that we do not need any decay in the m—direction since none of the three norms of X h involves differences in the m—direction. Definition 3.1 For a discrete function 77 defined on Z x Z, let ~ 877(77’1 m) : 711(7),, m)"(") m) Recall the definition of (Pm, 1412(1), m) = Hh(m)7?o(n) — éAhanmmn, m). Wlth 770 E 1,2;(Z). 87 Theorem 3.2 There exists ()0 > 0 and r > 0 such that, if B, = {w E X), : lelx. < r} and “none, 3 60, then i) S‘Dm : B, —+ B, continuously and ii) |IS,,O(V) — Sm(u)IIxh < AIIV — ,uIIXh, for some A < 1 and 14/1 6 B,, with A, 60, and 7‘ independent of h > 0. Proof(i): First, we would like to determine a bound for the operator norm of S on X h. Obviously, llgflllzgzfih S Il’lllzgzggh (3-2) and sup ||§17(-.m)llzg S sup ||n(-.m)|lz,2,- (3-3) mEZ mEZ Note that 0n,h(7hn)(n, m) = 7;,(n + 1, m)3n,h17(n, m) + 8n,h7h(n, m)n(n — 1, m). Then 1/2 llan,h(nn)llzgo)2 < llamhnllzsgzgn, + 8161121 (ha 2 (6.171(1). m)|2|n(n -1.m)l2) nah "- mEZ , 2/5 1/10 E IlanmnlI123)?“+ 51113013 Z Ian,h7h(n) m)|”) (ha 2 WW — 1) m)|10) (3-4) "6 mlSMh-3 mez 1/10 ch1/5M2/5 S Ilan,h77llz;,-o13n h + _—N—_ SUP ’13 Z W71 — limlllo (3-5) ' nEZ mez nEZ mEZ 1/2 1/5 S(1+5)||?7||xh, 88 with 5 depending on N and M. Note that (3.4) follows from H61der’s inequality and (3.5) follows from (3.1). By this estimate, (3.2), and (3.3), we have lls‘nllxh S (1+ é)||17||x,.- (3-6) Recall that from the proof of Theorem 1.21, we were able to show that I|,,077||Xh S C5,, + Cr5 for 17 E B,. By choosing r and (50 small enough, we have 1. c60+Cr5S . 1+ 6 Hence, (1511.47))...s(1+&)1+E =r for 77 6 B,, which implies 30,, :B, —) 3,. Proof(ii): Since S is linear, (3.6) implies that if [1,12 6 Xh, then ”59%)” — S‘I’nol/llx). S (1+ 5)||‘I’nou - ‘I’nonlxh- (3-7) From the proof of Theorem 1.21(ii), if 11,11 6 B,, then ll‘I’nou - (Paul/Hm. S C(“’llu - Vllxh- (3-8) 89 1 If we choose r small enough so that Cr4 < Till—5’ then (3.7) and (3.8) imply llgéflof‘ _ 3(1)?)onth 5 All” — VIIXM with A = Cr4(1 + E) < 1. Remark 3.3 In the definition of 7),, if we fix m E Z such that Imh3| S M, then 7h(n, m) = 0 only if InhI > 2N h‘1/5. This was necessary to insure (3.1), which led to (3.6) independent of h > 0. If we assume h is bounded below, then the h‘l/5 can be discarded, i.e., we can select 7,, with 7,,(n, m) = 0 only if |nh| > 2N. Hence, S,,0 is a contraction mapping on B, C X), for some 7‘ > 0. By the Contraction Mapping Principle, there exists a unique 17 E B, such that S (1)1107? = 77 which implies VhQnofl = 77 Let A = N/h E N. If (n, m) E Z x Z such that InI S A and |mh3| S M, then c11.401111) = 7701.171)- (3.9) By definition of the Operator (Pm, (3.9) implies that 77 solves (KdV): for —A + 3 S n S A — 3 and |mh3I S M. Since the definition of 62,,,w(n, m) includes w(n + 3, m) and w(n — 3, m), it is necessary to have n satisfy —A + 3 S n S A — 3 so that when the linear operator (as in (1.1)) is applied to both sides of (3.9), we maintain equality. 90 3.2 The Coefficients Q[n,m] In this section, we now turn to the actual implementation of the fixed point iteration. Fix h > 0 and let 770 E 1,2,(Z). Recall the definition of the operator Hh(m), Hh(m)no(n) = 1 [Wk e“"h9(fio)h(9) d0 2_7r —1r/h (1 + sgn(m)i sin3(h0))l’"| 6 1 [WM e—inhO -imarctan(sin3(h9)) thez 770(kleikh9 ‘ 2; ./). (1 + sin6(h0))ImI/? h w/1 . . . 3 : k _/ —i(n—k)h0 —imarctan(sm (h0)) (10 £0“ )21r —1r/h e e (1 + sin6(h0))l”"l/2 1 1f . . - 3 1 = k _/ -:(n—k)0—zmarctan(sm (0)) d 262’“ I27. .. e (1 + sin60)lml/2 1 7' cos((n - k)6 + m arctan(sin3(0))) = Z 770095;] (1 + sin6 6)ImI/2 kez ‘” Definition 3.4 For n, m E Z, let _ 1 1' cos(n0 + marctan(sin3(0))) an, ml - Er- _W (1 + sins 0)Im|/2 d6. It follows that Hh(m)770(n) = : 110(11)an — k, m). (3.10) kEZ We can represent Ahw(n, m) using these coefficients as well. Note that QIn, m] does not depend on h. These coefficients can be computed and stored for future use. To begin, we choose the zero function as our starting point. Independent of 1', this guarantees that we are starting in 8,. With Q[n, m] precomputed for n and m 91 such that Inhl S 2N h"1/ 5 and Imh3| S M, we compute the first iteration, (’01. Let wl (n, m) = (1),,0 (0) (n, m) = Hh(m)770(n) - 14.6....(oxn. m) 5 = Hh(m)770(”) = £7700“)an - k) m]. It follows that (01 solves the associated discrete linear equation. To obtain the first iteration, we need to multiply by the cutoff, i.e., 601 (n, m) = 7,,(n, m) -w1(n, m) = 3(1),“, (0)(n, m), which still solves the linear equation for n and m such that —A + 3 S n S A —— 3 and |mh3I S M. Since (721 is the solution to the associated linear equation, we check the computation’s accuracy by means of maximum relative error, max max InISA—3 0,,0(IJI(n, m) = SQmw1(n, m). 92 By induction, (I'Jk(n, m) = 7h(n, m) - ,,0(12k_1(n,m) = $0,012,-de for k 6 N. At each stage, we compute the maximum relative error by replacing (221 with a), in (3.12) for m > 0 and m < 0. Finally, we compute the Xh-norm of each iteration. Recall, this involves three size estimates, 1/5 1/2 Ilaklll=llakllzgqgh= h 2 (’13 Z l@k(n,m)llo) , In|S2Nh—3/5 ImISMh3 1/2 ll‘IJkll2 = llamhaklllg'flfim = lnlgmgm (’13 Z lan,h0k(n,m)l2) , lmlSMh" and 1/2 lldklls = SUP H5110 m)||,2 = max ’1 land" m)|2 - mEZ , h ImISMha Inls2§1—6/5 , 3.3 A Few Examples In this section, we give a few examples of the convergence and divergence of the operator ,,O. Using Mathematica, we have computed Q[n, m] for |n| S 32 and Iml S 32. Hence, our iterates will be defined on InI S 16 and ImI S 32. In the following examples, the step size h will be no smaller than 0.1. By Remark 3.3, we can set 1 lnl s 8. Iml _<. 32 3+2 —16SnS —8,|m| 332 7h(n,m)= —Lg-+2 8SnS16,ImIS32 0 else. 93 This implies that we only consider the points {(n, m) : In| S 5 and ImI S 32} for the solution of (K dV)51". Since we are limited in the number of coefiicients available, we will choose our initial data in each example to be compactly supported near the origin. Example 3.5 Suppose our initial data is given by the function 0.5x + 0.5 —0.5 S x S 0 110(1) = -0.5x + 0.5 0 g x g 0.5 0 IxI 2 0.5. Note that an is an even function. It is easy to see that if w(x, t) solves the associated linear problem (H) (see pg. 4), then so does w(—x, —t). Furthermore, if u(x, t) solves (KdV)4, then u(—x, —t) does as well. By uniqueness, u(—x, —t) = u(x,t). In fact, this symmetry prOperty should be true for all the iterates of the operator 5' in (3), since the right-hand side of (I H ) (see pg. 4), namely g(x, t), will have the property that —g(x, t) = g(—x, —t). Recall that we choose the zero function as our starting point. It follows that 5(")(0)(x, t) = SI")(0)(—x, —t) (3.13) for all x,t E R and n E N. Let h = 1. To begin, we need to discretize our initial data. We do so by choosing the values of uo at the integers. Thus, our discrete initial data is 0.5 n = 0 1M") = 0 n 75 0. In this example, we run four successive iterations. Recall that the first iteration solves the associated linear equation. To measure the accuracy of the computation, we compute the the maximum relative error of the first iteration with respect to the 94 linear equation using (3.11) with A = 8 and M = 32. The result is 8.62601 x 10‘“. All three norms of each iteration and the maximum relative error (M RE) of each iteration with respect to the (K dV)§I equation are listed in Figure 3.6. For the last three iterations, we include the distance of the current iteration from the previous one measured in the Xh-norm under the column with heading “distance”. As in the continuous setting, we have (3%)(I‘I(0)(n,m) = (§¢no)I’°I(0)(-n. -m) (3-14) for all n and m on our grid and k 6 N. One can see from the norms that these iterations are converging very rapidly. Since the relative error is becoming extremely small, we can conclude that the iterations are converging to a solution of (K dV)j,‘I. iteration M RE norm one norm two norm three distance 1 0.00094209 0.516778 0.409208 0.5 2 0.0000615788 0.516887 0.408461 0.5 0.03331 3 3.87118 x 10‘7 0.51688 0.408455 0.5 0.0000105811 4 2.43126 x 10‘9 0.51688 0.408455 0.5 5.08566 x 10‘8 Figure 3.6 Example 3.7 Let the initial data no be as in Example 3.5. However, in this example, we set h = 0.5. If no(n) = uo(nh), then 0.5 n = 0 77001) = 0 n 75 0, which is the same discrete initial data as in Example 3.5. One might think that we should get the same iterations as well. However, recall that the operator A), is 95 dependent upon h. This influences each iteration. As in the previous example, we compute the first four iterations (starting from the zero function) and their relevant data (see Figure 3.8). Once again, these iterations have the property mentioned in (3.14). By considering the distance estimates between successive iterations and the size of the maximum relative error, one can see that the iterations are converging exponentially to the solution of (K dV);d guaranteed by Theorem 3.2. See Figure 3.8. iteration M RE norm one norm two norm three distance 1 0.0188406 0.365417 0.289354 0.353553 2 0.0000305815 0.365431 0.289222 0.353553 0.0231647 3 4.78587 x 10’8 0.365431 0.289222 0.353553 4.6644 x 10‘7 4 6.90206 x 10"10 0.365431 0.289222 0.353553 7.68505 x 10-10 Figure 3.8 Example 3.9 Suppose our initial data is given by the function 1 —1SxS0 110(23): —1 0SxS1 0 IxIZl. Unlike the previous example, this function is discontinuous. However, it does have finite left and right-hand limits for all x 6 R. Assuming h = 1, let 1 170(7)) = 5 (5133» 110(1). + a) + 5133+ u0(n — 5)) , 96 i.e., 0 5 n = —1 0 n = 0 770(71) = —0.5 n = 1 0 other for n E Z. Also, note that uo is an odd function. In this case, if w(x, t) solves the associated linear equation (H), then so does —w(—x, —t) and if u(x, t) solves (K W )4, then so does —u(—x,—t). Hence, u(—x, —t) = u(x, t). By a similar argument as before, this is also true of all the iterates of the operator S. As in the previous example, our initial guess is the zero function. However, in this example, we compute the first seven iterations. This data illustrates the odd reflexivity noted above. Notice that for iterations five through seven, the only column which changes is the “distance” column which becomes extremely small (see Figure 3.10). This implies that by the fourth iteration, we are very close to the fixed point of our operator. Apparently, the best we can hope for in this example in terms of the relative error is 1.04304 x 10-8. iteration M RE norm one norm two norm three distance 1 0.0099275 0.587994 0.737931 0.707107 2 0.000044212 0.588019 0.738561 0.707107 0.036962 3 1.35837 x 10‘7 0.588019 0.738563 0.707107 5.61944 x 10‘6 4 1.04304 x 10‘8 0.588019 0.738563 0.707107 2.61774 x 10"8 5 1.04309 x 10‘8 0.588019 0.738563 0.707107 1.37671 x 10"10 6 1.04309 x 10“8 0.588019 0.738563 0.707107 6.09476 x 10’13 7 1.04309 x 10“8 0.588019 0.738563 0.707107 2.2745 x 10‘15 Figure 3.10 97 Example 3.11 Again consider the initial data from Example 3.9, but with h = 0.1. Hence, r 05 n = —10 1 ——9 S n S —1 0 = 0 770(71): I —1 1S n S 9 —05 n: 10 1 0 other. Again, we compute seven iterations. Unlike the previous examples, the maximum relative error of the first iteration is quite large, approximately 5.9. This is mainly due to the larger data which can be seen from the norms in Figure 3.12 and the smaller h. However, by the seventh iteration the maximum relative error is once again very small. iteration M RE norm one norm two norm three distance 1 5.90663 0.93733 0.834393 1.33463 2 2.00782 0.932334 0.838427 1.31646 0.186061 3 0.0562563 0.932087 0.838494 1.31646 0.00255853 4 0.00457199 0.932068 0.838484 1.31646 0.0000613827 5 0.000312503 0.932067 0.838483 1.31646 4.0051 x 10’6 6 9.41672 x 10‘6 0.932067 0.838483 1.31646 2.27996 x 10"7 7 2.4938 x 10‘6 0.932067 0.838483 1.31646 1.31978 x 10'8 Figure 3.12 Example 3.13 In the previous example, the initial data was larger than in previous examples. This seemed to lead to slower convergence of the iterations. Can the initial data be too large to yield convergence? With h = 1, we consider the following initial 98 data, 0 InI > 5 770(71) = n —5 S n S 5. After two iterations, it is evident that the iterations in this example are diverging in the space X), (see Figure 3.14). This suggests that the condition in Theorem 1.21 that the initial data be small is necessary. iteration M RE norm one norm two norm three distance 1 128.373 8.07398 6.46922 10.4881 2 1.44893 x 1012 850.296 1000.54 1071.73 1070.25 Figure 3.14 BIBLIOGRAPHY BIBLIOGRAPHY [1] J. Bona and R. Scott, Solutions of the Korteweg—deVries equation, in fractional order Sobolev spaces, Duke Math J. 43, 1976 87-99. [2] M. Frazier and B J awerth, Applications of the (,0 and wavelet transforms to the theory of function spaces, Wavelets and Their Applications, Ruskai et al., eds., Jones and Bartlett, Boston, 1992, 377-417. [3] I. I. Hirschman, Jr., A convexity theorem for certain groups of transformations, Journal d’Analyse, vol. 2, 1952, 209-218. [4] Y. Katznelson, An Introduction to Harmonic Analysis, Dover Publications Inc., New York, 1976. [5] C. Kenig, G. Ponce, and L. Vega, Well-posedness and scattering results for the generalized Korteweg-de Vries equation via the contraction mapping principle, Comm. in Pure and Appl. Math. 66, 1993, 527-620. [6] W. Rudin, Real and Complex Analysis, third ed., McGraw Hill, New York, 1987. [7] E. Stein, Interpolation of Linear Operators, Trans. Amer. Math. Soc. 83, 1956, 482-492. [8] E. Stein, Singular Integrals and Difierentiability Properties of Functions. Prince- ton Univ. Press, 1970. [9] E. Stein and C. Weiss, Introduction to Fourier Analysis on Euclidean Spaces. Princeton Univ. Press, 1971. [10] E. Stein, Harmonic Analysis: Real- Variable Methods, Orthogonality, and Oscil- latory Integrals, Princeton Univ. Press, Vol 42, 1993. [11] G. D. Smith, Numerical Solutions to Partial Differential Equations: Finite Dif- ference Methods, 2nd ed., Clarendon Press, Oxford, 1978. 99 IIlllllllll]llllllllllf