03 02 15

Page 1

Metodi Analitici e Numerici per l'Ingegneria II prova parziale 3 F ebbraio 2015 DOMANDA 1 Si consideri il seguente problema:    −u′′(x) + u′(x) = x; 0 < x < 1 u(0) = −1 u′(1) = 1 : 1. (T) Si ricavi la formulazione debole del problema giustificando i passaggi e motivando la scelta degli spazi funzionali. 2. (T) Si dimostri che la forma bilineare ottenuta al punto precedente ` e continua e coerciva. 3. (T) Si enunci il lemma (o teorema) di Lax-Milgram. 4. (E) Si disegnino le funzioni di base dello spazio di elementi finiti lineari definiti su una par- tizione unifome dell’intervallo [0; 1] di passo h = 1=4, che approssimano lo spazio funzionale scelto per la formulazione debole. 5. (E) Si scriva il sistema lineare corrispondente all’approssimazione ad elementi finiti. 6. (M) Si implementi in MATLAB il sistema, lo si risolva e si disegni il grafico della soluzione. DOMANDA 2 1. (T) Si determini una soluzione u = u(x; t) del problema    ut(x; t) − uxx(x; t) = 0 in (0 ; ) × (0; +∞); u(x; 0) = − 4 2 x2 + 2  x + 1 in (0 ; ); u(0; t) = 1 e u(; t) = −1; in (0 ; +∞): senza riportare tutta la separazione di variabili, ma solo i passaggi principali. In particolare esprimere i coefficienti sotto forma di integrale senza calcolarli esplicitamente. 2. (E) Qual ` e la soluzione stazionaria? 3. (T) La soluzione u ` e continua su [0 ; ] × [0; +∞)? 4. (T) Enunciare il principio del massimo per u(x; t). 5. (T) Dimostrare l’unicit` a della soluzione per il problema. 6. (M) Risolvere l’equazione con MATLAB (usando DCT o DST per calcolare i coefficienti) e disegnare il grafico della soluzione per t = 1. (Si ricorda che per una funzione f (x), utilizzando la formula del punto medio, ak ≃ 2 N N∑ j=1 f ( (2j − 1) 2N ) cos ( k(2j − 1) 2N ) e utilizzando la formula del trapezio bk ≃ 2…