ON LOCAL EXISTENCE OF SOLUTIONS OF THE FREE BOUNDARY PROBLEM FOR AN INCOMPRESSIBLE VISCOUS SELF-GRAVITATING FLUID MOTION

Size: px
Start display at page:

Download "ON LOCAL EXISTENCE OF SOLUTIONS OF THE FREE BOUNDARY PROBLEM FOR AN INCOMPRESSIBLE VISCOUS SELF-GRAVITATING FLUID MOTION"

Transcription

1 APPLICATIONES MATHEMATICAE 27,3(2), pp P.B.MUCHAandW.ZAJA CZKOWSKI(Waszawa) ON LOCAL EXISTENCE OF SOLUTIONS OF THE FREE BOUNDARY PROBLEM FOR AN INCOMPRESSIBLE VISCOUS SELF-GRAVITATING FLUID MOTION Abstact. The local-in-time existence of solutions of the fee bounday poblem fo an incompessible viscous self-gavitating fluid motion is poved. We show the existence of solutions with lowest possible egulaity fo this poblem such that u W ( T ) with > 3. The existence is poved by the method of successive appoximations whee the solvability of the Cauchy Neumann poblem fo the Stokes system is applied. We have to undeline that in the L p -appoach the Lagangian coodinates must be used. We ae looking fo solutions with lowest possible egulaity because this simplifies the poof and deceases the numbe of compatibility conditions. 1. Intoduction. In this pape we conside the motion of a viscous incompessible fluid in a bounded domain t R 3 with a fee bounday S t which is unde the self-gavitational foce. Let v = v(x,t) be the velocity of the fluid, p = p(x,t) the pessue, ν the constant viscosity coefficient and p the extenal pessue. Then the poblem is descibed by the following system: (1.1) v t +v v divt(v,p) = U in T, divv = in T, T(u,p) n = p n on S T, v t= = v in, t t= =, S t t= = S, v n = ϕ t / ϕ on S T, 2 Mathematics Subject Classification: 35Q3, 76D5. Key wods and phases: local existence, Navie Stokes equations, incompessible viscous baotopic self-gavitating fluid, shap egulaity, anisotopic Sobolev space. ReseachsuppotedbyPolishKBNGant2P3A3816. [319]

2 32 P.B.MuchaandW.Zaja czkowski whee T = t T t {t}, S T = t T S t {t}, ϕ(x,t) = descibes S t at least locally, n is the unit outwad vecto nomal to S t, n = ϕ/ ϕ, t is the domain at time t, S t = t, t T. Moeove, the dot denotes the scala poduct in R 3. By T = T(v,p) we denote the stess tenso of the fom (1.2) T(v,p) = {T ij } i,j=1,2,3 = { pδ ij +D ij (v)} i,j=1,2,3 whee (1.3) D(v) = {D ij (v)} i,j=1,2,3 = {ν(v i,xj +v j,xi )} i,j=1,2,3 is the velocity defomation tenso. Moeove, U( t,x,t) is the self-gavitational potential (1.4) U( t,x,t) = k t dy x y, whee k is the gavitation constant and some aguments of U ae omitted in evident cases. In view of the equation (1.1) 2 and the kinematic condition (1.1) 6 the total volume is conseved: (1.5) t = t dx = dx =. Let be given. Then we intoduce the Lagangian coodinates ξ as the initial data fo the following Cauchy poblem: x (1.6) t = v(x,t), x t= = ξ, ξ = (ξ 1,ξ 2,ξ 3 ). Integating (1.6), we obtain a tansfomation which connects the Euleian x and the Lagangian ξ coodinates, (1.7) x = x(ξ,t) ξ + t u(ξ,t ) x u (ξ,t), wheeu(ξ,t) = v(x u (ξ,t),t) and the index u in x u (ξ,t) will beomitted when no confusion can aise. Then fom (1.1) 6 we have t = {x R 3 : x = x(ξ,t), ξ } and S t = {x R 3 : x = x(ξ,t), ξ S = S = }. Ou aim is to pove the local-in-time existence of solutions to poblem (1.1) with lowest possible egulaity. Theefoe we apply the L p -appoach. The esult of the pape is the following theoem. Theoem 1.1. Let > 3, v W 2 2/ (), S W 2 1/. Then thee exists T > such that fo all T T thee exists a unique solution (u,p) of (1.1) such that u W ( T ), p W 1, ( T ) W ( S T )

3 Fee bounday poblem 321 and the following estimate holds: (1.8) u W ( T ) + p W 1, ( T ) + p W ( S T ) c(t) v W 2 2/ (). To pove Theoem 1.1 we need solvability of the Cauchy Neumann poblem fo the Stokes system fom [3]. To ecall the esult we fomulate the poblem u t divt(u,p) = F in T, (1.9) divu = G in T, n T(u,p) = H on S T, u t= = u in, whee T = [,T] and S T = S [,T]. Theoem 1.2 (see [3]). Let > 3, F L ( T ), G W 1, ( T ), G t divf = divb +A, A,B L ( T ), H W (S T ), u W 2 2/ (), S W 2 2/, and assume the compatibility conditions (1.1) divu = G(x,), n T(u,p ) S = H(x,), whee p = p t=. Then thee exists a unique solution (u,p) to poblem (1.9) such that u W ( T ), p W 1, ( T ) W (S T ) and the following estimate holds: (1.11) u W ( T ) + p W 1, ( T ) + p W (S T ) C(T)[ F L ( T )+ G W 1, ( T ) + B L ( T ) + A L ( T ) + H W (S T ) + u 2 2/ W () ], whee C(T) is a constant inceasing with T which does not depend on the solution (u, p). Poblem (1.1) without the self-gavitation foce is consideed in [4]. Moeove we ecall that the local existence of solutions to poblem (1.1) with suface tension is shown in [5]. 2. Notation. We neeheanisotopic Sobolev spaces W m,n (Q T ) whee m,n R + {}, 1 and Q T = Q (,T), with the nom

4 322 P.B.MuchaandW.Zaja czkowski (2.1) u W m,n (Q T ) = Q + u(x,t) dx m [ m ] Q m =[ m ] n [ n ] Q Q dx D m QQ D [n] D n x u(x,t) dx D m x u(x,t) Dx m u(x,t) dxdx x x s+( m [ m ]) t u(x,t) dx t u(x,t) D [n] t (x,t ), t t 1+(n [n]) whee s = dimq, [α] is the integal pat of α, D l x = l 1 x1... l s xs whee l = (l 1,...,l s ) is a multiindex. In the poof we will use the following esults. and then Poposition 2.1 (see [1]). Let u W m,n ( T ), m,n R +. If q κ = fo all ε (,1). 3 i=1 ( α i ) ( 1 q m + β ) 1 q n < 1 D β t Dα x u L q ( T ) ε 1 κ u W m,n ( T ) +cε κ u L ( T ) Poposition 2.2 (see [1, 2]). Let u W 2m,m ( T ), m R +. If 2m 1/ > then u = u ST is well defined as a function in W 2m 1/,m 1/(2) (S T ) and u 2m 1/,m 1/(2) W (S T ) c u W 2m,m ( T ). Poposition 2.3 (see [1, 2]). Let u W 2m 1/,m 1/(2) (S T ), m R +. If 2m 1/ > then thee exists a function ũ W 2m,m ( T ) such that ũ ST = u and the following estimate holds: ũ W 2m,m ( T ) c u W 2m 1/,m 1/(2) (S T ). In ou consideations we will use well known imbedding theoems fo Sobolev spaces. All constants ae denoted by the same lette c.

5 Fee bounday poblem Poof of Theoem 1.1. To pove local existence of solutions to poblem (1.1) we wite it in the Lagangian coodinates: (3.1) u t div u T u (u,q) = u U u in T, div u u = in T, n u T u (u,q) = p n u on S T, u t= = v on, whee u(ξ,t) = v(x(ξ,t),t), q(ξ,t) = p(x(ξ,t),t), u = ξ i,x ξi, T u (u,q) = D u (u) qi,d u (u) = ν{ξ k,xi u j,ξk +ξ k,xj u i,ξk } i,j=1,2,3, U u (ξ,t) = kj y(ξ,t)dξ x(ξ,t) y(ξ,t), whee J x(ξ,t) is the Jacobian of the tansfomation x = x(ξ,t), div u u = ξ k,xi u i,xk, n u (ξ,t) = n(x(ξ,t),t), I is the unit matix and the summation convention ove epeated indices is used. To pove the existence of solutions to (3.1) we use the following method of successive appoximations: (3.2) u m+1,t div um T um (u m+1,q m+1 ) = um U um in T, div um u m+1 = in T, n um T um (u m+1,q m+1 ) = p n um on S T, u m+1 t= = v on, whee m =,1,2,... and u m is teated as a given function. Assume u = and q =. To apply Theoem 1.2 we wite (3.2) in the fom u m+1,t divt(u m+1,q m+1 ) = div um T um (u m+1,q m+1 ) divt(u m+1,q m+1 )+ um U um in T, (3.3) divu m+1 = divu m+1 div um u m+1 in T, n T(u m+1,q m+1 ) = n T(u m+1,q m+1 ) n um T um (u m+1,q m+1 ) p n um on S T, u m+1 t= = v on, whee the opeatos without index contain deivatives with espect to ξ and n is the unit outwad vecto nomal to S.

