Well-posedness of hyperbolic Initial Boundary Value Problems

Size: px
Start display at page:

Download "Well-posedness of hyperbolic Initial Boundary Value Problems"

Transcription

1 Well-poseness of hyperbolic Initial Bounary Value Problems Jean-François Coulombel CNRS & Université Lille 1 Laboratoire e mathématiques Paul Painlevé Cité scientifique VILLENEUVE D ASCQ CEDEX, France jfcoulom@math.univ-lille1.fr March 16, 2004 Abstract Assuming that a hyperbolic initial bounary value problem satsifies an a priori energy estimate with a loss of one tangential erivative, we show a well-poseness result in the sense of Haamar. The coefficients are assume to have only finite smoothness in view of applications to nonlinear problems. This shows that the weak Lopatinskii conition is roughly sufficient to ensure well-poseness in appropriate functional spaces. AMS subject classification: 35L50, 35L40 1 Introuction In this paper, we consier hyperbolic Initial Bounary Value Problems (IBVPs) in several space imensions. Such problems typically rea: t U + A j(t, x) xj U + D(t, x)u = f(t, x), t ]0, T [, x R +, B(t, y) U x =0 = g(t, y), t ]0, T [, y R 1, (1) U t=0 = U 0 (x), x R +. The space variable x lies in the half-space R + := {x = (x 1,..., x ) R /x > 0}, y = (x 1,..., x 1 ) enotes a generic point of R 1, an t = x 0 is the time variable. The A j s an D are square n n matrices, while B is a p n matrix of maximal rank (the integer p is given below). For simplicity, we shall only eal with noncharacteristic problems, but we point out that the analysis can be reprouce with only minor changes for uniformly characteristic problems (we shall go back to this in our final remarks). To prove the well-poseness of (1), there are basically four steps (see e.g. [3] for a complete escription). One first proves a priori energy estimates for smooth solutions. Then one efines a ual problem an shows the existence of weak solutions (this works because the original an the ual problems usually share the same stability properties). The thir step is to show that weak solutions are strong solutions an thus satisfy the energy estimate. Eventually, one constructs solutions of the IBVP. The first step of this analysis is linke to the so-calle (uniform) Lopatinskii conition (or uniform Kreiss-Lopatinskii conition), see [10]. Namely, the uniform Lopatinskii conition yiels an energy estimate in L 2, with no loss of erivative from the source terms (f, g) to the solution U. The secon step relies on Hahn-Banach an Riesz theorems (see [3]). One obtains weak solutions for which it is not possible to apply the a priori energy estimate. Thus, in the thir step, one introuces a tangential mollifier, regularizes the weak 1

2 solution, applies the a priori estimate to the regularize sequence, an passes to the limit. This proceure was introuce in [11]. The fourth step is to take into account the initial atum U 0, an it was first achieve in [18]. In all the above mentionne results, it is crucial that the first step yiels an energy estimate without loss of erivatives. (We shall say that such problems are stable problems). Such an estimate hols either because the bounary conitions are maximally issipative (or strictly issipative, which is even better), either because the uniform Lopatinskii conition is satisfie. However, it is known that this stability conition is not met by some physically interesting problems. Examples of situations where the uniform Lopatinskii conition breaks own are provie by elastoynamics (with the well-known Rayleigh waves [22, 20]), shock waves or contact iscontinuities in compressible flui mechanics, see e.g. [12, 16]. For such nonstable problems, there is no L 2 estimate, but in some weakly stable situations, one can prove a priori energy estimates with a loss of one tangential erivative from the source terms to the solution. Without entering etails, these problems are those for which the so-calle Lopatinskii eterminant vanishes at orer 1 in the hyperbolic region of the cotangent of the bounary T R t,y R R. For noncharacteristic problems, such energy estimates with loss of one erivative have been erive by the author in [5], an for uniformly characteristic problems, similar energy estimates have been erive by P. Secchi an the author in [6]. (Note that for the Rayleigh waves problem, the Lopatinskii eterminant vanishes in the elliptic region of the cotangent of the bounary, an the situation is slightly better, as shown in [20]). In this paper, we show how to solve the IBVP for such weakly stable problems where losses of erivatives occur. More precisely, we show how to construct solutions of (1), with U 0 = 0, provie that we have an a priori estimate with a loss of one tangential erivative, both for the initial problem (1) an for a ual problem. The construction of a weak solution is quite classical, but still, it requires some attention. Then, we shall regularize our weak solution by using a tangential mollifier. Unlike in the case of stable problems, where any tangential mollifier is suitable, we shall show here that the choice of the mollifier is crucial in our context. Our result is that weak solutions are what we shall call semi-strong solutions. In the en, we shall prove a well-poseness result (in the sense of Haamar) for the IBVP (1), when U 0 = 0. The case of general initial ata is aresse in our final remarks. The paper is organize as follows. In view of possible applications to nonlinear problems, we have chosen to work with low regularity coefficients. Of course, this choice will introuce technical ifficulties, an we have foun it appropriate to give in section 2 all the notations an results on paraifferential calculus that will be use throughout this paper. In section 3, we state precisely our weak stability assumption, an give our main result. In section 4, we prove that (1) amits weak solutions, an that these weak solutions are semi-strong solutions. Up to a few technical etails, this ensures well-poseness for zero initial ata. In section 5, we give some extensions of our results, an make a few comments. 2 Paraifferential calculus with a parameter In this section, we collect some efinitions an results on paraifferential calculus. We refer to the original works by Bony an Meyer [1, 15] an also to [14, 17] for the introuction of the parameter. The reaer will fin etaile proofs in these references. We first introuce some norms on the usual Sobolev spaces. For all γ 1, an for all s R, we equip the space H s (R ) with the following norm: u 2 s,γ := 1 (2π) R λ 2s,γ (ξ) û(ξ) 2 ξ, λ s,γ (ξ) := (γ 2 + ξ 2 ) s/2. We shall write 0 rather than 0,γ for the (usual) L 2 norm. 2

3 The classification of paraifferential symbols (with a parameter) is the following: Definition 2.1. A paraifferential symbol of egree m R an regularity k (k N) is a function a(x, ξ, γ) : R R [1, + [ C q q such that a is C with respect to ξ an for all α N, there exists a constant C α verifying (ξ, γ), α ξ a(, ξ, γ) W k, (R ) C α λ m α,γ (ξ). The set of paraifferential symbols of egree m an regularity k is enote by Γ m k (R ). It is equippe with the obvious semi-norms. We enote by Σ m k (R ) the subset of paraifferential symbols a Γ m k (R ) such that for a suitable ε ]0, 1[ the partial Fourier transform of a satisfies (ξ, γ), Supp F x a(, ξ, γ) {ζ R / ζ ε (γ 2 + ξ 2 ) 1/2 }. Of course, the symbols in Σ m k (R ) are C functions with respect to both variables x an ξ, an for all a Σ m k (R ), we have the estimates (x, ξ, γ), β x α ξ a (x, ξ, γ) C α,β λ m α + β,γ (ξ). Thus any symbol a Σ m k (R ) belongs to Hörmaner s class S1,1 m [9] an efines an operator Op γ (a) on the Schwartz class S(R ) by the usual formula u S(R ), x R, Op γ (a)u(x) := 1 (2π) e ix ξ a(x, ξ, γ) û(ξ) ξ. R We shall use the following terminology: Definition 2.2. A family of operators {P γ } efine for γ 1 will be sai of orer m (m R) if the operators P γ are uniformly boune from H s+m (R ) to H s (R ) for all s, inepenently of γ: γ 1, u H s+m (R ), P γ u s,γ C s u s+m,γ. The following Theorem is crucial: Theorem 2.1. If a Σ m k (R ), k N an m R, the family {Op γ (a)} is of orer m. More precisely, for all s R, there exists a positive constant C such that γ 1, u H s+m (R ), Op γ (a)u s,γ C u s+m,γ. The constant C only epens on s, m, on the confinement parameter ε ]0, 1[, an on a finite number N of semi-norms of a (N only epens on s an m). The regularization of symbols in the class Γ m k (R ) is achieve by a convolution with amissible cut-off functions: Definition 2.3. Let ψ : R R [1, + [ [0, + [ be a C function such that the following estimates hol for all α, β N : (ζ, ξ, γ), α ζ β ξ ψ (ζ, ξ, γ) C α,β λ α β,γ (ξ). We shall say that ψ is an amissible cut-off function if there exist real numbers 0 < ε 1 < ε 2 < 1 satisfying ψ(ζ, ξ, γ) = 1 if ζ ε 1 (γ 2 + ξ 2 ) 1/2, ψ(ζ, ξ, γ) = 0 if ζ ε 2 (γ 2 + ξ 2 ) 1/2. 3

4 An example of cut-off function is the following: we choose a nonnegative C function χ 0 on R R such that ξ γ 2 1 ξ γ 2 2 = χ 0 (ξ 1, γ 1 ) χ 0 (ξ 2, γ 2 ), { χ 0 (ξ, γ) = 1 if ( γ 2 + ξ 2) 1/2 1/2, χ 0 (ξ, γ) = 0 if ( γ 2 + ξ 2) 1/2 1. We efine a function ϕ 0 (ξ, γ) := χ 0 (ξ/2, γ/2) χ 0 (ξ, γ). Then the function ψ 0 efine by ψ 0 (ζ, ξ, γ) := p 0 χ 0 (2 2 p ζ, 0) ϕ 0 (2 p ξ, 2 p γ) (2) is an amissible cut-off function (one can take ε 1 = 1/16 an ε 2 = 1/2). If ψ is an amissible cut-off function, the inverse Fourier transform K ψ of ψ(, ξ, γ) satisfies (ξ, γ), α ξ Kψ (, ξ, γ) L 1 (R ) C α λ α,γ (ξ). These L 1 bouns for the erivatives α ξ Kψ yiel the following result: Proposition 2.1. Let ψ be an amissible cut-off function. The mapping a σa ψ (x, ξ, γ) := K ψ (x y, ξ, γ) a(y, ξ, γ) y R is continuous from Γ m k (R ) to Σ m k (R ) for all m (the confinement parameter of σa ψ is ε 2 ). If a Γ m 1 (R ), then a σa ψ Γ m 1 0 (R ). In particular, if ψ 1 an ψ 2 are two amissible cut-off functions an a Γ m 1 (R ), then σ ψ 1 a σ ψ 2 a Σ m 1 0 (R ). Fixing an amissible cut-off function ψ, we efine the paraifferential operator T ψ,γ a formula T ψ,γ a := Op γ (σ ψ a ). by the If ψ 1 an ψ 2 are two amissible cut-off functions an a Γ m 1 (R ), then Proposition 2.1 an Theorem 2.1 show that the family {T ψ 1,γ a T ψ 2,γ a } is of orer (m 1). The symbolic calculus is base on the following Theorem: Theorem 2.2. Let a Γ m 1 (R ) an b Γ m 1 (R ). Then ab Γ m+m 1 (R ) an the family {T ψ,γ a T ψ,γ b T ψ,γ ab } γ 1 is of orer m + m 1 for all amissible cut-off function ψ. Let a Γ m 1 (R ). Then the family {(T ψ,γ a ) T ψ,γ a } γ 1 is of orer m 1 for all amissible cut-off function ψ. Let a Γ m 2 (R ) an b Γ m 2 (R ). Then ab Γ m+m 2 (R ) an the family {T ψ,γ a T ψ,γ b T ψ,γ ab T ψ,γ i j ξ j a xj b } γ 1 is of orer m + m 2 for all amissible cut-off function ψ. 4

5 Let a Γ m 2 (R ). Then the family {(Ta ψ,γ ) T ψ,γ a T ψ,γ i j ξ j xj a } γ 1 is of orer m 2 for all amissible cut-off function ψ. An easy consequence of Theorem 2.2 is that, for any symbols a Γ m 2 (R ) an b Γ m 2 (R ) that commute, the remainer T γ a T γ b T γ b T γ a T γ i{a,b} = [T γ a, T γ b ] T γ i{a,b} is of orer m + m 2. Here above, the notation {a, b} stans for the Poisson bracket of a an b: {a, b} := ξj a xj b xj a ξj b. j We now stuy the case of paraproucts: they are efine by the particular choice of ψ 0 as cut-off function, where ψ 0 is efine by (2). We shall write Ta γ instea of T ψ 0,γ a for the associate paraifferential operators. We have the following important result: Theorem 2.3. Let a W 1, (R ), u L 2 (R ) an γ 1. Then we have a u T γ a u 0 C γ a W 1, (R ) u 0, a xj u T γ a xj u 0 C a W 1, (R ) u 0, a u T γ a u 1,γ C a W 1, (R ) u 0, for a suitable constant C that is inepenent of (a, u, γ). If in aition a W 2, (R ), we have a u T γ a u 1,γ C γ a W 2, (R ) u 0, a xj u T γ a xj u 1,γ C a W 2, (R ) u 0, for a suitable constant C that is inepenent of (a, u, γ). We can exten the paraifferential calculus to symbols efine on a half-space in the following way: let Ω enote the half-space R ]0, + [= R The space L2 x (Ht,y) s is equippe with the norm u 2 s,γ := + 0 u(, x ) 2 s,γ x. Again, we shall write 0 rather than 0,γ when s = 0 (that is, for the usual norm in L 2 (Ω)). We enote by Γ m k (Ω) the set of symbols a(x 0,..., x, ξ, γ) efine on Ω R [1, + [ such that the mapping x a(, x, ) is boune into Γ m k (R ). We efine the paraifferential operator Ta γ by the formula u C 0 (Ω), x 0, (T γ a u)(, x ) := T γ a(x ) u(, x ). Using Theorem 2.3 an integrating with respect to x, we obtain for all symbol a W 1, (Ω) an all u L 2 (Ω) the estimates: a u T γ a u 0 C γ a W 1, (Ω) u 0, a xj u T γ a xj u 0 C a W 1, (Ω) u 0, j = 0,..., 1. When a W 2, (Ω), one obtains an estimate with a gain of two tangential erivatives: a u T γ a u 1,γ C γ a W 2, (Ω) u 0, a xj u T γ a xj u 1,γ C a W 2, (Ω) u 0, j = 0,..., 1. 5

6 3 Statement of the result Recall that Ω enotes the half-space R t,y R + x. We first make the following assumption on the coefficients of (1): Assumption 1. The A j s are efine on Ω an belong to W 2, (Ω). There exists δ > 0 such that for all (t, x) Ω one has et A (t, x) δ. The matrix B is efine on R an belongs to W 2, (R ). It has maximal rank p, where p equals the number of positive eigenvalues of A (that is, the number of incoming characteristics). The system is symmetric hyperbolic, that is, there exists a (real) matrix value mapping S W 2, (Ω) verifying (t, x) Ω, S(t, x) = S(t, x) T, S(t, x) δ I, S(t, x) A j (t, x) = A j (t, x) T S(t, x). We now make our first weak stability assumption on system (1): Assumption 2. For any D 1 W 1, (Ω), an for any symbol D 2 Γ 0 1 (Ω), there exists a constant C (that epens only on δ, A j W 2, (Ω), D 1 W 1, (Ω), B W 2, (R ) an on a finite number of seminorms of the symbol D 2 ) an there exists a constant γ 0 1 such that for all U C0 (Ω) an for all γ γ 0 one has γ U U x =0 2 0 C γ 3 f 2 1,γ + 1 ) γ 2 g 2 1,γ, where f := A 1 γu + x0 U + A j xj U + D 1 U + T γ D 2 U, g := B U. x =0 Before stating our last assumption, we make a couple of remarks. In the erivation of energy estimates, one usually replaces the linear operator U A 1 γu + x0 U + A j xj U + D 1 U + T γ D 2 U by its paraifferential version U T γ (γ+iξ 0 )A 1 1 U + T γ iξ j A 1 A j U + x U + T γ U + T γ A 1 D D 1 2 U, an treats the errors as source terms 1. These errors have the following form: ( ) γ A 1 U T γ U, or A 1 A 1 A j xj U T γ U, or A 1 iξ j A 1 A D 1 U T γ U. j A 1 D 1 To absorb these errors in an estimate with a loss of one tangential erivative, one nees the regularity state in assumptions 1 an 2 for the coefficients (see Theorem 2.3 in the preceeing section). 1 Recall that we use a tangential symbolic calculus for which x is seen as a parameter. 6

7 The crucial point in assumption 2 is that the energy estimate is inepenent of the lower orer term in the interior equation. More precisely, if the energy estimate hols with D 1 = D 2 = 0, it is not clear whether it also hols for arbitrary D 1 an D 2. This is a major ifference with the stable case where there is no loss of erivative (an therefore, one can treat lower orer terms as source terms in energy estimates). In the framework of weakly stable problems, the lower orer terms in the interior equations can not be neglecte an one nees to pay special attention. One way to rephrase assumption 2 is the following: energy estimates with loss of one tangential erivative hol, inepenently of the lower orer terms, an inepenently of their nature (meaning classical, or paraifferential, or a linear combination of the two). We now turn to our last assumption, that is the analogue of assumption 2 for a ual problem. First recall the following efinition: Definition 3.1. A ual problem for (1) is a linear problem that reas: { t V + AT j x j V + D V = f (t, x), t ]0, T [, x R +, M (t, y) V = g x =0 (t, y), t ]0, T [, y R 1, where M is a (n p) n matrix of maximal rank such that (t, y) R, B (t, y) T B(t, y) + M (t, y) T M(t, y) = A (t, y, 0), (3) for suitable p n an (n p) n matrices B an M, an such that B, M, M belong to W 2, (R ). Our final assumption is that the energy estimate with loss of one tangential erivative is also satisfie by one ual problem 2, when the parameter γ is change into γ: Assumption 3. There exists a ual problem (that is, a matrix M satisfying (3)) such that for any D 1 W 1, (Ω), an for any symbol D 2 Γ 0 1 (Ω), there exists a constant C (that epens only on δ, A j W 2, (Ω), D 1 W 1, (Ω), M W 2, (R ) an on a finite number of seminorms of the symbol D 2 ) an there exists a constant γ 0 1 such that for all V C0 (Ω) an for all γ γ 0 one has γ V V x =0 2 0 C γ 3 f 2 1,γ + 1 ) γ 2 g 2 1,γ, where f := (A T ) 1 γv x0 V A T j xj V + D 1 V + T γ D 2 V, g := M V. x =0 In terms of the Lopatinskii conition, assumption 3 means that for one ual problem, the backwar Lopatinskii conition egenerates at orer 1 in the hyperbolic region of the cotangent of the bounary. In practice, one can usually compute explicit ual bounary conitions for which the Lopatinskii eterminant equals that of the original problem (1). Thus the erivation of energy estimates for a ual problem is usually a irect consequence of energy estimates for the original problem. In all what follows, we always make assumptions 1, 2 an 3. The result is the following: Theorem 3.1. Let D W 1, (Ω), an let T > 0. Then, for all functions f(t, x) an g(t, y) verifying: f, t f, x1 f,, x 1 f L 2 (Ω T ), Ω T :=], T [ R +, g H 1 (ω T ), ω T :=], T [ R 1, 2 Note that, with our efinition, the ual problem is not uniquely efine. 7

8 an such that f an g vanish for t < 0, there exists a unique U L 2 (], T [ R +), whose trace on {x = 0} belongs to L 2 (], T [ R 1 ), that vanishes for t < 0, an that is a solution to { t U + A j(t, x) xj U + D(t, x)u = f(t, x), t ], T [, x R +, B(t, y) U = g(t, y), x =0 t ], T [, y R 1. In aition, U C([0, T ]; L 2 (R +)) an the following estimate hols for all t [0, T ] an all real number γ γ 0 : e 2γt U(t) 2 L 2 (R + ) + γ e γs U 2 L 2 (Ω + t) e γs U x =0 2 L 2 (ω t) C γ e γs f 2 L 2 (Ω + 1 t) γ 3 e γs t,y f 2 L 2 (Ω + t) e γs g 2 L 2 (ω + 1 ) t) γ 2 e γs g 2 L 2 (ω t). The constant C an the parameter γ 0 only epen on δ, A j W 2, (Ω), D W 1, (Ω), B W 2, (R ) an M W 2, (R ). 4 Proof of the main result In this section, we first show existence an uniqueness of solutions for the Bounary Value Problem, with source terms (f, g) in weighte spaces. For γ 1, we efine the spaces L 2 γ(ω) := exp(γt)l 2 (Ω), H 1 γ(ω) := exp(γt)h 1 (Ω). We also efine the spaces H(Ω) := {v L 2 (Ω) s.t. t v, x1 v,..., x 1 v L 2 (Ω)} = L 2 (R + x ; H 1 (R t,y)), H γ (Ω) := exp(γt)h(ω) = {v D (Ω) s.t. exp( γt)v H(Ω)}, H(Ω) := {v H(Ω) s.t. t v, x1 v,..., x v H(Ω)}. (4) The spaces L 2 γ(r ) an H 1 γ(r ) are efine in a similar way. The space L 2 γ(ω) is equippe with the obvious norm: v L 2 γ (Ω) := exp( γt)v 0, an the space H γ (Ω) is equippe with the norm v Hγ(Ω) := ṽ 1,γ with ṽ := exp( γt)v. Similarly, the space H 1 γ(r ) is equippe with the norm w H 1 γ (R ) := w 1,γ with w := exp( γt)w. Some elementary, though useful, properties of the spaces H(Ω) an H(Ω) are collecte in appenix A at the en of this paper. In particular, we show that elements of H(Ω) amit a trace in H 3/2 (R ) (though they o not necesarily belong to H 2 (Ω), since the efinition (4) oes not require x 2 f L 2 (Ω)). We consier a zero orer coefficient D W 1, (Ω), that we fix once an for all, an we wish to prove a well-poseness result for the following Bounary Value Problem: { LU := t U + A j(t, x) xj U + D(t, x)u = f(t, x), (t, x) Ω, B(t, y) U = g(t, y), (t, y) (5) x =0 R, when the source terms f an g belong to H γ (Ω) an Hγ(R 1 ), an γ is large. In view of assumption 2, we expect to obtain a unique solution U in L 2 γ(ω) whose trace on the bounary {x = 0} belongs to L 2 γ(r ). In the en, we shall localize this result on a finite time interval. 8

9 For later use, we efine the norm of the coefficients: N := A j W 2, (Ω) + D W 1, (Ω) + B W 2, (R ) + M W 2, (R ), (6) where M is given by assumption 3, an represents the ual bounary conitions. 4.1 Preliminary estimates We first show that the original problem, as well as the ual problem, satisfy an energy estimate in L 2 (H 1 ) when the source terms are in L 2 (that is, we can shift the inices of regularity). More precisely, we have the following result: Lemma 4.1. Let D 1 W 1, (Ω), an let D 2 Γ 0 1 (Ω). There exists a constant C (that epens only on δ, A j W 2, (Ω), D 1 W 1, (Ω), B W 2, (R ) an on a finite number of seminorms of the symbol D 2 ) an there exists a constant γ 1 1 such that for all U C0 (Ω) an for all γ γ 1 one has γ U U x =0 2 0 C γ 3 f 1 2 1,γ + 1 ) γ 2 g 1 2 1,γ, where f 1 := T γ (γ+iξ 0 )A 1 1 U + T γ iξ j A 1 A j U + x U + T γ A 1 D 1+D 2 U, g 1 := T γ B U x =0, an one also has γ U 2 1,γ + U x =0 2 1,γ C γ 3 f ) γ 2 g 2 2 0, where f 2 := A 1 γu + x0 U + A j xj U + D 1 U, g 2 := B U. x =0 (7) Proof. The first inequality is easily prove using the estimates given in Theorem 2.3: A 1 B U x =0 T γ B U x =0 1,γ C B W 1, (R ) U x =0 0, (γu + x 0 U) T γ (γ+iξ 0 )A 1 A 1 A j xj U T γ iξ j A 1 A j U 1,γ C A 1 W 2, (Ω) U 0, U 1,γ C A 1 A j W 2, (Ω) U 0, thanks to assumption 1. Consequently, using assumption 2, the triangle inequality an choosing γ large enough, one can absorb the error terms in the left han sie of the inequality. We now turn to the secon estimate. Let U C0 (Ω), an efine It is clear that we have f := T γ (γ+iξ 0 )A 1 g := T γ B U x =0, 1 U + T γ iξ j A 1 A j W := T γ (γ 2 + ξ 2 ) 1/2 U = T γ λ 1,γ U. U + x U + T γ U, A 1 D 1 W 0 = U 1,γ, an W x =0 0 = U x =0 1,γ, 9

10 an we are thus le to erive an energy estimate of W in L 2. Thanks to assumption 1 an to Theorem 2.2, we compute: T γ (γ+iξ 0 )A 1 1 W + T γ iξ j A 1 A j W + x W + T γ W A 1 D 1 = T γ f + T γ λ U + 1 T γ U + 1,γ i{(γ+iξ 0 )A 1,λ 1,γ } {ξ j A 1 A j,λ 1,γ } Rγ 2 U = T γ f + T γ λ W + 1 T γ W + 1,γ i{(γ+iξ 0 )A 1,λ 1,γ }λ 1,γ {ξ j A 1 A j,λ 1,γ }λ 1,γ Rγ 1 W, where R γ 2 is a family of orer 2, an Rγ 1 is a family of orer 1. In a similar way, we also compute T γ B W = T γ x g + R γ =0 λ 1,γ 1 W, x =0 where, once again, R γ with the symbol is a family of orer 1. Now, we apply the first estimate of Lemma We get D 2 := i{(γ + iξ 0 )A 1 γ W W x =0 2 0 C, 1 λ 1,γ }λ 1,γ {ξ j A 1 A j, λ 1,γ }λ 1,γ Γ 0 1(Ω). γ 3 T γ f + R γ λ 1,γ 1 W 2 1,γ + 1 ) γ 2 T γ g + R γ λ 1,γ 1 W x =0 2 1,γ. Going back to the efinition of W, an choosing γ large enough, we have alreay obtaine the L 2 (H 1 ) estimate for the paraifferential problem, namely: γ U 2 1,γ + U x =0 2 1,γ C γ 3 f ) γ 2 g 2 0. The result now follows from Lemma 4.2 that we give just below. This Lemma enables us to control the istance (in L 2 ) between f an f 2, an the istance between g an g 2 (f 2 an g 2 are efine by (7)). Lemma 4.2. Let γ 1, a W 1, (R ) an v H 1 (R ). Then one has (a T γ a ) v 0 C a W 1, (R ) v 1,γ, for a suitable constant C that oes not epen on γ, a, v. If, in aition, a W 2, (R ), one has (a T γ a ) v 0 C γ a W 2, (R ) v 1,γ, (a T γ a ) xj v 0 C a W 2, (R ) v 1,γ. Proof. We prove Lemma 4.2 for v in the Schwartz class S(R ). The conclusion follows from a ensity/continuity argument. Decompose v as v = v + j xj v j, with v 0 γ v 1,γ, v j 0 v 1,γ. 10

11 This ecomposition hols with We easily get v (ξ) := γ 2 γ 2 + ξ 2 v(ξ), v j(ξ) := iξ j γ 2 + ξ 2 v(ξ). a v T γ a v 0 (a T γ a )v 0 + j (a T γ a ) xj v j 0 C a W 1, (R ) v 1,γ, thanks to Theorem 2.3. The secon part of Lemma 4.2 is prove in the same way, an we omit the etails. When ealing with symbols efine on a half-space, one simply integrates the estimates of Lemma 4.2. The result is a gain of one or two tangential erivatives, epening on the regularity of the multiplicator. Then one can en the proof of Lemma 4.1. The etails are left to the reaer. Of course, similar a priori estimates hol true for the ual problem (with γ change into γ), since the assumptions an the regularity of the coefficients are exactly the same. We therefore have: Lemma 4.3. Let D W 1, (Ω). There exists a constant C (that epens only on δ, A j W 2, (Ω), D W 1, (Ω) an M W 2, (R )) an there exists a constant γ 1 1 such that for all V C 0 (Ω) an for all γ γ 1 one has γ V 2 1,γ + V x =0 2 1,γ C γ 3 f ) γ 2 g 2 0, where f := (A T ) 1 γv x0 V A T j xj V + D V, g := M V. x =0 Recall that M represents the bounary conitions for the ual problem. With the help of our L 2 (H 1 ) estimate, we are going to construct weak solutions of (5). 4.2 Existence of weak solutions This paragraph is evote to the proof of the following result: Proposition 4.1. There exists γ 2 (N, δ) 1 such that for γ γ 2, f H γ (Ω), an g H 1 γ(r ), there exists U L 2 γ(ω) satisfying U x =0 of istributions). H 1/2 γ (R ) an U is a solution to (5) (in the sense Proof. We first commute (5) with the weight exp( γt), an we are le to search a function Ũ L 2 (Ω) that is a solution to { L γ Ũ := γũ + LŨ = f(t, x), (t, x) Ω, B(t, y) Ũ = g(t, y), (t, y) (8) x =0 R, with ( f, g) := exp( γt)(f, g) H(Ω) H 1 (R ). The formal ajoint (L γ ) of the operator L γ is efine by (L γ ) V := γv t V A T j xj V + D T xj A T j V. 11

12 Using relation (3) an formally integrating by parts, (8) reas V C 0 (Ω), Ũ, (Lγ ) V L 2 (Ω) = f, V L 2 (Ω) + g, B V x =0 L 2 (R ) + MŨ x =0, M V x =0 L 2 (R ), where the scalar proucts are enote as follows: U 1, U 2 L 2 (Ω) := U 1 (x) U 2 (x) x, Ω U 1, U 2 L 2 (R ) := U 1 (t, y) U 2 (t, y) ty. R We efine a set of appropriate test functions: } F := {V C0 (Ω) s.t. M V = 0 x=0 C0 (Ω). Thanks to assumption 3 (with D 2 = 0 an D 1 = D T xj A T j ), we observe that the operator (L γ ) is one-to-one in the vector space F. Consequently, we may efine a linear form l with the following formula: V F, l[(l γ ) V ] := f, V L 2 (Ω) + g, B V x =0 L 2 (R ). (9) The following estimate is now a consequence of Lemma 4.3: l[(l γ ) V ] f 1,γ V 1,γ + C g 1,γ V x =0 1,γ ( ) C V 1,γ + V x =0 1,γ C γ 3/2 (Lγ ) V 0, provie that γ is large enough, say γ γ 2 (N, δ). Now we apply Hahn-Banach theorem an we can thus exten l as a (continuous) linear form over the whole space L 2 (Ω). Thanks to Riesz theorem, we conclue that there exists a function Ũ L2 (Ω) verifying: V F, l[(l γ ) V ] = Ũ, (Lγ ) V L 2 (Ω). In particular, it is clear that L γ Ũ = f in the sense of istributions. Using that A is invertible, the trace of Ũ on the bounary {x = 0} is well-efine an belongs to H 1/2 (R ), see [3, chapter 7] 3. Moreover, the following Green s formula hols: V C 0 (Ω), Ũ, (Lγ ) V L 2 (Ω) = f, V L 2 (Ω) + A Ũ x =0, V x =0 H 1/2 (R ),H 1/2 (R ). Combining with the efinition of l, see (9), we obtain: Using (3), we observe that the matrix V F, g BŨ x =0, B V x =0 H 1/2 (R ),H 1/2 (R ) = 0. ( B M ) W 2, (R ) is invertible. We can therefore conclue that BŨ x = g. This completes the proof of Proposition = The proof in [3] is one with C boune coefficients, but it extens to Lipschitzean coefficients by using paraifferential techniques to estimate commutators, see e.g. [4] an appenix B for such estimates. 12

13 We have constructe a solution Ũ L2 (Ω) of (8), whose trace belongs to H 1/2 (R ) (so that the bounary conitions have a clear meaning). We point out that the L 2 (H 1 ) estimate given by Lemma 4.3 is crucial in orer to obtain Ũ L2 (Ω). If we ha only use the L 2 estimate given by assumption 3, we woul have obtaine a solution Ũ L2 (H 1 ). In view of assumption 2, we expect the function Ũ to amit a trace in L2 an to satisfy an appropriate energy estimate, namely: γ Ũ Ũ x =0 2 0 C γ 3 f 2 1,γ + 1 ) γ 2 g 2 1,γ. In the next paragraph, we show that this property hols. In particular, there exists a unique solution of (5) in L 2 γ(ω), an its trace belongs to L 2 γ(r ) when γ is large. Of course, such an existence-uniqueness result will hol inepenently of the zero orer term D. Before showing this result, we first state an analogue of Proposition 4.1 when the source terms are in L 2 γ. As a matter of fact, we want to solve the BVP for source terms with tangential erivatives in L 2 γ, but in the analysis, we shall see that we also nee to solve BVPs with source terms that are only in L 2 γ. Proposition 4.2. There exists γ 2 (N, δ) 1 such that for γ γ 2, f L 2 γ(ω), an g L 2 γ(r ), there exists U L 2 (R + ; Hγ 1 (R )) satisfying U x =0 H 3/2 γ (R ) an U is a solution to (5) (in the sense of istributions). Proof. Most of the proof is similar to the proof of Proposition 4.1. Keeping the same notations for the linear form l, an for the vector space F, an using assumption 2, we easily obtain the existence of Ũ L2 (R + ; H 1 (R )) such that V F, l[(l γ ) V ] = Ũ, (Lγ ) V L 2 (H 1 ),L 2 (H 1 ). In particular, one has L γ Ũ = f in the sense of istributions. The problem is now to give a meaning to the bounary conitions. This is solve by a trace lemma, which we state in appenix C at the en of this paper. Using this result, we can conclue that the trace of Ũ on {x = 0} is well-efine an belongs to H 3/2 (R ). Moreover, the following Green s formula hols: V C 0 (Ω), Ũ, (Lγ ) V L 2 (H 1 ),L 2 (H 1 ) = f, V L 2 (Ω)+ A Ũ x =0, V x =0 H 3/2 (R ),H 3/2 (R ). Using a continuity/ensity argument, the equality hols for all functions V H 2 (Ω) such that M V = 0 on the bounary. As was one in the proof of Proposition 4.1, we obtain: g BŨ x =0, B V x =0 H 3/2 (R ),H 3/2 (R ) = 0, provie that V H 2 (Ω) an M V = 0 on the bounary. Because B an M belong to W 2, (R ), for all function µ H 3/2 (R ), there exists V H 2 (Ω) such that µ = B V, an x =0 M V x =0 = 0. We can therefore conclue that g = BŨ x = Weak=semi-strong The result is the following: Theorem 4.1. Let ( f, g) H(Ω) H 1 (R ), an let Ũ L2 (Ω) be a solution to (8) 4, for γ sufficiently large. Then there exist a sequence (U ν ) in H(Ω), a boune sequence ( ν ) in the set of symbols Γ 0 1 (Ω), an a boune sequence (bν ) in the set of symbols Γ 1 1 (R ), that satisfy the following properties: 4 Recall that the trace of Ũ is automatically in H 1/2 (R ) an the bounary conitions make sense. 13

14 U ν Ũ in L2 (Ω), U ν x =0 Ũ x =0 in H 1/2 (R ), L γ U ν + A T γ ν U ν f in H(Ω), B U ν x =0 + T γ b ν U ν x =0 g in H1 (R ). In particular, Ũ x =0 belongs to L2 (R ) an the following energy estimate hols: γ Ũ Ũ x =0 2 0 C γ 3 f 2 1,γ + 1 ) γ 2 g 2 1,γ. (10) Recall that the space H(Ω) is efine by (4) 5. There is a similar result for solutions of (5) with source terms in L 2 γ. One simply nees to shift the inices (the regularize sequence belongs to H γ (Ω) an so on). The a priori estimate in L 2 (Hγ 1 ) is the inequality (7) in Lemma 4.1. We omit the proof in this case, an focus on Theorem 4.1. As etaile in the introuction, we are going to introuce a tangential mollifier in orer to regularize Ũ. For all ε ]0, 1], we efine the following symbol ϑ ε: (ξ, γ) R [1, + [, ϑ ε (ξ, γ) := as well as the corresponing Fourier multiplier: Θ γ ε := T γ ϑ ε = (γ 2 ε t,y ) 1. 1 γ 2 + ε ξ 2, With slight abuse of notations, we let Θ γ ε act on functions efine over R an on functions efine over the half-space Ω (where we use symbolic calculus with respect to the tangential coorinates, an the Fourier transform has to be unerstoo as a partial Fourier transform). This mollifier is exactly the one use in [7] (after introucing the parameter γ). As we shall see later on, it has some particularly nice commutation properties with the operator L γ (these properties are expresse by relation (2.11) in [7]). Elementary properties of the mollifier Θ γ ε are liste below: Lemma 4.4. Let γ 1, ε ]0, 1] an s R. Then for all v L 2 (R + ; H s (R )), one has: Θ γ ε v s,γ 1 γ 2 v s,γ, Θ γ ε v s+1,γ 1 γ ε v s,γ, Θ γ ε v s+2,γ 1 ε v s,γ. If v L 2 (R + ; H s+2 (R )), one has Θ γ ε v v/γ 2 s,γ ε γ 4 v s+2,γ. In particular, one has Θ γ ε v v/γ 2 s,γ 0 when ε 0, for all v L 2 (R + ; H s (R )). The proof of Theorem 4.1 is base on several estimates of commutators. Before starting the proof, we recall a lemma by Frierichs: 5 Recall also that the trace on {x = 0} of any element v H(Ω) is well-efine an belongs to H 3/2 (R ). In particular, it belongs to H 1 (R ) an the thir point of the Theorem makes sense, see Theorem A.1 in appenix A. 14

15 Lemma 4.5 (Frierichs). Let a W 1, (Ω). There exists a constant C that epens only on a W 1, (Ω) such that, for all γ 1, for all ε ]0, 1], an for all v L 2 (Ω), one has [a, Θ γ ε ] v 1,γ C v 0. Furthermore, one has [a, Θ γ ε ] v 1,γ 0 when ε 0, for all v L 2 (Ω). Let a W 2, (Ω) an j {0,..., 1}. There exists a constant C that epens only on a W 2, (Ω) such that, for all γ 1, for all ε ]0, 1], an for all v L 2 (Ω), one has [a xj T γ a xj, Θ γ ε ] v 1,γ C v 0. Furthermore, one has [a xj T γ a xj, Θ γ ε ] v 1,γ 0 when ε 0, for all v L 2 (Ω). We postpone the proof of Lemma 4.5 to appenix B (the first part is well-known), an we now give the proof of Theorem 4.1. Proof. Define A irect computation yiels U ε := Θ γ ε Ũ L2 (R + x ; H 2 (R t,y)). A 1 Lγ U ε = Θ γ ε (A 1 f) + γ [A 1 1, Θγ ε ] Ũ + [A 1D, Θγ ε ] Ũ + j=0 [A 1 A j xj, Θ γ ε ] Ũ, where L γ is efine by (8), an where we use the convention A 0 = I. Thanks to Lemma 4.4 an to Lemma 4.5, we alreay have Θ γ ε (A 1 f) + γ [A 1, Θγ ε ] Ũ + [A 1D, Θγ ε ] Ũ = 1 γ 2 A 1 f + r ε, r ε 1,γ 0. This is because A 1 ecomposition W 2, (Ω), D W 1, (Ω), an A 1 f L 2 (R + ; H 1 (R )). Using the [A 1 A j xj, Θ γ ε ] Ũ = [A 1A j xj T γ ia 1 A jξ j, Θ γ ε ] Ũ + [T γ an using Lemma 4.5 (recall that A 1 A j W 2, (Ω)), we obtain A 1 Lγ U ε = 1 γ 2 A 1 1 f + r ε + j=0 ia 1 A jξ j, Θ γ ε ] Ũ, [T γ ia 1 A jξ j, Θ γ ε ] Ũ, (11) where r ε 1,γ tens to 0. The remaining commutators are zero orer terms in Ũ, uniformly with respect to ε. Therefore, one cannot neglect them an treat these commutators as source terms. What saves the ay is that these commutators can be ecompose in the following way 6 : [T γ ia 1 A jξ j, Θ γ ε ] Ũ = T γ j,ε U ε + r ε, where j,ε is a symbol in Γ 0 1 (Ω) that is uniformly boune with respect to ε. Inee, we use Theorem 2.2 to compute [T γ ia 1 A jξ j, Θ γ ε ] Ũ = T γ {A 1 A jξ j,ϑ ε} Ũ + R γ ε Ũ = T γ j,ε ϑ ε Ũ + R γ ε Ũ = T γ j,ε U ε + R γ ε Ũ, (12) 6 If ϑ ε ha compact support in ξ, such a ecomposition woul not hol. The choice of the mollifier ϑ ε is therefore crucial. 15

