Musica e Origini I Parte; presentazione didattica. Prof.ssa Serra Manuela
Presentazione parte cecio
-
Upload
andrea-ruggiu -
Category
Entertainment & Humor
-
view
35 -
download
0
description
Transcript of Presentazione parte cecio
![Page 1: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/1.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Il problema Phase Field Crystal bidimensionale
Cesare Molinari 782407Andrea Alessandro Ruggiu 782766
Politecnico di Milano
16/01/2014
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 1 / 20
![Page 2: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/2.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Sommario
1 Discretizzazione e Linearizzazione
2 Risultati Numerici
3 Ulteriori Approcci
4 Sviluppi Futuri
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 2 / 20
![Page 3: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/3.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Formulazione Debole: Metodo di Splitting
Obiettivo: utilizzare Elementi Finiti P1 in spazioMetodo: Splitting di Equazione in Sistema equivalente
∂ρ∂t = ∆σ
σ = ϕ (ρ) + 2Dk2ω + D∆ω
ω = ∆ρ
Formulazione Debole: dato ρ0 ∈ L2 (Ω) vogliamo trovareρ, σ, ω ∈ L2 (0,T ;H1
p (Ω))t.c. ∀u, v , z ∈ H1
p (Ω)(∂ρ∂t , u
)0
+ (∇σ,∇u)0 = 0,
(σ, v)0 −(2Dk2ω, v
)0 + (D∇ω,∇v)0 − (ϕ (ρ) , v)0 = 0,
(ω, z)0 + (∇ρ,∇z)0 = 0,
con ρ (t = 0) = ρ0Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 3 / 20
![Page 4: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/4.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Formulazione Debole: Metodo di Splitting
Obiettivo: utilizzare Elementi Finiti P1 in spazioMetodo: Splitting di Equazione in Sistema equivalente
∂ρ∂t = ∆σ
σ = ϕ (ρ) + 2Dk2ω + D∆ω
ω = ∆ρ
Formulazione Debole: dato ρ0 ∈ L2 (Ω) vogliamo trovareρ, σ, ω ∈ L2 (0,T ;H1
p (Ω))t.c. ∀u, v , z ∈ H1
p (Ω)(∂ρ∂t , u
)0
+ (∇σ,∇u)0 = 0,
(σ, v)0 −(2Dk2ω, v
)0 + (D∇ω,∇v)0 − (ϕ (ρ) , v)0 = 0,
(ω, z)0 + (∇ρ,∇z)0 = 0,
con ρ (t = 0) = ρ0Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 3 / 20
![Page 5: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/5.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Formulazione Debole: Metodo di Splitting
Obiettivo: utilizzare Elementi Finiti P1 in spazioMetodo: Splitting di Equazione in Sistema equivalente
∂ρ∂t = ∆σ
σ = ϕ (ρ) + 2Dk2ω + D∆ω
ω = ∆ρ
Formulazione Debole: dato ρ0 ∈ L2 (Ω) vogliamo trovareρ, σ, ω ∈ L2 (0,T ;H1
p (Ω))t.c. ∀u, v , z ∈ H1
p (Ω)(∂ρ∂t , u
)0
+ (∇σ,∇u)0 = 0,
(σ, v)0 −(2Dk2ω, v
)0 + (D∇ω,∇v)0 − (ϕ (ρ) , v)0 = 0,
(ω, z)0 + (∇ρ,∇z)0 = 0,
con ρ (t = 0) = ρ0Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 3 / 20
![Page 6: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/6.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Formulazione Debole: Metodo di Splitting
Obiettivo: utilizzare Elementi Finiti P1 in spazioMetodo: Splitting di Equazione in Sistema equivalente
∂ρ∂t = ∆σ
σ = ϕ (ρ) + 2Dk2ω + D∆ω
ω = ∆ρ
Formulazione Debole: dato ρ0 ∈ L2 (Ω) vogliamo trovareρ, σ, ω ∈ L2 (0,T ;H1
p (Ω))t.c. ∀u, v , z ∈ H1
p (Ω)(∂ρ∂t , u
)0
+ (∇σ,∇u)0 = 0,
(σ, v)0 −(2Dk2ω, v
)0 + (D∇ω,∇v)0 − (ϕ (ρ) , v)0 = 0,
(ω, z)0 + (∇ρ,∇z)0 = 0,
con ρ (t = 0) = ρ0Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 3 / 20
![Page 7: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/7.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Formulazione Debole: Metodo di Splitting
Obiettivo: utilizzare Elementi Finiti P1 in spazioMetodo: Splitting di Equazione in Sistema equivalente
∂ρ∂t = ∆σ
σ = ϕ (ρ) + 2Dk2ω + D∆ω
ω = ∆ρ
Formulazione Debole: dato ρ0 ∈ L2 (Ω) vogliamo trovareρ, σ, ω ∈ L2 (0,T ;H1
p (Ω))t.c. ∀u, v , z ∈ H1
p (Ω)(∂ρ∂t , u
)0
+ (∇σ,∇u)0 = 0,
(σ, v)0 −(2Dk2ω, v
)0 + (D∇ω,∇v)0 − (ϕ (ρ) , v)0 = 0,
(ω, z)0 + (∇ρ,∇z)0 = 0,
con ρ (t = 0) = ρ0Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 3 / 20
![Page 8: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/8.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Formulazione Debole: Metodo di Splitting
Obiettivo: utilizzare Elementi Finiti P1 in spazioMetodo: Splitting di Equazione in Sistema equivalente
∂ρ∂t = ∆σ
σ = ϕ (ρ) + 2Dk2ω + D∆ω
ω = ∆ρ
Formulazione Debole: dato ρ0 ∈ L2 (Ω) vogliamo trovareρ, σ, ω ∈ L2 (0,T ;H1
p (Ω))t.c. ∀u, v , z ∈ H1
p (Ω)(∂ρ∂t , u
)0
+ (∇σ,∇u)0 = 0,
(σ, v)0 −(2Dk2ω, v
)0 + (D∇ω,∇v)0 − (ϕ (ρ) , v)0 = 0,
(ω, z)0 + (∇ρ,∇z)0 = 0,
con ρ (t = 0) = ρ0Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 3 / 20
![Page 9: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/9.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Spazi di Semidiscretizzazione Elementi Finiti
Considerando una Triangolazione Th =⋃Nh
j=1 Kj di Nh elementisu Ω, definiamo
Vh =vh ∈ C 0 (Ω) : vh |Kj
∈ P1, ∀Kj ∈ Th
e
Vh =vh ∈ Vh : vh(x , y) periodica di periodo Ω
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 4 / 20
![Page 10: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/10.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Spazi di Semidiscretizzazione Elementi Finiti
Considerando una Triangolazione Th =⋃Nh
j=1 Kj di Nh elementisu Ω, definiamo
Vh =vh ∈ C 0 (Ω) : vh |Kj
∈ P1, ∀Kj ∈ Th
e
Vh =vh ∈ Vh : vh(x , y) periodica di periodo Ω
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 4 / 20
![Page 11: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/11.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Semidiscretizzazione in Spazio
Dato ρh0 ∈ Vh, ∀t > 0 fissato trovareρh (t) , σh (t) , ωh (t) ∈ Vh t.c. ∀uh, vh, zh ∈ Vh(∂ρh∂t , uh
)0
+ (∇σh,∇uh)0 = 0
(σh, vh)0 −(2Dk2ωh, vh
)0 + (D∇ωh,∇vh)0 − (ϕ (ρh) , vh)0 = 0
(ωh, zh)0 + (∇ρh,∇zh)0 = 0
con ρh (t = 0) = ρh0
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 5 / 20
![Page 12: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/12.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Semidiscretizzazione in Spazio
Dato ρh0 ∈ Vh, ∀t > 0 fissato trovareρh (t) , σh (t) , ωh (t) ∈ Vh t.c. ∀uh, vh, zh ∈ Vh(∂ρh∂t , uh
)0
+ (∇σh,∇uh)0 = 0
(σh, vh)0 −(2Dk2ωh, vh
)0 + (D∇ωh,∇vh)0 − (ϕ (ρh) , vh)0 = 0
(ωh, zh)0 + (∇ρh,∇zh)0 = 0
con ρh (t = 0) = ρh0
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 5 / 20
![Page 13: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/13.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Semidiscretizzazione in Spazio
Dato ρh0 ∈ Vh, ∀t > 0 fissato trovareρh (t) , σh (t) , ωh (t) ∈ Vh t.c. ∀uh, vh, zh ∈ Vh(∂ρh∂t , uh
)0
+ (∇σh,∇uh)0 = 0
(σh, vh)0 −(2Dk2ωh, vh
)0 + (D∇ωh,∇vh)0 − (ϕ (ρh) , vh)0 = 0
(ωh, zh)0 + (∇ρh,∇zh)0 = 0
con ρh (t = 0) = ρh0
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 5 / 20
![Page 14: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/14.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Discretizzazione Completa con Eulero all’Indietro inTempo
Dato ρh0 ∈ Vh, ∀n = 0, 1, . . . fissato trovareρn+1h , σn+1
h , ωn+1h ∈ Vh t.c. ∀uh, vh, zh ∈ Vh
(ρn+1h −ρnh
∆tn, uh
)0
+ (∇σh,∇uh)0 = 0
(σh, vh)0 −(2Dk2ωh, vh
)0 + (D∇ωh,∇vh)0 +
−(ϕ(ρn+1h
), vh)0 = 0
(ωh, zh)0 +(∇ρn+1
h ,∇zh)0 = 0
con ρ0h = ρh0
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 6 / 20
![Page 15: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/15.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Discretizzazione Completa con Eulero all’Indietro inTempo
Dato ρh0 ∈ Vh, ∀n = 0, 1, . . . fissato trovareρn+1h , σn+1
h , ωn+1h ∈ Vh t.c. ∀uh, vh, zh ∈ Vh
(ρn+1h −ρnh
∆tn, uh
)0
+ (∇σh,∇uh)0 = 0
(σh, vh)0 −(2Dk2ωh, vh
)0 + (D∇ωh,∇vh)0 +
−(ϕ(ρn+1h
), vh)0 = 0
(ωh, zh)0 +(∇ρn+1
h ,∇zh)0 = 0
con ρ0h = ρh0
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 6 / 20
![Page 16: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/16.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Discretizzazione Completa con Eulero all’Indietro inTempo
Dato ρh0 ∈ Vh, ∀n = 0, 1, . . . fissato trovareρn+1h , σn+1
h , ωn+1h ∈ Vh t.c. ∀uh, vh, zh ∈ Vh
(ρn+1h −ρnh
∆tn, uh
)0
+ (∇σh,∇uh)0 = 0
(σh, vh)0 −(2Dk2ωh, vh
)0 + (D∇ωh,∇vh)0 +
−(ϕ(ρn+1h
), vh)0 = 0
(ωh, zh)0 +(∇ρn+1
h ,∇zh)0 = 0
con ρ0h = ρh0
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 6 / 20
![Page 17: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/17.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Linearizzazione: Algoritmo Iterativo di Newton
Inizializzazione: ρn,0h = ρn−1h
Problema linearizzato: per n = 0, 1, . . . fissata, ∀k = 0, 1, . . .trovare ρn,k+1
h , σn,k+1h , ωn,k+1
h ∈ Vh t.c. ∀uh, vh, zh ∈ Vh
(ρn,k+1h −ρn,0h
∆tn, uh
)0
+ (∇σh,∇uh)0 = 0
(σh, vh)0 −(2Dk2ωh, vh
)0 + (D∇ωh,∇vh)0 −
(ϕ(ρn,kh
), vh
)0
−(ϕ′(ρn,kh
)ρn,k+1h , vh
)0
+(ϕ′(ρn,kh
)ρn,kh , vh
)0
= 0
(ωh, zh)0 +(∇ρn,k+1
h ,∇zh)
0= 0
con ρ0h = ρh0
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 7 / 20
![Page 18: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/18.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Linearizzazione: Algoritmo Iterativo di Newton
Inizializzazione: ρn,0h = ρn−1h
Problema linearizzato: per n = 0, 1, . . . fissata, ∀k = 0, 1, . . .trovare ρn,k+1
h , σn,k+1h , ωn,k+1
h ∈ Vh t.c. ∀uh, vh, zh ∈ Vh
(ρn,k+1h −ρn,0h
∆tn, uh
)0
+ (∇σh,∇uh)0 = 0
(σh, vh)0 −(2Dk2ωh, vh
)0 + (D∇ωh,∇vh)0 −
(ϕ(ρn,kh
), vh
)0
−(ϕ′(ρn,kh
)ρn,k+1h , vh
)0
+(ϕ′(ρn,kh
)ρn,kh , vh
)0
= 0
(ωh, zh)0 +(∇ρn,k+1
h ,∇zh)
0= 0
con ρ0h = ρh0
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 7 / 20
![Page 19: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/19.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Linearizzazione: Algoritmo Iterativo di Newton
Inizializzazione: ρn,0h = ρn−1h
Problema linearizzato: per n = 0, 1, . . . fissata, ∀k = 0, 1, . . .trovare ρn,k+1
h , σn,k+1h , ωn,k+1
h ∈ Vh t.c. ∀uh, vh, zh ∈ Vh
(ρn,k+1h −ρn,0h
∆tn, uh
)0
+ (∇σh,∇uh)0 = 0
(σh, vh)0 −(2Dk2ωh, vh
)0 + (D∇ωh,∇vh)0 −
(ϕ(ρn,kh
), vh
)0
−(ϕ′(ρn,kh
)ρn,k+1h , vh
)0
+(ϕ′(ρn,kh
)ρn,kh , vh
)0
= 0
(ωh, zh)0 +(∇ρn,k+1
h ,∇zh)
0= 0
con ρ0h = ρh0
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 7 / 20
![Page 20: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/20.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Linearizzazione: Algoritmo Iterativo di Newton
Inizializzazione: ρn,0h = ρn−1h
Problema linearizzato: per n = 0, 1, . . . fissata, ∀k = 0, 1, . . .trovare ρn,k+1
h , σn,k+1h , ωn,k+1
h ∈ Vh t.c. ∀uh, vh, zh ∈ Vh
(ρn,k+1h −ρn,0h
∆tn, uh
)0
+ (∇σh,∇uh)0 = 0
(σh, vh)0 −(2Dk2ωh, vh
)0 + (D∇ωh,∇vh)0 −
(ϕ(ρn,kh
), vh
)0
−(ϕ′(ρn,kh
)ρn,k+1h , vh
)0
+(ϕ′(ρn,kh
)ρn,kh , vh
)0
= 0
(ωh, zh)0 +(∇ρn,k+1
h ,∇zh)
0= 0
con ρ0h = ρh0
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 7 / 20
![Page 21: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/21.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Test di Arresto
I kn < toll (nelle Simulazioni, toll = 10−3)
dove I kn è l’Incremento, definito da
I kn =
√‖dk
n ρh‖2H1(Ω) + ‖dk
n σh‖2H1(Ω) + ‖dk
nωh‖2H1(Ω)
|Ω|.
con
dkn (·) = (·)n,k+1 − (·)n,k
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 8 / 20
![Page 22: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/22.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Test di Arresto
I kn < toll (nelle Simulazioni, toll = 10−3)
dove I kn è l’Incremento, definito da
I kn =
√‖dk
n ρh‖2H1(Ω) + ‖dk
n σh‖2H1(Ω) + ‖dk
nωh‖2H1(Ω)
|Ω|.
con
dkn (·) = (·)n,k+1 − (·)n,k
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 8 / 20
![Page 23: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/23.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Diagramma di Fase
Obiettivo: riproduzione del Diagramma di Fase (ricavatoempiricamente in [EKHG])
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 9 / 20
![Page 24: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/24.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Validazione Numerica del Diagramma di Fase
Metodo:
fissato il valore di ρ = 0.07, Simulazioni per differenti ε;fissato il valore di ε = 0.15, Simulazioni per differenti ρ
con Domino Ω = [0, 370]2 e Dato Iniziale:
ρ0 =
ρ+ C
(cos (qX ) cos
(qY√
3
)− 1
2 cos(
2qY√3
))se (X ,Y ) ∈ [220, 250]× [−70,−20] ,
ρ altrimenti
dove C = 0.446, q = 0.66 e (X ,Y ) è un sistema di coordinateruotato di −π
3Griglia e Timestep: griglia non strutturata con 170 nodi perlato; passo temporale ∆t = 4
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 10 / 20
![Page 25: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/25.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Validazione Numerica del Diagramma di Fase
Metodo:
fissato il valore di ρ = 0.07, Simulazioni per differenti ε;fissato il valore di ε = 0.15, Simulazioni per differenti ρ
con Domino Ω = [0, 370]2 e Dato Iniziale:
ρ0 =
ρ+ C
(cos (qX ) cos
(qY√
3
)− 1
2 cos(
2qY√3
))se (X ,Y ) ∈ [220, 250]× [−70,−20] ,
ρ altrimenti
dove C = 0.446, q = 0.66 e (X ,Y ) è un sistema di coordinateruotato di −π
3Griglia e Timestep: griglia non strutturata con 170 nodi perlato; passo temporale ∆t = 4
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 10 / 20
![Page 26: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/26.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Validazione Numerica del Diagramma di Fase
Metodo:
fissato il valore di ρ = 0.07, Simulazioni per differenti ε;fissato il valore di ε = 0.15, Simulazioni per differenti ρ
con Domino Ω = [0, 370]2 e Dato Iniziale:
ρ0 =
ρ+ C
(cos (qX ) cos
(qY√
3
)− 1
2 cos(
2qY√3
))se (X ,Y ) ∈ [220, 250]× [−70,−20] ,
ρ altrimenti
dove C = 0.446, q = 0.66 e (X ,Y ) è un sistema di coordinateruotato di −π
3Griglia e Timestep: griglia non strutturata con 170 nodi perlato; passo temporale ∆t = 4
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 10 / 20
![Page 27: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/27.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Validazione Numerica del Diagramma di Fase
Metodo:
fissato il valore di ρ = 0.07, Simulazioni per differenti ε;fissato il valore di ε = 0.15, Simulazioni per differenti ρ
con Domino Ω = [0, 370]2 e Dato Iniziale:
ρ0 =
ρ+ C
(cos (qX ) cos
(qY√
3
)− 1
2 cos(
2qY√3
))se (X ,Y ) ∈ [220, 250]× [−70,−20] ,
ρ altrimenti
dove C = 0.446, q = 0.66 e (X ,Y ) è un sistema di coordinateruotato di −π
3Griglia e Timestep: griglia non strutturata con 170 nodi perlato; passo temporale ∆t = 4
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 10 / 20
![Page 28: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/28.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Validazione Numerica del Diagramma di Fase
Metodo:
fissato il valore di ρ = 0.07, Simulazioni per differenti ε;fissato il valore di ε = 0.15, Simulazioni per differenti ρ
con Domino Ω = [0, 370]2 e Dato Iniziale:
ρ0 =
ρ+ C
(cos (qX ) cos
(qY√
3
)− 1
2 cos(
2qY√3
))se (X ,Y ) ∈ [220, 250]× [−70,−20] ,
ρ altrimenti
dove C = 0.446, q = 0.66 e (X ,Y ) è un sistema di coordinateruotato di −π
3Griglia e Timestep: griglia non strutturata con 170 nodi perlato; passo temporale ∆t = 4
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 10 / 20
![Page 29: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/29.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Diagramma di Fase per ε = 0.15 fissato al variare diρ
ρ = 0.3: Regime a Stato Stazionario Costante (Fase Liquida)
C:/Andrea/Tabella030.jpg
VIDEO
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 11 / 20
![Page 30: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/30.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Diagramma di Fase per ε = 0.15 fissato al variare diρ
ρ = 0.15: Regime a Stato Stazionario Esagonale (Fase SolidaCristallina)
C:/Andrea/Tabella015.jpg
VIDEOMolinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 12 / 20
![Page 31: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/31.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Diagramma di Fase per ε = 0.15 fissato al variare diρ
ρ = 0.03: Regime a Stato Stazionario a Striscie
C:/Andrea/Tabella003.jpg
VIDEO
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 13 / 20
![Page 32: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/32.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Crescita di Cristalli differentemente orientati
Valori delle Costanti: ρ = 0.285 e ε = 0.25 (RegimeEsagonale)Dato iniziale:
ρ0 =
ρ+ C(cos (qx) cos
(qy√
3
)− 1
2 cos(
2qy√3
))se (x , y) ∈ [140, 170]× [220, 260]
ρ+ C(cos (qX1) cos
(qY1√
3
)− 1
2 cos(
2qY1√3
))se (X1,Y1) ∈ [160, 180]× [−30, 20]
ρ+ C(cos (qX2) cos
(qY2√
3
)− 1
2 cos(
2qY2√3
))se (X2,Y2) ∈ [70, 100]× [280, 310]
ρ altrimenti
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 14 / 20
![Page 33: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/33.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Crescita di Cristalli differentemente orientati
Valori delle Costanti: ρ = 0.285 e ε = 0.25 (RegimeEsagonale)Dato iniziale:
ρ0 =
ρ+ C(cos (qx) cos
(qy√
3
)− 1
2 cos(
2qy√3
))se (x , y) ∈ [140, 170]× [220, 260]
ρ+ C(cos (qX1) cos
(qY1√
3
)− 1
2 cos(
2qY1√3
))se (X1,Y1) ∈ [160, 180]× [−30, 20]
ρ+ C(cos (qX2) cos
(qY2√
3
)− 1
2 cos(
2qY2√3
))se (X2,Y2) ∈ [70, 100]× [280, 310]
ρ altrimenti
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 14 / 20
![Page 34: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/34.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Crescita di Cristalli differentemente orientati
VIDEOIn Regime Stazionario, evidenza di Lacune (ovvero zone diconfine in cui avviene brusco il cambio di struttura cristallina)
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 15 / 20
![Page 35: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/35.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Crescita di Cristalli differentemente orientati
VIDEOIn Regime Stazionario, evidenza di Lacune (ovvero zone diconfine in cui avviene brusco il cambio di struttura cristallina)
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 15 / 20
![Page 36: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/36.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Test di Convergenza in Tempo
Settings:Dominio: Ω = [0, 370]2
Dato Iniziale: come in Esperimento di ValidazioneDiagramma di FaseValori Costanti: D = k = 1, ε = 0.14 e ρ = 0.07Griglia: non strutturata di 170 nodi per lato
Errore Relativo in H1 di ciascuna Soluzione ρh,∆t rispetto aquella con il passo di tempo più fine ρ∗h al Tempo FinaleT = 128 per diversi Passi di Tempo ∆t:
err (∆t) =1√|Ω|‖ρh,∆t − ρ∗h‖H1(Ω)
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 16 / 20
![Page 37: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/37.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Test di Convergenza in Tempo
Settings:Dominio: Ω = [0, 370]2
Dato Iniziale: come in Esperimento di ValidazioneDiagramma di FaseValori Costanti: D = k = 1, ε = 0.14 e ρ = 0.07Griglia: non strutturata di 170 nodi per lato
Errore Relativo in H1 di ciascuna Soluzione ρh,∆t rispetto aquella con il passo di tempo più fine ρ∗h al Tempo FinaleT = 128 per diversi Passi di Tempo ∆t:
err (∆t) =1√|Ω|‖ρh,∆t − ρ∗h‖H1(Ω)
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 16 / 20
![Page 38: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/38.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Test di Convergenza in Tempo
C:/Andrea/Tabel.jpg
0 0.5 1 1.5 2 2.5 3 3.5 40
0.05
0.1
0.15
0.2
0.25
Decrescita Lineare dell’Errore Relativo: è prova di Convergenzain Tempo Lineare
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 17 / 20
![Page 39: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/39.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Test di Convergenza in Tempo
C:/Andrea/Tabel.jpg
0 0.5 1 1.5 2 2.5 3 3.5 40
0.05
0.1
0.15
0.2
0.25
Decrescita Lineare dell’Errore Relativo: è prova di Convergenzain Tempo Lineare
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 17 / 20
![Page 40: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/40.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Approcci della Decomposizione dell’Energia
Problema: il Funzionale Energia Libera non è Convesso;Possibile approccio (proposto in [WWL] e seguito poi in[EW]): Decomposizione di Funzionale Energia Libera comedifferenza di due Funzionali Convessi
Infatti, perEc [ρ] =
∫Ω
[12 |∆ρ|
2 + 1−ε2 ρ2 + 1
4ρ4]dx ;
Ee [ρ] =∫
Ω |∇ρ|2 dx ,
si ha cheE = Ec − Ee
Risultato Teorico: Stabilità della Discretizzazione in TempoSemi-Implicita che deriva da
ρn+1 − ρn
∆t= ∆
(δEc
(ρn+1)δρ
)−∆
(δEe (ρn)
δρ
).
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 18 / 20
![Page 41: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/41.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Approcci della Decomposizione dell’Energia
Problema: il Funzionale Energia Libera non è Convesso;Possibile approccio (proposto in [WWL] e seguito poi in[EW]): Decomposizione di Funzionale Energia Libera comedifferenza di due Funzionali Convessi
Infatti, perEc [ρ] =
∫Ω
[12 |∆ρ|
2 + 1−ε2 ρ2 + 1
4ρ4]dx ;
Ee [ρ] =∫
Ω |∇ρ|2 dx ,
si ha cheE = Ec − Ee
Risultato Teorico: Stabilità della Discretizzazione in TempoSemi-Implicita che deriva da
ρn+1 − ρn
∆t= ∆
(δEc
(ρn+1)δρ
)−∆
(δEe (ρn)
δρ
).
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 18 / 20
![Page 42: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/42.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Approcci della Decomposizione dell’Energia
Problema: il Funzionale Energia Libera non è Convesso;Possibile approccio (proposto in [WWL] e seguito poi in[EW]): Decomposizione di Funzionale Energia Libera comedifferenza di due Funzionali Convessi
Infatti, perEc [ρ] =
∫Ω
[12 |∆ρ|
2 + 1−ε2 ρ2 + 1
4ρ4]dx ;
Ee [ρ] =∫
Ω |∇ρ|2 dx ,
si ha cheE = Ec − Ee
Risultato Teorico: Stabilità della Discretizzazione in TempoSemi-Implicita che deriva da
ρn+1 − ρn
∆t= ∆
(δEc
(ρn+1)δρ
)−∆
(δEe (ρn)
δρ
).
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 18 / 20
![Page 43: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/43.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Approcci della Decomposizione dell’Energia
Problema: il Funzionale Energia Libera non è Convesso;Possibile approccio (proposto in [WWL] e seguito poi in[EW]): Decomposizione di Funzionale Energia Libera comedifferenza di due Funzionali Convessi
Infatti, perEc [ρ] =
∫Ω
[12 |∆ρ|
2 + 1−ε2 ρ2 + 1
4ρ4]dx ;
Ee [ρ] =∫
Ω |∇ρ|2 dx ,
si ha cheE = Ec − Ee
Risultato Teorico: Stabilità della Discretizzazione in TempoSemi-Implicita che deriva da
ρn+1 − ρn
∆t= ∆
(δEc
(ρn+1)δρ
)−∆
(δEe (ρn)
δρ
).
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 18 / 20
![Page 44: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/44.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Approcci della Decomposizione dell’Energia
Problema: il Funzionale Energia Libera non è Convesso;Possibile approccio (proposto in [WWL] e seguito poi in[EW]): Decomposizione di Funzionale Energia Libera comedifferenza di due Funzionali Convessi
Infatti, perEc [ρ] =
∫Ω
[12 |∆ρ|
2 + 1−ε2 ρ2 + 1
4ρ4]dx ;
Ee [ρ] =∫
Ω |∇ρ|2 dx ,
si ha cheE = Ec − Ee
Risultato Teorico: Stabilità della Discretizzazione in TempoSemi-Implicita che deriva da
ρn+1 − ρn
∆t= ∆
(δEc
(ρn+1)δρ
)−∆
(δEe (ρn)
δρ
).
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 18 / 20
![Page 45: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/45.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Approcci della Decomposizione dell’Energia
Accenno all’approccio in [EW]: aggiungere e togliere untermine a entrambe le Componenti dell’Energia senzamodificarne la Convessità e in modo tale da ricondurre laNon-Linearità tutta nella Componente calcolata in esplicitoall’istante temporale precedente
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 19 / 20
![Page 46: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/46.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Sviluppi Futuri
Visto che i Lavori riguardo alla Decomposizionedell’Energia sono tutti in Differenze Finite, si potrebbeimplementare questa idea in Schemi a Elementi Finiti;Splitting a due equazioni tramite Spazi a Elementi FinitiNURBS (globalmente C1), che permettono il calcolodell’Energia Libera;Simulazioni Numeriche dell’Equazione PFC in TreDimensioni;
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 20 / 20
![Page 47: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/47.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Sviluppi Futuri
Visto che i Lavori riguardo alla Decomposizionedell’Energia sono tutti in Differenze Finite, si potrebbeimplementare questa idea in Schemi a Elementi Finiti;Splitting a due equazioni tramite Spazi a Elementi FinitiNURBS (globalmente C1), che permettono il calcolodell’Energia Libera;Simulazioni Numeriche dell’Equazione PFC in TreDimensioni;
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 20 / 20
![Page 48: Presentazione parte cecio](https://reader033.fdocumenti.com/reader033/viewer/2022051818/549798d5ac7959292e8b53af/html5/thumbnails/48.jpg)
ProblemaPFC 2D
Molinari -Ruggiu
Discretizzazionee Linearizza-zione
RisultatiNumerici
UlterioriApprocci
SviluppiFuturi
Sviluppi Futuri
Visto che i Lavori riguardo alla Decomposizionedell’Energia sono tutti in Differenze Finite, si potrebbeimplementare questa idea in Schemi a Elementi Finiti;Splitting a due equazioni tramite Spazi a Elementi FinitiNURBS (globalmente C1), che permettono il calcolodell’Energia Libera;Simulazioni Numeriche dell’Equazione PFC in TreDimensioni;
Molinari - Ruggiu (Politecnico di Milano) Problema PFC 2D 16/01/2014 20 / 20