Για την επίλυση προβλημάτων πολλαπλών πεδίων / πολλαπλών φυσικών έχουν προταθεί διάφορες τεχνικές στη βιβλιογραφία. Οι μέθοδοι χαλάρωσης στις διεπαφές είναι μια ενδιαφέρουσα προσέγγιση για την επίλυση αυτών των προβλημάτων. Υποθέτοντας κάποιες αρχικές τιμές στις διεπαφές του προβλήματος, οι μέθοδοι χαλάρωσης στις διεπαφές επιλύουν επαναληπτικά τα υποπροβλήματα και χαλαρώνουν τις τιμές στις διεπαφές έως ότου επιτευχθεί σύγκλιση. Τα κυριότερα πλεονεκτήματα των μεθόδων αυτών είναι ότι οι ρυθμοί σύγκλισής τους εξαρτώνται μόνο από τις παραμέτρους των ίδιων των προβλημάτων, τις παραμέτρους που σχετίζονται με την ανάλυση των προβλημάτων σε υποπροβλήματα και τις παραμέτρους που σχετίζονται με τον τελεστή που εφαρμόζεται στις διεπαφές. Στην παρούσα μεταπτυχιακή εργασία παρουσιάζεται μια νέα υλοποίηση μιας μεθόδου χαλάρωσης στις διεπαφές, της μεθόδου GEO. Η GEO βασίζεται σε έναν απλό μηχανισμό γεωμετρικής διόρθωσης και δρα επαναληπτικά, ώστε να χαλαρώσει τις τιμές στις διεπαφές. Πιο συγκεκριμένα, προσθέτει στις παλιές τιμές των διεπαφών έναν γεωμετρικά σταθμισμένο συνδυασμό των κανονικών παραγώγων των οριακών σημείων των γειτονικών πεδίων. Η υλοποίηση της μεθόδου GEO γίνεται στην πλατφόρμα FEniCS, η οποία είναι μια συλλογή ελεύθερου λογισμικού για την αυτοματοποιημένη, αποδοτική υλοποίηση διαφορικών εξισώσεων. Η υλοποίηση πραγματοποιείται με τρόπο τέτοιο που ο χρήστης να μπορεί εύκολα να ορίζει τις ιδιότητες των πεδίων του προβλήματος (όπως, γεωμετρία του προβλήματος, τελεστής μερικής διαφορικής εξίσωσης, οριακές συνθήκες, συνθήκες διεπαφών). Στο παρασκήνιο γίνεται η δημιουργία ή / και η βελτιστοποίηση των πλεγμάτων (τριγωνικά στοιχεία) για κάθε πεδίο του προβλήματος, η επίλυση των τοπικών προβλημάτων μερικών διαφορικών εξισώσεων και η εμφάνιση των υπολογισμένων τιμών για το πρόβλημα συνολικά και για τις διεπαφές. Η κυριότερη πρόκληση της υλοποίησης είναι η πρόσβαση στις τιμές των λύσεων στις διεπαφές και η εισαγωγή των χαλαρωμένων τιμών πίσω στα υποπροβλήματα για να αποτελέσουν τις οριακές συνθήκες των υποπροβλημάτων. Πραγματοποιούνται πειράματα σε προβλήματα-μοντέλα, ώστε να ελεγχθεί η εφαρμοσιμότητα και η σύγκλιση της μεθόδου. Στη συνέχεια, υλοποιείται μια παράλληλη έκδοση της μεθόδου και πραγματοποιούνται συγκριτικά πειράματα απόδοσης. Η παράλληλη υλοποίησή μας αποδεικνύεται πολύ χρήσιμη, ειδικά για μεγάλα προβλήματα. / Several approaches have been suggested in the literature for the solution of multiphysics / multidomain problems. Interface Relaxation (IR) methods are an interesting approach for the solution of these problems. Assuming initial guesses on the interfaces of the original problem, IR methods iteratively solve the subproblems and relax for new values on the interfaces until convergence is succeeded. Their main advantages are that their rates of convergence only depend on the parameters of the problem itself, the parameters related to its decomposition into subproblems and the parameters related to the operator imposed on the interfaces. In the present master thesis a new implementation of an IR method named GEO is presented. GEO is based on a simple geometric correction mechanism and acts iteratively so as to relax the values of the solution on the interfaces. In particular, it adds to the old interface values a geometrically weighted combination of the normal boundary derivatives of the adjacent subdomains. In this thesis GEO is implemented in FEniCS. The FEniCS project is a collection of free software for automated, efficient solution of differential equations. In order to evaluate the GEO implementation, it is applied on two different PDE problems with the same differential equation and boundary conditions and different domains. FEniCS methods are used to specify the problem's subdomains properties (i.e. geometry, PDE operator and boundary/interface conditions). They are also used to generate and/or refine meshes (triangular elements) for each subdomain, solve the local PDE problems and show the computed results in the global domain and on the interfaces. Getting values of the solutions on the interface (boundaries of the subproblems) and passing the new relaxed values back to the subproblems as updated values for the boundary conditions is the main challenge of the IR methodology implementation. Experiments are performed for 2-dimensional elliptic partial differential model problems with partitions in multiple subdomains and the results are examined in terms of the method's applicability and convergence. A parallel implementation of the GEO method using FEniCS is also presented, as well as its performance comparison to the serial implementation. The parallel implementation proves to be really useful, especially for large problems.