Posts by Călin Vamoş

Abstract

 

Authors

A. Georgescu

C. Vamoş

N. Suciu

Keywords

References

PDF

soon

About this paper

Cite this paper as:
Journal

St. Cerc. Mec. Apl.

Publisher Name
DOI

Not available yet.

Print ISSN

Not available yet.

Online ISSN

Not available yet.

Google Scholar Profile

Paper in HTML form

1990 Suciu-Georgescu-Vamos - Scenarii de turbulenta in cdrul haosului determinist

SCENARII DE TURBULENTĂ ÎN CADRUL HAOSULUI DETERMINIST *)

de ADELINA GEORGESCUInstitutul de MatematicăBucureştiC. VAMOŞCentrul de Prevedere a Vremii ClujN. SUCIUInstitutul de Meteorologie și Hidrologie București

In lucrare se analizează cîteva noi direcţii în matematica aplicată (haos determinist, turbulență și fractali).

FLUGTUATIII EXTERNE

Intre cele două regimuri principale (laminar și turbulent) de mișcare a fluidelor există un regim de trecere, numit tranziție, care prezintă atît proprietăți laminare, cît și turbulente. Curgerile laminare sînt descrise macroscopic cu ajutorul unor cîmpuri fundamentale (viteza, presiunea, densitatea, temperatura), care constituie soluția clasică a modelului Navier-Stokes-Fourier (N-S-F) diferențial sau integral. Să considerăm pentru simplitate cazul materialelor ne micropolare. Atunci acest model constă din patru ecuații de bilanț (pentru următoarele mărimi fizice: masă, impuls, energie și entropie) și condiții la limită și inițiale, în clase de, funcții Riemann integrabile. Forma diferențială a unei astfel de ecuații de bilanț este
d g d t = f + p + s d g d t = f + p + s (dg)/((d)t)=f+p+s\frac{\mathrm{d} g}{\mathrm{~d} t}=f+p+sdg dt=f+p+s
unde : g , f , p g , f , p g,f,pg, f, pg,f,p și s s sss sint respectiv densități ale mărimii fizice considerate ale fluxului, producției și aportului de astfel de mărime. Specificarea materialului se face postulînd pe bază experimentală sau deducînd din teoria constitutivă legături între f , p f , p f,pf, pf,p și s s sss pe de o parte şi aceşti g g ggg pe de alta, în legături (ecuații constitutive) apărînd presiunea și temperatura. Anumiți coeficienţi care apar datorită liniarității ecuațiilor constitutive se deter-
mină experimental macroscopic sau pe baze statistice. Deja de aici sînt posibile două interpretări ale mărimilor caracteristice fundamentale :
a) ele sînt medii statistice ale mărimilor aleatoare ce caracterizează mişearea la nivel microscopic. In acest caz mărimile nou apărute în ecuațiile constructive au interpretări pur macroscopice, coeficienții se determină experimental, iar caracteristicile fundamentale sînt ne aleatoare. Această interpretare conduce la existența unui prag (limita continuumului) de valabilitate a modelării ;
b) caracteristicile fundamentale sînt mărimi aleatoare, dar în modelare intră doar media lor statistică și anume : în membrul stîng al lụi (1) g g ggg reprezintă media mărimii g ~ g ~ tilde(g)\tilde{g}g~, iar f , p f , p f,pf, pf,p și s s sss sint legături macroscopice care modelează efectul fluctuațiilor g ~ g g ~ g tilde(g)-g\tilde{g}-gg~g microscopice la scară macroscopică. Această interpretare nu este obișnuită în cazul laminar.
În mişearea turbulentă, cîmpurile g g ggg nụ mai sînt fizic relevante, iar aparatura de măsură folosită în regim laminar, aici nu mai este adecvată. Turbulența este un regim de curgere principial distinct de cel laminar : în cazul turbulent modelul N-S-F își păstrează valabilitatea, dar cîmpurile. (instantanee) furnizate de el nu mai sînt folositoare, datorită variaţiilor deosebit de rapide și complicate, practic aleatoare. O variație asemănătoare are loc pentru functia f : ( 0 , ) R , f ( x ) = sin x 1 f : ( 0 , ) R , f ( x ) = sin x 1 f:(0,oo)rarrR,f(x)=sin x^(-1)f:(0, \infty) \rightarrow \mathbb{R}, f(x)=\sin x^{-1}f:(0,)R,f(x)=sinx1 care, pentru x x xxx mari, este lent variabilă în vreme ce variația sa devine din ce în ce mai rapidă, cînd x 0 x 0 x rarr0x \rightarrow 0x0. Cu toate acestea f f fff este infinit diferențiabilă dar stabilirea valorii sale într-un punct x 0 x 0 x_(0)x_{0}x0 apropiat de 0 este dificilă, cerînd o precizie din ce în ce mai mare. Foarte aproape de 0 este practic imposibilă o astfel de determinare, o mică abatere în specificarea lui x 0 x 0 x_(0)x_{0}x0 conducînd pentru f ( x 0 ) f x 0 f(x_(0))f\left(x_{0}\right)f(x0) la orice valoare cuprinsă între -1 și +1 (fig. 1). Literatura
Fig. 1
oferă multe exemple de comportamente complicate ale unor funcții ce descriu dinamici deterministe. Deci, caracterul complicat al traiectoriilor în mișcare turbulentă și rapida variație a cîmpurilor menționate, nu implică neapărat cerința ca aceste cîmpuri să fie aleatoare şi deci ca modelul N -S -F să fie considerat în clase de funcții aleatoare. Totuși, practic, acest model nu mai este adecvat pentru descrierea la scară macroscopică a turbulenței, tot aşa cum în regimul laminar nu era relevantă o descriere microscopică discretă. Atunci, avem două posibilități analoage lui a) şi b) : a') mărimile caracteristice mişeării turbulente sînt medii statistice ale celor laminare, iar f , p f , p f,pf, pf,p și s s sss sint legate de acestea prin postulări de modele de turbulență și fără să legăm pe f , p f , p f,pf, pf,p și s s sss de vreo caracteristică a mișcării laminare. In acest caz, mediile statistice din turbulență apar ca duble medieri ale mărimilor microscopice ; b') caracteristicile turbulente sînt
mărimile aleatoare (sau mai bine zis puternic variabile şi complete) întilnite în modelarea laminară in timp ce f , p f , p f,pf, pf,p și s s sss sînt rezultatul la scarai turbulenței a fluctuațiilor din caracteristicile laminare. Aceste interpretări țineau seama de modelul N-S-F valabil în cazul laminar. Există însă și o altă posibilitate și anume de a nu mai lega caracteristicile turbulente de cele laminare și a postula de la inceput că mărimile caracteristice tulbulenței sînt aleatoare, dar că acelea fizic relevante sînt mediile lor statistice legate de o scară molară, f , p f , p f,pf, pf,p şi s s sss incluzind efectul macroscopic al fluctuațiilor la aceea scară [1], [2]. Atunci, in loc de modelul N-S-F determinist, avem unul stocastic, in clase de astfel de functii aleatoare. In teoriile clasice ale turbulenței se postulează tacit că în regim turbulent este valabil un model stochastic in clase de funcții aleatoare, fără să se specifice legătura între mediile statistice ale acestor funcţii și mărimile din modelul N-S-F determinist (care la numere Reynold mici era laminar) și fără să se lege scara continuumului turbulent, continuu construit por-nindu-se de la un sistem discret de grupări de molecule (moli) care pot să joace rolul moleculelor din cazul laminar. Această nespecificare face imposibilă o comparare și o legătură între mişearea laminauă și cea turbulentă, deoarece cele două regimuri sînt caracterizate de modele distincte și pe de altă parte provoacă serioase dificultăți în interpretarea fizică a rezultatelor teoretice. Se pare că modelul lungimii de amestec al lui Prandtl, bazat pe experimentele lui Nikuradze, ia în considerație scara molară drept scara microscopică a turbulenței, ceea ce explică concordanța sa deosebită cu rezultatele măsurătorilor. Acest punct de vedere este actualmente dezvoltat de şcoala sovietică condusă de V. V. Struminskii [2].
Caracterul aleator al solutiei modelului N-S-F stocastic provine din caracterul aleator al datelor (condiții inițiale, condiții la limită, vîscozitate turbulentă, forțe exterioare aleatoare). De aceea in acest model fluctuatiile au natură exterioară (zgomot) [11]), sînt introduse in sistemul fluid și nu se dezvoltă din perturbaţii. In consecință, modelul N-S-F este adecvat descrierii turbulenței dezvoltate şi nu tranziției.
Modelele matematice efectiv folosite in investigarea turbulenţei dezvoltate se obţin simplificind modelul N-S-F stocastic intr-un mod diferit decit in cazul laminar, deoarece aici viscozitatea, nestationaritatea şi amestecul sînt condiții necesare de turbulenţa. Aceste simplificări se fac diferenţiat, tinîndu-se seama de proprietățile mediilor și corelațiilor și a legăturii dintre ele. Astfel, in teoriile statistice naive de tip Reynolds se postulează relaţia de închidere (modele de turbulență) de tip algebric, diferențial, integral sau mixt intre medii temporale și corelatii. În teoriile statistice propriu-zise se folosese medii probabilistice, dar si aici se introduc ipoteze simplificatoare : asupra mediilor și anume izotropia (in teoria lui Kolmogorov), sau se postulează relaţii între spectre și corelaţii (ca în teoria lui Heisenberg) sau se prescriu ordinele de mărime ale corelațiilor (în teoria lui Kraichnan) [3]. Teoria statistică a lui Hopf pare cea mai naturală și generală, singura care să permită o trecere spre interpretări laminare. Complexitatea acestei teorii face ca interesul actual pentru ea să fie pur teoretic. Pe baza ei s-a dezvoltat teoria solutiilor statistice bidimensionale ale lui Foias, de asemenea, deocamdata, numai de
interes teoretic, datorită gradului mare de abstracţie al aparatului matematic, folosit, dar care deja a zdruncinat credinţa de un secol că orice mişcare tubulentă este tridimensională [4].

