Kinetic modeling of polyatomic heat and mass transfer in rectangular microchannels

The present study aims at estimating the heat and the mass transfer coefficients in the case of the polyatomic gas flows through long rectangular microchannels driven by small and large pressure (Poiseuille flow) and temperature (Thermal creep flow) drops. The heat and mass transfer coefficients are presented for all gas flow regimes, from free molecular up to hydrodynamic ones, and for channels with different aspect ratios as well as for various values of translational and rotational Eucken factors. The applied values of the Eucken factors were extracted based on the Rayleigh-Brillouin experiments and the kinetic theory of gases. The numerical study has been performed on the basis of a kinetic model for linear and non-linear gas molecules considering the translational and rotational degrees of freedom. The solution of the obtained system of the kinetic equations is implemented on the Graphics Processing Units (GPUs), allowing the reduction of the computational time by two orders of magnitude. The results show that the Poiseuille mass transfer coefficient is not affected by the internal degrees of freedom and the non-dependence of the previous observed deviations with the experimental data on the molecular nature of the gas molecules is confirmed. However, the study shows that the deviation between monatomic and polyatomic values of the mass transfer coefficient in the thermal creep flow is increased as the gas rarefaction is decreased, and for several polyatomic gases met in practical applications in the temperature range from 300 to 900 K might reach 15%. In addition, the effect of the internal degrees of freedom on the heat transfer coefficient is found to be rather significant. The polyatomic heat transfer coefficients are obtained essentially higher than the monatomic ones, with the maximum difference reaching about 44% and 67% for linear and non-linear gas molecules. In view of the large differences between monatomic and polyatomic gases, the present results may be useful in the design of technological devices in which the thermal creep phenomenon plays a dominant role

Preliminary analysis of alternative divertors for DEMO

A physics and engineering analysis of alternative divertor configurations is carried out by examining benefits and problems by comparing the baseline single null solution with a Snowflake, an X- and a Super-X divertor. It is observed that alternative configurations can provide margin and resilience against large power fluctuations, but their engineering has intrinsic difficulties, especially in the balance between structural solidity and accessibility of the components and when the specific poloidal field coil positioning poses further constraints. A hybrid between the X- and Super-X divertor is proposed as a possible solution to the integration challenge

Flows through channels of various cross sections in the whole range of the knudsen number using deterministic and stochastic approaches