6 324 P.B.MuchaandW.Zaja czkowski Fist we obtain a unifom bound fo the sequence {u m } m= detemined by (3.3). Lemma 3.1. Assume that S W 2 1/, v W 2 2/ (). Then (3.4) u m W ( T ) + q m W 1, ( T ) + q m W (S T ) if T is small enough. c( v W 2 2/ Poof. Applying Theoem 1.2 to poblem (3.3) yields (), S W 2 2/ (3.5) u m+1 W ( T ) + q m+1 W 1, ( T ) + q m+1 W (S T ) c divt(u m+1,q m+1 ) div um T um (u m+1,q m+1 ) L ( T ) +c um U um L ( T ) +c divu m+1 div um u m+1 W 1, ( T ) +c n T(u m+1,q m+1 ) n um T um (u m+1,q m+1 ) W (S T ) +c n um W +c ((I A (u m ))u m+1 ) t L ( T ), (S T ) +c v 2 2/ W () whee ((I A (u Ì m ))u m+1 ) t is teated as B fom Theoem 1.2 (Ã = ) and t A ij (u m ) = δ ij + u mi,ξ j dτ, A kl (u m) = A 1 lk (u m). Hee we note that div um u m+1 = A 1 kl ξ l u l m+1 = div ξ (A u m+1 ), which follows fom 3 k=1 ξ k A lk (u m )(ξ,t) =. All the above elations hold unde the assumption that div um 1 u m =. To continue the induction we need to have div um u m+1 =, but this is given by (3.3) 2. Now we estimate the paticula tems fom the Ì.h.s. of (3.5). Define a m = T ( 1)/ u m W ( T ), α t m(t) = {α ij (u m )} = { u mi,ξ j dτ}. Fo > 3 we have α m L () ca m. To estimate the fist tem on the.h.s. of (3.5) we calculate div um T um (u m+1,q m+1 ) divt(u m+1,q m+1 ) = {ν(ξ lxj ξ kxj x s x sξl δ σi +ξ lxj ξ kxi x s x sξl δ σj )u m+1σ,ξk +ν(ξ lxi ξ kxj δ jk δ jl )u m+1i,ξl ξ k +ν(ξ lxj ξ kξi δ lj δ ki )u m+1j,ξl ξ k (ξ lxj δ lj q m+1,ξl )}, whee the matix ξ,x depends on u m. )

