Swirl-Bypass Nozzle for CO 2 Two-Phase Ejectors: Numerical Design Exploration

: In this work, a novel ejector design concept of a swirl-bypass nozzle is proposed to improve off-design performance of CO 2 two-phase ejectors. The swirl-bypass nozzle allows part of the ﬂow to bypass into the ejector mixing chamber to generate swirl. The design of such a device is investigated using a 3D multiphase CFD model. An extensive experimental test campaign is conducted to validate the baseline homogeneous equilibrium CFD model. The model’s prediction motive mass ﬂow rate within 2–12% error and suction mass ﬂow rate was predicted with 3–50% error. Based on the tested ejector geometry, simulations of different ejector swirl-bypass inlets are conducted. The results show that, for the current design, total entrainment of the ejector is reduced by 2–20% with the swirl-bypass inlet. The axial position of the bypass inlet plays a primary role in the bypass inlet ﬂow rate, and, consequently, in suction ﬂow reduction. This is found to be due to the bypass ﬂow blocking off the suction mass ﬂow rate, which has a net negative impact on performance. Finally, several design improvements to improve future designs are proposed.


Introduction 1.Literature Review
Environmentally friendly heating and cooling solutions have seen a large growth in recent years.With the ratification of the EU F-gas legislation and the Kigali amendment to the Montreal Protocol [1], increased efforts have been put toward the development of efficient and environmentally friendly cooling and heat-pumping systems.Due to the rapid transition to environmentally friendly refrigerants, much research has been put toward finding efficient system solutions with natural working fluids, such as ammonia, water, hydrocarbons, and CO 2 .
However, the high operating pressure of the CO 2 transcritical cycle is associated with increased expansion losses.These losses can be partially recovered by implementing work-recovery devices such as two-phase ejectors.Ejectors work by using the expansion energy of the high-pressure flow (motive flow) and expand the flow through a super-sonic converging-diverging nozzle.The high-velocity motive flow is then mixed with a secondary flow (suction flow) and the mixture is repressurized by slowing the flow down in a diffuser.This effectively pumps the suction flow to a higher pressure without additional work, improving system efficiency.Ejector performance is, therefore, commonly evaluated using the ratio of pumped suction mass flow rate to motive mass flow rate (Equation ( 1)), referred to as the entrainment ratio of the ejector.ω = MFR suction /MFR motive (1) Two-phase R744 (CO 2 ) ejectors are today implemented in several cooling and heatpump applications [18][19][20][21][22][23][24][25][26][27][28], and can improve system COP by 10-30% [19,29].
The promising benefits of using ejector-supported R744 systems have promoted much research into ejector design and operation.Due to the complex fluid dynamics within these components, much experimental and numerical research has been devoted to better the understanding of multiphase ejectors and their design.Accurate numerical models are of high value for R744 ejector design due to several reasons, such as high sensitivity to small changes in its many design parameters, the high cost of the experimental testing of ejectors, the large differences in optimal design for different applications, and the possibility of performing low-cost exploratory investigations into their design.This has prompted research into advanced models for R744 ejectors, such as the 0D and 1D models [30][31][32][33][34][35], which have been applied for system-level calculations and ejector design.A limitation of numerical modeling using thermodynamic or 1D approaches is their inability to generally predict phenomena such as flow separation or shock-wave patterns.This problem becomes important when exploring novel ejector designs.Furthermore, new geometric ejector designs that are not within conventional ejector design space are not possible to investigate with such approaches.In recent years, emphasis has been put toward developing accurate and fast multiphase computational fluid dynamics (CFD) models for R744 ejectors, where the full 2D or 3D ejector physics are resolved, allowing for more exploratory investigations of ejector designs.Accurate CFD modeling has the added benefit of allowing detailed insights into ejector physics as all flow variables are available.However, these models are, in general, more complex and more computationally costly.
The most prominent model used in R744 two-phase ejectors modeling is the homogeneous equilibrium models (HEM) [36][37][38][39][40][41][42][43][44].In general, for high-pressure and temperature motive conditions, the HEM provides good prediction accuracy at a low computational cost.However, as the motive pressure and temperature are reduced, non-equilibrium effects are no longer negligible, and delayed phase transition must be considered.Several models have been presented in the literature to model the transition to equilibrium, such as homogeneous relaxation models [40,45,46], and mixture models [47][48][49][50].These models treat transition by considering a delayed time scale to reach equilibrium, or by modeling the phase change mechanisms directly.Non-equilibrium models have better predictive capabilities than equilibrium-based approaches for motive operating conditions in the "off-critical" region (i.e., T m 30 • C, P m 75 bar) [51].Accurate prediction of R744 ejector flow is still an ongoing research area.Significant uncertainties are also associated with the prediction of the secondary ejector flow.This is considered a consequence of the lack of detailed understanding of the suction entrainment process, its dependency on accurate motive flow conditions, and the complexity of the multiphase turbulent mixing process in the ejector.For detailed reviews of the current status of R744 two-phase ejector modeling, see Ringstad et al. [51] and Nowak et al. [52].Extensive experimental tests are needed to validate and verify the numerical models.Detailed experimental observations of R744 ejectors are limited, even though some work has been carried out on visualization and local measurements [53][54][55], additional experimental results with varied and open-access ejector geometries are needed.
A primary reason for the large interest in ejector modeling is the necessity for a detailed understanding of ejector design and its impact on ejector performance under varied operating conditions, which has been highlighted as a key focus area for ejector developments [19,51,52].Several studies on R744 ejector design have been conducted to better understand the interdependence of different design parameters [29,44,56].Improved ejector-design processes have been explored using numerical CFD modeling in combination with optimization techniques.Palacz et al. [38,39] combined an evolutionary optimization algorithm with CFD modeling to optimize ejector design.Ringstad et al. [57] presented a machine-learning-based design methodology for the ejector design and performance mapping of R744 ejectors.
An efficient ejector design is highly dependent on the system operating conditions.When operated at off-design operating conditions, a passive (no active control) ejector system can dramatically reduce system COP.Nakagawa et al. [29] reported between a 34 to 82% drop in system COP for an ejector system operated in off-design conditions.Similarly, Lucas and Koehler [58] reported a 10 to 17% drop in system performance of an ejector-based system at off-design conditions.Therefore, ejector control is important for system performance.Different capacity control strategies using R744 ejectors have been presented in the literature; for a detailed review, see Gullo et al. [59].The main technologies for this application are the multiblock ejector [60], motive swirl control [61], adjustable needle control [33,62], and pulse-width modulation [63].
Alternatively, ejector performance can be improved by active control of the suction flow.One such design is the suction-bypass design, where a controllable secondary suction inlet downstream from along the mixing section is implemented for off-design performance improvement.This inlet splits the suction flow within the ejector so that parts of the flow bypass parts of the mixing section.This allows for the suction flow to enter the ejector in a low-pressure zone under off-design conditions.The bypass concept improves ejector performance by reducing entropy production in supersonic shocks and allows for more optimal suction inflow conditions.In such designs, the entrainment ratio is defined by the total mass flow rate pumped by the recovered work from the motive flow, Equation (2): The goal of such a design is to increase the total pumped mass flow beyond what can be achieved with only the suction port.The bypass concept for ejectors has been investigated with various working fluids, with steam [64,65], methane [66,67], air [68,69], and R744 [70,71].Tang et al. [65] investigated the bypass concept for steam ejectors using numerical CFD calculations.The study identified several low-pressure zones that could be utilized, and close to a 4% improvement in entrainment ratio was obtained by implementing the bypass nozzle.The bypass concept was implemented for pressure regulation in the steam ejector [65], and findings indicate an increase in entrainment ratio of up to 26% using this strategy.Chen et al. [68] presented a CFD investigation into the bypass ejector concept for air ejectors, finding that up to a 10% increased entrainment ratio could be generated.Later, Chen et al. [69] presented additional numerical results with different geometries for the bypass inlet.They found that significant improvements by up to 32.8% in entrainment ratio could be made by optimization of the bypass geometry parameters.Chen et al. [66] investigated an ejector study with methane as the working fluid where an ejector with an adjustable needle and bypass inlet was implemented, achieving up to 74.5% improved entrainment.Chen et al. [67] followed this with a detailed investigation of the methane ejector with two bypass inlets.With one bypass inlet, up to a 34.8% improvement was obtained in comparison to without a bypass.By the implementation of two bypass ports, a larger improvement of up to 48.9% was achieved.While the findings with single-phase ejectors are valuable for comparison, much of the knowledge is not directly transferable to two-phase ejectors such as those predominantly implemented in R744 cycles.
The bypass concept for R744 two-phase ejectors was first proposed and explored numerically by Bodys et al. [70].They explored different bypass geometries using a 2D CFD model.The exploration was performed with six bypass locations in the ejector diffuser and with two bypass inlet angles at different pressure lift conditions.Their findings were that the bypass inlet could improve suction mass flow rate at low-pressure lifts by up to 37%.Furthermore, they found that the optimal location for the bypass inlet was at a distance of 40% of the mixing chamber into the diffuser.Recently, Bodys et al. [71] followed up this study with an experimental investigation of a bypass suction inlet.In this work, the authors proved the applicability of the suction-bypass concept for ejector efficiency improvement and proposed a control strategy for such solutions.The resulting efficiency improvement was as high as 37% at low-pressure lifts when the bypass inlet was introduced.Furthermore, a full 3D CFD investigation was conducted to gain a better understanding of the ejector flow physics.The authors mention the potential for introducing a swirl component to the model for further improvements of the bypass concept.
The concept of adding the swirl component to the motive and suction flow at the ejector inlets has been proposed in previous works to intensify the mixing process in the mixing section of the ejector [72,73].The idea is that the surface area of the interface between the suction and motive flow is increased and more momentum is transferred.Additionally, the centrifugal force can encourage faster spreading and, therefore, mixing of the central flow core.Bodys et al. [73] investigated the effects of the swirl-flow component at the ejector motive and suction inlets of a two-phase R744 ejector.They concluded that some improvements could be gained in the entrainment ratio by the implementation of a motive inlet swirl component.However, the suction flow swirl component only limited the suction flow rate through the suction nozzle and had a negative impact on the entrainment ratio.In the mentioned study, as the swirl component was added before the suction inlet, the impact of the swirl component in the mixing chamber was reduced.

