Existence of equilibria in articulated bearings in presence of cavity

Size: px
Start display at page:

Download "Existence of equilibria in articulated bearings in presence of cavity"

Transcription

1 J. Math. Anal. Appl ) Existence of equilibria in articulate bearings in presence of cavity G. Buscaglia a, I. Ciuperca b, I. Hafii c,m.jai c, a Centro Atómico Bariloche an Instituto Balseiro, 8400 Bariloche, Argentina b CNRS UMR 5208 Université Lyon 1, MAPLY, Bât. 101, Villeurbanne Ceex, France c CNRS UMR 5208, INSA e Lyon, Bât. Léonar e Vinci, 69621, France Receive 20 July 2006 Available online 1 February 2007 Submitte by P.G.L. Leach Abstract The existence of equilibrium solutions for a lubricate system consisting of an articulate boy sliing over a flat plate is consiere. Here we consier the case when cavity can occur Elsevier Inc. All rights reserve. Keywors: Lubrication; Tilting-pa bearing; Equilibrium; Reynols variational inequality 1. Introuction This work is a continuation of a recent paper [1] treating the existence of equilibrium in a tilting-pa thrust bearings with cavitation effects isregare. Articulate sliers have two egrees of freeom: The first one is the vertical isplacement a uner the effect of a force F, applie at a position x 0, an of the pressure loa px where p = px) is the pressure of the flui between the surfaces of the slier. The secon egree of freeom is the tilt or pitch) angle θ see Fig. 1) see [1] for more explanations on the physical moel). * Corresponing author. aresses: gustavo@cab.cnea.gov.ar G. Buscaglia), ciuperca@maply.univ-lyon1.fr I. Ciuperca), ima.hafii@insa-lyon.fr I. Hafii), mohamme.jai@insa-lyon.fr M. Jai) X/$ see front matter 2007 Elsevier Inc. All rights reserve. oi: /j.jmaa

2 842 G. Buscaglia et al. / J. Math. Anal. Appl ) Fig. 1. Scheme of an articulate slier, containing efinitions of the ifferent components of the evice an some geometrical variables an non-imensional coorinates use in the text. In [1] we treate the case in which h, the non-imensional istance between the surfaces, is non-increasing with x. This guarantees positivity of the pressure, p, on all the omain. Thus p satisfies the Reynols equation [h 3 x) p ] = h, x, x 1 p = 0, x. 1.1) In the present work we consier that the system can have a situation with h non-increasing firstly an then non-ecreasing convergent ivergent case). In such situation the solution of 1.1) is not always non-negative an we must replace Eq. 1.1) by the corresponing variational inequation. In all this paper we restrict our stuy to the one-imensional case. Thus the equilibrium problem becomes: fin p,a,θ) K R R satisfying h 3 p ϕ p) p = F, xp = Fx 0, hx) = h 0 x) + a + θx, h ϕ p) ϕ p), ϕ K, 1.2) with =], 1[, K ={ϕ H0 1), ϕ 0}, x 0, F>0an h 0 : R a given function. We consier the following hypothesis on the function h 0 : h 0 C 2 ), h 0 > 0, x, h 0 is a non-ecreasing function, 1.6) h 0 0, x, h 0 0) = h 0 0) = ) 1.4) 1.5)

3 G. Buscaglia et al. / J. Math. Anal. Appl ) A typical example of such a function is h 0 x) = 1 x) 2. In [1] we prove the existence of at least a solution of 1.2) 1.5) with θ<0inwhich case the variational inequality 1.2) becomes the Reynols equation. The result prove in [1] says that for any F>0there exists a solution of 1.2) 1.5) with θ<0provie that the articulation point x 0 is situate not far from the right en sie x = 1) of. We then can hope to improve the result by enlarging the omain where we search for the unknown θ; it is clear that for θ h 0 ) the unique solution p of 1.2) is p 0 then 1.3) an 1.4) are impossible. In the present paper we prove that for any F>0an any x 0 there exists at least a solution of the problem 1.2) 1.5) with θ< h 0 ); this improves the result in [1] in the one-imensional case. Nevertheless, this stronger result is obtaine uner the stronger hypothesis 1.6) than that of [1]. In Section 2 we give some preliminary results on the behaviour of the solution of a variational inequality with respect to a parameter. In Section 3 we give the esire existence result by reucing our problem to an algebraic equation with an unknown θ. 2. Preliminaries In this section we shall give some preliminary results concerning the behaviour of the mapping: G : A ]0, + [ Px)x ]0, + [ with P = PA) is the unique solution of the Reynols variational inequality P K, H 0 + A) 3 p ϕ P) H 0 x ϕ P), ϕ K, 2.1) where H 0 : [0, + [ is a given boune non-negative function. We shall consier two principal situations: The function H 0 is assume ecreasing in. In this case the variational inequality becomes an equality Reynols equation 1.1) with h replace by H 0 + A). The function H 0 is assume ecreasing in [,] an increasing in [,1] for some number. In this case we must consier the general Reynols variational inequality. We give now the first general result on the behaviour of GA) when G goes to infinity. Proposition 2.1. Suppose that H 0 is boune an non-negative. Then lim GA) = 0. A + More precisely, we have 2 GA) 2 H 0 L 2 ) A 3.

4 844 G. Buscaglia et al. / J. Math. Anal. Appl ) Proof. Taking ϕ = 0 in 2.1) an since H 0 is non-negative we obtain ) P 2 A 3 H 0 x L 2) P x L 2 ) which gives easily the result by Poincaré inequality. In this section we are mainly intereste on the behaviour of GA) when A tens to 0 an on the monotony of G. These kin of problems have been consiere in [1,2] in the context of Reynols equation only. We recall these results in the next subsection an then we exten them to the case of Reynols variational inequality Equation case In this subsection the following hypothesis on H 0 will be consiere: H 0 W 1, ), min H 0 x) = 0, x H 0 x) 0 a.e.x, mes{ x : H 0 x) < 0} > 0 an we remark that in this case the problem 2.1) becomes [ H 0 + A) 3 P ] = H 0, x =], 1[, x P) = P1) = 0. We now recall a result prove in [1,2]. This result claims that GA) goes to infinity when A goes to 0 if H 0 x) is equivalent to 1 x)α with α>0 in a neighbourhoo of x = 1. Proposition 2.2. Assume, in aition to hypothesis 2.2), that there exist 0 <M 0 <M 1, α 0 an δ 0 ]0, 1[ such that 2.2) 2.3) M 0 1 x) α H 0 x) M 11 x) α a.e. x ]1 2δ 0, 1[. 2.4) Then there exists c = cm 0,M 1,δ 0,α)such that the solution P of 2.3) satisfies an 1 1 B 1 1 B ) ) B Px)x c log 1 + M 1 2, B ]0,δ 0 ], A>0, if α = 0, A Px)x ca 2 α 2, B ]0,δ 0 ], A<B α, if α>0. Corollary 2.3. Uner hypotheses of Proposition 2.2, we have lim GA) =+. A 0 We now give the following technical result.

