Η ηλεκτροεγκεφαλογραφία (EEG) και η μαγνητοεγκεφαλογραφία (MEG) είναι ιδιαίτερα χρήσιμες μέθοδοι εγκεφαλικών απεικονίσεων διότι έχουν πολύ καλή χρονική ανάλυση, της τάξεως του sec. Επειδή οι εγκεφαλικές διεργασίες εκτελούνται με μικρότερους ρυθμούς, οι εγκεφαλικές απεικονίσεις μέσω της EEG και MEG δίνουν τη δυνατότητα να παρακολουθούμε το λειτουργικό εγκέφαλο. Στο ευθύ πρόβλημα EEG γνωρίζουμε τη νευρωνική διέγερση που αναπτύσσεται εσωτερικά του εγκεφάλου και υπολογίζουμε το παραγόμενο ηλεκτρικό δυναμικό σε κάθε σημείο στο εξωτερικό ή στο εσωτερικό του εγκεφάλου. Αντίστοιχα στο ευθύ πρόβλημα MEG υπολογίζουμε το μαγνητικό δυναμικό εξωτερικά του εγκεφάλου. Στο αντίστροφο πρόβλημα της EEG γνωρίζουμε το ηλεκτρικό δυναμικό, από μετρήσεις, εξωτερικά του εγκεφάλου, στην επιφάνεια του κρανίου, και ζητάμε να προσδιορίσουμε τη νευρωνική διέγερση που το προκάλεσε. Για τη MEG γνωρίζουμε το μαγνητικό δυναμικό που καταγράφεται εξωτερικά του κρανίου και ζητάμε τη νευρωνική διέγερση που αναπτύχθηκε εσωτερικά του εγκεφάλου. Στην παρούσα εργασία επιλύονται οκτώ προβλήματα. Το ευθύ πρόβλημα και το αντίστροφο πρόβλημα της EEG και της MEG σε δύο περιπτώσεις. Στην πρώτη περίπτωση η νευρωνική διέγερση εντοπίζεται σε ένα μικρό ευθύγραμμο τμήμα, δηλαδή το ρεύμα που δημιουργήθηκε εσωτερικά του εγκεφάλου λόγω μιας εγκεφαλικής διεργασίας αναπαρίσταται με δίπολα που κατανέμονται κατά μήκος ενός μικρού ευθύγραμμου τμήματος. Επιλύουμε αρχικά το ευθύ πρόβλημα της EEG και στη συνέχεια το αντίστροφο πρόβλημα. Καταλήγουμε σε ένα μη γραμμικό σύστημα που, στη γενική περίπτωση, επιλύεται αριθμητικά για να υπολογίσει τη θέση και τη ροπή της διπολικής πηγής, τον προσανατολισμό και το μήκος του ευθύγραμμου τμήματος. Αναλυτικά επιλύουμε δύο ειδικές περιπτώσεις και βρίσκουμε μοναδική λύση. Στην πρώτη ειδική περίπτωση το ευθύγραμμο τμήμα είναι παράλληλο στον άξονα, ενώ στη δεύτερη στο άξονα και το κέντρο του, και στις δύο περιπτώσεις, είναι πάνω στο άξονα . Το ευθύ και το αντίστροφο πρόβλημα το επιλύουμε και για την περίπτωση της MEG δίνοντας αναλυτικά τη λύση στην περίπτωση που το ευθύγραμμο τμήμα είναι παράλληλο στο άξονα και το κέντρο του είναι επάνω στο άξονα. Στη δεύτερη περίπτωση μελετάμε ακριβώς τα ίδια προβλήματα όταν η νευρωνική διέγερση εντοπίζεται σε ένα μικρό κυκλικό δίσκο που το επίπεδό του είναι κάθετο στο διάνυσμα θέσης του κέντρου του δίσκου. Επιλύουμε το πρόβλημα σε συγκριμένη θέση του δίσκου για να απλοποιηθούν οι υπολογισμοί. Εφαρμόζουμε κατάλληλες στροφές Euler ώστε το επίπεδό του να βρεθεί σε θέση παράλληλη στο επίπεδο και το κέντρο του δίσκου να βρίσκεται πάνω στον άξονα. Στη συνέχεια εφαρμόζουμε στροφές Euler και το επαναφέρουμε στην αρχική θέση. Αφού επιλύσουμε τα ευθύ προβλήματα της EEG και MEG ξεχωριστά, προσδιορίζουμε, επιλύοντας το αντίστροφο, τη θέση του κυκλικού δίσκου. Καταλήγουμε σε ένα σύστημα μη γραμμικό που απαιτεί αριθμητική επίλυση τόσο για την EEG όσο και για τη MEG Αναλυτικά, επιλύουμε το αντίστροφο για την EEG σε μια ειδική περίπτωση κατά την οποία ο δίσκος είναι παράλληλος στο επίπεδο και το κέντρο του βρίσκεται πάνω στον άξονα και βρίσκουμε τη μοναδική λύση η οποία προσδιορίζει τη θέση του δίσκου, την ακτίνα του και τη ροπή της διπολικής πηγής. Με την επίλυση των προβλημάτων αυτών επιβεβαιώνουμε και τα αποτελέσματα της εργασίας των Albanese και Monk. Συγκεκριμένα, έδειξαν ότι δε μπορεί να προσδιοριστεί ο φορέας του ρεύματος που εντοπίζεται σε χώρο τριών διαστάσεων. Στην παρούσα διατριβή καθορίζουμε την έκτασης της νευρωνικής διέγερσης όταν η διάσταση του φορέα της είναι μικρότερη του τρία. / Electroencephalography (EEG) and Magnetoencephalography (MEG) are the two brain imaging modalities which have the necessary temporal resolution, sec for the study of the functional brain. Albanese and Monk have demonstrated that it is impossible to identify the extent of a localized three-dimensional current distribution lying inside a three-dimensional conductive medium. The purpose of the present work is to show that, as already predicted by Albanese and Monk, this result is not true if the current distribution is restricted on a one or two- dimensional set. The calculation of the values of the electric potential on the surface of the head defines the forward problem of EEG, while the calculation of the magnetic flux density a few centimeters outside the head defines the forward problem of MEG. The inverse EEG problem seeks to identify the neuronal current within the brain from the knowledge of the electric potential on the surface of the head. The corresponding inverse MEG problem seeks this neuronal current from the knowledge of the magnetic flux outside the head. In the present dissertation we study eight particular problems. They concern the forward and the inverse problem of EEG and MEG in two special geometric cases. In one geometrical case the neuronal current is supported on a small line segment and the neuronal current is represented by a dipole distribution along this line segment. First we solve the forward EEG problem and then we solve the inverse problem of identifying the location, the orientation, the size and the average dipolar moment over the line segment. We arrive at a nonlinear algebraic system which we solve analytically in two special cases. Next we solve the corresponding forward and the inverse MEG problems for the same structure. A second case concerns the relative EEG and MEG problems when the current is supported on a small disc normal to a radius of the conducting sphere. As before, we solve the EEG and MEG problems separately and then we solve the inverse EEG and MEG problems which determine the position, the orientation and the size of the disk supporting the primary neuronal current.