Swirl-Bypass Concept
To the best of the authors' knowledge, the concept of ejector improvement through combined bypass flow and swirl generation has not been considered.In this paper, such a solution is presented and discussed in terms of its design and flow characteristics.The novel swirl-bypass concept is presented in this paper to potentially cover two applications: (A) improving entrainment ratio at off-design operating conditions and (B) serving as an additive motive-capacity control technique.These applications are illustrated in Figure 1A,B.In this work, concept (A) is investigated.The novel idea is to add an additional angled inlet to the ejector downstream of the motive nozzle to generate swirl flow.At off-design conditions, the inlet can be used to bypass a part of the suction flow, similar to the standard bypass inlet, and the flow can be used to generate swirl.The swirl-inducing suction flow could then potentially improve mixing efficiency, indirectly and directly improving the entrainment ratio.This corresponds to concept (A).
Alternatively, the swirl-bypass inlet can be used to control the motive mass flow rate.Similar to the motive swirl control mechanism presented by Zhu and Elbel [61], the high-pressure motive flow can be split into two streams.One stream enters the motive nozzle normally, while, unlike the motive swirl control discussed by Zhu and Elbel [61], the second stream enters the swirl-bypass entrance inside the ejector mixing chamber.The second motive flow is then used to generate swirling flow and can be controlled by a valve upstream of the bypass inlet, corresponding with concept B. Such a design allows for the control of the motive flow not only to reduce the motive flow below the choked design value but also to increase the motive mass flow by bypassing the motive nozzle to generate swirl.This concept will perform better than a parallel expansion valve as long as the swirl flow has any positive impact on the flow.

Knowledge Gap
In this work, the novel concept of a swirl-inducing bypass inlet is explored numerically using 3D CFD simulation.As limited comparable ejector designs are available, the design space of such a bypass inlet will be tested with varied bypass geometries.Different inclination angles and locations of the bypass chamber are explored.The flow structure and its implications on ejector operation and design are discussed in terms of swirl and mixing.The results indicate that the design of this type of ejector is complex, and must be considered in conjunction with the rest of the ejector design.This study will, in this way, enable future designs of a swirl-bypass inlet and several design considerations and suggestions are presented.Limited knowledge regarding the influence of swirl on two-phase ejectors is available; however, the potential of utilizing swirl in ejectors has been highlighted by several authors in the literature [61,71].This article further fills this knowledge gap by supplying detailed simulation results and analysis of the influence and decay of swirl in the ejector.Furthermore, the baseline HEM CFD model without a bypass is validated against new experimental results with an ejector geometry designed for high entrainment at low-pressure lift.

