Coupled Thermodynamic And CFD Approaches Applied To A Supersonic Air Ejector

This paper presents a systematic comparison of ejector performance predictions by a thermodynamic and a CFD model for different operating conditions. The thermodynamic model developed by Galanis and Sorin (2015) assumes the primary flow is always choked, and irreversibilities due to viscous dissipation are taken into account through polytropic efficiencies. The CFD model developed by Croquer et al. (2015) using the software ANSYS Fluent v15.0 has already been validated for supersonic ejectors working with R134a. A standard high Reynolds number k- ω SST turbulence model coupled with the perfect gas law is used to model the turbulent air flows. The dimensions of the ejector were first determined by the thermodynamic model and then used in the CFD model. The thermodynamic model predicts higher entrainment ratios for double choking operation and somewhat different values of the critical and limiting pressure ratios. The CFD model validates the similarity solutions characteristic of ejectors using perfect gases. The present results show in particular that identical inlet pressure and temperature ratios induce the same entrainment ratio as well as the same critical and limiting pressure ratios. Both models confirm also that similar diameter ratios between the primary nozzle throat and the constant area section lead to the same values of the entrainment ratio. Thus, for double-choking operations, the entrainment ratio depends on the inlet pressure and temperature ratios rather than on the individual values of these four properties as is the case for ejectors with real fluids. It also shows that the position of the shock varies linearly with the compression ratio in qualitative agreement with the assumption used in the thermodynamic model. Finally, the main assumptions made to build the thermodynamic model have been checked and discussed a posteriori using the CFD results

Analyse dynamique (CFD) et thermodynamique combinée dans un éjecteur supersonique en présence de gouttelettes