HAOSUL DETERUINIST

În turbulenţă, simplificările nu privese mărimile instantanee, ci mediile, corelaţiile sau spectrele lor ; rezultă că în interacțiunile model fizic - model matematic - experiment, modelul fizic este cel mai important, deşi el este mereu corectat de experiment şi teorie (fig. 2). A se compara cu situația mai simplă din cazul laminar (fig. 3). Să remarcăm con-
Fig. 2
Fig. 3
diționarea experimentului de către ideea fizică, obligatorie, de altfel, în toată fizica [5], [6], [7]; de aceea, in laboratoarele de prestigiu, experimentatorii sînt aleşi dintre cei mai buni cunoscători ai teoriilor fizice. Pe de altă parte, se constată că dacă în cazul general matematicienii depind de ideiile fizice și de datele experimentale, în turbulență ei furnizează in plus, sugestii fizicienilor și experimentatorilor, apărînd astfel, concepte fizice, gîndite inițial matematic. In acest fel, a luat naştere matematica fizică, care apare în toate situaţiile in care ideile fizice sînt complexe, ne intuitive, ci sugerate de modelul matematic, deci ea apare mult și in turbulență.
Problema principială, majoră a mecanicii fluidelor este apariția comportamentului turbulent aleator, haotic dintr-un comportament regulat, laminar; cum ajung perturbațiile interne, caracterizate de mărimi nealeatoare, să se transforme în fluctuaţii (aleatoare) fără aport de stocasticitate din exterior. Cunoaşterea diferitelor tipuri particulare de turbulențe date de teoriile clasice menţionate este importantă practic dar stăpînirea mecanismului originii turbulenței ar fi covîrşitoare, deoarece ea ar permite adevăratul control al turbulenței, adică folosirea în profunzime, dirijata, a caracteristicilor turbulenţei şi nu în ultimul rînd a energiei sale mari. In plus, observarea structurilor coerente laminare dezvoltate în mişcări turbulente sporesc impasul în care se găsesc fundamentele teoriei turbulenței. Dificultaţile semnalate se întîlnesc in forme specifice In studiul celor mai diferite sisteme fizice şi ele se formulează astfel : de ce evoluţiile ordonate devin haotice şi de ce pe un fond de dezordine poate apare ordinea? Un răspuns la aceste întrebări îl oferă teoriile bazate pe conceptul de haos determinist : evoluţiile complet deterministe conțin în ele însele germenii dezordinii, fără aport de zgomot din exterior. Din punct de vedere matematic haosul determinist este o dinamică parti-
culară în spațiul fazelor (ataşat evoluției sistemului fizic, evoluție ce depinde de anumiţi parametri) in care portretul de fază are proprietăți probabiliste. De aceea, haosul deterministe este numit și haos de fază sau haos din spațiul fazelor. Interpretarea și caracterizarea sa matematică aparține matematicii fizice în care fenomenul fizic este interpretat pornin-du-se de la concepte matematice ataşate dinamicii din spațiul fazelor. Astfel, evoluția unui sistem fizic se numeşte haotică, dacă sistemul dinamic ataşat are un comportament haotic. Dinamica haotică de fază are diferite definiții. De exemplu, după Yorke, dinamica este haotică dacă mulțimea atractivă din spațiul fazelor conține o infinitate de orbite periodice, o infinitate numărabilă de orbite neperiodice și majoritatea orbitelor sînt hiperbolice, iar după Mandelbrot, dacă sistemul dinamic are un atractor fractal. Indiferent de definiția matematică, haosul determinist implică o dezordine macroscopică, constatată experimental sau observată direct. Trecerea de la dinamică regulată la una haotică, rămînind în cadrul determinist, are loc conform unor "scenarii", „drumuri spre turbulență" tipice. Prin turbulență, aici se înțelege o dinamică analoagă turbulenței dezvoltate a fluidelor și caracterizată de mărimi aleatoare. Haosul determinist este capătul "drumului", pe cînd generarea sa are loc pe parcursul unei întregi zone de variație a parametrilor fizici, corespunzătoare tranziţiei laminar - turbulent.

EXEMPLE MODEL

In marea lor majoritate, fenomenele lumii vii sau nevii sînt modelate de probleme pentru ecuaţii cu derivate partiale al căror comportament este departe de a fi înțeles. In anumite situații (e.g. cazul disipativ) unui astfel de fenomen i se ataşează o problemă Cauchy pentru un sistem de ecuații diferențiale ordinare, numită exemplu model, care păstrează caracteristicile esențiale de fază ale fenomenului fizic. Ideea reducerii studiului calitativ al unei ecuații cu derivate partiale la acela al unei ecuatii diferentiale ordinare este una dintre cele mai profunde din ultimii ani. Ea a apărut în urma prelucrării unei cantități imense de rezultate obținute la calculator, cu care ocazie s-a observat că mulțimea soluţiilor anumitor ecuaţii cu derivate parțiale se „aplatisează" ținînzînd, pe măsura depărtării de momentul inițial ( t = 0 ) ( t = 0 ) (t=0)(t=0)(t=0), către o mulţime de dimensiune finită (varietate inerțială) de soluţii ale exemplului model [8]. Teoretic, Foiaş a ajuns să definească această varietate stabilind margini superioare ale numărului modurilor necesare să descrie (pentru t t t rarr oot \rightarrow \inftyt ) mulțimea soluțiilor către care tinde mulţimea tuturor soluţiilor modelului NavierStokes.
Uneori, exemplele model conțin tot ecuaţii cu derivate parţiale dar mult mai simple. In alte situații, echivalența ecuației cu derivate parțiale cu cea diferențială ordinară se referă nu la proprietăți asimptotice pentru t t t rarr oot \rightarrow \inftyt (ca în cazul problemelor de stabilitate și atractivitate menționate) ci λ 0 λ 0 lambda rarr0\lambda \rightarrow 0λ0 sau λ λ lambda rarr oo\lambda \rightarrow \inftyλ unde λ λ lambda\lambdaλ este un parametru fizic.

MATEMATICA EXPERIMENTALĂ ŞI MATEMATICA APLICATĂ