Multiphase Model
The CFD model used in this work is a homogeneous equilibrium model (HEM) based on the formulation by Smolka et al. [36].The HEM assumes that both phases are at full mechanical, thermal, and thermodynamic equilibrium; hence, the flow can be described using a single velocity-( u), temperature-(T), and pressure-field (P).The HEM is, therefore, classified as a pseudo-fluid or single-fluid approach where a single set of transport equations, Equations (3)-( 5), are solved.Here, the energy equation, Equation (5), is reformulated to an equivalent enthalpy formulation.
Here, the subscript mix indicates the pseudo-fluid mixture properties.ρ, p, h, q refer to the density, pressure, enthalpy, and heat flux, respectively.The effective stress tensor τ ij,eff is the laminar (Newtonian) and turbulent stress tensors combined, τ eff = τ laminar + τ turb .In this equation, h is the specific enthalpy, and Γ eff is the effective diffusion coefficient.The terms Ṡh1,2,3 describe the mechanical energy, the irreversible dissipation of the kinetic energy variations, and the dissipation of the turbulent kinetic energy, respectively, explained in detail in the work by Smolka et al. [36].
The enthalpy and pressure can then be used to uniquely identify the thermodynamic equilibrium state, and, from this, the thermodynamic properties can be calculated.Properties are typically divided into thermophysical (Equation ( 6)) and transport (Equation ( 7)) properties.
where ρ, α, µ, λ, c p , s are the pseudo-fluid density, vapor volume fraction, kinematic viscosity, thermal conductivity, heat capacity, and entropy, respectively.Mixture properties are evaluated using mass-and volume-weighted averages.
The properties of liquid and gas are evaluated by the pressure and enthalpy, interpolated from a look-up table.The CoolProp library [74] for R744 is based on the Span-Wagner equation of state, which is considered the most accurate equation of state (EoS) for CO 2 and is widely used for R744 ejector simulations [51].Here, CoolProp is used to generate the look-up tables that are imported in Fluent using UDFs.
The HEM model for R744 two-phase ejectors is considered accurate for motive flow inlet conditions at supercritical pressure and temperatures, referred to as "near-critical" conditions [51].In the present study, the operating conditions will primarily concern high pressure and temperature operating conditions, where the assumption of equilibrium flow is considered fair [46].Compared to other approaches presented in the literature, the HEM model is preferable in robustness and stability.Numerical stability is critical for modeling the complex geometry considered in this work.The HEM model has been extensively validated for R744 two-phase ejectors [51] in terms of mass flow rate predictions.Still, the prediction uncertainty of local variables, such as pressures and velocities, is not wellunderstood due to the lacking experimental validation data.Due to the high accuracy in the considered domain and superior stability characteristics, the HEM model is chosen for this work.

Turbulence Model
Turbulence modeling is a key feature of accurate CFD predictions.The mixing process inside the ejector is, to a large extent, governed by the turbulent viscosity, ν t , predicted by the turbulence model.For R744 two-phase ejectors, it is common to use pseudo-fluid twoequation turbulence models, such as the k − and k − ω models.The turbulent viscosity, ν t , can then be calculated based on two turbulence parameters.The models involve solving the set of transport equations for these turbulent properties.As studied by Bodys et al. [50], different turbulence models can yield significantly different predictions of turbulent flow characteristics.An initial study was conducted to compare the k − and k − ω models for the swirl bypass inlet geometry, discussed in Section 4.4.Due to the complexity of the geometry and sharp gradients near the bypass inlet, a fully wall-resolved simulation with turbulent mesh resolution of y + ≈ 1 was not achievable and unstable simulations would occur.Instead, the mesh was resolved to y + ≈ 30+.In conclusion, the k − ω SST model is used due to its better accuracy.

Numerical Approach
The CFD problem is solved using a 3D steady-state pressure-based coupled algorithm in ANSYS Fluent 2020 R2.The numerical schemes used were the second-order upwind scheme for the momentum, density, enthalpy (UDS), and turbulence transport equations.Second-order schemes were chosen to reduce the influence of numerical viscosity, while higher order discretizations were disregarded due to numerical stability concerns near the swirl bypass inlet.Higher order accuracy will improve gradient estimation, which has an impact on shock resolution.The pressure was calculated using the PRESTO!scheme, which is the preferred pressure discretization for high swirl flows [75].Gradients were evaluated with the least-square cell gradient calculation, which is associated with comparable or higher accuracy than the Green Gauss node-and cell-based approaches, respectively, for regions of high mesh non-orthogonality [75].This is primarily a concern in the vicinity of the swirl-bypass inlet, as the baseline mesh has otherwise high orthogonal quality.The boundary conditions were specified as pressure inlets and outlets.Constant pressure and enthalpy were specified at the inlets.These boundary conditions are standard for CFD modeling of two-phase R744 ejectors [51], and ensure that the modeling problem is well-posed.These boundary conditions are set according to the experimentally measured values; thus, measurement uncertainty limits the accuracy of the model prediction.The turbulence boundary conditions were set to 5% turbulence intensity and a turbulent viscosity ratio of 5.The performance impact of turbulent boundary conditions is commonly found to be negligible [50] due to the turbulence generated in the ejector typically being more significant than the inlet turbulence.

Calculation of Ejector Parameters
Different ejector parameters will be investigated in this study to better understand the ejector flow structure.To evaluate the effects on ejector mixing-layer thickness, the vorticity thickness, δ w , is used [76]: The calculation is based on defining the thickness by the largest gradient in the mixing layer compared to the velocity difference across the boundary.The velocity difference, ∆U ∞ , is calculated from: where U motive is calculated based on the maximum velocity of the motive flow, and U suction is calculated based on the velocity in the section between the upper boundary layer and the motive core.The characteristic velocity was, for this region, calculated based on the velocity at r s = 0.9 • R wall .

