Abstract Near-wall modelling is one of the most challenging aspects of CFD computations. In fact, integration-to-the-wall with low-Reynolds approach strongly affects accuracy of results, but strongly increases the computational resources required by the simulation. A compromise between accuracy and speed to solution is usually obtained through the use of wall functions, especially in RANS computations, which normally require that the first cell of the grid to fall inside the log-layer (50 < y + < 200) [1]. This approach can be generally considered as robust, however the derivation of wall functions from attached flow boundary layers can mislead to non-physical results in presence of specific flow topologies, e.g. recirculation, or whenever a detailed boundary layer representation is required (e.g. aeroacoustics studies) [2]. In this work, a preliminary attempt to create an alternative data-driven wall function is performed, exploiting artificial neural networks (ANNs). Whenever enough training examples are provided, ANNs have proven to be extremely powerful in solving complex non-linear problems [3]. The learner that is derived from the multi-layer perceptron ANN, is here used to obtain two-dimensional, turbulent production and dissipation values near the walls. Training examples of the dataset have been initially collected either from LES simulations of significant 2D test cases or have been found in open databases. Assessments on the morphology and the ANN training can be found in the paper. The ANN has been implemented in a Python environment, using scikit-learn and tensorflow libraries [4][5]. The derived wall function is implemented in OpenFOAM v-17.12 [6], embedding the forwarding algorithm in run-time computations exploiting Python3.6m C_Api library. The data-driven wall function is here applied to k-epsilon simulations of a 2D periodic hill with different computational grids and to a modified compressor cascade NACA aerofoil with sinusoidal leading edge. A comparison between ANN enhanced simulations, available data and standard modelization is here performed and reported.

Abstract The multidisciplinary performance of a bypass jet propulsion engine is investigated and simulated in this paper by using the pseudo bond graph approach. This theory is applied in the in-house software tool ASTOR (AircraftEngine Simulation of Transient Operation Research) to simulate the different operating points and the system dynamics in transient load cases with high accuracy. In this paper, the modulation of the V2500-A1 turbofan jet engine in the pseudo bond graph notation is described in detail. The thermodynamic cycle of a steady-state operating point is calculated and a transient load case is investigated to simulate the system performance. The results correspond with the data from a commercial software and with literature. However, some significant discrepancies are detected, which can be explained by the current model simplification. In a nutshell, ASTOR is able to simulate the dynamic performance of the overall system and can be used to improve the process of the construction or regeneration of a complex capital good such as a jet engine.

Abstract This study presents a CFD comparison of a piccolo tube inlet anti-icing system with benchmark flight test data and a shielded swirl design. Aircraft engine inlets must prevent ice accretion that could cause damage or inhibit safe operation. A common way to meet this requirement is to heat the inlet with hot air distributed with a piccolo tube or a swirl anti-icing system. Piccolo tube systems generally use less engine bleed air but are more complex and heavier than swirl systems. The CFD model for this study includes all the inlet key characteristics such as piccolo tube holes, material properties, and thickness. The model is evaluated in a two-step process. Water collection rates for flight in icing conditions are calculated. These results are then mapped onto a conjugate heat transfer film runback model. The CFD results include surface temperatures and the amount of liquid film runback that flows off the heated surface. Benchmark test data is presented illustrating good inlet surface temperature agreement with the piccolo tube CFD model for flight in both dry air and icing conditions. This CFD modeling approach is also applied to a unique shielded swirl anti-icing system. The results of this study confirm that a shielded swirl anti-icing system can be designed to provide the same level of anti-icing protection that is achieved with a piccolo tube configuration without overheating the inlet or requiring more engine bleed flow.