16 where the symbol j,ε is given by j,ε = 1 k=0 an where R γ ε is of orer 1, uniformly with respect to ε: 2ε ξ j ξ k γ 2 + ε ξ 2 x k (A 1 A j) Γ 0 1(Ω), (13) R γ ε V 1,γ C V 0, V L 2 (Ω), ε ]0, 1]. The symbols j,ε efine by (13) are uniformly boune in Γ 0 1 (Ω) with respect to ε. We observe that ϑ ε = γ 2 +εσ ε, with σ ε boune in Γ 2 1 (Ω). We also observe that j,ε = εα j,ε, with α j,ε boune in Γ 2 1 (Ω). Consequently, the remainer Rγ ε in (12) satisfies lim ε 0 Rγ ε V 1,γ = 0, V C0 (Ω). Thanks to the uniform boun on R γ ε, we may conclue that R γ ε Ũ 1,γ tens to zero as ε tens to zero. Using (12) in equation (11), an efining 1 ε := j,ε Γ 0 1(Ω), j=0 we obtain A 1 Lγ U ε + T γ U ε = 1 ε γ 2 A 1 f + r ε, r ε 1,γ 0, (14) an the symbols ε are boune in Γ 0 1 (Ω). Note that (14) also reas x U ε = A 1 1 γu ε + x0 U ε + A j xj U ε + DU ε T γ U ε + 1 ε γ 2 A 1 f + r ε H(Ω), j=0 an we thus have U ε H(Ω), with H(Ω) efine by (4). For the bounary conitions, one procees in an entirely similar way, an gets BU ε x =0 + T γ b ε U ε x =0 = 1 γ 2 g + r ε, r ε 1,γ 0, (15) with symbols b ε boune in Γ 1 1 (R ). Using (14) an (15), the first part of the Theorem is prove, provie that we efine (with slight abuse of notations): U ν := γ 2 U εν H(Ω), ν := εν, b ν := b εν, ε ν := 1 ν + 1. Now we show that Ũ x =0 L2 (R ) an that (10) hols. First note that the operators V A 1 Lγ V + T γ ν V an V BV x =0 are continuous from the space H(Ω) into H(Ω) an from H(Ω) into H 1 (R ) (use Theorem A.1 in appenix A for the bounary operator). Since C0 (Ω) is ense in H(Ω) (see Proposition A.1 in appenix A), it is clear that the a priori energy estimate given by assumption 2 still hols when U H(Ω) (an not only when U C0 (Ω)). 16