Matematica experimentală nu este o știință, ci o îndeletnicire care furnizează rezultate numerice deduse din rezolvarea la calculator a unor probleme pentru care nu se dispune de suficiente informații calitative riguroase, privind în special, existența și unicitatea soluției. Ea a fost impusă de cerințe practice stringente și s-a dezvoltat mult în ultimul timp, interferînd cu arta și matematica [9]. Descrierea completă a dinamicii de fază, chiar în cazul exemplelor model, este imposibilă în absența calculatorului. Dar, pe de altă parte, interpretarea corectă a rezultatelor furnizate de calculator devine chiar principial imposibilă în cazul unor ecuaţii care depind sensibil de data inițială. In decelarea erorii calculatorului de micile variații ale soluției, un rol din ce în ce mai mare revine matematicienilor care posedă cunoștințe de teoria modelării. În anii noștri, periodic au loc congrese asupra modelării, se țin cursuri de modelare la marile instituţii de învățămînt superior, apar cărți dedicate acestui subiect [10], [11]. Matematica aplicată, prin știința modelării, este cea care sesizează mutațiile importante ale științei contemporane. In particular ea orientează studiul turbulenței în direcții principial noi în care s-a angajat cercetarea mondială: legătura între laminar și turbulent via haos determinist al sistemelor dinamice.

SISTEME DINAMICE

