Στην παρούσα εργασία παρουσιάζεται μια μέθοδος για την αριθμητική προσομοίωση δισδιάστατης, μη συνεκτικής ροής με ελεύθερη επιφάνεια, που προκύπτει από τη διάδοση κυμάτων βαρύτητας πάνω από πυθμένα με τυχαία μορφολογία. Η μέθοδος βασίζεται στην αριθμητική επίλυση των εξισώσεων Euler, που υπόκεινται σε πλήρως μη γραμμικές οριακές συνθήκες ελεύθερης επιφάνειας και κατάλληλες οριακές συνθήκες πυθμένα, εισόδου και εξόδου, χρησιμοποιώντας ένα υβριδικό σχήμα πεπερασμένων διαφορών και ψευδό-φασματικής μεθόδου. Οι εξισώσεις ροής μετασχηματίζονται έτσι ώστε τα όρια του υπολογιστικού πεδίου να είναι ανεξάρτητα του χρόνου. Η επαλήθευση της μεθόδου επίλυσης γίνεται με την εφαρμογή της στο πρόβλημα της κατανομής θερμοκρασίας σε λεπτή ορθογωνική πλάκα, υπό σταθερές συνθήκες. Για την ελαχιστοποίηση της ανάκλασης χρησιμοποιείται ζώνη απορρόφησης στην περιοχή απορροής. Η αποτελεσματικότητα της ζώνης απορρόφησης τεκμηριώνεται με την παρουσίαση αποτελεσμάτων προσομοίωσης διάδοσης γραμμικών κυματισμών σε πυθμένα σταθερού βάθους. Προκύπτει ότι η ζώνη απορρόφησης που βασίζεται στην επιβολή εξωτερικής δυναμικής πίεσης στην ελεύθερη επιφάνεια εμφανίζει την ελάχιστη ανάκλαση. Αποτελέσματα παρουσιάζονται για την προσομοίωση ροής με ελεύθερη επιφάνεια πάνω από πυθμένα σταθερής κλίσης 1:10 και 1:50, για διαφορετικά μήκη και ύψη κυμάτων εισόδου. Ο μετασχηματισμός των γραμμικών κυμάτων πάνω από την περιοχή σταθερής κλίσης συμφωνεί με τη θεωρία γραμμικής διασποράς για ροή με μικρά ύψη κύματος. Για μη γραμμικούς κυματισμούς, το μήκος κύματος μειώνεται πάνω από την περιοχή σταθερής κλίσης, ενώ η ανύψωση της ελεύθερης επιφάνειας αποκλίνει από την αρχική ημιτονοειδή μορφή και το ύψος κύματος αυξάνει λόγω της ρήχωσης. / A method for the numerical simulation of two-dimensional, inviscid, free-surface flow resulting from the propagation of regular gravity water waves over topography with arbitrary bottom shape is presented. The method is based on the numerical solution of the Euler equations subject to the fully nonlinear free-surface boundary conditions and the appropriate bottom, inflow and outflow conditions using a hybrid scheme of finite-differences and pseudo-spectral method. The formulation includes a boundary-fitted transformation. The validation of the pressure solver is accomplished by applying it to the two-dimensional, temperature, steady problem, under differing boundary conditions. For the free-surface flow, a wave absorption zone is attached at the outflow domain in order to minimize reflection effects. The absorption zone effectiveness is validated by the simulation of linear waves propagation over constant-depth bottom. Minimal reflection occurs when an appropriate external dynamic pressure is imposed on the free surface of the absorption zone. Results are presented for cases of wave propagation over constant slope bottom, with slopes 1:10 and 1:50, for various incoming wavelengths and wave heights. Over the bottom slope, lengths of waves in the linear regime are modified according to linear theory dispersion. For waves in the nonlinear regime, wave lengths are becoming shorter, while the free surface elevation deviates from its initial sinusoidal shape and the wave height increases due to shoaling.