Large-time behaviour of Hele-Shaw flow with injection or suction for perturbed balls in R N
|
|
- Sandra Lawson
- 6 years ago
- Views:
Transcription
1 Large-time behaviour of Hele-Shaw flow with injection or suction for perturbed balls in R N E. Vondenhoff Department of Mathematics and Computer Science Technische Universiteit Eindhoven P.O. Box 53, 5600 MB Eindhoven, The Netherlands Abstract Large-time behaviour of Hele-Shaw flow with surface tension and with injection or suction in one point is discussed. We consider domains which are initially small perturbations of balls. Radially symmetric solutions are stationary after a suitable time-dependent rescaling. The evolution of perturbations can be described by a non-local nonlinear parabolic evolution equation. Global existence results and decay properties are derived using energy estimates in Sobolev spaces. AMS subject classifications: 35R35, 35K55, 76D7 Key words: moving-boundary problem, non-local parabolic equation Introduction Hele-Shaw flow with surface tension and with injection or suction at a single point is described by a family of domains t Ωt in R N, parameterised by time t, and two functions v, t : Ωt R N and p, t : Ωt R that satisfy the following movingboundary problem: div v = µδ in Ωt, v = p in Ωt, p = γκ on Γt := Ωt. 3 Here, κ, t : Γt R is the mean curvature of the moving boundary t Γt of the domain taken negative if Ωt is convex, µ is the injection rate if µ > 0 or the suction rate if µ < 0, γ is a positive constant, and δ is the Dirac delta distribution. The evolution of the boundary Γt is specified by the requirement that its normal velocity v n is given by v n = v n. 4
2 The equations,, and 3 can be combined to give p = µδ in Ωt, 5 p = γκ on Γt = Ωt. 6 Hence, p solves a Dirichlet problem for any time t and on Γt we get from 4 v n = p n. Besides liquid flow in a Hele-Shaw cell [], the model and variations of it describe the growth of tumours [] and porous media flow [3, 4]. Short-time existence of solutions for a closely related problem has been proved by Escher and Simonett [5]. For the three-dimensional injection problem, global existence of classical solutions and largetime behaviour have been found in [0] for domains that are initially small perturbations of balls. The three-dimensional suction case requires some extra conditions in order to get global existence results see [0] as well. The ratio of suction rate to γ needs to be small and the position of the geometric centre of the initial domain has to be at the origin, where the suction point is located. In this paper we prove global existence results of classical solutions for N = and N 4. Let σ N be the area of the unit sphere S N. It is easy to check that if the initial domain has the same volume as the unit ball σ N N, then the volume V of the domain as a function of time is given by V t = σ N + µt. 7 N The suction problem only makes sense for t T := σ N µn. If Ω0 = B N := {x R N : x < }, then the moving domain will be a ball around the origin with radius αt given by αt = N µnt σ N +. In order to prove stability of this solution we rescale the domain by a factor αt such that B N becomes a stationary solution. We introduce a function r, t : S N R that describes a small star-shaped perturbation of this stationary solution. In [0] an evolution equation for r in Hölder spaces is derived and linearised. This equation is autonomous only for N = 3. Therefore we cannot use the techniques that are used in [0] in other dimensions. Instead, we consider the problem in Sobolev spaces and work with energy estimates. In Section we introduce the same evolution equation on Sobolev spaces. In Section 3 we use the linearisation of the evolution operator to prove global existence and decay properties for perturbations. For the injection problems, we find existence for all t > 0. In the two-dimensional suction problem, we need the extra condition that the geometric centre of the initial domain and the suction point coincide, in order to prove that all liquid can be removed. The domain vanishes as a round point. In contrast to the three-dimensional case, there is no restriction on the suction rate in the two-dimensional case. For N 4, we do not know whether or not a global existence result for the suction problem can be found.
3 dimension µ > 0 µ < 0 global existence and decay global existence for correct geometric centre and decay see Theorem 3.5 see Theorem 3.3 almost global existence see Theorem 4. 3 global existence and decay global existence for correct geometric centre when µ and decay γ < 3π 5 see [0] see [0] almost global existence see [0] 4 global existence almost global existence and decay see Theorem 4. see Theorem 3.6 Table : Large-time behaviour of Hele-Shaw flow with surface tension and injection or suction in one point Both for N = and for N 4, µ < 0, any portion of liquid, smaller than the entire domain, can be removed if the initial domain is close enough to a ball. We shall call this result almost global existence. It will be proved in Section 4. The evolution equation and its linearisation In [0], an evolution equation for the motion of the domain Ωt in Hölder spaces is derived and linearised. In this section we introduce the same equation in a Sobolev space setting. Let H N k be the vector space of harmonic homogeneous polynomials of degree k in N variables. Define the spherical harmonics as restrictions of these polynomials to the unit sphere S N : S N k = { q S N : q H N } k. For each k N 0, choose a basis s k,j νn,k j= of S N k that is orthonormal in L S N : S N k = s k,, s k,,..., s k,νn,k, where νn, k is the dimension of S N k. From [7] Lemma, we know that the spherical harmonics { } sk,, s k,,..., s k,νn,k, k=0 form an orthonormal basis for L S N. Let, 0 be the usual inner product on L S N. For each r L S N define r k,j by r k,j := r, s k,j 0. 3
4 For all s > 0, equip the Sobolev space H s S N with the inner product r, r s = k,j k + s r k,j r k,j. In the sequel we will use the Sobolev embedding theorem: If k N 0, β 0,, and s > N + k + β, then H s S N C k,β S N and H s+ B N C k,β B N. We will also use the fact that for s > N, Hs+ B N and H s S N are Banach algebras. We restrict ourselves to domain evolutions t Ωt that can be described by a continuous function R, t : S N, satisfying { } x Ωt = Ω R,t = x R N \ {0} : x < + R {0}. x Besides R, t we introduce r, t such that which is equivalent to Ω r,t = αt Ω R,t, + rt = + Rt. 8 αt We will often write rt instead of r, t and if we consider a fixed domain, then the argument t will be suppressed. From now on we assume that r H s S N where Define Γ r := Ω r. Introduce zr, : S N Γ r by s > N zr, ξ = + rξ ξ, nr, as the function that maps an element ξ S N to the exterior unit normal vector on Γ r at the point zr, ξ, κr, as the function that maps an element ξ S N to the mean curvature of Γ r at zr, ξ. We will often use the notations zr, nr, and κr instead of zr,, nr,, and κr,. In [8] Chapter 3 Lemma 6, it is proved that on a neighbourhood U of zero in H s S N the operators n : U H s S N N and κ : U H s S N are analytic. By [9] Theorem 6.08, there exists an extension operator E LH s S N, H s+ B N, such that Er S N = r. 0 4
5 N Define z : H s S N H s+ B N by where zr, = zr. zr, x := + Erx x, Lemma.. There exists a δ > 0 such that if r s < δ then zr : B N Ω r is bijective and zr C Ω r N. Proof. From the Sobolev embedding theorem we get H s S N C 4 S N. The bijectivity follows from [] Lemma.. Using [] Lemma.3, we can prove the other statement as well. On a neighbourhood U of zero in H s S N we define the following operators: N A : U LH s 3 B N, H s 7 B N and Q : U L H s 3 B N, H s 5 B N by Aru := u zr zr = i,k,l j i,l r j k,l r u x i x k and Qru := u zr zr = i,k j k,i r u x k e i, where e i is the i-th unit vector in R N and j k,i r are the coefficients of the inverse of the matrix J r = zr N N x H s B N. The elements j k,i r are in H s B N for r s small because of continuity of inversion near the identity J 0 in Banach algebras. Because of 9 the space H s 7 B N is a Banach algebra. Therefore the operators A and Q are well-defined. S : U LH s 3 B N, H s 7 B N H s S N by Aru Sru :=. Tru ϕ : U H s S N by ϕr, x := Φ + rxx, where Φ : R N R is the fundamental solution of the Laplacian: Φx := We often write φr instead of φr,. ln x N =, π N σ N x N N σ N N
6 Lemma.. For r s small, Sr is an isomorphism between H s 3 B N and H s 7 B N H s S N. Proof. S0 is an isomorphism between H s 3 B N and H s 7 B N H s S N. Because S is smooth [] Lemma.6, Sr is an isomorphism for small r. Suppose that 9 holds. By [0] we have for r in a suitable neighbourhood U of zero in H s S N r t = Fr, t, 4 αt 3 where Fr, t = γf r + µαt 3 N F r, 5 for a third order operator F : U H s 3 S N and a first order operator F : U H s S N given by F r = Erκr, F r = Erϕr + lr, where E : U LH s S N, H s 3 S N and l : U H s S N are defined by [ [ ]] 0 Tr Qr Sr ξ nr, ξ ψ Erψ := ξ nr, ξ ξ and lr := σ N + r N + r σ N. In fact, F also maps a neighbourhood of zero in H s S N to H s 3 S N. Introduce a new time variable τ = τt such that τ0 = 0 and For N 3 this implies τt = σ N µ N 3 dτ dt = αt 3. 6 µnt 3 N +. 7 σ N The original time interval was infinite for the injection problems and finite for the suction problems. Via 7, this does not change for the new time variable τ in the case N 4. For N = however, the new injection problem is defined on a finite time interval 0, τ max while the suction problem is defined on 0,. We have lim τt = π, for N =, µ > 0, t µ lim τt =, for N =, µ < 0, t T lim τt =, for N 4, µ > 0, t lim τt = σ N, for N 4, µ < 0. t T µ N 3 6
7 We shall denote these limit values for τ by τ max. We get r τ = Fr, tτ = γf r + µαtτ 3 N F r. Here tτ is the value of t that corresponds to τ. To simplify notation, we will write from now on ατ instead of αtτ and Fr, τ instead of Fr, tτ, so r τ = Fr, τ = γf r + µατ 3 N F r. 8 Note that the evolution operator F is independent of τ in the three dimensional case. In all other dimensions, time dependence only occurs in front of F. Lemma.3. Analyticity of the evolution operator Suppose that s > N+7. The mapping F is analytic from a neighbourhood U of zero in H s S N to H s 3 S N. The mapping F H s 3 S N. is analytic from a neighbourhood U of zero in H s S N to Proof. These statements can be proved in a similar way as was done in [8] Chapter 3 for Hele-Shaw flow without injection or suction and in [] and [0] for F and F in Hölder spaces. The linearisations of F and F around zero are polynomials of the Dirichlet-to- 0 Neumann mapping φ N φ := Tr S0 n0. This is a mapping of order φ that satisfies for all s S N k. From [0] we get and where p and p are the polynomials N s = ks, F 0r = p N r 9 F 0r = p N r 0 p k = k 3 + N k N k and p k = σ N k + N σ N. 7
8 3 Energy estimates and existence results In this section we find estimates for r, Fr, τ s, in order to prove stability of the stationary solution r 0. We use estimates for r, F 0r s and r, F 0r s that are easily obtained from 9 and 0. For the injection problems we prove that for small r0 there exists a global solution rt of 4 that converges to zero if t goes to infinity. For the two dimensional suction problem we need to restrict ourselves to domains with certain geometric properties, in order to get an existence result on the time interval [0, T ]. Let ω be a bijection between {l, m N : l < m N} and {,,..., N }. We define the following differential operators on functions on S N : D ωl,m := x l x m. x m x l For each i {,,..., N }, Di is the infinitesimal generator of a semigroup of operators h R h : R h f = e hdi f = f g h, f L S N, where g h : S N S N are rotations of the unit sphere. Lemma 3.. Let s > N+7. For r Hs+ S N with r s small, we have the generalised chain rule of differentiation: D i F k r = F kr[d i r], k =,. If in addition r H s+ S N, then we have the second order generalised chain rule of differentiation: D i D j F k r = F kr[d i D j r] + F k r[d i r, D j r], k =,. Proof. Because F k k =, commutes with rotations, i.e. we get F k r g h = F k r g h, D i F k r = lim h 0 h R h IF k r = lim h 0 h F kr g h F k r = lim h F kr g h F k r = lim h F kr[r g h r] [ ] = F kr r g h r lim = F h 0 h kr[d i r], h 0 h 0 where I is the identity. From this we can derive the second order generalised chain rule as well see also [6] Section 4 and [8] Chapter 5. Let for σ > 0, σ, be the norm on H σ S N induced by the inner product r, r σ, := r, r σ + i,j D i D j r, D i D j r σ. This norm is equivalent to the norm σ that we introduced earlier. 8
9 Lemma 3.. If r S N k then D ir S N k. Proof. The spaces S N k are invariant under rotations. The lemma follows from this and the fact that D i generates a semigroup of rotations. Theorem 3.3. Let N =, µ > 0, and λ 0 0, µ π. Suppose that s > 5. There exists a δ > 0 such that if r 0 H s S with r 0 s < δ then the problem r τ = Fr, τ, with r0 = r 0, has a solution r C[0, τ max, H s S C [0, τ max, H s 3 S. Furthermore, ξ, τ rτξ C S 0, τ max. If we regard r as a function of the original time variable t, then πλ 0 µt rt s, π + µ r0 s,, t [0,. 3 Proof. The theorem follows from the inequality r, Fr, τ s, λ 0 ατ r s,, 4 for all r H s+3 S with r s small. First we find a similar estimate for the Fréchet derivatives F 0 and F 0. Perturbation arguments and the chain rule lead to 4. Combining 6 and 4 we get algebraic decay of r as a function of t, given by 3. Throughout the proof we will assume that r H s+3 S with r s < δ, where δ is sufficiently small. The symbol C always denotes a sufficiently large constant, which is independent of r.. Take η > 0 such that λ 0 < ηµ π. Define c := γp k + ηµατp k γp k + ηµp k inf = inf > 0, k N 0,τ 0 k + 3 k N 0 k + 3 and define ε := min{c, ηµ π λ 0 }.. Let F be the Frechet derivative of F with respect to the first argument. From 9 and 0 we have the following estimate for the linear part of Fr, τ: r, F 0, τr s = γr, F 0r s + µατr, F 0r s = γr, F 0r s + ηµατr, F 0r s + ηµατr, F 0r s = k + s + 3 γp k ηµατp k k k,j ηατ k + s + µp k k k,j + c r s r k,j r k,j ηµ ατ r π s
10 In the last step we used pk k + π. 3. Now we find an estimate for the remaining nonlinear part. From the analyticity of F and F near zero and the fact that F 0 = F 0 = 0, we have for r in a neighbourhood of zero in H s S F r F 0r s 7 C r s, F r F 0r s 5 C r s. 3 Here the demand s > 5 is crucial. Now we get γr, F r F 0r s + µατr, F r F 0r s C r 3 s + ατ r 3 s From the chain rule we get r, Fr, τ s, = γr, F r s + µατr, F r s +γ D i D j r, F r[d i D j r] s + µατ D i D j r, F r[d i D j r] s i,j i,j +γ i,j D i D j r, F r[d i r, D j r] s + µατ D i D j r, F r[d i r, D j r] s. i,j We divide the right-hand side into three parts and we estimate these parts separately. 5. Adding 5 and 6 we get for the first part of 7 γr, F r s + µατr, F r s c r s c r s ηµ ατ r π s + C r 3 3 s + ατ r 3 s 3 ηµ ατ r π s + Cδ 3 7 r s + ατ r s
11 6. For the second part of 7 we get from similar arguments γd i D j r, F r[d i D j r] s + µατd i D j r, F r[d i D j r] s = γd i D j r, F 0[D i D j r] s + µατd i D j r, F 0[D i D j r] s +γd i D j r, {F r F 0} [D i D j r] s +µατd i D j r, {F r F 0} [D i D j r] s ηµ ατ D i D j r π s 3 c D i D j r s +Cδ D i D j r s + ατ D i D j r s 3. 9 In the last step we used analyticity of F near zero in H s S N and analyticity of F near zero in H s 3 S N. 7. Because of Lemma.3, there exists a C > 0 such that for r near the origin in H s S we have F r X C, for X = L H s S H s S, H s 7 S and F r X C, for X = L H s 3 S H s 3 S, H s 5 S. Therefore, the third part of 7 can be estimated as follows: γd i D j r, F r[d i r, D j r] s + µατd i D j r, F r[d i r, D j r] s C r s+ 3 r s+ + ατ r s+ r s C r s r s+ + ατ r 3 s 3 r s+ Cδ r s+ + ατ r 3 s+ Here we used the following interpolation inequalities:. 30 r s+ r s C r s r s+ 3, C r s 3 r s+. 8. Adding the results of 8, 9, and 30 and using equivalence of the norms σ and σ, we get r, Fr, τ s, c r ηµ s, ατ r π s 3, +Cδ r s, + ατ r s 3,. If we choose δ < ε C, then we get r, Fr, τ s, c ε r s, ηµ π ε ατ r s 3, λ 0 ατ r s 3, λ 0ατ r s,. 3
12 9. In [8] Chapter 6 Proposition 9 and 0, local existence results for Stokes flow with injection or suction are proved. Hele-Shaw flow can be treated in a similar way. Combining this local existence result and 3 we get global existence of a solution r C[0, τ max, H s S C [0, τ max, H s 3 S such that r C 0, τ max, C S. From this it follows that ξ, τ rτξ C S 0, τ max. Furthermore, for small ϑ > 0 and τ > ϑ, we have rτ s, yτ where y : [ϑ, τ max R satisfies dy dτ = λ 0ατy, with yϑ = rϑ s,. Solving this ODE we get for τ > ϑ yτ = e λ0 R τ ϑ α τd τ rϑ s,. After reintroducing the original time variable by 6 we get R t y = e λ0 tϑ α t d t rϑ s, µt tϑ = π πλ 0 µ + rϑ s,. Here, tϑ denotes the value of the original time variable t for τ equal to ϑ. Letting ϑ go to zero we find 3. Define M N = { r C 0 S N : dx = σ } N Ω r N, xdx = 0. Ω r Note that M N is the set of continuous functions on S N for which the corresponding domains Ω r have the volume of the unit ball and a geometric centre that coincides with the origin. It is proved in [0] Lemma 3.7, that if r is a solution of 8 with r0 M N, then rt M N for all t. Introduce the Hilbert spaces H σ S N as H σ S N = {r H σ S N : r, s 0 = 0, s S N 0 S N }. Define f : H s S N R R N by Ω f r := r dx σ N N Ω r xdx Let P : H s S N H s S N be the orthogonal projection onto H s S N with respect to the L S N -inner product and let φ : H s S N R R N H s S N be defined by f r φr :=. P r By the Implicit Function Theorem φ is an analytic bijection between a neighbourhood of zero in H s S N and a neighbourhood of zero in R R N H s S N. This is proved for different function spaces in []. On a suitable neighbourhood U of zero in H s S N we define ψ : U M N by ψ r := φ 0, r..
13 Lemma 3.4. Let s > N+7. For r Hs+ S N with r s small, we have D i D j ψ r = ψ r[d i D j r] + ψ r[d i r, D j r]. 3 Proof. It is sufficient to show that ψ commutes with rotations. If r H s+ S N, then ψ r M N and for a rotation g : R N R N we have r g H s+ S N and ψ r g M N. Because rotations and P commute we have P ψ r g = P ψ r g = r g. Therefore ψ r g = φ 0, P ψ r g = φ 0, r g = ψ r g. This proves the lemma. Now we derive a global existence result for the suction case for N =. Like in the proof of Theorem 3.3, we get this result from energy estimates. The suction case is more complicated then the injection case, first of all because we need to restrict ourselves to the subset M, which is not a vector space. This problem is solved by considering an equivalent problem on H s S given by equation 35 and using the bijection ψ between H s S and M. Once we have found existence of solutions r of this equation, we have existence of corresponding solutions r = ψ r of equation 33. The second complication here is that we need to split up the time interval [0, in two parts, [0, ˆT ] and [ ˆT,. On the first interval, the norm of solutions r that we find might grow up to a value δ. On the second interval we need an energy estimate, that is sharper than the one that we found on the first interval, in order to obtain exponential decay for solutions of 35. For any ratio of µ to γ, a suitable ˆT exists, because for large time surface tension dominates suction. For the three dimensional problem see [0] this is not the case, because eigenvalues of the linearisations of the evolution operators do not change in time. 6γ Theorem 3.5. Let N =, µ < 0, and take λ 0 0, 5. Suppose that s > 5. There 5 exists a δ > 0 and a M > 0 such that if r 0 H s S M with r 0 s < δ then the problem r = Fr, τ, 33 τ with r0 = r 0, has a solution r C[0,, H s S C [0,, H s 3 S with ξ, τ rτξ C S 0,. We have rτ s, Me λ0τ r 0 s,. 34 Proof. Again, the symbol C is used for a constant that may vary throughout the proof.. Note that pk k + 3 Furthermore p k k + 3 such that for τ ˆT and k decreases in k. Therefore we have for k γp k k + 3 6γ 5 5 < λ 0. is bounded and lim τ ατ = 0. Therefore there exists a ˆT γp k + µ ατp k < λ k
14 Choose K N such that for k > K we have γp k + µ p k < 0 and let P K : L S L S be the orthogonal projection with respect to the L S - inner product onto the orthoplement of K k=0 S k. Define c > 0 and c > 0 by and c := γp k µ ατp k γp k µ α inf = inf ˆT p k > λ k,τ ˆT k + 3 k k c := γp k µ ατp k inf = inf k>k,τ 0 k + 3 k>k The positivity of c follows from the fact that γpk µ pk goes to infinity. Define ε := min{c λ 0, c }. γp k µ p k. k + 3 k + 3. Assume for the moment that r satisfies 33. Then r := P r satisfies converges to γ if k r τ = P Fψ r, τ. 35 First we will prove solvability of this equation, finding estimates for r, P Fψ r, τ s, for r H s+3 S and r s < δ with δ small enough. 3. Introduce on a suitable neighbourhood U of zero in H s S the operators F : U H s 3 S and F : U H s S as F = P F ψ and F = P F ψ. These operators are compositions of analytic operators, so they are analytic themselves. From a simple calculation we see that ψ 0 is the identity on H s S. Therefore, F k0 r = F k0 r, for k =, and we have the following estimate for the linear part of γ F +µατ F : γ r, F 0 r s + µατ r, F 0 r s = k Kk + s + 3 γp k + µ ατp k k k>kk + s + 3 γp k + µ ατp k k + 3 r k,j r k,j C r 0 c P K r s = C r 0 + c I P K r s c r s C r 0 c r s. 36 Here we used the fact that I P K : L S H s S is bounded. 4
15 4. For the suction problem we have ατ. Therefore, the nonlinear parts can be estimated in the following way: γ r, F r F 0 r s + µατ r, F r F 0 r s C r 3 s. 37 Here we used the analyticity of F and F, as we did for F and F in the proof of Theorem From the chain rule we get γ r, F r s, + µατ r, F r s, = γ r, F r s + µατ r, F r s +γ i,j D i D j r, F r[d i D j r] s + µατ i,j D i D j r, F r[d i D j r] s +γ i,j D i D j r, F r[d i r, D j r] s + µατ D i D j r, F r[d i r, D j r] s. i,j The right-hand side consists of three parts that will be estimated separately on both intervals [0, ˆT ] and [ ˆT,. We start with [0, ˆT ]. 6. From 36 and 37, we have for the first part γ r, F r s + µατ r, F r s C r 0 c r s + Cδ r s The second part can be treated in a way similar to the proof of Theorem 3.3. We use 36, the boundedness of ατ, and the analyticity of F and F to get γd i D j r, F r[d i D j r] s + µατd i D j r, F r[d i D j r] s = γd i D j r, F 0[D i D j r] s + µατd i D j r, F 0[D i D j r] s +γd i D j r, { F r F 0 } [D i D j r] s +µατd i D j r, { F r F 0 } [D i D j r] s C D i D j r 0 c D i D j r s + Cδ D i D j r s. 39 Here, we also used the fact that D i D j r H s S if r H s+3 S. This follows from Lemma For the third part, we refer to the proof of Theorem 3.3 as well and we use the boundedness of ατ and the analyticity of F and F. In this way we find γd i D j r, F r[d i r, D j r] s + µατd i D j r, F r[d i r, D j r] s Cδ r s
16 9. Combining 38, 39, and 40 and using equivalence of the norms s+ 3 and s, we get on the interval [0, ˆT ] If we take δ < ε C, then we get γ r, F r s, + µατ r, F r s, C r 0, c r s, + Cδ r s,. γ r, F r s, + µατ r, F r s, C r 0, C r s,. 4 Define r 0 := P r 0. Using local existence results as in Theorem 3.3 and diminishing δ if necessary we find an S > 0 such that if r 0 s, δ then 35 has a solution r on [0, S], with r0 = r 0. Take δ < e C ˆT δ and assume r 0 s, δ. Then one can show by induction over K N 0 that 35 has a solution on [KS, K + S] [0, ˆT ] that satisfies rτ s, e Cτ r 0 s, there. Therefore r ˆT s, e C ˆT δ < δ. 0. Now we treat the interval [ ˆT,. Again we consider the chain rule and distinguish between three parts. Because of the boundedness of α we have for the first part γ r, F r s + µατ r, F r s γ r, F 0 r s + µατ r, F 0 r s + C r 3 s = k k + s + 3 γp k + µ ατp k r k + 3 k,j + C r 3 s c r s + Cδ r s. 4 Note that in the summation we start counting from k =, because r H s S.. For the second part we use the same strategy as for the first time interval, to obtain γd i D j r, F r[d i D j r] s + µατd i D j r, F r[d i D j r] s c D i D j r s + Cδ D i D j r s. 43. For the third part, we get exactly the same result as for the first time interval, cf Adding 4, 43, and 40 and using equivalence of the norms s+ 3 and s, we get γ r, F r s, + µατ r, F r s, c r s, + Cδ r s,. 6
17 Taking δ < ε C we find γ r, F r s, + µατ r, F r s, λ 0 r s, λ 0 r s,. Using [8] Chapter 6 Proposition 9 and 0 again, we can extend the solution r that we found on [0, ˆT ] to [ ˆT,, such that for τ [ ˆT, rτ s, e λ0τ ˆT r ˆT s,. Combining the results on both intervals, we get existence of a M > 0 independent of r0, such that for any τ [0, Define rτ s, M e λ0τ r0 s,. r = ψ r. From the smoothness of ψ and the fact that ψ0 = 0 we see that there exists a M > 0 such that if r 0 is small enough, then r is a solution to 33 that satisfies 34. In contrast to the three-dimensional suction problem that is discussed in [0], we do not need any restriction on the ratio of suction rate to γ. Now we derive a theorem for global existence for the higher dimensional case. For injection, we have to deal with the problem that eigenvalues of the linearisation corresponding to spherical harmonics of order zero and one go to zero for large time. The nonlinear part cannot be controlled anymore. In order to deal with this, we use the bijection φ near the origin between H s S N and R R N H s S N and consider the evolution of P r on H s S N and the evolution of f r, the zeroth and first Richardson moments, separately. We write down an equation for P r and find an energy estimate for its evolution operator. This equation differs from the one that we found in the proof of Theorem 3.5 because the evolution may not be in M N. We use the fact that the zeroth and first Richardson moments as function of time are known beforehand. For the suction problem for N 4 we do not get any global existence result because if t approaches T, more and more eigenvalues of the linearised evolution equation become positive. In other words: For large t, F dominates F. Suction can no longer be controlled by surface tension. Theorem 3.6. Let N 4 and µ > 0. Suppose that s > N+8. There exists a δ > 0 and a M > 0 such that if r 0 H s S N with r 0 s < δ then the problem r τ = Fr, τ, 44 with r0 = r 0, has a solution r C[0,, H s S N C [0,, H s 3 S N. Furthermore, ξ, τ rτξ C S N 0,. If we regard r as a function of t conform 7, then M rt s µnt σ N + r 0 s. 7
18 Proof.. Introduce the number c > 0 as c := inf k Choose λ 0 0, c 3 and define ε := c 3 λ 0. γp k. k + 3. It is easily checked that the geometric centre of a moving domain Ω Rt satisfying -4 does not change in time. From this and 7 it follows that solutions r of 44 satisfy where f rτ = V0 V 0 ατ N ατ N+ m 0 =: := f m r 0. 0 Vτ m τ, For notational convenience we introduce q τ := V τ, m τ T. Assume for the moment that r satisfies 44. Then r := P r satisfies r τ = P F φ q τ, r, τ. 45 First we prove solvability of this equation, finding estimates for r, P F φ q τ, r, τ s,, for q 0 small, r H s+3 S N and r s < δ, with δ small enough. 3. Because r H s S N we have γ r, F 0 r s + µατ 3 N r, F 0 r s γ r, F 0 r s c r s Because of Lipschitz continuity of F k φ, for k =,, and because ψ = φ 0, we have P F k φ q τ, r P F k ψ r s = P F k φ q τ, r P F k φ 0, r s C q τ. 47 From a simple calculation it follows that ψ 0 is the identity on H s S N. Therefore the restriction of F k 0 to Hs S N is the Fréchet derivative around zero of the analytic mapping P F k ψ on H s S N and we have P F k ψ r F k0 r s 7 C r s. 48 8
19 Combining 47 and 48 we get the following estimate: { γ r, P F φ q τ, r } s r, F 0 r s +µατ 3 N { r, P F φ q τ, r s r, F 0 r s } C q τ r s + ατ 3 N q τ r s + r 3 s + ατ 3 N r 3 s C q τ r s + r 3 s. 49 Here we used the fact that ατ 3 N. 5. From the chain rule we get r, P F φ q τ, r, τ s, where for k =, = γ F + G + H + µατ 3 N F + G + H, 50 F k = r, P F k φ q τ, r s, G k = i,j D i D j r, P F k φ q τ, r [ D i D j φ q τ, r ] s, H k = i,j D i D j r, P F k φ q τ, r [ D i φ q τ, r, D j φ q τ, r ] s. We will estimate the terms containing F k, G k, and H k separately. 6. Using 46 and 49, we have γf + µατ 3 N F c r s + C q τ r s + C r 3 s Now we find an estimate for the terms of G. Because by definition ψ r = φ 0, r, we have G = i,j D i D j r, I ij + J ij + K ij s, 5 where I ij = P F φ q τ, r [ D i D j φ q τ, r ] P F φ 0, r [ D i D j φ 0, r ], J ij = P F ψ r [ψ r[d i D j r]], K ij = P F ψ r [ψ r[d i r, D j r]]. 9
20 Here we used Lemma 3.4. Because P F φ [D i D j φ ] is Lipschitz continuous from R R N H s+3 S N to H s S N, we get I ij s C q τ. By the Cauchy-Schwarz inequality we have D i D j r, I ij s C q τ D i D j r s. 53 Using 46, the analyticity of F ψ, the fact that ψ 0 is the identity and ψ0 = 0, we get D i D j r, J ij s D i D j r, P F 0 [D i D j r] s + C r s D id j r s c D i D j r s + C r s r s There exists a C > 0, such that for r near the origin in H s S N we have P F ψ r ψ r X C for X = L H s S N H s S N, H s 7 S N. Therefore we have By an interpolation inequality we have K ij s 7 C r s+. D i D j r, K ij s C r s+ 3 r s+ C r s+ r 3 s. 55 Adding 53, 54, and 55 we get for 5 D i D j r, I ij + J ij + K ij s c D i D j r s + C r s+ r 3 s + q τ D i D j r s. For the terms of G we get from similar arguments D i D j r, P F φ q τ, r [ D i D j φ q τ, r ] s C r s+ r s 3 + q τ D i D j r s. Here, we used the estimate D i D j r, P F 0D i D j r s 0. This is a consequence of 0. Because ατ 3 N, we have γg + µατ 3 N G c D i D j r s + C r s+ r 3 s + q τ D i D j r s. 56 0
21 8. From arguments that we used in Theorem 3.3 step 7 we obtain γh + µατ 3 N H C r s+ 3 φ q τ, r s+ Again we used an interpolation inequality. 9. Adding 5, 56, and 57 we get for 50 r, P F φ q τ, r, τ s, C q τ r s+ 3 + r s+ r s+ 3 C q τ r s+ 3 + r s r s+ 3 c r s, + C q τ r s, + δ r s, + q τ r s,. 57 c + Cδ r s, + C q τ + c 3 r s, + C q τ 4 + c 3 r s, c 3 + Cδ r s, + C q τ. Here we used Cauchy s inequality and the fact that q τ 4 q τ for small q 0. If we choose δ < ε C, then r, P F φ q τ, r, τ s, λ 0 r s, + C q τ λ 0 r s, + C q 0 ατ N. 0. Arguing as in the proof of Theorem 3.5 we get global existence of a solution r of 45 for fixed q 0 and for r0 = P r 0 small enough. Furthermore, we have rτ s yτ where y : [0, R satisfies dy dτ = λ 0y + C q 0 ατ N, with y0 = P r 0 s. This ODE can be solved using the variation of constants formula: We have τ 0 yτ = e λ0τ y0 + C q 0 τ e λ0 τ τ τ α τ N d τ e λ0 τ τ d τ + 0 λ 0 0 e λ0 τ τ α τ N d τ. α τ N e λ0τ e λ0τ + C α τ C N ατ N. τ τ α τ N e λ0 τ τ d τ
22 We omitted the exponential terms because they are dominated by the algebraic terms. The result is rτ s Ce λ0τ P r 0 s + C ατ N q 0.. Now we construct a solution r of the original problem by setting rτ := φ q τ, rτ. From the boundedness of φ near the origin we get or rτ s Ce λ0τ P r 0 s + rt s Ce λ0τt P r 0 s + C q0 αt N + P r 0 s C αt N r 0 s. C ατ N q 0 58 C αt N q 0 Note that because of 58, if we restrict ourselves to the case r 0 M N, which means q 0 = 0, then we have faster convergence. We do not treat the three dimensional problems here. Finding energy estimates for N = 3 is easier because the evolution operator does not depend on time. The results are similar to those that are found in [0] in a Hölder space setting. 4 Almost global existence results for the suction problems In this section we find almost global existence results for the suction problems. Both the cases N = and N 4 will be treated. For almost global existence we do not need to restrict ourselves to the case r0 M N. Remember that { for N =, τ max := σ N µ N 3 for N 4. Theorem 4.. Let N = or N 4 and µ < 0. Let T + 0, τ max and s > N+8. There exists a δ > 0 such that if r 0 H s S N with r 0 s < δ, then there exists a solution r C[0, T +, H s S N C [0, T +, H s 3 S N of r τ = Fr, τ, 59 with r0 = r 0. Furthermore, ξ, τ rτξ C S N 0, T +.
23 Proof. For the case N =, we argue as in the proof of Theorem 3.5. There, we split up the time interval in two parts. A different approach for the second time interval was necessary there, because we wanted to show global existence and exponential decay in τ assuming that r M. Here we only consider the first time interval, choosing ˆT T +. In the estimates in steps 3-8 of the proof of Theorem 3.5 we replace r by r and F k by F k. All estimates that are found for the evolution operators F k k =, on the first interval hold for the operators F k as well, because up to equation 4 we did not use the fact that r M. In this way we derive that if r s < δ, for δ small, then we have r, Fr, τ s, C r s,. We choose δ < δ e CT+ and use local existence results as before to prove the theorem for N =. For N 4, ατ 3 N goes to infinity if τ approaches τ max. However, on the time interval [0, T +, ατ 3 N is bounded. Therefore we can use the same strategy as in the proof of Theorem 3.5 on the first of the two intervals to prove the theorem. References [] Chen, X., and Friedman, A. A free boundary problem for an elliptic-hyperbolic system: an application to tumor growth. SIAM J. Math. Anal. 35, 4 003, electronic. [] Elliott, C. M., and Ockendon, J. R. Weak and variational methods for moving boundary problems, vol. 59 of Research Notes in Mathematics. Pitman Advanced Publishing Program, Boston, Mass., 98. [3] Escher, J., and Simonett, G. Maximal regularity for a free boundary problem. NoDEA Nonlinear Differential Equations Appl., 4 995, [4] Escher, J., and Simonett, G. Analyticity of the interface in a free boundary problem. Math. Ann. 305, 3 996, [5] Escher, J., and Simonett, G. Classical solutions for Hele-Shaw models with surface tension. Adv. Differential Equations, 4 997, [6] Günther, M., and Prokert, G. Existence results for the quasistationary motion of a free capillary liquid drop. Z. Anal. Anwendungen 6, 997, [7] Müller, C. Spherical harmonics, vol. 7 of Lecture Notes in Mathematics. Springer- Verlag, Berlin, 966. [8] Prokert, G. Parabolic evolution equations for quasistationary free boundary problems in capillary fluid mechanics. Technische Universiteit Eindhoven, Eindhoven, 997. Dissertation, Technische Universiteit Eindhoven, Eindhoven, 997. [9] Renardy, M., and Rogers, R. C. An introduction to partial differential equations, vol. 3 of Texts in Applied Mathematics. Springer-Verlag, New York,
24 [0] Vondenhoff, E. Asymptotic behaviour of injection and suction for Hele-Shaw flow in R 3 with surface tension near balls. CASA-Report Technische Universiteit Eindhoven, Eindhoven, 006. [] Vondenhoff, E. Long-time behaviour of classical Hele-Shaw flows with injection near expanding balls. CASA-Report Technische Universiteit Eindhoven, Eindhoven,
Long-time behaviour of classical Hele-Shaw flows with injection near expanding balls
Long-time behaviour of classical Hele-Shaw flows with injection near expanding balls E Vondenhoff Dep of Mathematics and Computer Science Technische Universiteit Eindhoven PO Box 513, 5600 MB Eindhoven,
More informationTakens embedding theorem for infinite-dimensional dynamical systems
Takens embedding theorem for infinite-dimensional dynamical systems James C. Robinson Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. E-mail: jcr@maths.warwick.ac.uk Abstract. Takens
More informationClassical solutions for the quasi-stationary Stefan problem with surface tension
Classical solutions for the quasi-stationary Stefan problem with surface tension Joachim Escher, Gieri Simonett We show that the quasi-stationary two-phase Stefan problem with surface tension has a unique
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationRigidity and Non-rigidity Results on the Sphere
Rigidity and Non-rigidity Results on the Sphere Fengbo Hang Xiaodong Wang Department of Mathematics Michigan State University Oct., 00 1 Introduction It is a simple consequence of the maximum principle
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 informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More information08a. Operators on Hilbert spaces. 1. Boundedness, continuity, operator norms
(February 24, 2017) 08a. Operators on Hilbert spaces Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ [This document is http://www.math.umn.edu/ garrett/m/real/notes 2016-17/08a-ops
More informationBIHARMONIC WAVE MAPS INTO SPHERES
BIHARMONIC WAVE MAPS INTO SPHERES SEBASTIAN HERR, TOBIAS LAMM, AND ROLAND SCHNAUBELT Abstract. A global weak solution of the biharmonic wave map equation in the energy space for spherical targets is constructed.
More informationNon-degeneracy of perturbed solutions of semilinear partial differential equations
Non-degeneracy of perturbed solutions of semilinear partial differential equations Robert Magnus, Olivier Moschetta Abstract The equation u + F(V (εx, u = 0 is considered in R n. For small ε > 0 it is
More informationTHE L 2 -HODGE THEORY AND REPRESENTATION ON R n
THE L 2 -HODGE THEORY AND REPRESENTATION ON R n BAISHENG YAN Abstract. We present an elementary L 2 -Hodge theory on whole R n based on the minimization principle of the calculus of variations and some
More informationStability of an abstract wave equation with delay and a Kelvin Voigt damping
Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability
More informationOn Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma. Ben Schweizer 1
On Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma Ben Schweizer 1 January 16, 2017 Abstract: We study connections between four different types of results that
More information2 A Model, Harmonic Map, Problem
ELLIPTIC SYSTEMS JOHN E. HUTCHINSON Department of Mathematics School of Mathematical Sciences, A.N.U. 1 Introduction Elliptic equations model the behaviour of scalar quantities u, such as temperature or
More informationCONNECTIONS BETWEEN A CONJECTURE OF SCHIFFER S AND INCOMPRESSIBLE FLUID MECHANICS
CONNECTIONS BETWEEN A CONJECTURE OF SCHIFFER S AND INCOMPRESSIBLE FLUID MECHANICS JAMES P. KELLIHER Abstract. We demonstrate connections that exists between a conjecture of Schiffer s (which is equivalent
More informationComplex geometrical optics solutions for Lipschitz conductivities
Rev. Mat. Iberoamericana 19 (2003), 57 72 Complex geometrical optics solutions for Lipschitz conductivities Lassi Päivärinta, Alexander Panchenko and Gunther Uhlmann Abstract We prove the existence of
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationBLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED
BLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED TAOUFIK HMIDI AND SAHBI KERAANI Abstract. In this note we prove a refined version of compactness lemma adapted to the blowup analysis
More informationAnalysis in weighted spaces : preliminary version
Analysis in weighted spaces : preliminary version Frank Pacard To cite this version: Frank Pacard. Analysis in weighted spaces : preliminary version. 3rd cycle. Téhéran (Iran, 2006, pp.75.
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationGeometry and the Kato square root problem
Geometry and the Kato square root problem Lashi Bandara Centre for Mathematics and its Applications Australian National University 7 June 2013 Geometric Analysis Seminar University of Wollongong Lashi
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 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 informationNon-degeneracy of perturbed solutions of semilinear partial differential equations
Non-degeneracy of perturbed solutions of semilinear partial differential equations Robert Magnus, Olivier Moschetta Abstract The equation u + FV εx, u = 0 is considered in R n. For small ε > 0 it is shown
More informationSYMMETRY OF POSITIVE SOLUTIONS OF SOME NONLINEAR EQUATIONS. M. Grossi S. Kesavan F. Pacella M. Ramaswamy. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 12, 1998, 47 59 SYMMETRY OF POSITIVE SOLUTIONS OF SOME NONLINEAR EQUATIONS M. Grossi S. Kesavan F. Pacella M. Ramaswamy
More informationLIFE SPAN OF BLOW-UP SOLUTIONS FOR HIGHER-ORDER SEMILINEAR PARABOLIC EQUATIONS
Electronic Journal of Differential Equations, Vol. 21(21), No. 17, pp. 1 9. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu LIFE SPAN OF BLOW-UP
More informationMEAN CURVATURE FLOW OF ENTIRE GRAPHS EVOLVING AWAY FROM THE HEAT FLOW
MEAN CURVATURE FLOW OF ENTIRE GRAPHS EVOLVING AWAY FROM THE HEAT FLOW GREGORY DRUGAN AND XUAN HIEN NGUYEN Abstract. We present two initial graphs over the entire R n, n 2 for which the mean curvature flow
More informationA note on the σ-algebra of cylinder sets and all that
A note on the σ-algebra of cylinder sets and all that José Luis Silva CCM, Univ. da Madeira, P-9000 Funchal Madeira BiBoS, Univ. of Bielefeld, Germany (luis@dragoeiro.uma.pt) September 1999 Abstract In
More informationThe continuity method
The continuity method The method of continuity is used in conjunction with a priori estimates to prove the existence of suitably regular solutions to elliptic partial differential equations. One crucial
More informationStrong uniqueness for second order elliptic operators with Gevrey coefficients
Strong uniqueness for second order elliptic operators with Gevrey coefficients Ferruccio Colombini, Cataldo Grammatico, Daniel Tataru Abstract We consider here the problem of strong unique continuation
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 informationGeometry and the Kato square root problem
Geometry and the Kato square root problem Lashi Bandara Centre for Mathematics and its Applications Australian National University 7 June 2013 Geometric Analysis Seminar University of Wollongong Lashi
More informationFunctional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...
Functional Analysis Franck Sueur 2018-2019 Contents 1 Metric spaces 1 1.1 Definitions........................................ 1 1.2 Completeness...................................... 3 1.3 Compactness......................................
More informationSome Applications to Lebesgue Points in Variable Exponent Lebesgue Spaces
Çankaya University Journal of Science and Engineering Volume 7 (200), No. 2, 05 3 Some Applications to Lebesgue Points in Variable Exponent Lebesgue Spaces Rabil A. Mashiyev Dicle University, Department
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 informationSolving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels
Solving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels Y.C. Hon and R. Schaback April 9, Abstract This paper solves the Laplace equation u = on domains Ω R 3 by meshless collocation
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationTD M1 EDP 2018 no 2 Elliptic equations: regularity, maximum principle
TD M EDP 08 no Elliptic equations: regularity, maximum principle Estimates in the sup-norm I Let be an open bounded subset of R d of class C. Let A = (a ij ) be a symmetric matrix of functions of class
More informationOn uniqueness in the inverse conductivity problem with local data
On uniqueness in the inverse conductivity problem with local data Victor Isakov June 21, 2006 1 Introduction The inverse condictivity problem with many boundary measurements consists of recovery of conductivity
More informationNumerical 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 informationAsymptotic behavior of the degenerate p Laplacian equation on bounded domains
Asymptotic behavior of the degenerate p Laplacian equation on bounded domains Diana Stan Instituto de Ciencias Matematicas (CSIC), Madrid, Spain UAM, September 19, 2011 Diana Stan (ICMAT & UAM) Nonlinear
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationDeviation Measures and Normals of Convex Bodies
Beiträge zur Algebra und Geometrie Contributions to Algebra Geometry Volume 45 (2004), No. 1, 155-167. Deviation Measures Normals of Convex Bodies Dedicated to Professor August Florian on the occasion
More informationChapter One. The Calderón-Zygmund Theory I: Ellipticity
Chapter One The Calderón-Zygmund Theory I: Ellipticity Our story begins with a classical situation: convolution with homogeneous, Calderón- Zygmund ( kernels on R n. Let S n 1 R n denote the unit sphere
More informationSYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS. Massimo Grosi Filomena Pacella S. L. Yadava. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 21, 2003, 211 226 SYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS Massimo Grosi Filomena Pacella S.
More informationMATH Final Project Mean Curvature Flows
MATH 581 - Final Project Mean Curvature Flows Olivier Mercier April 30, 2012 1 Introduction The mean curvature flow is part of the bigger family of geometric flows, which are flows on a manifold associated
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationRecall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
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 informationRobustness for a Liouville type theorem in exterior domains
Robustness for a Liouville type theorem in exterior domains Juliette Bouhours 1 arxiv:1207.0329v3 [math.ap] 24 Oct 2014 1 UPMC Univ Paris 06, UMR 7598, Laboratoire Jacques-Louis Lions, F-75005, Paris,
More informationOrientation transport
Orientation transport Liviu I. Nicolaescu Dept. of Mathematics University of Notre Dame Notre Dame, IN 46556-4618 nicolaescu.1@nd.edu June 2004 1 S 1 -bundles over 3-manifolds: homological properties Let
More informationSquare roots of perturbed sub-elliptic operators on Lie groups
Square roots of perturbed sub-elliptic operators on Lie groups Lashi Bandara (Joint work with Tom ter Elst, Auckland and Alan McIntosh, ANU) Centre for Mathematics and its Applications Australian National
More informationS chauder Theory. x 2. = log( x 1 + x 2 ) + 1 ( x 1 + x 2 ) 2. ( 5) x 1 + x 2 x 1 + x 2. 2 = 2 x 1. x 1 x 2. 1 x 1.
Sep. 1 9 Intuitively, the solution u to the Poisson equation S chauder Theory u = f 1 should have better regularity than the right hand side f. In particular one expects u to be twice more differentiable
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. Note on the fast decay property of steady Navier-Stokes flows in the whole space
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES Note on the fast decay property of stea Navier-Stokes flows in the whole space Tomoyuki Nakatsuka Preprint No. 15-017 PRAHA 017 Note on the fast
More informationCompressible Euler equations with second sound: asymptotics of discontinuous solutions. Beixiang Fang Reinhard Racke
Universität Konstanz Compressible Euler equations with second sound: asymptotics of discontinuous solutions Beixiang Fang Reinhard Racke Konstanzer Schriften in Mathematik Nr. 306, August 0 ISSN 430-3558
More informationKernel Method: Data Analysis with Positive Definite Kernels
Kernel Method: Data Analysis with Positive Definite Kernels 2. Positive Definite Kernel and Reproducing Kernel Hilbert Space Kenji Fukumizu The Institute of Statistical Mathematics. Graduate University
More informationHARDY INEQUALITIES WITH BOUNDARY TERMS. x 2 dx u 2 dx. (1.2) u 2 = u 2 dx.
Electronic Journal of Differential Equations, Vol. 003(003), No. 3, pp. 1 8. ISSN: 107-6691. UL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) HADY INEQUALITIES
More informationTravelling bubbles in a moving boundary problem of Hele-Shaw type
Travelling bubbles in a moving boundary problem of Hele-Shaw type G. Prokert, TU Eindhoven Center for Analysis, Scientific Computing, and Applications (CASA) g.prokert@tue.nl Joint work with M. Günther,
More information4 Hilbert spaces. The proof of the Hilbert basis theorem is not mathematics, it is theology. Camille Jordan
The proof of the Hilbert basis theorem is not mathematics, it is theology. Camille Jordan Wir müssen wissen, wir werden wissen. David Hilbert We now continue to study a special class of Banach spaces,
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationL19: Fredholm theory. where E u = u T X and J u = Formally, J-holomorphic curves are just 1
L19: Fredholm theory We want to understand how to make moduli spaces of J-holomorphic curves, and once we have them, how to make them smooth and compute their dimensions. Fix (X, ω, J) and, for definiteness,
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
More informationMEASURE VALUED DIRECTIONAL SPARSITY FOR PARABOLIC OPTIMAL CONTROL PROBLEMS
MEASURE VALUED DIRECTIONAL SPARSITY FOR PARABOLIC OPTIMAL CONTROL PROBLEMS KARL KUNISCH, KONSTANTIN PIEPER, AND BORIS VEXLER Abstract. A directional sparsity framework allowing for measure valued controls
More informationChapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems
Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter
More informationLinear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems
Linear ODEs p. 1 Linear ODEs Existence of solutions to linear IVPs Resolvent matrix Autonomous linear systems Linear ODEs Definition (Linear ODE) A linear ODE is a differential equation taking the form
More informationb i (µ, x, s) ei ϕ(x) µ s (dx) ds (2) i=1
NONLINEAR EVOLTION EQATIONS FOR MEASRES ON INFINITE DIMENSIONAL SPACES V.I. Bogachev 1, G. Da Prato 2, M. Röckner 3, S.V. Shaposhnikov 1 The goal of this work is to prove the existence of a solution to
More informationLinear Ordinary Differential Equations
MTH.B402; Sect. 1 20180703) 2 Linear Ordinary Differential Equations Preliminaries: Matrix Norms. Denote by M n R) the set of n n matrix with real components, which can be identified the vector space R
More informationPartial Density Functions and Hele-Shaw Flow
Partial Density Functions and Hele-Shaw Flow Jackson Van Dyke University of California, Berkeley jayteeveedee@berkeley.edu August 14, 2017 Jackson Van Dyke (UCB) Partial Density Functions August 14, 2017
More informationNotation. General. Notation Description See. Sets, Functions, and Spaces. a b & a b The minimum and the maximum of a and b
Notation General Notation Description See a b & a b The minimum and the maximum of a and b a + & a f S u The non-negative part, a 0, and non-positive part, (a 0) of a R The restriction of the function
More informationFréchet differentiability of the norm of L p -spaces associated with arbitrary von Neumann algebras
Fréchet differentiability of the norm of L p -spaces associated with arbitrary von Neumann algebras (Part I) Fedor Sukochev (joint work with D. Potapov, A. Tomskova and D. Zanin) University of NSW, AUSTRALIA
More informationASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS. Tian Ma. Shouhong Wang
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 11, Number 1, July 004 pp. 189 04 ASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS Tian Ma Department of
More informationand finally, any second order divergence form elliptic operator
Supporting Information: Mathematical proofs Preliminaries Let be an arbitrary bounded open set in R n and let L be any elliptic differential operator associated to a symmetric positive bilinear form B
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 informationError estimates for the Raviart-Thomas interpolation under the maximum angle condition
Error estimates for the Raviart-Thomas interpolation under the maximum angle condition Ricardo G. Durán and Ariel L. Lombardi Abstract. The classical error analysis for the Raviart-Thomas interpolation
More informationCUTOFF RESOLVENT ESTIMATES AND THE SEMILINEAR SCHRÖDINGER EQUATION
CUTOFF RESOLVENT ESTIMATES AND THE SEMILINEAR SCHRÖDINGER EQUATION HANS CHRISTIANSON Abstract. This paper shows how abstract resolvent estimates imply local smoothing for solutions to the Schrödinger equation.
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
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 informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationMemoirs on Differential Equations and Mathematical Physics
Memoirs on Differential Equations and Mathematical Physics Volume 51, 010, 93 108 Said Kouachi and Belgacem Rebiai INVARIANT REGIONS AND THE GLOBAL EXISTENCE FOR REACTION-DIFFUSION SYSTEMS WITH A TRIDIAGONAL
More informationCONVOLUTION OPERATORS IN INFINITE DIMENSION
PORTUGALIAE MATHEMATICA Vol. 51 Fasc. 4 1994 CONVOLUTION OPERATORS IN INFINITE DIMENSION Nguyen Van Khue and Nguyen Dinh Sang 1 Introduction Let E be a complete convex bornological vector space (denoted
More informationResearch Article Existence and Uniqueness of the Solution for a Time-Fractional Diffusion Equation with Robin Boundary Condition
Abstract and Applied Analysis Volume 211, Article ID 32193, 11 pages doi:1.1155/211/32193 Research Article Existence and Uniqueness of the Solution for a Time-Fractional Diffusion Equation with Robin Boundary
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More informationFourier Transform & Sobolev Spaces
Fourier Transform & Sobolev Spaces Michael Reiter, Arthur Schuster Summer Term 2008 Abstract We introduce the concept of weak derivative that allows us to define new interesting Hilbert spaces the Sobolev
More informationWeak Solutions to Nonlinear Parabolic Problems with Variable Exponent
International Journal of Mathematical Analysis Vol. 1, 216, no. 12, 553-564 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ijma.216.6223 Weak Solutions to Nonlinear Parabolic Problems with Variable
More informationSplines which are piecewise solutions of polyharmonic equation
Splines which are piecewise solutions of polyharmonic equation Ognyan Kounchev March 25, 2006 Abstract This paper appeared in Proceedings of the Conference Curves and Surfaces, Chamonix, 1993 1 Introduction
More informationEquilibria with a nontrivial nodal set and the dynamics of parabolic equations on symmetric domains
Equilibria with a nontrivial nodal set and the dynamics of parabolic equations on symmetric domains J. Földes Department of Mathematics, Univerité Libre de Bruxelles 1050 Brussels, Belgium P. Poláčik School
More informationSome lecture notes for Math 6050E: PDEs, Fall 2016
Some lecture notes for Math 65E: PDEs, Fall 216 Tianling Jin December 1, 216 1 Variational methods We discuss an example of the use of variational methods in obtaining existence of solutions. Theorem 1.1.
More informationCOMPLEX SPHERICAL WAVES AND INVERSE PROBLEMS IN UNBOUNDED DOMAINS
COMPLEX SPHERICAL WAVES AND INVERSE PROBLEMS IN UNBOUNDED DOMAINS MIKKO SALO AND JENN-NAN WANG Abstract. This work is motivated by the inverse conductivity problem of identifying an embedded object in
More informationERRATUM TO AFFINE MANIFOLDS, SYZ GEOMETRY AND THE Y VERTEX
ERRATUM TO AFFINE MANIFOLDS, SYZ GEOMETRY AND THE Y VERTEX JOHN LOFTIN, SHING-TUNG YAU, AND ERIC ZASLOW 1. Main result The purpose of this erratum is to correct an error in the proof of the main result
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationOn Shalom Tao s Non-Quantitative Proof of Gromov s Polynomial Growth Theorem
On Shalom Tao s Non-Quantitative Proof of Gromov s Polynomial Growth Theorem Carlos A. De la Cruz Mengual Geometric Group Theory Seminar, HS 2013, ETH Zürich 13.11.2013 1 Towards the statement of Gromov
More informationu( x) = g( y) ds y ( 1 ) U solves u = 0 in U; u = 0 on U. ( 3)
M ath 5 2 7 Fall 2 0 0 9 L ecture 4 ( S ep. 6, 2 0 0 9 ) Properties and Estimates of Laplace s and Poisson s Equations In our last lecture we derived the formulas for the solutions of Poisson s equation
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 informationON QUALITATIVE PROPERTIES OF SOLUTIONS OF QUASILINEAR ELLIPTIC EQUATIONS WITH STRONG DEPENDENCE ON THE GRADIENT
GLASNIK MATEMATIČKI Vol. 49(69)(2014), 369 375 ON QUALITATIVE PROPERTIES OF SOLUTIONS OF QUASILINEAR ELLIPTIC EQUATIONS WITH STRONG DEPENDENCE ON THE GRADIENT Jadranka Kraljević University of Zagreb, Croatia
More informationBlow-up on manifolds with symmetry for the nonlinear Schröding
Blow-up on manifolds with symmetry for the nonlinear Schrödinger equation March, 27 2013 Université de Nice Euclidean L 2 -critical theory Consider the one dimensional equation i t u + u = u 4 u, t > 0,
More informationKadath: a spectral solver and its application to black hole spacetimes
Kadath: a spectral solver and its application to black hole spacetimes Philippe Grandclément Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris F-92195 Meudon, France philippe.grandclement@obspm.fr
More informationA PROJECTED HESSIAN GAUSS-NEWTON ALGORITHM FOR SOLVING SYSTEMS OF NONLINEAR EQUATIONS AND INEQUALITIES
IJMMS 25:6 2001) 397 409 PII. S0161171201002290 http://ijmms.hindawi.com Hindawi Publishing Corp. A PROJECTED HESSIAN GAUSS-NEWTON ALGORITHM FOR SOLVING SYSTEMS OF NONLINEAR EQUATIONS AND INEQUALITIES
More information