Swirl Injection Geometry
The swirl-bypass inlet is an additional inlet entering the mixing chamber of the ejector.As opposed to the standard bypass inlet, the novel swirl-bypass nozzle is angled with a tangential component to the ejector.The effect of such a design is investigated in detail in this study.Several studies have shown that positive efficiency effects can be gained by increasing mixing by the use of swirl flow.To the best of the author's knowledge, the bypass injection for swirl generation in ejectors has not yet been explored in the literature.For this design, a compromise has to be made between moving the swirl-bypass inlet downstream for better off-design performance or upstream in the mixing chamber for optimal mixing performance.This is studied further in this work.
Different bypass swirl designs are possible, such as implementing a swirl chamber where the swirl is generated before being sucked into the ejector.However, in this work, a design where the bypass channel is directly inserted tangentially into the mixing chamber is explored.
The swirl-bypass port is illustrated with the main parameters shown in Figure 2.These parameters are: the relative bypass location measured from the nozzle exit position, x b = (x bypass inlet − x motive outlet )/D mix ; the bypass inlet width W and length L; the bypass tangential angle α t and the bypass inclination angle α i .The tangential angle, α t , is the angle formed in the cross-sectional plane of the mixing chamber (y-z-cross section in Figure 2) between the center line of the bypass inlet and the radial vector from the mixingchamber center at their intersecting point.This angle can vary between 0 • and 90 • , where 0 • corresponds to a bypass with no tangential component and 90 • corresponds to a fully tangential inlet.The bypass inclination angle is the angle between the mixing chamber centerline and the bypass inlet in the axial direction.In the suction nozzle, the suction flow is slightly expanded through a converging channel.This converts a part of the suction pressure into velocity.Similarly, a compression ratio, CR = A in /A out , is used to correct the bypass-inlet area.This area must also be corrected for the larger area on the cylindrical mixing chamber surface.The area correction is then added to the inlet of the suction bypass nozzle.Due to the tangential angle of bypass entrance, the problem is inherently three dimensional and will, in this study, be investigated with 3D CFD modeling.Here, the integration of suction and bypass ports is not considered and is left for future studies.The bypass inlet dimension was defined with a length L = D mix/2 and a width W = D mix (1 − cos(30 • ))/2, such that 30 • of the mixing-chamber circumference would be open.These numbers were chosen based on an initial calculation, where it was found that a larger width would limit the swirl production as a significant amount of flow would be angled toward the ejector center line with a minimal tangential component.The ejector geometry that this inlet is tested with is the geometry presented in Table 1 and was experimentally tested and numerically validated in this article.This was chosen to have a fair model validation for the baseline geometry and to compare any performance improvement or reduction with an already optimized design.However, it is noted that performance improvements are believed to be better for an ejector designed for higher pressure lift and lower entrainment ratio, discussed further in Section 6.

Model Validation
The HEM has been previously extensively studied for two-phase R744 ejectors [36,46,77].In general, the accuracy is considered to be within 5% for motive mass flow rate and 10-20% for suction mass flow rate [51] when near super-critical motive conditions are considered.Further validation and mesh independence verification are conducted with new experimental data in the following sections.

Experimental Test Campaign
To validate the 3D HEM model results, an experimental test campaign was set up at the NTNU/SINTEF Energy Research laboratory in Trondheim-Norway.The experimental rig is an R744 system with one base-load compressor and two parallel compressors, two R744 gas-coolers for heat rejection to a glycol-and a water loop, two evaporators, and an expansion device test section where a novel ejector design is installed.The detailed system description can be found in the work by Banasiak et al. [78].The system is equipped with pressure-, temperature-, and mass-flow sensors, shown in Table 2.The system has been installed with a novel R744 two-phase ejector designed for lowpressure lift and high entrainment ratio intended for heat pump conditions.The main ejector-shape design parameters are listed in Table 1.These ejector dimensions were based on a design optimization using an advanced 1D ejector model [79].The ejector was produced in stainless steel with a provided wall-roughness estimate of Ra0.8.To evaluate ejector performance additional pressure, temperature and mass-flow rate sensors were installed as displayed in Figure 3.The full range of experimental test points is shown in Table A1 in the Appendix A.

Comparison with Experimental Results
Four data points will be used for validation of the CFD model corresponding to different motive conditions at 90 bar, 80 bar, 75 bar, and 83 bar and pressure lifts of 2.1 bar, 3.6 bar, 1.5 bar, and 2.5 bar, respectively.These cases will be referred to as cases I-IV and correspond to the experimental data points 49, 59, 62, and 79 in Table A1, respectively.Case I is a measurement at a higher operating pressure and intermediate pressure lift, case II is a measurement at an intermediate operating pressure and high-pressure lift, case III is a measurement at low operating pressure and low-pressure lift, and case IV is a measurement at an intermediate operating pressure and intermediate pressure lift.
The experimental results and the calculated CFD results at those operating conditions are shown in Table 3.Here, the subscripts m, s and o refer to motive, suction and outlet conditions, respectively.The experimental measurements were calculated from a time series where close to a steady state was achieved.The standard deviation over the measurement period is, therefore, calculated for the mass flow rates and are shown along with the experimental validation data.The results from the simulations compared to the experimental data are shown in Table 3.As found previously in the literature [51], a good agreement between experimental measurements and computed solution for the motive MFR is found for the high motive pressure operating conditions, with less than 3%, 6% and 6% difference in predicted MFR for cases I, II, and IV, respectively.As the pressure and temperature are reduced, the motive MFR is underpredicted by −12% by the model in comparison to experimental results, shown in Figure 4.This is in agreement with previous findings with the HEM and is in the literature attributed to non-equilibrium phase change [46].The suction MFR prediction is, in comparison, more challenging to accurately reproduce and the accuracy is, in general, lower than for the motive MFR.Here, the suction MFR is over-predicted by 3%, 31%, 14%, and 51% for cases I, II, III, and IV, respectively.The prediction is thought to be better for the high-pressure cases (I, II) due to the better prediction of the motive flow structure.For the lower pressure case (III), the motive flow structure is less accurately reproduced locally in the ejector, which reduces the prediction accuracy of the mixing process and, therefore, the suction MFR.The suction MFR prediction error is the largest for case IV.Accurate prediction of the suction flow is a common challenge in modeling R744 ejectors, and errors up to 100% are not uncommon.These errors are typically attributed to turbulence modeling; however, currently, there exists no consensus on this [51].

Mesh Independence Study
A mesh study was conducted to assure that the 3D physics effects are well-presented by the mesh.The four validation points, I, II, III, and IV are calculated using meshes with varied mesh resolutions.Using the ANSYS ICEM software, a 3D-structured mesh was generated with high orthogonal quality and low skewness in the main dimensions.The meshes were generated based on a high-quality in-house meshing algorithm for ejectors, implemented into the ejector CFD automation toolset presented by the author in previous openly available work [57].This will allow for machine-learning-based optimization of the ejector bypass inlet in future works.An example of the meshes generated is shown in Figure 5. Here, the impact of mesh resolution on the motive flow shock train and the performance parameters is analyzed.Different meshes are compared and referred to as meshes A-E.The meshes A, B, and C were generated with 700,000, 1,000,000, and 2,000,000 cells, respectively.Two additional meshes, D and E, were generated with higher resolution in the premixing chamber and along the ejector walls, respectively.The cell counts for these meshes were 1,400,000 cells for mesh D and 2,500,000 cells for mesh E.