Abstract The paper provides design and performance data for two envisaged year-2050 state-of-the-art engines: a geared high bypass turbofan for intercontinental missions and a contra-rotating pusher open rotor targeting short to medium range aircraft. It defines component performance and cycle parameters, general powerplant arrangements, sizes and weights. Reduced thrust requirements for future aircraft reflect expected improvements in engine and airframe technologies. Advanced simulation platforms have been developed, using the software PROOSIS, to model the engines and details of individual components, including custom elements for the open rotor engine. The engines are optimised and compared with ‘baseline’ year-2000 turbofans and an anticipated year-2025 entry-into-service open rotor to quantify the relative fuel-burn benefits. A preliminary scaling with non-optimised year-2050 ‘reference’ engines, based on Top-of-Climb (TOC) thrust and bypass ratio, highlights the trade-offs between reduced specific fuel consumption (SFC) and increased weight and engine diameter. These parameters are then converted into mission fuel burn using linear and non-linear trade factors from aircraft models. The final turbofan has an optimised design-point bypass ratio (BPR) of 16.8, and a maximum overall pressure ratio (OPR) of 75.4 for a 31.5% TOC thrust reduction and a 46% mission fuel burn reduction per passenger kilometre compared to the respective year-2000 baseline engine and aircraft combination. The final open rotor SFC is 9.5% less than the year-2025 open rotor and 39% less than the year-2000 turbofan, while the TOC thrust increases by 8% versus the 2025 open rotor, due to assumed increase in aircraft passenger capacity. Combined with airframe improvements, the final open rotor-powered aircraft has a 59% fuel-burn reduction per passenger kilometre relative to its year-2000 baseline.

Abstract The development of Ultra-High Bypass Ratio (up to 20) engines with the aim of improving the propulsive efficiency, introduces new challenges for the transmission system in terms of heat management and power losses, since the amount of power transferred through the gearbox is greatly increased. In this respect the accurate estimate of losses at the various flow regimes realized during a typical aeroengine mission within a Power Gear Box (PGB) is essential for the correct design and operation of the engine itself. This paper proposes a computational methodology to estimate all fluid-dynamic (load-independent) losses, which become a major source of dissipation at the high rotational speeds typical of aeroengines, developing within an epicyclic PGB. The overall procedure is based on the superposition principle and approaches the three fluid-dynamic losses, namely injection, windage and meshing losses, with different numerical techniques. The simulations of windage effects, which consider the actual PGB geometry including the carrier disk, the lubricant spray-bar and the external casing, are based on steady-state computations. In order to introduce such simplification, a scaling procedure that avoids interference of the stationary and rotating interfaces was implemented following the outcomes of a previous analysis. The computation of meshing losses employs a fully unsteady dynamic mesh approach and considers a portion of the meshing gears only. Both the sun-planet and planet-ring meshing were considered showing that the latter introduces a much lower level of losses. Finally the injection losses are calculated considering the oil jet momentum variation with simplified methods based on 0D modelling. The proposed procedure, based on the superposition principle and applied to a planetary power gear train, is tested against experimental results described in a previous paper focused on a meshing gear pair.

Abstract Body force modelling in numerical simulations of axial compressors and fans has been used extensively in the literature to assess aerodynamic performance in uniform and non-uniform flow at relatively low computational cost. Existing approaches require calibration or, in the case of purely analytical models, are unable to accurately predict losses. It can also be challenging to capture the chordwise loading distribution with existing analytical models. This paper introduces a new body force modelling approach in which blade loading is computed using an isolated-airfoil, analytical model supplemented by a trained artificial-neural-network-based correction factor for finite pitch effects. The loading model is derived from potential flow theory and accounts for both camber and thickness effects; it captures both the overall and local loading. The approach is currently implemented for low-Mach number (incompressible) flows. The body forces causing flow turning derive directly from the corrected loading model. Forces arising from viscous losses are modeled by solving the integral boundary layer equations along streamlines within blade rows based on the fictitious edge velocities computed by the loading model. The viscous loss force is a function of the local dissipation coefficients. The approach is implemented within a traditional finite-volume computational fluid dynamics solver. In this paper, the application is limited to 2D cascades. To assess the approach, results from the body force model are compared to blade-to-blade solutions from MISES. The key findings are (1) that a relatively modest set of training data for the neural network produces a robust finite pitch correction, and (2) that the modelling approach is able to successfully capture the flow turning and losses associated with a variety of low-speed compressor cascades without any calibration specific to the blade row(s) being modeled.

Abstract Reynolds-averaged Navier-Stokes (RANS) equations are employed for aerodynamic and aeroelastic modelling in axial compressors. Their solutions are highly dependent on the turbulence models for closure. The main objective of this work is to assess the widely used Spalart-Allmaras model’s suitability for compressor flows. For this purpose, an extensive investigation of the sources of uncertainties in a high-speed multi-stage compressor rig was carried out. The grid resolution near the casing end wall, which affects the tip leakage flow and casing boundary layer, was found to have a major effect on the stability limit prediction. Refinements in this region led to a stall margin loss prediction. It was found that this loss was exclusively due to the destruction term in the SA model.