5 G. Buscaglia et al. / J. Math. Anal. Appl ) Lemma 2.4. Uner the hypothesis 2.2) an in aition Θ + 2M 0 1 x) H 0 x) Θ + 2M 11 x), x [1 2δ 0, 1], 2.5) with 0 <Θ 1, 0 <M 0 <M 1 an δ 0 ]0, 1[, then there exists C = CM 0,M 1,δ 0 )>0 such that Px C A + Θ 2 for all A>0, such that A + Θ 2 <δ0 2. Proof. We follow the lines of the proof of Proposition 2.2. Step 1: We take φ C 2 ) with φ 1on[1 δ 0, 1] an φ1 2δ 0 ) = 0. We first show that there exists C 1 > 0 epening only on M 0,M 1 an δ 0 such that Px) C 1 q 1 x)φx), x [1 2δ 0, 1], with q 1 x) = Θ1 x) 2 + M 0 1 x) 3 [Θ1 x) + M 1 1 x) 2 + A] 3. Proof of Step 1. We apply the maximum principle on [1 2δ 0, 1]. Since q 1 1) = 0 an φ1 2δ 0 ) = 0, it suffices to prove that [ C 1 A + H 0 ) 3 q ] 1φ) H 0 x x x), x [1 2δ 0, 1]. Accoring to 2.5), it is enough to show the existence of C 2 > 0 epenent on M 0,M 1 an δ 0 such that [ A + H 0 ) 3 q ] 1φ) C 2 Θ + 2M0 1 x) ). x x By eveloping, it is enough to show with I k C 2 Θ + 2M0 1 x) ), k= 1, 2,...,5, 2.6) 5 I 1 = 3 A + H 0 ) 2 H 0 q 1 φ, I 2 = 3 A + H 0 ) 2 H 0 φ q 1, I 3 = A + H0 ) 3 q 1, I 4 = 2 A + H0 ) 3 q 1 φ, I 5 = A + H0 ) 3 q 1 φ. Some elementary calculations yiels q 1 x) Θ1 x) + M 0 1 x) 2 C3 [Θ1 x) + M 1 1 x) 2 + A] 3, x [1 δ 0, 1[, 2.7) q 1 x) Θ + 1 x C4 [Θ1 x) + M 1 1 x) 2 + A] 3, x [1 δ 0, 1[, 2.8) with C 3,C 4 > 0 epening only on M 0 an M 1.

6 846 G. Buscaglia et al. / J. Math. Anal. Appl ) Now, integrating 2.5) between x an 1 with x ], 1[ an aing A we obtain A + Θ1 x) + M 0 1 x) 2 A + H 0 x) A + Θ1 x) + M 1 1 x) 2, x [1 2δ 0, 1[. Using now the expression of q 1, the inequalities 2.7) 2.9) an the hypothesis 2.5), we easily euce 2.6) which ens the proof of Step 1. Step 2: We prove the esire inequality. Using the elementary inequality Θ1 x) 1 x) 2 + Θ 2, we euce Θ1 x) + M 1 1 x) 2 + A M 1 + 1)1 x) 2 + A + Θ 2. From the non-negativity of P an Step 1 we euce 1 1 P P C 1 q 1 x) x, 1 δ 0 1 δ 0 then P C 1 M δ 0 1 x) 3 x M 1 + 1)1 x) 2 + A + Θ 2 ) 3. Let x ]1 δ 0, 1[ such that 1 x = A + Θ 2. Then P C 1 M 0 M 1 + 2) 3 A + Θ 2 ) 3 1 x 1 x) 3 x C 1 M 0 = A + Θ 2 ) 2, 4M 1 + 2) 3 A + Θ 2 ) 3 which ens the proof. We finally give the regularity an monotonicity result which is prove in [1]. Proposition 2.5. The function G is of class C an G A) < 0, A> Inequation case In all this subsection the following hypothesis on the given function H 0 is consiere H 0 C 2 ) an H 0 x) > 0, x, an there exists ], 1[ such that H 0 ) = H ) ) = 0, H 0 x) < 0 for x< an H 0 x) > 0 for x>. The hypothesis 2.10) implies H 0 x) > 0 for any x {}. We recall that we consier P epening on A) the solution of 2.1). The analoguous of Proposition 2.2 is the following 2.9)

7 G. Buscaglia et al. / J. Math. Anal. Appl ) Proposition 2.6. Let 0 ], 1[. Then there exists C>0epening only on 0, min x H 0 x) an max x H 0 x) such that for any > 0, we have Px)x C ) A for any A with 0 <A< Proof. From results of [3], we have that P>0on],[, an P satisfies H 0 + A) 3 P ) = H 0 x, x ],[. Since P) = 0 an P) 0, we have by maximum principle P q, x [,], with q the solution of the problem H 0 + A) 3 q ) = H 0 x, q) = q) = 0. x ],[, Using now the Taylor evelopment of H 0 x) aroun we easily obtain 2.11) 1 2 min H 0 y) x) H 0 x) 1 y ],[ 2 min H 0 y) x) for any x ],[. y ],[ 2.12) By an affine a change of variables: x ], 1[ x ],[, x x = x + 1 2, 2.13) the problem 2.11) can be written in ], 1[ an we can easily see that the hypothesis 2.4) is satisfie with α = 1. Applying Proposition 2.2 we euce the result since P 0on. Corollary 2.7. lim GA) =+. A 0 We now stuy the monotony of GA). From the results of [3] an assumptions 2.10), we have two possibilities: Px)>0, x, an P satisfies H 0 + A) 3 P ) = H 0, x, x P) = P1) = ) There exists β ],1[ such that Px)>0, x ],β[, Pβ)= P β) = 0 an Px)= 0, x [β,1].