Flow-Rate Prediction
The mesh independence of the mass flow rate prediction must be verified.The motive and suction MFR predicted at different resolutions of mesh quality for the different cases are shown in Table 4.For all cases, the motive mass flow rate prediction is nearly independent of the mesh resolution, and only a 1% overprediction error is found between the coarsest and finest meshes.This is in agreement with previous findings with the HEM [51] and is ascribed to the choked motive flow in the ejector.For all cases, the suction MFR is reduced by increased mesh resolution until reaching around 2000k cells.For case I, the difference between the coarsest and finest mesh in suction MFR prediction is 6%.For cases II and IV the refinement has a large impact on the suction MFR as the mesh is refined to 2,500,000; cells.For mesh E (2,500,000) the accuracy of the suction MFR is, for cases I, II, and III, close to the experimentally obtained data points, and this mesh is concluded to be sufficiently refined for further computation.The results with meshes A, B, C, and E are compared for the resolution of the motive flow shock train.The pressure distribution along the center ejector axis for case I is shown in Figure 6.Meshes A, B, and C underestimate the shock strength in comparison to mesh E. A more resolved mesh in this region lowers the minimum pressure reached in the shock.
Figure 7 shows the Mach number distribution in the ejector for case IV for two mesh resolutions (1,400,000 and 2,500,000) corresponding to mesh D and E. The Mach number lines for Mach 1 are illustrated with a red contour line.In the figure, it is clear that mesh refinement to 2,500,000 cells improves the detailed resolution of the shock train.The results with mesh E show that the high-velocity core is extended by refining the mesh.This is attributed to the reduction in numerical viscosity from mesh refinement.The difference in predicted mass flow rates is very small, so the effect of this additional resolution does not significantly affect the entrainment of the secondary flow.In total, it was decided to use Mesh E with 2,500,000 cells for the following calculations.

Comparison of Turbulence Models
The choice of turbulence model was based on an initial comparison of the k − realizable model and the k − ω SST model models.A calculation using mesh D for experimental test case I with the two turbulence models was conducted and the results are shown in Table 5.Both models predict the same motive MFR as this is primarily governed by supersonic choking.For prediction of the suction MFR, it is found that the k − ω model predicts the experimentally measured MFR accurately with an over-prediction error of about 2%.The k − model, on the other hand, severely under-predicts the suction MFR with an error of approximately −38%.This difference in predictive ability could be influenced by many parameters, such as mesh resolution, geometry, and model parameters.In addition, these turbulence models are not tuned for multiphase flows, so these findings are not generalizable outside two-phase R744 ejectors.Still, the findings agree with the detailed turbulence model study of Bodys et al. [50], where the k − ω model predictions were, overall, more accurate and generally over-predicted the suction MFR, and the k − realizable model predictions were overall less accurate and generally under-predicted the suction MFR.Based on these calculations, the authors decided to continue using the k − ω SST model.

Bypass Inlet
To verify the importance of meshing and simulating the bypass inlet section, simulations with the bypass inlet section geometrically resolved were compared to simulations with the bypass specified as a boundary on the ejector wall.Introducing the bypass entrance poses a meshing challenge due to the high gradients and large cell skewness near the interception of the mixing chamber and the bypass inlet.Still, simulating the entrance was numerically preferable to specifying the boundary conditions on the ejector wall, as it yielded numerically stable results.Meshing and solution of the bypass entrance region add computational cost to the simulations; however, it also improves the accuracy of the model as a more physically reasonable boundary condition could be specified at the start of the bypass inlet.Furthermore, this removes the need for specifying boundary conditions with specified swirl and axial components, which introduces additional model uncertainties as these components are determined by the bypass inlet geometry.

Results
In this work, the swirl-bypass concept is explored by investigating the influence of its main geometry parameters, and their influence on the ejector flow field.First, different swirl-bypass ejector geometries are investigated and discussed regarding optimal design.Secondly, the flow structure and swirl characteristics are explored and analyzed in detail.The simulations are conducted with operating conditions corresponding to case IV and case I.The simulations are set up with boundary conditions for the suction and bypass inlet at the same pressure.The overall performance of the ejector is, therefore, calculated with the total entrainment of the ejector, Equation (2).
In general, the simulations with the bypass inlet showed that gas suction could be achieved; however, the swirl flow would significantly reduce the suction mass flow rate, with a net negative impact on entrainment.This is found in all configurations of geometric design and operating conditions.
The CFD simulations were computed on the NTNU IDUN computing cluster [80].The cluster has more than 70 nodes and 90 GPUs.Each node contains two Intel Xeon E5-2630 v2 CPUs, at least 128 GB of main memory, and is connected to an Infiniband network.Storage is provided by two storage arrays and a Lustre parallel-distributed file system.Computation time using one 20-core node for a single data point was 1-2 weeks depending on the specific geometry.

Bypass Inlet Design
A study of the different geometric parameters is conducted.The results are compared to a baseline simulation without a bypass inlet.Due to the large computational cost of each simulation, the study is limited to investigating the swirl-bypass position, the swirl inclination-and tangential angle, and the pressure lift at these two operating conditions.The results are shown in Table 6.
Table 6.Parameter study of the different geometric features of the bypass inlet on bypass performance.Simulations are identified by the simulation ID Sim, where the prefix N indicates no-bypass.The column named Diff is defined as the difference in percentage between the total entrainment (bypass and suction) compared to the same operating condition without bypass.