17 Thanks to assumption 2, an to the bouneness of ν in Γ 0 1 (Ω), we know that there exists a constant C = C(N, δ) an a positive number γ 3 (N, δ) such that ) γ U ν U ν x =0 2 0 C for all ν N, an for all γ γ 3. Decomposing γ 3 A 1 Lγ U ν + T γ U ν 2 ν 1,γ + 1 γ 2 BU ν 2 x =0 1,γ BU ν x =0 = (BU ν x =0 + T γ b ν U ν x =0 ) T γ b ν U ν x =0, an using (14)-(15), we get (for γ large enough): where γ U ν U ν x =0 2 0 C γ 3 A 1 f + rν int 2 1,γ + 1 γ 2 g + rb ν 2 1,γ r int ν 1,γ 0, r b ν 1,γ 0., ), (16) The sequence (U ν ) is thus boune in L2 x (R ), an therefore, up to extracting a subsequence, =0 it converges weakly in L 2 (R ) towar some function u L 2 (R ). Since the whole sequence (U ν ) converges towar Ũ in x =0 x =0 H 1/2 (R ), this implies Ũ x =0 L2 (R ). Moreover, we know that the sequence (U ν ) converges strongly towar Ũ in L2 (Ω), an (16) yiels γ Ũ Ũ x =0 2 0 C lim inf γ 3 A 1 f + rν int 2 1,γ + 1 ) γ 2 g + rb ν 2 1,γ, C γ 3 f 2 1,γ + 1 ) γ 2 g 2 1,γ. This completes the proof. We summarize Proposition 4.1 an Theorem 4.1 by the following well-poseness result for the Bounary Value Problem (5): Theorem 4.2. Let D W 1, (Ω). There exists γ 3 (N, δ) such that for γ γ 3, f H γ (Ω) an g H 1 γ(r ), there exists a unique solution U L 2 γ(ω) of the following system: { LU = t U + A j(t, x) xj U + D(t, x)u = f(t, x), (t, x) Ω, B(t, y) U x =0 = g(t, y), (t, y) R, This solution satisfies U x =0 L2 γ(r ) an the following estimate hols: γ U 2 L 2 γ (Ω) + U x =0 2 L 2 γ (R ) C γ 3 f 2 H + 1 ) γ(ω) γ 2 g 2 Hγ 1(R ). In aition, there exists a sequence (U ν ) in H 1 γ(ω) that satisfies the following properties: U ν U in L 2 γ(ω), U ν x =0 U x =0 in L 2 γ(r ), LU ν f in L 2 γ(ω), B U ν x =0 g in L2 γ(r ). The last part of the Theorem is prove in [17], using Frierichs lemma (lemma 4.5). We make the following important comments: even though the source terms f an g have tangential erivatives in L 2 γ, there is no hope to prove, for instance, that LU ν converges towar f in H γ (Ω), see Theorem 4.1. Eventually, the convergence of the traces in L 2 γ can be obtaine because we alreay know (thanks to Theorem 4.1) that the trace of U belongs to L 2 γ(r ). When the source terms are only in L 2 γ, there is an analogous result: 17