8 848 G. Buscaglia et al. / J. Math. Anal. Appl ) In this case P satisfies the problem H 0 + A) 3 P ) = H 0 x, x ],β[, P) = Pβ)= 0, P β) = 0, Px)= 0, x [β,1]. 2.15) We remark that β is here an unknown. By integrating 2.15), we obtain H 0 + A) 3 P x = H 0 C an P β) = 0givesC = H 0 β). We then euce P 1 x) = H 0 x) + A) 2 H 0β) + A H 0 x) + A) 3. Due to P) = 0 we obtain x Px)= s H 0 s) + A) 2 H 0 β) + A ) x s H 0 s) + A) ) Finally the conition Pβ)= 0, allows us to obtain a scalar equation satisfie by β. It is then natural to introuce the following application: ψ : [, 1] ]0, + [ R, efine by β ψβ,a)= s H 0 s) + A) 2 H 0 β) + A ) β s H 0 s) + A) ) Then for any A, the problem 2.15) reuces to the existence of β ],1[ such that ψβ,a)= 0. It is clear that ψ C 2 ], 1[ ]0, + [) an we have ψ β β β, A) = H 0 β) s H 0 s) + A) ) Hypotheses 2.10) imply that for any A,max β [,1] ψβ,a)= ψ,a), ψ β > 0forβ ],[ an ψ β < 0forβ ],1[. Since ψ,a) = 0 an from the above consierations, for any A>0 the existence an uniqueness of a solution of 2.15) is equivalent to the conition ψ1,a)<0. Then we have prove: Proposition 2.8. If ψ1,a)<0, then there exists β = βa) ],1[ such that P satisfies 2.15). If ψ1,a) 0, then P satisfies 2.14).

9 G. Buscaglia et al. / J. Math. Anal. Appl ) The next lemma is useful in the stuy of the monotony of G. Lemma 2.9. Assume that H 0 is non-ecreasing. Let ν ],1] such that H 0)>H 0 ν) an φ : [0, + [ ]0, + [ an arbitrary continuous function. Let q be the solution of the problem φh 0 ) q q) = qν) = 0. Then ν qx)x > 0. ) = H 0 x, x ],ν[, 2.19) Remark This result can be seen as a kin of maximum principle in the sense: if 1 H 0 x < 0, then 1 qx>0. Proof. By integrating 2.19) there exists c R such that φh 0 )q = H 0 c. 2.20) Due to the assumptions 2.10), three ifferent cases are possible: Case 1: q x) 0, x ],ν[ or q x) 0, x ],ν[ with mes{x: q x) = 0}=0; this is impossible because q) = qν) = 0. Case 2: There exists ν 1 ],ν[ such that q x) > 0 for x ],ν 1 [, q x) < 0 for x ]ν 1,ν[, q ν 1 ) = 0. It is clear that qx) > 0, x ],ν[ which proves the result. Case 3: There exist ν 1 ],[ an ν 2 ],ν[ such that q x) > 0 if x ],ν 1 [ ]ν 2,ν[, q x) < 0 if x ]ν 1,ν 2 [, q ν 1 ) = q ν 2 ) = 0. Let ν ],ν[ such that H 0 ν ) = H 0 ν). It is clear that ν <ν 1. We have: ν ν 0 x)qx) x = H 0 q x = H ν ν H 0 q x From 2.20) an using also H 0 ν ) = H 0 ν) we euce ν ν H 0 q x = 0. ν H 0 q x. 2.21)

10 850 G. Buscaglia et al. / J. Math. Anal. Appl ) We then easily obtain from 2.21) ν H 0 qx)x > 0. Let ν 3 ]ν 1,ν 2 [ such that qν 3 ) = 0, qx) > 0ifx ],ν 3 [ an qx) < 0ifx ]ν 3, 1[. Since H 0 is non-ecreasing, we have ν qx)x > 1 H 0 ν 3) ν 3 H 0 x)qx) x + 1 H Together with 2.22) this proves the claime result. 0 ν 3) We have a first result concerning the monotony of G. ν ν 3 H 0 x)qx) x. 2.22) Proposition Suppose that H 0 is non-ecreasing an A>0 such that ψ1,a) 0. Then G is erivable in A an we have G A) < 0. A Proof. Due to Proposition 2.8, we istinguish two cases following the sign of ψ1,a). Case 1: ψ1,a)>0. Uner this hypothesis P satisfies 2.14) equation case). The continuity of ψ implies ψ1,a )>0 for any A in a neighborhoo of A. From the Implicit Function Theorem, the solution P is erivable in A. Denoting q 1 = P A,wehaveG A) = q 1x) x where q 1 is the solution of the problem H 0 + A) 3 q ) 1 = 3 H 0 + A) 2 P ), x, x x q 1 ) = q 1 1) = ) By integrating 2.23) there exists a constant C 1 A) R such that H 0 + A) 3 q 1 x = 3H 0 + A) 2 P x + C 1A). 2.24) Multiplying 2.24) by H 0 + A an ifferentiating in x we obtain H 0 + A) 4 q ) 1 = 3 H 0 + A) 3 P ) + C 1 A) H 0 x x x. Using 2.14) we euce H 0 + A) 4 q ) 1 = C 1 A) 3 ) H 0 x x x. By linearity we can write q 1 = C 1 A) 3 ) ˆq ) with ˆq 1 solution of