During the last decade research in rarefied gas dynamics has attracted a lot of attention. This refreshed interest is due to applications in the emerging field of nano- and micro-fluidics, as well as to the more traditional fields of vacuum technology under low, medium and high vacuum conditions and high altitude aerodynamics. Some of these applications may include mass flow controllers in gas metering, micro propulsion nozzles in high altitude and space gas dynamics, pipe networks, pumps and distillation towers in vacuum systems, membranes and porous media in filtering, fabrication processes in microelectronics and gaseous devices in micro-electromechanical systems. It is noted that gas flows through channels of various cross sections may be basic elements of all these equipments. The gas rarefaction is specified by the Knudsen number (Kn), which is defined as the ratio of the mean free path over a characteristic macroscopic length of the problem. In general, when the flow is considered as far from local equilibrium, then the well known hydrodynamic equations cannot be applied since the continuum assumption and the associated constitutive laws are not valid anymore. In these cases two main approaches can be implemented in order to solve problems in the whole range of the Knudsen number. The first approach is based on the kinetic theory as expressed by the Boltzmann equation or its associated kinetic models. The main unknown is the distribution function, which obeys a kinetic equation, while all the macroscopic quantities of practical interest are obtained by the moments of the distribution function. In many occasions kinetic model equations are very reliable and can be used as an alternative to the Boltzmann equation, producing accurate results with much less computational effort. Typical examples are the BGK and Shakhov models for isothermal and non-isothermal flows respectively. The second approach is the Direct Simulation Monte Carlo (DSMC) method, which actually circumvents the direct solution of a kinetic equation and over the years it has been found to be very efficient in solving high speed flows. It is noted that the DSMC method recovers the solution of Boltzmann equation in the limit of vanishing discretization. Both approaches may be used in the case of rarefied flows through channels. When the flow is very slow or when the channel is adequate long, then the flow is fully developed (no variation in the flow direction) and can be handled by applying linear kinetic theory. When the channel is short, end effects cannot be neglected and the flow becomes non linear. Then, the most efficient approach is the DSMC method. Within the above framework, the present Ph.D. thesis is focused on the investigation of rarefied gas flow through long and short channels of various cross sections (cylindrical, annular cylindrical, orthogonal, triangular and trapezoidal), using kinetic equations and he DSMC method. In the former approach, the linear BGK kinetic model is solved using the Discrete Velocity Method (DVM) is applied, which has been extensively used over the years with considerable success in the solution of the Boltzmann equation or of reliable kinetic models. It is characterized by the deterministic discretization of the kinetic equations in the physical and molecular velocity spaces. In the latter approach, an upgraded version of the DSMC algorithm based on the Non Time Counter Scheme is used. This is a probabilistic scheme. The fundamental idea of the method is to track a large number of statistically representative particles. The primary approximation of DSMC is to uncouple the molecular motion and the intermolecular collisions over small time intervals. Particle motions are modeled deterministically, while the collisions are treated statistically. The core of the DSMC algorithm consists of four primary processes: free particle motion, index and cross-reference of the particles, simulation of collisions, and sampling of the flow field.Την τελευταία δεκαετία έντονο ερευνητικό ενδιαφέρον έχει εκδηλωθεί στο πεδίο της δυναμικής των αραιοποιημένων αερίων. Η στροφή της επιστημονικής κοινότητας προς την κατεύθυνση αυτή οφείλεται κυρίως στην εμφάνιση νέων τεχνολογιών και εφαρμογών που περιλαμβάνουν μεταξύ άλλων την τεχνολογία χαμηλού, μέσου και υψηλού κενού, την αεροδιαστημική μεγάλων υψομέτρων και τέλος την νανοτεχνολογία. Η αραιοποίηση ενός αερίου καθορίζεται από τον αριθμό Knudsen (Kn), ο οποίος ορίζεται ως το πηλίκο της μέσης ελεύθερης διαδρομής των μορίων του αερίου προς ένα χαρακτηριστικό μήκος του εκάστοτε προβλήματος. Εν γένει, σε ροές μακριά από την θερμοδυναμική ισορροπία, η θεωρία του συνεχούς μέσου καταρρέει και οι γνωστές υδροδυναμικές και καταστατικές εξισώσεις παύουν να ισχύουν. Στις περιπτώσεις αυτές η αντιμετώπιση του προβλήματος βασίζεται στην Κινητική Θεωρία των αερίων, η οποία περιλαμβάνει δύο βασικές προσεγγίσεις για την επίλυση σε όλο το εύρος του αριιθμού Knudsen. Η πρώτη προσέγγιση βασίζεται στην κινητική θεωρία, όπως αυτή εκφράζεται μέσω της εξίσωσης Boltzmann ή των αντίστοιχων κινητικών μοντέλων. Ο βασικός άγνωστος είναι η συνάρτηση κατανομής, ενώ οι μακροσκοπικές ποσότητες που ενδιαφέρουν προκύπτουν από διάφορες ροπές της συνάρτησης κατανομής. Τονίζεται ότι σε αρκετές εφαρμογές τα κινητικά μοντέλα αποδεικνύονται αρκετά αξιόπιστα και είναι δυνατόν να χρησιμοποιηθούν ως εναλλακτική επιλογή της εξίσωσης Boltzmann, έχοντας ακριβή αποτελέσματα με λιγότερο υπολογιστικό κόστος σε όλο το εύρος του αριθμού Kn. Χαρακτηριστικά παραδείγματα κινητικών μοντέλων είναι το BGK και το Shakhov για ισοθερμοκρασιακές και μη ισοθερμοκρασιακές ροές αντίστοιχα. Η δεύτερη προσέγγιση είναι η μέθοδος Direct Simulation Monte Carlo (DSMC), η οποία παρακάμπτει την επίλυση της εξίσωσης Boltzmann και έχει αποδειχθεί ότι είναι πολύ αποτελεσματική σε περιπτώσεις γρήγορων ροών. Οι δύο παραπάνω προσεγγίσεις χρησιμοποιούνται για την επίλυση αραιοποιημένων ροών διαμέσου αγωγών. Στην περίπτωση αργών ροών ή όταν το μήκος των αγωγών είναι αρκετά μεγάλο, τότε η ροή είναι πλήρως ανεπτυγμένη και περιγράφεται από τη γραμμική κινητική θεωρία. Στην περίπτωση αγωγών πεπερασμένου μήκους, τα φαινόμενα των άκρων δεν μπορούν να αγνοηθούν και η ροή είναι μη γραμμική. Τότε για την επίλυση τέτοιων ροών η πιο αποδοτική μέθοδος είναι η μέθοδος DSMC. Στο πλαίσιο αυτό η παρούσα διδακτορική διατριβή εστιάζεται στην μελέτη ροών αραιοποιημένων αερίων σε αγωγούς μικρού και μεγάλου μήκους με ποικίλες διατομές (κυλινδρική, δακτυλοειδή, ορθογώνια, τριγωνική και τραπεζοειδή) εφαρμόζοντας ντετερμινιστικές και στοχαστικές μεθοδολογίες αντίστοιχα. Η πρώτη προσέγγιση βασίζεται στην επίλυση του γραμμικού κινητικού μοντέλου BGK με την χρήση της μεθόδου των διακριτών ταχυτήτων (Discrete Velocity Method - DVM), η οποία έχει εφαρμοστεί με επιτυχία στην επίλυση της εξίσωσης Boltzmann ή των αντίστοιχων κινητικών μοντέλων. Η μέθοδος DVM χαρακτηρίζεται από την ντετερμινιστική διακριτοποίηση των κινητικών εξισώσεων στον φυσικό χώρο και στον χώρο των μοριακών ταχυτήτων. Στην τελευταία προσέγγιση γίνεται χρήση ενός αναβαθμισμένου αλγορίθμου της μεθόδου DSMC, ο οποίος βασίζεται στο μοντέλο Non Time Counter για την μοντελοποίηση των μοριακών συγκρούσεων. Η μέθοδος DSMC αποτελεί μια στοχαστική μεθοδολογία. Η θεμελιώδης ιδέα της μεθόδου έγκειται στην μοντελοποίηση ενός στατιστικά μεγάλου αριθμού σωματιδίων. Η βασική προσέγγιση της μεθόδου DSMC είναι ο διαχωρισμός της ελεύθερης κίνησης των σωματιδίων μοντελοποιείται ντετερμινιστικά, ενώ οι μοριακές συγκρούσεις μοντελοποιούνται στοχαστικά. Ο πυρήνας του αλγορίθμου της μεθόδου DSMC αποτελείται από τέσσερεις βασικές διαδικασίες: α) την ελεύθερη κίνηση των σωματιδίων, β) την μετονομασία των σωματιδίων, γ) την μοντελοποίηση των μοριακών συγκρούσεων και δ) τον υπολογισμό των μακροσκοπικών ποσοτήτων. Τα τέσσερα παραπάνω βήματα πραγματοποιούνται ανεξάρτητα σε κάθε χρονικό βήμα