We study the flow of a dilute polymer solution in a wavy channel under steady-state flow conditions by employing the non-equilibrium thermodynamics two-fluid model (Mavrantzas-Beris, 1992), allowing for the coupling between polymer concentration and polymer stresses. The resulting highly complex system of partial differential equations describing inhomogeneous transport phenomena in the fluid are solved with an efficient implementation of the mixed finite-element method. We present numerical results for polymer concentration, stress, velocity and fluxes of polymer as a function of the non-dimensional parameters of the problem (the Deborah number , the Peclet number , the Reynolds number , the ratio of the solvent viscosity to the total fluid viscosity , and the constriction ratio of the channel width ). We find that the constricted part of the wall is depleted of polymer, when the polymer diffusion length-scale, expressed by the ratio of / , increases. The migration is more pronounced for macromolecules characterized by longer relaxation times, and takes place towards the expanding part of the channel or towards the centerplane. Migration is also enhanced by the width variability of the channel: the more corrugated the channel, the stronger the transfer of polymer to the centerplane. This increases the spatial extent of polymer depletion near the wall or induces a zone of sharp variation in polymer stress and concentration, which moves away from the channel wall, especially in lower polymer concentration. The development of a polymer-depleted layer smooths out the boundary layer which is known to arise with Boger fluids at the walls of such corrugated channels or tubes and gives rise to an “apparent” slip in the constricted section of the wall and to a very low value of the drag force on the wall. When and where boundary layers arise, they scale as (1/De) for the stresses and as (De⁄Pe)^(1⁄3) for the concentration. / --

- / -

- / -

- / -

- / -

- / -

O σκοπός της παρούσας διδακτορικής διατριβής είναι διττός και, εμπεριέχει δραστηριότητα τόσο στο κομμάτι της εφαρμοσμένης όσο και της βασικής έρευνας. Πιο συγκεκριμένα, περιλαμβάνει την εφαρμογή σύγχρονων υπολογιστικών μεθόδων (Μέθοδος Πεπερασμένων Στοιχείων, Μέθοδος Πεπερασμένων Όγκων) στη μελέτη της ροής του αίματος, καθώς και την ανάπτυξη σύγχρονων υπολογιστικών μεθοδολογιών που δε στηρίζονται στη χρήση πλέγματος. Ταυτόχρονα μελετάται η αποτελεσματική εφαρμογή των μεθόδων της Υπολογιστικής Ρευστομηχανικής στην ιατρική πρακτική (Simulation Treatment Planning), για την αποτελεσματική πρόληψη, διάγνωση και θεραπευτική αντιμετώπιση των νόσων του καρδιακού και περιφερικού αγγειακού συστήματος. Η διαδικασία υλοποίησης των προσομοιώσεων έχουν σκοπό την υποβοήθηση του θεράποντα ιατρού στη λήψη ιατρικής απόφασης σχετικά με τη θεραπευτική αγωγή. Παράλληλα, στο κομμάτι της βασικής έρευνας αναπτύσσονται σύγχρονες υπολογιστικές μέθοδοι οι οποίες πρόκειται να καλύψουν τις αδυναμίες που παρουσιάζουν οι διαδεδομένες υπολογιστικές μέθοδοι. Η υλοποίησή τους βρίσκει ανταπόκριση και πεδίο εφαρμογής σε διάφορους τομείς της Επιστήμης και της Μηχανικής. Έτσι, θα παρουσιαστούν αποτελέσματα μόνο στο τομέα της Υπολογιστικής Ρευστομηχανικής. Τα κριτήρια αυτά θα στηριχθούν σε ποσοτικές συνδυαστικές αναλύσεις, οι οποίες ενσωματώνουν το state-of-the-art της μορφολογικής απεικόνισης των αγγείων στο state-of-the-art των μεθοδολογιών της υπολογιστικής ρευστομηχανικής, οι οποίες είναι άμεσα σχετιζόμενες με την εκτίμηση των αιμοδυναμικών παραγόντων και τάσεων που εφαρμόζονται σε πάσχουσες περιοχές του αγγειακού συστήματος, όπως στενώσεις, θρόμβοι και ανευρύσματα. / The aim of present doctoral thesis is double fold and, includes research activity both in the fields of applied and basic research. More precisely, it includes the application of modern numerical methods (Finite Element Method, Finite Volume Method) for the study of blood flow, as well as the development of modern numerical methodologies, which do not rely on the use of a computational mesh, that is the so-called meshless or Meshfree methods. Furthermore, the effective application of sophisticated numerical methods in the medical practice (Simulation Treatment Planning) has been studied, since there is a great necessity for effective prevention, diagnosis and therapeutic confrontation of illnesses the cardiac and vascular system. The simulation conducted they aim to assist the doctor in the decision-making. At the same time, regarding of area of the basic research, sophisticated numerical methods were developed and applied to various applications of science and engineering. More precisely, results will be presented for Computational Fluid Dynamics problems.