Η παρούσα διδακτορική διατριβή πραγματεύεται τη λειτουργική συμπεριφορά μιας σύγχρονης μηχανής με έκτυπους πόλους κατά τη διάρκεια δύο ειδών σφαλμάτων, τα οποία παρουσιάζονται παρακάτω, για τις δύο συνθήκες σύνδεσης μίας σύγχρονης γεννήτριας με το δίκτυο. Ειδικότερα, μελετήθηκε η περίπτωση βραχυκυκλώματος στην τροφοδοσία του τυλίγματος διέγερσης της σύγχρονης μηχανής όταν είναι συνδεδεμένη σε ισχυρό δίκτυο και είτε ο αριθμός στροφών αυτής διατηρείται απόλυτα σταθερός, είτε ένας ελεγκτής τη συγκρατεί στο σύγχρονο αριθμό στροφών. Ακόμη διερευνήθηκε η περίπτωση εσωτερικού σφάλματος στο τύλιγμα του στάτη για τις δύο προαναφερθείσες περιπτώσεις. Υπολογίστηκαν και μελετήθηκαν η ηλεκτρομαγνητική ροπή και η μαγνητική επαγωγή κατά τη διάρκεια κάθε φαινομένου, καθώς επίσης τα ρεύματα σε δρομέα και στάτη συμπεριλαμβανομένων των ρευμάτων στον κλωβό απόσβεσης και του ρεύματος βραχυκύκλωσης σε τμήματα του τυλίγματος του στάτη. Αρχικά, περιγράφεται αναλυτικά ο τρόπος με τον οποίο μοντελοποιήθηκε η σύγχρονη γεννήτρια με έκτυπους πόλους, ενώ παράλληλα αναλύεται η μέθοδος με την οποία μοντελοποιήθηκαν τα σφάλματα σε στάτη και δρομέα και ο τρόπος με τον οποίο προσομοιώθηκε στο πρόγραμμα πεπερασμένων στοιχείων ο παραλληλισμός της μηχανής στο δίκτυο παροχής ηλεκτρικής ενέργειας. Στη συνέχεια, γίνεται αναφορά στον τρόπο με τον οποίο ορίζονται οι περιοχές του μοντέλου, οι εξισώσεις που επιλύει το πρόγραμμα πεπερασμένων στοιχείων για την εξαγωγή των αποτελεσμάτων, ο ορισμός των οριακών συνθηκών και τέλος περιγράφεται η εφαρμογή της μεθόδου πεπερασμένου στοιχείων που εφαρμόσθηκε στο συγκεκριμένο μοντέλο. Έπειτα προσομοιώθηκε η περίπτωση βραχυκυκλώματος στην τροφοδοσία του τυλίγματος διέγερσης, καθώς ο στάτης της σύγχρονης μηχανής είναι συνδεδεμένος σε ισχυρό δίκτυο με το δρομέα να στρέφεται με σταθερό αριθμό στροφών. Κατά τη διάρκεια του σφάλματος πραγματοποιείται λεπτομερής καταγραφή της μαγνητικής επαγωγής και της ηλεκτρομαγνητικής ροπής, καθώς και όλων των ρευμάτων στο τύλιγμα στάτη και δρομέα. Προκύπτουν χρήσιμα συμπεράσματα για τη συμπεριφορά της μηχανής σε όλη τη διάρκεια του σφάλματος και καταγράφονται τα ηλεκτρομαγνητικά μεγέθη και γίνεται αξιολόγηση της συμπεριφοράς της κατά τη διάρκεια αυτού του μεταβατικού φαινομένου. Ακόμη, εξετάζεται η ίδια περίπτωση σφάλματος στο τύλιγμα διέγερσης, αλλά ο αριθμός στροφών του δρομέα διατηρείται σταθερός μέσω ενός ελεγκτή στροφών. Παρατηρείται ότι η συμπεριφορά της μηχανής αλλά και όλα τα ηλεκτρομαγνητικά μεγέθη αυτής είναι τελείως διαφορετικά συγκρινόμενα με την προηγούμενη περίπτωση. Στη συνέχεια αυτής της διδακτορικής διατριβής εξετάζεται η συμπεριφορά της σύγχρονης μηχανής στην περίπτωση εσωτερικού βραχυκυκλώματος στο τύλιγμα του στάτη, καθώς είναι συνδεδεμένη σε ισχυρό δίκτυο με το δρομέα να στρέφεται με σταθερό αριθμό στροφών. Αναλυτικότερα, μελετώνται τα ρεύματα στο τύλιγμα του δρομέα καθώς και στο τύλιγμα του στάτη για βραχυκύκλωμα μεταξύ σπειρών που ανήκουν σε ίδια ή διαφορετική φάση. Υπολογίζεται το ρεύμα βραχυκύκλωσης και μελετάται ο τρόπος που επηρεάζει τα φασικά ρεύματα του στάτη τόσο ποσοτικά όσο και ποιοτικά. Επίσης, αναλύονται τα ρεύματα στον κλωβό απόσβεσης και μελετάται η συμπεριφορά τους κατά τη διάρκεια του βραχυκυκλώματος. Η εισαγωγή ελέγχου στροφών μεταβάλει τη συμπεριφορά της σύγχρονης γεννήτριας και υπολογίζονται αναλυτικά τα ηλεκτρομαγνητικά μεγέθη της σύγχρονης μηχανής προκύπτοντας συμπεράσματα για τον τρόπο με τον οποίο το βραχυκυκλωμένο τύλιγμα επηρεάζει τα μεγέθη αυτά, ενώ παράλληλα καθορίζεται ο ρόλος που έχουν οι φάσεις που συμμετέχουν στο βραχυκύκλωμα. Τέλος, γίνεται μια σύντομη σύγκριση για το πόσο επηρεάζει ο αριθμός σπειρών τη συμπεριφορά της σύγχρονης μηχανής στην περίπτωση σφάλματος στο στάτη κατά την περίπτωση κατά την οποία η μηχανή είναι συνδεδεμένη σε ισχυρό δίκτυο με σταθερό τον αριθμό στροφών του δρομέα. Πιο συγκεκριμένα, αναλύονται τα ρεύματα σε στάτη και δρομέα καθώς και η ηλεκτρομαγνητική ροπή, για τις περιπτώσεις που οι βραχυκυκλωμένες σπείρες ανήκουν στην ίδια και διαφορετική φάση αλλά με διαφορετικό αριθμό βραχυκυκλωμένων σπειρών για κάθε περίπτωση. Προκύπτει το συμπέρασμα ότι καθοριστικό ρόλο στη διαμόρφωση των ηλεκτρομαγνητικών μεγεθών κατά τη διάρκεια του σφάλματος έχει ο αριθμός των βραχυκυκλωμένων σπειρών, ανεξάρτητα από τον αριθμό των φάσεων που συμμετέχουν στο βραχυκύκλωμα. / This thesis deals with the functional behavior of a salient pole synchronous generator during two kinds of short-circuits. In particular, we studied the case of short circuit in the supply of the excitation winding of the synchronous machine when it is connected to an infinite bus and either the rotor speed was absolutely constant, or a simple PI- Controller maintained the synchronous speed equal to the synchronous. Additionally, the case of an internal fault in the stator winding for the two cases mentioned above was investigated. The electromagnetic torque and the magnetic flux density in each case were calculated and studied, as well as the stator and rotor currents, including the damper cage, and the short-circuit current in the faulty loop of the stator winding. Firstly, it is described in detail the way in which the salient pole synchronous generator was modeled and it is analyzed the method by which the faults are modeled in both stator and rotor and the way these faults were simulated, in the finite element program. Additionally, it is presented the way in how the areas of the model are defined, the equations that were solved through finite element software, in order to extract the results, the definition of the boundary conditions and finally it is described the finite element method, which was applied to this specific model. The case of a short circuit in the supply of the field winding while the stator of the synchronous machine is connected to the grid and the rotor speed is held constant and equal to the synchronous one, is examined. During this fault the magnetic flux, the electromagnetic torque and all the stator and rotor currents are measured in detail. Useful conclusions about the behavior of the machine throughout this kind of short-circuit were derived, all the electromagnetic magnitudes were recorded and an assessment of the generator behavior during this transient phenomenon is made. Similarly, the same type of fault is analyzed, but the speed of the rotor is maintained constant through a speed controller. It is observed that the behavior of the machine and all the electromagnetic magnitudes are quite different compared to the previous case. In this dissertation is examined the behavior of the hydrogenerator in the case of an inter-turn short circuit in the stator winding, while it is connected to the grid with a constant rotor speed. Specifically, it is examined the currents in the rotor and the stator winding for a short circuit between turns that belong to the same or to different phases. The short circuit current is calculated and it is presented the way that it affects quantitative and qualitative the stator phase currents. It is also analyzed the damper currents and it is studied their behavior during the short circuit. The speed controller alters the behavior of the synchronous generator and all the electromagnetic magnitudes of this machine are analytically calculated, resulting significant conclusions on how the faulty loop affects these quantities, while it is set out the role of the participating phases in this short-circuit. Finally, a brief comparison of the way that the number of the shorted turns affects the behavior of the simulated machine in the case of an inter-turn stator fault, while it is connected to the grid with a fixed number of the rotor revolutions. Specifically, it is analyzed the stator and rotor currents and the electromagnetic torque, for the cases that the short-circuited turns belong either to the same or to different phases, but with different number of shorted turns. It is concluded that a key role in determining the electromagnetic magnitudes during this fault has the number of the short-circuited turns and not the number of the phases that are involved in the short circuit.