Case
Sim The first main parameter for a bypass inlet design is the downstream location of the bypass inlet.The bypass locations are defined along the mixing section in increments of mixing chamber diameters.The bypass inlet position was varied between 2 × D mix and 10 × D mix .
The influence of bypass inlet location is visualized in Figure 8.The figure shows that the choice of inlet bypass position has a major impact on the total pumped flow rate and total entrainment ratio, Equation (2).For a bypass location close to the ejector motive nozzle, the suction effect is more significant, producing a relatively larger bypass flow rate, approximately 1.6 g/s, and 2.1 g/s, for case IV and case IV, respectively.In comparison, the suction nozzle produces 29.5 g/s and 32.8 g/s, for case IV and case IV, respectively.As the bypass port is moved downstream in the mixing chamber, the bypass flow quickly drops to nearly zero.On the other hand, the entrained suction flow increases as the bypass inlet is moved toward the end of the mixing chamber.This is attributed to the bypass flow blocking and disturbing the flow path of the suction flow entering the mixing chamber.Compared to the ejector without a bypass inlet, the effect of the bypass is to reduce the total suction flow as the total flow rate accounted for by the suction nozzle is approximately 10 times larger than the bypass nozzle.This is seen for both operating conditions and all configurations of the bypass port.This blockage effect is illustrated in Figure 9. Just after the bypass inlet, colored in blue, the axial velocity through the ejector is reduced to small or even locally negative values.This forces the suction flow that would pass through this area to deviate its flow path, which drastically reduces suction flow.Additionally, the path taken by the suction flow around the bypass flow creates a vortex, illustrated in Figure 10, that disturbs the main motive flow and is moved off-center.The effect of the flow inclination and tangential angles are evaluated by varying these parameters for the bypass geometry at bypass position x b = 4.The inclination angle is varied between 0 • to 45 • .The variation in this parameter has only a minor effect on the suction and bypass flow rates.Increasing the flow inclination angle to 20 • has a positive impact on the suction flow, increasing it by approximately 1-2% for cases IV and I, and had only a minor negative impact of 4% on the bypass flow for case IV and a negligible impact for case I. Further increasing the inclination angle to 45 • further increases the suction mass flow rate by approximately 1%.However, this reduces the bypass flow rate by approximately the same amount.This is in agreement with the findings of Bodys et al. [70], where the impact of the bypass inlet angle was found to be negligible.
Varying the tangential angle of the bypass inlet has a significant impact on both suction and bypass flow.Reducing this angle directs the bypass flow more directly toward the central flow inside the mixing chamber.The central low pressure is, therefore, able to pump more bypass flow into the ejector.However, this comes at the cost of further blocking and disturbing the main suction flow.By decreasing the tangential angle to 60 • , the bypass flow rate would increase by 85% for case IV and 66% for case I. Further decreasing this angle to 40 • from 60 • increased the flow rate by a further 8%.Similarly to the bypass position, the increase in bypass flow rate reduces the suction flow rate.By reducing the tangential bypass angle, the swirl component of the bypass velocity is reduced as the inlet direction is changed from tangential to radial.

Pressure Lift
The effect of varying pressure lift depends on the bypass geometry.For the bypass entrance near the end of the mixing chamber (x b = 10, simulation 4 and 10) the bypass is significantly increased by 730% when the pressure lift is reduced from P lift = 2.5 bar to P lift = 1 bar.While for a bypass entrance close to the start of the mixing chamber (x b = 4, simulation 14 and 16), the influence of the pressure lift on bypass flow rate is smaller, producing only an increase by 20%.Consequently, when the bypass flow is increased at a lower pressure lift, the swirl blockage effect is increased and overall performance is reduced in comparison to the ejector design without a bypass.This is seen for the simulations with a bypass inlet near the end of the ejector mixing chamber, where the difference in total entrainment between an ejector with (simulation 4 and 10) and without (simulation N3) bypass is increased from −1.6% to −14%.However, for bypass inlets near the start of the ejector mixing chamber, this difference is reduced from −20.5% to −18.8%.This increase in performance could be attributed to the higher bypass performance of the bypass inlet for lower pressure lifts.It is noted that the ejector configuration with an inlet further toward the end of the mixing chamber at a low pressure lift corresponds most closely with the conventional bypass ejector design, for instance, discussed by Bodys et al. [70], where the bypass is introduced in the ejector diffuser.However, this configuration is believed to gain less value from increased mixing by introducing a swirl component as the remaining distance to exploit swirl mixing is shorter.Figure 11 shows the pressure distribution along the mixing-chamber center line.It can be seen that, for higher pressure lift, a significant increase in pressure of approximately 2 bar is seen along the mixing chamber, while only increasing by 0.1 bar in the low-pressure lift case.This implies that the position of the bypass flow rate is more sensitive to the bypass port position for high-pressure-lift operating conditions than for low-pressure lifts.This is important to design for when considering an ejector that will operate at varied pressure lifts.

Flow Structure of Swirl Mixing
Plotting the axial velocity distribution (x-direction in figures) along the radial direction at different axial positions in the mixing chamber for simulation N1, the barrel shock flow structure is seen, see Figure 12.The flow structure is symmetric around the ejector center axis.This structure is dampened and smoothed out as the flow moves toward the end of the mixing chamber.At approximately x b = 4, the velocity profile reaches the classical parabolic shape with one velocity maxima along the flow center.These profiles are then used to analyze the boundary layer thickness along the mixing chamber in an ejector without a bypass inlet according to Equation (8), shown in Figure 13.The mixing layer thickness starts very thin with a thickness below 0.1 mm due to the two high-velocity peaks of the barrel shock.As the velocities are decreased, the mixing layer increases to a maximum of 0.9 mm at x b 8 [-], near halfway through the mixing chamber.The mixing layer then decreases further down the mixing chamber as the velocities are largely evened out.
Referring to the direction towards the bypass entrance (negative z-direction) as angle θ bp = 0, the velocity profiles for the different directions in a swirl-bypass flow is shown in Figure 14 and 15 for simulation 5. Comparing this to the flow structure of the mixing layer without a swirl bypass, it is clear that symmetry is lost.There is a clear low-velocity zone right after the bypass entrance with close to negative axial velocity values, which, as the flow follows the swirl direction path, creates a high-velocity zone on the opposite side of the mixing chamber (θ bp = 180 • ).
Figure 16 shows the mixing-layer thickness for different directions of velocity profiles calculated based on Equation (8).The mixing-layer thickness has a similar profile to the result without bypass.The figure shows that the mixing layer is for all angles increasing up until a maximum near x b = 6, with a local minimum near x b = 4, the location of the bypass inlet.Near the bypass inlet (x b = 4), as the flow enters and reaches θ bp = 90-180 • , a lowvelocity zone is produced.This low-velocity zone increases the velocity difference between the flow core and the wall, which makes the mixing-layer thickness smaller.The flow that would go through this zone is instead diverted to the region on the opposite side of the bypass inlet, θ bp = 270 • , which produces a corresponding high-velocity zone.This highvelocity zone in the suction flow region evens out the flow profile causing the wider mixing layer.As the flow develops down the mixing chamber, these high and low-velocity regions are shifted with the swirling of the flow.Eventually, the high and low flow regions are inverted and the mixing thickness is largest for angle θ bp = 90 • and smallest for the angle θ bp = 270 • .The non-symmetric swirl structure then continues further into the diffuser.It is observed that, both for the swirl bypass and the standard ejector geometry, Figures 13 and 16, the mixing-layer thickness is evolving along the mixing chamber.Compared to theoretical estimates of this spreading rate of the mixing layer of a jet [81,82], the spreading rate of the mixing layer in the ejector varies significantly over the ejector mixing chamber.In this figure, no region of steady mixing-thickness growth is available for comparison.This is due to the barrel-shock structure in this ejector, which makes comparison challenging.Figures 17 and 18 show the development of the tangential velocity along the different orthogonal radial directions.Figure 17 clearly shows a spike of positive tangential flow near the wall at the starting location of the bypass inlet (positive radial direction).Towards the center, the flow travels in the opposite direction, setting up a vortex that dies out further downstream in the mixing chamber.On the opposing side, the flow tangential-flow direction has evened out over the cross section.In the y direction, Figure 18 shows that, right before the bypass entrance, a significant increase in tangential velocity occurs in the opposing direction to the bypass induced flow.These tangential flows in both positive and negative y directions are both moving away from the bypass inlet and towards the other side of the ejector mixing chamber.This is believed to be caused by the space taken up by the bypass flow displacing the suction flow.Downstream from the bypass inlet, the tangential velocity profile from the bypass flow is established.