Abstract : This research project has as main objective to study in detail the internal flow features of single-phase supersonic ejectors for refrigeration applications, and the potential effects of injecting droplets on the performance of the device. To this end, a numerical approach is proposed which has been separated into two parts: First, a RANS modelling strategy for supersonic ejectors has been outlined combining the NIST real gas equations database [NIST, 2010] and the k − ω SST turbulence model in its low-Reynolds number formulation. The proposed approach agrees within 5% (resp. 2%) to the experimental entrainment ratio (resp. compression ratio) data of García del Valle et al. [2014], properly captures the main internal flow features and has a reasonable computational cost. This RANS model has been applied in the analysis of a supersonic R134a ejector for refrigeration purposes, showing in particular that the secondary flow is entrained by momentum transfer through the mixing shear layer, that the distance between the primary nozzle exit and the shock-waves in the constant area section varies between 9 and 16 times the primary nozzle exit diameter and that the important axial character of the flow limits mixing of both inlet flows until after the shock train. Furthermore, an exergy analysis through the device shows that the mixing and the oblique shock waves are responsible for between 50% and 70% of the generated losses, the latter might be attenuated through droplet injection in the constant area section. Moreover, it has been shown that drop-in replacement of the working fluid with HFOs R1234yf and R1234ze(E) leads to mild changes in the ejector performance but reduces the HDRC system COP (resp. cooling capacity) in average by 7.1% (resp. 23.3%). Lastly, a comparison of the model predictions with the thermodynamic model of Galanis and Sorin [2016] for an air ejector, shows that as the working fluid approaches the ideal gas behaviour, the flow can be adimensionalized in terms of the secondary inlet temperature and pressure, the motive nozzle throat and the entrainment and compression ratios. In the second part, the influence of droplets has been studied from a local perspective by extending the RANS model to include a discrete phase, which affects the main flow through exchanges of momentum and thermal energy, and from a global perspective by building a thermodynamic model, which predicts the entrainment and limiting compression ratio given a fixed geometry and operating conditions. Both approaches present very good agreement in terms of p, T and M a internal profiles. Results for a supersonic ejector with R134a as baseline working fluid and droplets injected at the constant area section show that the flow structure has perceptible changes only at the highest injection fraction considered 10%, which induces boundary layer detachment, reduces the shock intensity by 8% and diminishes the superheat at the ejector outlet by 15 ◦C. Nonetheless, ejector performance metrics are severely affected as the limiting compression ratio, Elbel efficiency and exergy performance reduce respectively by 5%, 11% and 15%, due mainly to the additional entropy generated through droplet injection and mixing with the main flow.Ce projet de recherche a pour objectif principal d’étudier en détail les caractéristiques de l’écoulement interne dans des éjecteurs supersoniques monophasiques pour des applications en réfrigération, et les effets potentiels de l’injection de gouttelettes sur les performances de l’appareil. A cette fin, une approche numérique est proposée et a été séparée en deux parties. Tout d’abord, une stratégie de modélisation RANS pour les éjecteurs supersoniques a été décrite en combinant la base de données pour les gaz réels NIST [NIST, 2010] et le modèle de turbulence k − ω SST dans sa formulation à bas nombre de Reynolds. L’approche proposée prédit avec un accord d’environ 5% (resp. 2%) le rapport d’entraînement (resp. rapport de compression) avec les données expérimentales de García del Valle et al. [2014]. Il capte également correctement les principales caractéristiques de l’écoulement interne et a un coût de calcul raisonnable. Ce modèle RANS a été appliqué à l’analyse d’un éjecteur supersonique au R134a utilisé à des fins de réfrigération, montrant en particulier que le flux secondaire est entraîné par un transfert d’impulsion à travers la couche de cisaillement, que la position de départ des ondes de choc dans la section constante se situe dans une plage de 9 à 16 fois le diamètre de sortie de la buse primaire et que l’important caractère axial du flux limite le mélange des deux écoulements d’entrée au-delà du train d’ondes de choc. De plus, une analyse exergétique à travers le dispositif montre que le mélange et les ondes de choc obliques sont responsables de 50% et 70% des pertes générées, ces dernières pouvant être atténuées par injection de gouttelettes dans la section à zone constante. De plus, il a été démontré que le remplacement direct du fluide de travail par les HFO R1234yf et R1234ze(E) entraîne de légers changements dans la performance de l’éjecteur mais réduit en moyenne le COP du système HDRC (resp. la capacité de refroidissement) de 7.1% (resp. 23.3%). Enfin, une comparaison des prédictions du modèle avec le modèle thermodynamique de Galanis and Sorin [2016] pour un éjecteur à air montre que lorsque le fluide de travail se rapproche du comportement de gaz idéal, l’écoulement peut être normalisé en fonction de la température et de la pression à l’entrée secondaire, la gorge de la tuyère principale et les rapports d’entraînement et de compression. Dans la seconde partie, l’influence des gouttelettes a été étudiée d’un point de vue local en étendant le modèle RANS à une phase discrète qui affecte le flux principal par des échanges de quantité de mouvement et d’énergie thermique, et d’un point de vue global en construisant un modèle thermodynamique qui prédit l’entraînement et le rapport de compression limitant étant donné une géométrie fixe et les conditions de fonctionnement. Les deux approches présentent un très bon accord en termes de profils internes de p, T et Ma. Les résultats pour un éjecteur supersonique au R134a comme fluide de base, avec des gouttelettes injectées à mi-chemin dans la section de la zone constante, montrent que la structure d’écoulement dans cette région présente des changements perceptibles seulement à la fraction d’injection la plus élevée, 10%, en diminuant l’intensité du choc de 8% et la surchauffe à la sortie de l’éjecteur de 15 ◦C. Néanmoins, la performance de l’éjecteur est sévèrement affectée vu que le rapport de compression, l’efficacité d’Elbel et le performance exergétique sont réduites respectivement de 5%, 11% et 15%, principalement en raison de l’entropie supplémentaire générée par l’injection de gouttelettes et le mélange avec le flux principal

Exergetic Optimisation of Vortex Tubes using a Thermodynamic Model

This article identifies sources of exergy losses in a vortex tube working with air by using a recently developed thermodynamic model and a reference experiment from the literature. Exergetic efficiency considering transiting exergy is used as the efficiency metrics in this work. When both the cold and hot outlets are useful, the exergetic efficiency reaches its maximum value for a cold mass fraction equal to 0.7. Interestingly, up to 45% of the inlet exergy is lost downstream of the vortex tube under this condition because of pressure losses in the cold tube and through measuring instruments. These losses do not contribute to the energy separation mechanism. Inside the vortex tube, the exergy irreversibly is mainly caused by the dissipation of kinetic exergy. The thermodynamic model is also used to identify the working conditions, which maximize the vortex tube efficiency. The efficiency is always at its maximum value when the inlet Mach number is equal to one. The optimum value of the cold outlet diameter, the mass fraction and the cold outlet axial Mach number changes depending on whether thermal exergy from both outlets can be used or not. Increasing the cold outlet pressure increases the exergetic efficiency as well as changing the optimal condition for all variables except the inlet Mach number. At the end, the optimal vortex tube is twice as efficient as the reference vortex tube. Finally, the model is employed to identify the best vortex tubes’ arrangement to maximize the exergetic efficiency for an open cycle with a fixed inlet pressure of six bar. This analysis demonstrates that the best arrangement is a cascade of vortex tubes, where a vortex tube unit with the maximum efficiency is placed first. Two other vortex tubes are two other vortex tubes are placed to recover waste pressure on the cold and hot streams from the first unit

