1. TEMEL TANIM VE EŞİTSİZLİKLER
Fizik, mühendislik ve bir çok bilim dallarında karşılaşılan problemleri çözmek için önce problemler matematiksel ifadelerle formüle edilir. Matematiksel modelleri kurulan problemlerin daha sonra uygun teknik ve yöntemlerle çözümleri araştırılır. Bu kısımda, genel tanımını vereceğimiz diferensiyel denklemler teorik ve pratik açıdan oldukça önemli olup tüm bilim dallarında çok geniş uygulama alanlarına sahiptir [1-4]. Diferensiyel denklemler tanımını vermeden önce analiz derslerinden bildiğimiz adi ve kısmi türevleri kısaca verelim.
x
bağımsız değişken, y de bu bağımsız değişkene bağlı bir değişken yani bağımlı değişken olsun. Bu durum y = f(x) olarak yazılır. y ifadesinin bağımsız değişkenlere göre çeşitli mertebeden adi türevleridx dy y ≡′ , 2 2 dx y d y ≡′′ , … ile gösterilir.
x
ve ybağımsız değişkenler, z de bu bağımsız değişkenlere bağlı bir değişken olsun. Bu durum z = f(x,y) olarak yazılır. z ifadesinin bağımsız değişkenlere göre çeşitli mertebeden kısmi türevlerix z zx ∂ ∂ ≡ , y z zy ∂ ∂ ≡ , 2 2 x z zxx ∂ ∂ ≡ , …
şeklinde ifade edilir.
Tanım 1.1.
Bir bağımlı değişkenin, bir veya daha fazla bağımsız değişkene göre çeşitli mertebeden türevlerini ihtiva eden bir denkleme diferensiyel denklem adı verilir.
Tanım 1.2.
Bir bağımlı değişkenin, bir tek bağımsız değişkene göre çeşitli mertebeden (adi) türevlerini içeren diferensiyel denkleme adi türevli diferensiyel denklem denir. Bir diferensiyel denklemin mertebesi, denklemde görülen en yüksek mertebeden türevin mertebesidir. n. mertebeden adi diferensiyel bir denklem genel olarak,
F(x,y,y′,...,y(n))=0 (1.1) kapalı şekilde gösterilebilir [5]. Bir a<x<b aralığında tanımlı
φ
fonksiyonu a<x<bF(x,
φ
(x),φ
′(x),...,φ
(n)(x))=0 (1.2)ise
φ
fonksiyonuna (1.1) denkleminin çözümüdür denir. Bir adi diferensiyel denklemin genel çözümü, diferensiyel denklemin mertebesi kadar sabit değeri parametre olarak kabul eden bir eğri ailesi olarak ortaya çıkar. Genel çözümdeki keyfi sabit veya sabitlerin belirli değerleri için elde edilen çözüme adi diferensiyel denklemin özel çözümü denir [6].Tanım 1.3.
Bir bağımlı değişkenin, en az iki bağımsız değişkene göre kısmi türevlerini içeren denkleme kısmi türevli diferensiyel denklem denir.
) , , , (x y z t u
u = fonksiyonu x,y,z ve t bağımsız değişkenleriyle bir D bölgesinde
tanımlı bir fonksiyon olsun. u fonksiyonunun x,y,z ve t değişkenlerine göre kısmi türevlerini içeren .n mertebeden bir kısmi diferensiyel denklem genel olarak,
{
)
0
,...,
,
,
,
,
,
,
,
,
,
,
(
tan ... , ,=
− e n t ttt xy tt zz yy xx t z y xu
u
u
u
u
u
u
u
u
u
t
z
y
x
F
(1.3) şeklinde gösterilir.Bir kısmi diferensiyel denklemdeki değişkenler Γ sınırına sahip bir D bölgesinde tanımlanır. t zaman değişkeni olmak üzere
t
1<
t
<
t
2 aralığında ve D bölgesinde (x,y,z)noktasında
ϕ
fonksiyonu n mertebeye kadar türevleri sürekli ise .ϕ
=ϕ
(x,y,z,t)fonksiyonuna n mertebeden kısmi diferensiyel denklemin çözümüdür denir [7]. Kısmi . diferensiyel denklemin genel çözümü, denklemin mertebesi kadar keyfi fonksiyon içerir ve çözümleri bir yüzey ailesi olarak ortaya çıkar. Bundan dolayı, adi türevli diferensiyel denklemlere kıyasla kısmi diferensiyel denklemlerin çözümlerini bulmak daha zordur. Bir kısmi diferensiyel denklemin özel çözümünün bulunması için bazı yardımcı şartlar gerekir. Bu şartlar genel olarak iki başlık altında toplanabilir.
)
(i Sınır şartları: Sınırlı bir D tanım bölgesinde, bir fiziksel olayın matematiksel durumunu kontrol eden bir kısmi diferensiyel denklemde u bağımlı değişkeni D tanım bölgesinin sınırında yazılır. Yani sınır şartları kısmi diferensiyel denklemin sağlandığı D bölgesinin Γ sınırı boyunca sağlanması gereken şartlardır. Sınır şartlarını aşağıda tanımlanan üç farklı şekil ile verebiliriz:
a) Dirichlet Sınır Şartları:
Bu durumda, u fonksiyonu D bölgesinin Γ sınırı üzerinde yazılır. 0< x <l olmak üzere l uzunluğundaki bir çubuk için sınır şartları
α
veβ
olmak üzere u(0)=α
, u(l)=β
şeklinde tanımlanır.
1
0
<
x <
l
ve 0< y <l2 olmak üzere bir dikdörtgensel tabak için sınır şartları1
) , 0
( y =α
u , u(l1,y)=α2,
u
(
x
,
0
)
=
β
1 ve u(x,l2)=β2 olarak yazılır. Eğer u bağımlı değişkeni sınırlar üzerindeki herhangi bir noktada sıfır ise sınır şartları homojen olarak adlandırılır. Aksi halde homojen değildir.b) Neumann Sınır Şartları:
Bu durumda, sınır şartları u nun normalleri boyunca yani u nun türevinde yazılır. l
uzunluğundaki bir çubuk için Neumann sınır şartları; ux( t0, )=
α
, ux( tl, )=β
şeklindedir.c) Karışık (mixed) veya Robin Sınır Şartları:
Bu durumda, u fonksiyonu ve aynı zamanda u nun .n türevi de sınırlar üzerinde belirlenir. Ve genellikle 2-boyutlu Laplace denklemlerinde kullanılır.
)
(ii Başlangıç şartları: Başlangıç şartları sistemin başlangıcında D bölgesi boyunca sağlanması gereken şartlardır. Yani bağımlı değişkenin kendisi veya türevleri cinsinden bağımsız değişkenlerin bir tek değeri için verilmiş ise bu şartlara başlangıç şartları denir.
Başlangıç şartlarıyla birlikte verilmiş kısmi diferensiyel denkleme ‘Cauchy problemi’ denir. Örneğin, 0 2 = − xx tt c u u , t>0; ) ( ) 0 , (x f x u = , ut(x,0)= g(x), ikinci mertebeden bir Cauchy problemidir.
İkinci mertebeden, iki bağımsız değişkenli bir kısmi diferensiyel denklem,
Auxx +Buxy +Cuyy +Dux +Euy +Fu+G=0 (1.4) genel şekliyle verilebilir. Burada A,B,C,D,E,F katsayı fonksiyonları ve G fonksiyonu da sabit veya değişken içeren fonksiyondur. Diğer taraftan ∆=B2 −4AC diskriminantını tanımlayalım [8].
Diskriminant Denklem Tipi Örnek İsimlendirme
∆>0 Hiperbolik utt −c2uxx =0 Dalga Denklemi
∆=0 Parabolik ut −kuxx =0 Isı Denklemi
<
Herhangi bir tipteki problemin çözümü, klasik Hadamard testi gereğince aşağıdaki üç şartı sağlarsa problem, “iyi durumlu”, en az bir şartı sağlamaz ise “kötü durumlu” olarak adlandırılır. Bu şartlar;
i. Varlık: En az bir çözüm vardır, ii. Teklik: En çok bir çözüm vardır,
iii. Kararlılık: Çözüm verilere sürekli bağlıdır.
Pratikte bir denklemin çözümünün varlığını tarif etmenin en iyi yolu problemdeki bütün şartları sağlayan ve problemde yerine konduğunda denklemi sağlayan bir çözüm yapılandırmaktır. Eğer çözümün tekliği gösterilirse denklemin çözümü bulunmuş demektir.
Birincisi açık bir mantıksal şart olup, ancak sadece fiziksel problemin bir çözüme sahip olması nedeniyle, matematiksel problemin bir çözüme sahip olacağını kolayca ifade edemeyeceğimiz akılda tutulmalıdır. Fiziksel problem bir tek çözüme sahip olabilir, ancak matematiksel problemin birden çok çözümü de olabilmektedir.
Son şart bir gerek şarttır. Uygulamada, ölçüm yönteminde küçük hatalar oluşabilir. Dolayısıyla, fiziksel olayı gösteren matematiksel problem için, verilerdeki küçük bir değişme, çözümde en fazla küçük bir değişikliğe yol açabilir.
Örnek olarak, 0 = + yy xx u u (1.5)
Laplace denklemini n>0 olmak üzere,
u(0,y)=0, y
n y
ux(0, )= 1sin (1.6) Cauchy verileriyle göz önüne alalım. Bu problemin değişkenlere ayrılma yöntemi ile elde edilen çözümü, nx y n y x u1( , )= 12 sinh sin (1.7)
şeklindedir. Başlangıç verileri u(0,y)=0 ve ux(0,y)=0 iken problemin çözümü
u
2=
0
aşikar çözümüdür. İki başlangıç verisi arasındaki fark
n
→
∞
iken0 sin lim −1 = ∞ → n y n
olur. Yani, başlangıç verisinde küçük bir değişiklik olmuştur. Bu başlangıç verisine karşılık gelen çözümler farkının
2
π
=y noktasında,
n
tek bir pozitif sayı olmak üzeren
→
∞
iken limit değerine bakalım;∞
=
−
=
=
−
− ∞ → ∞ → ∞ → 1 2 2 22
lim
sinh
1
lim
)
2
,
(
)
2
,
(
lim
n
e
e
nx
n
x
u
x
u
nx nx n n nπ
π
olur. Yani, başlangıç verilerinde yapılan küçük bir değişiklik çözümde büyük bir değişikliğe yol açmıştır. Böylece (1.5) ve (1.6) Cauchy problemi iyi durumlu olmadığı sonucuna varılır [10].
Tanım 1.4.
Ω, IRn de bir bölge ve ppozitif bir reel sayı olsun. Ω üzerinde tanımlı, ölçülebilen ve
∫
Ω ∞ < dx x u( ) p , x∈Ω (1.8)özelliğini sağlayan tüm u(x) fonksiyonlarının uzayına Lp(Ω) uzayı denir. Lp(Ω) uzayı bir Banach uzayıdır ve bu uzayda norm,
p p L u x dx u p 1 ) ( ( ) =
∫
Ω Ω (1.9) şeklindedir. Tanım 1.5.Ω bölgesinde ölçülebilen bir
u
fonksiyonunu göz önüne alalım. Hemen hemen herΩ ∈
x için
u
(
x
)
≤
K
olacak şekilde bir K sabiti bulunabiliyorsau
fonksiyonuna hemen hemen her yerde sınırlı denir. Böyle K ların en büyük alt sınırlarına dau
fonksiyonunun Ω bölgesindeki esas supremumu denir veess
sup
x Ω∈u
(
x
)
ile gösterilir. Böylece L∞(Ω) ile Ω bölgesinde tanımlı bütün hemen hemen her yerde sınırlı fonksiyonların sınıfı gösterilebilir ve⋅ normu aşağıdaki gibi tanımlanır:
sup ( ) ) ( ess u x u x L∞ Ω = ∈Ω . (1.10) Tanım 1.6. n IR ⊂ Ω
Ω, 1 iki bölge ve
Ω
1⊂
Ω
olsun. Eğer Ω1 ⊂Ω iseΩ
1 bölgesine, Ω bölgesinin kesin alt bölgesi denir veΩ
1⊂⊂
Ω
şeklinde ifade edilir.Tanım 1.7.
) (Ω p loc
L aşağıdaki gibi tanımlanır:
{
: ( ),}
. ) (Ω = ∈ p Ω1 ∀Ω1⊂⊂Ω p loc u u L LTanım 1.8.
Ω, IRn de bir bölge ve
u
da Ω bölgesinde tanımlı bir fonksiyon olsun.u
fonksiyonunun desteği aşağıdaki biçimde tanımlanır:
{
x : u(x) 0}
u
supp = ∈Ω ≠ .
Eğer suppu⊂⊂Ω ise u fonksiyonuna Ω bölgesinde kompakt dayanaklı fonksiyon denir.
Tanım 1.9.
Ω, IRn de açık bir bölge ve m pozitif bir tamsayı olsun. Cm(Ω) ile Ω bölgesinde tanımlı ve
α
≤
m
mertebesinden sürekli kısmi türevleri olan fonksiyonların uzayı gösterilir.) (Ω ∞
C ile Ω bölgesinde tanımlı olan ve her mertebeden sürekli kısmi türevleri olan
fonksiyonların sınıfı gösterilir.
) (
0 Ω
∞
C ile de, C∞(Ω) sınıfından olan kompakt destekli fonksiyonların sınıfı belirtilir.
Tanım 1.10.
i
α
, i=1,2,...,n pozitif tamsayılar olsun.α
=(α
1,...,α
n) sıralı n-lisine çoklu indis denir.α
ileα
’nın boyu yani∑
= = n i i 1α
α
ile gösterilir. 1≤i≤n olmak üzere;
n n x x x D α α α α ∂ ∂ ∂ ∂ ∂ ∂ = ... 2 1 2 1
α
mertebesinden diferansiyel operatörüdür ve u∈Cα(Ω) ise,n n x x x u x u D α α α α ∂ ∂ ∂ = ... ) ( ) ( 1 1
şeklinde ifade edilir.
Tanım 1.11.
) ( ,v∈L1loc Ω
u olsun. Eğer ∀
φ
∈C0∞(Ω) fonksiyonu için,∫
∫
Ω Ω − = v x x dx dx x D x u( ) αφ( ) ( 1)α ( )φ( )oluyorsa
v
fonksiyonunau
fonksiyonun Ω bölgesindeα
mertebesinden Dαu tipinden genelleştirilmiş türevi denir.Tanım 1.12.
Ω, IRn de bir bölge, m pozitif bir sayı ve 1≤ p<∞ olsun. Dαu,
u
fonksiyonunun genelleştirilmiş türevi olmak üzere,{
u
L
D
u
L
m
}
W
m,p(
Ω
)
=
∈
p(
Ω
)
:
α
∈
p(
Ω
),
∀
0
≤
α
≤
ile tanımlanan Wm,p(Ω) uzayına Sobolev uzayı denir. Bu uzayda norm,
p m p L Wmp D u p u 1 ) ( ) ( , =
∑
≤ Ω Ω α α , 1≤ p<∞şeklindedir. Sobolev uzayı bir Banach uzayıdır. Eğer p=2 ise bir Hilbert uzaydır.
1.1. Bazı Önemli Eşitsizlikler
Tanım 1.13. (Cauchy Eşitsizliği):
) ( ,g∈ L2 Ω f ise, fg∈ L1(Ω) ’dır ve ) ( ) ( ) ( 2 2 1Ω ≤ Ω Ω L L L f g g f dır.
Tanım 1.14. (Hölder Eşitsizliği):
1 > p ve 1 +1 =1 q p olsun. ∈ (Ω), p L f g∈Lq(Ω) ise fg∈ L1(Ω)’dır ve ) ( ) ( ) ( 1Ω ≤ LpΩ LqΩ L f g g f dır.
Tanım 1.15. (
ε
-Young Eşitsizliği):0 >
q b p a ab q q p p ε ε + 1 ≤ eşitsizliği sağlanır.
Tanım 1.16. (Genelleştirilmiş Young Eşitsizliği):
0 > i p , i=1,2,...,n ve
∑
==
n i 1p
i1
1
olduğunu kabul edelim. Bu durumda verilen
, IR ai∈ ∀ i=1,2,...,n sayıları için
∑
= ≤ n i i p i n p a a a a i 1 2 1 ... eşitsizliği sağlanır.Tanım 1.17. (Poincare Eşitsizliği):
Ω kümesi IRn ’de sınırlı bir bölge ve ∂Ω da C1 sınıfından olsun. 1≤ p<∞ olduğunu kabul edelim. Bu durumda sadece
n
ve p sayıları ile Ω bölgesine bağlı bir Csayısı vardır öyle ki, her u∈W1,p(Ω) fonksiyonu için
) (Ω Ω ≤ +
∫
p L Du C udx u eşitsizliği sağlanır.Tanım 1.18. (Friedrichs Eşitsizliği):
Ω sınırlı bir bölge olsun. Öyle bir c(Ω)>0 sabiti vardır ki ∀u∈W01,2(Ω) için
) ( ) ( 2 2 Ω ≤ (Ω) Ω L L c Du u eşitsizliği sağlanır.
Tanım 1.19. (Wirtinger Eşitsizliği): ] , 0 [ 2 l L
w∈ olsun. λ=l2 π2 olmak üzere
∫
≤∫
l l x x dx w dx x w 0 0 2 2( ) λ ( ) eşitsizliği sağlanır.2. Dördüncü Mertebeden Doğrusal Olmayan Bir Dalga Denkleminin Çözümlerinin Davranışı
Bu bölümde, sondaj borusu dinamiğini belirten dördüncü mertebeden doğrusal olmayan bir dalga denklemi ele alınmıştır. Sondaj borusu, bir gemi ya da platformdan deniz tabanına uzanıp delme operasyonlarında kullanılan ve matkap kordonu ile delme işleminde ortaya çıkan çamuru zapteden ince, uzun, dikey bir borudur. Amacımız aşağıdaki problemin sıfır çözümünün global asimptotik olarak kararlı olduğunu göstermektir.
[
]
∈ = = = = = + + − + + (2.2) . , 0 ) , ( ) , ( ) , 0 ( ) , 0 ( (2.1) , 0 ) ( IR t t l u t l u t u t u u bu u u x a ku mu xx xx p t t tx x x xxxx tt γBu problemde m, k, b ve pverilen pozitif sayılardır.
γ
ise bir reel sayıdır. a(.) fonksiyonu ise alttan sınırlı ve C1[ ]
0,l sınıfındandır [11-15].2.1. Sıfır Çözümünün Global Asimptotik Kararlığı:
Eğer sıfır çözümü kararlı ve tüm çözümler t→∞ iken sıfıra gidiyor ise, sıfır çözümü global asimptotik kararlı olarak adlandırılır. (2.1) ve (2.2) probleminin tüm çözümleri klasik çözümlerdir. Esas sonucumuz aşağıdaki teoremdir [16].
Teorem 2.1.1
Aşağıdaki şartların sağlandığını varsayalım:
i. m, k, b ve pverilen pozitif sayılar ve
γ
ise bir reel sayıdır.ii. a(.) fonksiyonu C1
[ ]
0,l sınıfındadır ve∀
x
∈
[ ]
0
,
l
için ( ) , 0 0 0 ≥ − ≥ c c x a ve 0 0 2 2 0 = > −c l d kπ
’dır.Bu durumda (2.1) probleminin (2.2) sınır şartını sağlayan sıfır çözümü,
2 1 0 2 2 1 ( ) ) , ( + =
∫
l xx t t u u dx u unormuna göre global asimptotik kararlıdır. Ayrıca, (2.1) denkleminin (2.2) sınır şartlarını sağlayan her çözümü aşağıdaki eşitsizliği sağlar:
( )
[ )
∞ ∈ ≥ ∈ ≤ + − + + + −∫
∫
( , ) , , 1, 0,11,, . 2 ) , ( 2 2( 2) ) 2 ( ) 1 ( 0 2 0 0 2 t p At p At dx t x u d dx t x u m p p p l xx l t (2.1.1) İspat: ) , (x tu fonksiyonu (2.1) denkleminin (2.2) sınır şartını sağlayan klasik çözümü olsun.
(2.1) denklemini ut ile çarparsak,
[
( )]
+ + 2 =0 − + p+ t tx t x x t xxxx t tt tu kuu u a x u u u bu mu γeşitliğinielde ederiz. Bu eşitlikte,
), ( 2 2 t tt t u t m u mu ∂ ∂ = ( ), 2 2 t tx t u x u u ∂ ∂ =γ γ ∂ ∂ + − ∂ ∂ = 2 2 1 ) ( t xxx tx xx xx xxxx t u t k u u u u x k u ku ,
[
]
[
]
[
]
[
2]
) ( 2 1 ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( x x t xx t xx t x tx x t xx t x t x x t u x a t u u x a x u u x a u u x a u u x a u u x a x u u x a u u x a x u x a u ∂ ∂ + ∂ ∂ − = − + + ∂ ∂ − = − ∂ ∂ − = − olduğundan, eşitliğini[
]
[
]
( ) 0 (2.1.2) 2 ) ( 2 1 ) ( 2 1 ) ( ) ( 2 2 2 2 2 2 = + ∂ ∂ + ∂ ∂ + ∂ ∂ − ∂ ∂ + − ∂ ∂ + ∂ ∂ + p t t x x t xx xx tx xxx t t u b u x u x a t u u x a x u t k u u u u x k u t m γ[
]
( ) 0 2 ) ( 2 1 ) ( 2 ) ( ) ( 2 0 2 0 2 0 2 0 0 2 0 0 0 2 = + + + − + − +∫
∫
∫
∫
+ = = = = = = = = l p t l x x t l x l x x x t l xx l x x xx tx l x x xxx t l t dx u b u dx u x a dt d u u x a dx u k dt d u u k u u k dx u m dt d γeşitliğini elde ederiz. Burada (2.2) sınır şartını kullanırsak, ( ) 0 2 1 2 2 0 2 0 2 2 2 + = + +
∫
∫
+ l p t l x xx t u a x u dx b u dx k u m dt d (2.1.3)eşitliğini elde ederiz. (2.1) denklemini bu kez
u
ile çarparak,[
( )]
+ + =0 − + p t t tx x x xxxx tt kuu ua x u uu buu u muuγ
eşitliğini elde ederiz. Bu eşitlikte, − ∂ ∂ = 2 ) ( t t tt uu u t m muu , − ∂ ∂ = t x t tx uu u u x uu
γ
( )γ
, 2 ) ( xxx x xx xx xxxx uu u u ku x k kuu − + ∂ ∂ = ,[
]
[
]
2 ) ( ) ( ) ( ) ( ) ( x x xx x x x u x a uu x a x uu x a x a x uu u x a u + ∂ ∂ − = − ∂ ∂ − = −olup bunları yukarıdaki yerlerine yazarsak
[
]
(2.1.4) 0 ) ( ) ( ) ( ) ( ) ( 2 2 2 = − − ∂ ∂ + + ∂ ∂ − − + − ∂ ∂ + − ∂ ∂ p t t t x t x x xx xx x xxx t t u buu u u uu x u x a uu x a x ku u u uu x k u uu t m γ[
( )]
( ) ( ) 0 ) ( ) ( m 0 0 0 0 2 0 0 2 0 0 0 2 0 = + − + + − + − + − + ∂ ∂∫
∫
∫
∫
∫
∫
= = = = = = = = l p t t l t x l x x t l x l x x x l xx l x x xx x l x x xxx l t l t dx u uu b dx u u uu dx u x a uu x a dx u k u u k u ku dx u m dx uu t γ γeşitliğini buluruz. Bu eşitlikte (2.2) sınır şartı kullanılarak, (2.1.5) 0 ) , ( ) ( )] , ( [ 0 0 2 2 2 = + − + + −
∫
∫
l p t t t x l x xx t t dx u uu b u u dx u x a u k u m u u m dt d γeşitliği elde edilir.
Keyfi bir δ1>0 sayısı seçerek (2.1.5) eşitliğini
δ
1 ile çarpalım ve sonra da (2.1.3) eşitliği ile toplarsak,(2.1.6) 0 ) , ( ) ( ) , ( ) ( 2 1 2 2 0 2 0 1 0 1 2 1 2 1 2 1 0 1 2 2 2 = + + − + + − + + +
∫
∫
∫
∫
+ l p t l p t t l t x x xx t l t x xx t dx u b dx u uu b u u dx u x a u k u m u u m dx u x a u k u m dt d δ γ δ δ δ δ δeşitliğini elde ederiz. Diğer taraftan,
=
+
+
∫
+
l t x xx tu
a
x
u
dx
m
u
u
k
u
m
t
E
0 1 2 2 2 1(
)
(
,
)
2
1
2
2
)
(
δ
(2.1.7) olsun. (2.1.7) tanımı ve∫
∫
+ − ≥ l p t l p t t u dx b u u dx uu b 0 1 1 0 1 δ δ(2.1.8) 0 ) , ( ) ( ) ( 0 2 0 1 1 0 1 2 1 2 1 2 1 1 ≤ + − − + + −
∫
∫
∫
+ + l p t l p t l t x x xx t dx u b dx u u b u u dx u x a u k u m t E dt d δ γ δ δ δ δeşitsizliği elde edilir.
(2.1.8) eşitsizliğinde
δ
1γ
(ux,ut) terimi için bir kestirim elde etmeye çalışalım. Bunun için önce aşağıdaki eşitsizliğin sağlandığını gösterelim:2 2. 2 2 xx x u l u π ≤ (2.1.9)
( )
uux x =ux +uuxx 2eşitliğini göz önüne alalım. Bu eşitliği
x
değişkenine göre 0 ’danl ’ye integre eder ve (2.2) sınır şartını göz önüne alırsak,
∫
=−∫
l l xx xdx uu dx u 0 0 2eşitliğini buluruz. Yukarıdaki ifadeye Cauchy eşitsizliğini uygularsak,
2 1 0 2 0 2 1 0 2 0 2 ≤ ≤
∫
∫
∫
∫
l xx l l l xx xdx u u dx u dx u dx ueşitsizliğini elde ederiz. Eşitsizliğin sağ tarafındaki ilk terime Wirtinger eşitsizliğini uygulayarak, 2 1 0 2 1 0 2 1 0 2 2
≤
∫
∫
∫
l xx l x xu
dx
u
dx
l
dx
u
π
eşitsizliğini buluruz. Son eşitsizlikte her iki tarafı
2 1 0 2
∫
luxdx terimine bölersek, 2 1 0 2 2 1 0 2
≤
∫
∫
l xx l xu
dx
l
dx
u
π
eşitsizliği elde edilir. Bu eşitsizlikte her iki tarafın karesi alınarak (2.1.9) eşitsizliğini elde ederiz.
Şimdi
δ
1γ
(ux,ut) terimi için bir kestirim elde etmeye çalışalım.γ
’nın bir reel sayı olduğu ve Cauchy eşitsizliği göz önüne alınırsa,t x l t x t x
u
u
u
dx
u
u
u
,
)
(
0 1 1 1γ
≤
δ
γ
∫
≤
δ
γ
δ
eşitsizliğini yazabiliriz. Burada Young eşitsizliğini kullanarak ve sonra da (2.1.9) eşitsizliğini uygulayarak, (2.1.10) 2 2 1 2 2 1 ) , ( 2 2 2 2 2 1 2 2 2 2 1 2 1 xx t x t t x u l u u u u u π γ δ γ δ γ δ + ≤ + ≤
eşitsizliğini buluruz. Elde ettiğimiz bu kestirimi (2.1.8) eşitsizliğinde kullanarak,
(2.1.11) 0 ) ( ) 2 1 (-) 2 ( ) ( 0 2 0 1 1 0 2 1 2 1 2 2 2 2 2 1 1 1 ≤ + − + − + − +
∫
∫
∫
+ l p+ t l p t l x t xx dx u b dx u u b dx u x a u m u l k t E dt d δ δ δ π γ δ δeşitsizliğini yazabiliriz. Şimdi
δ
1 sayısı için, 2 2 0 2 1 2 0 l dγ
π
δ
< < (2.1.12)şartını ekleyelim. Teoremin .ii şartındaki 2 0 0
2
0 = >
−c l d
k
π
eşitsizliğini göz önüne alarak ve(2.1.12) şartı kullanılırsa,
0
)
(
2
2
2 2 0 1 0 1 2 2 2 1 1 2 2 2 2 1 1
≥
−
=
>
−
=
−
=
π
δ
δ
π
γ
δ
δ
π
γ
δ
δ
k
l
k
l
k
d
c
l
L
olduğunu görürüz. (2.1.11) eşitsizliğini göz önüne alalım.
(2.1.13) ) ( ) 2 1 ( ) ( 0 2 0 1 1 0 2 1 2 2 2 2 1 1
∫
∫
∫
+ + − + − − − + + ≤ l p t l p t l x xx t t t dx u b dx u u b dx u x a u L u m u m u m t E dt d δ δ δ(2.1.13) eşitsizliğini kullanarak (2.1) denkleminin sıfır çözümünün kararlı olduğunu şöyle gösterebiliriz. Teoremin .ii şartını, yani; a(x)≥−c0 şartını kullanarak
∫
∫
∫
+
∫
≥
−
l l x xx l l x xxdx
a
x
u
dx
k
u
dx
c
u
dx
u
k
0 0 2 0 2 0 0 2 2)
(
eşitsizliğini yazabiliriz. Burada (2.1.9) eşitsizliğini kullanarak,
(2.1.14) ) ( 0 2 0 2 2 2 0 2 2 2 0 2 2 0 2 0 0 2 2
∫
∫
∫
= − ≥ − ≥ − ≥ + l xx xx xx xx x xx l l x xx dx u d u l c k u l c u k u c u k dx u x a dx u k π π olduğunu görürüz. E(t) fonksiyonunu∫
+
+
=
l x xx tu
a
x
u
dx
k
u
m
t
E
0 2 2 2)
(
2
1
2
2
)
(
şeklinde tanımlayalım. (2.1.14) eşitsizliği ve E(t) fonksiyonu birlikte göz önüne alınırsa,
2 0 2 2 2 ) (t m ut d uxx E ≥ + eşitsizliğini buluruz.
E
t
u
u
xxdx
l t)
(
)
(
2 0 20
=
∫
+
olduğu göz önüne alınırsamin
{
,}
( ) 2 1 ) (t m d0 E0 t E ≥ (2.1.15)eşitsizliğini yazabiliriz. Diğer taraftan a(.)∈C1
[ ]
0,l sınıfından olduğu içina
fonksiyonu[ ]
0
,
l
aralığında birA
1 maksimum değeri vardır. Bu özelliği, E(t) fonksiyonu tanımında kullanarak,∫
∫
+ + ≤ + + ≤ l x xx t l x xx t dx u A u k u m dx u x a u k u m t E 0 2 1 2 2 0 2 2 2 2 2 2 ) ( 2 1 2 2 ) (eşitsizliğini buluruz. Sağ tarafta (2.1.9) eşitsizliğini kullanarak ta
)
(
,
max
2
1
)
(
21 0 2t
E
k
A
l
m
t
E
+
≤
π
eşitsizliğini buluruz. Ayrıca (2.1.3) eşitsizliğinden ve E(t) fonksiyonun tanımından,
( ) 0 0 2 ≤ − =
∫
+ l p t dx u b t E dt d (2.1.16)eşitsizliği sağlanır. Böylece (2.1) denkleminin (2.2) şartını sağlayan sıfır çözümünün kararlı olduğu neticesini çıkarırız.
Şimdi de (2.1) denkleminin (2.2) sınır şartını sağlayan sıfır çözümünün global asimptotik kararlı olduğunu gösterelim.
(2.1.16) eşitsizliğini 0 ’dan t ’ye integre edersek,
∫ ∫
= + − + t l p t dxds u b E t E 0 0 2 0 ) 0 ( ) (eşitliğini buluruz. E(t)≥0 olduğu göz önüne alınırsa yukarıdaki eşitlikten,
∫∫
+ ≤ t l p t b E dxds u 0 0 2 (0) (2.1.17) eşitsizliği yazılır.Şimdi (2.1.13) eşitsizliğine tekrar dönelim ve
∫
−
−
−
l x xx tL
u
a
x
u
dx
u
m
0 2 1 2 2)
(
δ
için bir kestirim elde etmeye çalışalım. Eğer
E
t
a
x
u
xdx
l2 0
0
(
)
=
∫
(
)
terimi negatif değilse,{
1
,
,
}
min
L
k
δ
(2.1.18) ) ( 2 ) ( , , 1 min ) ( 1 0 2 2 2 0 1 2 1 2 2 t E D dx u x a u k u m k L dx u x a u L u m l x xx t l x xx t ≥ + + ≥ + + δ
∫
δ∫
eşitsizliği sağlanır. − = 1 222 0 1 2 2 , 1 minπ
γ
δ
δ
l d k D olsun. Eğera
x
u
xdx
l 2 0)
(
∫
negatif ise, herx
için0 ) (x c a ≥− olduğundan,
∫
∫
≥
+
−
+
+
l x xx t l x xx tL
u
a
x
u
dx
m
u
L
u
c
u
dx
u
m
0 2 0 1 2 2 0 2 1 2 2)
(
δ
δ
eşitsizliğini yazabiliriz. Son eşitsizliğin sağ tarafında (2.1.9) eşitsizliğini kullanarak,
2 2 2 0 1 2 2 0 2 1 2 2
)
(
t xx xx l x xx tu
l
c
u
L
u
m
dx
u
x
a
u
L
u
m
π
δ
δ
≥
+
−
+
+
∫
eşitsizliğini yazabiliriz. Burada L sayısının tanımı ve c0l2 π2 =k −d0 olduğu göz önüne alınırsa uxx 2 teriminin katsayısı,
− = − − − = − 2 2 2 1 0 1 0 1 2 2 2 2 1 1 2 2 0 1 2 ) ( 2
π
γ
δ
δ
δ
π
γ
δ
δ
π
δ
l d d k l k l c Lolur. Bu eşitliği ve
D
2 sayısının tanımı göz önüne alınarak,(
)
+ + ≥ + − ≥ − + ≥ + +∫
∫
l x xx t xx t xx t l x xx t dx u x a u k u m D u k u m l d k u l d u m dx u x a u L u m 0 2 2 2 2 2 2 2 2 2 1 0 1 2 2 2 2 1 0 1 2 0 2 1 2 2 ) ( 2 , 1 min 2 ) ( π γ δ δ π γ δ δ δ( ) 2 ( ). 0 2 2 1 2 2
∫
≥ + + l x xx t Lu a x u dx D E t u m δ (2.1.19)Şimdi D =3 min
{
D1,D2}
dersek, (2.1.18) ve (2.1.19) eşitsizliklerinden,
+
+
∫
≥
l x xx tL
u
a
x
u
dx
D
E
t
u
m
0 3 2 1 2 2)
(
2
)
(
δ
(2.1.20)eşitsizliğini elde ederiz. Eğer (2.1.20) eşitsizliğini (2.1.13) eşitsizliğinde kullanılırsa,
m m B= + + 2 1 1
δ
olmak üzere, () 2 ( ) 0 2 0 1 1 3 2 1∫
∫
+ + − + − ≤ l p t l p t t D E t b uu dx b u dx u B t E dt dδ
(2.1.21)eşitsizliğini buluruz. (2.1.21) eşitsizliğini 0’dan t ’ye integre edersek,
[
]
ds dx u b dxds u u b dxds u B t E E ds s E D t p t l t t l t l p t t∫∫
∫
∫ ∫
∫ ∫
+ + − + + − ≤ 0 2 0 0 0 0 0 0 1 1 2 1 1 3 ) ( ) 0 ( ) ( 2δ
eşitsizliğini elde ederiz. (2.1.16) eşitsizliğine göre E(t) artmayan bir fonksiyon olduğundan
∫
≤
tds
s
E
D
t
tE
D
0 3 3(
)
2
(
)
2
yazılabilir. Bunu yukarıdaki son eşitsizlikte kullanırsak,[
]
(2.1.22) ) ( ) 0 ( ) ( 2 0 2 0 0 0 1 1 0 0 0 0 1 1 2 1 1 3 ds dx u b dxds u u b dxds u u b dxds u B t E E t tE D t p t l t l p t t l t l p t t∫∫
∫∫
∫∫
∫ ∫
+ + + − + + + − ≤ δ δeşitsizliğini buluruz. Bu eşitsizlikte, sağ taraftaki son terim atılarak ifade daha da büyür:
2 ( )
[
(0) ( )]
. 0 0 0 0 1 1 2 1 1 3∫ ∫
∫ ∫
+ + + − ≤ t l t l p t tdxds b u u dxds u B t E E t tE D δ (2.1.23)Şimdi (2.1.23) eşitsizliğinin sağ tarafındaki terimler için birer kestirim elde etmeye çalışacağız. Bunu önce E1(t) fonksiyonunu alttan sınırlandırabileceğimizi göstermekle başlayalım. Poincare-Friedrichs eşitsizliğini ve (2.1.9) eşitsizliğini göz önüne alırsak,
2 4 4 2 2 2 2 xx x u l u l u π π ≤ ≤ (2.1.24) olduğunu görürüz. E1(t)’nin tanımından, a(x)≥−c0 şartını kullanılarak,
) , ( 2 2 2 ) ( 1 0 2 0 2 2 1 t l x xx t u dx m u u c u k u m t E ≥ + −
∫
−δolduğunu elde ederiz. Bu eşitsizlikte, (2.1.9) eşitsizliğini kullanırsak,
) , ( 2 2 ) , ( 2 2 2 ) ( 1 2 0 2 1 2 2 2 0 2 2 1 t xx t t xx xx t u u m u d u m u u m u l c u k u m t E δ δ π − + = − − + ≥
olduğunu görürüz. Bu eşitsizliğin sağ tarafındaki son terime p= q=2 eşlenik sayıları için Young eşitsizliği uygulanırsa,
2 1 2 1 2 0 2 1 2 2 2 2 ) (t m ut d uxx m u m ut E ≥ + −δ −δ
ifadesini buluruz. Şimdi bu eşitsizlikte (2.1.24) eşitsizliğini kullanarak,
(
)
4 2 4 1 0 2 1 1 2 1 1 2 ) (t m ut d ml uxx E − + − ≥ π δ δ (2.1.25)eşitsizliğini elde ederiz. δ 1 sayısı üzerine,
< 4 4 0 2 2 0 2 1 , 2 min ml d l d π γ π δ
şartını eklersek son eşitsizliğin sağ tarafındaki katsayılar negatif değer alamaz. Buna göre
(
)
− = 4 4 1 1 1 min 1 , π δ δ ml m d olmak üzere E1(t)≥d1(
ut 2 + uxx 2)
≥0 (2.1.26) eşitsizliğini elde ederiz. Bu ifade E1(t) fonksiyonun alttan pozitif bir fonksiyonla sınırlı olduğunu gösterir. Buna göre,E1(0)−E1(t)≤E1(0) (2.1.27) ifadesini yazabiliriz.
(2.1.23) eşitsizliğindeki diğer terimleri üstten sınırlandırmak için önce
∫
l q dx u 0 , q>1 integraline bir kestirim bulmaya çalışacağız. u(0,t)=0 olduğundan, Newton-Leibniz formülüne göre,
∫
= − x x s t ds u t u t x u 0 ) , ( ) , 0 ( ) , ( veya∫
= x x s t ds u t x u 0 ) , ( ) , ( eşitliği bulunur. Bu eşitlikten,∫
∫
≤ ≤ l x x x s t ds u x t dx u t x u 0 0 ) , ( 1 ) , ( ) , ( , ∀x∈[ ]
0,lolduğunu kolayca görebiliriz. Burada Cauchy eşitsizliğini kullanırsak,
2 1 0 2( , ) ) , ( ≤
∫
l x x t dx u l t x u , ∀x∈[ ]
0,leşitsizliğini buluruz. Yukarıdaki eşitlik ∀x∈
[ ]
0,l için doğru olduğundan, [ ] 2 1 0 2 , 0 ( , ) ( , ) max ≤∫
∈ l x l x u x t l u x t dx (2.1.28)eşitsizliğini yazabiliriz. (2.1.28) eşitsizliğine göre, q≥1 olmak üzere
[ ] 2 1 0 2 ) 2 1 ( ) 1 ( , 0 1 1 0 ) , ( ) , ( max ) , ( ≤ ≤
∫
∫
+ ∈ l x q l x q q l q dx t x u l t x u l dx t x ueşitsizliği sağlanır. Son eşitsizliği aşağıdaki şekilde yazmak mümkündür:
2 0 2 2 ) 2 ( 0 ) , ( ) , ( q l x q l q dx t x u l dx t x u + ≤
∫
∫
+ . (2.1.29)Şimdi (2.1.29) eşitsizliğini ve Hölder eşitsizliğini kullanarak, (2.1.23) eşitsizliğindeki
∫∫
t l utdxds B 0 0 2 ve∫∫
+ t l p t dxds u u b 0 0 1 1p p p′= +2 ve 2 2 + = ′ p
q eşlenik sayıları için Hölder eşitsizliğini iki defa kullanarak,
) 2 ( 2 0 0 2 ) 2 ( 0 0 ) 2 ( 2 0 2 ) 2 ( 0 0 0 0 2 1 1 1 + + + + + + ≤ ≤ ≡
∫∫
∫∫
∫
∫ ∫
∫∫
p t l p t p p t l p l p t p p t l t l t dxds u dxds B ds dx u dx B dxds u B Ieşitsizliğini elde ederiz. Burada (2.1.17) eşitsizliğini göz önüne alırsak,
(0) , ) 2 ( 2 ) 2 ( 1 + + = p p p b E Bl C I1 ≤C1tp(p+2) (2.1.30)
olduğunu görürüz. (2.1.23) eşitsizliğindeki son terimi ise şöyle kestirebiliriz:
1 2 + + = ′ p p p ve 2 + = ′ p
q eşlenik sayıları için iki defa Hölder eşitsizliğini kullanarak,
) 2 ( 1 0 0 2 ) 2 ( ) 1 ( 0 0 2 1 ) 2 ( 1 0 2 ) 2 ( ) 1 ( 0 0 2 1 0 0 1 1 2 + + + + + + + + + + + ≤ ≤ ≡
∫∫
∫∫
∫
∫ ∫
∫∫
p t l p p p t l p t p l p p p t l p t t l p t dxds u dxds u b ds dx u dx u b dxds u u b I δ δ δeşitsizliğini bulabiliriz. (2.1.29) ve (2.1.9) eşitsizliklerini kullanarak yukarıdaki eşitsizliğin sağ tarafındaki son terim için
2 2 0 2 2 2 4 0 2 4 0 2 2 4 2 + + + + + + ≤ ≤
∫
∫
∫
p l xx p p l p l x p p dx u l l dx u l dx u πeşitsizliğini elde ederiz. Ayrıca (2.1.15) eşitsizliğinden,
∫
≥ l xxdx u t E d 0 2 0 ) ( 2yazabiliriz. E(t) azalan fonksiyon olduğundan E(t)≤E(0), ∀t≥0 ve 0
0 0 2 2 E d dx u l xx ≤
∫
(2.1.31) ) 0 ( ) 2 ( 1 2 ) 2 ( 1 0 2 ) 2 ( 0 2 2 2 ) 4 ( ) 2 ( ) 1 ( 1 2 + + + + + + + ≤ ≤
∫ ∫
p p t l p xx p p p p t C ds dx u l l b E b I π δeşitsizliği bulunur. Burada,
[
]
12 0 ) 2 ( 2 ) 8 3 ( ) 2 ( 2 ) 4 3 ( ) 2 ( 1 1 2 2 ) 0 ( = + + + + + d l E b C p p p p p π δ dir.Şimdi, elde edilen (2.1.27), (2.1.30) ve (2.1.33) kestirimlerini, (2.1.23) eşitsizliğinde kullanırsak aşağıdaki eşitsizliği
) 2 ( 1 2 ) 2 ( 1 1 3 ( ) (0) 2D tE t ≤E +Ctp p+ +C t p+ buluruz. Bu eşitsizliği,
(
)
[
( 1)( 2)]
2 ) 2 ( 2 1 1 1 1 3 (0) 2 ) (t ≤ D − E t− +C t− p+ +C t− p+ p+ Ebiçiminde yazabiliriz. A=
(
2D3)
−1[
E1(0)+C1 +C2]
dersek aşağıdaki ifade her t≥1 için doğrudur:( )
[ )
∞ ∈ ≥ ∈ ≤ + − + + − . , 1 , 1 , , 1 , 0 , ) ( ) 2 ( 2 ) 2 ( ) 1 ( t p At p At t E p p pYukarıdaki eşitsizlikte, (2.1.15) eşitsizliğini göz önüne alırsak,
{
}
( )
[ )
∞ ∈ ≥ ∈ ≤ + − + + − , , 1 , 1 , , 1 , 0 , ) ( , min 2 1 ) 2 ( 2 ) 2 ( ) 1 ( 0 0 t p At p At t E d m p p pelde edilir. Bu eşitsizlik bize (2.1) ve (2.2) probleminin sıfır çözümünün global asimptotik olarak kararlı olduğunu gösterir.
Teorem 2.1.2
IR
p=0, k>0, b>0, γ∈ ve a(x) fonksiyonu, Teoremin 2.1.1 ’in .ii şartını sağlıyorsa,
lineer denkleminin (2.2) sınır şartını sağlayan her çözümü için 2 2 2t, xx t u Ke u + ≤ −δ K,δ2 >0 (2.1.33) eşitsizliği sağlanır [15]. İspat: ), , ( tx
u (2.1.33) denkleminin (2.2) sınır şartını sağlayan bir çözümü olduğunu kabul edelim. p=0 olduğu göz önüne alınırsa (2.1.8) eşitsizliğinden aşağıdaki ifadeyi
(2.1.34) 0 ) , ( ) ( ) ( 0 2 0 1 0 1 2 1 2 1 2 1 1 ≤ + − − + + −
∫
∫
∫
l t l p t l t x x xx t dx u b dx u u b u u dx u x a u k u m t E dt d δ γ δ δ δ δyazabiliriz. δ2∈(0,δ1) keyfi bir sayı olsun. (2.1.7) de E1(t) ’nin tanımı göz önüne alınırsa ve (2.1.34) eşitsizliğine δ2E1(t) ’yi ekleyip çıkarırsak aşağıdaki ifadeyi
(2.1.35) 0 ) , ( ) ( ) , ( ) ( 2 2 2 ) ( ) ( 0 2 0 1 1 0 2 1 2 1 2 1 1 2 0 2 2 2 2 2 2 1 2 1 ≤ + − − + + − − − − − +
∫
∫
∫
∫
l t l t t x l x xx t t l x xx t dx u b dx u u b u u dx u x a u k u m u u m dx u x a u k u m t E t E dt d δ γ δ δ δ δ δ δ δ δ δ δelde ederiz. Burada bazı kestirimlere ihtiyaç duyacağız.
∫
l xdx u x a 0 2 )
( terimi için Teorem 2.1.1 ’in
.
ii şartı ve (2.1.9) eşitsizliği göz önüne alınırsa,
(2.1.36) 2 2 ) ( 2 0 2 2 1 2 2 0 0 2 2 1 0 0 2 2 1
∫
∫
∫
− − ≥ − − ≥ − l xx l x l x dx u l c dx u c dx u x a δ δ π δ δ δ δyukarıdaki eşitsizlik doğrudur. δ1γ(ux,ut) ve
∫
l t dx u u b 0 1δ terimleri için
ε
-Young ve (2.1.9)eşitsizliklerini kullanarak, (2.1.37) 4 2 ) , ( 2 2 2 2 2 1 2 2 2 2 1 2 1 xx t x t t x u b l u b u b u b u u π γ δ γ δ γ δ + ≤ + ≤ (2.1.38) 4 4 2 4 4 2 1 2 2 2 1 2 0 1 xx t t l t u bl u b u b u b dx u u b π δ δ δ + ≤ + ≤
∫
eşitsizlikleri yazılabilir ve benzer şekilde δ2δ1m(u,ut) terimi için de
ε
-Young ve (2.1.24) eşitsizliklerini kullanarak, (2.1.39) 2 2 2 2 ) , ( 2 4 4 1 2 2 1 2 2 1 2 2 1 2 1 2 xx t t t u ml u m u m u m u u m π δ δ δ δ δ δ δ δ δ δ + ≤ + ≤eşitsizliğini yazabiliriz. Elde ettiğimiz (2.1.36), (2.1.37), (2.1.38) ve (2.1.39) eşitsizliklerini (2.1.35) eşitsizliğinde göz önüne alırsak,
(2.1.40) 0 2 2 2 2 2 ) ( ) ( 2 4 4 2 1 2 2 2 2 1 4 4 1 2 0 2 0 1 2 2 1 2 1 1 2 1 ≤ − − − − + − − − + + xx t u bl b l ml d d u m m m b t E t E dt d π δ π γ δ π δ δ δ δ δ δ δ δ δ
eşitsizliği elde edilir. Burada (2.1.25) eşitsizliğini de sağlayan δ1 sayısı ile δ2 ≤δ1 sayısını,
0 ) ( ) ( 2 1 1 t + E t ≤ E dt d δ (2.1.41)
eşitsizliğini doğrulayacak kadar küçük seçebiliriz. Bu eşitsizliğin her iki tarafını eδ2t ile
çarparak,
[
e 2E1(t)]
≤0 dtd δt
eşitsizliğini elde ederiz. Buradan
E (t) E (0)e 2t 1 1 δ − ≤ (2.1.42)
elde edilir. Buna göre (2.1.26) eşitsizliği de göz önüne alınırsa istenen
, 2 2 2 t xx t u Ke u + ≤ −δ K,δ2 >0
eşitsizlik elde edilir.
3. Genel Boru Ve Dönen Mil Denklemlerinin Sıfır Çözümlerinin Asimptotik Davranışı
Bu bölümde amacımız, x ∈(0,L), t∈ IR+ için aşağıdaki problemin sıfır çözümünün global asimptotik kararlı olduğunu göstermektir. Bunun için aşağıdaki problemi alalım.
wtt +awxxxx +bwtxxxx +cwt +α(t)wxx +l(x)β(t)wxx +γ(t)wxt =0 (3.1) w(0,t)=w(L,t)=wx(0,t)=wx(L,t)=0 (3.2) veya
w(0,t)=w(L,t)=wxx(0,t)=wxx(L,t)=0. (3.3)
Burada a, b, c verilen pozitif sayılar ve α(t), β(t), γ(t), l(x) pozitif sınırlı fonksiyonlar olsunlar. vtt +α1vxxxx −a(t)wt −b(t)v−c(t)w+β1vtxxxx+γvt −d(t)w=0 (3.4) wtt +α2wxxxx −a(t)vt −b(t)w+c(t)v+β2wtxxxx +γwt +d(t)v=0 (3.5) + ∈ = = = = = = = = IR t , 0 ) , ( ) , 0 ( ) , ( ) , 0 ( ) , ( ) , 0 ( ) , ( ) , 0 ( t L w t w t L v t v t L w t w t L v t v x x x x (3.6) veya + ∈ = = = = = = = = IR t , 0 ) , ( ) , 0 ( ) , ( ) , 0 ( ) , ( ) , 0 ( ) , ( ) , 0 ( t L w t w t L v t v t L w t w t L v t v xx xx xx xx (3.7)
sistemini göz önüne alalım. Burada verilen α1, α2, β1, β2, γ pozitif sayılar ve ) ( ), ( ), ( ), (t b t c t d t
a pozitif sınırlı fonksiyonlar olsunlar.
(3.2) veya (3.3) sınır şartları altında,
0=Elwxxxx+diElwtxxxx+dewt +MU2wxx +(L−x)MUwxx +2MUwxt +(M +m)wtt (3.8) denklemini ve (3.6) veya (3.7) sınır şartları altında,
El1vxxxx+mvtt −2mΩwt −mΩ2v−mΩw+diEl1vtxxxx+dem(vt −Ωw)=0 (3.9) El2wxxxx +mwtt −2mΩvt −mΩ2w−mΩv+diEl2wtxxxx +dem(wt −Ωv)=0 (3.10) sistemlerini düşünebiliriz [17].
(3.8) denklemi, L uzunluğundaki düzgün bir borunun dinamiğini belirtir. Burada m her birim uzunluktaki kütleyi, El eğilim katsayısını, M her birim uzunluktaki bir sıvı kütlesinin taşınmasını, U(t) pozitif x yönünde zaman değişim hızını, d ve i d sırasıyla iç ve e dış bükülme katsayılarını gösterir. (3.9) ve (3.10) eşitlikleri zaman değişim açısal hızı olan Ω(t) ile dönen bir milin dinamiğini tanımlar. L, m, di ve de (3.8) denkleminde olduğu gibi aynıdır.
1
El ve El eğilim katsayılarıdır. (3.8), (3.2) veya (3.8), (3.3) probleminin sıfır çözümü, 2
(
)
12 0 2 2 +∫
L wxx wt dxnormuna göre belirli sınır şartları altında kararlıdır [17]. Benzer şekilde, (3.9), (3.10), (3.6) veya (3.9), (3.10), (3.7) probleminin sıfır çözümü de
(
)
12 0 2 2 2 2 + + +∫
L vxx wxx vt wt dxnormuna göre belirli sınır şartları altında kararlıdır [17].
Burada amaç; (3.1), (3.2) (veya (3.1), (3.3)) ve (3.4), (3.5), (3.6) (veya (3.4), (3.5), (3.7)) problemlerinin sıfır çözümlerinin global asimptotik kararlı olduğunu yeterli şartlar altında garanti etmektir. Sonuçlarımızı elde etmek için Lyapunov fonksiyon tekniğini kullanacağız. Bu bölümde kullanılan bu metot, lineer olmayan sistemlere de uygundur.
İfade edilen vektörü u ile gösterelim. Bu durumda denge ifadesi de u=0 olur. Sınır şartlarını sağlayan u elemanını içeren uzayı da U ile gösterelim. t=0 için u nun ilk ifadesi
0
u olsun ve t anında ise u(t,u0) olur. Bir özel normda U uzayı üzerinde tanımlansın. Aşağıdaki özellikleri sağlayan U uzayında tanımlanan bir w fonksiyonunu elde ederek Lyapunov’un metodunu genişletebiliriz. w fonksiyonu için, aşağıdaki üç özellik u=0 denge ifadesinin kararlılığı için yeterli bir koşulu gösterir , bu koşullar;
i. ∀u ∈U için w(u)≤c1 u 2 şartını sağlayan bir c1>0 vardır, ii. ∀u ∈U için w(u)≥c2 u 2 şartını sağlayan bir c2 >0 vardır, iii. w
(
u(t,u0))
≤0,dt d
dır [18-20].
Eğer sıfır çözümü kararlı ve tüm çözümler t→∞ iken sıfıra yaklaşıyorsa, sıfır çözümü global asimptotik kararlı olarak adlandırılır.
Burada, 2 1 0 2( ) =
∫
L dx x w w ve =∫
L dx x v x w v w 0 ) ( ) ( ) , ( anlamındadır. 2 2 π λ=L olmak üzere,∫
≤∫
L x L dx x w dx x w 0 2 0 2( ) λ ( ) (3.11)Wirtinger eşitsizliğini göz önüne alalım. Bu durumda aşağıdaki eşitsizliğin
∫
≤∫
L xx L x x dx w x dx w 0 2 0 2( ) λ ( ) (3.12)sağlandığını gösterelim.
(
wwx)
x =wx2 +wwxx eşitliğini göz önüne alalım. Bu eşitliğix
değişkenine göre 0 ’dan L ’ye integre eder ve (3.2) sınır şartını göz önüne alarak,
(
)
(
)
dx ww dx w dx ww dx w t w t w t L w t L w dx ww w dx ww L xx L x L xx L x x x L xx x L x x∫
∫
∫
∫
∫
∫
= − ⇒ + = − ⇒ + = 0 0 2 0 0 2 0 2 0 ) , 0 ( ) , 0 ( ) , ( ) , (eşitliği elde edilir. Son eşitliğin sağ tarafına Cauchy eşitsizliğini uygulayıp gerekli düzenlemeler yapılarak ≤
∫
∫
∫
L xx L L xdx w dx w dx w 0 2 0 2 2 0 2eşitsizliği elde edilir. Burada (3.12) Wirtinger eşitsizliği göz önüne alınarak
≤
∫
∫
∫
L xx L x L xdx w dx w dx w 0 2 0 2 2 0 2 λeşitsizliği elde edilir. Eşitsizliğin her iki tarafını
∫
L xdx
w
0
2 terimine bölersek, (3.12) eşitsizliği elde
edilir.
∫
≤∫
L xx L dx x w dx x w 0 2 2 0 2 ) ( ) ( λ (3.13) eşitsizliği elde edilir.3.1. Global Asimptotik Kararlılık:
Teorem 3.1.1
Aşağıdaki şartların sağlandığını kabul edelim:
i. a, b, c verilen pozitif sayılardır.
ii. α(.), β(.) aşağıdaki şartları sağlayan C1
[
0,∞)
sınıfından pozitif fonksiyonlardır;0 2 1 2 ) ( ) ( ) ( ) (t +β t L+β t L λ +γ t ≤a α , (3.1.1) − ≤ + ′ + ′(t) (t)L 2(t)L1 n0 a a0 λ β β α , (3.1.2) burada a , 0 λ a a <0 ve − + − − = 3 2 , , 1 min 1 2 2 0 0 L c b b a a n λ λ λ
yukarıdaki şartı sağlayan keyfi bir pozitif sayı olsun. iii. γ(.), C1
[
0,∞)
sınıfından pozitif bir fonksiyondur. iv. l(.)∈C1[ ]
0,L ve 0≤l(x)≤L, (3.1.3) l′(x) ≤L1 , ∀x∈[ ]
0,L (3.1.4) olup burada L , 1 + < b c L1 2 2 λ (3.1.5)eşitsizliği sağlayan pozitif bir sayıdır. (3.1), (3.2) veya (3.1), (3.3) probleminin sıfır çözümü,
(
)
12 0 2 2 +∫
L wxx wt dxnormuna göre global asimptotik olarak kararlıdır. Bununla birlikte (3.2) veya (3.3) sınır şartlarını sağlayan, (3.1) denkleminin her çözümü için
wxx 2 + wt 2≤K1 e−δt; K1,δ >0, (3.1.6) eşitsizliği sağlanır. İspat 3.1.1 ) , ( tx
w , (3.2) veya (3.3) sınır şartlarını sağlayan (3.1) denkleminin bir çözümü olduğunu kabul edelim ve n ise daha sonra tanımlanacak bir parametre olsun. (3.1) denklemini
t w ile çarparsak, 0 ) ( ) ( ) ( ) ( 2 + + + = + + + t xxxx t txxxx t t xx t xx t xt tt tw aww bw w cw t ww l x t ww t w w w α β γ
eşitliğini elde ederiz. Bu eşitlikte,
( )
2 2 1 t tt t w t w w ∂ ∂ = ,( )
2 2 ) ( ) ( t xt wt x t w w t ∂ ∂ =γ γ ,(
)
( )
2 2 xx xx tx xxx t xxxx t w t a w w w w x a w aw ∂ ∂ + − ∂ ∂ = ,(
)
(
)
2 txx txx tx txxx t txxxx t w w bw x b w w x b w bw + ∂ ∂ − ∂ ∂ = ,(
2)
(
)
( ) 2 2 1 ) ( ) ( 2 1 ) ( t xx x t wtwx t wx x w t t w w t α α α α + ′ ∂ ∂ + ∂ ∂ − = ,(
)
(
t x)
t x x x xx t w w x l t w w x l x t w x l t w x l t t w w t x l ) ( ) ( ) ( ) ( ) ( ) ( 2 1 ) ( ) ( 2 1 ) ( ) ( 2 2 ′ − ∂ ∂ + ′ + ∂ ∂ − = β β β β βolduğundan, yukarıdaki eşitliği
( )
2(
)
( )
2 2 2 1 xx xx tx xxx t t w t a w w w w x a w t ∂ ∂ + − ∂ ∂ + ∂ ∂(
t txxx)
(
tx txx)
txx t(
x)
(
t wtwx)
x w t t cw bw w w x b w w x b ( ) ( ) 2 1 2 2 2 α α ∂ ∂ + ∂ ∂ − + + ∂ ∂ − ∂ ∂ +(
)
( ) ( ) ( )(
( ))
( ) ( ) 0 2 1 ) ( ) ( 2 1 ) ( 2 1 2 2 2 − ′ = ∂ ∂ + ′ + ∂ ∂ − ′ + x x x l x wtwx t l x wtwx x t w x l t w x l t t w t β β β β αbiçiminde yazabiliriz. Bu eşitliği, x değişkenine göre integre edersek ve sonra da (3.2) veya (3.3) sınır şartlarını göz önüne alarak,
− − +
∫
L x x xx t w t w t l x w dx a w dt d 0 2 2 2 2 ) ( ) ( 2 1 ) ( 2 1 2 2 1 β α ( ) ( ) ( ) ( ) 0 2 1 ) ( 2 1 0 0 2 2 2 2 = ′ − ′ + ′ + + +∫
∫
L x t L x x t txx c w t w t l x w dx t l x ww dx w b α β β (3.1.7)eşitliği elde edilir.
(3.1) denklemini bu kez nw ile çarparak,
0 ) ( ) ( ) ( ) ( + + = + + + + xxxx txxxx t xx xx xt tt anww bnww cnww n t ww nl x t ww n t ww nww α β γ
eşitliğini elde ederiz. Bu eşitlikte,
(
)
2 t t tt ww nw t n nww − ∂ ∂ = ,( )
2 2 t w cn cnwwt ∂ ∂ = ,(
)
(
)
2 xx xx x xxx xxxx w w anw x an ww x an anww + ∂ ∂ − ∂ ∂ = ,( )
xx(
txxx xxt x)
txxxx w w w w x bn w t bn bnww − ∂ ∂ + ∂ ∂ = 2 2 ,(
)
( ) 2 ) ( ) ( xx wwx n t wx x t n ww t nα α − α ∂ ∂ = ,(
x)
x x xx l x ww nl x t w nl x t ww x t n ww t x nl( )β( ) β( ) ( ) − ( )β( ) 2 − ′( )β( ) ∂ ∂ = ,(
t)
x t xt ww n t w w x t n ww t nγ( ) γ( ) − γ( ) ∂ ∂ =olup bunları yukarıda elde ettiğimiz denklemde yerlerine yazarsak
(
)
2(
)
(
)
2 xx xx x xxx t t w w anw x an ww x an nw ww t n + ∂ ∂ − ∂ ∂ + − ∂ ∂( )
2(
)
( )
2 ( )(
)
( ) 2 2 2 xx txxx xxt x t w n t x wwx n t wx cn w w w w x bn w t bn α α − ∂ ∂ + ∂ ∂ + − ∂ ∂ + ∂ ∂ +(
( ))
( ) ( ) ( ) ( ) ( )(
)
( ) 0 ) ( 2 − = ∂ ∂ + ′ − − ∂ ∂ + x x x wwt n t wxwt x t n ww t x l n w t x nl ww x l x t nβ β β γ γeşitliğini elde ederiz. Bu eşitliği, x değişkenine göre integre edersek ve sonra da (3.2) veya (3.3) sınır şartlarını göz önüne alarak,
2 2 2 2 2 2 ) , (w wt bn wxx cn w nwt an wxx n dt d + − + + ( ) ( ) ( ) ( ) ( ) ( )( , ) 0 0 0 2 2 = − ′ − − −
∫
∫
x t L x L x x n t l x w dx n t l x ww dx n t w w w t nα β β γ (3.1.8)eşitliğini elde ederiz.
(3.1.7) ve (3.1.8) eşitliklerini taraf tarafa toplarsak,
− − +
∫
L x x xx t w t w t l x w dx a w dt d 0 2 2 2 2 ) ( ) ( 2 1 ) ( 2 1 2 2 1 β α 2 2 2 2 2 ) ( 2 1 2 2 ) , (w wt bn wxx cn w b wtxx c wt t wx n + + + α′ + + +∫
∫
− ′ − + ′ + L xx t x t L xdx t l x ww dx n w anw w x l t 0 2 2 0 2 ( ) ( ) ) ( ) ( 2 1 β β ( ) ( ) ( ) ( ) ( ) ( )( , ) 0 0 0 2 2 − − ′ − = −∫
∫
x t L x L x x n t l x w dx n t l x ww dx n t w w w t nα β β γ (3.1.9)eşitliğini elde ederiz.
2 2 0 2 2 2 2 2 2 ) , ( ) ( ) ( 2 1 ) ( 2 1 2 2 1 ) , ( w cn w bn w w n dx w x l t w t w a w w w E xx t L x x xx t t n + + + − − + = α β
∫
(3.1.10) ve ) , )( ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( 2 1 ) ( 2 1 ) , ( 0 0 2 2 0 2 2 0 2 2 2 2 t x L x L x x L xx t x t L x x t txx t n w w t n dx ww x l t n dx w x l t n w t n w an w n dx w w x l t dx w x l t w t w c w b w w H γ β β α β β α − ′ − − − + − ′ − ′ + ′ + + =∫
∫
∫
∫
(3.1.11)olsun. (3.1.9) eşitliğinde, (3.1.10) ve (3.1.11) eşitliklerini göz önüne alarak, 0 ) , ( ) , ( t + n t = n w w H w w E dt d (3.1.12) eşitliğini elde ederiz.
) , (w wt
t t t w w n w w n w w n ≤ ≤ ) , ( ) , (
eşitsizliğini yazabiliriz. Burada Young eşitsizliğini kullanarak ve sonra da (3.13) eşitsizliğini uygulayarak, 2 2 2 2 2 2 2 2 1 2 1 ) , ( xx t t t w n w n w w n w w n λ + ≤ + ≤ (3.1.13) eşitsizliğini buluruz. (3.1.10) da (3.11), (3.12) ve (3.1.13) eşitsizliklerini kullanarak,
[
2]
2 (1 ) 2 2 1 ) ) ( ) ( ( ) ( 2 1 ) , ( t xx t n w w a n b t t L w n w E ≥ − λ − −λ α +β + − (3.1.14)eşitsizliğini elde ederiz.
Eğer λ2 − b<0 ise, n<1 olmak üzere uygun bir k sabiti için (3.1.1) eşitsizliğinden, 0
En(w,wt)≥k0
(
wxx 2 + wt 2)
(3.1.15) eşitsizliğini elde ederiz.Eğer λ2 − b>0 ise, − − > b a a n 2 0 , 1 min λ λ (3.1.16) olmak üzere uygun bir k sabiti için (14) eşitsizliğinden, 1
En(w,wt)≥k1
(
wxx 2 + wt 2)
(3.1.17) eşitsizliğini elde ederiz.(3.1.8) de, (3.1.13) eşitsizliği kullanılarak,
(
2 2 0)
2 2 2 1 2 1 ) , ( t xx t n w w n nb nc a a w w E ≤ λ + + λ + + λ + (3.1.18)eşitsizliği elde edilir.
{
n nb nc a a n}
k2 =max λ2 + + λ2 + + 0λ,1+
olsun. Bu durumda, (3.1.18) eşitsizliğinden,
En(w,wt)≤k2
(
wxx 2 + wt 2)
(3.1.19) eşitsizliği elde edilir.Diğer taraftan γ(t)(wx,wt) terimi için bir kestirim elde etmek üzere Cauchy eşitsizliğini kullanarak, t x t t x w w t w w t w w t ) ( ) , ( ) ( ) , )( ( γ γ γ ≤ ≤
yukarıdaki eşitsizliği elde ederiz. Burada Young eşitsizliğini kullanarak,
2 2 2 2 1 2 ) ( ) , ( ) (t wx wt ≤γ t wx + wt γ (3.1.20)
eşitsizliği elde edilir.
Şimdi
∫
′ L xdx ww x l t 0 ) ( ) (β terimi için bir kestirim elde edelim. Bunun için sırasıyla
Cauchy eşitsizliğini ve (3.1.4) eşitsizliğini kullanarak,
w w L t w w x l t w w x l t w w x l t dx ww x l t x x x x L x 1 0 ) ( ) ( ) ( ) ( ) ( ) , ) ( ( ) ( ) ( ) ( β β β β β ≤ ′ ≤ ′ ≤ ′ ≤ ′
∫
eşitsizliğini elde ederiz. Son eşitsizlikte, (3.11) eşitsizliği göz önüne alınarak,
1 2 0 ) ( ) ( ) ( x L xdx t L w ww x l t β λ β
∫
′ ≤ (3.1.21)eşitsizliği elde edilir.
(3.1.11) eşitliğindeki
∫
′ L t xwdx w x l t 0 ) ( ) (β terimine sırasıyla Cauchy eşitsizliğini ve
(3.1.4) eşitsizliğini kullanarak, t x t x t x x L t x w w L t w w x l t w w x l t w w x l t dx w w x l t 1 0 ) ( ) ( ) ( ) ( ) ( ) , ) ( ( ) ( ) ( ) ( β β β β β ≤ ′ ≤ ′ ≤ ′ ≤ ′
∫
2 1 2 1 2 0 2 2 ) ( ) ( ) ( x t L t x w L w L t dx w w x l t
∫
′ ≤ β + β (3.1.22)eşitsizliği elde edilir.
(3.1.11) de (3.12), (3.13), (3.1.3), (3.1.20), (3.1.21) ve (3.1.22) ifadeleri göz önüne alınarak,
(
)
2 1 2 2 1 2 2 1 2 3 2 ) ( ) ( ) ( 2 1 2 ) ( ) ( ) ( ) ( ) , ( t x t n w n L c b v L t L t t t L t L t t a n w w H + − − + + ′ + ′ − − − − − ≥ λ β β α γ λ β β α λ (3.1.23)eşitsizliği elde edilir. Eğer (3.1.1) ve (3.1.3) şartlarını (3.1.23) eşitsizliğinde kullanarak,
− + − − = 3 2 , , 1 min 1 2 2 0 0 L c b b a a n λ λ λ olmak üzere 0 0 2 n n n ≤ ≤ için Hn(w,wt)≥0 ve En(w,wt)≤0 dt d
elde edilir. Böylece, )
,
( t
n w w
E ’nin (3.1), (3.2) veya (3.1), (3.3) problemi için Lyapunov fonksiyonu olduğunu elde
ederiz. Bu yüzden,(3.1), (3.2) veya (3.1), (3.3) probleminin sıfır çözümü
(
)
12 0 2 2 +∫
L wxx wt dxnormuna göre kararlıdır. δ >0 olmak üzere (3.1.12) den,
) , )( ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( 2 1 ) ( 2 1 2 2 ) , ( ) ( ) ( 2 ) ( 2 2 2 ) , ( ) , ( 0 0 2 2 2 2 0 0 2 2 2 2 2 2 0 2 2 2 2 t x L x L x x xx t L x t L x x t txx xx t L x x xx t t n t n w w t n dx ww x l t n dx w x l t n w t n w an w n dx w w x l t dx w x l t w t w c w b w cn w bn w w n dx w x l t w t w a w w w E w w E dt d γ β β α β β α δ δ δ β δ α δ δ δ δ + ′ + + + − + ′ − ′ − ′ − − − + + + − − + = +
∫
∫
∫
∫
∫
(3.1.24)eşitliği elde edilir. (3.1.24) de (3.11), (3.12), (3.13) ve (3.1.1) ifadelerini kullanarak,
(
)
2 2 1 2 1 2 2 1 2 2 0 2 2 3 2 2 n ) ( ) ( ) ( 2 na 2 ) ( ) ( ) ( ) ( 2 2 2 2 ) ( ) , ( ) , ( t xx t n t n w b c L n w L t L t t t L t L t t n nc b n n a a w w E w w E dt d + + + − − + + ′ + ′ + + + + + + + + + = + λ δ δ β β α λ γ λ β β α λ λ δ δ λ δ λ δ δ (3.1.25)eşitliğini elde ederiz.
0 0
2 n n
n ≤
≤ olmak üzere (3.1.25) de (3.1.1) ve (3.1.2) eşitsizlikleri göz önüne alınarak,
(
)
2 2 1 2 0 0 2 2 0 2 2 3 2 2 n 2 2 2 2 2 ) ( ) , ( ) , ( t xx t n t n w b c L n w a a n n nc b n n a a w w E w w E dt d + + + − − + − − + + + + ≤ + λ δ δ λ λ δ δ λ δ λ δ δ (3.1.26)eşitsizliği elde edilir.
(3.1.26) eşitsizliğinde δ >0 sayısını 0 ) , ( ) , ( t + n t ≤ n w w E w w E dt d δ (3.1.27) eşitsizliğini sağlayacak şekilde seçebiliriz.
Bu eşitsizliğin her iki tarafı eδt ile çarpılırsa,
[
]
(
)
(
( ,0), ( ,0))
(3.1.28) ) , ( 0 ) 0 , ( ), 0 , ( ) , ( 0 ) , ( 0 ) , ( ) , ( t x n t n x n t n t t n t t n t t n t e x w x w E w w E x w x w E w w E e w w E e dt d w w E e w w E dt d e δ δ δ δ δ δ ≤ ⇒ ≤ − ⇒ ≤ ⇒ ≤ +eşitsizliği elde edilir. (3.1.28) eşitsizliğinde (3.1.15) veya (3.1.17) eşitsizlikleri göz önüne alınarak,
t t
xx w K e
w 2 + 2 ≤ 1 −δ