Abstract As a consequence of the increasing share of renewable energy sources in present-day electrical grid systems, time variations of the power demand for fossil fuel plants can become more sudden. Therefore, an ability to respond to sudden load changes becomes an important issue for power generation gas turbines. This paper describes a real-time model for predicting the transient performance of gas turbines. The method includes basic transient phenomena, such as volume packing and the heat transfer between the working fluid and the structural elements. The dynamics of components are quantified by solving ordinary differential equations with appropriate initial and boundary conditions. Compressor and turbine operating points are determined from corresponding performance maps previously calculated using sophisticated aerodynamic, through-flow codes. This includes a sufficient number of such characteristics to account for the variations in speed and machine geometry. The developed dynamic model was verified by comparison of simulation results with experimentally recorded operating parameters for a real engine. This includes the start-up sequence and the change of load. Additional simulation covers the system response to a step increase in fuel flow. The simulation is carried out faster than the real process.

Abstract Modeling of the engine behavior in windmilling conditions is an important engineering objective. The relight capability of the engine is mainly associated with the air mass flow rate that passes through the engine high pressure core in those conditions. This is one of the parameters that drive the combustion chamber volume. Predicting the engine behavior is challenging, especially early in the development process. The pressure losses along the core are distributed between the different stages of the compressors and turbines, which are operated extremely far from their design point. Engine manufacturers must anticipate with sufficient margins to ensure that the specifications are met when the engine is finally qualified in flight. This article focuses on the behavior of compressor cascades in such conditions, corresponding to high negative incidences. A recently designed high pressure compressor is studied in windmilling conditions using three complementary approaches. First, engine tests are used to obtain validated 0D data for two flight conditions. Then, state of the art Lattice-Boltzmann Method (LBM) simulations are carried out to have insight in the flow characteristics inside the compressor. They are compared to the available experimental data. Finally, a 1D model using stage by stage Euler equation for turbomachinery is used. This kind of modeling is of particular interest because it can be used early in the design process. The correlations for losses and deviation angle from the literature are modified to account for the particularity of the flow in those conditions. One shows that the three approaches give consistent results.

Abstract Despite significant advancements in computational power and various numerical modeling in past decades, flow simulation of a multi-stage axial-flow compressor is still one of the most active areas of research, for it is the critical component in engine performance and operability, and there are so many elements that need to be looked into to predicting correct matching of the stages and accurate flow distribution inside the machine. Modeling unsteadiness, both deterministic and random types, and real geometries are among the most important features to be considered in such prediction. The authors have conducted in their previous studies a series of unsteady RANS (URANS) simulations of a 6.5-stage high-speed highly-loaded axial-flow compressor, and explored many unsteady effects as well as effects of real geometries such as Variable Stator Vane (VSV) clearance and inter-stage seal leakage flow on the compressor performance. However, all the analyses failed to predict correct stage matching, total pressure and temperature radial profiles, or mass-flow with adequate accuracies. In the present study, an Improved Delayed Detached Eddy Simulation (IDDES) with SST k-omega model is applied to the simulation of the same compressor configuration at aerodynamic design point. Fifth-order WENO scheme is employed for improved spatial accuracy to suppress significant increase in mesh size. Total number of mesh points are over 400 million for 1/10th sector model. Computations are ensemble averaged for 20 sector passage. Computed overall performance and flow field are compared with the compressor rig test data. The predictions of inter-stage total temperature radial profiles are noticeably improved over the URANS with the same mesh, discretization scheme and eddy turbulence model. Good comparison with the rig data indicates the current simulation is properly capturing the span-wise mixing phenomena. Unsteady flow field are compared between IDDES and URANS to locate the cause for the enhanced mixing. It is shown that components of Reynolds stress responsible for radial diffusion and anisotropic features are intensified in the tip leakage vortex at the rotor exit for the IDDES.

Abstract Designing a turbofan to operate in distorted inlet conditions is an issue of growing interest. In such conditions however, fan design can be computationally challenging. Indeed, subject to neither axi-symmetrical nor periodic inlet conditions, computations must be carried out throughout the whole circumferential domain i.e. 360°. Besides, the classical CFD approach implies the use of URANS computations so as to capture the distortion transfer across the fan stage. Eventually, computations are too onerous to be used in design loops. In this context, this paper presents a methodology to effectively assess a fan blade design domain in distorted conditions. This methodology is based on a body-force source term approach formulated in order to accurately recreate deviations, loads and losses across the fan stage. It notably enables to gain two orders in terms of restitution time and thus the possibility to use optimization tools. The design domain of this study is based on variations of profile chord, blade leading and trailing edges angles applied at two different relative heights. A Latin Hypercube Sampling (LHS) associated with a meta-model based on Radial Basis Functions (RBF) enables to assess the impact of geometric variations on performance and operability. Although this study emphasizes that some modeling improvements are still necessary, it also demonstrates the potential of the body-force methodology to conduct fan design when subject to inlet distortion.