Discussion
In this work, the bypass-swirl-generation inlet has been numerically investigated with CFD.The applied HEM CFD model was validated with new experimental results with a new ejector geometry.The results show the highest accuracy for high pressure and temperature motive operating conditions where predictions for motive mass flow rate are within an error 6%, which agree with previous findings in the literature with similar HEM approaches [51].The validation work presented in this article has highlighted that significant differences are still observed between CFD modeling and measurements.This is especially true for the suction flow, where between 3 and 50% prediction error is observed.This is believed to be due to the local prediction accuracy of motive flow, the lack of accurate turbulence models that account for two-phase effects, and the lack of appropriate multiphase models for this application.In general, the numerical models are limited for high accuracy prediction of the swirl-bypass and suction flow rate.These will depend on the accurate prediction of local flow variables such as pressure and velocities, especially in the vicinity of the bypass inlet.As of yet, very limited local experimental validation data is available for R744 ejectors, primarily due to the small size, high pressures, high velocities, and multiphase flow inside the device.Still, comparison between different models can give valuable insight into the impact of different ejector design variables.
An ejector-bypass concept intended to generate swirl in the ejector mixing chamber is tested for various geometric configurations.The results show that the key impact of the swirl-bypass inlet for this specific ejector geometry is to disturb the suction flow, creating a non-symmetric mixing that reduces ejector performance.This effect is dominant regardless of the bypass angles and position.The main cause of the decrease in performance is attributed to the blocking effect of the bypass flow, limiting the available flow area of the suction flow, the turbulent energy loss induced by the non-optimal collision between these two flows, the lack of swirl component velocity recovery at the outlet, and non-symmetric effects induced by the swirl inlet.Another significant feature of the bypass swirl flow is the non-symmetric flow structure caused by having only one bypass inlet.The non-symmetry caused significant disruption to the flow path of the suction and motive flow that carried on into the diffuser.Based on these simulations, it is found that 3D effects, such as induced non-symmetric flow, are highly important to accurately predict ejector performance with swirl, and that, in general, 2D flow solutions are inadequate for designing such devices.

Open Access
To promote open access to research into environmentally sustainable solutions with ejectors, the research produced will be made available openly.The ejector geometry, code for the generation of high-quality ejector swirl bypass mesh in ANSYS ICEM, experimental data, and velocity and mixing layer calculations are available at: https://github.com/knutringstad/Swirl_bypass_ejector (accessed on 10 September 2022).

Further Work and Future Design Considerations
Based on the results of this study, several improvements and suggestions for future study of swirl bypass ejectors have been found.The primary reason for the poor performance of the presented ejector is the blockage effect that is caused by the bypass flow obstructing the path of the suction flow.Increasing the diameter of the ejector mixing chamber is believed to allow for the bypass flow to flow primarily in the outer section of the mixing chamber, near the wall.This would reduce the swirl component's radial extension into the suction flow, which, in turn, will reduce the efficiency of the swirl mixing.However, the bigger mixing chamber would allow the suction flow to be less disturbed by the bypass inlet, improving the suction flow rate.Secondly, varying the size and shape of the bypass opening could reduce the impact of the suction flow blocking.Reducing the size of the opening will increase the velocity of the flow from the bypass inlet.This could be applied to intensify mixing while potentially reducing the blocked area of the suction flow.Optimization of the shape of the opening is also of interest; however, for simplicity of component manufacturing rectangular or circular openings are preferred.
Another significant feature of the bypass swirl flow is the non-symmetric flow structure caused by having only one bypass inlet.The non-symmetry caused significant disruption to the flow path of the suction and motive flow that carried on into the diffuser.By increasing the bypass inlet velocity (reducing bypass port flow area) the flow would reach further around the perimeter of the mixing chamber.This could reduce non-symmetry as the bypass flow is more evenly distributed around the mixing chamber axial direction.Another potential design improvement would be to introduce multiple inlets that were offset from each other.Such a solution would reduce non-symmetric effects that push the flow off center.This solution will, however, require additional study to find a good design in terms of manufacturability and port connections.
In this study, the swirl-velocity component is quickly smoothed out by viscosity as the flow progresses through the ejector.However, this effect may be significant for an ejector that produces a larger swirl component.In that case, the diffuser needs to be designed for swirl-component recovery.
Future studies should investigate using a holistic approach to device design.This work has only considered the swirl-bypass inlet as an addition to an already high-performance ejector design.It is, therefore, difficult to gain the largest benefit from the swirl mixing without integrating the bypass port with the other ejector design parameters.This integration could be considered by the use of optimization techniques such as those used in the authors' previous work [57].The intended application of the swirl bypass inlet is to generate a swirl flow that increases the mixing of motive and suction flow.Each ejector is designed to have this occur inside the mixing chamber.When the pressure lift is lower than design conditions, a higher potential for work recovery is present; however, the mixing chamber will be under-dimensioned for the entrainment ratio.Here, the potential for ejector performance improvement by swirl mixing is the largest.In this work, an ejector designed for already low-pressure lifts was used, so limited improvements could be gained.Future studies should investigate the optimal ranges of swirl mixing operation and these should be integrated into the ejector design.This work has not considered the integration of the bypass and suction port, which should, in general, be designed practically and compactly with limited flow losses.The swirl-bypass concept can also be combined with other ejector concepts.Combining a swirl-bypass and a diffuser-bypass concept could potentially improve ejector entrainment by improving the performance of the ejector under different off-design conditions.Combining motive-inlet swirl and bypass-swirl generation can potentially generate synergistic effects if optimized appropriately to generate additional swirl flow.Advanced ejector design concepts, especially in combination, require appropriate control strategies for optimal performance, which need to be explored in future work.
The authors also highlight the potential of the second application of a swirl-bypass inlet, referred to in this paper as concept (B), for motive flow control.Such a device has not yet been explored and is left for further work.Most of the design suggestions mentioned are also applicable to such a design.The bypass-design considerations regarding positioning of the bypass inlet for the optimal flow rate are then not applicable, as these are primarily decided by the motive pressure.