11 H 0 + A) 4 ˆq ) 1 x x ˆq 1 ) =ˆq 1 1) = 0. G. Buscaglia et al. / J. Math. Anal. Appl ) = H 0, x, x We now nee to know the sign of C 1 A) 3 an of 1 ˆq 1x) x. Due to 2.24) we obtain q 1 x = 3 P H 0 + A x + C 1A) H 0 + A) ) an ue to q 1 ) = q 1 1) = 0, we have 1 P x C 1 A) = 3 H 0 + A) 1 H. 2.27) 0 + A) 3 Integrating 2.14) there exists C 2 A) R such that H 0 + A) 3 P x = H 0 + A C 2 A). 2.28) Diviing by H 0 + A) 4 an integrating we have 1 P x H 0 + A) x = 1 H 0 + A) 3 C 2 A) Substituting in 2.27) we obtain 1 C 1 A) 3 = 3C 2 A) H 0 + A) 4 1 H 0 + A). 3 1 H 0 + A) 4. Diviing 2.28) by H 0 + A) 3, integrating in an using P) = P1) = 0 we euce C 2 A) > 0 which gives C 1 A) < 3. Let us now prove by absur that 2.29) H 0 1)<H 0 ). Assuming that max s ],1[ H 0 s) = H 0 1) we euce H0 1) + A ) 1 s 1 H 0 s) + A) 3 s H 0 s) + A) ) which will contraict hypothesis ψ1,a)>0. This implies 2.30). We now apply Lemma 2.9 with ν = 1 an φz)= z + A) 4 an we euce that the solution ˆq 1 of 2.26) satisfies ˆq 1 x > 0. With 2.25) an 2.29) we obtain the result in this case.

12 852 G. Buscaglia et al. / J. Math. Anal. Appl ) Case 2: ψ1,a)<0. We have ψ1,a )<0 for any A in a neighborhoo of A. We euce that P = PA ) satisfies 2.15) with A replace by A an β = βa ) ],1[. From 2.18) we obtain ψ β βa), A) < 0 because βa) >. Then the Implicit Function Theorem implies that the application β is of class C 1 in A. An elementary calculus gives an ψ β β, A) = 3 A s H 0 s) + A) H 0 β) + A ) β β A) = 3 H 0β) + A) β H 0s) + A) 4 s β On the other han we have GA) = β Px)x H 0 β) β H 0s) + A) 3 s s H 0 s) + A) ) H 0s) + A) 3 s. 2.32) with P = PA)given by 2.16) an β = βa). It is clear that G is erivable in A an since Pβ)= 0, an elementary calculus gives G A) = β q 2 x) x with q 2 : [,β] R given by q 2 x) = 3 + H 0 β)β A) ) x H0 s) + A ) 3 s + 3 H 0 β) + A ) x H0 s) + A ) 4 s. Differentiating in x an multiplying by H 0 s) + A) 4 we obtain that q 2 satisfies [ H 0 + A) 4 q ] 2 = [ 3 + H 0 x x β)β A) ] H 0 x, q 2 ) = q 2 β) = 0, where q 2 β) = 0 is a irect consequence of 2.32). By linearity we have as in Case ) 2.34) q 2 = [ 3 + H 0 β)β A) ] ˆq 2, 2.35) where ˆq 2 satisfies [ H 0 + A) 4 ˆq ] 2 = H 0 x, ˆq 2 ) =ˆq 2 β) = 0. x ],β[, 2.36)

13 G. Buscaglia et al. / J. Math. Anal. Appl ) On the other han, from 2.32) we obtain 3 + H 0 β)β A) = 3 H 0 β) + A ) β H 0 + A) 4 β H 0 + A) 3 which implies 3 + H 0 β)β A) > 0. We now use Lemma 2.9 in orer to fin the sign of 1 ˆq 2x) x. Let us first prove by absur that H 0 β) < H 0 ). Assume that max ],β[ H 0 = H 0 β). This implies H0 β) + A ) β s β H 0 s) + A) 3 > s H 0 s) + A) 2 which contraicts the equality that ψβ,a)= 0, which shows 2.38). Applying Lemma 2.9 to problem 2.36) with ν = β an φz)= z + A) 4, we obtain β ˆq 2 x > 0. Combining 2.39), 2.35) an 2.37) the claime result is prove. 2.37) 2.38) 2.39) Proposition If A>0 satisfies ψ1,a)= 0, then ψ A 1,A)>0. Proof. From the efinition 2.17) of ψ we have by elementary calculus ψ 1 A 1,A)= 3 H 0 1) H 0 x) H 0 x) + A) 4 x. Exactly as in 2.30), we prove that H 0 1)<H 0 ), which implies the existence of ν ],[ such that H 0 ν) = H 0 1). We now write ψ ν A 1,A)= 3 φx) 1 H 0 x) + A x + 3 with φx)= H 01) H 0 x) H 0 x)+a) 3. It is clear that ψ1,a)= 0 means 1 Since we have φx)x = 0. ν φx) H 0 x) + A x 2.40)

14 854 G. Buscaglia et al. / J. Math. Anal. Appl ) φx)<0 for x ],ν[ an φx)>0 forx ]ν,1[, we euce ψ A 1,A)> 3 H 0 1) + A ν which gives the result by using 2.40). φx)x + 3 H 0 1) + A 1 ν φx)x We are now in a position to state the main result of this section. Proposition Assume in aition to hypothesis 2.10) that H 0 is non-ecreasing. Then if the function G is continuous, then it is strictly ecreasing. Proof. From Proposition 2.12, we euce that it can exists at most a single point A>0 such that ψ1,a)= 0. Then Proposition 2.11 gives the claime result. 3. The existence result We point out that we want to show the existence of a solution p,a,θ) of the system 1.2) 1.5). We will search the unknown a, θ) in the biimensional set { Q = a, θ) R 2 : θ< h 0 ), a > min h0 x) + θx )}. x Let us introuce the following applications g 1,g 2 : Q R efine by g 1 a, θ) = px F, g 2 a, θ) = xp x Fx 0 with p = pa,θ) the unique solution of the variational inequality 1.2) where h is given by 1.5). It is clear that the original problem 1.2) 1.4) is reuce to fin a, θ) Q solution of the system { g1 a, θ) = 0, 3.1) g 2 a, θ) = 0. The regularity of g 1 an g 2 are not obvious since we cannot use there the Implicit Function Theorem. In the next proposition, we show the continuity of g 1 an g 2 by a irect metho. Proposition 3.1. The functions g 1 an g 2 are continuous on Q. Proof. Let a, θ) in Q an a n,θ n ) a, θ) such that there exists b 0 > 0 with h 0 x) + θ n x + a n b 0, x. 3.2) Let p n = pa n,θ n ). Taking φ = 0 in 1.2) an using 3.2) we euce p n H 1 0 ) C with C inepenent of n.