7 Fee bounday poblem 325 Since x iξj = δ ij + Ì t u iξ j (τ)dτ = δ ij +α ij and ξ jxi is the invese matix to x iξj, we have ξ jxi = δ ij +φ ij (α), whee φ ij is a polynomial matix-valued function which contains tems of α and α 2 (α = {α ij }). Then ξ j,xi x k = φ ij,αs α s,ξσ ξ σ,xk, whee α s,ξσ = Ì t u,ξ s ξ σ (τ)dτ. Then we wite the fist tem of the.h.s. of (3.5) in the fom ψ 1 (α m )(I A(u m ))(u m+1,ξξ +q m+1,ξ )+ψ 2 (α m )A(u m ),ξ u m+1,ξ L ( T ) φ(a m )a m ( u m+1 W ( T ) + q m+1 W 1, ( T ) ), whee ψ i ae some functions with ψ i () and φ always denotes an inceasing positive function. We estimate the thid tem by the same quantity. The fouth tem can be expessed in the fom whee ψ 3 (α m )α m u m+1,ξ +ψ 4 (α m )α m q m+1 W (S T ) I Next we have L ( ( ψ 3 (α m )α m u m+1,ξ W (S T ) + ψ 4 (α m )α m q m+1 W (S T ) I +J, ψ 3 (α m )α m u m+1,ξ W 1()dτ ( + ψ 3 (α m )α m u m+1,ξ W 1/2 (,T) L+K. ψ 3 (α m )α m u m+1,ξ L () dτ ( ( ( ψ 3,αm (α m ) ψ 3 (α m ) t t u m,ξξ dτ α m u m+1,ξ u m,ξξ dτ u m+1,ξ L () L () ψ 3 (α m )α m u m+1,ξξ L () dτ L1 +L 2 +L 3 +L 4.

8 326 P.B.MuchaandW.Zaja czkowski Continuing, we have L 1 φ(a m (T))a m (T) u m+1 W L 2 +L 3 φ(a m (T))a m (T) Next we examine ( K dξ ( t ( T ), u m,ξξ dτ u m+1,ξ φ(a m (T))a m (T)a 2 m(t) u m+1 W ( T ), L 4 φ(a m (T))a m (T) u m+1 W ( T ). ( + dξ + ( dξ L () ψ 3(α m (t)) ψ 3 (α m (t )) α m (t) u m+1,ξ (t) t t 1+/2 K 1 +K 2 +K 3. Using the fomula (3.6) ψ 3 (α m (t)) ψ 3 (α m (α m (t )) we obtain K 1 +K 2 ( φ(a m (T)) dξ ψ 3(α m (t )) α m (t) α m (t ) u m+1,ξ (t) t t 1+/2 ψ 3(α m (t)) α m (t) u m+1,ξ (t) u m+1,ξ (t ) t t 1+/2 = (α m (t) α m (t )) ( φ(a m (T)) dξ 1 ( φ(a m (T)) u m W ( T ) dξ ψ 3,αm (t )(α m (t )+s(α m (t) α m (t ))ds Ì t u t m,ξ dτ u m+1,ξ (t) t t 1+/2 t t t /2 2 Integating with espect to t we get (/2 > 1) t u m,ξ dτ u m+1,ξ (t) t t /2 2 u m+1,ξ K4. K 4 φ(a m (T))T /2 1 u m W ( T ) u m+1 W ( T ).

9 Finally Fee bounday poblem 327 K 3 φ(a m (T))a m (T) u m+1 W ( T ). Summaizing the above consideations we obtain Similaly, we obtain I φ(a m (T))a m (T) u m+1 W ( T ) +φ(a m (T))T /2 1 u m W ( T ) u m+1 W ( T ). J φ(a m (T))a m (T) q m+1 W 1,1/2 ( T ) +φ(a m (T))T /2 1 u m W ( T ) q m+1 W 1,1/2 ( T ), whee q m+1 is an extension of q m+1 W (S T ). The fifth tem on the.h.s. of (3.5) is estimated by ψ 5 (α m ) W whee M 1 φ(a m (T)) ( M 2 dξ S ( using (3.6) we have ( φ(a m (T)) dξ S ( dξ S (S T ) ( t ψ 5 (α m (t)) 1 1/ W (S) ( + ψ 5 (α m (t)) 1/2 1/(2) W (,T) dξ S M 1 +M 2, u m,ξ dτ W 1 () φ(am (T))T 1/ a m (T), ψ 5(α m (t)) ψ 5 (α m (t )) t t 1+(1/2 1/(2)) Ì t u t m,ξ (ξ,τ)dτ t t 1+(1/2 1/(2)) u m,ξ (ξ,τ) dτ ( φ(a m (T))T 1/2+/2 u m W ( T ). ; t t /2 3/2 The seventh tem of the.h.s. of (3.5) will be consideed in the fom ((I A (u m ))u m+1 ) t L ( T ) ψ 6 (α m )α m u m+1,t L ( T ) + ψ 7 (α m )u m,ξ u m+1 L ( T ) N 1 +N 2

10 328 P.B.MuchaandW.Zaja czkowski and we have (3.7) N 1 φ(a m (T)) α m u m+1 W ( T ), ( N 2 φ(a m (T)) (u m,ξ v,ξ +v,ξ ) v + φ(a m (T))T 1/ v 2 W 2 2/ () +φ(a m (T)) u m W ( T ) a m+1 +φ(a m (T)) v 2 2/ W () Tβ ( u m W ( T ) + v W 2 2/ () ). t u m+1,t ) L(T) In the last tem of the.h.s. of (3.7) 2 we have applied the imbedding W 1,1/2 ( T ) C β (,T;L ()) with < β < 1/2 1/. This enables us to get u m,ξ v,ξ L T β ( u m W ( T ) + v W 2 2/ () ). Finally we conside the second tem of the.h.s. of (3.5). We have J xu m um U um = (ξ,t) um x um (ξ,t) x um (ξ,t) dξ = = Assuming that um x um (ξ,t) (x um (ξ,t) x um (ξ,t)) x um (ξ,t) x um (ξ,t) 3 J xu m (ξ,t)dξ um x um (ξ,t) (ξ ξ Ì Ì 1 )(1+ dsì t su m (ξ +s(ξ ξ ),τ)dτ) ξ ξ dsì t su m (ξ +s(ξ ξ ),τ)dτ 3 J xu m (ξ,t)dξ. t u mξ (,τ) L ()dτ < 1, we obtain um U um L ( T ) φ(a m (T))T 1/ dξ ξ ξ 2 φ(a m (T))T 1/. L () Fo simplicity we intoduce X k = u k W ( T ) + q k W 1, ( T ) + q m W (S T ).

11 Fee bounday poblem 329 Summing up the estimates fo all tems of the.h.s. of (3.5) we get X m+1 a m φ(a m )X m+1 +φ(a m )T /2 1 X m X m+1 Putting we have (3.8) +φ(a m )a m T 1/ +φ(a m )T 1/2+/2 X m +φ(a m )X m a m+1 (T) +(T 1/ +T β )φ(a m )+φ(a m (T))T β X m. { 1 a = min, 2 2, 1, } 2,β X m+1 T a φ(a m )X m X m+1 +T a φ(a m )X m +φ(a m )T a X m X m+1 +T a φ(a m ). By induction we pove that X k 1 (X = ). Taking T such that T 1 and T a φ(1) 1/4, inseting X m 1 in (3.8), we obtain X m X m X m , which gives X m+1 1. The poof of the lemma is complete. Lemma 3.2. Assume that S W 2 1/ exist u W ( T ) and p W 1,, v W 2 2/ ( T ) W (). Then thee (S T ) such that (S T ) u m u in W ( T ) and q m p in W 1, ( T ) W as m fo T small enough. Poof. We show that {(u m,q m )} n=1 is convegent. Fo this pupose we conside v m = u m+1 u m, m = q m+1 q m which satisfy the system (3.9) v m,t divt(v m, m ) = div um T um (u m+1,q m+1 ) div um 1 T um 1 (u m,q m ) divt(v m,q m ) + um U um um 1 U um 1 I in T, divv m = divv m div um u m+1 +div um 1 u m J in T, n T(v m, m ) = n T(v m, m ) n um T(u m+1,q m+1 ) +n um 1 (u m,q m ) p n um +p n um 1 K on S T, v m t= = on. By Theoem 1.2 we obtain an estimate on solutions of (3.9): (3.1) v m W ( T ) + m W 1, ( T ) + m W (S T ) c( I L ( T ) + J W 1, ( T ) + K + B W (S T ) L ( )), T whee B is defined by the elation J t = divb.

12 33 P.B.MuchaandW.Zaja czkowski Fist we estimate the tems of the.h.s. of (3.9) 1 in L ( T ). Let I = I 1 +I 2. We examine I 1 = div um T um (u m+1,q m+1 ) div um 1 T um 1 (u m,q m ) divt(v m,q m ) = (div um T um divt)(v m, m ) +(div um T um div um 1 T um 1 )(u m,q m ) I 11 +I 12. I 11 is estimated in the same way as the fist tem of the.h.s. of (3.5): I 11 L ( T ) φ(a m )a m ( v W ( T ) + m W 1, ( T ) ). Fo the second tem we have I 12 L ( T ) = (div um (T um T um 1 ) Next we conside +(div um div um 1 )T um 1 )(u m,q m ) L ( T ) φ(a m )T ( 1)/ v m 1 W ( T ). I 2 = um U um um 1 U um 1 = um (U um U um 1 )+( um um 1 )U um 1 I 21 +I 22. The fist tem is I 21 = A(u m ) ( xm (ξ,t) y m (ξ,t) x m (ξ,t) y m (ξ,t) 3J y m (ξ,t) x m 1(ξ,t) y m 1 (ξ,t) x m 1 (ξ,t) y m 1 (ξ,t) 3 ) dξ x m 1(ξ,t) y m 1 (ξ,t) x m 1 (ξ,t) y m 1 (ξ,t) 3J y m 1 (ξ,t) x m (ξ,t) y m (ξ,t) = A(u m ) x m (ξ,t) y m (ξ,t) 3(J y m (ξ,t) J ym 1 (ξ,t))dξ ( xm (ξ,t) y m (ξ,t) A(u m ) x m (ξ,t) y m (ξ,t) 3 ) J ym 1 (ξ,t)dξ I 211 +I 212, whee x k (ξ,t) = ξ + Ì t u k(ξ,t ) and y k (ξ,t) = ξ + Ì t u k(ξ,t ). Since we have J ym (ξ,t) J ym 1 (ξ,t) c u m u m 1 W ( T ), The same holds fo I 212. Since I 211 L ( T ) φ(a m )T 1/ v m 1 W ( T ).

13 Fee bounday poblem 331 x m (ξ,t) y m (ξ,t) x m (ξ,t) y m (ξ,t) 3 x m 1(ξ,t) y m 1 (ξ,t) x m 1 (ξ,t) y m 1 (ξ,t) 3 dξ we obtain and c u m u m 1 W ( T ), I 212 L ( T ) φ(a m )T 1/ v m 1 W ( T ). We estimate J in W 1, ( T ) by the same quantity. Let K = K 1 +K 2, whee K 1 = n T(v m, m ) n um T(u m+1,q m+1 )+n um 1 (u m,q m ) = (n T(v m, m ) n um T(v m, m )) +(n um T um n um 1 T um 1 )(u m,q m ) K 11 +K 12 K 2 = p (n um n um 1 ). The tem K 11 is estimated just as the fouth tem of the.h.s. of (3.5): K 11 W (S T ) φ(a m(t))a m (T) v m W ( T ) The second tem is +φ(a m (T))T /2 1 u m W ( T ) v m W ( T ). K 12 = (n um (T um T um 1 )+(n um n um 1)T um 1 )(u m,q m ) K 121 +K 122. Fo K 121 we have K 121 W (S T ) φ(a m(t))t ( 1)/ v m 1 W ( T ). Since n um n um 1 W (S T ) φ(a m(t))t ( 1)/ v m 1 W ( T ), we conclude that K 122 W (S T ) + K 2 W (S T ) φ(a m (T))T ( 1)/ v m 1 W ( T ). Finally we have to examine B which is defined by J t = divb. We have J = divv m div um u m 1 +div um 1 u m = (divv m div um v m )+(div um 1 u m div um u m ) J 1 +J 2. To examine J 1 we poceed as in the case of the seventh tem of the.h.s. of (3.5). By the same agument as in Lemma 3.1 we have J 2 = (div um 1 div um )u m = div ((A m 1 A m )u m). Hence we put B 2 = ((A m 1 A m)u m ) t and we get, just as fo N 2 in the poof of Lemma 3.1, the following estimate:

14 332 P.B.MuchaandW.Zaja czkowski B 2 L ( T ) φ(a m (T))(T ( 1)/ +T β ) v m 1 W ( T ), whee < β < 1/2 1/. Define Y m = v m W ( T ) + m W 1, ( T ) + m W (S T ). Summing up the estimates fo all tems of the.h.s. of (3.1) we obtain Y m φ(a m (T))(T ( 1)/ +T /2 1 +T β )Y m +φ(a m (T))(T ( 1)/ +T β )Y m 1. Taking T so small that φ(a m (T))(T ( 1)/ +T /2 1 +T β ) 1/2 we get Y m φ(a m (T))(T ( 1)/ +T β )Y m 1. Thus if φ(a m (T))(T ( 1)/ +T β ) < 1 we have a contaction, hence Y m as m. This yields the existence of u W ( T ) and p W 1, ( T ) (S T ) such that W u m u in W ( T ), q m p in W 1, ( T ) W (S T ). The poof of the lemma is complete. By Lemma 3.2 we see that system (3.1) has a unique solution (u,q) in W ( T ) W 1, ( T ) W (S T ). By Lemma 3.1 we get the estimate (3.11) u W ( T ) + q W 1, ( T ) + q W (S T ) c( v W 2 2/ (), S W 2 2/ Since fo u W ( T ) the tansfomation (1.7) is invetible, fom (3.11) we obtain estimate (1.8). Theoem 1.1 is poved. ). Refeences [1] O.V.Besov,V.P.Il inands.m.nikol skiĭ,integalrepesentationsoffunctions and Imbedding Theoems, Nauka, Moscow, 1975(in Russian). [2] O.A.Ladyzhenskaya,V.A.SolonnikovandN.N.Ual tseva,lineaand Quasilinea Equations of Paabolic Type, Ame. Math. Soc., Povidence, RI, [3] P.B.MuchaandW.M.Zaj aczkowski,ontheexistencefothecauchy Neumann poblemfothestokessysteminthel p-famewok,studiamath.,toappea. [4] V.A.Solonnikov,Onnonstationaymotionofanisolatedvolumeofaviscous incompessible fluid, Izv. Akad. Nauk SSSR 51(1987), (in Russian).

15 Fee bounday poblem 333 [5] V.A.Solonnikov,Solvabilityonafinitetimeintevalofthepoblemofevolution of a viscous incompessible fluid bounded by a fee suface, Algeba Anal. 3(1991), (in Russian). Piot Bogus law Mucha Institute of Applied Mathematics and Mechanics Wasaw Univesity Banacha Waszawa, Poland mucha@hyda.mimuw.edu.pl Wojciech Zaja czkowski Institute of Mathematics Polish Academy of Sciences Śniadeckich 8-95 Waszawa, Poland and Institute of Mathematics and Opeations Reseach Militay Univesity of Technology Kaliskiego Waszawa, Poland wz@impan.gov.pl Received on ; evised vesion on

As is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3.

As is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3. Appendix A Vecto Algeba As is natual, ou Aeospace Stuctues will be descibed in a Euclidean thee-dimensional space R 3. A.1 Vectos A vecto is used to epesent quantities that have both magnitude and diection.

More information

Green s Identities and Green s Functions

Green s Identities and Green s Functions LECTURE 7 Geen s Identities and Geen s Functions Let us ecall The ivegence Theoem in n-dimensions Theoem 7 Let F : R n R n be a vecto field ove R n that is of class C on some closed, connected, simply

More information

PROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.

PROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr. POBLM S # SOLUIONS by obet A. DiStasio J. Q. he Bon-Oppenheime appoximation is the standad way of appoximating the gound state of a molecula system. Wite down the conditions that detemine the tonic and

More information

Stress, Cauchy s equation and the Navier-Stokes equations

Stress, Cauchy s equation and the Navier-Stokes equations Chapte 3 Stess, Cauchy s equation and the Navie-Stokes equations 3. The concept of taction/stess Conside the volume of fluid shown in the left half of Fig. 3.. The volume of fluid is subjected to distibuted

More information

J. N. R E DDY ENERGY PRINCIPLES AND VARIATIONAL METHODS APPLIED MECHANICS

J. N. R E DDY ENERGY PRINCIPLES AND VARIATIONAL METHODS APPLIED MECHANICS J. N. E DDY ENEGY PINCIPLES AND VAIATIONAL METHODS IN APPLIED MECHANICS T H I D E DI T IO N JN eddy - 1 MEEN 618: ENEGY AND VAIATIONAL METHODS A EVIEW OF VECTOS AND TENSOS ead: Chapte 2 CONTENTS Physical

More information

KOEBE DOMAINS FOR THE CLASSES OF FUNCTIONS WITH RANGES INCLUDED IN GIVEN SETS

KOEBE DOMAINS FOR THE CLASSES OF FUNCTIONS WITH RANGES INCLUDED IN GIVEN SETS Jounal of Applied Analysis Vol. 14, No. 1 2008), pp. 43 52 KOEBE DOMAINS FOR THE CLASSES OF FUNCTIONS WITH RANGES INCLUDED IN GIVEN SETS L. KOCZAN and P. ZAPRAWA Received Mach 12, 2007 and, in evised fom,

More information

Do not turn over until you are told to do so by the Invigilator.

Do not turn over until you are told to do so by the Invigilator. UNIVERSITY OF EAST ANGLIA School of Mathematics Main Seies UG Examination 2015 16 FLUID DYNAMICS WITH ADVANCED TOPICS MTH-MD59 Time allowed: 3 Hous Attempt QUESTIONS 1 and 2, and THREE othe questions.

More information

Dymore User s Manual Two- and three dimensional dynamic inflow models

Dymore User s Manual Two- and three dimensional dynamic inflow models Dymoe Use s Manual Two- and thee dimensional dynamic inflow models Contents 1 Two-dimensional finite-state genealized dynamic wake theoy 1 Thee-dimensional finite-state genealized dynamic wake theoy 1

More information

ON LACUNARY INVARIANT SEQUENCE SPACES DEFINED BY A SEQUENCE OF MODULUS FUNCTIONS

ON LACUNARY INVARIANT SEQUENCE SPACES DEFINED BY A SEQUENCE OF MODULUS FUNCTIONS STUDIA UNIV BABEŞ BOLYAI, MATHEMATICA, Volume XLVIII, Numbe 4, Decembe 2003 ON LACUNARY INVARIANT SEQUENCE SPACES DEFINED BY A SEQUENCE OF MODULUS FUNCTIONS VATAN KARAKAYA AND NECIP SIMSEK Abstact The

More information

GROWTH ESTIMATES THROUGH SCALING FOR QUASILINEAR PARTIAL DIFFERENTIAL EQUATIONS

GROWTH ESTIMATES THROUGH SCALING FOR QUASILINEAR PARTIAL DIFFERENTIAL EQUATIONS Annales Academiæ Scientiaum Fennicæ Mathematica Volumen 32, 2007, 595 599 GROWTH ESTIMATES THROUGH SCALING FOR QUASILINEAR PARTIAL DIFFERENTIAL EQUATIONS Teo Kilpeläinen, Henik Shahgholian and Xiao Zhong

More information

Transformation of the Navier-Stokes Equations in Curvilinear Coordinate Systems with Maple

Transformation of the Navier-Stokes Equations in Curvilinear Coordinate Systems with Maple Global Jounal of Pue and Applied Mathematics. ISSN 0973-1768 Volume 12, Numbe 4 2016, pp. 3315 3325 Reseach India Publications http://www.ipublication.com/gjpam.htm Tansfomation of the Navie-Stokes Equations

More information

SOME SOLVABILITY THEOREMS FOR NONLINEAR EQUATIONS

SOME SOLVABILITY THEOREMS FOR NONLINEAR EQUATIONS Fixed Point Theoy, Volume 5, No. 1, 2004, 71-80 http://www.math.ubbcluj.o/ nodeacj/sfptcj.htm SOME SOLVABILITY THEOREMS FOR NONLINEAR EQUATIONS G. ISAC 1 AND C. AVRAMESCU 2 1 Depatment of Mathematics Royal

More information

Review: Electrostatics and Magnetostatics

Review: Electrostatics and Magnetostatics Review: Electostatics and Magnetostatics In the static egime, electomagnetic quantities do not vay as a function of time. We have two main cases: ELECTROSTATICS The electic chages do not change postion

More information

EM Boundary Value Problems

EM Boundary Value Problems EM Bounday Value Poblems 10/ 9 11/ By Ilekta chistidi & Lee, Seung-Hyun A. Geneal Desciption : Maxwell Equations & Loentz Foce We want to find the equations of motion of chaged paticles. The way to do

More information

THE CONE THEOREM JOEL A. TROPP. Abstract. We prove a fixed point theorem for functions which are positive with respect to a cone in a Banach space.

THE CONE THEOREM JOEL A. TROPP. Abstract. We prove a fixed point theorem for functions which are positive with respect to a cone in a Banach space. THE ONE THEOEM JOEL A. TOPP Abstact. We pove a fixed point theoem fo functions which ae positive with espect to a cone in a Banach space. 1. Definitions Definition 1. Let X be a eal Banach space. A subset

More information

Asymptotically Lacunary Statistical Equivalent Sequence Spaces Defined by Ideal Convergence and an Orlicz Function

Asymptotically Lacunary Statistical Equivalent Sequence Spaces Defined by Ideal Convergence and an Orlicz Function "Science Stays Tue Hee" Jounal of Mathematics and Statistical Science, 335-35 Science Signpost Publishing Asymptotically Lacunay Statistical Equivalent Sequence Spaces Defined by Ideal Convegence and an

More information

arxiv: v1 [math.na] 8 Feb 2013

arxiv: v1 [math.na] 8 Feb 2013 A mixed method fo Diichlet poblems with adial basis functions axiv:1302.2079v1 [math.na] 8 Feb 2013 Nobet Heue Thanh Tan Abstact We pesent a simple discetization by adial basis functions fo the Poisson

More information

A NOTE ON VERY WEAK SOLUTIONS FOR A CLASS OF NONLINEAR ELLIPTIC EQUATIONS

A NOTE ON VERY WEAK SOLUTIONS FOR A CLASS OF NONLINEAR ELLIPTIC EQUATIONS SARAJEVO JOURNAL OF MATHEMATICS Vol3 15 2007, 41 45 A NOTE ON VERY WEAK SOLUTIONS FOR A CLASS OF NONLINEAR ELLIPTIC EQUATIONS LI JULING AND GAO HONGYA Abstact We pove a new a pioi estimate fo vey weak

More information

SOME GENERAL NUMERICAL RADIUS INEQUALITIES FOR THE OFF-DIAGONAL PARTS OF 2 2 OPERATOR MATRICES

SOME GENERAL NUMERICAL RADIUS INEQUALITIES FOR THE OFF-DIAGONAL PARTS OF 2 2 OPERATOR MATRICES italian jounal of pue and applied mathematics n. 35 015 (433 44) 433 SOME GENERAL NUMERICAL RADIUS INEQUALITIES FOR THE OFF-DIAGONAL PARTS OF OPERATOR MATRICES Watheq Bani-Domi Depatment of Mathematics

More information

Math 124B February 02, 2012

Math 124B February 02, 2012 Math 24B Febuay 02, 202 Vikto Gigoyan 8 Laplace s equation: popeties We have aleady encounteed Laplace s equation in the context of stationay heat conduction and wave phenomena. Recall that in two spatial

More information

working pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50

working pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50 woking pages fo Paul Richads class notes; do not copy o ciculate without pemission fom PGR 2004/11/3 10:50 CHAPTER7 Solid angle, 3D integals, Gauss s Theoem, and a Delta Function We define the solid angle,

More information

TheWaveandHelmholtzEquations

TheWaveandHelmholtzEquations TheWaveandHelmholtzEquations Ramani Duaiswami The Univesity of Mayland, College Pak Febuay 3, 2006 Abstact CMSC828D notes (adapted fom mateial witten with Nail Gumeov). Wok in pogess 1 Acoustic Waves 1.1

More information

ON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION. 1. Introduction. 1 r r. r k for every set E A, E \ {0},

ON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION. 1. Introduction. 1 r r. r k for every set E A, E \ {0}, ON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION E. J. IONASCU and A. A. STANCU Abstact. We ae inteested in constucting concete independent events in puely atomic pobability

More information

Appendix A. Appendices. A.1 ɛ ijk and cross products. Vector Operations: δ ij and ɛ ijk

Appendix A. Appendices. A.1 ɛ ijk and cross products. Vector Operations: δ ij and ɛ ijk Appendix A Appendices A1 ɛ and coss poducts A11 Vecto Opeations: δ ij and ɛ These ae some notes on the use of the antisymmetic symbol ɛ fo expessing coss poducts This is an extemely poweful tool fo manipulating

More information

Journal of Inequalities in Pure and Applied Mathematics

Journal of Inequalities in Pure and Applied Mathematics Jounal of Inequalities in Pue and Applied Mathematics COEFFICIENT INEQUALITY FOR A FUNCTION WHOSE DERIVATIVE HAS A POSITIVE REAL PART S. ABRAMOVICH, M. KLARIČIĆ BAKULA AND S. BANIĆ Depatment of Mathematics

More information

A STABILITY RESULT FOR p-harmonic SYSTEMS WITH DISCONTINUOUS COEFFICIENTS. Bianca Stroffolini. 0. Introduction

A STABILITY RESULT FOR p-harmonic SYSTEMS WITH DISCONTINUOUS COEFFICIENTS. Bianca Stroffolini. 0. Introduction Electonic Jounal of Diffeential Equations, Vol. 2001(2001), No. 02, pp. 1 7. ISSN: 1072-6691. URL: http://ejde.math.swt.edu o http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) A STABILITY RESULT

More information

Numerical approximation to ζ(2n+1)

Numerical approximation to ζ(2n+1) Illinois Wesleyan Univesity Fom the SelectedWoks of Tian-Xiao He 6 Numeical appoximation to ζ(n+1) Tian-Xiao He, Illinois Wesleyan Univesity Michael J. Dancs Available at: https://woks.bepess.com/tian_xiao_he/6/

More information

2. Electrostatics. Dr. Rakhesh Singh Kshetrimayum 8/11/ Electromagnetic Field Theory by R. S. Kshetrimayum

2. Electrostatics. Dr. Rakhesh Singh Kshetrimayum 8/11/ Electromagnetic Field Theory by R. S. Kshetrimayum 2. Electostatics D. Rakhesh Singh Kshetimayum 1 2.1 Intoduction In this chapte, we will study how to find the electostatic fields fo vaious cases? fo symmetic known chage distibution fo un-symmetic known

More information

On absence of solutions of a semi-linear elliptic equation with biharmonic operator in the exterior of a ball

On absence of solutions of a semi-linear elliptic equation with biharmonic operator in the exterior of a ball Tansactions of NAS of Azebaijan, Issue Mathematics, 36, 63-69 016. Seies of Physical-Technical and Mathematical Sciences. On absence of solutions of a semi-linea elliptic euation with bihamonic opeato

More information

A dual-reciprocity boundary element method for axisymmetric thermoelastodynamic deformations in functionally graded solids

A dual-reciprocity boundary element method for axisymmetric thermoelastodynamic deformations in functionally graded solids APCOM & ISCM 11-14 th Decembe, 013, Singapoe A dual-ecipocity bounday element method fo axisymmetic themoelastodynamic defomations in functionally gaded solids *W. T. Ang and B. I. Yun Division of Engineeing

More information

New problems in universal algebraic geometry illustrated by boolean equations

New problems in universal algebraic geometry illustrated by boolean equations New poblems in univesal algebaic geomety illustated by boolean equations axiv:1611.00152v2 [math.ra] 25 Nov 2016 Atem N. Shevlyakov Novembe 28, 2016 Abstact We discuss new poblems in univesal algebaic

More information

TUTORIAL 9. Static magnetic field

TUTORIAL 9. Static magnetic field TUTOIAL 9 Static magnetic field Vecto magnetic potential Null Identity % & %$ A # Fist postulation # " B such that: Vecto magnetic potential Vecto Poisson s equation The solution is: " Substitute it into

More information

Hydroelastic Analysis of a 1900 TEU Container Ship Using Finite Element and Boundary Element Methods

Hydroelastic Analysis of a 1900 TEU Container Ship Using Finite Element and Boundary Element Methods TEAM 2007, Sept. 10-13, 2007,Yokohama, Japan Hydoelastic Analysis of a 1900 TEU Containe Ship Using Finite Element and Bounday Element Methods Ahmet Egin 1)*, Levent Kaydıhan 2) and Bahadı Uğulu 3) 1)

More information

Mean Curvature and Shape Operator of Slant Immersions in a Sasakian Space Form

Mean Curvature and Shape Operator of Slant Immersions in a Sasakian Space Form Mean Cuvatue and Shape Opeato of Slant Immesions in a Sasakian Space Fom Muck Main Tipathi, Jean-Sic Kim and Son-Be Kim Abstact Fo submanifolds, in a Sasakian space fom, which ae tangential to the stuctue

More information

What Form of Gravitation Ensures Weakened Kepler s Third Law?

What Form of Gravitation Ensures Weakened Kepler s Third Law? Bulletin of Aichi Univ. of Education, 6(Natual Sciences, pp. - 6, Mach, 03 What Fom of Gavitation Ensues Weakened Keple s Thid Law? Kenzi ODANI Depatment of Mathematics Education, Aichi Univesity of Education,

More information

DonnishJournals

DonnishJournals DonnishJounals 041-1189 Donnish Jounal of Educational Reseach and Reviews. Vol 1(1) pp. 01-017 Novembe, 014. http:///dje Copyight 014 Donnish Jounals Oiginal Reseach Pape Vecto Analysis Using MAXIMA Savaş

More information

Vectors, Vector Calculus, and Coordinate Systems

Vectors, Vector Calculus, and Coordinate Systems ! Revised Apil 11, 2017 1:48 PM! 1 Vectos, Vecto Calculus, and Coodinate Systems David Randall Physical laws and coodinate systems Fo the pesent discussion, we define a coodinate system as a tool fo descibing

More information

Semicanonical basis generators of the cluster algebra of type A (1)

Semicanonical basis generators of the cluster algebra of type A (1) Semicanonical basis geneatos of the cluste algeba of type A (1 1 Andei Zelevinsky Depatment of Mathematics Notheasten Univesity, Boston, USA andei@neu.edu Submitted: Jul 7, 006; Accepted: Dec 3, 006; Published:

More information

Math Notes on Kepler s first law 1. r(t) kp(t)

Math Notes on Kepler s first law 1. r(t) kp(t) Math 7 - Notes on Keple s fist law Planetay motion and Keple s Laws We conside the motion of a single planet about the sun; fo simplicity, we assign coodinates in R 3 so that the position of the sun is

More information

A PROOF OF THE INF-SUP CONDITION FOR THE STOKES EQUATIONS ON LIPSCHITZ DOMAINS

A PROOF OF THE INF-SUP CONDITION FOR THE STOKES EQUATIONS ON LIPSCHITZ DOMAINS A PROOF OF THE INF-SUP CONDITION FOR THE STOKES EQUATIONS ON LIPSCHITZ DOMAINS JAMES H. BRAMBLE Abstact. The pupose of this pape is to pesent a athe simple poof of an inequality of Nečas [9] which is equivalent

More information

(read nabla or del) is defined by, k. (9.7.1*)

(read nabla or del) is defined by, k. (9.7.1*) 9.7 Gadient of a scala field. Diectional deivative Some of the vecto fields in applications can be obtained fom scala fields. This is vey advantageous because scala fields can be handled moe easily. The

More information

Supplementary material for the paper Platonic Scattering Cancellation for Bending Waves on a Thin Plate. Abstract

Supplementary material for the paper Platonic Scattering Cancellation for Bending Waves on a Thin Plate. Abstract Supplementay mateial fo the pape Platonic Scatteing Cancellation fo Bending Waves on a Thin Plate M. Fahat, 1 P.-Y. Chen, 2 H. Bağcı, 1 S. Enoch, 3 S. Guenneau, 3 and A. Alù 2 1 Division of Compute, Electical,

More information

A CHARACTERIZATION ON BREAKDOWN OF SMOOTH SPHERICALLY SYMMETRIC SOLUTIONS OF THE ISENTROPIC SYSTEM OF COMPRESSIBLE NAVIER STOKES EQUATIONS

A CHARACTERIZATION ON BREAKDOWN OF SMOOTH SPHERICALLY SYMMETRIC SOLUTIONS OF THE ISENTROPIC SYSTEM OF COMPRESSIBLE NAVIER STOKES EQUATIONS Huang, X. and Matsumua, A. Osaka J. Math. 52 (215), 271 283 A CHARACTRIZATION ON BRAKDOWN OF SMOOTH SPHRICALLY SYMMTRIC SOLUTIONS OF TH ISNTROPIC SYSTM OF COMPRSSIBL NAVIR STOKS QUATIONS XIANGDI HUANG

More information

Perturbation to Symmetries and Adiabatic Invariants of Nonholonomic Dynamical System of Relative Motion

Perturbation to Symmetries and Adiabatic Invariants of Nonholonomic Dynamical System of Relative Motion Commun. Theo. Phys. Beijing, China) 43 25) pp. 577 581 c Intenational Academic Publishes Vol. 43, No. 4, Apil 15, 25 Petubation to Symmeties and Adiabatic Invaiants of Nonholonomic Dynamical System of

More information

arxiv: v1 [math.ca] 31 Aug 2009

arxiv: v1 [math.ca] 31 Aug 2009 axiv:98.4578v [math.ca] 3 Aug 9 On L-convegence of tigonometic seies Bogdan Szal Univesity of Zielona Góa, Faculty of Mathematics, Compute Science and Econometics, 65-56 Zielona Góa, ul. Szafana 4a, Poland

More information

Absorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere

Absorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere Applied Mathematics, 06, 7, 709-70 Published Online Apil 06 in SciRes. http://www.scip.og/jounal/am http://dx.doi.og/0.46/am.06.77065 Absoption Rate into a Small Sphee fo a Diffusing Paticle Confined in

More information

Geometry of the homogeneous and isotropic spaces

Geometry of the homogeneous and isotropic spaces Geomety of the homogeneous and isotopic spaces H. Sonoda Septembe 2000; last evised Octobe 2009 Abstact We summaize the aspects of the geomety of the homogeneous and isotopic spaces which ae most elevant

More information

COORDINATE TRANSFORMATIONS - THE JACOBIAN DETERMINANT

COORDINATE TRANSFORMATIONS - THE JACOBIAN DETERMINANT COORDINATE TRANSFORMATIONS - THE JACOBIAN DETERMINANT Link to: phsicspages home page. To leave a comment o epot an eo, please use the auilia blog. Refeence: d Inveno, Ra, Intoducing Einstein s Relativit

More information

Lacunary I-Convergent Sequences

Lacunary I-Convergent Sequences KYUNGPOOK Math. J. 52(2012), 473-482 http://dx.doi.og/10.5666/kmj.2012.52.4.473 Lacunay I-Convegent Sequences Binod Chanda Tipathy Mathematical Sciences Division, Institute of Advanced Study in Science

More information

On a quantity that is analogous to potential and a theorem that relates to it

On a quantity that is analogous to potential and a theorem that relates to it Su une quantité analogue au potential et su un théoème y elatif C R Acad Sci 7 (87) 34-39 On a quantity that is analogous to potential and a theoem that elates to it By R CLAUSIUS Tanslated by D H Delphenich

More information

Chapter 5 Force and Motion

Chapter 5 Force and Motion Chapte 5 Foce and Motion In Chaptes 2 and 4 we have studied kinematics, i.e., we descibed the motion of objects using paametes such as the position vecto, velocity, and acceleation without any insights

More information

Solution of a Spherically Symmetric Static Problem of General Relativity for an Elastic Solid Sphere

Solution of a Spherically Symmetric Static Problem of General Relativity for an Elastic Solid Sphere Applied Physics eseach; Vol. 9, No. 6; 7 ISSN 96-969 E-ISSN 96-9647 Published by Canadian Cente of Science and Education Solution of a Spheically Symmetic Static Poblem of Geneal elativity fo an Elastic

More information

Chapter 5 Force and Motion

Chapter 5 Force and Motion Chapte 5 Foce and Motion In chaptes 2 and 4 we have studied kinematics i.e. descibed the motion of objects using paametes such as the position vecto, velocity and acceleation without any insights as to

More information

Final Review of AerE 243 Class

Final Review of AerE 243 Class Final Review of AeE 4 Class Content of Aeodynamics I I Chapte : Review of Multivaiable Calculus Chapte : Review of Vectos Chapte : Review of Fluid Mechanics Chapte 4: Consevation Equations Chapte 5: Simplifications

More information

Numerical Integration

Numerical Integration MCEN 473/573 Chapte 0 Numeical Integation Fall, 2006 Textbook, 0.4 and 0.5 Isopaametic Fomula Numeical Integation [] e [ ] T k = h B [ D][ B] e B Jdsdt In pactice, the element stiffness is calculated numeically.

More information

Mathematical Model of Magnetometric Resistivity. Sounding for a Conductive Host. with a Bulge Overburden

Mathematical Model of Magnetometric Resistivity. Sounding for a Conductive Host. with a Bulge Overburden Applied Mathematical Sciences, Vol. 7, 13, no. 7, 335-348 Mathematical Model of Magnetometic Resistivity Sounding fo a Conductive Host with a Bulge Ovebuden Teeasak Chaladgan Depatment of Mathematics Faculty

More information

SOLVING THE VISCOUS COMPOSITE CYLINDER PROBLEM BY SOKOLOV S METHOD

SOLVING THE VISCOUS COMPOSITE CYLINDER PROBLEM BY SOKOLOV S METHOD 1 SOLVING THE VISCOUS COMPOSITE CYLINDER PROBLEM BY SOKOLOV S METHOD NGO By CO. H. TRAN, PHONG. T. Faculty of Mathematics & Infomatics, Univesity of Natual Sciences VNU-HCM Abstact : The pape pesents some

More information

SPECTRAL SEQUENCES. im(er

SPECTRAL SEQUENCES. im(er SPECTRAL SEQUENCES MATTHEW GREENBERG. Intoduction Definition. Let a. An a-th stage spectal (cohomological) sequence consists of the following data: bigaded objects E = p,q Z Ep,q, a diffeentials d : E

More information

Doubling property for the Laplacian and its applications (Course Chengdu 2007)

Doubling property for the Laplacian and its applications (Course Chengdu 2007) Doubling popety fo the Laplacian and its applications Couse Chengdu 007) K.-D. PHUNG The oiginal appoach of N. Gaofalo and F.H. Lin Fo simplicity, we epoduce the poof of N. Gaofalo and F.H. Lin in the

More information

Computational Methods of Solid Mechanics. Project report

Computational Methods of Solid Mechanics. Project report Computational Methods of Solid Mechanics Poject epot Due on Dec. 6, 25 Pof. Allan F. Bowe Weilin Deng Simulation of adhesive contact with molecula potential Poject desciption In the poject, we will investigate

More information

ON EXISTENCE, UNIQUENESS AND L r -EXPONENTIAL STABILITY FOR STATIONARY SOLUTIONS TO THE MHD EQUATIONS IN THREE-DIMENSIONAL DOMAINS

ON EXISTENCE, UNIQUENESS AND L r -EXPONENTIAL STABILITY FOR STATIONARY SOLUTIONS TO THE MHD EQUATIONS IN THREE-DIMENSIONAL DOMAINS ANZIAM J. 46004, 95 09 ON EXISTENCE, UNIQUENESS AND L -EXPONENTIAL STABILITY FOR STATIONARY SOLUTIONS TO THE MHD EQUATIONS IN THREE-DIMENSIONAL DOMAINS CHUNSHAN ZHAO and KAITAI LI Received 4 August, 000;

More information

12th WSEAS Int. Conf. on APPLIED MATHEMATICS, Cairo, Egypt, December 29-31,

12th WSEAS Int. Conf. on APPLIED MATHEMATICS, Cairo, Egypt, December 29-31, th WSEAS Int. Conf. on APPLIED MATHEMATICS, Caio, Egypt, Decembe 9-3, 7 5 Magnetostatic Field calculations associated with thick Solenoids in the Pesence of Ion using a Powe Seies expansion and the Complete

More information

Vectors, Vector Calculus, and Coordinate Systems

Vectors, Vector Calculus, and Coordinate Systems Apil 5, 997 A Quick Intoduction to Vectos, Vecto Calculus, and Coodinate Systems David A. Randall Depatment of Atmospheic Science Coloado State Univesity Fot Collins, Coloado 80523. Scalas and vectos Any

More information

THE LAPLACE EQUATION. The Laplace (or potential) equation is the equation. u = 0. = 2 x 2. x y 2 in R 2

THE LAPLACE EQUATION. The Laplace (or potential) equation is the equation. u = 0. = 2 x 2. x y 2 in R 2 THE LAPLACE EQUATION The Laplace (o potential) equation is the equation whee is the Laplace opeato = 2 x 2 u = 0. in R = 2 x 2 + 2 y 2 in R 2 = 2 x 2 + 2 y 2 + 2 z 2 in R 3 The solutions u of the Laplace

More information

The inviscid limit of incompressible fluid flow in an annulus

The inviscid limit of incompressible fluid flow in an annulus The inviscid limit of incompessible fluid flow in an annulus Saa Fietze, Robet Geity, and Tiago Picon ugust 8, 2006 bstact Incompessible, ciculaly symmetic fluid flow in a two-dimensional annulus = {x

More information

Method for Approximating Irrational Numbers

Method for Approximating Irrational Numbers Method fo Appoximating Iational Numbes Eic Reichwein Depatment of Physics Univesity of Califonia, Santa Cuz June 6, 0 Abstact I will put foth an algoithm fo poducing inceasingly accuate ational appoximations

More information

FREE TRANSVERSE VIBRATIONS OF NON-UNIFORM BEAMS

FREE TRANSVERSE VIBRATIONS OF NON-UNIFORM BEAMS Please cite this aticle as: Izabela Zamosa Fee tansvese vibations of non-unifom beams Scientific Reseach of the Institute of Mathematics and Compute Science Volume 9 Issue pages 3-9. The website: http://www.amcm.pcz.pl/

More information

Fluid flow in curved geometries: Mathematical Modeling and Applications

Fluid flow in curved geometries: Mathematical Modeling and Applications Fluid flow in cuved geometies: Mathematical Modeling and Applications D. Muhammad Sajid Theoetical Plasma Physics Division PINSTECH, P.O. Niloe, PAEC, Islamabad Mach 01-06, 010 Islamabad, Paistan Pesentation

More information

AP-C WEP. h. Students should be able to recognize and solve problems that call for application both of conservation of energy and Newton s Laws.

AP-C WEP. h. Students should be able to recognize and solve problems that call for application both of conservation of energy and Newton s Laws. AP-C WEP 1. Wok a. Calculate the wok done by a specified constant foce on an object that undegoes a specified displacement. b. Relate the wok done by a foce to the aea unde a gaph of foce as a function

More information

A New Approach to General Relativity

A New Approach to General Relativity Apeion, Vol. 14, No. 3, July 7 7 A New Appoach to Geneal Relativity Ali Rıza Şahin Gaziosmanpaşa, Istanbul Tukey E-mail: aizasahin@gmail.com Hee we pesent a new point of view fo geneal elativity and/o

More information

arxiv: v1 [math.co] 4 May 2017

arxiv: v1 [math.co] 4 May 2017 On The Numbe Of Unlabeled Bipatite Gaphs Abdullah Atmaca and A Yavuz Ouç axiv:7050800v [mathco] 4 May 207 Abstact This pape solves a poblem that was stated by M A Haison in 973 [] This poblem, that has

More information

MATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE. We consider second order constant coefficient scalar linear PDEs on R n. These have the form

MATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE. We consider second order constant coefficient scalar linear PDEs on R n. These have the form MATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE ANDRAS VASY We conside second ode constant coefficient scala linea PDEs on R n. These have the fom Lu = f L = a ij xi xj + b i xi + c i whee a ij b i and

More information

Vanishing lines in generalized Adams spectral sequences are generic

Vanishing lines in generalized Adams spectral sequences are generic ISSN 364-0380 (on line) 465-3060 (pinted) 55 Geomety & Topology Volume 3 (999) 55 65 Published: 2 July 999 G G G G T T T G T T T G T G T GG TT G G G G GG T T T TT Vanishing lines in genealized Adams spectal

More information

AE301 Aerodynamics I UNIT B: Theory of Aerodynamics

AE301 Aerodynamics I UNIT B: Theory of Aerodynamics AE301 Aeodynamics I UNIT B: Theoy of Aeodynamics ROAD MAP... B-1: Mathematics fo Aeodynamics B-2: Flow Field Repesentations B-3: Potential Flow Analysis B-4: Applications of Potential Flow Analysis AE301

More information

But for simplicity, we ll define significant as the time it takes a star to lose all memory of its original trajectory, i.e.,

But for simplicity, we ll define significant as the time it takes a star to lose all memory of its original trajectory, i.e., Stella elaxation Time [Chandasekha 1960, Pinciples of Stella Dynamics, Chap II] [Ostike & Davidson 1968, Ap.J., 151, 679] Do stas eve collide? Ae inteactions between stas (as opposed to the geneal system

More information

A Multivariate Normal Law for Turing s Formulae

A Multivariate Normal Law for Turing s Formulae A Multivaiate Nomal Law fo Tuing s Fomulae Zhiyi Zhang Depatment of Mathematics and Statistics Univesity of Noth Caolina at Chalotte Chalotte, NC 28223 Abstact This pape establishes a sufficient condition

More information

( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.

( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx. 9. LAGRANGIAN OF THE ELECTROMAGNETIC FIELD In the pevious section the Lagangian and Hamiltonian of an ensemble of point paticles was developed. This appoach is based on a qt. This discete fomulation can

More information

Research Article On Alzer and Qiu s Conjecture for Complete Elliptic Integral and Inverse Hyperbolic Tangent Function

Research Article On Alzer and Qiu s Conjecture for Complete Elliptic Integral and Inverse Hyperbolic Tangent Function Abstact and Applied Analysis Volume 011, Aticle ID 697547, 7 pages doi:10.1155/011/697547 Reseach Aticle On Alze and Qiu s Conjectue fo Complete Elliptic Integal and Invese Hypebolic Tangent Function Yu-Ming

More information

Syntactical content of nite approximations of partial algebras 1 Wiktor Bartol Inst. Matematyki, Uniw. Warszawski, Warszawa (Poland)

Syntactical content of nite approximations of partial algebras 1 Wiktor Bartol Inst. Matematyki, Uniw. Warszawski, Warszawa (Poland) Syntactical content of nite appoximations of patial algebas 1 Wikto Batol Inst. Matematyki, Uniw. Waszawski, 02-097 Waszawa (Poland) batol@mimuw.edu.pl Xavie Caicedo Dep. Matematicas, Univ. de los Andes,

More information

A Bijective Approach to the Permutational Power of a Priority Queue

A Bijective Approach to the Permutational Power of a Priority Queue A Bijective Appoach to the Pemutational Powe of a Pioity Queue Ia M. Gessel Kuang-Yeh Wang Depatment of Mathematics Bandeis Univesity Waltham, MA 02254-9110 Abstact A pioity queue tansfoms an input pemutation

More information

A THREE CRITICAL POINTS THEOREM AND ITS APPLICATIONS TO THE ORDINARY DIRICHLET PROBLEM

A THREE CRITICAL POINTS THEOREM AND ITS APPLICATIONS TO THE ORDINARY DIRICHLET PROBLEM A THREE CRITICAL POINTS THEOREM AND ITS APPLICATIONS TO THE ORDINARY DIRICHLET PROBLEM DIEGO AVERNA AND GABRIELE BONANNO Abstact. The aim of this pape is twofold. On one hand we establish a thee citical

More information

ECE 3318 Applied Electricity and Magnetism. Spring Prof. David R. Jackson ECE Dept. Notes 13

ECE 3318 Applied Electricity and Magnetism. Spring Prof. David R. Jackson ECE Dept. Notes 13 ECE 338 Applied Electicity and Magnetism ping 07 Pof. David R. Jackson ECE Dept. Notes 3 Divegence The Physical Concept Find the flux going outwad though a sphee of adius. x ρ v0 z a y ψ = D nˆ d = D ˆ

More information

arxiv: v1 [physics.pop-ph] 3 Jun 2013

arxiv: v1 [physics.pop-ph] 3 Jun 2013 A note on the electostatic enegy of two point chages axiv:1306.0401v1 [physics.pop-ph] 3 Jun 013 A C Tot Instituto de Física Univesidade Fedeal do io de Janeio Caixa Postal 68.58; CEP 1941-97 io de Janeio,

More information

2. Plane Elasticity Problems

2. Plane Elasticity Problems S0 Solid Mechanics Fall 009. Plane lasticity Poblems Main Refeence: Theoy of lasticity by S.P. Timoshenko and J.N. Goodie McGaw-Hill New Yok. Chaptes 3..1 The plane-stess poblem A thin sheet of an isotopic

More information

d 2 x 0a d d =0. Relative to an arbitrary (accelerating frame) specified by x a = x a (x 0b ), the latter becomes: d 2 x a d 2 + a dx b dx c

d 2 x 0a d d =0. Relative to an arbitrary (accelerating frame) specified by x a = x a (x 0b ), the latter becomes: d 2 x a d 2 + a dx b dx c Chapte 6 Geneal Relativity 6.1 Towads the Einstein equations Thee ae seveal ways of motivating the Einstein equations. The most natual is pehaps though consideations involving the Equivalence Pinciple.

More information

Classical Mechanics Homework set 7, due Nov 8th: Solutions

Classical Mechanics Homework set 7, due Nov 8th: Solutions Classical Mechanics Homewok set 7, due Nov 8th: Solutions 1. Do deivation 8.. It has been asked what effect does a total deivative as a function of q i, t have on the Hamiltonian. Thus, lets us begin with

More information

Semi-Lagrangian simulations of the diocotron instability

Semi-Lagrangian simulations of the diocotron instability Semi-Lagangian simulations of the diocoton instability Eic Madaule, Seve Adian Histoaga, Michel Mehenbege, Jéôme Péti To cite this vesion: Eic Madaule, Seve Adian Histoaga, Michel Mehenbege, Jéôme Péti.

More information

ON THE INVERSE SIGNED TOTAL DOMINATION NUMBER IN GRAPHS. D.A. Mojdeh and B. Samadi

ON THE INVERSE SIGNED TOTAL DOMINATION NUMBER IN GRAPHS. D.A. Mojdeh and B. Samadi Opuscula Math. 37, no. 3 (017), 447 456 http://dx.doi.og/10.7494/opmath.017.37.3.447 Opuscula Mathematica ON THE INVERSE SIGNED TOTAL DOMINATION NUMBER IN GRAPHS D.A. Mojdeh and B. Samadi Communicated

More information

1 Equations of linear elasticity

1 Equations of linear elasticity C5.2 Elasticity & Plasticity Daft date: 17 Januay 2017 1 1 1 Equations of linea elasticity 1.1 Hooke s law Robet Hooke (1678 wote... it is... evident that the ule o law of natue in evey spinging body is

More information

Suggested Solutions to Homework #4 Econ 511b (Part I), Spring 2004

Suggested Solutions to Homework #4 Econ 511b (Part I), Spring 2004 Suggested Solutions to Homewok #4 Econ 5b (Pat I), Sping 2004. Conside a neoclassical gowth model with valued leisue. The (epesentative) consume values steams of consumption and leisue accoding to P t=0

More information

Regularity Criteria for the Magneto-micropolar Fluid Equations in Terms of Direction of the Velocity

Regularity Criteria for the Magneto-micropolar Fluid Equations in Terms of Direction of the Velocity Applied Mathematical Sciences, Vol. 8, 01, no. 71, 51-59 HIKARI td, www.m-hikai.com http://dx.doi.og/10.1988/ams.01.05 Regulaity Citeia fo the Magneto-micopola Fluid Equations in Tems of Diection of the

More information

arxiv: v1 [math.co] 6 Mar 2008

arxiv: v1 [math.co] 6 Mar 2008 An uppe bound fo the numbe of pefect matchings in gaphs Shmuel Fiedland axiv:0803.0864v [math.co] 6 Ma 2008 Depatment of Mathematics, Statistics, and Compute Science, Univesity of Illinois at Chicago Chicago,

More information

RADIALLY SYMMETRIC SOLUTIONS TO THE GRAPHIC WILLMORE SURFACE EQUATION

RADIALLY SYMMETRIC SOLUTIONS TO THE GRAPHIC WILLMORE SURFACE EQUATION RADIALLY SYMMETRIC SOLUTIONS TO THE GRAPHIC WILLMORE SURFACE EQUATION JINGYI CHEN AND YUXIANG LI Abstact. We show that a smooth adially symmetic solution u to the gaphic Willmoe suface equation is eithe

More information

Internet Appendix for A Bayesian Approach to Real Options: The Case of Distinguishing Between Temporary and Permanent Shocks

Internet Appendix for A Bayesian Approach to Real Options: The Case of Distinguishing Between Temporary and Permanent Shocks Intenet Appendix fo A Bayesian Appoach to Real Options: The Case of Distinguishing Between Tempoay and Pemanent Shocks Steven R. Genadie Gaduate School of Business, Stanfod Univesity Andey Malenko Gaduate

More information

PHYS 705: Classical Mechanics. Small Oscillations

PHYS 705: Classical Mechanics. Small Oscillations PHYS 705: Classical Mechanics Small Oscillations Fomulation of the Poblem Assumptions: V q - A consevative system with depending on position only - The tansfomation equation defining does not dep on time

More information

Liquid gas interface under hydrostatic pressure

Liquid gas interface under hydrostatic pressure Advances in Fluid Mechanics IX 5 Liquid gas inteface unde hydostatic pessue A. Gajewski Bialystok Univesity of Technology, Faculty of Civil Engineeing and Envionmental Engineeing, Depatment of Heat Engineeing,

More information

Application of homotopy perturbation method to the Navier-Stokes equations in cylindrical coordinates

Application of homotopy perturbation method to the Navier-Stokes equations in cylindrical coordinates Computational Ecology and Softwae 5 5(): 9-5 Aticle Application of homotopy petubation method to the Navie-Stokes equations in cylindical coodinates H. A. Wahab Anwa Jamal Saia Bhatti Muhammad Naeem Muhammad

More information

On the Poisson Approximation to the Negative Hypergeometric Distribution

On the Poisson Approximation to the Negative Hypergeometric Distribution BULLETIN of the Malaysian Mathematical Sciences Society http://mathusmmy/bulletin Bull Malays Math Sci Soc (2) 34(2) (2011), 331 336 On the Poisson Appoximation to the Negative Hypegeometic Distibution

More information

A 1. EN2210: Continuum Mechanics. Homework 7: Fluid Mechanics Solutions

A 1. EN2210: Continuum Mechanics. Homework 7: Fluid Mechanics Solutions EN10: Continuum Mechanics Homewok 7: Fluid Mechanics Solutions School of Engineeing Bown Univesity 1. An ideal fluid with mass density ρ flows with velocity v 0 though a cylindical tube with cosssectional

More information

Question Bank. Section A. is skew-hermitian matrix. is diagonalizable. (, ) , Evaluate (, ) 12 about = 1 and = Find, if

Question Bank. Section A. is skew-hermitian matrix. is diagonalizable. (, ) , Evaluate (, ) 12 about = 1 and = Find, if Subject: Mathematics-I Question Bank Section A T T. Find the value of fo which the matix A = T T has ank one. T T i. Is the matix A = i is skew-hemitian matix. i. alculate the invese of the matix = 5 7

More information