Στην παρούσα διπλωματική εργασία εξετάζονται τα χαρακτηριστικά των υποθαλάσσιων αποθέσεων της ερυθράς ιλύος που καλύπτουν τεράστια έκταση του πυθμένα του κεντρικού Κορινθιακού κόλπου και προκύπτουν ως μεταλλευτικό κατάλοιπο της επεξεργασίας βωξιτών για την παραγωγή αλουμινίου. / In the present diplomatic work are examined the characteristics of submarine red mud tailings that cover enormous extent of seabed of central corinthian gulf and result as mining residue of treatment of bauxites for the production of aluminium.

Στην παρούσα διατριβή διερευνάται η επίδραση παράκτιων μη-θραυόμενων κυματισμών στη μορφολογική ισορροπία αμμώδους πυθμένα με πτυχώσεις, η θραύση εκχείλισης κυμάτων πάνω από πυθμένα σταθερής κλίσης, καθώς και τα συνεπαγόμενα κυματογενή ρεύματα στη ζώνη απόσβεσης. Για το σκοπό αυτό αναπτύσσονται μοντέλα αριθμητικής προσομοίωσης, τα οποία πραγματοποιούν επίλυση των δισδιάστατων και τρισδιάστατων εξισώσεων ασυμπίεστης, συνεκτικής ροής με ελεύθερη επιφάνεια. Η αριθμητική επίλυση των εξισώσεων ροής, Navier-Stokes, επιτυγχάνεται με τη χρήση κλασματικής μεθόδου για τη χρονική ολοκλήρωση, ενώ η χωρική διακριτοποίηση πραγματοποιείται μέσω ενός υβριδικού σχήματος πεπερασμένων διαφορών και ψευδο-φασματικών μεθόδων προσέγγισης. Στις προσομοιώσεις της θραύσης εκχείλισης κύματος γίνεται χρήση της μεθόδου προσομοίωσης μεγάλων κυμάτων LWS, σύμφωνα με την οποία επιλύονται μόνο οι μεγάλες χωρικές διακυμάνσεις της ταχύτητας και της ελεύθερης επιφάνειας, ενώ η επίδραση των μικρότερων διακυμάνσεων περιγράφεται μέσω ενός μοντέλου διατμητικών τάσεων υποκλίμακας (SGS), ανάλογα με ότι ισχύει στη μέθοδο προσομοίωσης μεγάλων δινών, LES. Ένα ανεξάρτητο μοντέλο για την προσομοίωση της μεταβολής μορφολογίας πυθμένα, μέσω μεταφοράς φορτίου πυθμένα, αναπτύσσεται και χρησιμοποιείται σε σύζευξη με τα μοντέλα προσομοίωσης δισδιάστατης ροής. H παροχή του φορτίου πυθμένα υπολογίζεται μέσω τροποποίησης γνωστών εμπειρικών σχέσεων, σε συνδυασμό με τη στιγμιαία διατμητική τάση πυθμένα από τη μονάδα προσομοίωσης της ροής. Από τις προσομοιώσεις ροής πάνω από πυθμένα με πτυχώσεις, προκύπτει ότι η παρουσία των πτυχώσεων επηρεάζει σημαντικά το κυματογενές οριακό στρώμα, ενώ οι μορφολογικές προσομοιώσεις οδηγούν στο συμπέρασμα ότι, η μακροπρόθεσμη ισορροπία των πτυχώσεων επέρχεται για συγκεκριμένη τιμή της γωνίας/συντελεστή δυναμικής τριβής, η οποία συσχετίζεται με τις διαστάσεις των πτυχώσεων και τα χαρακτηριστικά του κύματος. Για τη θραύση εκχείλισης εξετάζονται οι περιπτώσεις κάθετης αλλά και υπό γωνία, ως προς την ακτογραμμή, διάδοσης κυμάτων πάνω από πυθμένα σταθερής κλίσης 1/35. Τα αποτελέσματα για τα χαρακτηριστικά της κάθετης θραύσης (ύψος και βάθος θραύσης, Ηb και db, αντίστοιχα) και του συνεπαγόμενου υποβρύχιου ρεύματος, συγκρίνονται με δημοσιευμένες πειραματικές μετρήσεις και η συμφωνία είναι ικανοποιητική. Το μοντέλο είναι σε θέση να προσομοιώσει την ανάπτυξη του επιφανειακού στροβίλου στο μέτωπο του θραυόμενου κύματος, η οποία συνοδεύεται από αύξηση της ισχύος των SGS τάσεων (μέχρι βάθους d/db ≈ 0.75) και διαδοχική μείωσή τους, μέχρι μηδενισμού, στα ρηχά της ζώνης απόσβεσης. Από τα αποτελέσματα για το πεδίο στροβιλότητας και τις SGS τάσεις, κατά την προσομοίωση της υπό γωνία θραύσης, παρατηρείται η σταδιακή θραύση του κύματος κατά μήκος της κορυφογραμμής, ενώ προκύπτει ότι οι τελευταίες παραμένουν ενεργές για περίπου δύο μήκη κύματος. Επίσης, η μέση ταχύτητα του παράλληλου ρεύματος προκύπτει πιο ενισχυμένη σε ρηχά βάθη στη ζώνη απόσβεσης (d/db < 0.5), ενώ η κατακόρυφη κατανομή του παρουσιάζεται σαφώς επηρεασμένη από την παρουσία του υποβρύχιου ρεύματος κοντά στον πυθμένα. / In the present thesis, the impact of nearshore, non-breaking waves on the morphological equilibrium of small scale patterns that appear in sandy beds, well-known as ripples, the spilling wave breaking over a bed of constant slope and the wave-induced currents developing in the surf zone, are investigated. Numerical models are developed for the simulation of the aforementioned phenomena, based on the numerical solution of the two/three-dimensional, incompressible, viscous, free-surface flow. The numerical solution of the flow equations, i.e. the Navier-Stokes equations, is accomplished by means of a time-splitting scheme of three stages for the temporal discretization and a hybrid scheme for the spatial discretization, with central finite differences along the streamwise direction and pseudo-spectral approximations, with Fourier modes and Chebyshev polynomials along the spanwise and vertical directions, respectively. The simulations of spilling wave breaking are performed employing the so-called large-wave simulation (LWS) method, according to which, large velocity and free-surface scales are fully resolved, while the effect of subgrid scales is modeled by eddy-viscosity stresses, similar to large-eddy simulation (LES) methodology. The bed morphology evolution, driven by the bed load sediment transport, is simulated by a morphology model, which performs the numerical solution of the sediment mass conservation equation, utilized coupled with the two-dimensional flow model. The bed load transport rate, is computed inserting bed shear stress timeseries (by the flow model) into published empirical formulas that have been modified to fit the characteristics of the investigated cases. For the case of rippled bed, it was found that the structure of the wave boundary layer is substantially influenced by the presence of the ripples, and that for a certain value of the dynamic friction angle/coefficient, which is correlated to the ripple dimensions and the wave characteristics, the ripples remain in quasi-steady equilibrium after each wave period. Wave breaking is investigated by the simulation of normal and oblique to the shoreline propagation, transformation and spilling breaking of incoming Stokes waves, over a bed of constant slope, tanβ = 1/35. Our numerical results are compared to published experimental measurements, and it is found that the LWS model predicts adequately the wave breaking parameters - breaking height, Ηb, and depth, db- and the distribution of the undertow current in the surf zone. The development of the surface roller in the breaking wavefront is also captured, while is connected to the increase of the strength of the sub-grid (SGS) stresses in the outer surf zone (up to d/db ≈ 0.75) and their successive decrease at shallower depths close to the shoreline. For the case of oblique wave breaking, the vorticity and the SGS stresses distribution in the surf zone clearly indicate the gradual breaking along the wave crestline, while the latter (SGS stresses) remain active for about two wavelenghts. Finally, the magnitude of the longshore current is found to be enhanced at shallower depths in the surf zone (d/db < 0.5), while its vertical distribution is affected by the interaction with the undertow current.