1D modelling of the liquid-gas jet pump performance

Abstract: The performance of a liquid-gas jet pump, which uses a high-velocity liquid flow to compress and entrain a gas flow, can be divided into two modes, on-design and off-design. The present paper investigates the mentioned modes of performance using a 1D model based on the conservation equations of mass, momentum and energy. Comparisons in terms of compression ratio and efficiency between the present model and experimental data show that the 1D model is capable of predicting the behavior of the liquid-gas jet pump for both modes. The effects of the primary flow velocity head, and areas of the mixing throat and diffuser on the performance of the jet pump are also investigated through a sensitivity study.Communication présentée lors du congrès international tenu conjointement par Canadian Society for Mechanical Engineering (CSME) et Computational Fluid Dynamics Society of Canada (CFD Canada), à l’Université de Sherbrooke (Québec), du 28 au 31 mai 2023

Numerical Simulation of the Heat Transfer in a Refrigerated Trailer Equipped with Eutectic Plates for Frozen Food Delivery

The present work reports the Computational Fluid Dynamics simulation and analysis of the heat transfer inside a refrigerated truck trailer equipped with three eutectic plates and fans. The numerical model solves the conjugated heat transfer inside the trailer in 2D using the �� −�� Shear Stress Transport (SST) turbulence model. It has been already favorably validated against the numerical and experimental data of Lafaye de Micheaux el al. (2015) by Croquer et al. (2019). These simulations are used to improve the configuration of the refrigeration system with the eutectic plates as well as to investigate the feasibility of the eutectic plates for the transport of frozen food products under different operating loads and transport temperature requirements. Three eutectic plates having an optimal inter-plate distance of 6 cm to maximize the air flow between the plates (Croquer et al., 2019) are either placed in series on the roof of the trailer or vertically at its back. For both configurations, fans are blowing the air from the eutectic plates to the inside of the trailer and modeled by adding a source term into the momentum equations. During the door opening period, the configuration with the plates placed on the roof of the trailer without the cargo has noticeably lower area-averaged temperature inside the trailer than the configuration with the plates placed on the back of the trailer due to the presence of the circulation zones and the cold plates located near the doorway. However, introduction of the cargo into the simulations eliminates the formation of the circulation zones that prevents the infiltration of the atmospheric air. Also, the configuration with the plates placed on the roof of the trailer allows the atmospheric air to infiltrate earlier, therefore resulting in an overall higher temperature observed in the cargo

Conséquences du réchauffement climatique sur les structures de protection côtière

Parmi les principales conséquences de l'élévation du niveau de la mer associée à la variabilité du climat, se trouvent les impacts sur les structures côtières. Jusqu'à présent, l'attention s'est principalement concentrée sur les zones côtières de faible altitude, car ce sont celles qui subiront de graves inondations. Pourtant, peu d'études ont mis l’emphase sur la performance des structures qui protège ces zones côtières. Dans ce travail, nous avons étudié l'impact de l'élévation du niveau de la mer sur la réponse des structures côtières et, en particulier sur les principaux modes de défaillance des brise-lames. Pour cela, nous avons réalisé une série de tests expérimentaux en laboratoire et des tests numériques simulant différents scénarios d'élévation du niveau de la mer. Nous pouvons conclure que le débit de débordement des vagues et les forces appliquées sur le mur augmentent avec l'élévation du niveau de la mer. Cependant, les tests expérimentaux indiquent que la prédiction d'une relation entre ces deux valeurs n'est pas simple et des programmes expérimentaux minutieux doivent être réalisés à l'avenir