15 G. Buscaglia et al. / J. Math. Anal. Appl ) We euce the existence of p K an of a subsequence of p n enote also by p n such that p n pin H0 1)-weak. We also have for any ϕ K, Using h 0 + a n + θ n x) 3 p n x ϕ ) 2 h 0 + a n + θ n x) 3 pn + h 0 + a n + θ n x) ϕ p n) x. 3.3) x x h 0 + a n + θ n x h 0 + a + θx in L ) strongly an lim inf n + ) 2 ) h 0 + a n + θ n x) 3 pn p 2 h 0 + a + θx) 3 x x an passing to the limit in 3.3), we obtain p = pa,θ). Since p is unique, we euce that the entire sequence p n converges to p in H0 1 ) weakly which easily gives the result. Now in orer to reuce 3.1) to a single equation, we nee the following result: Proposition 3.2. For any θ< h 0 ) there exists an unique solution a of the equation g 1 a, θ) = 0 such that a, θ) Q. Proof. We istinguish the two cases. Case 1: θ 0. We apply the results of Section 2.1 which H 0 x) = h 0 x) + θx θ, which satisfies hypothesis 2.2) an we classically euce the result as a consequence of Propositions 2.1, 2.5 an Corollary 2.3 with α = 0forθ<0an α = 1forθ = 0. Case 2: θ ]0, h 0 )[. Since h 0 is increasing because h 0 > 0) there exists a unique solution θ ], 1[ of the equation h 0 θ ) = θ 3.4) an we have that the function x h 0 x) + a + θx is ecreasing on ], θ [ an increasing on ] θ, 1[. It is convenient to write hx) = h 0 x) + a + θx in the form hx) = b + rx) with b = h 0 θ ) + a + θ θ 3.5) an rx) = h 0 x) h 0 θ ) h 0 x θ )x θ ). 3.6) It is clear that

16 856 G. Buscaglia et al. / J. Math. Anal. Appl ) r C 2 ), r θ ) = r θ ) = 0, r x) < 0 forx< θ, r x) > 0 forx> θ, r x) = h 0 x) > 0, x. It is also clear that the require problem becomes: show the existence of a unique b>0 an p K, satisfying: b + r) 3 p ϕ p) ϕ p) b + r), x ϕ K, 3.7) p = F. 3.8) It is natural to introuce the application b ]0, + [ G 1 b) = px ]0, + [ with p the solution of 3.7) an the problem 3.7) 3.8) is equivalent to: fin b>0 such that G 1 b) = F. 3.9) In fact G 1 b) = g 1 b h 0 θ ) θ θ,θ).) From 3.5) an Proposition 3.1 the function G 1 is continue since θ is fixe. We apply the results of Section 2.2 with H 0 = r an A = b. Now the result is a consequence of Proposition 2.1, Corollary 2.7 an Proposition Now for any θ< h 0 ), we enote by a = aθ) the unique solution of g 1 a, θ) = 0. We introuce the application S : ], h 0 )[ R efine by ) Sθ) = g 2 aθ),θ an it is clear that the problem 1.2) 1.5) or 3.1) is reuce to search for at least a zero of S. We remark that the restriction of S on ], 0[ is exactly the function S introuce in [1] since the variational inequality becomes equation in this case). We prove in [1] an the result is of course vali here): lim Sθ) > 0. θ The next result proves that S takes a negative value in at least a point. We obtain the following result. 3.10) Proposition 3.3. There exists θ 0 ]0, h 0 )[ such that Sθ 0 )<0. Proof. Using the notations of the proof of Case 2 of Proposition 3.2, we have Sθ) = px x 0 ) with px),b > 0 the unique solution of 3.7) 3.8), with r given by 3.6). 3.11)

17 G. Buscaglia et al. / J. Math. Anal. Appl ) It is clear that if θ h 0 ), then θ an we obtain r) 0. Since there exists θ in a neighborhoo of h 0 such that r)<r1) 3.12) we euce the existence of β θ ], 1[ such that the solution p of 3.7) satisfies r + b) 3 p ) = r x on ],β θ [, 3.13) p) = pβ θ ) = 0, 3.14) p β θ ) = 0, 3.15) px) = 0 on]β θ, 1[ 3.16) if not one woul euce exactly as for 2.38) that r1)<r) which contraict 3.12)). We then have exactly as in 2.38) that rβ θ )<r) which implies that rβ θ ) 0 when θ h 0 ). On the other han, from 3.6), we easily euce 3.17) rx) 1 2 inf h 0 x)x θ) 2 x an the positivity of h 0 implies β θ θ 0. We then euce β when θ h 0 ). Then for any x 0 ], 1[, we can choose θ in a neighborhoo of h 0 ) such that β θ <x 0. Since p = 0forx β θ,wehave 3.18) β θ Sθ) = px x 0 ) an finally the positivity of p on ],β[ an 3.18) give the claime result. It now remains to show the continuity of the function S. Proposition 3.4. The function S is continuous on ], h 0 )[. Proof. The ifficulty of the proof of the continuity of S comes from the change in the nature of problem satisfie by a, p) when θ varies in the mentione interval. We fixe θ ], h 0 )[ an we consier three cases. Case 1: θ 0 < 0. The continuity of S in θ 0 was prove in [1] see Proposition 3.3). Case 2: θ 0 ]0, h 0 )[. Letθ n θ 0 an we can suppose that θ n [0,θ 1 ] with 0 <θ 1 < h 0 ) fixe. We enote by θn the solution of 3.4) with θ = θ n an r n x) = h 0 x) h 0 θn ) h 0 θ n )x θn ). We also enote b n > 0,p n x) the unique solution of 3.7) 3.8) with r = r n.

18 858 G. Buscaglia et al. / J. Math. Anal. Appl ) We prove by absur that there exists a compact ˆK ]0, + [ inepenent of n such that b n ˆK, n N. Let us suppose the contrary, then we have two cases: 3.19) Case 2.1: There exists a subsequence enote by b n such that b n +. Then by Proposition 2.1 with H 0 = r n an A = b n, we euce that p n x 0 which contraicts p n x = F. Case 2.2: There exits a subsequence enote by b n such that b n 0. We apply Proposition 2.6 with H 0 = r n, 0 = θ1 an A = b n notice that r n x) = h 0 x) is inepenent of n) an we euce that p n x + which also contraicts p n = F. Then we prove 3.19). Then there is a subsequence enote also by b n an an element b ˆK such that b n b. By a similar proof to that of Proposition 3.1 we euce the existence of a subsequence enote by p n an of an element p K such that p n p weakly in H 1 ) with p an b satisfying 3.7) 3.8) with rx) = h 0 x) h 0 θ0 ) h 0 θ 0 )x θ0 ) since r n r strongly in L )). By uniqueness of the solution of the 3.7) 3.8), we euce that all the sequence b n,p n ) converges to the limit. Then Sθ n ) Sθ 0 ). Case 3: θ 0 = 0. We also take a sequence θ n 0 an we can consier two cases: Case 3.1: θ n > 0. We show exactly as in Case 2 that Sθ n ) S0). Case 3.2: θ n < 0. We are in the hypothesis of Section 2.1, so Sθ n ) = p n x x 0 )x with b n > 0,p n x) solution of [ bn + h 0 x) θ n 1 x) ) ] 3 p n = [ h0 x) θ n 1 x) ], x p n ) = p n 1) = 0, p n x = F. Exactly as in Case 2.1 we prove that no subsequence b n can exist. Let us now suppose that a subsequence of b n satisfies b n 0. We apply Lemma 2.4 with H 0 x) = h 0 x) θ n 1 x) an A = b n.wehave H 0 x) = θ n h 0 x) an the assumptions of Lemma 2.4 are satisfie with M 0 = min x h 0, M 1 = max x h 0, δ 0 = 1 an Θ = θ n. Then we obtain p n x + which contraicts p n x = F. We then prove that there is no subsequence of b n which tens to 0. We euce the existence of a compact ˆK ]0, + [ such that b n ˆK, n N. The en of the proof is as in Case 2.

19 G. Buscaglia et al. / J. Math. Anal. Appl ) Now applying the intermeiate value theorem, we obtain from 3.10), Propositions 3.3 an 3.4 an the main result of this paper. Theorem 3.5. For any x 0 ], 1[ an F>0 there exists at least a solution p,a,θ) of 1.2) 1.5). Acknowlegments This work was partially supporte by grant A03E01 from ECOS-SECYT an by the CMIMF Action Intérée N o MA/04/94. References [1] G. Buscaglia, I. Ciuperca, I. Hafii, M. Jai, Existence of equilibria in articulate bearings, J. Math. Anal. Appl ) [2] I. Ciuperca, I. Hafii, M. Jai, Singular perturbation problem for the incompressible Reynols equation, Electron. J. Differential Equations ) 2006) [3] G. Bayaa, M. Chambat, Sur quelques moélisations e la zone e cavitation en lubrification hyroynamique Some moelling of the cavitation region in hyroynamic lubrication), J. Méc. Théor. Appl )

Equilibrium analysis for a mass-conserving model in presence of cavitation

Equilibrium analysis for a mass-conserving model in presence of cavitation Equilibrium analysis for a mass-conserving model in presence of cavitation Ionel S. Ciuperca CNRS-UMR 58 Université Lyon, MAPLY,, 696 Villeurbanne Cedex, France. e-mail: ciuperca@maply.univ-lyon.fr Mohammed

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

Topological Sensitivity Analysis for Three-dimensional Linear Elasticity Problem

Topological Sensitivity Analysis for Three-dimensional Linear Elasticity Problem Topological Sensitivity Analysis for Three-imensional Linear Elasticity Problem A.A. Novotny, R.A. Feijóo, E. Taroco Laboratório Nacional e Computação Científica LNCC/MCT, Av. Getúlio Vargas 333, 25651-075

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

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

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

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

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

Uniqueness of limit cycles of the predator prey system with Beddington DeAngelis functional response

Uniqueness of limit cycles of the predator prey system with Beddington DeAngelis functional response J. Math. Anal. Appl. 290 2004 113 122 www.elsevier.com/locate/jmaa Uniqueness of limit cycles of the preator prey system with Beington DeAngelis functional response Tzy-Wei Hwang 1 Department of Mathematics,

More information

SINGULAR PERTURBATION PROBLEM FOR THE INCOMPRESSIBLE REYNOLDS EQUATION

SINGULAR PERTURBATION PROBLEM FOR THE INCOMPRESSIBLE REYNOLDS EQUATION Electronic Journal of Differential Equations, Vol. 20062006), No. 83, pp. 1 19. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu login: ftp) SINGULAR

More information

Exponential asymptotic property of a parallel repairable system with warm standby under common-cause failure

Exponential asymptotic property of a parallel repairable system with warm standby under common-cause failure J. Math. Anal. Appl. 341 (28) 457 466 www.elsevier.com/locate/jmaa Exponential asymptotic property of a parallel repairable system with warm stanby uner common-cause failure Zifei Shen, Xiaoxiao Hu, Weifeng

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

Linear and quadratic approximation

Linear and quadratic approximation Linear an quaratic approximation November 11, 2013 Definition: Suppose f is a function that is ifferentiable on an interval I containing the point a. The linear approximation to f at a is the linear function

More information

State-Space Model for a Multi-Machine System

State-Space Model for a Multi-Machine System State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal

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

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

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

Journal of Algebra. A class of projectively full ideals in two-dimensional Muhly local domains

Journal of Algebra. A class of projectively full ideals in two-dimensional Muhly local domains Journal of Algebra 32 2009 903 9 Contents lists available at ScienceDirect Journal of Algebra wwwelseviercom/locate/jalgebra A class of projectively full ieals in two-imensional Muhly local omains aymon

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

Control Volume Derivations for Thermodynamics

Control Volume Derivations for Thermodynamics Control olume Derivations for Thermoynamics J. M. Powers University of Notre Dame AME 327 Fall 2003 This ocument will give a summary of the necessary mathematical operations necessary to cast the conservation

More information

Convergence rates of moment-sum-of-squares hierarchies for optimal control problems

