YOMEDIA
ADSENSE
Dirichlet cauchy problem for second order parabolic equations in domains with edges
4
lượt xem 1
download
lượt xem 1
download
Download
Vui lòng tải xuống để xem tài liệu đầy đủ
In this paper, we study the first initial boundary value problem for second-order parabolic equations in cylinders with piecewise smooth base. Some results on the unique existence and on the smoothness with respect to time of the solution are given.
AMBIENT/
Chủ đề:
Bình luận(0) Đăng nhập để gửi bình luận!
Nội dung Text: Dirichlet cauchy problem for second order parabolic equations in domains with edges
- JOURNAL OF SCIENCE OF HNUE Natural Sci., 2011, Vol. 56, No. 3, pp. 13-24 DIRICHLET-CAUCHY PROBLEM FOR SECOND-ORDER PARABOLIC EQUATIONS IN DOMAINS WITH EDGES Do Van Loi Hong Duc University E-mail: 37loilinh@gmail.com Abstract. In this paper, we study the first initial boundary value prob- lem for second-order parabolic equations in cylinders with piecewise smooth base. Some results on the unique existence and on the smoothness with respect to time of the solution are given. Keywords: Domains with edges, generalized solution, regularity. 1. Introduction Initial boundary value problems for non-stationary systems and equations in cylinders with non-smooth bases have been investigated in [2-8], where some impor- tant results on the unique existence, smoothness and asymptotic of the solution for the problems in Sobolev spaces were given. The initial boundary value problems for parabolic equations in a cylinder with base containing conical points [5] or in a dihe- dral angle [8] have been studied in Sobolev spaces with weights. In the present paper, we consider the first initial boundary value problem for second-order parabolic equa- tions in domains with edges. We study the existence, uniqueness and smoothness with respect to time of the generalized solution for this particular problem. The paper is organized as follows. In the second section we define Dirichlet- Cauchy problem for second-order parabolic equations in domains with edges. In the 3rd section we study the solvability of the problem. The last section is intended to regularity with respect to time of the Generani Zaisedzed solution. 2. Formulation of the problem Let Ω be a bounded domain in Rn , n ≥ 2. Its boundary ∂Ω is a piecewise smooth surface consisting of finitely many (n − 1)-dimensional smooth surfaces Γi . We assume that the surface Γi intersects only the surfaces Γi−1 , Γi+1 along smooth (n − 2)-Dimensional manifolds li−1 , li+1 . Without loss of generality we shall deal 13
- Do Van Loi explicitly with the case when ∂Ω consists two surfaces Γ1 , Γ2 which intersects along a manifold l0 . Assume that in a neighbourhood of each point of l0 the set Ω is diffeo- morphic to a dihedral angle. For any point P ∈ l0 two half-spaces T1 (P ), T2 (P ) tangent to Ω and a two-dimensional plane π(P ) normal to l0 are defined. We denote by ν(P ) the angle in the plane π(P ) (on the side of Ω) bounded by the rays T1 (P ) ∩ π(P ), T2 (P ) ∩ π(P ) and by β(P ) the aperture of this angle. Denote QT = Ω × (0, T ), ST = ∂Ω × (0, T ) for each T : 0 < T ≤ ∞. For each multi-index p = (p1 , . . . , pn ) ∈ Nn , |p| = p1 + · · · + pn , the symbol D p u = ∂ |p| u/∂xp11 ...∂xpnn denotes the generalized derivative of order α with respect to x = (x1 , ..., xn ); utk ∂ k u/∂tk is the generalized derivative of order k with respect to t. Throughout the paper we need the following functional spaces (see [1]): H l (Ω) is the space consisting of all functions u(x) defined on Ω such that 21 X Z kukH m (Ω) = |D p u|2 dx < +∞. |p|≤m Ω o ˚l (Ω) is the completion of C ∞ (Ω) in the norm of the space H l (Ω). H Let γ be a positive number. We denote the following spaces: H l,k (QT , γ) is the space consisting of all functions u(x, t) defined on QT such that Z X k X kuk2H m,k (QT ,γ) = |D p u|2 + |utj |2 e−γt dxdt < +∞, QT |p|≤m j=1 Z X kuk2H m,0 (QT ,γ) = |D p u|2 e−γt dxdt < +∞. QT |p|≤m ˚l,k (QT , γ) is the closure in H l,k (QT , γ) of the set of all infinitely differentiable H functions on QT which vanish near ST . L2 (QT , γ) is the space of functions u(x, t) defined on QT with the normal Z 12 2 −γt kukL2 (QT ,γ) = |u| e dxdt . QT Now we introduce the differential operator 14
- Dirichlet-Cauchy problem for second-order parabolic equations in domains with edges ∂ ∂u X X n n ∂u L(x, t, D) = − aij (x, t) + bi (x, t) + c(x, t)u, i,j=1 ∂xi ∂xj i=1 ∂xi where aij (x, t), bi (x, t) and c(x, t) are bounded real-valued functions on QT and aij (., t) belongs to C 1 (Ω), t ∈ (0, T ). Moreover, suppose that aij , i, j = 1, ..., n, are continuous in x ∈ Ω uniformly with respect to t ∈ (0, T ) and n X aij (x, t)ξi ξj ≥ µ0 |ξ|2 (2.1) i,j=1 for all ξ ∈ Rn \{0} and (x, t) ∈ QT , where µ0 = const > 0. We consider the following problem in the cylinder QT : ut + L(x, t, D)u = f (2.2) with the initial conditions u|t=0 = 0, (2.3) and boundary conditions
- u
- = 0. (2.4) ST Denote n Z X X n Z Z ∂u ∂v ∂u B(u, v; t) = aij (., t) dx + bi (., t) vdx + c(., t)uvdx. i,j=1 Ω ∂xj ∂xi i=1 ∂xi Ω Ω Applying condition (2.1) and similar arguments as the proof of Guarding in- equality it follows that B(u, u; t) ≥ µ0 kuk2H 1(Ω) − λ0 kuk2L2(Ω) (2.5) o for all u(x, t) ∈ H 1,1 (QT , γ), where µ0 = const > 0, λ0 = const ≥ 0. Without loss of generality we shall deal explicitly with the case when λ0 = 0, since by a substitution v = eλ0 t u problem (2.2) - (2.4) can be transformed to a problem with constant λ0 = 0. The function u(x, t) is called a generalized solution in the space H 1,1 (γ, QT ) of o problem (2.2) - (2.4) if and only if u(x, t) ∈ H 1,1(QT , γ), u(x, 0) = 0 and the equality hut , vi + B(u, v; t) = hf, vi (2.6) o holds for all v ∈ H 1 (Ω). 15
- Do Van Loi 3. The well-posedness Theorem 3.1. Let f ∈ L2 (QT , γ0 ), γ0 > 0, and suppose that the coefficients of the operator L satisfy sup{|aij |, |bi |, |c| : i, j = 1, . . . , n; (x, t) ∈ QT } ≤ µ. Then for each γ > γ0 , problem (2.2) - (2.4) has unique generalized solution u in the ˚1,1 (QT , γ) and the following estimate holds space H kuk2H 1,1 (QT ,γ) ≤ Ckf k2L2 (QT ,γ0 ), (3.1) where C is a constant independent of u and f . This solution depends continuously on f . Proof. Firstly, we will prove the existence by Galerkin’s approximating method. Let {ωk (x)}∞ ˚1 k=1 be an orthogonal basis of H (Ω) which is orthonormal in L2 (Ω). Puting XN N u (x, t) = CkN (t)ωk (x), k=1 where CkN (t), t ∈ [0, T ), k = 1, ..., N, are the solution of the system of the following ordinary differential equations: (uN N t , ωk ) + B(u , ωk ; t) = (f, ωk ), t ∈ [0, T ), k = 1, ..., N, (3.2) with the initial conditions CkN (0) = 0, k = 1, ..., N, (3.3) here (., .) is the inner products in the space L2 (Ω). Let us multiply (3.2) by CkN (t), then take the sum with respect to k from 1 to N, we arrive at (uN N N N N t , u ) + B(u , u ; t) = (f, u ), t ∈ [0, T ). Now adding this equality to its complex conjugate, we get d kuN k2L2 (Ω) + 2ReB(uN , uN ; t) = 2Re(f, uN ). (3.4) dt From (2.5) with λ0 = 0, we obtain ReB(uN , uN ; t) ≥ µ0 kuN k2H 1 (Ω) . 16
- Dirichlet-Cauchy problem for second-order parabolic equations in domains with edges On the other hand, by the Cauchy inequality, for an arbitrary positive number ε, we have 2|(f, uN )| ≤ 2kf kL2 (Ω) kuN kL2 (Ω) ≤ Ckf k2L2 (Ω) + εkuN k2L2 (Ω) , where C = C(ε) is a constant independent of uN , f and t. Combining the estimations above, we get from (3.4) that d kuN (., t)k2L2 (Ω) + 2µkuN (., t)k2H 1 (Ω) ≤ Ckf (., t)k2L2 (Ω) + εkuN (., t)k2L2 (Ω) (3.5) dt for a.e. t ∈ [0, T ). Now write η(t) := kuN (., t)k2L2 (Ω) ; ξ(t) := kf (., t)k2L2 (Ω) , t ∈ [0, T ). Then (3.5) implies η ′ (t) ≤ ε.η(t) + ξ(t), for a.e. t ∈ [0, T ). Thus the differential form of Gronwall-Belmann’s inequality yields the estimate Zt εt η(t) ≤ Ce ξ(s)ds), t ∈ [0, T ). (3.6) 0 We obtain from (3.6) the estimate Zt kuN (., t)k2L2 (Ω) ≤ Ce(ε+γ0 )t e−γ0 s kf k2L2 (Ω) ds ≤ Ce(γ0 +ε)t kf k2L2 (QT ,γ0 ). 0 Now multiplying both sides of the inequality above by e−γt , γ > γ0 + ε, then inte- grating them with respect to t from 0 to T , we obtain kuN k2L2 (QT ,γ) ≤ Ckf k2L2 (QT ,γ0 ) . (3.7) Returning once more to inequality (3.5), we multiply both sides of the inequality above by e−γt , then integrating them with respect to t from 0 to τ, τ ∈ (0, T ), we obtain Zτ Zτ −γt d N 2 e ku kL2 (Ω) dt + 2µ e−γt kuN k2H 1 (Ω) dt dt 0 0 ≤ C kf k2L2 (QT ,γ0 ) + kuN k2L2 (QT ,γ) . 17
- Do Van Loi Noting that Zτ Zτ Zτ −γt d N 2 d −γt N 2 e ku kL2 (Ω) dt = e ku kL2 (Ω) dt + γ e−γt kuN k2L2 (Ω) dt dt dt 0 0 0 Zτ = e−γτ kuN (x, τ )k2L2 (Ω) + γ e−γt kuN k2L2 (Ω) dt ≥ 0. 0 We employ the inequalities above to find Zτ 2µ e−γt kuN k2H 1 (Ω) dt ≤ Ckf k2L2 (QT ,γ0 ) , ∀τ ∈ (0, T ). (3.8) 0 Since the right-hand side of (3.8) is independent of τ , we get kuN k2H 1,0 (QT ,γ) ≤ Ckf k2L2(QT ,γ0 ), (3.9) where C is a constant independent of u, f and N. ˚1 (Ω), with kvk2 1 ≤ 1 and write v = v 1 +v 2 where v 1 ∈ span{ωk }N Fix any v ∈ H H (Ω) k=1 ⊥ and (v 2 , ωk ) = 0, k = 1, ..., N, (v 2 ∈ span{ωk }N 1 k=1 ). We have kv kH 1 (Ω) ≤ kvkH 1 (Ω) ≤ 1. From (3.2), we get 1 1 1 (uN N t , v ) + B(u , v ; t) = (f, v ), for a.e. t ∈ [0, T ). P N From uN (x, t) = CkN (t)ωk , we can see that k=1 1 1 1 (uN N N t , v) = (ut , v ) = (f, v ) − B(u , v ; t), consequently, 2 N 2 |(uN t , v)| ≤ C kf k L2 (Ω) + ku k H (Ω) . 1 Since the inequality above holds for all v ∈ H ˚1 (Ω), kvkH 1 (Ω) ≤ 1, we receive the following inequality 2 2 N 2 kuN k t L2 (Ω) ≤ C kf k L2 (Ω) + ku k 1 H (Ω) (3.10) Multiplying (3.10) by e−γt , γ > γ0 + ε, then integrating them with respect to t from 0 to T , and by using (3.9), we obtain 2 2 kuN t kL2 (QT ,γ) ≤ Ckf kL2 (QT ,γ0 ) . (3.11) 18
- Dirichlet-Cauchy problem for second-order parabolic equations in domains with edges Combining (3.9) and (3.11), we arrive at kuN k2H 1,1 (QT ,γ) ≤ Ckf k2L2 (QT ,γ0 ), (3.12) here C is a absolute constant. According to inequality (3.12), by standard weakly convergent arguments, we can conclude that the sequence {uN }∞ N =1 possesses a subsequence convergent to a func- ˚ 1,1 tion u ∈ H (QT , γ), which is a generalized solution of the problem (2.2)-(2.4). Moreover, it follows from (3.12) that inequality (3.1) holds. Finally, we will prove the uniqueness of generalized solution. It suffices to check that the only the generalized solution of problem (2.2)-(2.4) with f ≡ 0 is u ≡ 0. (3.13) By setting v = u(., t) in identity (2.6) (for f ≡ 0) and adding it to its complex conjugate, we get d (ku(., t)k2) + 2ReB(u, u; t) = 0. dt Since (2.5), we have d (kuk2L2 (Ω) ) + 2µkuk2H 1(Ω) ≤ 0, for a.e. t ∈ [0, T ). dt From this inequality and Gronwall-Belmann’s inequality, it implies (3.13). By esti- mate (3.12), we also get that the solution u depends continuously on f . 4. The smoothness with respect to time variable of the generalized solution In this section, we will show that u is smooth with respect to time variable, provided the coefficients of the operator L and the right-hand side f of system (2.2), are smooth with respect to time variable. We also show that the smoothness with respect to time variable of u is independent of the base Ω of cylinders. To simplify notation, we write ∂ k aij ∂ k bi ∂k c aijtk = , b k = , c k = ∂tk it ∂tk t ∂tk and n Z X X n Z Z ∂u ∂v ∂u Btk (u, v; t) = aijtk (., t) dx + bitk (., t) vdx + c(., t)tk uvdx. i,j=1 Ω ∂xj ∂xi i=1 ∂xi Ω Ω 19
- Do Van Loi Theorem 4.1. Let h ∈ N∗ , and we assume that n o (i) sup |aijtk |, |bitk |, |ctk | : i, j = 1, . . . , n; (x, t) ∈ QT , k ≤ h ≤ µ, (ii) ftk ∈ L2 (QT , γk ), k ≤ h; ftk (x, 0) = 0, 0 ≤ k ≤ h − 1. Then for an arbitrary real number γ satisfying γ > γ0 , the generalized solution u∈H ˚1,1 (QT , γ) of problem (2.2)-(2.4) has derivatives with respect to t up to order h with utk ∈ H˚1,1(QT , γ + ε), k = 0, ..., h, and we have the estimate h X kuth k2H 1,1 (QT ,γ+ε) ≤ C kftj k2L2 (QT ,γ0 ), (4.1) j=0 where C is a constant independent of u, f and ε > 0 is small. Proof. From the assumptions on the coefficients of operator L and the function f , it implies that the solution {CkN }N k=1 of problem (3.2)-(3.3) has derivatives with respect to t up to order h + 1. We will prove by induction that h X 2 kuN th kH 1,0 (QT ,γ+ε) ≤ C kftj k2L2 (QT ,γ0 ) . (4.2) j=0 Firstly, we differential h times both sides of (3.2) with respect to t to finding the following equality: h X N h (uth+1 , ωk ) + Bth−l (uN tl , ωk ; t) = (fth , ωk ), k = 1, ..., N (4.3) l l=0 From these equalities together with the initial (3.3) and the assumption (ii), we can show by induction on h that uN tk |t=0 = 0 for k = 0, . . . , h. (4.4) dh CkN (t) We multiply equality (4.3) by and sum k = 1, . . . , N, to discover dth Xh N N h (uth+1 , uth ) + Bth−j (uN N N tj , uth ; t) = (fth , uth ). j=0 j Adding this equality to its complex conjugate, we get Xh d N 2 h ku (., t)kL2 (Ω) + 2 Re Bth−j (uN N N tj , uth ; t) = 2 Re(fth , uth ). (4.5) dt j=0 j 20
- Dirichlet-Cauchy problem for second-order parabolic equations in domains with edges According to Theorem 3.1, it implies that inequalities (4.2) holds for h = 0. Assume now inequalities (4.2) is valid for k = h − 1, we will show it to be true for k = h. Returning to equality (4.5), the second term in the left-hand side of (4.5) is written in the following form: Xh h 2Re Bth−j (uN N tj , uth ; t) = j=0 j Xh−1 N N h = 2ReB(uth , uth ; t) + 2Re Bth−j (uN N tj , uth ; t). j=0 j Hence, from (4.5) we have h−1 X d N 2 N N N h ku (., t)kL2 (Ω) +2ReB(uth , uth ; t) = 2 Re(fth , uth )−2 Re Bth−j (uN N tj , uth ; t). dt j=0 j (4.6) For shortly, we denote by I, II, the terms from the first to the second, respectively, of the right-hand side of (4.6). By using the asumption (i) and the Cauchy inequality, we obtain the following estimates: (I) ≤C(ε1 )kfth k2L2 (Ω) + ε1 kuth k2L2 (Ω) . h−1 X 2 N 2 (II) ≤C(ε2 ) kuN tj kH 1 (Ω) + ε2 µnhkuth kH 1 (Ω) . j=0 Employing the estimates above, we get from (4.6) that d kuN (., t)k2L2 (Ω) + 2ReB(uN N 2 2 th , uth ; t) ≤ C(ε1 )kfth kL2 (Ω) + ε1 kuth kL2 (Ω) dt h−1 X 2 N 2 + C(ε2 ) kuNtj kH 1 (Ω) + ε2 µnhkuth kH 1 (Ω) . j=0 (4.7) By using again (2.5), we obtain from (4.7) the estimate d kuN (., t)k2L2 (Ω) + (2µ0 − ε2 µnh)kuN 2 2 2 th kH 1 (Ω) ≤ ε1 kuth kL2 (Ω) + C(ε1 )kfth kL2 (Ω) dt h−1 X 2 + C(ε2 ) kuNtj kH 1 (Ω) . j=0 (4.8) 21
- Do Van Loi ε1 By let ε2 > 0 such that (2µ0 − ε2 µnh) < µ0 , and let ε = min{µ0 ;1} . From (4.8), we get X h−1 d kuN (., t)k2L2 (Ω) + kuNh k 2 t H (Ω)1 ≤ εku h k 2 t L2 (Ω) + C kf h k 2 t L2 (Ω) + kuN 2 tj kH 1 (Ω) . dt j=0 (4.9) Now write h−1 X 2 2 2 η(t) := kuN th (., t)kL2 (Ω) ; ξ(t) := kfth kL2 (Ω) + kuN tj kH 1 (Ω) , t ∈ [0, T ). j=0 Then (4.9) implies η ′ (t) ≤ ε.η(t) + Cξ(t), for a.e. t ∈ [0, T ). Thus the differential form of Gronwall-Belmann’s inequality yields the estimate Zt εt η(t) ≤ Ce ξ(s)ds), t ∈ [0, T ). (4.10) 0 We obtain from (4.10) the estimate Zt h−1 X 2 2 2 kuN th (., t)kL2 (Ω) ≤ Ce εt kfth (., s)kL2 (Ω) + kuN tj (., s)kH 1 (Ω) ds 0 j=0 h−1 X (ε+γ0 )t ≤ Ce kfth k2L2 (QT ,γ0 ) + Ce (ε+γ)t kuN 2 tj kH 1,0 (QT ,γ+ε) j=0 Multiplying both sides of the inequality above by e−(γ+ε)t , then integrating them with respect to t from 0 to T and using the induction assumptions, we arrive at h X 2 kuN th kL2 (QT ,γ+ε) ≤C kftj k2L2 (QT ,γ0 ) . (4.11) j=0 Returning once more to inequality (4.9), we multiply both sides of the inequality above by e−(γ+ε)t , then integrating them with respect to t from 0 to τ, τ ∈ (0, T ), we obtain Zτ Zτ d e−(γ+ε)t kuNh k2 dt + 2µ0 e−(γ+ε)t kuN 2 th kH 1 (Ω) dt dt t L2 (Ω) 0 0 h X ≤C kftj k2L2 (QT ,γ0 ) . j=0 22
- Dirichlet-Cauchy problem for second-order parabolic equations in domains with edges Noting that Zτ Zτ Zτ −(γ+ε)t d N 2 d −(γ+ε)t N 2 e kuth kL2 (Ω) dt = e kuth kL2 (Ω) dt + (γ + ε) e−γt kuN 2 th kL2 (Ω) dt dt dt 0 0 0 Zτ = e−(γ+ε)τ kuN 2 th (x, τ )kL2 (Ω) + (γ + ε) e−γt kuN 2 th kL2 (Ω) dt ≥ 0. 0 We employ the inequalities above to find Zτ h X −(γ+ε)t 2 e kuN th kH 1 (Ω) dt ≤C kftj k2L2 (QT ,γ0 ) , ∀τ ∈ (0, T ). (4.12) 0 j=0 Since the right-hand side of (4.12) is independent of τ , we get h X 2 kuN th kH 1,0 (QT ,γ+ε) ≤C kftj k2L2 (QT ,γ0 ), (4.13) j=0 where C is a constant independent of u, f and N. By the arguments analogous to the proof of Theorem 3.1, and using the asumption (i) again, we obtain the estimate h X 2 kuN tk kL2 (QT ,γ+ε) ≤ C kftj k2L2 (QT ,γ0 ) , (4.14) j=0 for k = 1, . . . , h + 1. Combining (4.13) and (4.14), we conclude that h X 2 kuN th kH 1,1 (QT ,γ+ε) ≤C kftj k2L2 (QT ,γ0 ) . (4.15) j=0 From this inequality, by again standard weakly convergent arguments, we can con- ∞ clude that the sequence {uN tk }N =1 possesses a subsequence convergent to a function ˚1,1 (QT , γ + ε); moreover, u(k) is k th generalized derivative in t of generalized u(k) ∈ H solution u of problem (2.2) - (2.4). The estimate (4.1) follows from (4.15) by passing the weak convergences. Acknowledgement. This work was supported by National Foundation for Science and Technology Development (NAFOSTED), Vietnam, under project No. 101.01.58.09. 23
- Do Van Loi REFERENCES [1] R. A. Adams. Sobolev Spaces, 1975. Academic Press. [2] A. Kokatov and B. A. Plamenevssky, 2005. On the asymptotic on solutions to the Neumann problem for hyperbolic systems in domain with conical point. English transl., St. Peterburg Math. J, 16, No. 3, pp. 477-506. [3] Nguyen Manh Hung and Cung The Anh, 2010. Asymtotic expansions of solutions of the first initial boundary value problem for the Schr¨odinger system near conical points of the boundary. Differentsial’nye Uravneniya, Vol. 46, No. 2, pp. 285-289. [4] N. M. Hung, 1999. Asymptotic behaviour of solutions of the first boundary-value problem for strongly hyperbolic systems near a conical point at the boundary of the domain. Math. Sbornik, 19, pp. 103-126. [5] Nguyen Manh Hung and Nguyen Thanh Anh, 2008. Regularity of solutions of initial-boundary value problems for parabolic equations in domains with conical points. Journal of Differential Equations, Vol. 245, Issue 7, pp. 1801-1818. [6] Nguyen Manh Hung, Jen-Chih Yao, 2009. On the asymptotics of solutions of the first initial boundary value problem for hyperbolic systems in infinite cylinders with base containing conical points. Nonlinear Anal. Vol. 71, Issue 5-6, pp. 1620-1635. [7] Hung N. M.,Yao J. C., 2010. Cauchy-Neumann problem for second-order hyper- bolic equations in cylinder with base containing conical points. Applicable Analy- sis, Vol. 89, pp. 1293-1311. [8] Solonnikov V. A., 1983. On the solvability of classical initial-boundary value prob- lem for the heat equation in a dihedral angle. Zap. Nachn. Sem. Leningr. Otd. Math. Inst., 127, pp. 7-48. 24
Thêm tài liệu vào bộ sưu tập có sẵn:
Báo xấu
LAVA
AANETWORK
TRỢ GIÚP
HỖ TRỢ KHÁCH HÀNG
Chịu trách nhiệm nội dung:
Nguyễn Công Hà - Giám đốc Công ty TNHH TÀI LIỆU TRỰC TUYẾN VI NA
LIÊN HỆ
Địa chỉ: P402, 54A Nơ Trang Long, Phường 14, Q.Bình Thạnh, TP.HCM
Hotline: 093 303 0098
Email: support@tailieu.vn