Abstract Turbochargers are a vital component for aiding engine manufacturers in meeting the latest emissions standards. However, their range of operation is limited for low mass flows by compressor surge. Operation in surge results in pressure and mass flow oscillations that are often damaging to the compressor and its installation. Since surge is a highly complex flow regime, full unsteady 3D models are generally too computationally expensive to run. The majority of current low-dimensional surge models use a cubic compressor characteristic that needs to be fitted to experimental data. Therefore, each time a compressor is studied using these models, costly experimental testing is required. In this paper, a new technique for obtaining an axisymmetric centrifugal compressor characteristic is presented. This characteristic is built using the equations of mass, momentum and energy from first principles in order to provide a more complete model than those currently obtained via experimental data. This approach enables us to explain the resulting cubic-like shape of the characteristic and hence to identify impeller inlet stall as a route into surge. The characteristic is used within a quasi-steady, map-based surge model in order to demonstrate its ability to predict the onset of surge while only providing geometric data as input. Validation is provided for this model by discussion of the qualitative flow dynamics and a good fit to experimental data, especially for low impeller speeds and pressure ratios.

Abstract The adoption of very high bypass ratios is an effective strategy to improve the performance of turbofans for civil aviation. Very high-bypass ratio turbofans are characterized by large-diameter fans and large intake diameters, whereas the overall length of the engine and its installation have to be contained to prevent excessive weight growth. As a result, the main components of the low-pressure compression system (LPC) are closer to each other in modern and future engines than they were in less recent designs, if distances are measured in terms of engine diameter. As the axial decay of potential distortions takes place over distances comparable to their tangential wavelength, modern and future engines also display stronger interactions between the intake, the fan, the OGVs and the pylons, causing performance and integrity issues that need to be addressed at the design stage. The direct computation of whole LPC systems through CFD is expensive on account of the large extent of the domain and the wide range of wavelengths/frequencies that need to be resolved. This makes CFD unsuitable for bypass design. However, as the interaction between intake, fan, OGV and bypass takes place mainly through flow potential over long wavelengths, drastic simplifications can be introduced. Correspondingly, considerable computational savings are possible by resorting to semi-analytical models. In this Part I of a two-part paper, a model for the intake/fan/OGV bypass interaction problem based on potential flow is presented. The model can describe the potential component of the LPC flow field, as well as vorticity distortions transmitted through the fan and those generated by the fan as a result of the non-uniform work done around the circumference. The method is validated against CFD results on a configuration typical of modern machines for civil aviation service and is demonstrated as a useful preliminary design tool.

Abstract The growth in diameter of turbofan engines exacerbates problems related to the interaction of the Outlet Guide Vanes (OGV), pylon and intake because it reduces the ratio between components gaps and disturbance wavelength. The main components of this interaction are the potential fields generated by the intake and by structural components in the bypass, the pylon and the Radial Drive Fairing (RDF). The OGV bladerow and the fan are immersed in these potential fields and suffer performance degradation as well as integrity issues as a result. Simple actuator-disc analysis shows that a uniform OGV cascade amplifies the effect of the pylon potential flow. Therefore, a number of methods have been proposed over the years to compute OGV exit flow angle patterns that result in an approximately circumferentially uniform static pressure field at fan exit. Within actuator disc approximations, the determination of the optimal exit flow angle pattern can be accomplished analytically but little information is obtained on how the geometry of the vanes ought to be modified. Consequently, it is not difficult to generate by this method OGV cascades that stall or choke locally. More recent contributions use CFD computations coupled to optimization methods to determine OGV patterns that reduce the distortion at the fan exit, while minimising some measure of OGV loss. Whilst in principle more rational, these methods encounter practical difficulties due the computational power needed to obtain reliable loss estimates while exploring large design spaces. In this paper the view is taken that the performance of the OGV bladerow can be preserved during the optimization process if the loading distribution of each vane is made to match the loading distribution of the nominal vane (i.e. the aerodynamic design intent with axisymmetric inlet and exit flow). As loading distributions are readily available from inviscid-type analysis, the generation of optimal OGV patterns can be accomplished with very reasonable computational expense using a method based on the model described in part I of this paper.