Convergence rates of moment-sum-of-squares hierarchies for optimal control problems Convergence rates of moment-sum-of-squares hierarchies for optimal control problems Milan Kora 1, Diier Henrion 2,3,4, Colin N. Jones 1 Draft of September 8, 2016 Abstract We stuy the convergence rate

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

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

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

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

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

Students need encouragement. So if a student gets an answer right, tell them it was a lucky guess. That way, they develop a good, lucky feeling.

Students need encouragement. So if a student gets an answer right, tell them it was a lucky guess. That way, they develop a good, lucky feeling. Chapter 8 Analytic Functions Stuents nee encouragement. So if a stuent gets an answer right, tell them it was a lucky guess. That way, they evelop a goo, lucky feeling. 1 8.1 Complex Derivatives -Jack

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

Lecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b.

Lecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b. b. Partial erivatives Lecture b Differential operators an orthogonal coorinates Recall from our calculus courses that the erivative of a function can be efine as f ()=lim 0 or using the central ifference

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

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

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

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

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

FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS. 1. Introduction

FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS. 1. Introduction FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS ALINA BUCUR, CHANTAL DAVID, BROOKE FEIGON, MATILDE LALÍN 1 Introuction In this note, we stuy the fluctuations in the number

More information

x f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0.

x f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0. Engineering Mathematics 2 26 February 2014 Limits of functions Consier the function 1 f() = 1. The omain of this function is R + \ {1}. The function is not efine at 1. What happens when is close to 1?

More information

Math 251 Notes. Part I.

Math 251 Notes. Part I. Math 251 Notes. Part I. F. Patricia Meina May 6, 2013 Growth Moel.Consumer price inex. [Problem 20, page 172] The U.S. consumer price inex (CPI) measures the cost of living base on a value of 100 in the

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

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

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

QF101: Quantitative Finance September 5, Week 3: Derivatives. Facilitator: Christopher Ting AY 2017/2018. f ( x + ) f(x) f(x) = lim

QF101: Quantitative Finance September 5, Week 3: Derivatives. Facilitator: Christopher Ting AY 2017/2018. f ( x + ) f(x) f(x) = lim QF101: Quantitative Finance September 5, 2017 Week 3: Derivatives Facilitator: Christopher Ting AY 2017/2018 I recoil with ismay an horror at this lamentable plague of functions which o not have erivatives.

More information

Qubit channels that achieve capacity with two states

Qubit channels that achieve capacity with two states Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March

More information

II. First variation of functionals

II. First variation of functionals II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent

More information

ELEC3114 Control Systems 1

ELEC3114 Control Systems 1 ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.

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

A Weak First Digit Law for a Class of Sequences

A Weak First Digit Law for a Class of Sequences International Mathematical Forum, Vol. 11, 2016, no. 15, 67-702 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1288/imf.2016.6562 A Weak First Digit Law for a Class of Sequences M. A. Nyblom School of

More information

Lecture 7: Interchange of integration and limit

Lecture 7: Interchange of integration and limit Lecture 7: Interchange of integration an limit Differentiating uner an integral sign To stuy the properties of a chf, we nee some technical result. When can we switch the ifferentiation an integration?

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

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

A global Implicit Function Theorem without initial point and its applications to control of non-affine systems of high dimensions

A global Implicit Function Theorem without initial point and its applications to control of non-affine systems of high dimensions J. Math. Anal. Appl. 313 (2006) 251 261 www.elsevier.com/locate/jmaa A global Implicit Function Theorem without initial point an its applications to control of non-affine systems of high imensions Weinian

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

Slide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13)

Slide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13) Slie10 Haykin Chapter 14: Neuroynamics (3r E. Chapter 13) CPSC 636-600 Instructor: Yoonsuck Choe Spring 2012 Neural Networks with Temporal Behavior Inclusion of feeback gives temporal characteristics to

More information

WUCHEN LI AND STANLEY OSHER

WUCHEN LI AND STANLEY OSHER CONSTRAINED DYNAMICAL OPTIMAL TRANSPORT AND ITS LAGRANGIAN FORMULATION WUCHEN LI AND STANLEY OSHER Abstract. We propose ynamical optimal transport (OT) problems constraine in a parameterize probability

More information

Polynomial Inclusion Functions

Polynomial Inclusion Functions Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl

More information

On the enumeration of partitions with summands in arithmetic progression

On the enumeration of partitions with summands in arithmetic progression AUSTRALASIAN JOURNAL OF COMBINATORICS Volume 8 (003), Pages 149 159 On the enumeration of partitions with summans in arithmetic progression M. A. Nyblom C. Evans Department of Mathematics an Statistics

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

12.11 Laplace s Equation in Cylindrical and

12.11 Laplace s Equation in Cylindrical and SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering

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

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

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite

More information

The Ehrenfest Theorems

The Ehrenfest Theorems The Ehrenfest Theorems Robert Gilmore Classical Preliminaries A classical system with n egrees of freeom is escribe by n secon orer orinary ifferential equations on the configuration space (n inepenent

More information

Extreme Values by Resnick

Extreme Values by Resnick 1 Extreme Values by Resnick 1 Preliminaries 1.1 Uniform Convergence We will evelop the iea of something calle continuous convergence which will be useful to us later on. Denition 1. Let X an Y be metric

More information

REPRESENTATIONS FOR THE GENERALIZED DRAZIN INVERSE IN A BANACH ALGEBRA (COMMUNICATED BY FUAD KITTANEH)

REPRESENTATIONS FOR THE GENERALIZED DRAZIN INVERSE IN A BANACH ALGEBRA (COMMUNICATED BY FUAD KITTANEH) Bulletin of Mathematical Analysis an Applications ISSN: 1821-1291, UL: http://www.bmathaa.org Volume 5 Issue 1 (2013), ages 53-64 EESENTATIONS FO THE GENEALIZED DAZIN INVESE IN A BANACH ALGEBA (COMMUNICATED

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

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

Notes on Lie Groups, Lie algebras, and the Exponentiation Map Mitchell Faulk

Notes on Lie Groups, Lie algebras, and the Exponentiation Map Mitchell Faulk Notes on Lie Groups, Lie algebras, an the Exponentiation Map Mitchell Faulk 1. Preliminaries. In these notes, we concern ourselves with special objects calle matrix Lie groups an their corresponing Lie

More information

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

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

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

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

model considered before, but the prey obey logistic growth in the absence of predators. In

model considered before, but the prey obey logistic growth in the absence of predators. In 5.2. First Orer Systems of Differential Equations. Phase Portraits an Linearity. Section Objective(s): Moifie Preator-Prey Moel. Graphical Representations of Solutions. Phase Portraits. Vector Fiels an

More information

Monotonicity for excited random walk in high dimensions

Monotonicity for excited random walk in high dimensions Monotonicity for excite ranom walk in high imensions Remco van er Hofsta Mark Holmes March, 2009 Abstract We prove that the rift θ, β) for excite ranom walk in imension is monotone in the excitement parameter