Η ανάπτυξη που σημειώθηκε τις τελευταίες δεκαετίες στον τομέα της εξόρυξης υδρογονανθράκων οδήγησε στην υιοθέτηση όλο και πιο σύγχρονων τεχνικών διαχωρισμού των προϊόντων στην πηγή. Οι φυγοκεντρικοί διαχωριστές είναι η πιο επιτυχημένη και διαδεδομένη εφαρμογή. Σε αυτή την γκάμα, έρχεται να προστεθεί και ο πρότυπος φυγοκεντρικός διαχωριστής που παρουσιάζεται στην παρούσα εργασία. Η ανάλυση προσπαθεί να ισορροπήσει μεταξύ μίας «αυστηρά επιστημονικής» ανάλυσης του στροβιλώδους πεδίου διφασικής ροής μέσα στο διαχωριστή και μιας πιο «πρακτικής» ανάλυσης της απόδοσης της συγκεκριμένης συσκευής. Η πρώτη είναι απαραίτητη λόγω της πρωτοτυπίας που παρουσιάζει το συγκεκριμένο επιστημονικό πεδίο και της έλλειψης προγενέστερων μελετών που να καλύπτουν το θέμα, ενώ η δεύτερη είναι επιβεβλημένη λόγω της «πρακτικότητας» που προσδοκούμε να έχει το τελικό αποτέλεσμα. Ο προσδιορισμός του προβλήματος που θέτει προς επίλυση η παρούσα εργασία και η ανασκόπηση βιβλιογραφίας που καθορίζει σε μεγάλο βαθμό την κατεύθυνση που θα ακολουθήσει η έρευνα, γίνεται στο 1ο κεφάλαιο. Για την αντιμετώπιση του προβλήματος ακολουθείται η τυπική διαδικασία εξέλιξης ενός βιομηχανικού προϊόντος σε τέσσερα στάδια. Στο 1ο στάδιο αρχικά αναλύεται η ιδέα –concept- λειτουργίας και η θεωρία πάνω στην οποία βασίζεται. Στο 2ο στάδιο γίνεται αρχική διερεύνηση λειτουργίας του διαχωριστή, με προσομοίωση της ροής με Υπολογιστή Ρευστομηχανική. Στο 3ο στάδιο γίνεται πειραματική – παραμετρική ανάλυση της λειτουργίας του διαχωριστή. Τέλος, στο 4ο στάδιο γίνεται επιβεβαίωση – πιστοποίηση της τεχνικής προσομοίωσης. Στο 2ο κεφάλαιο γίνεται μία θεωρητική ανάλυση των εξισώσεων που διέπουν τη ροή των ρευστών. Οι εξισώσεις κίνησης αναπτύσσονται και εξηγείται πως διαμορφώνονται μετά την ανάπτυξη της τύρβης. Στο σημείο αυτό γίνεται αντιληπτό ότι δεν είναι δυνατή η επίλυσή τους με αναλυτικές μεθόδους αλλά με αριθμητικές. Αυτή είναι και η απαρχή της εισόδου των αριθμητικών μεθόδων στη μηχανική ρευστών. Στο 3ο κεφάλαιο παρατίθεται το θεωρητικό υπόβαθρο πάνω στο οποίο βασίζεται η ανάλυση των διφασικών ροών αλλά και των στροβιλωδών ροών, έννοιες που συναντώνται στα επόμενα κεφάλαια. Η θεωρητική ανάλυση ολοκληρώνεται στο 4ο κεφάλαιο με την εμβάθυνση στην αριθμητική επίλυση των εξισώσεων του 2ου κεφαλαίου. Αρχίζουμε με μία σύντομη παρουσίαση των μεθόδων που αναπτύχθηκαν τις τελευταίες δεκαετίες, επικεντρώνοντας στη μέθοδο που χρησιμοποιούμε στην παρούσα εργασία. Στη συνέχεια παρατίθενται διάφορα μοντέλα προσομοίωσης ροών επικεντρώνοντας στις διφασικές ροές. Επειδή το αντικείμενο αυτό περιλαμβάνει εκτεταμένες μαθηματικές αναλύσεις, μεταφέρθηκαν αρκετά τμήματα στο παράρτημα 1, ώστε το κυρίως κείμενο της εργασίας να παρουσιάζει ομοιογένεια και να διατηρεί την ισορροπία μεταξύ της θεωρητικής ανάλυσης και της πρακτικής – ερευνητικής προσέγγισης που αναφέρθηκε προηγουμένως. Από το 5ο κεφάλαιο ξεκινά η παρουσίαση του καθαρά ερευνητικού έργου και ακολουθείται η λογική και χρονολογική σειρά που αναλύθηκε πριν, με τα τέσσερα στάδια εξέλιξης του «προϊόντος». Σε αυτό λοιπόν το κεφάλαιο αναλύεται η ιδέα –concept- λειτουργίας και η θεωρία πάνω στην οποία βασίζεται και παρουσιάζονται σχέδια από την Η ανάπτυξη που σημειώθηκε τις τελευταίες δεκαετίες στον τομέα της εξόρυξης υδρογονανθράκων οδήγησε στην υιοθέτηση όλο και πιο σύγχρονων τεχνικών διαχωρισμού των προϊόντων στην πηγή. Οι φυγοκεντρικοί διαχωριστές είναι η πιο επιτυχημένη και διαδεδομένη εφαρμογή. Σε αυτή την γκάμα, έρχεται να προστεθεί και ο πρότυπος φυγοκεντρικός διαχωριστής που παρουσιάζεται στην παρούσα εργασία. Η ανάλυση προσπαθεί να ισορροπήσει μεταξύ μίας «αυστηρά επιστημονικής» ανάλυσης του στροβιλώδους πεδίου διφασικής ροής μέσα στο διαχωριστή και μιας πιο «πρακτικής» ανάλυσης της απόδοσης της συγκεκριμένης συσκευής. Η πρώτη είναι απαραίτητη λόγω της πρωτοτυπίας που παρουσιάζει το συγκεκριμένο επιστημονικό πεδίο και της έλλειψης προγενέστερων μελετών που να καλύπτουν το θέμα, ενώ η δεύτερη είναι επιβεβλημένη λόγω της «πρακτικότητας» που προσδοκούμε να έχει το τελικό αποτέλεσμα. Ο προσδιορισμός του προβλήματος που θέτει προς επίλυση η παρούσα εργασία και η ανασκόπηση βιβλιογραφίας που καθορίζει σε μεγάλο βαθμό την κατεύθυνση που θα ακολουθήσει η έρευνα, γίνεται στο 1ο κεφάλαιο. Για την αντιμετώπιση του προβλήματος ακολουθείται η τυπική διαδικασία εξέλιξης ενός βιομηχανικού προϊόντος σε τέσσερα στάδια. Στο 1ο στάδιο αρχικά αναλύεται η ιδέα –concept- λειτουργίας και η θεωρία πάνω στην οποία βασίζεται. Στο 2ο στάδιο γίνεται αρχική διερεύνηση λειτουργίας του διαχωριστή, με προσομοίωση της ροής με Υπολογιστή Ρευστομηχανική. Στο 3ο στάδιο γίνεται πειραματική – παραμετρική ανάλυση της λειτουργίας του διαχωριστή. Τέλος, στο 4ο στάδιο γίνεται επιβεβαίωση – πιστοποίηση της τεχνικής προσομοίωσης. Στο 2ο κεφάλαιο γίνεται μία θεωρητική ανάλυση των εξισώσεων που διέπουν τη ροή των ρευστών. Οι εξισώσεις κίνησης αναπτύσσονται και εξηγείται πως διαμορφώνονται μετά την ανάπτυξη της τύρβης. Στο σημείο αυτό γίνεται αντιληπτό ότι δεν είναι δυνατή η επίλυσή τους με αναλυτικές μεθόδους αλλά με αριθμητικές. Αυτή είναι και η απαρχή της εισόδου των αριθμητικών μεθόδων στη μηχανική ρευστών. Στο 3ο κεφάλαιο παρατίθεται το θεωρητικό υπόβαθρο πάνω στο οποίο βασίζεται η ανάλυση των διφασικών ροών αλλά και των στροβιλωδών ροών, έννοιες που συναντώνται στα επόμενα κεφάλαια. Η θεωρητική ανάλυση ολοκληρώνεται στο 4ο κεφάλαιο με την εμβάθυνση στην αριθμητική επίλυση των εξισώσεων του 2ου κεφαλαίου. Αρχίζουμε με μία σύντομη παρουσίαση των μεθόδων που αναπτύχθηκαν τις τελευταίες δεκαετίες, επικεντρώνοντας στη μέθοδο που χρησιμοποιούμε στην παρούσα εργασία. Στη συνέχεια παρατίθενται διάφορα μοντέλα προσομοίωσης ροών επικεντρώνοντας στις διφασικές ροές. Επειδή το αντικείμενο αυτό περιλαμβάνει εκτεταμένες μαθηματικές αναλύσεις, μεταφέρθηκαν αρκετά τμήματα στο παράρτημα 1, ώστε το κυρίως κείμενο της εργασίας να παρουσιάζει ομοιογένεια και να διατηρεί την ισορροπία μεταξύ της θεωρητικής ανάλυσης και της πρακτικής – ερευνητικής προσέγγισης που αναφέρθηκε προηγουμένως. Από το 5ο κεφάλαιο ξεκινά η παρουσίαση του καθαρά ερευνητικού έργου και ακολουθείται η λογική και χρονολογική σειρά που αναλύθηκε πριν, με τα τέσσερα στάδια εξέλιξης του «προϊόντος». Σε αυτό λοιπόν το κεφάλαιο αναλύεται η ιδέα –concept- λειτουργίας και η θεωρία πάνω στην οποία βασίζεται και παρουσιάζονται σχέδια από την σχεδίασή του στο AUTOCAD. Εκτός από αυτό όμως, κρίθηκε σκόπιμο να προχωρήσουμε ένα ακόμη βήμα, με τη δόμηση του πλέγματος του πεδίου ροής του διαχωριστή. Ο λόγος είναι καθαρά τεχνικός και έγκειται στην παραπλήσια μορφή των εξής δύο εργασιών: της σχεδίασης στο AUTOCAD και της σχεδίασης της γεωμετρίας στο GAMBIT. Το πλέγμα ακολούθησε αμέσως μετά, για να μην επανερχόμαστε στη χρήση του GAMBIT στο επόμενο κεφάλαιο. Στο 6ο κεφάλαιο γίνεται η αρχική διερεύνηση λειτουργίας του διαχωριστή, με προσομοίωση της ροής με Υπολογιστική Ρευστομηχανική. Ουσιαστικά ολοκληρώνεται η φάση της μοντελοποίησης, με την επεξεργασία του πλέγματος στο FLUENT σε επίπεδο post processing αλλά και με το «στήσιμο» της προσομοίωσης. Τέλος, με το «στήσιμο» του μοντέλου, είμαστε σε θέση να διεξάγουμε την προμοίωση και να συλλέξουμε αποτελέσματα στη συνέχεια. Στην ανάλυση που γίνεται, ακολουθείται η ροή: θεωρητική τεκμηρίωση, υλοποίηση προσομοίωσης ροής και συλλογή – επεξεργασία αποτελεσμάτων. Από τα αποτελέσματα της ανάλυσης προκύπτει ότι η απόδοση αυξάνει με την ταχύτητα εισόδου του μίγματος κάτι που είναι επιστημονικά σωστό – αποδεκτό. Επιτεύχθηκε σύγκλιση της λύσης και παρήχθησαν εικόνες της ροής μέσα στο διαχωριστή που μας δίνουν μια ποιοτική αποτίμηση του πεδίου ροής. Επειδή τα εξαγόμενα της εργασίας εδώ είναι κυρίως εικόνες από το FLUENT που περιγράφουν τη ροή, αυτές παρατίθενται στο παράρτημα 2 για να διευκολυνθεί η ροή της διατριβής. Στο 7ο κεφάλαιο γίνεται η πειραματική – παραμετρική ανάλυση της λειτουργίας του διαχωριστή. Στην ανάλυση αυτή, ακολουθείται και πάλι η ροή: θεωρητική τεκμηρίωση, διεξαγωγή πειράματος και συλλογή – επεξεργασία αποτελεσμάτων. Διερευνήθηκαν οι επιδράσεις των ανεξαρτήτων παραμέτρων συνδυασμένες σαν κλάσμα όγκου αέρα, στην απόδοση του διαχωριστή. Η συσκευή παρουσίασε καλή απόδοση σε συγκεκριμένες περιοχές λειτουργίας κάτι που υποδεικνύει ότι υπάρχει δυναμική για περαιτέρω έρευνα και βελτίωση. Σε αυτό το κεφάλαιο, δίνεται ιδιαίτερη έμφαση στην περιγραφή της πειραματικής εγκατάστασης, η οποία δεν προϋπήρχε, αλλά σχεδιάστηκε και κατασκευάστηκε από το μηδέν για την παρούσα εργασία. Σχέδια από τη φάση της σχεδίασης της εγκατάστασης στο AUTOCAD αλλά και λεπτομέρειές της, παρουσιάζονται αναλυτικά στο παράρτημα 3 μαζί με χαρακτηριστικές εικόνες από τη λειτουργία της. Το πρώτο εξαγόμενο από το πείραμα, είναι η ανάλυση της απόδοσης της συσκευής όπως παρουσιάστηκε στο 7ο κεφάλαιο. Το δεύτερο εξαγόμενο, -εξίσου σημαντικό- είναι η πιστοποίηση της μοντελοποίησης που έγινε στο 6ο κεφάλαιο. Αυτή η διαδικασία λαμβάνει χώρα στο 8ο κεφάλαιο. Το εξαγόμενο εδώ, είναι η τυποποίηση των ρυθμίσεων της υπολογιστικής διαδικασίας με τέτοιο τρόπο, ώστε στον επόμενο κύκλο (loop) πιστοποίησης του προϊόντος να έχουμε περισσότερα και πιο αξιόπιστα αποτελέσματα από το πρώτο στάδιο. Η ακρίβεια της πρόβλεψης με το υπολογιστικό πρόγραμμα, προέκυψε από την καλή αντιστοιχία που επιτεύχθηκε μεταξύ των αποτελεσμάτων της προσομοίωσης με το βέλτιστο μοντέλο, με τα πειραματικά αποτελέσματα τόσο σε σχέση με την παροχή αέρα στην έξοδο, όσο και με την απόδοση διαχωρισμού της συσκευής. Τέλος, στο 9ο κεφάλαιο κατατίθεται σύντομα –και με την εμπειρία που αποκτήθηκε- μία σειρά από δυνατότητες βελτίωσης του διαχωριστή. Είμαστε στο σημείο που αφήνουμε την ακαδημαϊκή ανάλυση του πεδίου ροής που ολοκληρώθηκε με επιτυχία και επικεντρωνόμαστε στην πιο πρακτική πτυχή, της βελτίωσης της απόδοσης της συσκευής. / The growth experienced in recent decades in the extraction of hydrocarbons, led to the adoption of increasingly modern technical products towards the phase separation near the source. The centrifugal separator is the most successful and widespread application. In this range of products, belongs the novel centrifugal separator presented in the present study. The analysis tries to balance between a "strictly scientific" analysis of the turbulent two-phase flow field in the separator and a more "practical" analysis of the performance of this device. The first is necessary because of the originality that has the particular scientific field and the lack of previous studies covering the issue; while the latter is essential because of the "practicality" we expect the final result to have. The problem that the present work is trying to solve, and the literature review that determines the direction to be followed by the research team, are presented in the 1st chapter. To tackle this problem, we follow a standard procedure of an industrial product development in four stages. In the first stage, the concept of function and the theory on which it is based are discussed. In the second stage, an initial investigation of the separator running a CFD simulation is performed. In the third step, an experimental - parametric analysis of the separator’s performance is being held. Finally, the fourth step is a confirmation - verification of CFD simulation. The 2nd chapter, is a theoretical analysis on the development of the governing equations of fluid flow. The equations of motion are developed and the transformation after the turbulence development is explained. At this point it is understood that it is not possible to solve them by analytical methods but with numerical. This is the beginning of the entry of numerical methods in fluid mechanics. The 3rd chapter, outlines the theoretical background upon which the analysis of two-phase flows and of turbulent flows is based. The theoretical analysis is completed in the 4th chapter by deepening in the numerical solution of equations of the second chapter. We begin with a brief presentation of the methods developed in recent decades, focusing on the method used in this work. Several simulation models focusing on the two-phase flow are presented here. Because this topic includes extensive mathematics, several sections were carried in Annex 1 to allow for consistency to the body of work, and maintain a balance between theoretical analysis and practice - research approach mentioned earlier. The presentation of pure research project starts in the 5th chapter, following the logical and chronological sequence that was analyzed before, with four stages in the “product development”. In this chapter, we discussed the concept of function and the theory on which it is based, and presented drawings of the design procedure in AUTOCAD. Apart from this however, it was decided to go a step further by building the grid of the flow field inside separator geometry. The reason is purely technical and is the similarity of the following two tasks: the design in AUTOCAD and geometry and mesh generation in GAMBIT. So the grid was generated in order not to revert to the use of GAMBIT in the next chapter. In the 6th chapter, we investigate the operation of the separator, by conducting an initial CFD simulation. An investigation of air-water separation process has been conducted using Υπολογιστική και Πειραματική Ανάλυση Διφασικής Ροής Αέρα – Νερού σε Φυγοκεντρικούς Διαχωριστές xi numerical simulation. Separator performance has been analyzed in terms of separation efficiency and pressure drop along the water flow path. The results of the analysis show that the performance increases with the inflow speed of the mixture which is scientifically correct - acceptable. Having achieved a converged solution and having produced images of the flow inside the separator we give a qualitative assessment of the flow field. Because the extracted work here is mostly pictures of the FLUENT describing the flow, these are listed in Annex 2 to facilitate the flow of the thesis. In the 7th chapter, the experimental - parametric analysis of the separator’s performance is conducted. New experimental data are reported, on the onsets of gas (air) and liquid (water) entrainment, mass flow rate and separation performance for the two-phase flow through the novel compact phase separator under various conditions. These data series cover the ranges of hmmhmwater/4/1303≤≤− and %90%4≤≤airvf. Influences of the independent parameters airm−0 and waterm−0 combined as air volume fraction vfair on the performance n, were considered. In this chapter, special emphasis is placed on describing the experimental setup, which did not exist, but was designed and built from scratch for this work. Designs from the design stage of the installation in AUTOCAD and its details are presented in detail in Annex 3 along with typical images of the operation. The first outcome of the experiment is the analysis of device’s performance as presented in chapter 7. The second outcome, -equally important- is the certification of modeling technique applied in chapter 6. This process takes place in the 8th chapter. The outcome, is to standardize the configuration of the computational process in such a way, that in the next development loop of the product, we have more and more reliable results than in the first stage. Finally, in chapter 9 (deposited soon after the experience gained), a series of possible improvements in the separator are presented. We leave it to the point that the academic analysis of the flow field was completed successfully and focus on more the practical aspect of improving the device’s performance.

It is well-known that, increasing the flow rate in polymer extrusion, the flow becomes unstable and the smooth extrudate surface becomes wavy and disordered to an increasing degree. In order to investigate the mechanisms responsible for these instabilities we perform a linear stability analysis of the steady extrusion of a viscoelastic fluid flowing through a planar die under creeping flow conditions. We consider the Phan-Thien-Tanner (PTT) model to account for the viscoelasticity of the material. We employ the mixed finite element method combined with an elliptic grid generator to account for the deformable shape of the interface. The generalized eigenvalue problem is solved using Arnoldi’s algorithm. We perform a thorough parametric study in order to determine the effects of all material properties and rheological parameters. We investigate in detail the effect of interfacial tension and the presence of a deformable interface. It is found that the presence of a finite surface tension destabilizes the flow as compared to the case of the stick-slip flow. We recognize two modes which are found to become unstable beyond a critical value of the Weissenberg number and perform an energy analysis to examine the mechanisms responsible for the destabilization of the flow and compare against the mechanisms that have been suggested in the literature. / --

- / -