Abstract 1D modelling is crucial for turbomachinery unsteady performance prediction and system response assessment. The purpose of the paper is to describe a newly developed 1D modelling (TURBODYNA) for turbomachinery. Different from classic 1D modelling, in TURBODYNA, rotor has been meshed and its unsteadiness due to flow field time scale is considered. Instead of direct using of performances maps, source terms are added in Euler equation set to simulate the rotor. By comparing 1D modelling with 3D CFD results, It shows that rotor unsteadiness is indispensable for a better prediction. In addition, different variables response to pulse differently. In the rotor, mass flow is close to quasi-steady while entropy is significantly unsteady. TURBODYNA can capture these features correctly and provide an accurate prediction on pressure wave transportation.

Abstract Reducing the uncertainties in the prediction of turbine inlet conditions is a crucial aspect to improve aero engine designs and further increase engine efficiencies. To meet constantly stricter emission regulations, lean burn combustion could play a key role for future engine designs. However, these combustion systems are characterized by significant swirl for flame stabilization and reduced cooling air mass flows. As a result, substantial spatial and transient variations of the turbine inlet conditions are encountered. To investigate the effect of the combustor on the high pressure turbine, a rotating cooled transonic high-pressure configuration has been designed and investigated experimentally at the DLR turbine test facility ‘NG-Turb’ in Göttingen, Germany. It is a rotating full annular 1.5 stage turbine configuration which is coupled to a combustor simulator. The combustor simulator is designed to create turbine inlet conditions which are hydrodynamically representative for a lean-burn aero engine. A detailed description of the test rig and its instrumentation as well as a discussion of the measurement results is presented in part I of this paper. Part II focuses on numerical modeling of the test rig to further extend the understanding of the measurement results. Integrated simulations of the configuration including combustor simulator and nozzle guide vanes are performed for leading edge and passage clocking position and the effect on the hot streak migration is discussed. The simulation and experimental results at the combustor-turbine interface are compared showing a good overall agreement. The relevant flow features are correctly predicted in the simulations, proving the suitability of the numerical model for application to integrated combustor-turbine interaction analysis.

Abstract The performance of radial inflow turbines, and specifically of turboexpanders for oil & gas applications, has been traditionally described in terms of efficiency versus velocity speed ratio (U/C) and discharge flow coefficient (Q/N). Especially in the testing phase, this latter parameter has been often preferred to the angle setting of moveable inlet guide vanes (IGV), which are standard equipment for most turboexpanders. In practice, the expander U/C has been often considered to give the performance backbone, while the Q/N ratio has been used for secondary corrections. Moreover, although the role of pressure ratio (PR) is recognized, its impact has been experimentally unexplored in those cases where testing facilities had capacity limitations. Eventually, in case of variable nozzles, the inlet flow capacity curve has been rarely included among the output performance variables, being the attention mainly focused on efficiency. In the present paper, beside an overview and an explanation of the physical meaning of traditional performance parameters, an alternative approach based on torque mapping versus U/C is introduced and discussed in detail. As a matter of fact, numerical and experimental data show smooth and regular trends when torque coefficient is used instead of adiabatic efficiency. Moreover, performance based on torque coefficient can be more conveniently extrapolated at extreme off-design conditions such as start-up (locked rotor condition) or full speed no load. The ease of extrapolation is particularly important for machine operability, which often requires accurate modeling of transient missions at very partial loads (as for instance during start-up or shut-down). Examples will be offered to show the advantages of torque coefficient representation and how sensitive this is to IGV setting and pressure.

Abstract The meanline modeling of asymmetrical twin-scroll turbocharger turbine performance has been in continuous development since its first introduction in 2008 by Daimler AG. Asymmetric turbine features scrolls of different sizes to enable high-pressure exhaust gas recirculation (EGR), which has been proven to be effective. However, the asymmetric feature of the scrolls has been identified as the major difficulty in the meanline modeling. This is particularly true at unequal-admission operating conditions, where the two volute scrolls are subjected to different magnitude of pressure and mass flow rate. Tuning a mean-line model for one particular turbine unit is usually possible, depending on the amount of flexibility in its loss sub-models. However the transferability of loss coefficients for the use of different turbine sizes has not been studied in the case of asymmetrical twin-scroll turbine. The feasibility in doing so, and the consequences, difficulties and outcomes, will be presented in this paper. The modeling work in this study will be carried out using commercial meanline tool — RITAL™ from Concepts NREC. The loss model is first tuned for the largest asymmetrical twin-scroll turbine in Daimler AG’s product line, at different admission conditions. After the model prediction has been validated with coldflow experimental data, the same loss coefficient settings are then used for performance prediction of a smaller asymmetrical twin-scroll turbine. The finding of this paper suggests that the fraction of the volute flow exit area varies with turbine admission state. The variation trend exhibits a unique function of the turbine geometry, and found to be consistent across the turbine sizes. On the other hand, the rotor loss coefficient may remain constant regardless of the turbine operation. The goal of this study is to establish the confidence level for blind performance prediction in the future.