18 Theorem 4.3. Let D W 1, (Ω). There exists γ 3 (N, δ) such that for γ γ 3, f L 2 γ(ω) an g L 2 γ(r ), there exists a unique solution U L 2 (R + x : Hγ 1 (R t,y)) of the following system: { t U + A j(t, x) xj U + D(t, x)u = f(t, x), (t, x) Ω, B(t, y) U = g(t, y), (t, y) x =0 R, This solution satisfies U x =0 H 1 γ (R ) an the following estimate hols: ( ) 1 γ 3 U 2 L 2 (R + ;Hγ 1 (R )) + γ2 U x =0 2 Hγ 1 (R ) C γ f 2 L 2 γ(ω) + g 2 L 2 γ(r ). One shoul also keep in min that the ual problem amits similar well-poseness results. 4.4 Well-poseness with zero initial ata. En of the proof Now, we show that Theorem 4.2 an Theorem 4.3 yiel a well-poseness result for the Initial Bounary Value Problem with zero initial ata. We first prove that the classical support lemma extens to weakly stable problems: Lemma 4.6. There exists γ 4 (N, δ) such that, if γ γ 4, (f, g) H γ (Ω) H 1 γ(r ) vanish for t < T 0, then the solution U L 2 γ(ω) of (5) vanishes for t < T 0. Moreover, if γ γ 4, f L 2 γ(ω), an g L 2 γ(r ), then the solution U L 2 (R + ; H 1 γ (R )) of (5) also vanishes for t < T 0. Proof. We give the proof when the source terms are in H γ (Ω) Hγ(R 1 ), but the proof is similar when the source terms are in L 2 γ. There is no loss of generality in assuming T 0 = 0. (Otherwise, use a translation t t T 0 ). We fix a function χ C (R) such that χ oes not vanish an { 1, if x 0, χ(t) = exp( t), if x 1. The function (t, y, x ) Ω χ (t)/χ(t) belongs to W 1, (Ω). enough, the only solution in L 2 γ(ω) to the linear problem LV χ (t) V = 0, (t, x) Ω, χ(t) B(t, y) V = 0, (t, y) x =0 R, Consequently, for all γ large (17) is the trivial solution V = 0, thanks to Theorem 4.2 (we use the essential fact that Theorem 4.2 hols for any zero orer term in W 1, (Ω)). Consier some ata (f, g) H γ (Ω) Hγ(R 1 ) that vanish for t < 0. Then we have (f, g) H γ+j (Ω) Hγ+j 1 (R ) for all integer j. Thanks to Theorem 4.2, we know that there exists a unique U j L 2 γ+j (Ω) satisfying { LUj = f(t, x), (t, x) Ω, B(t, y) U j x =0 = g(t, y), (t, y) R. The function χ(u j+1 U j ) belongs to L 2 γ+j (Ω) an one checks that it is a solution to (17). Therefore it equals zero, an U j+1 = U j = = U 0. Furthermore, we know that sup j N 1 γ + j f H γ+j (Ω) < + an sup j N 1 γ + j g H 1 γ+j (R ) < +, 18

19 because f an g vanish for t < 0. Thus Theorem 4.2 yiels sup j This implies that U 0 vanishes for t < 0. U j L 2 γ+j (Ω) = sup U 0 L 2 j γ+j (Ω) < +. We introuce a few notations: for T > 0, let Ω T := Ω {t < T } =], T [ R +, an let ω T :=], T [ R 1. The spaces L 2 γ(ω T ), L 2 γ(ω T ), an H γ (Ω T ) are efine similarly as L 2 γ(ω) etc. The efinition of the norms in L 2 γ(ω T ) an L 2 γ(ω T ) is clear. As regars the norm in H γ (Ω T ), it is efine by 1 f 2 H γ(ω T ) := γ2 f 2 L 2 γ(ω T ) + xj f 2 L 2 γ(ω T ). The norm of H 1 γ(ω T ) is efine in a similar way. We are now able to en the proof of Theorem 3.1. Proof. We consier source terms f H(Ω T ), an g H 1 (ω T ), that vanish in the past. We note that f an g belong to H γ (Ω T ) an to H 1 γ(ω T ) for all γ 1. We exten f an g as functions f H(Ω) an g H 1 (R ). Because f an g vanish for t < 0, it is also straightforwar that f H γ (Ω) an g H 1 γ(r ) for all γ 1. Consequently, for γ large enough, there exists a unique U L 2 γ(ω) such that j=0 { LU = f (t, x), (t, x) Ω, B(U ) x =0 = g (t, y), (t, y) R, an U satisfies the corresponing energy estimate, see Theorem 4.2. Furthermore, U vanishes in the past, thanks to Lemma 4.6. We also have γ U 2 L 2 γ(ω T ) + (U ) x =0 2 L 2 γ(ω T ) γ U 2 L 2 γ(ω) + (U ) x =0 2 L 2 γ(r ) C γ 3 f 2 H + 1 ) γ(ω) γ 2 g 2 Hγ 1(R ) C γ 3 f 2 H γ(ω T ) + 1 ) γ 2 g 2 Hγ 1(ω T ). The restriction U of U to Ω T belongs to L 2 (Ω T ) because U vanishes in the past an U L 2 γ(ω T ) when γ is large. We have thus constructe a solution in L 2 (Ω T ) to the localize problem: { LU = f(t, x), (t, x) ΩT, B(t, y) U = g(t, y), (t, y) ω (18) x =0 T. We now show uniqueness of such a solution. Let U L 2 (Ω T ) have a trace in L 2 (ω T ), vanish in the past, an satisfy { LU = 0, (t, x) ΩT, BU x =0 = 0, (t, y) ω T. Let ε > 0, an consier a function χ C (R) such that χ(t) = 1 if t T 2ε, an χ(t) = 0 if t T ε. Define U χ := χu. Then U χ L 2 (Ω) an U χ vanishes in the past, so U χ L 2 γ(ω) for all γ 1. Moreover, we compute: { LUχ = χ (t)u, (t, x) Ω, BU χ x =0 = 0, (t, y) R. 19

20 Observe that χ U L 2 γ(ω) for all γ 1, an χ U vanishes for t < T 2ε. We can thus apply Lemma 4.6 (with source terms in L 2 γ) 7, an conclue that U χ vanishes for t < T 2ε. Consequently, U vanishes for t < T, that is, U = 0. To en the proof of Theorem 3.1, we show the continuity with respect to the time variable. We consier source terms f H(Ω T ) an g H 1 (ω T ) that vanish in the past, an we continue these functions as f H(Ω) an g H 1 (R ). We alreay know that the unique solution U L 2 (Ω T ) of (18) that vanishes in the past, is the restriction to Ω T of the solution U L 2 γ(ω) to the global problem { LU = f, (t, x) Ω, BU x =0 = g, (t, y) R, Moreover, the trace of U on {x = 0} belongs to L 2 (R ), an vanishes in the past. To prove Theorem 3.1, it is sufficient to show the continuity of U with respect to the time variable. Thanks to Theorem 4.2, we know that there exists a sequence (U ν ) in H 1 γ(ω) verifying U ν U in L 2 γ(ω), U ν x =0 U x =0 in L 2 γ(r ), LU ν LU = f in L 2 γ(ω), B U ν x =0 g in L 2 γ(r ). Using the Frierichs symmetrizer S (see assumption 1), the classical energy estimate for symmetric hyperbolic systems in a half-space reas (see e.g. [11]): e 2γt U ν (t) 2 L 2 (R + ) + γ U ν 2 L 2 γ(ω t) C γ LU ν 2 L 2 γ (Ωt) + U ν x =0 2 L 2 γ (ωt) ), an we also have e 2γt U ν (t) U ν (t) 2 L 2 (R + ) + γ U ν U ν 2 L 2 γ(ω t) ( ) 1 C γ LU ν LU ν 2 L 2 γ(ω + U ν t) U ν x =0 x =0 2 L 2 γ(ω t). Passing to the limit, we obtain the continuity of U with respect to the time variable. previous estimates for the trace of U on ω t yiels the estimate state in Theorem Some remarks 5.1 The IBVP with general initial ata Using Theorem 3.1, one woul like to show well-poseness of the IBVP (1) with general initial ata. Assume first that the coefficients A j s an D, as well as the Frierichs symmetrizer S, are C, boune, an with boune erivatives. Exten those coefficients to the whole space R +1, so that the system remains symmetric hyperbolic. Then for all f L 1 (]0, T [; H 2 (R )), an for all U 0 H 2 (R ), one can construct a solution U (1) C 0 ([0, T ]; H 2 (R )) C 1 ([0, T ]; H 1 (R )) to the Cauchy problem { t U (1) + A j(t, x) xj U (1) + D(t, x)u (1) = f(t, x), t ]0, T [, x R, U (1) t=0 = U 0 (x), x R. For the IBVP (1), one seeks the solution uner the form U = U (1) + U (2), with U (2) solution to t U (2) + A j(t, x) xj U (2) + D(t, x)u (2) = 0, t ]0, T [, x R +, B(t, y)u (2) (1) x = g BU =0, t ]0, T [, y R 1 x, =0 U (2) t=0 = 0, x R. 7 It is crucial here to have a well-poseness result for source terms in L 2 γ. The 20

21 Consequently, if the source term g belongs to H 1 (]0, T [ R 1 ), an if the initial ata U 0 H 2 (R +) satisfy the compatibility conition g t=0 = B(0, y)(u 0 ) x =0, then one can solve the IBVP (1) with a source term f L 1 (]0, T [; H 2 (R +)), thanks to Theorem 3.1. However, this strategy harly applies when the A j s are in W 2, (Ω) an D is only in W 1, (Ω). The problem is to solve the Cauchy problem with initial ata, for instance in H 2 (R ), an to obtain a solution on ]0, T [ R, such that its trace on ]0, T [ R 1 belongs to H 1. This oes not seem possible with a zero orer coefficient in W 1,. We therefore prefer not to pursue this issue, which is a little beyon the scope of this paper. However, for C boune coefficients, an with ata that satisfy the above mentione compatibility conition, the techniques of [18] shoul yiel a well-poseness result for the IBVP (1). The result of [19] even suggests that, uner this compatibility conition, the IBVP is well-pose with initial ata in H 1 (R +). When ealing with nonlinear problems, one usually solves the nonlinear equations by a sequence of linearize problems with zero initial ata an source terms that vanish in the past, see e.g. [12, 14, 17]. This is another reason why we o not pursue the stuy of general initial ata. 5.2 Uniformly characteristic IBVP In applications, it often happens that the bounary is characteristic, that is, the eterminant of the matrix A vanishes on the bounary. In many of these cases (at least in many physically relevant problems), the rank of A is constant only on the bounary, an the bounary conitions are maximally issipative. In such situations, the corresponing IBVP has been stuie in great etails by many authors, even at the level of quasilinear equations, see e.g. [8, 21] an the references cite therein. When the bounary conitions satisfy the uniform Lopatinskii conition, an when the rank of A is constant in a neighborhoo of the bounary, the linear IBVP was stuie in [13]. When losses of tangential erivatives occur, an when the bounary is uniformly characteristic (this happens for instance in the stuy of contact iscontinuities, see [6]), one can reprouce the analysis evelope above. More precisely, assume that there exist two invertible matrices Q 1,2 (t, x) W 2, (Ω), such that 0 n0 (t, x) Ω, Q 1 (t, x)a (t, x)q 2 (t, x) = I n+ I n, where I k is the ientity matrix in R k, an n 0, n ± are fixe integers. Then the problem { t U + A j(t, x) xj U + D(t, x)u = f(t, x), (t, x) Ω T, B(t, y) U = g(t, y), x =0 (t, y) ω T, with source terms f H(Ω T ) an g H 1 (ω T ) that vanish in the past, satisfies a well-poseness result that is entirely analogous to Theorem 3.1, provie that the analogues of assumptions 2 an 3 for characteristic problems are satisfie. The only ifference is that we can control only the noncharacteristic part of the trace of the solution U on the bounary. With the help of our analysis, the verification of the well-poseness of the linearize equations for the vortex sheets problem (as stuie in [6]) is thus essentially reuce to the calculation of the Lopatinskii eterminant for a suitable ual problem. 21

22 Acknowlegments I thank Patrick Gérar for bringing the reference [7] to my attention. I also thank Guy Métivier for enlightening numerous technical etails. Research of the author was supporte by the EU finance network HYKE, HPRN-CT A Some properties of anisotropic Sobolev spaces In R +1, a generic point is enote by x = (x 0,..., x ). We use the notation x = (x, x ) with x R an x R. We also use the notation Ω = R +1 + = {x R+1 s.t. x > 0}. We efine the following spaces H(R +1 ) := {f L 2 (R +1 ) s.t. x0 f,..., x 1 f L 2 (R +1 )}, H(R +1 ) := {f H(R +1 ) s.t. x0 f,..., x f H(R +1 )}, H(Ω) := {f L 2 (Ω) s.t. x0 f,..., x 1 f L 2 (Ω)}, H(Ω) := {f H(Ω) s.t. x0 f,..., x f H(Ω)}. The spaces H(R +1 ) an H(R +1 ) are equippe with the following norms 8 : f 2 H(R +1 ) := 1 (2π) +1 (1 + ξ 2 ) f(ξ) 2 ξ, R +1 f 2 H(R +1 ) := 1 (2π) +1 (1 + ξ 2 )(1 + ξ 2 ) f(ξ) 2 ξ, R +1 where we have ecompose ξ = (ξ, ξ ) R R. The spaces H(Ω) an H(Ω) are equippe with the norms: f 2 H(Ω) := f 2 L 2 (Ω) + x 0 f 2 L 2 (Ω) + + x 1 f 2 L 2 (Ω), f 2 H(Ω) := f 2 H(Ω) + x 0 f 2 H(Ω) + + x f 2 H(Ω). The following ensity result is stanar, an is prove by a truncation/regularization argument (see [2] for etails): Proposition A.1. The space C0 (R+1 ) is ense in both H(R +1 ) an H(R +1 ). The space C0 (Ω) is ense in both H(Ω) an H(Ω). The following result is also very classical: Proposition A.2. There exist two continuous linear mappings E : H(Ω) H(R +1 ) an E : H(Ω) H(R +1 ) such that for all u H(Ω) (resp. u H(Ω)), Eu = u (resp. Eu = u) almost everywhere in Ω. Observe that for the extension operator E, it is sufficient to consier a continuation by 0 outsie of Ω (this is because there is no normal erivative x in the efinition of H(Ω)). We en this short appenix with the following result: Theorem A.1. The mapping Γ : u C 0 (Ω) u(x, 0) C 0 (R ) can be continue in a unique way as a continuous linear mapping Γ : H(Ω) H 3/2 (R ). 8 Here we take γ = 1 for the sake of simplicity, but it is clear that introucing the parameter γ in the efinition of the norms oes not change the results state below. 22

PDE Notes, Lecture #11

PDE Notes, Lecture #11 PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =

More information

The Generalized Incompressible Navier-Stokes Equations in Besov Spaces

The Generalized Incompressible Navier-Stokes Equations in Besov Spaces Dynamics of PDE, Vol1, No4, 381-400, 2004 The Generalize Incompressible Navier-Stokes Equations in Besov Spaces Jiahong Wu Communicate by Charles Li, receive July 21, 2004 Abstract This paper is concerne

More information

Calculus of Variations

Calculus of Variations 16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

REGULARITY OF WEAKLY WELL POSED HYPERBOLIC MIXED PROBLEMS WITH CHARACTERISTIC BOUNDARY

REGULARITY OF WEAKLY WELL POSED HYPERBOLIC MIXED PROBLEMS WITH CHARACTERISTIC BOUNDARY REGULARITY OF WEAKLY WELL POSED HYPERBOLIC MIXED PROBLEMS WITH CHARACTERISTIC BOUNDARY ALESSANDRO MORANDO AND PAOLO SECCHI Abstract. We study the mixed initial-boundary value problem for a linear hyperbolic

More information

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

3 The variational formulation of elliptic PDEs

3 The variational formulation of elliptic PDEs Chapter 3 The variational formulation of elliptic PDEs We now begin the theoretical stuy of elliptic partial ifferential equations an bounary value problems. We will focus on one approach, which is calle

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy, NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which

More information

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics

More information

Agmon Kolmogorov Inequalities on l 2 (Z d )

Agmon Kolmogorov Inequalities on l 2 (Z d ) Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,

More information

Tractability results for weighted Banach spaces of smooth functions

Tractability results for weighted Banach spaces of smooth functions Tractability results for weighte Banach spaces of smooth functions Markus Weimar Mathematisches Institut, Universität Jena Ernst-Abbe-Platz 2, 07740 Jena, Germany email: markus.weimar@uni-jena.e March

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

ON ISENTROPIC APPROXIMATIONS FOR COMPRESSIBLE EULER EQUATIONS

ON ISENTROPIC APPROXIMATIONS FOR COMPRESSIBLE EULER EQUATIONS ON ISENTROPIC APPROXIMATIONS FOR COMPRESSILE EULER EQUATIONS JUNXIONG JIA AND RONGHUA PAN Abstract. In this paper, we first generalize the classical results on Cauchy problem for positive symmetric quasilinear

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES

SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES Communications on Stochastic Analysis Vol. 2, No. 2 (28) 289-36 Serials Publications www.serialspublications.com SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES

More information

Lower bounds on Locality Sensitive Hashing

Lower bounds on Locality Sensitive Hashing Lower bouns on Locality Sensitive Hashing Rajeev Motwani Assaf Naor Rina Panigrahy Abstract Given a metric space (X, X ), c 1, r > 0, an p, q [0, 1], a istribution over mappings H : X N is calle a (r,

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Linear Algebra- Review And Beyond. Lecture 3

Linear Algebra- Review And Beyond. Lecture 3 Linear Algebra- Review An Beyon Lecture 3 This lecture gives a wie range of materials relate to matrix. Matrix is the core of linear algebra, an it s useful in many other fiels. 1 Matrix Matrix is the

More information

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine

More information

The continuity equation

The continuity equation Chapter 6 The continuity equation 61 The equation of continuity It is evient that in a certain region of space the matter entering it must be equal to the matter leaving it Let us consier an infinitesimal

More information

LOCAL AND GLOBAL MINIMALITY RESULTS FOR A NONLOCAL ISOPERIMETRIC PROBLEM ON R N

LOCAL AND GLOBAL MINIMALITY RESULTS FOR A NONLOCAL ISOPERIMETRIC PROBLEM ON R N LOCAL AND GLOBAL MINIMALITY RSULTS FOR A NONLOCAL ISOPRIMTRIC PROBLM ON R N M. BONACINI AND R. CRISTOFRI Abstract. We consier a nonlocal isoperimetric problem efine in the whole space R N, whose nonlocal

More information

A nonlinear inverse problem of the Korteweg-de Vries equation

A nonlinear inverse problem of the Korteweg-de Vries equation Bull. Math. Sci. https://oi.org/0.007/s3373-08-025- A nonlinear inverse problem of the Korteweg-e Vries equation Shengqi Lu Miaochao Chen 2 Qilin Liu 3 Receive: 0 March 207 / Revise: 30 April 208 / Accepte:

More information

INVERSE PROBLEM OF A HYPERBOLIC EQUATION WITH AN INTEGRAL OVERDETERMINATION CONDITION

INVERSE PROBLEM OF A HYPERBOLIC EQUATION WITH AN INTEGRAL OVERDETERMINATION CONDITION Electronic Journal of Differential Equations, Vol. 216 (216), No. 138, pp. 1 7. ISSN: 172-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu INVERSE PROBLEM OF A HYPERBOLIC EQUATION WITH AN

More information

Abstract A nonlinear partial differential equation of the following form is considered:

Abstract A nonlinear partial differential equation of the following form is considered: M P E J Mathematical Physics Electronic Journal ISSN 86-6655 Volume 2, 26 Paper 5 Receive: May 3, 25, Revise: Sep, 26, Accepte: Oct 6, 26 Eitor: C.E. Wayne A Nonlinear Heat Equation with Temperature-Depenent

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

DECOMPOSITION OF POLYNOMIALS AND APPROXIMATE ROOTS

DECOMPOSITION OF POLYNOMIALS AND APPROXIMATE ROOTS DECOMPOSITION OF POLYNOMIALS AND APPROXIMATE ROOTS ARNAUD BODIN Abstract. We state a kin of Eucliian ivision theorem: given a polynomial P (x) an a ivisor of the egree of P, there exist polynomials h(x),

More information

Chaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena

Chaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena Chaos, Solitons an Fractals (7 64 73 Contents lists available at ScienceDirect Chaos, Solitons an Fractals onlinear Science, an onequilibrium an Complex Phenomena journal homepage: www.elsevier.com/locate/chaos

More information

Discrete Operators in Canonical Domains

Discrete Operators in Canonical Domains Discrete Operators in Canonical Domains VLADIMIR VASILYEV Belgoro National Research University Chair of Differential Equations Stuencheskaya 14/1, 308007 Belgoro RUSSIA vlaimir.b.vasilyev@gmail.com Abstract:

More information

The effect of dissipation on solutions of the complex KdV equation

The effect of dissipation on solutions of the complex KdV equation Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

Local well-posedness of nonlocal Burgers equations

Local well-posedness of nonlocal Burgers equations Local well-poseness of nonlocal Burgers equations Sylvie Benzoni-Gavage To cite this version: Sylvie Benzoni-Gavage. Local well-poseness of nonlocal Burgers equations. Differential Integral Equations,

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

LOCAL WELL-POSEDNESS OF NONLINEAR DISPERSIVE EQUATIONS ON MODULATION SPACES

LOCAL WELL-POSEDNESS OF NONLINEAR DISPERSIVE EQUATIONS ON MODULATION SPACES LOCAL WELL-POSEDNESS OF NONLINEAR DISPERSIVE EQUATIONS ON MODULATION SPACES ÁRPÁD BÉNYI AND KASSO A. OKOUDJOU Abstract. By using tools of time-frequency analysis, we obtain some improve local well-poseness

More information

Generalized Tractability for Multivariate Problems

Generalized Tractability for Multivariate Problems Generalize Tractability for Multivariate Problems Part II: Linear Tensor Prouct Problems, Linear Information, an Unrestricte Tractability Michael Gnewuch Department of Computer Science, University of Kiel,

More information

On the Cauchy Problem for Von Neumann-Landau Wave Equation

On the Cauchy Problem for Von Neumann-Landau Wave Equation Journal of Applie Mathematics an Physics 4 4-3 Publishe Online December 4 in SciRes http://wwwscirporg/journal/jamp http://xoiorg/436/jamp4343 On the Cauchy Problem for Von Neumann-anau Wave Equation Chuangye

More information

Function Spaces. 1 Hilbert Spaces

Function Spaces. 1 Hilbert Spaces Function Spaces A function space is a set of functions F that has some structure. Often a nonparametric regression function or classifier is chosen to lie in some function space, where the assume structure

More information

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012 CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration

More information

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix

More information

ORDINARY DIFFERENTIAL EQUATIONS AND SINGULAR INTEGRALS. Gianluca Crippa

ORDINARY DIFFERENTIAL EQUATIONS AND SINGULAR INTEGRALS. Gianluca Crippa Manuscript submitte to AIMS Journals Volume X, Number 0X, XX 200X Website: http://aimsciences.org pp. X XX ORDINARY DIFFERENTIAL EQUATIONS AND SINGULAR INTEGRALS Gianluca Crippa Departement Mathematik

More information

Global Solutions to the Coupled Chemotaxis-Fluid Equations

Global Solutions to the Coupled Chemotaxis-Fluid Equations Global Solutions to the Couple Chemotaxis-Flui Equations Renjun Duan Johann Raon Institute for Computational an Applie Mathematics Austrian Acaemy of Sciences Altenbergerstrasse 69, A-44 Linz, Austria

More information

A new proof of the sharpness of the phase transition for Bernoulli percolation on Z d

A new proof of the sharpness of the phase transition for Bernoulli percolation on Z d A new proof of the sharpness of the phase transition for Bernoulli percolation on Z Hugo Duminil-Copin an Vincent Tassion October 8, 205 Abstract We provie a new proof of the sharpness of the phase transition

More information

Logarithmic spurious regressions

Logarithmic spurious regressions Logarithmic spurious regressions Robert M. e Jong Michigan State University February 5, 22 Abstract Spurious regressions, i.e. regressions in which an integrate process is regresse on another integrate

More information

ON THE SHIFT DIFFERENTIABILITY OF THE FLOW GENERATED BY A HYPERBOLIC SYSTEM OF CONSERVATION LAWS. Stefano Bianchini S.I.S.S.A. (Trieste) May 1999

ON THE SHIFT DIFFERENTIABILITY OF THE FLOW GENERATED BY A HYPERBOLIC SYSTEM OF CONSERVATION LAWS. Stefano Bianchini S.I.S.S.A. (Trieste) May 1999 ON THE SHIFT DIFFEENTIABILITY OF THE FLOW GENEATED BY A HYPEBOLIC SYSTEM OF CONSEVATION LAWS Stefano Bianchini S.I.S.S.A. (Trieste) May 999 Abstract. We consier the notion of shift tangent vector introuce

More information

The Exact Form and General Integrating Factors

The Exact Form and General Integrating Factors 7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily

More information

Interconnected Systems of Fliess Operators

Interconnected Systems of Fliess Operators Interconnecte Systems of Fliess Operators W. Steven Gray Yaqin Li Department of Electrical an Computer Engineering Ol Dominion University Norfolk, Virginia 23529 USA Abstract Given two analytic nonlinear

More information

Chapter 6: Energy-Momentum Tensors

Chapter 6: Energy-Momentum Tensors 49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.

More information

7.1 Support Vector Machine

7.1 Support Vector Machine 67577 Intro. to Machine Learning Fall semester, 006/7 Lecture 7: Support Vector Machines an Kernel Functions II Lecturer: Amnon Shashua Scribe: Amnon Shashua 7. Support Vector Machine We return now to

More information

Research Article Global and Blow-Up Solutions for Nonlinear Hyperbolic Equations with Initial-Boundary Conditions

Research Article Global and Blow-Up Solutions for Nonlinear Hyperbolic Equations with Initial-Boundary Conditions International Differential Equations Volume 24, Article ID 724837, 5 pages http://x.oi.org/.55/24/724837 Research Article Global an Blow-Up Solutions for Nonlinear Hyperbolic Equations with Initial-Bounary

More information

THE KELLER-SEGEL MODEL FOR CHEMOTAXIS WITH PREVENTION OF OVERCROWDING: LINEAR VS. NONLINEAR DIFFUSION

THE KELLER-SEGEL MODEL FOR CHEMOTAXIS WITH PREVENTION OF OVERCROWDING: LINEAR VS. NONLINEAR DIFFUSION THE KELLER-SEGEL MODEL FOR CHEMOTAXIS WITH PREVENTION OF OVERCROWDING: LINEAR VS. NONLINEAR DIFFUSION MARTIN BURGER, MARCO DI FRANCESCO, AND YASMIN DOLAK Abstract. The aim of this paper is to iscuss the

More information

An extension of Alexandrov s theorem on second derivatives of convex functions

An extension of Alexandrov s theorem on second derivatives of convex functions Avances in Mathematics 228 (211 2258 2267 www.elsevier.com/locate/aim An extension of Alexanrov s theorem on secon erivatives of convex functions Joseph H.G. Fu 1 Department of Mathematics, University

More information

Suppressing Chemotactic Blow-Up Through a Fast Splitting Scenario on the Plane

Suppressing Chemotactic Blow-Up Through a Fast Splitting Scenario on the Plane Arch. Rational Mech. Anal. Digital Object Ientifier (DOI) https://oi.org/1.17/s25-18-1336-z Suppressing Chemotactic Blow-Up Through a Fast Splitting Scenario on the Plane Siming He & Eitan Tamor Communicate

More information

SEMILINEAR DEGENERATE PARABOLIC SYSTEMS AND DISTRIBUTED CAPACITANCE MODELS. Brooke L. Hollingsworth and R.E. Showalter

SEMILINEAR DEGENERATE PARABOLIC SYSTEMS AND DISTRIBUTED CAPACITANCE MODELS. Brooke L. Hollingsworth and R.E. Showalter DISCRETE AND CONTINUOUS DYNAMICAL SYSTEMS Volume, Number, January 995 pp. 59 76 SEMILINEAR DEGENERATE PARABOLIC SYSTEMS AND DISTRIBUTED CAPACITANCE MODELS Brooke L. Hollingsworth an R.E. Showalter Department

More information

Brooke L. Hollingsworth and R. E. Showalter Department of Mathematics The University of Texas at Austin Austin, TX USA

Brooke L. Hollingsworth and R. E. Showalter Department of Mathematics The University of Texas at Austin Austin, TX USA SEMILINEAR DEGENERATE PARABOLIC SYSTEMS AND DISTRIBUTED CAPACITANCE MODELS Brooke L. Hollingsworth an R. E. Showalter Department of Mathematics The University of Texas at Austin Austin, TX 7872 USA Abstract.

More information

SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. where L is some constant, usually called the Lipschitz constant. An example is

SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. where L is some constant, usually called the Lipschitz constant. An example is SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. Uniqueness for solutions of ifferential equations. We consier the system of ifferential equations given by x = v( x), () t with a given initial conition

More information

A COMBUSTION MODEL WITH UNBOUNDED THERMAL CONDUCTIVITY AND REACTANT DIFFUSIVITY IN NON-SMOOTH DOMAINS

A COMBUSTION MODEL WITH UNBOUNDED THERMAL CONDUCTIVITY AND REACTANT DIFFUSIVITY IN NON-SMOOTH DOMAINS Electronic Journal of Differential Equations, Vol. 2929, No. 6, pp. 1 14. ISSN: 172-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu A COMBUSTION MODEL WITH UNBOUNDED

More information

REAL ANALYSIS I HOMEWORK 5

REAL ANALYSIS I HOMEWORK 5 REAL ANALYSIS I HOMEWORK 5 CİHAN BAHRAN The questions are from Stein an Shakarchi s text, Chapter 3. 1. Suppose ϕ is an integrable function on R with R ϕ(x)x = 1. Let K δ(x) = δ ϕ(x/δ), δ > 0. (a) Prove

More information

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.

More information

Stable and compact finite difference schemes

Stable and compact finite difference schemes Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long

More information

Dissipative numerical methods for the Hunter-Saxton equation

Dissipative numerical methods for the Hunter-Saxton equation Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary

More information

LECTURE NOTES ON DVORETZKY S THEOREM

LECTURE NOTES ON DVORETZKY S THEOREM LECTURE NOTES ON DVORETZKY S THEOREM STEVEN HEILMAN Abstract. We present the first half of the paper [S]. In particular, the results below, unless otherwise state, shoul be attribute to G. Schechtman.

More information

On some parabolic systems arising from a nuclear reactor model

On some parabolic systems arising from a nuclear reactor model On some parabolic systems arising from a nuclear reactor moel Kosuke Kita Grauate School of Avance Science an Engineering, Wasea University Introuction NR We stuy the following initial-bounary value problem

More information

Sturm-Liouville Theory

Sturm-Liouville Theory LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

Acute sets in Euclidean spaces

Acute sets in Euclidean spaces Acute sets in Eucliean spaces Viktor Harangi April, 011 Abstract A finite set H in R is calle an acute set if any angle etermine by three points of H is acute. We examine the maximal carinality α() of

More information

ON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM

ON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract

More information

Second order differentiation formula on RCD(K, N) spaces

Second order differentiation formula on RCD(K, N) spaces Secon orer ifferentiation formula on RCD(K, N) spaces Nicola Gigli Luca Tamanini February 8, 018 Abstract We prove the secon orer ifferentiation formula along geoesics in finite-imensional RCD(K, N) spaces.

More information

Kinematics of Rotations: A Summary

Kinematics of Rotations: A Summary A Kinematics of Rotations: A Summary The purpose of this appenix is to outline proofs of some results in the realm of kinematics of rotations that were invoke in the preceing chapters. Further etails are

More information

Euler Equations: derivation, basic invariants and formulae

Euler Equations: derivation, basic invariants and formulae Euler Equations: erivation, basic invariants an formulae Mat 529, Lesson 1. 1 Derivation The incompressible Euler equations are couple with t u + u u + p = 0, (1) u = 0. (2) The unknown variable is the

More information

Iterated Point-Line Configurations Grow Doubly-Exponentially

Iterated Point-Line Configurations Grow Doubly-Exponentially Iterate Point-Line Configurations Grow Doubly-Exponentially Joshua Cooper an Mark Walters July 9, 008 Abstract Begin with a set of four points in the real plane in general position. A to this collection

More information

Lower Bounds for an Integral Involving Fractional Laplacians and the Generalized Navier-Stokes Equations in Besov Spaces

Lower Bounds for an Integral Involving Fractional Laplacians and the Generalized Navier-Stokes Equations in Besov Spaces Commun. Math. Phys. 263, 83 831 (25) Digital Obect Ientifier (DOI) 1.17/s22-5-1483-6 Communications in Mathematical Physics Lower Bouns for an Integral Involving Fractional Laplacians an the Generalize

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

Analysis of a penalty method

Analysis of a penalty method Analysis of a penalty metho Qingshan Chen Department of Scientific Computing Floria State University Tallahassee, FL 3236 Email: qchen3@fsu.eu Youngjoon Hong Institute for Scientific Computing an Applie

More information

A Sketch of Menshikov s Theorem

A Sketch of Menshikov s Theorem A Sketch of Menshikov s Theorem Thomas Bao March 14, 2010 Abstract Let Λ be an infinite, locally finite oriente multi-graph with C Λ finite an strongly connecte, an let p

More information

WELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL PRESSURE IN SOBOLEV SPACES

WELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL PRESSURE IN SOBOLEV SPACES Electronic Journal of Differential Equations, Vol. 017 (017), No. 38, pp. 1 7. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu WELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL

More information

HYPOCOERCIVITY WITHOUT CONFINEMENT. 1. Introduction

HYPOCOERCIVITY WITHOUT CONFINEMENT. 1. Introduction HYPOCOERCIVITY WITHOUT CONFINEMENT EMERIC BOUIN, JEAN DOLBEAULT, STÉPHANE MISCHLER, CLÉMENT MOUHOT, AND CHRISTIAN SCHMEISER Abstract. Hypocoercivity methos are extene to linear kinetic equations with mass

More information

Lower Bounds for the Smoothed Number of Pareto optimal Solutions

Lower Bounds for the Smoothed Number of Pareto optimal Solutions Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.

More information

Robustness and Perturbations of Minimal Bases

Robustness and Perturbations of Minimal Bases Robustness an Perturbations of Minimal Bases Paul Van Dooren an Froilán M Dopico December 9, 2016 Abstract Polynomial minimal bases of rational vector subspaces are a classical concept that plays an important

More information

GLOBAL SOLUTIONS FOR 2D COUPLED BURGERS-COMPLEX-GINZBURG-LANDAU EQUATIONS

GLOBAL SOLUTIONS FOR 2D COUPLED BURGERS-COMPLEX-GINZBURG-LANDAU EQUATIONS Electronic Journal of Differential Equations, Vol. 015 015), No. 99, pp. 1 14. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu GLOBAL SOLUTIONS FOR D COUPLED

More information

On the number of isolated eigenvalues of a pair of particles in a quantum wire

On the number of isolated eigenvalues of a pair of particles in a quantum wire On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in

More information

. ISSN (print), (online) International Journal of Nonlinear Science Vol.6(2008) No.3,pp

. ISSN (print), (online) International Journal of Nonlinear Science Vol.6(2008) No.3,pp . ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.6(8) No.3,pp.195-1 A Bouneness Criterion for Fourth Orer Nonlinear Orinary Differential Equations with Delay

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

Gradient flow of the Chapman-Rubinstein-Schatzman model for signed vortices

Gradient flow of the Chapman-Rubinstein-Schatzman model for signed vortices Graient flow of the Chapman-Rubinstein-Schatzman moel for signe vortices Luigi Ambrosio, Eoaro Mainini an Sylvia Serfaty Deicate to the memory of Michelle Schatzman (1949-2010) Abstract We continue the

More information

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS MARK SCHACHNER Abstract. When consiere as an algebraic space, the set of arithmetic functions equippe with the operations of pointwise aition an

More information

International Journal of Pure and Applied Mathematics Volume 35 No , ON PYTHAGOREAN QUADRUPLES Edray Goins 1, Alain Togbé 2

International Journal of Pure and Applied Mathematics Volume 35 No , ON PYTHAGOREAN QUADRUPLES Edray Goins 1, Alain Togbé 2 International Journal of Pure an Applie Mathematics Volume 35 No. 3 007, 365-374 ON PYTHAGOREAN QUADRUPLES Eray Goins 1, Alain Togbé 1 Department of Mathematics Purue University 150 North University Street,

More information

Calculus and optimization

Calculus and optimization Calculus an optimization These notes essentially correspon to mathematical appenix 2 in the text. 1 Functions of a single variable Now that we have e ne functions we turn our attention to calculus. A function

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

arxiv: v1 [physics.flu-dyn] 8 May 2014

arxiv: v1 [physics.flu-dyn] 8 May 2014 Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa

More information

NONLINEAR MARKOV SEMIGROUPS AND INTERACTING LÉVY TYPE PROCESSES

NONLINEAR MARKOV SEMIGROUPS AND INTERACTING LÉVY TYPE PROCESSES NONLINEAR MARKOV SEMIGROUPS AND INTERACTING LÉVY TYPE PROCESSES Vassili N. Kolokoltsov November 21, 2006 To appear in Journal of Statistical Physics. Abstract Semigroups of positivity preserving linear

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

CLARK-OCONE FORMULA BY THE S-TRANSFORM ON THE POISSON WHITE NOISE SPACE

CLARK-OCONE FORMULA BY THE S-TRANSFORM ON THE POISSON WHITE NOISE SPACE CLARK-OCONE FORMULA BY THE S-TRANSFORM ON THE POISSON WHITE NOISE SPACE YUH-JIA LEE*, NICOLAS PRIVAULT, AND HSIN-HUNG SHIH* Abstract. Given ϕ a square-integrable Poisson white noise functionals we show

More information

Calculus of Variations

Calculus of Variations Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,

More information

Introduction to variational calculus: Lecture notes 1

Introduction to variational calculus: Lecture notes 1 October 10, 2006 Introuction to variational calculus: Lecture notes 1 Ewin Langmann Mathematical Physics, KTH Physics, AlbaNova, SE-106 91 Stockholm, Sween Abstract I give an informal summary of variational

More information

Characterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations

Characterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations Characterizing Real-Value Multivariate Complex Polynomials an Their Symmetric Tensor Representations Bo JIANG Zhening LI Shuzhong ZHANG December 31, 2014 Abstract In this paper we stuy multivariate polynomial

More information