Energy and exergy analysis of possible alternatives to R134a in a vapour compression refrigeration cycle of a water cooler unit

Abstract: In this paper, the energy and exergy performance of a vapour-compression refrigeration cycle of a water dispenser unit has been analyzed theoretically using different refrigerants as possible alternative substitutes to R134a. The selected low ?Global Warming Potential (GWP) ? refrigerants are: HydroFluoroOlefins (HFO) R1234yf and R1234ze(E), HydroCarbons (HC) R290 and R600a, and HydroFluoroCarbon (HFC) R152a. The process was evaluated based on the evaporator and condenser temperatures, which range between -10 and 5?C and between 30 and 45?C, respectively. The theoretical model based on the first and second laws of thermodynamics has been developed using the Matlab environment. The performances of the vapor compression cycle are discussed in terms of coefficient of performance, exergy destruction and exergy efficiency. The results show that the maximum COP achieved is 5.12 and 5.10 for R152a and R600a respectively, the highest total exergy destruction is about 82.82 W for R290, the highest exergy efficiency is about 55.12% for R600a.Communication présentée lors du congrès international tenu conjointement par Canadian Society for Mechanical Engineering (CSME) et Computational Fluid Dynamics Society of Canada (CFD Canada), à l’Université de Sherbrooke (Québec), du 28 au 31 mai 2023

Analyse dynamique (CFD) et thermodynamique combinée dans un éjecteur supersonique en présence de gouttelettes