Sistemul dinamic este un concept geometric (de topologie diferențială) care modelează mișcarea (evoluția) unui sistem material cu ajutorul unei mulțimi infinite de transformări indexate după un număr natural n n nnn (în care caz sistemul se numeşte discret) sau după un parametru real t t ttt numit timp (în care caz sistemul dinamic este continuu). Deci, un sistem dinamic discret este o aplicative T : N C T : N C T:NrarrCT: \mathbb{N} \rightarrow \mathcal{C}T:NC, unde C C C\mathcal{C}C este mulțimea transformărilor T n : M M T n : M M T_(n):M rarr MT_{n}: M \rightarrow MTn:MM ale spatiului fazelor M M MMM care au proprietățile T n + m == T n T m m , n N T n + m == T n T m m , n N T_(n+m)==T_(n)@T_(m)m,n inNT_{n+m}= =T_{n} \circ T_{m} m, n \in \mathbb{N}Tn+m==TnTmm,nN iar T 0 T 0 T_(0)T_{0}T0 este transformarea identitate. T T TTT duce un număr natural n n nnn într-o funcţie T n T n T_(n)T_{n}Tn și anume T ( n ) = T n T ( n ) = T n T(n)=T_(n)T(n)=T_{n}T(n)=Tn. La rîndul ei funcția T n T n T_(n)T_{n}Tn duce orice element din M M MMM intr-un element din M M MMM. Aplicaţiile de interval g : [ a , b ] [ a , b ] , [ a , b ] R g : [ a , b ] [ a , b ] , [ a , b ] R g:[a,b]rarr[a,b],[a,b]subRg:[a, b] \rightarrow[a, b],[a, b] \subset \mathbb{R}g:[a,b][a,b],[a,b]R, generează sisteme dinamice discrete asociate unei rezolvări iterative a ecuației x = g ( x ) x = g ( x ) x=g(x)x=g(x)x=g(x) și anume T n ( x ) = g ( x n ) T n ( x ) = g x n T_(n)(x)=g(x_(n))T_{n}(x)=g\left(x_{n}\right)Tn(x)=g(xn), deci T n ( x ) = g n ( x 0 ) T n ( x ) = g n x 0 T_(n)(x)=g^(n)(x_(0))T_{n}(x)=g^{n}\left(x_{0}\right)Tn(x)=gn(x0) unde x 0 x 0 x_(0)x_{0}x0 este valoarea inițială din procesul iterativ, iar M = [ a , b ] M = [ a , b ] M=[a,b]M=[a, b]M=[a,b]. Aflarea soluției acestei ecuații revine la găsirea punctului fix ** al aplicației g g ggg. Dacă g g ggg este contracție, atunci un procedeu de aproximaţii succesive în care x n = g n ( x 0 ) x n = g n x 0 x_(n)=g^(n)(x_(0))x_{n}=g^{n}\left(x_{0}\right)xn=gn(x0) converge către x x x^(**)x^{*}x deci lim n x n = x lim n x n = x lim_(n rarr oo)x_(n)=x^(**)\lim _{n \rightarrow \infty} x_{n}=x^{*}limnxn=x pentru orice alegere a lui x 0 x 0 x_(0)x_{0}x0. In acest procedeu apare şirul x 0 , x 1 , x 2 , x 0 , x 1 , x 2 , x_(0),x_(1),x_(2),dotsx_{0}, x_{1}, x_{2}, \ldotsx0,x1,x2,, numit orbită sau traiectoria lui x 0 x 0 x_(0)x_{0}x0, definit deci ca T 0 ( x 0 ) , T 1 ( x 0 ) , T 2 ( x 0 ) , T 0 x 0 , T 1 x 0 , T 2 x 0 , T_(0)(x_(0)),T_(1)(x_(0)),T_(2)(x_(0)),dotsT_{0}\left(x_{0}\right), T_{1}\left(x_{0}\right), T_{2}\left(x_{0}\right), \ldotsT0(x0),T1(x0),T2(x0),, reprezentînd locul pozițiilor lui x 0 x 0 x_(0)x_{0}x0, cînd timpul discret n n nnn crește de la zero la oo\infty și deci descriind dinamica discretă a lui x 0 x 0 x_(0)x_{0}x0 în R R R\mathbb{R}R. Considerații geometrice simple arată că grafic x n + 1 x n + 1 x_(n+1)x_{n+1}xn+1 se construieşte ducînd din x n x n x_(n)x_{n}xn o perpendiculară pe axa x x xxx care înțeapă graficul lui g g ggg in P 1 P 1 P_(1)P_{1}P1. Din P 1 P 1 P_(1)P_{1}P1 se duce o paralelă
la axa x x xxx care ințeapă prima bisectoare în P 2 P 2 P_(2)P_{2}P2; proiectia lui P 2 P 2 P_(2)P_{2}P2 pe axa x x xxx este tocmai x n + 1 x n + 1 x_(n+1)x_{n+1}xn+1 (fig. 4). Orbita lui x 0 x 0 x_(0)x_{0}x0 este un şir de puncte astfel construite de pe axa x x xxx. Ea este descrisă cu ajutorul valorilor în x 0 x 0 x_(0)x_{0}x0 ale tuturor transfor-
Fig. 4
mărilor T n T n T_(n)T_{n}Tn cu n 0 n 0 n >= 0n \geqslant 0n0 deci folosind întregul sistem dinamic T T TTT dar referitor numai la un singur element x 0 x 0 x_(0)x_{0}x0 al spatiului fazelor. Cu ajutorul lui T T TTT se poate descrie dinamica tuturor punctelor lui [ a , b ] [ a , b ] [a,b][a, b][a,b] deci se poate găsi mulțimea orbitelor tuturor punctelor din spatiul fazelor, mulțime numită portretul de fază al lui T T TTT. Produsul N × M N × M NxxM\mathbb{N} \times \mathbb{M}N×M se numește spațiul fazelor extins, iar șirul perechilor ( 0 , x 0 ) , ( 1 , x 1 ) , ( 2 , x 2 ) , 0 , x 0 , 1 , x 1 , 2 , x 2 , (0,x_(0)),(1,x_(1)),(2,x_(2)),dots\left(0, x_{0}\right),\left(1, x_{1}\right),\left(2, x_{2}\right), \ldots(0,x0),(1,x1),(2,x2), se numeşte orbită din spațiul fazelor extins a lui x 0 x 0 x_(0)x_{0}x0. Deoarece T n T n T_(n)T_{n}Tn sint funcții, în spațiul fazelor extins orbitele nu se intersectează în schimb se pot interesecta în spațiul fazelor. Aspectul portretului de fază al unui sistem dinamic este dictat de prezența în cadrul acestui portret a unor orbite speciale, atractive. Orbita care pleacă din x 0 [ a , b ] x 0 [ a , b ] x_(0)in[a,b]x_{0} \in[a, b]x0[a,b] se numește partial atractivă pentru orbitele care pleacă din x 0 [ c , d ] [ a , b ] x 0 [ c , d ] [ a , b ] x_(0)^(')in[c,d]sub[a,b]x_{0}^{\prime} \in[c, d] \subset[a, b]x0[c,d][a,b] (unde [ c , d ] x 0 [ c , d ] x 0 [c,d]∋x_(0)[c, d] \ni x_{0}[c,d]x0 ) dacă distanța dintre ea și aceste orbite tinde la zero, deci dacă lim n | T n ( x 0 ) T n ( x 0 ) | = 0 lim n T n x 0 T n x 0 = 0 lim_(n rarr oo)|T_(n)(x_(0))-T_(n)(x_(0)^('))|=0\lim _{n \rightarrow \infty}\left|T_{n}\left(x_{0}\right)-T_{n}\left(x_{0}^{\prime}\right)\right|=0limn|Tn(x0)Tn(x0)|=0. [ c , d ] [ c , d ] [c,d][c, d][c,d] se numeşte bazin de atractie. Dacă bazinul [ c , d ] [ c , d ] [c,d][c, d][c,d] coincide cu întregul spaţiu al fazelor [ a , b a , b a,ba, ba,b ], atunci orbita ce porneşte din x 0 x 0 x_(0)x_{0}x0 se numeşte total atractivă sau, simplu, atractivă. Există o mulțime de alte tipuri de atractivitate sau de comportamente asimptotice pentru n n n rarr oon \rightarrow \inftyn ale orbitelor. Dintre orbitele tipice care pot fi atractive menţionăm orbitele periodice (sau ciclurile) de perioadă n n nnn sau n n nnn-ciclurile în care T m ( n 1 ) + k ( x 0 ) = T k ( x 0 ) T m ( n 1 ) + k x 0 = T k x 0 T_(m(n-1)+k)(x_(0))=T_(k)(x_(0))T_{m(n-1)+k}\left(x_{0}\right)=T_{k}\left(x_{0}\right)Tm(n1)+k(x0)=Tk(x0) deci de forma x 0 , x 1 , x 2 , x n , x 0 , x 1 , x 2 , , x n , x 0 , x 1 , x 2 , x n , x 0 , x 1 , x 2 , , x n , x_(0),x_(1),x_(2),dotsx_(n),x_(0),x_(1),x_(2),dots,x_(n),dotsx_{0}, x_{1}, x_{2}, \ldots x_{n}, x_{0}, x_{1}, x_{2}, \ldots, x_{n}, \ldotsx0,x1,x2,xn,x0,x1,x2,,xn, unde 0 k n 1 0 k n 1 0 <= k <= n-10 \leqslant k \leqslant n-10kn1, n N , m N .1 n N , m N .1 n inN^(**),m inN.1n \in \mathbb{N}^{*}, m \in \mathbb{N} .1nN,mN.1 - ciclul este orbita constantă x 0 , x 0 , x 0 , x 0 , x 0 , x 0 , x_(0),x_(0),x_(0),dotsx_{0}, x_{0}, x_{0}, \ldotsx0,x0,x0, in care, deci x 0 x 0 x_(0)x_{0}x0 este punct fix sau punct de echilibru al aplicatiei g g ggg care a generat pe T T TTT. Imaginea în spațiul fazelor a 1 -ciclului coincide cu imaginea punctului fix x 0 x 0 x_(0)x_{0}x0. Dacă un n n nnn-ciclu descrie mişcarea lui x 0 x 0 x_(0)x_{0}x0, atunci x 0 x 0 x_(0)x_{0}x0 este punct fix al lui g n g n g^(n)g^{n}gn și se numeşte punct periodic de perioadă n n nnn. In spațiul fazelor un n n nnn-ciclu se reprezintă prin n n nnn puncte. Punctul fix este un punct periodic de perioadă 1 . In cazul cînd g g ggg este o contracție, atunci sistemul dinamic T T TTT generat de g g ggg are orbita constantă x , x , x , x , x , x , x^(**),x^(**),x^(**),dotsx^{*}, x^{*}, x^{*}, \ldotsx,x,x, (unde x x x^(**)x^{*}x este punctul fix al lui g g ggg ) drept orbită total atractivă.
Mulțimea formată dintr-o orbită total atractivă sau din mulțimea tuturor orbitelor parțial atractive ale lui T T TTT se numeşte mulțime atractivă a lui T T TTT. Dacă mulțimea atractivă este local o varietate, atunci ea se numeşte atractor. Cînd T T TTT este generat de aplicația g g ggg iar aplicația g g ggg este o contracție, atunci mulțimea atractivă a lui T T TTT este un atractor compus dintr-un singur element (orbita constantă corespunzătoare punctului fix). Dacă g g ggg
nu este contracție, atunci atractorul poate fi format dintr-un număr finit sau infinit de orbite parțial áractive. O mulțime atractivă care nu este varietate se numeşte atractor straniu sau atractor haotic; ea conține o infinitate de orbite parțial atractive (i.e. al căror bazin de atracţie este strict conținut în M M MMM ) între care orbitele neatractive rătăcesc haotic. Caracterul determinant al dinamicii de fază rezultă din aceea că T n T n T_(n)T_{n}Tn sînt funcții, deci duc un x 0 x 0 x_(0)x_{0}x0 într-un singur T n ( x 0 ) T n x 0 T_(n)(x_(0))T_{n}\left(x_{0}\right)Tn(x0). Dar cu tot acest caracter determinist, cînd T T TTT posedă un atractor straniu dinamica de fază este haotică. Prin dinamică haotică se înțelege un portret de fază format din orbite complicat de descris (din cauza ordinei aproape întimplătoare a termenilor săi), orbite care tind indecis cînd către o orbită a atractorului straniu, cînd către alta, nu părăsesc vecinătatea acestuia, dar nici nu tind către el după o regulă simplă, deşi această tindere este complet specificată de valorile T n ( x 0 ) T n x 0 T_(n)(x_(0))T_{n}\left(x_{0}\right)Tn(x0). Deci, mişcarea unui punct x 0 M x 0 M x_(0)in Mx_{0} \in Mx0M este haotică dacă, deși deterministă, ea are proprietăți nepredictibile (stocastice). Această dualitate: determinism pe de o parte (deoarece T n T n T_(n)T_{n}Tn sint funcții) și impredictibilitatea practică (datorită complexității lui T n ( x 0 ) = g n ( x 0 ) T n x 0 = g n x 0 T_(n)(x_(0))=g^(n)(x_(0))T_{n}\left(x_{0}\right)=g^{n}\left(x_{0}\right)Tn(x0)=gn(x0) ) se numește haos determinist [12].
Se remarcă că acelaşi tip de haos determinist, deşi într-o măsură mult mai mică, era cuprins și în funcția sin x 1 x 1 x^(-1)x^{-1}x1. Dar acest haos nu este încă turbulență, deoarece el este teoretic predictibil, determinist, în timp ce prin definiție dinamică turbulentă este caracterizată de mărimi aleatoare. Cu toate acestea, haosul determinist este un capăt de drum spre turbulență, deoarece dinamica de fază este impredictibilă practic, dar are anumite proprietăți statistice comune cu turbulența. In acest fel, dacă sistemul dinamic T T TTT depinde de un parametru real λ λ lambda\lambdaλ și, pe măsură ce λ λ lambda\lambdaλ creşte, mulțimea atratică a lui T T TTT se modifică, din atractor devenind atractor straniu, rezultă că odată cu creşterea lui λ λ lambda\lambdaλ se parcurge un drum spre turbulență și acest drum este un scenariu care depinde de sistemul dinamic T T TTT concret. S-a constatat că există mai multe scenarii tipice pe care le vom expune la sfîrşitul lucrării după ce vom fi introdus şi alte noțiuni.
Sistemele dinamice discrete pot fi generate și de funcţi complexe, de variabilă complexă. Mult studiate sînt funcțiile raționale (cît de două polinoame) care posedă atractori bogați și dinamici haotice.
In cazul sistemelor dinamice continue numite și curgeri, toate definițiile și rezultatele de mai sus se transcriu corespunzător.
Dimensiunea unti sistem dinamic este dată de dimensiunea spatiului fazelor. Sistemele dinamice continue finit dimensionale sînt generate de soluții ale problemei Cauchy, pentru care soluţia există şi este unică,
(1) x ( 0 ) = x 0 , (1) x ( 0 ) = x 0 , {:(1) vec(x)(0)= vec(x)_(0)",":}\begin{equation*} \vec{x}(0)=\vec{x}_{0}, \tag{1} \end{equation*}(1)x(0)=x0,
pentru o ecuație (vectorială) diferențială ordinară, autonomă
(2) d x d t = f ( x ) (2) d x d t = f ( x ) {:(2)(dx)/((d)t)= vec(f)( vec(x)):}\begin{equation*} \frac{\mathrm{d} x}{\mathrm{~d} t}=\vec{f}(\vec{x}) \tag{2} \end{equation*}(2)dx dt=f(x)
( f f vec(f)\vec{f}f se numeşte cimp de vectori), dacă (1), (2) are soluție unică şi anume
(3) T l ( x 0 ) = x ( t , x 0 ) (3) T l x 0 = x t , x 0 {:(3)T_(l)( vec(x)_(0))= vec(x)(t, vec(x)_(0)):}\begin{equation*} T_{l}\left(\vec{x}_{0}\right)=\vec{x}\left(t, \vec{x}_{0}\right) \tag{3} \end{equation*}(3)Tl(x0)=x(t,x0)
unde x ( t , x 0 ) x t , x 0 vec(x)(t, vec(x)_(0))\vec{x}\left(t, \vec{x}_{0}\right)x(t,x0) este soluţia acelei probleme. De fapt, sistemul dinamic este generat de mulțimea acestor probleme corespunzătoare tuturor datelor x 0 M x 0 M vec(x)_(0)in M\vec{x}_{0} \in Mx0M sau, echivalent, de problema Cauchy a cărei dată inițială descrie pe M M MMM. Invers, oricărui sistem dinamic T T TTT, ale cărui valori T t T t T_(t)T_{t}Tt sînt definite de (3), i se poate ataşa o problemă Cauchy (1), (2) a cărei soluție să genereze pe T T TTT. Specificarea lui x 0 x 0 x_(0)x_{0}x0 în scrierea soluției lui (1), (2) arată că acea soluție corespunde datei initiale x 0 x 0 vec(x)_(0)\vec{x}_{0}x0. Aceeaşi scriere uşurează și o interpretare dinamică bazată pe definiția (3) : orbita (traiectoria) lui x 0 x 0 vec(x)_(0)\vec{x}_{0}x0 este imaginea funcţiei x ( t , x 0 ) x t , x 0 vec(x)(t, vec(x)_(0))\vec{x}\left(t, \vec{x}_{0}\right)x(t,x0) de t t ttt (cu x 0 x 0 vec(x)_(0)\vec{x}_{0}x0 fixat). Cînd variază atît t t ttt cît și x 0 x 0 vec(x)_(0)\vec{x}_{0}x0 se obține portretul de fază al lui T T TTT și, corespunzător, mulțimea soluțiilor tuturor problemelor Cauchy ale căror date inițiale sînt în spaţiul fazelor M M MMM, care este izomorf cu R n R n R^(n)\mathbb{R}^{n}Rn și deci îl alegem drept R n R n R^(n)\mathbb{R}^{n}Rn. Traiectoriile de fază ale sistemelor dinamice continue sînt deci curbe din R n R n R^(n)\mathbb{R}^{n}Rn, variabila independentă t t ttt a problemei (1), (2) jucînd rol de parametru pe acele traiectorii. Caracterizarea de continuu atribuită lui T T TTT se datoreşte faptului că traiectoriile sale sînt continuu muri matematice (fiind curbe). 1 -ciclurilor din cazul discret le corespund aici soluțiile staţionare ale lui (2) deci vectorii constanți x R x R x^(**)inRx^{*} \in \mathrm{R}xR, numiți puncte de echilibru, sau, simplu, echilibre, care satisfac ecuaţia
(4) f ( x ) = 0 (4) f x = 0 {:(4) vec(f)( vec(x)^(**))=0:}\begin{equation*} \vec{f}\left(\vec{x}^{*}\right)=0 \tag{4} \end{equation*}(4)f(x)=0
un analog al k k kkk-ciclurilor fiind torurile T k 1 T k 1 T^(k-1)T^{k-1}Tk1 din R n R n R^(n)\mathbb{R}^{n}Rn. Curbele de pe T k T k T^(k)T^{k}Tk sînt k k kkk-periodice în raport cu t t ttt. In particular, T 1 T 1 T^(1)T^{1}T1 este orice curbă închisă topologic echivalentă cu cercul și deseori este numit ciclu limită.
Deşi, Poincaré a fost primul care a observat că sistemele dinamice neliniare (i.e. c u f c u f cu vec(f)c u \vec{f}cuf neliniar) pot avea mulțimi atractive bizare (i.e. atractorii stranii), totuși, incepind cu el, o atenție deosebită a fost acordată unor orbite atractive speciale: punctele de echilibru și ciclurile limită. Indepărtarea atenției de la atractorii stranii a avut un efect negativ asupra dezvoltării teoriei sistemelor dinamice; probabil însă că aceasta nu putea fi evitat, deoarece aceşti atractori intrau în categoria mulţimilor ciudate, ocolite de toţi matematicienii. Insuși Cantor cînd a dedus mulțimea care-i poartă numele a exclamat: „văd dar nu cred". Abia după anii 1975, odată cu introducerea fractalilor de către Mandelbrot, mulțimile neobișnuite au început să fie studiate sistematic în cadrul teoriei geometrice a măsurii.
Spre deosebire de cazul discret în care traiectoriile de fază sînt şiruri de puncte discrete, in sistemele dinamice continue traiectoriile de fază sînt curbe. Acestea se pot intersecta în spațiul fazelor dar nu in spațiul fazelor extins, care acum este R × M R × M Rxx M\mathbb{R} \times MR×M, deoarece T t T t T_(t)T_{t}Tt sint funcţii și, corespunzător, problema Cauchy (1), (2) are soluție unică. Mişeării pe traiectorie a unui singur punct x 0 x 0 vec(x)_(0)\vec{x}_{0}x0 din spațiul fazelor î corespunde variația în timp a soluţiei x ( t , x 0 ) x t , x 0 vec(x)(t, vec(x)_(0))\vec{x}\left(t, \vec{x}_{0}\right)x(t,x0) a problemei Cauchy ataşată in care x 0 x 0 vec(x)_(0)\vec{x}_{0}x0 este fixat. O problemă Cauchy pentru o ecuație diferențială ordinară neautonomă de tipul (3) în care f = f ( t , x ) f = f ( t , x ) vec(f)= vec(f)(t, vec(x))\vec{f}=\vec{f}(t, \vec{x})f=f(t,x) generează, dacă soluția sa există și este unică un sistem dinamic n + 1 n + 1 n+1n+1n+1 dimensional, deoarece in acest caz t t ttt se ia drept
Ca orice funcție, un sistem dinamic T T TTT este cunoscut cînd se cunose valorile sale, adică transformările T t T t T_(t)T_{t}Tt pentru t R t R AA t inR\forall t \in \mathbb{R}tR, sau, ținînd seama de (3), cînd se cunoaşte explicit dependența soluției x x vec(x)\vec{x}x de x 0 x 0 vec(x)_(0)\vec{x}_{0}x0. Aceasta este posibil doar cînd f f vec(f)\vec{f}f este funcție liniară sau afină de x x vec(x)\vec{x}x. De exemplu, în cazul afin f ( x ) = A x + B f ( x ) = A x + B vec(f)( vec(x))=A vec(x)+B\vec{f}(\vec{x})=A \vec{x}+Bf(x)=Ax+B, unde A A AAA este o matrice n × n n × n n xx nn \times nn×n, iar B B BBB o matrice n × 1 , A n × 1 , A n xx1,An \times 1, An×1,A și B B BBB depinzînd de t t ttt, soluţia x ( t , t 0 , x 0 ) x t , t 0 , x 0 vec(x)(t,t_(0), vec(x)_(0))\vec{x}\left(t, t_{0}, \vec{x}_{0}\right)x(t,t0,x0) a problemei (1), (2) corespunzătoare în care x ( t 0 ) = x 0 x t 0 = x 0 vec(x)(t_(0))= vec(x)_(0)\vec{x}\left(t_{0}\right)=\vec{x}_{0}x(t0)=x0 este dată de formula de variație a constantelor
x ( t , t 0 , x 0 ) = C ( t , t 0 ) x 0 + t 0 t C ( t , s ) B ( s ) d s x t , t 0 , x 0 = C t , t 0 x 0 + t 0 t C ( t , s ) B ( s ) d s vec(x)(t,t_(0), vec(x)_(0))=C(t,t_(0)) vec(x)_(0)+int_(t_(0))^(t)C(t,s)B(s)ds\vec{x}\left(t, t_{0}, \vec{x}_{0}\right)=C\left(t, t_{0}\right) \vec{x}_{0}+\int_{t_{0}}^{t} C(t, s) B(s) \mathrm{d} sx(t,t0,x0)=C(t,t0)x0+t0tC(t,s)B(s)ds
unde C ( t , t 0 ) C t , t 0 C(t,t_(0))C\left(t, t_{0}\right)C(t,t0) este matricea fundamentală a lui (2) și are drept elemente C i j = x i ( t , t 0 , e j ) C i j = x i t , t 0 , e j C_(ij)=x_(i)(t,t_(0), vec(e)_(j))C_{i j}=x_{i}\left(t, t_{0}, \vec{e}_{j}\right)Cij=xi(t,t0,ej). In acest caz, sistemul dinamic continuu T T TTT generat de
(1), (2) are valorile T t = C ( t , t 0 ) + t 0 t C ( t , s ) B ( s ) d s T t = C t , t 0 + t 0 t C ( t , s ) B ( s ) d s T_(t)=C(t,t_(0))+int_(t_(0))^(t)C(t,s)B(s)dsT_{t}=C\left(t, t_{0}\right)+\int_{t_{0}}^{t} C(t, s) B(s) \mathrm{d} sTt=C(t,t0)+t0tC(t,s)B(s)ds.
Următoarea problemă Cauchy, clasică, pentru un sistem de ecuații cu derivate partiale
(5) lim t 0 u ( t , x ) = u 0 ( x ) , d u d t = N ( x , u ) (5) lim t 0 u ( t , x ) = u 0 ( x ) , d u d t = N ( x , u ) {:(5)lim_(t rarr0) vec(u)(t"," vec(x))= vec(u)_(0)( vec(x))","(d( vec(u)))/((d)t)= vec(N)( vec(x)"," vec(u)):}\begin{equation*} \lim _{t \rightarrow 0} \vec{u}(t, \vec{x})=\vec{u}_{0}(\vec{x}), \frac{\mathrm{d} \vec{u}}{\mathrm{~d} t}=\vec{N}(\vec{x}, \vec{u}) \tag{5} \end{equation*}(5)limt0u(t,x)=u0(x),du dt=N(x,u)
în care N N vec(N)\vec{N}N este un operator în general neliniar, generează un sistem dinamic continuu infinit dimensional T T TTT unde T t ( u 0 ) = u ( t , x ; u 0 ) , T t T t u 0 = u t , x ; u 0 , T t T_(t)( vec(u)_(0))= vec(u)(t,( vec(x)); vec(u)_(0)),T_(t)T_{t}\left(\vec{u}_{0}\right)=\vec{u}\left(t, \vec{x} ; \vec{u}_{0}\right), T_{t}Tt(u0)=u(t,x;u0),Tt sint definite pe spatiul fazelor infinit dimanional M M MMM al functiilor netede u 0 ( x ) u 0 ( x ) vec(u)_(0)( vec(x))\vec{u}_{0}(\vec{x})u0(x) care satisfac anumite condiții la limită pe Ω Ω del Omega\partial \OmegaΩ. Cu toate că explicitarea dependenței lui T t T t T_(t)T_{t}Tt de u 0 u 0 vec(u)_(0)\vec{u}_{0}u0 este în general imposibilă, simpla existență a acestor funcții și anumite proprietăți ale lor le fac deosebit de folositoare in teoria calitativă. Problemelor fizicii matematice li se ataşează sisteme dinamice mai ales pentru acest studiu calitativ, care în afara existenței, unicității și stabilității privește și atractivitatea, bifurcația și comportamentul asimptotic al soluțiilor acelor probleme. Sistemele dinamice fiind obiecte geometrice conduc la denumirea de teorie geometrică a ecuațiilor diferențiale pentru teoria calitativă.
Punerea problemelor Cauchy sub formă de sistem dinamic este o problemă matematică deosebit de dificilă.
Problema bifurcației mulțimii soluţilor unei probleme pentru o ecuație diferențială apare cînd acele probleme depind de un parametru fizic λ λ lambda\lambdaλ sau de mai mulți. Bifurcația mulțimii soluțiilor problemelor stationare se numeşte bifurcație statică și analizează dependența acelei mulțimi de parametru (dependență care nu este funcție de λ λ lambda\lambdaλ, ci funcție multiformă de λ λ lambda\lambdaλ ). Teoriile de bifurcație statică furnizează poziția și numărul punctelor de bifurcație, numărul și forma varietăților bifurcate, mărginirea acestora, comportamentul lor pentru valori deosebite ale lui λ λ lambda\lambdaλ. Studiul în jurul punctelor de bifurcaţie se numeşte local; cel al varietăţilor
bifurcate este un studiu global. Bifurcația mulțimii S a soluțiilor problemelor nestaționare se numeşte bifurcație dinamică. λ 0 λ 0 lambda_(0)\lambda_{0}λ0 este punct de bifurcație pentru S S SSS dacă la traversarea lui λ 0 , S λ 0 , S lambda_(0),S\lambda_{0}, Sλ0,S își schimbă proprietățile topologice, deci dacă protretul de fază pentru λ < λ 0 λ < λ 0 lambda < lambda_(0)\lambda<\lambda_{0}λ<λ0 nu este topologic echivalent cu cel pentru λ > λ 0 λ > λ 0 lambda > lambda_(0)\lambda>\lambda_{0}λ>λ0. Spre deosebire de cazul static o bifurcație dinamică poate avea loc cînd se modifică doar proprietățile atractive ale atractorului sistemului dinamic ataşat, dar nu și componența acelui atractor.