Abstract To create an adequate computational model of oil behaviour in an aeroengine bearing chamber previous work at the Gas Turbine and Transmissions Research Centre (G2TRC) suggests it is necessary to be able to model oil shedding from bearings, breaking up into droplets/ligaments and forming thin and thick films driven by gravity and shear. Our previously published work using Fluent successfully coupled volume of fluid with the Eulerian thin film model (ETFM) and identified the challenges coupling the ETFM with the discrete phase modelling (DPM). For this latter work comparison was made to published experimental and modelling data in which a jet is injected into a duct breaking up into droplets before forming a wall film. In this paper the use of the open-source CFD code OpenFOAM is investigated for this application recognising that such an approach eliminates some of the restrictions in a commercial product. A transient solver for spray particle cloud modelling and thin liquid film transport ( sprayParcelFilmFoam ) has been developed and incorporated within OpenFOAM. Fully coupled DPM-ETFM is presented, capable of modelling both primary atomization and secondary breakup. In addition two new film sub-models have been implemented for film stripping and edge separation. In order to achieve accurate statistical representation of droplets, modifications to the DPM particle injector code were implemented. CFD results are validated against published high speed imaging and phase Doppler experimental data and in addition there is a comparison to computational results obtained using ANSYS Fluent. The fidelity of both the solver and the novel surface film sub-models are evaluated against average film thickness measurements along the duct centreline. With the inclusion of both film stripping and edge separation, a normalized root mean squared deviation of 5.1 % was achieved when compared to film thickness measurements, improving significantly on the results obtained with Fluent. A comparison with experimental data of particle diameters and velocities downstream of the expansion edge gives good qualitative agreement. Future work is recommended to provide a better formulation for the edge-separated droplet diameters. Analysis of film momentum source terms highlights the necessity for including both the gas and hydrostatic pressure source terms within the film momentum transport equation. This CFD investigation has successfully established a fully coupled two-way DPM-ETFM approach. This work illustrates an advance in bearing chamber modelling capability and has established a necessary foundation for future aeroengine bearing chamber film modelling.

Abstract Modern aero-engines are characterized by compact components (fan, compressor, combustor, and turbine). Such proximity creates a complex interaction between the components and poses a modeling challenge due to the difficulties in identifying a clear interface between components since they are usually modeled separately. From a numerical point of view, the simulation of a complex compact aero-engine system requires interaction between these individual components, especially the combustor-turbine interaction. The combustor is characterized by a subsonic chemically reacting and swirling flow while the high-pressure turbine (HPT) stage has flow which is transonic. Furthermore, the simulation of combustor-turbine interactions is more challenging due to aggressive flow conditions such as non-uniform temperature, non-uniform total-pressure, strong swirl, and high turbulence intensity. The simulation of aero-engines, where combustor-turbine interactions are important, requires a methodology that can be used in a real engine framework while ensuring numerical requirements of accuracy and stability. Conventionally, such a simulation is carried out using one of the two approaches: a combined simulation (or joint-simulation) of the combustor and the HPT geometries, or a co-simulation between the combustor and the turbine with the exchange of boundary conditions between these two separate domains. The primary objective of this paper is to assess the effectiveness of the joint simulation versus the co-simulation and propose a more practical approach for modeling combustor and turbine interactions. First, a detailed grid independence study with hexahedral and polyhedral meshes is performed to select the required polyhedral mesh. Then, an optimal location of the interface between the combustor and the nozzle guide vane (NGV) is identified. Co-simulations are then performed by exchanging information between the combustor and the NGV at the interface, wherein the combustor is solved using LES while the NGV is solved using RANS. The joint combustor-NGV simulations are solved using LES. The effect of the combustor-NGV interaction on the flow field and hot streak migration is analyzed. The results suggest that the joint simulation is computationally efficient and more accurate since both components are modelled together.