CHARACTERIZING THE STABILIZATION SIZE FOR SEMI-IMPLICIT FOURIER-SPECTRAL METHOD TO PHASE FIELD EQUATIONS
|
|
- Merry Robertson
- 5 years ago
- Views:
Transcription
1 CHARACTERIZING THE STABILIZATION SIZE FOR SEMI-IMPLICIT FOURIER-SPECTRAL METHOD TO PHASE FIELD EQUATIONS DONG LI, ZHONGHUA QIAO, AND TAO TANG Abstract. Recent results in the literature provide computational evidence that stabilized semi-implicit time-stepping method can efficiently simulate phase field problems involving fourth-order nonlinear diffusion, with typical examples like the Cahn-Hilliard equation and the thin film type equation. The up-to-date theoretical explanation of the numerical stability relies on the assumption that the derivative of the nonlinear potential function satisfies a Lipschitz type condition, which in a rigorous sense, implies the boundedness of the numerical solution. In this work we remove the Lipschitz assumption on the nonlinearity and prove unconditional energy stability for the stabilized semi-implicit time-stepping methods. It is shown that the size of stabilization term depends on the initial energy and the perturbation parameter but is independent of the time step. The corresponding error analysis is also established under minimal nonlinearity and regularity assumptions.. Introduction In this work we consider two phase field models: the Cahn-Hilliard equation and the molecular beam epitaxy equation MBE with slope selection. The Cahn-Hilliard CH equation was originally developed in [5] to describe phase separation in a two-component system such as metal alloy. It typically takes the form t u = u + fu, x,t Ω,,. u = u, t= where u = ux, t is a real-valued function which represents the difference between two concentrations. Due to this fact the equation. is invariant under the sign change u u. Another common form for CH is { t u = w.2 w = ǫ u + ǫ fu. As ǫ the chemical potential w tends to a limit which solves the two-phase Hele-Shaw Mullins- Sekerka problem see [2] for a heuristic derivation, [] for a convergence proof under the assumption that classical solution to the limiting Hele-Shaw problem exists. In. the spatial domain Ω is taken to be the usual 2π-periodic torus T 2 = R 2 /2πZ 2. For simplicity we only consider the periodic case but our analysis can be generalized to other settings such as bounded domain with Neumann boundary conditions. The free energy term fu is given by fu = F u = u 3 u, Fu = 4 u The parameter > is often called diffusion coefficient. Usually one is interested in the physical regime < in which the dynamics of. is close to the limiting Hele-Shaw problem after some transient time. 99 Mathematics Subject Classification. 35Q35, 65M5, 65M7. Key words and phrases. Cahn-Hilliard, energy stable, large time stepping, epitaxy, thin film.
2 2 D. LI, Z. QIAO, AND T. TANG For smooth solutions to., the total mass is conserved: d Mt, Mt = ux, tdx..4 dt In particular Mt if M =. Throughout this work we will only consider initial data u with mean zero. At the Fourier side this implies the zero th mode û =. One can then define fractional Laplacian s u for s <. The energy functional associated with. is Eu = 2 u 2 + Fu dx..5 Ω As is well known, Eq.. can be regarded as a gradient flow of Eu in H. The basic energy identity takes the form d dt Eut + t u 2 2 =..6 Note that t u has mean zero and t u is well-defined. Alternatively to avoid using, one can write.6 as d dt Eut + u + fu 2 dx =..7 Ω It follows from the energy identity that Ω Eut Eus, t s..8 This gives a priori control of H -norm of the solution. The global wellposedness of. is not an issue thanks to this fact. There is by now an extensive literature on the numerical simulation of the CH equation and related phase field models, see, e.g., [4, 7, 9,, 5, 6, 26, 3] and the references therein. On the analysis side, it is noted that Feng and Prohl [2] gave the error analysis of a semi-discrete in time and fully discrete finite element method for CH. Under a certain spectral assumption on the linearized CH operator more precisely, one has to assume the existence of classical solutions to the corresponding Hele-Shaw problem, they proved an error bound which depends on / polynomially. It is known that explicit schemes usually suffer severe time step restrictions and generally do not obey energy conservation. To enforce the energy decay property and increase the time step, a good alternative is to use implicit-explicit semi-implicit schemes in which the linear part is treated implicitly such as backward differentiation in time and the nonlinear part is evaluated explicitly. For example, in [7] Chen and Shen considered the semi-implicit Fourier-spectral scheme for. set = û n+ k û n k t = k 4 û n+ k k 2 fun k,.9 where û n denotes the Fourier coefficient of u at time step. On the other hand, the semi-implicit schemes can generate large truncation errors. As a result smaller time steps are usually required to guarantee accuracy and energy stability. To resolve this issue, a class of large time-stepping methods were proposed and analyzed in [3, 6, 26, 29, 3]. The basic idea is to add an O t stabilizing term to the numerical scheme to alleviate the time step constraint whilst keeping energy stability. The choice of the O t term is quite flexible. For example, in [3] the authors considered the Fourier spectral approximation of the modified Cahn-Hilliard-Cook equation t C = ac 2 C 3 C κ 2 C.. The explicit Fourier spectral scheme is see equation 6 therein Ĉ n+ k,t Ĉn k,t t = ik { ac 2 [ik { C + C 3 } n k + κ k 2 Ĉ n k,t] r }k..
3 PHASE-FIELD EQUATIONS 3 The time step for the above scheme has a severe constraint t κ K 4,.2 where K is the number of Fourier modes in each coordinate direction. To increase the allowed time step, the authors of [3] added a term Ak 4 Ĉn+ Ĉn to the RHS of.. Note that on the real side, this term corresponds to the fourth order dissipation, i.e. A 2 C n+ C n which roughly is of order O t. In [6], a stabilized semi-implicit scheme was considered for the CH model, with the use of an order O t stabilization term A u n+ u n. Under a condition on A of the form: { A max x Ω 2 un x un+ x + u n x 2}, 2 n,.3 one can obtain energy stability.8. Note that the condition.3 depends nonlinearity on the numerical solution. In other words, it implicitly uses the L -bound assumption on u n in order to make A a controllable constant. In [26], Shen and Yang proved energy stability of semi-implicit schemes for the Allen-Cahn and the CH equations with truncated nonlinear term. More precisely it is assumed that max f u L.4 u R which is what we referred to as the Lipschitz assumption on the nonlinearity in the abstract. The same assumption was adopted recently in [3] to analyze stabilized Crank-Nicolson or Adams-Bashforth scheme for both the CH equations. In a recent work of [4], Bertozzi et al. considered a nonlinear diffusion model of the form t u = fu u + gu u, where gu = fuφ u, and f, φ are given smooth functions. In addition f is assumed to be nonnonnegative. The numerical scheme considered in [4] takes the form u n+ u n = A 2 u n+ u n fu n u n + gu n u n,.5 t where A > is a parameter to be taken large. One should note the striking similarity between this scheme and the one introduced in [3]. In particular in both papers the biharmonic stabilization of the form A 2 u n+ u n was used. The analysis in [4] is carried out under the additional assumption that sup fu n A <..6 n This is reminiscent of the L bound on u n. Roughly speaking, all prior analytical developments are conditional in the sense that either one makes a Lipschitz assumption on the nonlinearity, or one assumes certain a priori L bounds on the numerical solution. It is very desirable to remove these technical restrictions and establish a more reasonable stability theory. Thus Problem: prove unconditional energy stability of large time-stepping semi-impliciumerical schemes for general phase field models. Here unconditional means thao restrictive assumptions should be imposed on the time step. Of course one should also develop the corresponding error analysis under minimal regularity conditions.
4 4 D. LI, Z. QIAO, AND T. TANG The purpose of this work is to settle this problem for the spectral Galerkin case. In a forthcoming work [2], we shall analyze the finite difference schemes for the CH model by using a completely different approach. We now state our main results. We first consider a stabilized semi-implicit scheme introduced in [6] following the earlier work [29]. It takes the form u n+ u n = 2 u n+ + A u n+ u n + Π N fu n, n,.7 u = Π N u. where > is the time step, and A > is the coefficient for the O regularization term. For each integer N 2, define { X N = span cosk x,sink x : k N, k Z 2}. Note that the space X N includes the constant function by taking k =. The L 2 projection operator Π N : L 2 Ω X N is defined by Π N u u,φ =, φ X N,.8 where, denotes the usual L 2 inner product on Ω. In yet other words, the operator Π N is simply the truncation of Fourier modes of L 2 functions to k N. Since π N u X N, by induction it is easy to check that u n X N for all n. Note that one can recast.7 into the usual weak formulation, for example: d t u n+,v + A u n+ u n, v + fu n, v + u n+, v =, v X N, where d t u n+ = u n+ u n /. However in our analysis it is more convenient to work with.7. Note that u n has mean zero for all n since we assume u has mean zero. Theorem. Unconditional energy stability for CH. Consider.7 with > and assume u H Ω L Ω with mean zero. Denote E = Eu the initial energy. There exists a constant β c > depending only on E such that if then where E is defined by.5. A β u 2 + log 2 +, β β c,.9 Eu n+ Eu n, n, Remark.. We stress that the above stability result works for any time step >. In particular the condition on the parameter A is independent of. In order to keep the argument simple, we do not try to optimize the dependence of A on the diffusion coefficient. This can certainly be pushed further. Remark.2. One should note that in.9, the lower bound log 2 is formally consistent with the predicted bound.3. In terms of the PDE solution ut,x, the bound.3 roughly asserts that A O ut 2. For the PDE solution, there is no L conservation and one has to trade it with the Ḣ T 2 bound with some logarithmic correction. The energy conservation gives ut H 2, and the log-correction gives log. Thus we need A log 2 from this heuristic. There is an analogue of Theorem. for the MBE equation. The MBE equation has the form t h = 2 h + g h, x,t Ω,, h = h, t=.2
5 PHASE-FIELD EQUATIONS 5 where h = hx,t : Ω R R represents the scaled height function of a thin film equation, and gz = z 2 z for z R 2. The domain Ω is again assumed to be the periodic torus T 2. Eq..2 can be regarded as an L 2 gradient flow of the energy functional Eh = 2 h G hdx,.2 where Gz = 4 z 2 2 for z R 2. Note the striking similarity between the MBE energy.2 and the CH energy.5. Roughly speaking, h is the correct scaling analogue of u in.. In fact it is well known that in D the MBE equation can be transformed into the CH equation through the change of variable u = x h. In recent [9] we obtained new upper and lower gradient bounds for the MBE equation in dimensions d 3. A refined wellposedness theory is also worked out there. Some of these results will be used in the H error analysis in this work. We refer to the introduction of [9] and also [, 2, 3, 7, 8, 27, 3] for some background material and related wellposedness/ill-posedness results. Consider the following semi-implicit scheme for MBE: h n+ h n = 2 h n+ + A h n+ h n + Π N g h n, n,.22 h = Π N h. This scheme was introduced and analyzed in [29] see also [22]. The authors of [29] first introduced the stabilized O t term of the form A h n+ h n as given in.22. They also proved that the energy stability.8 under the condition A 2 hn hn+ + h n 2, n Again, it is seen that A depends implicitly on the L bound on the numerical solution h n. The result below will provide a clean description on the size of the constant A, in the sense A is independent of the L bound on the numerical solution. Theorem.2 Unconditional energy stability for MBE. Consider.22 with >. Assume h H 2 Ω with mean zero. Assume also h L Ω. There exists a constant β c > depending only on E such that if then where E is defined by.2. A β h 2 + log 2 +, β β c, Eu n+ Eu n, n, We now state the results for error estimates. We start with the CH equation. Theorem.3 L 2 error estimate for CH. Let >. Let u H s, s 4 with mean zero. Let ut be the solution to. with initial data u. Let u n be defined according to.7 with initial data Π N u. Assume A satisfies the same condition in Theorem.. Define t m = m, m. Then ut m u m 2 A e Ctm C 2 N s +. Here C > depends only on u,, C 2 > depends on u,,s. For the MBE equation, we have the following H error estimate. Note that due to the use of H space the error bound below involves N s instead of N s. Theorem.4 H error estimate for MBE. Let > and h H s, s 5 with mean zero. Let ht be the solution to the MBE equation with initial data h. Let h n be defined according to.22 with initial data Π N h. Assume A satisfies the same condition as in Theorem.2. Define t m = m, m. Then ht m h m 2 A e Ctm C 2 N s +, Ω
6 6 D. LI, Z. QIAO, AND T. TANG where C > depends on h,, C 2 > depends on,h,s. We close this section by introducing some notation and preliminaries used in this paper. We shall use X+ to denote X + ǫ for arbitrarily small ǫ >. Similarly we can define X. We denote by T d = R d /2πZ d the 2π-periodic torus. Let Ω = T d. For any function f : Ω R, we use f L p = f Lp Ω or sometimes f p to denote the usual Lebesgue L p norm for p. If f = fx,y : Ω Ω 2 R, we shall denote by f p L x L p 2 to y denote the mixed-norm: f p L x L p 2 = fx,y p y L 2 y Ω 2 p L x Ω In a similar way one can define other mixed-norms such as f C t H etc. x m For any two quantities X and Y, we denote X Y if X CY for some constant C >. Similarly X Y if X CY for some C >. We denote X Y if X Y and Y X. The dependence of the constant C on other parameters or constants are usually clear from the context and we will often suppress this dependence. We denote X Z,,Z m Y if X CY where the constant C depends on the parameters Z,,Z m. We use the following convention for Fourier expansion on Ω = T d : Ffk = ˆfk = fxe ix k dx, fx = ˆfke ik x Ω 2π d. k Z d For f : T d R and s, we define the H s -norm and Ḣs -norm of f as f H s = + k 2s ˆfk 2 2, f Ḣs = k 2s ˆfk 2 2. k Z d k Z d provided of course the above sums are finite. Note that for s = f Ḣ = f 2. If f has mean zero, then ˆf = and in this case f H s k 2s ˆfk 2 2. k Z d For mean zero functions, we can define the fractional Laplacian s, s R via the relation s fk = k s ˆfk, k Z d. The mean zero condition is only needed for s <. Occasionally we will need to use the Littlewood Paley LP frequency projection operators. To fix the notation, let φ C c R d and satisfy φ, φ ξ = for ξ, φ ξ = for ξ 2. Let φξ := φ ξ φ 2ξ which is supported in /2 ξ 2. For any f S R d, j Z, define j fξ = φ2 j ξ ˆfξ, Ŝ j fξ = φ 2 j ξ ˆfξ, ξ R d. Let f : T d R be a smooth function. Note that f can be regarded as a tempered distribution on R d for which j f can be defined as above. For any p q, we recall the following Bernstein inequalities see [9] for a standard proof s j f Lp T d 2 js j f Lp T d, s R;.24 j f L q T d 2 jd p q f L p T d, j Z;.25 S j f Lq T d 2 jd p q f Lp T d, j
7 PHASE-FIELD EQUATIONS 7 In later sections, we will use sometimes without explicit mentioning the following interpolation inequality on T 2 : for s > and any f H s T 2 with mean zero, we have f L T 2 s f Ḣ T 2 log3 + f H s T Remark. The mean-zero condition is certainly needed in view of the f Ḣ term on the RHS. If it is replaced by f H then the inequality holds for any f noecessarily with mean zero. We include a proof of.27 for the sake of completeness. Since f has mean zero we have j f = for j < 2. Let j Z whose value will be chosen later. By using the Bernstein inequality, we have f L T 2 2 j j f L2 T j 2 js f Hs T 2, 2 j j j>j s j + 3 f Ḣ + 2 js f H s. Choosing j = const log3 + f H s then yields Proof of Stability results In this section, we will provide rigorous proofs for the stability results, i.e., Theorems. and Proof of Theorem.. Rewrite.7 as u n+ A Π N = + 2 A un A fun. 2. Lemma 2.. There is an absolute constants c > such that for any n, A + u n+ 3 c H 2 + E n +, 2.2 T2 A u n+ Ḣ T 2 + A + 3 A un 2 u n Ḣ T 2, 2.3 where E n = Eu n. Proof. In this proof for any two quantities X and Y, we shall use the notation X Y to denote X CY where C > is an absolute constant. Also for any f : T 2 R, we denote f as its average on T 2. Below one should note that fu n = u n 3 u n generally does not have mean zero even though u n has mean zero. A u n+ 3 H 2 + u n 2 + A 2 fu n fu n 2 A + u n 2 + A un 3 u n u n A + + E n +. A For u n+ Ḣ, we have This completes the proof of Lemma 2.. u n+ Ḣ u n Ḣ + A un 3 u n Ḣ + A + 3 A un 2 u n Ḣ. Lemma 2.2. For any n, E n+ E n + A u n+ u n 2 2 u n+ u n 2 2 u n un
8 8 D. LI, Z. QIAO, AND T. TANG Proof. In this proof we denote by, the usual L 2 inner product. Recall u n+ u n = 2 u n+ + A u n+ u n + Π N fu n. Taking the inner product with u n+ u n on both sides and using the identity we get b b a = 2 b 2 a 2 + b a 2, a,b R d, 2.5 u n+ u n un+ 2 2 u n u n+ u n A u n+ u n 2 2 = Π N fu n, u n+ u n. 2.6 Since all u n have Fourier modes supported in k N, we have Π N fu n, u n+ u n = fu n,u n+ u n. 2.7 By the Fundamental Theorem of Calculus, we have recall f = F Thus Fu n+ Fu n = fu n u n+ u n + = fu n u n+ u n + u n+ u n u n+ u n = fu n u n+ u n + un+ u n 2 f su n+ sds 3s 2 u n+ sds 4 3u n 2 + u n u n u n+ 2. u n+ u n E n+ E n + 2 un+ u n A + 2 un+ u n 2 2 = 4 un+ u n 2,3u n 2 + u n u n u n+ u n+ u n u n 2 + u n u n u n+ 4 u n+ u n 2 u n un Finally observe u n+ u n un+ u n u n+ u n 2 u n+ u n 2 un+ u n 2 2. The desired inequality then follows easily. Remark. By using the auxiliary function gs = Fu n + su n+ u n and the Taylor expansion g = g + g + we get g s sds, Fu n+ = Fu n + fu n u n+ u n 2 un+ u n 2 + u n+ u n 2 f u n + su n+ u n sds,
9 PHASE-FIELD EQUATIONS 9 where fz = z 3 and f z = 3z 2 for z R. From this it is easy to see that This bound will also suffice. LHS of 2.8 u n+ u n max{ un 2, u n+ 2 }. Proof of Theorem.. We inductively prove for all n, E n E, 2.9 A + u n 3 H 2 c + E +, A 2. where c > is the same absolute constant in Lemma 2.. We proceed in two steps. In Step below, we first verify that if the statement holds for some n, then it holds for n +. In Step 2, we check the base case, namely for n = the statement holds. We organize our whole argument in this reverse order rather than checking the base case n = first and then perform induction because the verification for the base case n = can be viewed as more or less a special case of the proof in Step. Step : the induction step n n +. Assume the induction holds for some n. We now verify the statement for n +. By Lemma 2., we have A + u n+ 3 H 2 c + A + E n + c + E +. A A Thus we only need to check E n+ E. In fact we shall show E n+ E n. By Lemma 2.2, we only need to show the inequality A un un We shall use the log-interpolation inequality see.27 and choose s = 3 2 for any f with mean zero: f L T 2 d f Ḣ T 2 log f 3 H 2, T 2 where d > is an absolute constant. In the rest of this proof, to ease the notation we shall use X E Y to denote X C E Y where C E is a constant depending only on E. Clearly u n d u n Ḣ log d u n 3 H A + 2E log 3 + c + E + A E 2 + log A + log }{{} + 2 log. 2.3 =:m Here in the above inequality, if then it is not difficult to check that the log term is not present. In the rest of this proof we shall just assume < without loss of generality. The case is similar and even easier. Now u n 2 E m 2 + log 2.
10 D. LI, Z. QIAO, AND T. TANG By 2.2 and Lemma 2., we have below in the third inequality we drop /A since A u n+ u n+ Ḣ log u n+ 3 H Therefore + A + un 2 A un Ḣ log u n+ 3 H un 2 A un Ḣ log E + m2 + log 2 A u n+ 3 H m + 2 log E m + 2 log + m log 3 A E m + m3 A log u n 2 + u n+ 2 E m + m log 6. A Therefore to show the inequality 2., it suffices to prove A + C E m + m log 6, 2.5 A where m = 2 + log A + log. Now we discuss two cases. Case : C E 3 log 6. In this case we choose A such that A E m 2 = + log A + log 2. Clearly for, we juseed to choose A E. On the other hand, for <, it suffices to take A = β log 2, with β sufficiently large depending only on E. Thus in both cases if we take A = β max{ log 2,}, with β E, then 2.5 holds. Case 2: C E 3 log 6. In this case we have log E + log. In this case we will not prove 2.5 but prove 2. directly. We first go back to the bound on u n. Easy to check that u n E m, u n+ E + m2 m. A The needed inequality on A then takes the form 2 A C E + m + m3. A Again we only need to choose A such that A E m 2. The same choice of A as in Case with β larger if necessary works. Concluding from both cases, we have proved the inequality 2. holds. This completes the induction step for n n +.
11 PHASE-FIELD EQUATIONS Step 2: verification of the base step n =. By Lemma 2. we have A + u 3 H 2 c + E +. A Therefore we only need to check E E. This amounts to checking the inequality A Π Nu u 2. Now clearly the estimate of u proceeds in the same way as in Step. On the other hand since we juseed to choose A such that This completes the proof of Theorem.. Π N u u, A E u 2 + log Proof of Theorem.2. This is similar to the proof of Theorem.. Therefore we only sketch the needed modifications. In terms of scaling it is useful to think of h n as u n in Theorem.. Write.22 as h n+ A Π N = + 2 A hn A g h n. 2.6 In place of Lemma 2. we have the following lemma. We omit the proof since it is quite similar. Lemma 2.3. There is an absolute constants c > such that Here E n = Eh n. Lemma 2.4. For any n, h n+ H 5 2 T2 c A + + A E n +, h n+ Ḣ2 T 2 + A + 3 A hn 2 h n Ḣ2 T 2. E n+ E n + A hn+ h n 2 2 h n+ h n max{ hn 2, h n+ 2 }. 2.7 Proof. Taking inner product with h n+ h n on both sides of.22, we get hn+ h n hn+ 2 2 h n h n+ h n A h n+ h n 2 2 = g h n, h n+ h n. Recall gz = z 2 z = G and Gz = 4 z 2 2. Introduce By using the expansion we get Hs = G h n + s h n+ h n. H = H + H + G h n+ G h n = g h n h n+ h n 2 + i h n+ h n j h n+ h n i,j= H s sds, ij G h n + s h n+ h n sds,
12 2 D. LI, Z. QIAO, AND T. TANG Now denote Gz = 4 z 4. Then = E n+ E n + hn+ h n hn+ h n A + 2 hn+ h n i h n+ h n j h n+ h n ij G h n + s h n+ h n sds,, i,j= where represents the constant function with value. Now since ij Gz = z 2 δ ij + 2z j z i, we have the point-wise bound ij Gz 3 z 2. Thus ij G h n + s h n+ h n 3max{ h n 2, h n+ 2 }. The desired inequality now follows from this and the interpolation inequality h 2 h 2 2 h This completes the proof of the lemma. Proof of Theorem.2. We only need to check the induction hypothesis E n E, A + h n 5 H 2 c + E +, A for n +. Here c > is the same absolute constant in Lemma 2.3. By Lemma 2.3, we have A + h n+ 5 H 2 c + A + E n + c + E +. A A Thus we only need to check E n+ E n. By Lemma 2.4, this amounts to proving the inequality We shall again use the inequality A max{ hn 2, h n+ 2 }. 2.9 f L T 2 d f Ḣ T 2 log f H 3 2 T2 + 3, 2.2 where d > is an absolute constant, and f has mean zero. Clearly h n d h n Ḣ2 log h n 5 H E A + d log 3 + c + E A The rest of the argumenow is similar to that in the Proof of Theorem.. We omit further repetitive details. 3. Bounds on the PDE solution of CH Consider t w = 2 w + fw, w = w. t= 3. Recall that the corresponding energy E is defined by.5.
13 PHASE-FIELD EQUATIONS 3 Proposition 3.. Let <. Assume the initial data w H 2 T 2 with mean zero. Assume w. Then E sup wt t< log + log E +, 3.2 where E = Ew. Proof. First consider the regime < t. Write Then wt = e t 2 w + t t e t s 2 fwsds. wt w + 2 t s 2 fws ds w + 2 t 2 w 3 L s,x [,t] + w L s,x [,t]. 3.3 By using a continuity argument on the quantity w L s,x [,t], we get sup wt, 3.4 t ǫ where ǫ > is a sufficiently small absolute constant. Strictly speaking the value of ǫ depends on the implied constants hidden in the inequalities w and <. Next we consider the L bound in the time regime t ǫ. First observe that by using energy conservation, we have E wt 2. Set t = t 2 ǫ. Then We bound the Ḣ+ -norm of w as t wt = e t t 2 wt + e t s 2 fwsds. t wt Ḣ+ + e t t 2 wt e t s ws 3 ws 2 ds t t t t wt Ḣ + t s 3 4 ws 3 + ws Ḣds Ḣ t 2 E E Then recall that w has mean zero 2 + E t 3 E 2 + E This completes the proof of Proposition 3.. wt wt 2 log + wt H + E + log + log E. Remark. By using the method in [9], one can prove a wellposedness result for w L 2 T 2. However we shall noeed this refinement here.
14 4 D. LI, Z. QIAO, AND T. TANG Proposition 3.2. Assume the initial data w has mean zero and w H s T 2, s 4. Then T t w 2 2dt,w + T. 3.6 Proof. To ease the notation we shall write,w as throughout this proof. By using the smoothing effect, it is easy to show that From energy conservation, we have sup t w H. 3.7 t< Thus by interpolation, This implies Now we only need to show Observe T t w 2 2dt. 3.8 t w 4 2dt. 3.9 t w 2 2dt + T. 3. t w 2 2dt. 3. t w = 3 w + 2 fw. Multiplying both sides by t w and integrating by parts, we get t w 2 2 = d 2 dt 2 w fw t wdx. Thus Ω d 2 dt 2 w 2 2 t w fw 2 t w 2 2 t w const w 3 H + w 4 H From this and standard H 4 global wellposedness theory, we get The desired inequality then follows. t w 2 2dt Error estimate for CH In this section we give the estimate for CH in L 2.
15 PHASE-FIELD EQUATIONS Auxiliary L 2 error estimate for near solutions. Consider v n+ v n = 2 v n+ + A v n+ v n + Π N fv n + G n, n, ṽ n+ ṽ n = 2 ṽ n+ + A ṽ n+ ṽ n + Π N fṽ n + G n 2, n, v = v, ṽ = ṽ, where v and ṽ has mean zero. Denote G n = G n G n 2. We first state and prove a simple lemma. Lemma 4. Discrete Gronwall inequality. Let > and y n, α n, β n for n =,,2,. Suppose Then for any m, we have In particular Proof. Clearly y n+ y n m y m exp n= α n y n + β n, n. α n y + m y m exp n= m k= α n y + exp m k= m j=k+ y n+ + α n y n + β n e αn y n + β n, n. 4. α j βk. 4.2 β k. 4.3 Thus exp n j= Summing n from to m, we get Thus 4.2 is obtained. n α j y n+ exp m exp j= j= α j y m y + α j y n + exp m n= exp n α j βn. j= n α j βn. j= Proposition 4.. For solutions of 4., assume for some N >, N 2 >, Then for any m, sup ṽ n N, n sup v n 2 N 2, n sup ṽ n 2 N n v m ṽ m 2 2 exp m C + N 4 + N 4 2 where C > is an absolute constant. v ṽ A v ṽ m Remark. The same proposition holds if Π N is replaced by the identity operator. n= G n 2 2, 4.5
16 6 D. LI, Z. QIAO, AND T. TANG Proof of Proposition 4.. Denote e n = v n ṽ n. Then e n+ e n = 2 e n+ + A e n+ e n + Π N fv n fṽ n + G n. 4.6 Taking L 2 -inner product with e n+ on both sides, we get 2 en+ 2 2 e n e n+ e n e n A 2 en+ 2 2 e n e n+ e n 2 2 Obviously = G n, e n+ + fv n fṽ n, Π N e n G n, e n+ 2 G n 2 2 On the other hand, recalling f z = 3z 2, we get fv n fṽ n = f ṽ n + se n dse n + 8 en = a + a 2 ṽ n 2 e n + a 3 ṽ n e n 2 + a 4 e n 3, where a i, i =,,4 are constants which can be computed explicitly. We now estimate the contribution of each term. In the rest of this proof, to ease the notation, we shall denote by C an absolute constant whose value may change from line to line. Clearly a + a 2 ṽ n 2 e n, e n+ C + ṽ n 2 e n 2 e n+ C + N 4 2 e n en By using the interpolation inequality f 4 f 2 2 f 2 2, we get Similarly a 4 e n 3, e n+ a 3 ṽ n e n 2, e n+ C ṽ n e n 2 4 e n+ 2 C N2 e n en+ 2 2 C N2 en 2 2 e n en+ 2 2 C N2 N2 2 e n en C en en+ 2 2 C en 2 2 e n en+ 2 2 C N4 2 en en Collecting the estimates, we get e n+ 2 2 e n A e n+ 2 2 e n 2 2 Define 4 G n C + N4 + N 4 2 e n y n = e n A e n 2 2, α = C + N4 + N2 4, βn = 4 G n 2 2. Then obviously y n+ y n αy n + β n. The desired result then follows from Lemma 4..
17 PHASE-FIELD EQUATIONS L 2 error estimate for CH Proof of Theorem.3. In this proof to ease the notation, we shall denote by C a constant depending only on,u. The value of C may vary from line to line. For any two quantities X and Y, we shall write X Y if X CY. Note that we shall still keep track of the dependence on the parameter A and also the regularity index s. We need to compare u n+ u n = 2 u n+ + A u n+ u n + Π N fu n, t u = 2 u + fu, 4.2 ũ = Π N u, u = u. We first rewrite the PDE solution u in the discretized form. Note that for one-variable function h = ht, we have the formulae htdt = h + htdt = h+ + h t + tdt, 4.3 h t tdt. 4.4 By using the above formulae and integrating the PDE for u on the time interval [,+ ], we get u+ u = 2 u+ + A u+ u + Π N fu + Π >N fu + G n, 4.5 where Π >N = Id Π N Id is the identity operator and Now G n = t u tdt + t fu + tdt A t udt. 4.6 G tn+ n 2 t u 2 dt + t u 2 dt f u + A. 4.7 By Proposition 3., we have u. Since t u 2 t u 2, we get Therefore by Proposition 3.2, m n= It is easy to check that which implies This gives G n 2 + A G n A 2 m n= t u 2 dt + A t u 2 2dt tm 2. t u 2 2dt + A 2 + t m. sup ut H s s t sup fut H s s. 4.8 t Π >N fu 2 2 s N 2s t m /. 4.9
18 8 D. LI, Z. QIAO, AND T. TANG Therefore Note that m G n Π >N fu 2 2 s + t m 2 + N 2s + A n= u Π N u 2 s N s, u Π N u 2 s N s. By Proposition 4., we then get u m ut m 2 2 s + A 2 e N Ctm 2s + N 2s + + t m 2 + N 2s. Since by assumption we have s 4, clearly by Cauchy-Schwartz This implies N 2s 2 + N 4s 2 + N 2s. u m ut m 2 s + Ae Ctm N s +. Remark 4.. From the above analysis, it is clear that our regularity assumption H s, s 4 on the initial data comes from bounding the term t u 2 2dt which in turn arose from rewriting the diffusion term 2 u into the time-discretized form. Recall t u = 2 u + fu. For < t, the linear effect is dominant and one can roughly regard t u 2 P <t 4 u, where P <t 4 is the Littlewood-Paley projection to the frequency regime ξ t 4. Heuristically speaking t u 2 2 t 2 P<t 4 2 u 2 2 t P <t 4 2 u 2 2 which is barely non-integrable in t, provided we assume H 4 regularity on u. Of course a well-known technique in these situations is to use the maximal regularity estimates of the linear semi-group to get integrability in t. In the L 2 case the usual energy estimate suffices and this is why we need H 4 regularity on the initial data. 5. Error estimate for MBE 5.. Auxiliary H estimate for MBE. For MBE we need to consider q n+ q n = 2 q n+ + A q n+ q n + Π N g q n + G n, q n+ q n = 2 q n+ + A q n+ q n + Π N g q n + G n 2, q = q, q = q, where we recall gz = z 2 z for z R 2. As before q and q are assumed to have mean zero. Denote G n = G n G n 2. Proposition 5.. Assume for some N > q n + q n 2 + q n 2 N. 5.2 Then for any m, q n q n 2 2 e sup n C +N 4 m where C > is an absolute constant. q q A q q m n= G n , 5.3
19 PHASE-FIELD EQUATIONS 9 Proof. Denote e n = q n q n. Then Then Thus e n+ e n = 2 e n+ + A e n+ e n + Π N g q n g q n G n. Taking L 2 -inner product with e n+ on both sides, we get 2 en+ 2 2 e n e n+ e n e n A 2 en+ 2 2 e n e n+ e n 2 2 = G n, e n+ +g q n g q n, Π N e n }{{}}{{} I I 2 For the first term on the RHS of 5.4, we simply bound it as For the second term I 2, recalling gz = z 2 z, we have We then obtain I G n en g q n g q n = O e n + O q n 2 e n + O q n e n 2 + O e n 3. g q n g q n 2 + N 2 e n 2 + N e n e n N 2 e n 2 + N e n 2 e n 2 + e n 2 e n N 2 e n I 2 C + N4 e n+ 2 2 e n 2 2 C + N4 The desired result then follows from Lemma 4.. e n en A e n+ 2 2 e n 2 2 e n G n Proof of Theorem.4. Similar to the proof of Theorem.3, we need to compare h n+ h n = 2 h n+ + A h n+ h n + Π N g h n, t h = 2 h + g h, h = Π N h, h = h. On the time interval [,+ ], we have where h+ h = 2 h+ + A h+ h + Π N g h + Π >N g h + G n, 5.9 G n = t h tdt + t g ht + tdt A t hdt.
20 2 D. LI, Z. QIAO, AND T. TANG Now we only need to verify the estimates: Recall T T t h 2 2dt,h + T, 5. t g h 2 2dt,h + T. 5. t h = 2 h + g h. Multiplying both sides by 3 t h and integrating by parts, we get t h 2 2 = d 2 dt 2 h g h t hdx, and d 2 dt 2 h 2 2 t h g h 2 t h 2 2 t h const h 3 H + h 5 H This together with standard local wellposedness theory, cf. [9] for more refined results yields The smoothing effect gives control for t. Thus T For the term t g h 2, we note that Thus Finally we get The theorem is proved. t g h 2 t h 2 2dt,h. t h 2 2dt,h + T. t h 2 h h 2 h 2 + t h 2,h t h T t g h 2 2dt,h + T. 5.4 ht m h m 2 + Ae Ctm N s Concluding remarks In this work we considered a class of large time-stepping methods for the phase field models such as the Cahn-Hilliard equation and the thin film equation with fourth order dissipation. We analyzed the representative case see.7 and.22 which is first order in time and Fourier-spectral in space, with a stabilization O t term of the form A u n+ u n. 6. For A sufficiently large A O log 2, we proved unconditional energy stability independent of the time step. The corresponding error analysis is also carried out in full detail L 2 for CH and H for MBE. It is worth emphasizing that our analysis does not require any additional Lipschitz assumption on the nonlinearity, or any a priori bounds on the numerical solution. It is expected our theoretical framework can be extended in several directions. We discuss a few such possibilities below the fold.
21 PHASE-FIELD EQUATIONS 2 General stabilization techniques. There are a myriad of ways of introducing the stabilization term. Taking the first order in time methods as an example, instead of 6., one can consider a more general form ABu n+ u n, 6.2 where B is a general operator. One example is B = 2 which is already used in the aforementioned works [3, 4]. Similarly one can consider B = s s > is real or even a general pseudo-differential operator. It will be interesting to carry out a comparative study of these different stabilization techniques and identify the corresponding stability regions. Another issue is to investigate the lower bound on the parameter A. In typical numerical simulations the stability is observed to hold for relatively small values of A the threshold value exhibits a weak dependence on the time step and the diffusion coefficient, cf. the numerical simulation results in [6]. This certainly merits further study and probably one has to fine-tune our analysis with some numerically verifiable bounds. Higher order time stepping methods. In [29], Xu and Tang considered a second order scheme for MBE: 3h n+ 4h n + h n + 2 h n+ 2 = A h n+ 2h n + h n + Π N g 2h n h n, n, 6.3 where h is the initial condition and h is computed by the first order scheme.22. Here to keep some consistency with our setup we have added the projection operator Π N in front of the nonlinear term. This scheme is called BD2/EP2 since it is obtained by combining a second-order backward differentiation BD2 for the time derivative term and a second-order extrapolation EP2 for the explicit treatment of the nonlinear term. A similar higher order BD3/EP3 scheme is also presented in [29]. The stability analysis in [29] is conditional in the sense that the choice of A depends on the a priori gradient bound on the numerical solution. Moreover, quite different from the first order in time methods, the energy stability for higher order methods typically takes the form Eh n Eh + O, n T, 6.4 where the implied constant in the O term usually depends on the time interval [, T]. In yet other words one cannot achieve strict monotonic decay of energy as in the first order case. A very natural problem is to extend our analysis to cover these cases. By using our analysis it is also possible to refine the stability results in [26] and remove the Lipschitz assumption on the nonlinearity in the case of second order implicit scheme. For second order semi-implicit schemes it is expected that our method can be extended to prove an unconditional stability result at least for time steps which are moderately small. We plan to address these issues in a future publication. General phase field models possibly with higher order dissipations. In [9], the authors considered the sixth order scalar model t u = ǫ 2 W u + ǫ 2 ηǫ 2 u W u, 6.5 where Wu = 4 u2 2 and η > is a given constant. This equation arises in the modeling of pore formation in functionalized polymers [4]. The numerical experiments in [9] used implicit time stepping together with Newton s method at each time step. From our point of view it will be interesting to use the numerical schemes similar to.7 and establish the corresponding stability and error convergence results. In a similar vein one can also consider the volume-preserving vector CH model in the same paper see equation 7 in [9] and also the nonlinear diffusion model in [4]. Yet another possibility is to study the model with general fractional dissipation which is already mentioned in the introduction of [9]. Also one can extend our analysis to the phase-fields models
22 22 D. LI, Z. QIAO, AND T. TANG of two-phase complex fluids see [28] for a pioneering study in this direction. In any case a first step in the analysis is to establish similar results as in [9]. The above list is certainly not exhaustive. For example we did not include the analysis of the Allen-Cahn model which will be quite similar to the CH case from our point of view. To keep the presentation simple we leave out the case of dimensions d = and d = 3 which can be similarly handled. One can also consider generalizing the analysis herein to finite difference schemes and even some hybrid schemes. In [2] we will introduce a completely new approach to tackle some of these problems. Another direction is to consider the phase field models with stochastic noises. One can introduce similar numerical stabilization techniques as in the deterministic case and prove stability and convergence in these settings. We plan to investigate these problems in the future. Acknowledgments. D. Li was supported by an Nserc discovery grant. The research of Z. Qiao is partially supported by the Hong Kong Research Council GRF grants 222, and NSFC/RGC Joint Research Scheme N HKBU24/2. The research of T. Tang is mainly supported by Hong Kong Research Council GRF Grants and Hong Kong Baptist University FRG grants. References [] N.D. Alikakos, P.W. Bates and X.F. Chen. Convergence of the Cahn-Hilliard equation to the Hele-Shaw model. Arch. Ration. Mech. Anal , [2] J. Bourgain and D. Li. Strong ill-posedness of the incompressible Euler equation in borderline Sobolev spaces. To appear in Invent. Math. [3] J. Bourgain and D. Li. Strong illposedness of the incompressible Euler equation in integer C m spaces. To appear in Geom. Funct. Anal. [4] A. Bertozzi, N. Ju and H. Lu. A biharmonic-modified forward time stepping method for fourth order nonlinear diffusion equations. Disc. Conti. Dyn. Sys. 29 2, no. 4, [5] J.W. Cahn, J.E. Hilliard. Free energy of a nonuniform system. I. Interfacial energy free energy, J. Chem. Phys [6] H.D. Ceniceros, R. L. Nos and A.M. Roma. Three-dimensional, fully adaptive simulations of phase-field fluid models. J. Comput. Phys., 229 2, pp [7] L.Q. Chen, J. Shen. Applications of semi-implicit Fourier-spectral method to phase field equations. Comput. Phys. Comm., 8 998, pp [8] F. Chen and J. Shen. Efficient energy stable schemes with spectral discretization in space for anisotropic Cahn-Hilliard systems. Commun. Comput. Phys., 3 23, [9] A. Christlieb, J. Jones, K. Promislow, B. Wetton, M. Willoughby. High accuracy solutions to energy gradient flows from material science models. J. Comput. Phys , part A, [] W. M. Feng, P. Yu, S. Y. Hu, Z. K. Liu, Q. Du and L. Q. Chen A Fourier spectral moving mesh method for the Cahn-Hilliard equation with elasticity. Commun. Comput. Phys., 5 29, pp [] G. Ehrlich and F.G. Hudda. Atomic view of surface diffusion: tungsten on tungsten. J. Chem. Phys , 36. [2] X. B. Feng and A. Prohl. Error analysis of a mixed finite element method for the Cahn-Hilliard equation. Numer. Math., 99 24, pp [3] X. Feng, T. Tang and J. Yang. Stabilized Crank-Nicolson/Adams-Bashforth schemes for phase field models. East Asian J. Appl. Math. 3 23, no., [4] N. Gavish, J. Jones, Z. Xu, A. Christlieb, K. Promislow. Variational models of network formation and ion transport: applications to perfluorosulfonate ionomer membranes. Polymers 4 22, [5] H. Gomez and T.J.R. Hughes. Provably unconditionally stable, second-order time-accurate, mixed variational methods for phase-field models. J. Comput. Phys., 23 2, pp [6] Y. He, Y. Liu and T. Tang. On large time-stepping methods for the Cahn-Hilliard equation. Appl. Numer. Math., 57 27, [7] B. Li and J.G. Liu. Thin film epitaxy with or without slope selection. Euro. Jnl of Appl. Math., 4 23, pp [8] D. Li. On a frequency localized Bernstein inequality and some generalized Poincaré-type inequalities. Math. Res. Lett. 2 23, no. 5, [9] D. Li, Z. Qiao and T. Tang. Gradient bounds for a thin film epitaxy equation. Submitted, 24. [2] D. Li, Z. Qiao and T. Tang. In preparation. [2] R.L. Pego. Front migration in the nonlinear Cahn-Hilliard equation. Proc. Roy. Soc. London A
23 PHASE-FIELD EQUATIONS 23 [22] Z. Qiao, Z. Zhang and T. Tang. An adaptive time-stepping strategy for the molecular beam epitaxy models. SIAM J. Sci. Comput. 33 2, no. 3, [23] P. Rybka and K. Hoffmann. Convergence of solutions to Cahn-Hilliard equation. Comm. Partial Differential Equations , no. 5 6, [24] R.L. Schwoebel and E.J. Shipsey. Step motion on crystal surfaces. J. Appl. Phys , [25] R.L. Schwoebel. Step motion on crystal surfaces II. J. Appl. Phys , 64. [26] J. Shen and X. Yang. Numerical approximations of Allen-Cahn and Cahn-Hilliard equations. Discrete Contin. Dyn. Syst. A, 28 2, [27] J. Shen, C. Wang, X. Wang, S.M. Wise. Second-order convex splitting schemes for gradient flows with Ehrlich-Schwoebel type energy: application to thin film epitaxy. SIAM J. Numer. Anal. 5 22, no., [28] J. Shen and X. Yang. Decoupled energy stable schemes for phase-field models of two-phase complex fluids. SIAM J. Sci. Comput , no., B22 B45. [29] C. Xu and T. Tang. Stability analysis of large time-stepping methods for epitaxial growth models. SIAM J. Numer. Anal , no. 4, [3] J. Zhu, L.-Q. Chen, J. Shen, and V. Tikare. Coarsening kinetics from a variable-mobility Cahn-Hilliard equation: Application of a semi-implicit Fourier spectral method, Phys. Rev. E 3, 6 999, pp [3] C. Wang, S. Wang, and S.M. Wise. Unconditionally stable schemes for equations of thin film epitaxy. Disc. Contin. Dyn. Sys. Ser. A, 28 2, pp D. Li Department of Mathematics, University of British Columbia, 984 Mathematics Road, Vancouver, BC, Canada V6TZ2 address: dli@math.ubc.ca Z. Qiao Department of Applied Mathematics, The Hong Kong Polytechnic University, Hung Hom, Hong Kong address: zhonghua.qiao@polyu.edu.hk T. Tang Department of Mathematics & Institute for Computational and Theoretical Studies, Hong Kong Baptist University, Kowloon, Hong Kong address: ttang@math.hkbu.edu.hk
Numerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
More informationERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION. Xiaofeng Yang. (Communicated by Jie Shen)
DISCRETE AND CONTINUOUS doi:1.3934/dcdsb.29.11.157 DYNAMICAL SYSTEMS SERIES B Volume 11, Number 4, June 29 pp. 157 17 ERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION Xiaofeng Yang
More informationHIGHLY EFFICIENT AND ACCURATE NUMERICAL SCHEMES FOR THE EPITAXIAL THIN FILM GROWTH MODELS BY USING THE SAV APPROACH
HIGHLY EFFICIENT AND ACCURATE NUMERICAL SCHEMES FOR THE EPITAXIAL THIN FILM GROWTH MODELS BY USING THE SAV APPROACH QING CHENG, JIE SHEN, XIAOFENG YANG Abstract. We develop in this paper highly efficient,
More informationNUMERICAL APPROXIMATIONS OF ALLEN-CAHN AND CAHN-HILLIARD EQUATIONS
UMERICAL APPROXIMATIOS OF ALLE-CAH AD CAH-HILLIARD EQUATIOS JIE SHE AD XIAOFEG YAG Abstract. Stability analyses and error estimates are carried out for a number of commonly used numerical schemes for the
More informationEXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN EQUATION WITH DEGENERATE MOBILITY
Electronic Journal of Differential Equations, Vol. 216 216), No. 329, pp. 1 22. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN
More informationwith deterministic and noise terms for a general non-homogeneous Cahn-Hilliard equation Modeling and Asymptotics
12-3-2009 Modeling and Asymptotics for a general non-homogeneous Cahn-Hilliard equation with deterministic and noise terms D.C. Antonopoulou (Joint with G. Karali and G. Kossioris) Department of Applied
More informationOn the local existence for an active scalar equation in critical regularity setting
On the local existence for an active scalar equation in critical regularity setting Walter Rusin Department of Mathematics, Oklahoma State University, Stillwater, OK 7478 Fei Wang Department of Mathematics,
More informationA PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM
A PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM Hyun Geun LEE 1, Jeong-Whan CHOI 1 and Junseok KIM 1 1) Department of Mathematics, Korea University, Seoul
More informationNUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART I: ERROR ANALYSIS UNDER MINIMUM REGULARITIES
NUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM PART I: ERROR ANALYSIS UNDER MINIMUM REGULARITIES XIAOBING FENG AND ANDREAS PROHL Abstract. In this first part
More informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
More informationNon-radial solutions to a bi-harmonic equation with negative exponent
Non-radial solutions to a bi-harmonic equation with negative exponent Ali Hyder Department of Mathematics, University of British Columbia, Vancouver BC V6TZ2, Canada ali.hyder@math.ubc.ca Juncheng Wei
More informationDissipative quasi-geostrophic equations with L p data
Electronic Journal of Differential Equations, Vol. (), No. 56, pp. 3. ISSN: 7-669. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) Dissipative quasi-geostrophic
More informationResearch Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations
Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationThe 2D Magnetohydrodynamic Equations with Partial Dissipation. Oklahoma State University
The 2D Magnetohydrodynamic Equations with Partial Dissipation Jiahong Wu Oklahoma State University IPAM Workshop Mathematical Analysis of Turbulence IPAM, UCLA, September 29-October 3, 2014 1 / 112 Outline
More informationEnergy transfer model and large periodic boundary value problem for the quintic NLS
Energy transfer model and large periodic boundary value problem for the quintic NS Hideo Takaoka Department of Mathematics, Kobe University 1 ntroduction This note is based on a talk given at the conference
More informationDynamics and stochastic limit for phase transition problems. problems with noise. Dimitra Antonopoulou IACM-FORTH
Dynamics and stochastic limit for phase transition problems with noise IACM-FORTH International Conference on Applied Mathematics, September 16-20, 2013, ACMAC, Heraklion, Greece Table of contents 1 1.1
More informationOSGOOD TYPE REGULARITY CRITERION FOR THE 3D NEWTON-BOUSSINESQ EQUATION
Electronic Journal of Differential Equations, Vol. 013 (013), No. 3, pp. 1 6. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu OSGOOD TYPE REGULARITY
More informationOPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES
OPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES RENJUN DUAN Department of Mathematics, City University of Hong Kong 83 Tat Chee Avenue, Kowloon, Hong Kong,
More informationSECOND-ORDER FULLY DISCRETIZED PROJECTION METHOD FOR INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
Tenth MSU Conference on Differential Equations and Computational Simulations. Electronic Journal of Differential Equations, Conference 3 (06), pp. 9 0. ISSN: 07-669. URL: http://ejde.math.txstate.edu or
More informationA Second Order Energy Stable Scheme for the Cahn-Hilliard-Hele-Shaw Equations
A Second Order Energy Stable Scheme for the Cahn-Hilliard-Hele-Shaw Equations Wenbin Chen Wenqiang Feng Yuan Liu Cheng Wang Steven M. Wise August 9, 07 Abstract We present a second-order-in-time finite
More informationA LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE. 1.
A LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE THOMAS CHEN AND NATAŠA PAVLOVIĆ Abstract. We prove a Beale-Kato-Majda criterion
More informationStep Bunching in Epitaxial Growth with Elasticity Effects
Step Bunching in Epitaxial Growth with Elasticity Effects Tao Luo Department of Mathematics The Hong Kong University of Science and Technology joint work with Yang Xiang, Aaron Yip 05 Jan 2017 Tao Luo
More informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
More informationarxiv: v1 [math.na] 15 Aug 2007
A POSTERIORI ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS OF THE CAHN-HILLIARD EQUATION AND THE HELE-SHAW FLOW XIAOBING FENG AND HAIJUN WU arxiv:78.116v1 [math.na] 15 Aug 7 Abstract. This paper develops
More informationA COUNTEREXAMPLE TO AN ENDPOINT BILINEAR STRICHARTZ INEQUALITY TERENCE TAO. t L x (R R2 ) f L 2 x (R2 )
Electronic Journal of Differential Equations, Vol. 2006(2006), No. 5, pp. 6. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) A COUNTEREXAMPLE
More informationStability and Convergence Analysis of Fully Discrete Fourier Collocation Spectral Method for 3-D Viscous Burgers Equation
J Sci Comput (01) 53:10 18 DOI 10.1007/s10915-01-961-8 Stability and Convergence Analysis of Fully Discrete Fourier Collocation Spectral Method for 3-D Viscous Burgers Equation Sigal Gottlieb Cheng Wang
More informationON THE MAXIMUM PRINCIPLE PRESERVING SCHEMES FOR THE GENERALIZED ALLEN CAHN EQUATION
COMMUN. MATH. SCI. Vol., No., pp. 7 c International Press ON THE MAXIMUM PRINCIPLE PRESERVING SCHEMES FOR THE GENERALIZED ALLEN CAHN EQUATION JIE SHEN, TAO TANG, AND JIANG YANG Abstract. This paper is
More informationWeak solutions for the Cahn-Hilliard equation with degenerate mobility
Archive for Rational Mechanics and Analysis manuscript No. (will be inserted by the editor) Shibin Dai Qiang Du Weak solutions for the Cahn-Hilliard equation with degenerate mobility Abstract In this paper,
More informationA generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem
A generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem Dave McCormick joint work with James Robinson and José Rodrigo Mathematics and Statistics Centre for Doctoral Training University
More informationGlobal well-posedness of the primitive equations of oceanic and atmospheric dynamics
Global well-posedness of the primitive equations of oceanic and atmospheric dynamics Jinkai Li Department of Mathematics The Chinese University of Hong Kong Dynamics of Small Scales in Fluids ICERM, Feb
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationPart IV: Numerical schemes for the phase-filed model
Part IV: Numerical schemes for the phase-filed model Jie Shen Department of Mathematics Purdue University IMS, Singapore July 29-3, 29 The complete set of governing equations Find u, p, (φ, ξ) such that
More informationStabilization for the Wave Equation with Variable Coefficients and Balakrishnan-Taylor Damping. Tae Gab Ha
TAIWANEE JOURNAL OF MATHEMATIC Vol. xx, No. x, pp., xx 0xx DOI: 0.650/tjm/788 This paper is available online at http://journal.tms.org.tw tabilization for the Wave Equation with Variable Coefficients and
More informationON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction
Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist
More informationNumerical Simulations on Two Nonlinear Biharmonic Evolution Equations
Numerical Simulations on Two Nonlinear Biharmonic Evolution Equations Ming-Jun Lai, Chun Liu, and Paul Wenston Abstract We numerically simulate the following two nonlinear evolution equations with a fourth
More informationON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM
ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM OLEG ZUBELEVICH DEPARTMENT OF MATHEMATICS THE BUDGET AND TREASURY ACADEMY OF THE MINISTRY OF FINANCE OF THE RUSSIAN FEDERATION 7, ZLATOUSTINSKY MALIY PER.,
More informationNUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART II: ERROR ANALYSIS AND CONVERGENCE OF THE INTERFACE
NUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART II: ERROR ANALYSIS AND CONVERGENCE OF THE INTERFACE XIAOBING FENG AND ANDREAS PROHL Abstract. In this
More informationGLOBAL EXISTENCE AND ENERGY DECAY OF SOLUTIONS TO A PETROVSKY EQUATION WITH GENERAL NONLINEAR DISSIPATION AND SOURCE TERM
Georgian Mathematical Journal Volume 3 (26), Number 3, 397 4 GLOBAL EXITENCE AND ENERGY DECAY OF OLUTION TO A PETROVKY EQUATION WITH GENERAL NONLINEAR DIIPATION AND OURCE TERM NOUR-EDDINE AMROUN AND ABBE
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationModels for dynamic fracture based on Griffith s criterion
Models for dynamic fracture based on Griffith s criterion Christopher J. Larsen Abstract There has been much recent progress in extending Griffith s criterion for crack growth into mathematical models
More informationExistence of minimizers for the pure displacement problem in nonlinear elasticity
Existence of minimizers for the pure displacement problem in nonlinear elasticity Cristinel Mardare Université Pierre et Marie Curie - Paris 6, Laboratoire Jacques-Louis Lions, Paris, F-75005 France Abstract
More informationPARTIAL DIFFERENTIAL EQUATIONS. Lecturer: D.M.A. Stuart MT 2007
PARTIAL DIFFERENTIAL EQUATIONS Lecturer: D.M.A. Stuart MT 2007 In addition to the sets of lecture notes written by previous lecturers ([1, 2]) the books [4, 7] are very good for the PDE topics in the course.
More informationLOCAL WELL-POSEDNESS FOR AN ERICKSEN-LESLIE S PARABOLIC-HYPERBOLIC COMPRESSIBLE NON-ISOTHERMAL MODEL FOR LIQUID CRYSTALS
Electronic Journal of Differential Equations, Vol. 017 (017), No. 3, pp. 1 8. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu LOCAL WELL-POSEDNESS FOR AN ERICKSEN-LESLIE S
More informationTHE RADIUS OF ANALYTICITY FOR SOLUTIONS TO A PROBLEM IN EPITAXIAL GROWTH ON THE TORUS
THE RADIUS OF ANALYTICITY FOR SOLUTIONS TO A PROBLEM IN EPITAXIAL GROWTH ON THE TORUS DAVID M AMBROSE Abstract A certain model for epitaxial film growth has recently attracted attention, with the existence
More informationEXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018
EXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018 While these notes are under construction, I expect there will be many typos. The main reference for this is volume 1 of Hörmander, The analysis of liner
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationRegularity of the density for the stochastic heat equation
Regularity of the density for the stochastic heat equation Carl Mueller 1 Department of Mathematics University of Rochester Rochester, NY 15627 USA email: cmlr@math.rochester.edu David Nualart 2 Department
More informationScientiae Mathematicae Japonicae Online, Vol. 5, (2001), Ryo Ikehata Λ and Tokio Matsuyama y
Scientiae Mathematicae Japonicae Online, Vol. 5, (2), 7 26 7 L 2 -BEHAVIOUR OF SOLUTIONS TO THE LINEAR HEAT AND WAVE EQUATIONS IN EXTERIOR DOMAINS Ryo Ikehata Λ and Tokio Matsuyama y Received November
More informationRANDOM PROPERTIES BENOIT PAUSADER
RANDOM PROPERTIES BENOIT PAUSADER. Quasilinear problems In general, one consider the following trichotomy for nonlinear PDEs: A semilinear problem is a problem where the highest-order terms appears linearly
More informationThe Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:
Oct. 1 The Dirichlet s P rinciple In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: 1. Dirichlet s Principle. u = in, u = g on. ( 1 ) If we multiply
More informationA FOURTH-ORDER ACCURATE FINITE DIFFERENCE SCHEME FOR THE EXTENDED-FISHER-KOLMOGOROV EQUATION. Tlili Kadri and Khaled Omrani
Bull. Korean Math. Soc. 0 (0), No. 0, pp. 0 https://doi.org/0.434/bkms.b603 pissn: 05-8634 / eissn: 34-306 A FOURTH-ORDER ACCURATE FINITE DIFFERENCE SCHEME FOR THE EXTENDED-FISHER-KOLMOGOROV EQUATION Tlili
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationGlobal regularity of a modified Navier-Stokes equation
Global regularity of a modified Navier-Stokes equation Tobias Grafke, Rainer Grauer and Thomas C. Sideris Institut für Theoretische Physik I, Ruhr-Universität Bochum, Germany Department of Mathematics,
More informationThe Research Grants Council of Hong Kong NSFC/RGC Joint Research Scheme Joint Completion Report
RGC Ref.: N_HKBU204/12 NSFC Ref.: 11261160486 (please insert ref. above) The Research Grants Council of Hong Kong NSFC/RGC Joint Research Scheme Joint Completion Report (Please attach a copy of the completion
More informationSpectrum and Exact Controllability of a Hybrid System of Elasticity.
Spectrum and Exact Controllability of a Hybrid System of Elasticity. D. Mercier, January 16, 28 Abstract We consider the exact controllability of a hybrid system consisting of an elastic beam, clamped
More informationCONVERGENCE OF FINITE DIFFERENCE METHOD FOR THE GENERALIZED SOLUTIONS OF SOBOLEV EQUATIONS
J. Korean Math. Soc. 34 (1997), No. 3, pp. 515 531 CONVERGENCE OF FINITE DIFFERENCE METHOD FOR THE GENERALIZED SOLUTIONS OF SOBOLEV EQUATIONS S. K. CHUNG, A.K.PANI AND M. G. PARK ABSTRACT. In this paper,
More informationarxiv: v1 [math.ap] 11 Apr 2019
Global well-posedness of the 3D Cahn-Hilliard equations Zhenbang Li Caifeng Liu arxiv:194.6191v1 [math.ap] 11 Apr 219 April 15, 219 Abstract The Cauchy problem of the Cahn-Hilliard equations is studied
More informationERROR ESTIMATES FOR SEMI-DISCRETE GAUGE METHODS FOR THE NAVIER-STOKES EQUATIONS : FIRST-ORDER SCHEMES
MATHEMATICS OF COMPUTATION Volume, Number, Pages S 5-578XX- ERROR ESTIMATES FOR SEMI-DISCRETE GAUGE METHODS FOR THE NAVIER-STOKES EQUATIONS : FIRST-ORDER SCHEMES RICARDO H. NOCHETTO AND JAE-HONG PYO Abstract.
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationSharp blow-up criteria for the Davey-Stewartson system in R 3
Dynamics of PDE, Vol.8, No., 9-60, 011 Sharp blow-up criteria for the Davey-Stewartson system in R Jian Zhang Shihui Zhu Communicated by Y. Charles Li, received October 7, 010. Abstract. In this paper,
More informationMAXIMAL AVERAGE ALONG VARIABLE LINES. 1. Introduction
MAXIMAL AVERAGE ALONG VARIABLE LINES JOONIL KIM Abstract. We prove the L p boundedness of the maximal operator associated with a family of lines l x = {(x, x 2) t(, a(x )) : t [0, )} when a is a positive
More informationCOMPARISON OF DIFFERENT NUMERICAL SCHEMES FOR THE CAHN HILLIARD EQUATION
J KSIAM Vol17, No3, 197 207, 2013 http://dxdoiorg/1012941/jksiam201317197 COMPARISON OF DIFFERENT NUMERICAL SCHEMES FOR THE CAHN HILLIARD EQUATION SEUNGGYU LEE 1, CHAEYOUNG LEE 1, HYUN GEUN LEE 2, AND
More informationFourth Order Convergence of Compact Finite Difference Solver for 2D Incompressible Flow
Fourth Order Convergence of Compact Finite Difference Solver for D Incompressible Flow Cheng Wang 1 Institute for Scientific Computing and Applied Mathematics and Department of Mathematics, Indiana University,
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationFinite difference method for heat equation
Finite difference method for heat equation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationNew Identities for Weak KAM Theory
New Identities for Weak KAM Theory Lawrence C. Evans Department of Mathematics University of California, Berkeley Abstract This paper records for the Hamiltonian H = p + W (x) some old and new identities
More informationA POSTERIORI ERROR ESTIMATES BY RECOVERED GRADIENTS IN PARABOLIC FINITE ELEMENT EQUATIONS
A POSTERIORI ERROR ESTIMATES BY RECOVERED GRADIENTS IN PARABOLIC FINITE ELEMENT EQUATIONS D. LEYKEKHMAN AND L. B. WAHLBIN Abstract. This paper considers a posteriori error estimates by averaged gradients
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationarxiv: v1 [math.ap] 18 May 2017
Littlewood-Paley-Stein functions for Schrödinger operators arxiv:175.6794v1 [math.ap] 18 May 217 El Maati Ouhabaz Dedicated to the memory of Abdelghani Bellouquid (2/2/1966 8/31/215) Abstract We study
More informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
More informationA second-order asymptotic-preserving and positive-preserving discontinuous. Galerkin scheme for the Kerr-Debye model. Abstract
A second-order asymptotic-preserving and positive-preserving discontinuous Galerkin scheme for the Kerr-Debye model Juntao Huang 1 and Chi-Wang Shu Abstract In this paper, we develop a second-order asymptotic-preserving
More informationA New Approach for the Numerical Solution of Diffusion Equations with Variable and Degenerate Mobility
A New Approach for the Numerical Solution of Diffusion Equations with Variable and Degenerate Mobility Hector D. Ceniceros Department of Mathematics, University of California Santa Barbara, CA 93106 Carlos
More informationLECTURE 3: DISCRETE GRADIENT FLOWS
LECTURE 3: DISCRETE GRADIENT FLOWS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and
More informationOBSERVABILITY INEQUALITY AND DECAY RATE FOR WAVE EQUATIONS WITH NONLINEAR BOUNDARY CONDITIONS
Electronic Journal of Differential Equations, Vol. 27 (27, No. 6, pp. 2. ISSN: 72-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu OBSERVABILITY INEQUALITY AND DECAY RATE FOR WAVE EQUATIONS
More information1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods
J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods 1 Introduction Achieving high order time-accuracy in the approximation of the incompressible Navier Stokes equations by means of fractional-step
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationPointwise convergence rate for nonlinear conservation. Eitan Tadmor and Tao Tang
Pointwise convergence rate for nonlinear conservation laws Eitan Tadmor and Tao Tang Abstract. We introduce a new method to obtain pointwise error estimates for vanishing viscosity and nite dierence approximations
More informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More informationTADAHIRO OH 0, 3 8 (T R), (1.5) The result in [2] is in fact stated for time-periodic functions: 0, 1 3 (T 2 ). (1.4)
PERIODIC L 4 -STRICHARTZ ESTIMATE FOR KDV TADAHIRO OH 1. Introduction In [], Bourgain proved global well-posedness of the periodic KdV in L T): u t + u xxx + uu x 0, x, t) T R. 1.1) The key ingredient
More informationISABELLE GALLAGHER AND MARIUS PAICU
REMARKS ON THE BLOW-UP OF SOLUTIONS TO A TOY MODEL FOR THE NAVIER-STOKES EQUATIONS ISABELLE GALLAGHER AND MARIUS PAICU Abstract. In [14], S. Montgomery-Smith provides a one dimensional model for the three
More informationCoupled second order singular perturbations for phase transitions
Coupled second order singular perturbations for phase transitions CMU 06/09/11 Ana Cristina Barroso, Margarida Baía, Milena Chermisi, JM Introduction Let Ω R d with Lipschitz boundary ( container ) and
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
More informationNULL CONDITION FOR SEMILINEAR WAVE EQUATION WITH VARIABLE COEFFICIENTS. Fabio Catalano
Serdica Math J 25 (999), 32-34 NULL CONDITION FOR SEMILINEAR WAVE EQUATION WITH VARIABLE COEFFICIENTS Fabio Catalano Communicated by V Petkov Abstract In this work we analyse the nonlinear Cauchy problem
More informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More informationFinite difference methods for the diffusion equation
Finite difference methods for the diffusion equation D150, Tillämpade numeriska metoder II Olof Runborg May 0, 003 These notes summarize a part of the material in Chapter 13 of Iserles. They are based
More information30 crete maximum principle, which all imply the bound-preserving property. But most
3 4 7 8 9 3 4 7 A HIGH ORDER ACCURATE BOUND-PRESERVING COMPACT FINITE DIFFERENCE SCHEME FOR SCALAR CONVECTION DIFFUSION EQUATIONS HAO LI, SHUSEN XIE, AND XIANGXIONG ZHANG Abstract We show that the classical
More informationOPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS
PORTUGALIAE MATHEMATICA Vol. 59 Fasc. 2 2002 Nova Série OPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS J. Saint Jean Paulin and H. Zoubairi Abstract: We study a problem of
More informationThe 2-d isentropic compressible Euler equations may have infinitely many solutions which conserve energy
The -d isentropic compressible Euler equations may have infinitely many solutions which conserve energy Simon Markfelder Christian Klingenberg September 15, 017 Dept. of Mathematics, Würzburg University,
More informationTime-dependent variational forms
Time-dependent variational forms Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Oct 30, 2015 PRELIMINARY VERSION
More informationApplications of the compensated compactness method on hyperbolic conservation systems
Applications of the compensated compactness method on hyperbolic conservation systems Yunguang Lu Department of Mathematics National University of Colombia e-mail:ylu@unal.edu.co ALAMMI 2009 In this talk,
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationSCATTERING FOR THE TWO-DIMENSIONAL NLS WITH EXPONENTIAL NONLINEARITY
SCATTERING FOR THE TWO-DIMENSIONAL NLS WITH EXPONENTIAL NONLINEARITY S. IBRAHIM, M. MAJDOUB, N. MASMOUDI, AND K. NAKANISHI Abstract. We investigate existence and asymptotic completeness of the wave operators
More informationNonlinear Analysis. A regularity criterion for the 3D magneto-micropolar fluid equations in Triebel Lizorkin spaces
Nonlinear Analysis 74 (11) 5 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na A regularity criterion for the 3D magneto-micropolar fluid equations
More information#A69 INTEGERS 13 (2013) OPTIMAL PRIMITIVE SETS WITH RESTRICTED PRIMES
#A69 INTEGERS 3 (203) OPTIMAL PRIMITIVE SETS WITH RESTRICTED PRIMES William D. Banks Department of Mathematics, University of Missouri, Columbia, Missouri bankswd@missouri.edu Greg Martin Department of
More informationSYMMETRY IN REARRANGEMENT OPTIMIZATION PROBLEMS
Electronic Journal of Differential Equations, Vol. 2009(2009), No. 149, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu SYMMETRY IN REARRANGEMENT
More informationON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume, Number, Pages S -9939(XX- ON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS N. S. HOANG AND A. G. RAMM (Communicated
More informationLOCALISATION AND COMPACTNESS PROPERTIES OF THE NAVIER-STOKES GLOBAL REGULARITY PROBLEM
LOCALISATION AND COMPACTNESS PROPERTIES OF THE NAVIER-STOKES GLOBAL REGULARITY PROBLEM TERENCE TAO Abstract. In this paper we establish a number of implications between various qualitative and quantitative
More informationEULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS
Qiao, H. Osaka J. Math. 51 (14), 47 66 EULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS HUIJIE QIAO (Received May 6, 11, revised May 1, 1) Abstract In this paper we show
More information