Αντικείμενο της διατριβής είναι η προσομοίωση της ροής και της μεταφοράς θερμότητας σε ηλιακά και ενεργειακά συστήματα. Η έμφαση δόθηκε στις δεξαμενές αποθήκευσης της θερμότητας που παράγεται στα συγκεκριμένα συστήματα, με στόχο τον χαρακτηρισμό των ενεργειακών απωλειών και την βελτιστοποίηση του σχεδιασμού τους. Κύριες δραστηριότητες της διατριβής θα είναι η περαιτέρω ανάπτυξη διαθέσιμων εργαλείων προσομοίωσης ροών φυσικής και μικτής συναγωγής, με διερεύνηση των νεώτερων εξελίξεων στην μοντελοποίηση με τη μέθοδο Προσομοίωσης Μεγάλων Δινών (LES). Αρχικά γίνεται εκτεταμένη επικύρωση με πειραματικά αποτελέσματα σε απλές γεωμετρικές διατάξεις (π.χ. ορθογωνικά κανάλια ή κοιλώματα με βαθμίδα θερμοκρασίας) από την βιβλιογραφία. Στη συνέχεια η μεθοδολογία εφαρμόζεται στον υπολογισμό ροών σε πιο ρεαλιστικές γεωμετρίες, επιλεγμένες από πρακτικές εφαρμογές, όπως οι δεξαμενές αποθήκευσης νερού. Αναλύονται σε βάθος οι δυναμικές διεργασίες και τα ροϊκά φαινόμενα τόσο κατά την προσαγωγή της θερμότητας στη δεξαμενή (φόρτιση) όσο και κατά την απαγωγή της (εκφόρτιση) και η επίδραση που έχουν αυτά στην αποδοτικότητα της αποθήκευσης με βάση κατάλληλους ποσοτικούς δείκτες. Από τα αποτελέσματα αναδεικνύεται η σημασία της μοντελοποίησης σε τέτοιου είδους συστήματα ως ένα σημαντικό εργαλείο στη διερεύνηση της απόδοσης τους, του ενεργειακού χαρακτηρισμού τους και ακολούθως στην προσπάθεια επίτευξης του βέλτιστου σχεδιασμού τους. / The subject of the thesis is the Simulation of Turbulent Flow and Heat Transfer in Solar and Energy Systems. Emphasis is given in the thermal storage component of these systems, with the aim of characterizing their energy losses and improve their design. Main activities of the thesis will be the further development of available computational tools for the simulation of flows in natural and mixed convection, incorporating some of the most recent developments in modeling, particularly in the Large Eddy Simulation (LES) method. Initially, an extensive validation with experimental results in simple geometric configurations is carried out (e.g. channels or differentially heated cavities). Subsequently, the methodology is applied in the calculation of flows for more realistic geometries selected from practical applications, such as various hot water storage tanks. Analysis is conducted of the dynamic processes and relevant physical phenomena during the heat supply (charging) to and removal (discharging) from the tank and their influence on the storage effectiveness using appropriate thermodynamic indices. From the simulation results, the significance of the flow and heat transfer modeling in these systems as a practical tool for studying their performance is demonstrated, by characterizing their energy content and significantly contributing to the process of optimizing their design.