More information

arxiv: v2 [math.cv] 2 Mar 2018

arxiv: v2 [math.cv] 2 Mar 2018 The quaternionic Gauss-Lucas Theorem Riccaro Ghiloni Alessanro Perotti Department of Mathematics, University of Trento Via Sommarive 14, I-38123 Povo Trento, Italy riccaro.ghiloni@unitn.it, alessanro.perotti@unitn.it

More information

Permanent vs. Determinant

Permanent vs. Determinant Permanent vs. Determinant Frank Ban Introuction A major problem in theoretical computer science is the Permanent vs. Determinant problem. It asks: given an n by n matrix of ineterminates A = (a i,j ) an

More information

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1 Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems

More information

Discrete Mathematics

Discrete Mathematics Discrete Mathematics 309 (009) 86 869 Contents lists available at ScienceDirect Discrete Mathematics journal homepage: wwwelseviercom/locate/isc Profile vectors in the lattice of subspaces Dániel Gerbner

More information

A LIMIT THEOREM FOR RANDOM FIELDS WITH A SINGULARITY IN THE SPECTRUM

A LIMIT THEOREM FOR RANDOM FIELDS WITH A SINGULARITY IN THE SPECTRUM Teor Imov r. ta Matem. Statist. Theor. Probability an Math. Statist. Vip. 81, 1 No. 81, 1, Pages 147 158 S 94-911)816- Article electronically publishe on January, 11 UDC 519.1 A LIMIT THEOREM FOR RANDOM

More information

Partial Differential Equations

Partial Differential Equations Chapter Partial Differential Equations. Introuction Have solve orinary ifferential equations, i.e. ones where there is one inepenent an one epenent variable. Only orinary ifferentiation is therefore involve.

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

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

Two formulas for the Euler ϕ-function

Two formulas for the Euler ϕ-function Two formulas for the Euler ϕ-function Robert Frieman A multiplication formula for ϕ(n) The first formula we want to prove is the following: Theorem 1. If n 1 an n 2 are relatively prime positive integers,

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

Sharp Thresholds. Zachary Hamaker. March 15, 2010

Sharp Thresholds. Zachary Hamaker. March 15, 2010 Sharp Threshols Zachary Hamaker March 15, 2010 Abstract The Kolmogorov Zero-One law states that for tail events on infinite-imensional probability spaces, the probability must be either zero or one. Behavior

More information

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7. Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency

More information

arxiv: v2 [math.dg] 16 Dec 2014

arxiv: v2 [math.dg] 16 Dec 2014 A ONOTONICITY FORULA AND TYPE-II SINGULARITIES FOR THE EAN CURVATURE FLOW arxiv:1312.4775v2 [math.dg] 16 Dec 2014 YONGBING ZHANG Abstract. In this paper, we introuce a monotonicity formula for the mean

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

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

Darboux s theorem and symplectic geometry

Darboux s theorem and symplectic geometry Darboux s theorem an symplectic geometry Liang, Feng May 9, 2014 Abstract Symplectic geometry is a very important branch of ifferential geometry, it is a special case of poisson geometry, an coul also

More information

(a 1 m. a n m = < a 1/N n

(a 1 m. a n m = < a 1/N n Notes on a an log a Mat 9 Fall 2004 Here is an approac to te eponential an logaritmic functions wic avois any use of integral calculus We use witout proof te eistence of certain limits an assume tat certain

More information

LEGENDRE TYPE FORMULA FOR PRIMES GENERATED BY QUADRATIC POLYNOMIALS

LEGENDRE TYPE FORMULA FOR PRIMES GENERATED BY QUADRATIC POLYNOMIALS Ann. Sci. Math. Québec 33 (2009), no 2, 115 123 LEGENDRE TYPE FORMULA FOR PRIMES GENERATED BY QUADRATIC POLYNOMIALS TAKASHI AGOH Deicate to Paulo Ribenboim on the occasion of his 80th birthay. RÉSUMÉ.

More information

arxiv: v1 [math.co] 15 Sep 2015

arxiv: v1 [math.co] 15 Sep 2015 Circular coloring of signe graphs Yingli Kang, Eckhar Steffen arxiv:1509.04488v1 [math.co] 15 Sep 015 Abstract Let k, ( k) be two positive integers. We generalize the well stuie notions of (k, )-colorings

More information

Spectral properties of a near-periodic row-stochastic Leslie matrix

Spectral properties of a near-periodic row-stochastic Leslie matrix Linear Algebra an its Applications 409 2005) 66 86 wwwelseviercom/locate/laa Spectral properties of a near-perioic row-stochastic Leslie matrix Mei-Qin Chen a Xiezhang Li b a Department of Mathematics

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

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

On state representations of time-varying nonlinear systems

On state representations of time-varying nonlinear systems On state representations of time-varying nonlinear systems Paulo Sérgio Pereira a Silva a, Simone Batista a, a University of São Paulo, Escola Politécnicca PTC Av. Luciano Gualberto trav. 03, 158, 05508-900

More information

On colour-blind distinguishing colour pallets in regular graphs

On colour-blind distinguishing colour pallets in regular graphs J Comb Optim (2014 28:348 357 DOI 10.1007/s10878-012-9556-x On colour-blin istinguishing colour pallets in regular graphs Jakub Przybyło Publishe online: 25 October 2012 The Author(s 2012. This article

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

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

Unit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule

Unit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule Unit # - Families of Functions, Taylor Polynomials, l Hopital s Rule Some problems an solutions selecte or aapte from Hughes-Hallett Calculus. Critical Points. Consier the function f) = 54 +. b) a) Fin

More information