PROPRIETĂȚI STATISTICE ALE SISTEMELOR DINAMCE

In afara sistemelor dinamice definite există unele mai generale studiate de teoria sistemelor dinamice, parte a topologiei diferentiale. Topologia diferenţială alături de analiza pe varietăţi și alte ramuri matematice alcătuiesc analiza globală în cadrul cărora se studiază azi evoluțiile complexe ale sistemelor. Topologia studiază aplicațiile dintre spații; cînd aceste aplicații sînt diferențiale atunci ele formează obiectul topologiei diferențiale. Ó ramură importantă a analizei globale este teoria ergodică care studiază mişcarea în spații, mai generale, cu măsură M ~ = ( M , Σ , μ ) M ~ = ( M , Σ , μ ) widetilde(M)=(M,Sigma,mu)\widetilde{M}=(\mathrm{M}, \Sigma, \mu)M~=(M,Σ,μ) unde M M MMM este spațiul fazelor sistemului dinamic T , Σ T , Σ T,SigmaT, \SigmaT,Σ este o σ σ sigma\sigmaσ-algebră de submulţimi ale lui M M MMM, iar μ μ mu\muμ este o măsură definită pe Σ Σ Sigma\SigmaΣ. Spațiile cu măsură sînt spaţii ale fazelor pentru obiecte T T TTT mai generale decît T T TTT, numite de asemenea, în teoria ergodică, sisteme dinamice. Din punctul de vedere al teoriei probabilităților M M MMM poate fi interpretat ca spațiu al realizărilor unui proces aleator, iar măsura μ μ mu\muμ pe M M MMM ca distribuție de probabilitate pentru acest proces. Un sistem dinamic se zice ergodic dacă μ ( A ) = 0 μ ( A ) = 0 mu(A)=0\mu(A)=0μ(A)=0 sau 1 pentru orice mulțime invariantă A A AAA. Dacă T ~ T ~ tilde(T)\tilde{T}T~ este un sistem dinamic ergodic atunci pentru orice x M x M x in Mx \in MxM, media în timp este egală cu media spațială, rezultat cunoscut sub numele de teorema ergodică și folosit în teoria lui Reynolds, deşi această teoremă n-a fost demonstrată pentru modelul Navier-Stokes. Mai demult, teoria ergodică privea în principal teoremele de mediere și problemele calitative; se credea că spectrul este principalul invariant al sistemelor dinamice. Ulterior, s-a constatat că entropia este principalul invariant al sistemelor dinamice, teoria ergodică devenind in ultimii 25 de ani, o știință care studiază proprietățile statistice ale sistemelor dinamice [13], cu ajutorul notiunilor fundamentale de ergodicitate, amestec, spectru, entropie. În cadrul său se studiază, în mod natural, proprietățile probabiliste ale sistemelor deterministe [10].
Studiul portretului de fază presupune determinarea : 1) dimensiunii atractorului (dimensiume care este de diferite tipuri : capacitatea termenilor, dimensiunea Hausdorff, dimensiunea fractală, informația, dimensiunea lui Renyi, dimensiunea lui Liapunov) ; 2) entropia (topologică, metrică, de tip Kolmogorov-Sinai) ; 3) spectrul exponenților Liapunov ; 4) funcții de distribuție de probabilitate pentru secțiunile și aplicațiilePoincaré.
De exemplu, în spații metrice o dimensiune metrică este capacitatea d c d c d_(c)d_{c}dc : fie A A AAA o mulţime mărginită din R n ; R n ; R^(n);\mathbb{R}^{n} ;Rn; atunci d c ( A ) = lim ε 0 ln N ( ε ) ln ( 1 / ε ) d c ( A ) = lim ε 0 ln N ( ε ) ln ( 1 / ε ) d_(c)(A)=lim_(epsi rarr0)(ln N(epsi))/(ln(1//epsi))d_{c}(A)=\lim _{\varepsilon \rightarrow 0} \frac{\ln N(\varepsilon)}{\ln (1 / \varepsilon)}dc(A)=limε0lnN(ε)ln(1/ε), unde N ( ε ) N ( ε ) N(epsi)N(\varepsilon)N(ε) este numărul minim de curburi n n nnn-dimensionale de latură ε ε epsi\varepsilonε care pot acoperi pe A A AAA. Dacă A A AAA este mulţimea lui Cantor, atunci
d c ( A ) = ln 2 ln 3 0 , 63 d c ( A ) = ln 2 ln 3 0 , 63 d_(c)(A)=(ln 2)/(ln 3)~~0,63d_{c}(A)=\frac{\ln 2}{\ln 3} \approx 0,63dc(A)=ln2ln30,63. Pentru diferite tipuri de mulţimi se pot defini diferite tipuri de dimensiuni. Atractorii pot fi caracterizaţi cu ajutorul dimensiunii topologice care este un număr intreg egal cu dimensiunea spatitului euclidian local izomorf cu acea varietate. Tipurile celelalte de dimensiuni sînt necesare pentru atractorii haotici pentru care nu se poate defini dimensiunea topologică, deci acele dimensiuni ale atractorilor haotici sînt numere neîntregi. De exemplu, atractorul straniu al lui Lorenz este o mulțime de puncte din R 3 R 3 R^(3)\mathbb{R}^{3}R3 care seamănă cu o suprafaţă, dar care are o anumită grosime mică, dimensiunea sa Hausdorff fiind 2,06 .
Faptul că un atractor este straniu rezultă din aceea că, cel mai mare exponent Liapunov λ = 0 1 P ( X ) [ ln | d g d X | ] d X λ = 0 1 P ( X ) ln d g d X d X lambda=int_(0)^(1)P(X)[ln|(dg)/((d)X)|]dX\lambda=\int_{0}^{1} P(X)\left[\ln \left|\frac{\mathrm{d} g}{\mathrm{~d} X}\right|\right] \mathrm{d} Xλ=01P(X)[ln|dg dX|]dX, unde P ( X ) P ( X ) P(X)P(X)P(X) este probabilitatea de a găsi o aplicație g g ggg în intervalul ( X , X + d X X , X + d X X,X+dXX, X+\mathrm{d} XX,X+dX ), este pozitiv. Pentru caracterizarea stocasticitătii atractorilor haotici cu ajutorul exponenţilor lui Liapunov recomandăm [14].
O proprietate importantă a atractorilor stranii este aceea că mişcarea asimptotică spre ei este haotică (aparent) și sensibilă de punctul de plecare (data inițială), această haoticitate (exprimată prin proprietățile probabilistice menţionate) fiind intrinsecă, fără aport de zgomot. Atractorii stranii sînt structural stabili (adică orbitele lor nu sînt afectate semnificativ de mici perturbaţii). Posibilitatea apariției atractorilor stranii este importantă, deoarece prin aceasta, comportarea pe timp lung a soluțiilor este de natură aparent aleatoare, ceea ce nu este de asteptat din punctul de vedere traditional al sistemelor dinamice. Această impredictibilitate pe termen lung este o consecință a proprietății de hiperbolicitate a sistemului dinamic, deci de aplatisare a spatitului fazelor in general in infinit de multe direcţii. De această proprietate este legată noțiunea de varietate inerţială.

MULTIMI FRAGTALE

Mulțimile fractale sau, pe scurt, fractalii, sînt mulțimi care nu pot fi caracterizate relevant de o dimensiune topologica a lor sau a oricărei părți conexe a lor, ci de o dimensiune Hausdorff fracționară. Ele sînt neobişnuite, parțial sau total discrete, pot fi sub formă de curbe dense, suprafețe subțiri, pot avea arie nulă dar conțin linii în orice directie și multe alte ciudățenii. Cu toate acestea ele sînt intilnite peste tot, formează regula și nu excepția. Neregularitatea lor caracteristică îndepărtează de ele pe cel obişnuit cu simetriile şi formele simple ale geometriei euclidiene : punct, dreaptă, plan, spațiu euclidian. Noțiunea de fractal a fost introdusă în 1975 de B. Mandelbrot drept mulțime care are dimensiunea Hausdorff mai mare ca dimensiunea topologică (a oricărei părţi conexe) a sa. Cei mai mulți fractali nu sînt varietăți. Cel mai simplu fractal este mulțimea triadică a lui Cantor. Alte exemple : prafurile dense, traiectoriile browniene, coastele mărilor și oceanelor, crestele munților, frontierele bazinelor de atracție ale anumitor aplicații. Unele mulțimi fractale sînt autosimilare
(asemenea cu părți ale lor) ; fractalii autosimilari pot fi generaţi iterativ. Deoarece, atractorii stranii cunoscuţi nu sînt varietăți şi posedă dimensiune Hausdorff fracționară, ei sînt fractali, de aceea atractorii haotici se mai numese si atractori fractali.

SCENARII DE TURBULENŢĂ

In experimentele numerice și în multe studii teoretice au fost puse in evidenţă drumuri tipice de trecere spre turbulență cînd parametrul λ λ lambda\lambdaλ al sistemului dinamic creşte. De obicei în orice evoluție se întîlnesc mai multe astfel de scenarii.
Scenariul lui Landau şi Hopf (L-H) de bifurcații statice repetate (fig. 5) : există un număr infinit de puncte de bifurcație ale mulțimii soluțiilor unei probleme staţionare, în fiecare astfel de puncte aparind noi
Fig. 5
soluții liniar stabile pe seama destabilizării soluțiilor din care au provenit. La turbulenţă se ajunge cu o mişcare cvaziperiodică. In particular, acest scenariu este observat in mişcarea Couette și in convecțiile termice. Acest scenariu a fost sugerat de modelul lui Burgers [15], apărînd și în [16].
Scenariul lui Ruelle si Takens (R-T): trecerea spre turbulență implică numai trei bifturcații. După prima apare o orbită periodică printr-o bifurcație Hopi, după cea de a doua apare o orbită cvaziperiodică (tor), iar după cea de a treia apare un atractor haotic. La fiecare bifurcație are loc și o pierdere de stabilitate [17]. Experimental acest scenariu a fost pus în evidență în particular în mişearea lui Couette. Acest scenariu a fost găsit prin experimente numerice, incă nefundamentate, asupra sistemului dinamic tridimensional al lui Lorenz. Remarcăm număral finit de bifurcații necesare pentru turbulentizarea (haoticizarea) orbitelor sistemului dinamic prevăzut de acest scenariu.
Dintre inconvenientele scenariului L H L H L-H\mathrm{L}-\mathrm{H}LH semnalăm: conduce la o statistică gaussiană care nu se observă în realitate ; complexitatea dinamicii este gradată spre deosebire de observația experimentală că turbulența apare la un anumit număr Reynolds bine definit; starea turbulentă prezisă de acest scenariu nu depinde sensibil de datele initiale, contrar observațiilor. Se constată prin comparație că aceste neajunsuri sînt inlăturate de scenariul R-T. Conform acestui scenariu, trecerea spre turbulență se face prin mişeări aperiodice (atractorii) ce posedă anumite proprietăți statistice. Scenariile L-H și R-T se referă la sisteme dinamice continue.
Scenariul lui Feigenbaum (F) al dublării perioadei se referă la sisteme dinamice discrete și a fost pus in evidență in 1978 prin experimente numerice [18] pentru cazul functiei logistice g ( x ) = 4 λ × ( 1 x ) g ( x ) = 4 λ × ( 1 x ) g(x)=4lambda xx(1-x)g(x)=4 \lambda \times(1-x)g(x)=4λ×(1x), unde λ [ 0 , 1 ] λ [ 0 , 1 ] lambda in[0,1]\lambda \in[0,1]λ[0,1] este un parametru real, iar x [ 0 , 1 ] x [ 0 , 1 ] x in[0,1]x \in[0,1]x[0,1]. Scenariul F F FFF constă dintr-o acumulare de bifurcații de mulțimi de orbite a căror perioadă se dublează : un punct de echilibru care în punctul de bifurcație își pierde atractivitatea este înlocuit cu un 2 -ciclu atractiv. Acesta își pierde atractivitatea în alt punct de bifurcație fiind înlocuit de un 2 2 2 2 2^(2)2^{2}22-ciclu atractiv și aşa mai departe, din orbite periodice de ordin 2 k 1 2 k 1 2^(k-1)2^{k-1}2k1, luînd naștere orbite periodice de ordin 2 k 2 k 2^(k)2^{k}2k, iar k k k rarr ook \rightarrow \inftyk. Scenariul F a fost pus in evidență in anumite experimente de convecție. După unii autori, acest scenariu este al lui Grossmann, May și Feigenbaum sau al lui Verhulst.
În experimente, în oscilatori electrici neliniari, în reacția BelousovJabotinski, în convecția Rayleigh-Bénard intr-un cîmp magnetic apare scenariul ( U U UUU ) universal și anume dincolo de punctul de acumulare ( 2 2 2^(oo)2^{\infty}2 ciclul) din scenariul F. Cu cîțiva ani înainte ca proprietățile de scalare universală ale aplicațiilor unidimensionale să fi fost descoperite de Feigenbaum, Metropolis și Stein au găsit că aplicațiile unidimensionale cu un singur punct de extrem posedă o dinamică universală, funcție de parametrul de bifurcație : după șirul dublărilor de perioadă a unei stări periodice cu o oscilație pe perioadă apar stări periodice cu K K KKK oscilații pe perioadă pentru toți K numere naturale și fiecare dintre aceste K K KKK-cicluri suferă şirul său propriu infinit de dublări de perioadă 2 n K 2 n K 2^(n)K2^{n} K2nK.
Scenariul lui Mandelbrot (M) a rezultat din experimente numerice [19] asupra sistemului dinamic generat de aplicația z z 2 C z z 2 C z rarrz^(2)-Cz \rightarrow z^{2}-\mathbb{C}zz2C unde z z zzz, c C c C c inCc \in \mathbb{C}cC, c fiind parametrul : în planul complex C C C\mathbb{C}C există o mulțime M M MMM conexă, formată dintr-un corp principal (o cardioidă) terminat cu muguri, aceștia reproducînd aproape asemenea cardioida și mugurii, această aproape similaritate repetîndu-se pînă la infinit. După poziția lui s s sss în diferitele porţiuni centrale sau laterale ale lui M M MMM, sistemul dinamic discret generat de aplicația menționată are un portret de fază sau altul, incluzînd mulțimi atractive mărginite de frontiere fractale conexe. Dacă c c ccc este în exteriorul lui M M MMM atunci această frontieră nu este conexă. De aceea, cînd c c ccc traversează frontiera spre exteriorul lui M M MMM frontiera bazinelor de atracție explodează, descompunîndu-se într-o infinitate de puncte discrete. Se conjecturează că această traversare modelează trecerea spre turbulență [20].
Scenariul Manneville-Pomeau a fost găsit in 1980 și precizează o tranziție intermitentă spre turbulență în care faza laminară se menține ca 1 / Re Re cr 1 / Re  Re  cr  1//sqrt(Re-" Re "_("cr "))1 / \sqrt{\operatorname{Re}-\text { Re }_{\text {cr }}}1/Re Re cr . Acest scenariu a fost găsit in experimente din convecții Rayleigh-Bénard și în circuite electronice [21].
Nici unul dintre scenariile de turbulență n-a fost demonstrat riguros pentru modelul Navier-Stokes.

BIBLIOGRAFIE

  1. ŞT. SĂVULESCU, Tranziția de là scurgerea laminară la cea turbulentă, Edit. Academiei Române, București, 1968.
  2. В.В. СТРУМИНСКИЙ, (ред.), Проблемы турбулентных течений, Наука, Москва, 1987.
  3. M. M. STANIS゙IČ, The mathematical theory of turbulence, Springer-Verlag, New York, 1958.
  4. C. FOIAS, Statistical study of Navier-Slokes equations. Rend. Sem. Mat. Padova, I: 48 (1973) ; II : 49 (1973).
  5. P. L. KAPIŢA, Experiment, teorie, practică, Edit. politică, București, 1978.
  6. D. L. DWOYER, M. Y. HUSSAINI, R. G. VOIGT (eds.), Theoretical approaches to turbulence. Applied Mathematical Sciences, 58, Springer-Verlag, New York, 1985.
  7. G. K. BATGHELOR, An introduction to fluid dynamics, Cambridge University Press, London, 1967.
  8. P. CONSTANTIN, C. FOIAS, O. MANLEY, R. TEMAM, Determining modes and dimension of turbulent flows, J. Fluid Mech., 150 (1985).
  9. H. O. PEITGEN, P. H. RICHTER, The beauly of fractals. Images of complex dynamical systems, Springer-Verlag, Berlin, 1986.
  10. A. LASOTA, M. C. MACKEY, Probabilistic properties of deterministic systems, Cambridge University Press, Cambridge, 1985.
  11. W. HORSTHEMKE, R. LEFEVER, Noise-induced transitions, Springer Series in Synergetics, 15, Springer-Verlag, Berlin, 1984.
  12. M. GHIL (ed.), Turbulence and predictability in geophysical fluid dynamics and climate dynamics, Proceedings of the International School of Physics „Enrico Fermi", cursul LXXXVIII, North-Holland, Amsterdam, 1984.
  13. I. P. CORNFELD, S. V. FOMIN, YA. G. SINAI, Ergodic Theory, seria Grundlehren der mathematischen Wisseschaften, 245 , Springer-Verlag, Berlin, 1982.
  14. K. IKEDA, K. MATSUMOTO, A very high-dimensional chaotic attractor, in „Dynamical systems and nonlinear oscillations (ed. G. Kegami), vol. I, World Scientifié Advanced Series in Dynamical Sytems, 1986, p. 247-253.
  15. E. HOPF, A mathematical example displaying features of turbulence, Comm. Apl. Mat., 1, 4 (1948).
  16. L. LANDAU, On the problem of turbulence, C.R. Acad. Sc. USSR, 4 45 (1944).
  17. D. RUELLE, F. TAKENS, On the nature of turbulence, Comm. Math. Phys., 20 (1971); 23 (1971).
  18. M. J. FEIGNBAUM, Quantitative universalily for a class of nonlinear transformations, J. Statist. Phys., 19 (1978).
  19. B. B. MANDELBROT, Les objets fractals, Flamarion, Paris, 1975.
  20. ADELINA GEORGESCU, Fractalii şi unele aplicații ale lor, St. Cerc. Fiz., 41, 3 (1989);
  21. G. PFISTER, Deterministic chaos in rotational Taylor-Couette flow, LNPh 235, SpringerVerlag, Berlin, 1985, p. 199-210.

  1. *) Lucrare prezentată pe data de 14 III 1989 in cadrul Seminarului de Mecanica Fluidelor condus de Prof. V. N. Constantinescu.
    ST. CERC. MEC. APL., TOM 49, NR. 4, 1990, P. 413-427
  2. a n + 1 a n + 1 a n+1-an+1-an+1a componentă a lui x x vec(x)\vec{x}x.

Related Posts

Turbulence scenarios in determinist chaos

Abstract   AuthorsA. Georgescu C. Vamoş N. SuciuKeywordsReferencesPDFsoon About this paperCite this paper as:JournalSt. Cerc. Mec. Apl. Publisher NameDOINot available…

Global Random Walk simulations of diffusion

AbstractRandom walk methods are suitable to build up convergent solutions for reaction-diffusion problems and were successfully applied to simulations of transport…