Conclusions
In this work, a novel swirl-bypass ejector design concept is explored using a highfidelity 3D CFD model.To verify model accuracy, the multiphase CFD model is validated against new experimental data collected at a range of operating conditions with a novel ejector design.A comparison of experimental and simulation data showed good agreement (within a 6% prediction error) for the motive flow rate prediction for high operating pressures, and poorer agreement (within a 12% prediction error) for lower operating conditions, as is often seen in the literature.The ejector swirl bypass concept performed poorly with a single inlet port.The simulations showed that the generated mixing effect of the swirl flow was minimal in comparison with the reduced efficiency due to the effect of the bypass flow blocking the suction flow area.The reduction in performance of the total entrained flow ranged from −2% to −20%, depending primarily on bypass position.This effect persisted for all different configurations of the bypass port; however, it was lower for bypass positions toward the end of the mixing chamber where less bypass flow rate is produced.The bypass flow rate was found to be highly sensitive to pressure lift for inlets near the end of the mixing chamber, and less sensitive when positioned closer to the motive nozzle.The bypass flow also introduces a non-symmetric effect that pushes the primary and secondary flow out of the ejector center line.This is believed to be causing large losses due to the onset of recirculation zones in the diffuser.The mixing-layer thickness was analyzed with and without the swirl bypass and the off-center motive flow was found to cause the mixing layer between motive and suction flow to be thicker and thinner on either side of the ejector mixing chamber.Lastly, different possibilities for improvements to this concept are suggested that could bring up the efficiency of this ejector design concept.These design improvements include:

•
Increase in mixing-chamber diameter.

Figure 1 .
Figure 1.Illustration of the two swirl-bypass concepts (A,B).Concept (A) is a suction swirl bypass ejector.Concept (B) is a motive swirl bypass ejector.

Figure 2 .
Figure 2. Generic ejector geometry with geometry parameters.Gray color signifies solid parts.

Figure 3 .
Figure 3. Simplified illustration of the system layout with new ejector and sensors installed.PT-pressure/temperature sensor, DP-differential pressure sensor, M-mass flow meter.

Figure 4 .
Figure 4. Comparison of experimental and numerical results in terms of motive (upper) and suction (lower) mass flow rate prediction.

Figure 5 .
Figure 5. Generic ejector geometry with geometry parameters.The red color signifies solid parts.

Figure 6 .
Figure 6.Pressure distribution for case I along the x-axis for meshes A, B, C, and E.

Figure 7 .
Figure 7. Mach number distribution for meshes E (top) and D (bottom) for case IV.Red line illustrates Mach 1.

Figure 8 .
Figure 8. Plot of bypass (blue) and suction (green) mass flow rate as a function of bypass inlet position.The total pumped flow rate is shown in red.

Figure 9 .
Figure 9. Top: Contour plot of the axial velocity in gray.The bypass inlet is colored blue.Bottom: Cross-sectional vector plots of the tangential velocity shown for locations along the axial direction after the bypass inlet.

Figure 10 .
Figure 10.Flowpath illustration of the flow coming from the bypass inlet.Different cross-sectional contour plots of the tangential velocity are shown for locations along the axial direction after the bypass inlet.5.1.2.Flow Inclination Angle, α i , and Flow Tangential Angle, α t

Figure 11 .
Figure 11.Pressure variation along the mixing-chamber center line for simulations N4 and N5.

Figure 12 .
Figure 12.Velocity plot in the y-direction for the ejector without a bypass inlet.The flow is symmetric around the center line in this flow.

Figure 13 .
Figure 13.Mixing layer thickness along the mixing-chamber center line without a bypass inlet.

Figure 14 .
Figure 14.Velocity plot in the z-direction.Positive r direction is in the direction towards the bypass inlet, θ bp = 0 • .

Figure 15 .
Figure 15.Velocity plot in the y direction.Positive r direction is in the direction furthest offset from the flow path of the bypass inlet, θ bp = 270 • .

Figure 16 .
Figure 16.Plot of the mixing-layer thickness for an ejector with a bypass inlet calculated towards the wall in different radial directions plotted along the mixing-chamber center line.

Figure 17 .
Figure 17.Tangential velocity plot in the z direction.Positive r direction is in the direction towards the bypass inlet, θ bp = 0 • .

Figure 18 .
Figure 18.Tangential velocity plot in the y direction.Positive r direction is in the direction furthest offset from the flow path of the bypass inlet, θ bp = 270 • .

Table 1 .
[57] dimensions of the ejector geometry with dimensions as defined in Ringstad et al.[57].

Table 2 .
Characteristics of the measuring devices.

Table 3 .
Comparison of experimental and CFD results at different operating conditions.δ = MFR CFD − MFR Exp .The MFR measurements are shown with the standard deviation over the experimental measurement period.Subscript exp refers to an experimentally obtained value.

Table 4 .
Mesh convergence study with different mesh sizes.

Table 5 .
Comparison of turbulence the k − ω SST and k − realizable models for Case I with mesh D.
• Reduction in bypass inlet size.•Optimization of bypass inlet opening shape.•Addingadditionalbypass inlets with a constant offset.Additionally, several directions for future research have been presented that can potentially improve future generations of the swirl-bypass ejector design: