Ghost-point methods for Elliptic and Hyperbolic Equations

In this thesis we have presented a finite-difference ghost-point method to solve elliptic and hyperbolic equations on arbitrary domains. The equations are discretized on a uniform Cartesian grid. At first we applied the Coco-Russo method, which represents a generalization of the finite-difference method for the elliptic equations on arbitrary domains, at the resolution of the Poisson equation. This method proposes a polynomial interpolation technique to impose boundary conditions and therefore the interpolation error can influence the accuracy order of the method itself. We have proposed linear and bilinear interpolation techniques. These conditions are imposed on the projections of the ghost points on the border of the domain. The numerical tests performed on the behaviors of the inverse matrix of the method, of the error and of the consistency error confirm the stability and convergence of the Coco-Russo method in 1D, 2D and 3D, in the case of Dirichlet problems and in the case of mixed problems. We have also presented a rigorous proof of the stability and convergence of the numerical method in the one-dimensional case. Once we certain of the convergence and stability of the Coco-Russo method, our interest it has moved to the study of the Euler equations of the gas dynamic. The Coco-Russo method was applied for the development of a semi-implicit method for Euler equations on domains that have obstacles, to impose boundary conditions in a manner similar to elliptic equations. This method being semi-implicit overcomes the problem of spatial restriction to guarantee the stability of the method typical of explicit methods.

In questa tesi abbiamo presentato un metodo ghost-point alle differenze finite per la risoluzione numerica delle equazioni ellittiche e iperboliche su domini arbitrari. Le equazioni sono discretizzate su una griglia Cartesiana regolare. In un primo momento abbiamo applicato il metodo Coco-Russo, che consiste in un metodo alle differenze finite applicato al caso di domini arbitrari, alla risoluzione dell'equazione di Poisson. Questo metodo propone una tecnica di interpolazione polinomiale per imporre le condizioni al contorno pertanto l'errore di interpolazione può influenzare l'ordine di accuratezza del metodo stesso. In particolare, noi abbiamo proposto tecniche di interpolazione lineare e bilineare. Tali condizioni sono imposte sulle proiezioni dei punti fantasma sul bordo del dominio. I test numerici eseguiti sui comportamenti della matrice inversa del metodo, dell'errore e dell'errore di consistenza confermano la stabilità e la convergenza del metodo Coco-Russo in 1D, 2D e 3D, sia nel caso di problemi di Dirichlet che nel caso di problemi misti. Abbiamo anche presentato una dimostrazione rigorosa della stabilità e della convergenza del metodo numerico nel caso monodimensionale. Una volta certi della convergenza e della stabilità del metodo Coco-Russo, il nostro interesse si è spostato allo studio delle equazioni di Eulero della gasdinamica. Il metodo Coco-Russo è stato quindi applicato per lo sviluppo di un metodo semi-implicito per equazioni di Eulero su domini che presentano ostacoli, per imporre le condizioni al contorno in modo simile a quanto fatto nelle equazioni ellittiche. Questo metodo essendo semi-implicito supera il problema della restrizione spaziale per garantire la stabilità del metodo tipico dei metodi espliciti.