| Docente: | Andrea Moiola |
| https://euler.unipv.it/moiola | |
| Email: | andrea.moiola@unipv.it |
| Telefono: | +39 0382 985656 |
| Ufficio: | E15, Dipartimento di Matematica |
| Pagina del corso: | https://euler.unipv.it/moiola/T/MN2024/MN2024.html |
| Pagina ufficiale: | https://unipv.coursecatalogue.cineca.it/insegnamenti/2024/8312/2009/9999/10057 |
| Semestre: | Autunno 2024 |
| Corso di laurea: | Laurea triennale in matematica |
|
Ingegneria computazionale e modellistica per materiali, strutture e tecnologie sostenibili
| |
| Lezioni: | Lunedì 14-16 |
| Mercoledì 14-17 Aula E10 Andremo nel laboratorio informatico del dipartimento di matematica tipicamente il mercoledì dopo le 15. | |
Le registrazioni delle lezioni dell'anno 2020 sono disponibili su Google Drive, il link per accedere si trova sulla pagina Kiro del corso: https://elearning.unipv.it/course/view.php?id=7605 | |
| Ricevimento: | Dopo le lezioni, oppure su appuntamento (accordandosi per email o a lezione) |
| Crediti formativi: | 6 |
| Ore di lezione: | 56, inclusi i laboratori |
| Esame: | Orale + relazione Matlab. |
| La relazione deve essere lunga al massimo 4 pagine ed essere inviata per email in pdf al massimo 24 ore prima dell'appello. | |
| Modello relazione: | .tex, .pdf |
A. Iserles, A First Course in the Numerical Analysis of Differential Equations, Cambridge University Press, 2a edizione, 2009.
R.J. LeVeque, Finite Difference Methods for Ordinary and Partial Differential Equations. Steady-state and Time-dependent Problems, SIAM 2007.
A. Quarteroni, R. Sacco, F. Saleri, P. Gervasio, Matematica Numerica, Springer, 4a edizione, 2014. Capitoli 11-12.
G. Strang, G. Fix, An Analysis of the Finite Element Method, Wellesey–Cambridge press, 2008.
E. Süli, An introduction to the Numerical Analysis of Partial Differential Equations, 2005, dispense disponibili su http://people.maths.ox.ac.uk/suli/nspde.ps.
E. Süli, D.F. Mayers, An introduction to Numerical Analysis, Cambridge University Press, 2003. Capitoli 13-14.
L.N. Trefethen, A. Birkisson, T.A. Driscoll, Exploring ODEs, SIAM, 2018. Pdf e files Matlab disponibili su http://people.maths.ox.ac.uk/trefethen/ExplODE/
A. Tveito, R. Winther, Introduction to Partial Differential Equations. A Computational Approach, Springer 2005.
| Lunedì 30.09.2024 | 2h | E10 | Introduzione al corso. Dettagli pratici.
Problemi ai valori iniziali e problemi ai limiti per equazioni alle derivate ordinarie. Esempio di problema ai limiti lineare ben posto e mal posto a seconda delle condizioni al bordo. Introduzione al metodo di shooting. |
1.1 |
| Mercoledì 2.10.2024 | 1h + 2h | E10 Lab.Inf. |
Il metodo di shooting combinato con il metodo di bisezione. Il metodo di shooting combinato con il metodo di Newton. Esercizio 1.2: uso di ode45 in Matlab. Esercizi 1.3, 1.6, 1.8: implementazione del metodo di shooting con bisezione e Newton. |
1.2 1.2.1 |
| Lunedì 7.10.2024 | 2h | E10 | Differenziazione numerica. Differenze finite in avanti, all'indietro e centrate per l'approssimazione di \(f'(x)\). Errori di troncamento. Differenze finite centrate per la derivata seconda. Interpretazione geometrica e attraverso l'interpolazione polinomiale. L'errore di arrotondamento: cancellazione, crescita dell'errore quando il passo \(h\) è molto piccolo. Stima dell'ordine di convergenza di un metodo numerico con polyfit. |
3.1 3.2 |
| Mercoledì 9.10.2024 | 1h + 2h | E10 Lab.Inf. |
Derivazione informale delle equazioni di continuità, di Fick, del calore, di Poisson, Laplace, e più generali equazioni di diffusione-trasporto-reazione lineari paraboliche. Esercizio 3.6: errore relativo delle differenze finite in avanti e centrate. Errore di troncamento e di roundoff. Stima dell'ordine di convergenza. Esercizi 3.7, 3.8, 3.9: differenze finite del quarto ordine, complex-step, per \(\sin(kx)\). |
2.1.1 2.1.2 2.1.3 |
| Lunedì 14.10.2024 | 2h | E10 |
Equazioni di diffusione-trasporto-reazione lineari: generalizzazioni, equazioni ellittiche. Problemi al bordo lineari in una dimensione. Analogia con problemi algebrici lineari. Esistenza e unicità della soluzione del problema di Dirichlet non omogeneo data l'unicità del problema omogeneo. Dimostrazione dell'unicità della soluzione del problema omogeneo con il metodo dell'energia. |
2.1.4 2.2.3 2.2.1 |
| Mercoledì 16.10.2024 | 2h + 1h | E10 Lab.Inf. |
Problemi al bordo lineari in una dimensione. Principio del massimo e unicità della soluzione del problema di Dirichlet omogeneo. La funzione di Green per \(-u''=f\) con condizioni di Dirichlet. Esercizi 3.6-3.9. |
2.2.2 2.2.4 |
| Lunedì 21.10.2024 | 2h | E10 |
La funzione di Green e la dipendenza continua della soluzione dai dati. Problemi di diffusione-reazione con condizioni al bordo di Neumann: mal posti per \(q=0\), ben posti per \(q>0\). Condizioni al bordo di Robin e periodiche. Il metodo delle differenze finite in una dimensione. Derivazione del metodo per un problema di diffusione-reazione di Dirichlet. Formulazione matriciale. |
2.2.4 2.2.5 4.1 |
| Mercoledì 23.10.2024 | 1h + 2h | E10 Lab.Inf. |
Ripasso: il metodo delle differenze finite in una dimensione. Quali domande ci porremo su questo metodo? Buona positura, convergenza, ordini di convergenza, implementazione efficiente, generalizzazioni. Principio del massimo discreto. Invertibilità della matrice delle differenze finite. Matrici monotone e loro proprietà. Esercizio 4.1: implementazione del metodo delle differenze finite per problemi di Dirichlet. Extra: Esercizio 4.14: convergenza del metodo delle differenze finite. |
4.2 |
| Lunedì 28.10.2024 | 2h | E10 |
Analisi dell'errore per il metodo delle differenze finite. Studio dell'errore di troncamento. Ripasso delle norme vettoriali e matriciali. Studio della stabilità del metodo delle differenze finite: stima della norma dell'inversa della matrice. Stima di convergenza del metodo. Commento sulla convergenza in norme diverse da quella infinito. Discretizzazione del problema di Neumann. Metodo delle differenze finite: approssimazione delle condizioni al bordo con differenze finite in avanti/indietro. |
4.3 |
| Mercoledì 30.10.2024 | 1h + 2h | E10 Lab.Inf. |
Metodo delle differenze finite per il problema di Neumann: approssimazione delle condizioni al bordo con differenze finite in avanti/indietro e con differenze finite centrate. Teorema dei cerchi di Gershgorin; dimostrazione. Matrici a predominanza diagonale e a predominanza diagonale stretta. Per \(q\) positivo, la matrice del metodo delle differenze finite per il problema di Neumann è invertibile, simmetrica, a predominanza diagonale stretta e definita positiva. Stima di stabilità. Esercizio 4.14: convergenza del metodo delle differenze finite. Esercizio 4.19: implementazione di due metodi alle differenze finite per il problema di Neumann, confronto degli ordini di convergenza. |
4.4 |
| Lunedì 4.11.2024 | 2h | E10 |
Implementazione del metodo delle differenze finite: importanza della sparsità. Decomposizione LU di una matrice tridiagonale e soluzione del sistema lineare corrispondente con costo lineare in \(n\). Metodo delle differenze finite per il problema al bordo peridioco. Soluzione di un sistema lineare che è perturbazione di rango uno di un sistema risolvibile in modo efficiente. Applicazione al metodo delle differenze finite per il problema periodico. |
4.5.1 4.5.2 |
| Mercoledì 6.11.2024 | 1h + 2h | E10 Lab.Inf. |
Problemi di diffusione e trasporto. Modello minimo di problema di diffusione e trasporto: \(-\epsilon u''+p u'=0\) in \((0,1)\), \(u(0)=0,u(1)=1\). Due casi limite: diffusione dominante \(0\lt p\ll\epsilon\) (perturbazione regolare) e trasporto dominante \(0\lt\epsilon\ll p\) (perturbazione singolare). Discretizzazione di problemi di diffusione e trasporto con il metodo delle differenze finite centrate. Oscillazioni spurie del metodo delle differenze finite centrate, numero di Péclet locale. Calcolo della soluzione esatta del metodo delle differenze finite centrate per il problema modello; verifica della monotonia della soluzione discreta solo per numero di Péclet locale minore di uno. Esercizio 4.28: metodo delle differenze finite per problemi periodici. Esercizio 4.32: funzione x=TridiagSolver(a,b,c,y) per la risoluzione efficiente di un sistema tridiagonale, test del codice, confronto dei tempi computazionali contro backslash e spdiags. Uso di TridiagSolver con il metodo delle differenze finite. Esercizio 4.36: risoluzione efficiente della perturbazione di rango 1 di un sistema tridiagonale; applicazione al metodo delle differenze finite per problemi periodici. Esercizio 4.40: metodo delle differenze finite centrate per il problema modello di diffusione e trasporto. |
4.6.1 4.6.2 |
| Lunedì 11.11.2024 | 2h | E10 |
Discretizzazione di problemi di diffusione e trasporto con il metodo upwind: definizione attraverso differenze finite in avanti e all'indietro. Viscosità numerica: il metodo upwind coincide con il metodo alle differenze centrate per un problema perturbato, con numero di Péclet locale minore di uno. Autovalori e autofunzioni di \(-u''=\lambda u\) con condizioni di Dirichlet. Approssimazione con il metodo delle differenze finite: autovettori e autovalori discreti. Stima dell'errore commesso dal metodo delle differenze finite. |
4.6.3 4.7 |
| Mercoledì 13.11.2024 | 1h + 2h | E10 Lab.Inf. |
Estensione dell'idea di matrice autoaggiunta (e delle sue proprietà spettrali) agli operatori differenziali: problemi di Sturm-Liouville. Proprietà degli operatori differenziali di Sturm-Liouville, dei loro autovalori e autofunzioni. Perché i problemi differenziali agli autovalori sono interessanti: motivi dalla matematica (espansione in autofunzioni, stabilità), dalla fisica (meccanica quantistica), dall'ingegneria (risonanze). Esercizio 4.47: confronto tra metodo upwind e differenze centrate. Esercizio 4.54: autovalori e autofunzioni di \(-u''=\lambda u\) con condizioni di Dirichlet. Esercizio 4.62: autovalori e autofunzioni di \(-u''+qu=\lambda u\) con potenziali \(q\) a singolo e doppio pozzo. |
4.7.1 |
| Lunedì 18.11.2024 | 2h | E10 |
Differenze finite per problemi non lineari. Il metodo delle differenze finite combinato con il metodo di Newton per l'equazione del pendolo. Differenza tra il residuo \(\|\mathbf G(\mathbf U^k)\|_\infty\) del metodo di Newton e l'errore \(\|\mathbf u-\mathbf U^k\|_\infty\) del metodo delle differenze finite. Aspetti positivi e negativi del metodo delle differenze finite; nelle prossime lezioni cercheremo un metodo che converga più velocemente e la cui incognita sia una funzione, siamo disposti a rinunciare alla sparsità della matrice. Introduzione generale al metodo di collocazione per il problema di Dirichlet. |
4.8 5.1 |
| Mercoledì 20.11.2024 | 3h | E10 Lab.Inf. |
Ripasso: metodo di collocazione per il problema di Dirichlet. Il metodo di collocazione spettrale polinomiale, i polinomi di Legendre integrati. Esercizio 4.66: differenze finite per un problema al bordo per l'equazione (non lineare) del pendolo. Esercizio 4.68: differenze finite per un problema non lineare che coinvolge anche la derivata prima della soluzione. |
5.2 |
| Lunedì 25.11.2024 | 2h | E10 | ||
| Mercoledì 27.11.2024 | 3h | E10 Lab.Inf. |
... ... |
|
| Lunedì 2.12.2024 | 2h | E10 | ||
| Mercoledì 4.12.2024 | 3h | E10 Lab.Inf. |
||
| Mercoledì 11.12.2024 | 3h | E10 Lab.Inf. |
||
| Lunedì 16.12.2024 | 2h | E10 | ||
| Mercoledì 18.12.2024 | 3h | E10 Lab.Inf. |