Abstract : This research project has as main objective to study in detail the internal flow features of single-phase supersonic ejectors for refrigeration applications, and the potential effects of injecting droplets on the performance of the device. To this end, a numerical approach is proposed which has been separated into two parts: First, a RANS modelling strategy for supersonic ejectors has been outlined combining the NIST real gas equations database [NIST, 2010] and the k − ω SST turbulence model in its low-Reynolds number formulation. The proposed approach agrees within 5% (resp. 2%) to the experimental entrainment ratio (resp. compression ratio) data of García del Valle et al. [2014], properly captures the main internal flow features and has a reasonable computational cost. This RANS model has been applied in the analysis of a supersonic R134a ejector for refrigeration purposes, showing in particular that the secondary flow is entrained by momentum transfer through the mixing shear layer, that the distance between the primary nozzle exit and the shock-waves in the constant area section varies between 9 and 16 times the primary nozzle exit diameter and that the important axial character of the flow limits mixing of both inlet flows until after the shock train. Furthermore, an exergy analysis through the device shows that the mixing and the oblique shock waves are responsible for between 50% and 70% of the generated losses, the latter might be attenuated through droplet injection in the constant area section. Moreover, it has been shown that drop-in replacement of the working fluid with HFOs R1234yf and R1234ze(E) leads to mild changes in the ejector performance but reduces the HDRC system COP (resp. cooling capacity) in average by 7.1% (resp. 23.3%). Lastly, a comparison of the model predictions with the thermodynamic model of Galanis and Sorin [2016] for an air ejector, shows that as the working fluid approaches the ideal gas behaviour, the flow can be adimensionalized in terms of the secondary inlet temperature and pressure, the motive nozzle throat and the entrainment and compression ratios. In the second part, the influence of droplets has been studied from a local perspective by extending the RANS model to include a discrete phase, which affects the main flow through exchanges of momentum and thermal energy, and from a global perspective by building a thermodynamic model, which predicts the entrainment and limiting compression ratio given a fixed geometry and operating conditions. Both approaches present very good agreement in terms of p, T and M a internal profiles. Results for a supersonic ejector with R134a as baseline working fluid and droplets injected at the constant area section show that the flow structure has perceptible changes only at the highest injection fraction considered 10%, which induces boundary layer detachment, reduces the shock intensity by 8% and diminishes the superheat at the ejector outlet by 15 ◦C. Nonetheless, ejector performance metrics are severely affected as the limiting compression ratio, Elbel efficiency and exergy performance reduce respectively by 5%, 11% and 15%, due mainly to the additional entropy generated through droplet injection and mixing with the main flow.Ce projet de recherche a pour objectif principal d’étudier en détail les caractéristiques de l’écoulement interne dans des éjecteurs supersoniques monophasiques pour des applications en réfrigération, et les effets potentiels de l’injection de gouttelettes sur les performances de l’appareil. A cette fin, une approche numérique est proposée et a été séparée en deux parties. Tout d’abord, une stratégie de modélisation RANS pour les éjecteurs supersoniques a été décrite en combinant la base de données pour les gaz réels NIST [NIST, 2010] et le modèle de turbulence k − ω SST dans sa formulation à bas nombre de Reynolds. L’approche proposée prédit avec un accord d’environ 5% (resp. 2%) le rapport d’entraînement (resp. rapport de compression) avec les données expérimentales de García del Valle et al. [2014]. Il capte également correctement les principales caractéristiques de l’écoulement interne et a un coût de calcul raisonnable. Ce modèle RANS a été appliqué à l’analyse d’un éjecteur supersonique au R134a utilisé à des fins de réfrigération, montrant en particulier que le flux secondaire est entraîné par un transfert d’impulsion à travers la couche de cisaillement, que la position de départ des ondes de choc dans la section constante se situe dans une plage de 9 à 16 fois le diamètre de sortie de la buse primaire et que l’important caractère axial du flux limite le mélange des deux écoulements d’entrée au-delà du train d’ondes de choc. De plus, une analyse exergétique à travers le dispositif montre que le mélange et les ondes de choc obliques sont responsables de 50% et 70% des pertes générées, ces dernières pouvant être atténuées par injection de gouttelettes dans la section à zone constante. De plus, il a été démontré que le remplacement direct du fluide de travail par les HFO R1234yf et R1234ze(E) entraîne de légers changements dans la performance de l’éjecteur mais réduit en moyenne le COP du système HDRC (resp. la capacité de refroidissement) de 7.1% (resp. 23.3%). Enfin, une comparaison des prédictions du modèle avec le modèle thermodynamique de Galanis and Sorin [2016] pour un éjecteur à air montre que lorsque le fluide de travail se rapproche du comportement de gaz idéal, l’écoulement peut être normalisé en fonction de la température et de la pression à l’entrée secondaire, la gorge de la tuyère principale et les rapports d’entraînement et de compression. Dans la seconde partie, l’influence des gouttelettes a été étudiée d’un point de vue local en étendant le modèle RANS à une phase discrète qui affecte le flux principal par des échanges de quantité de mouvement et d’énergie thermique, et d’un point de vue global en construisant un modèle thermodynamique qui prédit l’entraînement et le rapport de compression limitant étant donné une géométrie fixe et les conditions de fonctionnement. Les deux approches présentent un très bon accord en termes de profils internes de p, T et Ma. Les résultats pour un éjecteur supersonique au R134a comme fluide de base, avec des gouttelettes injectées à mi-chemin dans la section de la zone constante, montrent que la structure d’écoulement dans cette région présente des changements perceptibles seulement à la fraction d’injection la plus élevée, 10%, en diminuant l’intensité du choc de 8% et la surchauffe à la sortie de l’éjecteur de 15 ◦C. Néanmoins, la performance de l’éjecteur est sévèrement affectée vu que le rapport de compression, l’efficacité d’Elbel et le performance exergétique sont réduites respectivement de 5%, 11% et 15%, principalement en raison de l’entropie supplémentaire générée par l’injection de gouttelettes et le mélange avec le flux principal

Thermodynamic Modelling of Supersonic Gas Ejector with Droplets

This study presents a thermodynamic model for determining the entrainment ratio and double choke limiting pressure of supersonic ejectors within the context of heat driven refrigeration cycles, with and without droplet injection, at the constant area section of the device. Input data include the inlet operating conditions and key geometry parameters (primary throat, mixing section and diffuser outlet diameter), whereas output information includes the ejector entrainment ratio, maximum double choke compression ratio, ejector efficiency, exergy efficiency and exergy destruction index. In single-phase operation, the ejector entrainment ratio and double choke limiting pressure are determined with a mean accuracy of 18 % and 2.5 % , respectively. In two-phase operation, the choked mass flow rate across convergent-divergent nozzles is estimated with a deviation of 10 % . An analysis on the effect of droplet injection confirms the hypothesis that droplet injection reduces by 8 % the pressure and Mach number jumps associated with shock waves occuring at the end of the constant area section. Nonetheless, other factors such as the mixing of the droplets with the main flow are introduced, resulting in an overall reduction by 11 % of the ejector efficiency and by 15 % of the exergy efficiency