Bogdan A. Danciu, George K. Giannakopoulos, Mathis Bode, Christos E. Frouzakis
{"title":"Multi-cycle Direct Numerical Simulations of a Laboratory Scale Engine: Evolution of Boundary Layers and Wall Heat Flux","authors":"Bogdan A. Danciu, George K. Giannakopoulos, Mathis Bode, Christos E. Frouzakis","doi":"10.1007/s10494-024-00576-w","DOIUrl":"https://doi.org/10.1007/s10494-024-00576-w","url":null,"abstract":"<p>Multi-cycle direct numerical simulations (DNS) of a laboratory-scale engine at technically relevant engine speeds (1500 and 2500 rpm) are performed to investigate the transient velocity and thermal boundary layers (BL) as well as the wall heat flux during the compression stroke under motored operation. The time-varying wall-bounded flow is characterized by a large-scale tumble vortex, which generates vortical structures as the flow rolls off the cylinder wall. The bulk flow is found to strongly affect the development of the BL profiles, especially at higher engine speeds. As a result, the large-scale flow structures lead to alternating pressure gradients near the wall, invalidating the flow equilibrium assumptions used in typical wall modeling approaches. The thickness of the velocity BL and of the viscous sublayer was found to scale inversely with engine speed and crank angle. The thermal BL thickness also scales inversely with engine speed but increases with in-cylinder temperature. In contrast, thermal displacement thickness, which is sometimes used as a proxy for thermal BL thickness, was found to decrease with increasing temperature in the bulk. Examination of the heat flux distribution revealed areas of increased heat flux, particularly at places characterized by strong flow directed towards the wall. In addition, significant cyclic variations in the surface-averaged wall heat flux were observed for both engine speeds. An analysis of the cyclic tumble ratio revealed that the cycles with lower tumble ratio values near top dead center (TDC), indicative of an earlier tumble breakdown, also exhibit higher surface averaged wall heat fluxes. These findings extend previous numerical and experimental results for the evolution of BL structure during the compression stroke and serve as an important step for future engine simulations under realistic operating conditions.</p>","PeriodicalId":559,"journal":{"name":"Flow, Turbulence and Combustion","volume":"58 1","pages":""},"PeriodicalIF":2.4,"publicationDate":"2024-08-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141937209","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Experimental and Full-Annulus Simulation Analysis of the Rotating Stall in a Centrifugal Compressor Stage with a Vaned Diffuser","authors":"Yufang Zhang, Shuai Li, Hechun Yu, Linlin Cui","doi":"10.1007/s10494-024-00578-8","DOIUrl":"https://doi.org/10.1007/s10494-024-00578-8","url":null,"abstract":"<p>Flow instability such as rotating stall and even surge occurs when the centrifugal compressor stage operates under low flow conditions. This phenomenon is an extremely complex dynamic process, and it is closely related to the aerodynamic performance and internal flow of the stage. Therefore, it is necessary to study the flow development characteristics in the stage. This paper employs experimental measurement and full-annulus numerical simulation to investigate the effects of diffuser stall on the aerodynamic performance of the compressor and the internal flow of the impeller. The propagation direction, speed, evolution characteristics, and the number of the stall cell were obtained by experimental measurement, and the numerical simulation method was verified. The numerical results that there is a stall limit cycle with counter-clockwise rotation between the flow rate and total pressure ratio of the compressor when the diffuser stalls. Meanwhile, it is found that the stall limit cycle is closely related to the separation strength of the internal flow in the compressor. Finally, the coherent flow structure near the vane shroud side is identified by the modal decomposition methods when the diffuser stalls. The research results in this paper promote an in-depth understanding of the stall mechanism of centrifugal compressors.</p>","PeriodicalId":559,"journal":{"name":"Flow, Turbulence and Combustion","volume":"22 1","pages":""},"PeriodicalIF":2.4,"publicationDate":"2024-08-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141937212","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Spatiotemporal Surface Temperature Measurements Resolving Flame-Wall Interactions of Lean H2-Air and CH4-Air Flames Using Phosphor Thermometry","authors":"Anthony O. Ojo, Abhijit Padhiary, Brian Peterson","doi":"10.1007/s10494-024-00571-1","DOIUrl":"10.1007/s10494-024-00571-1","url":null,"abstract":"<div><p>Spatiotemporal wall temperature (T<sub>wall</sub>) distributions resulting from flame-wall interactions of lean H<sub>2</sub>-air and CH<sub>4</sub>-air flames are measured using phosphor thermometry. Such measurements are important to understand transient heat transfer and wall heat flux associated with various flame features. This is particularly true for hydrogen, which can exhibit a range of unique flame features associated with combustion instabilities. Experiments are performed within a two-wall passage, in an optically accessible chamber. The phosphor ScVO<sub>4</sub>:Bi<sup>3+</sup> is used to measure T<sub>wall</sub> in a 22 × 22 mm<sup>2</sup> region with 180 µm/pixel resolution and repetition rate of 1 kHz. Chemiluminescence imaging is combined with phosphor thermometry to correlate the spatiotemporal dynamics of the flame with the heat signatures imposed on the wall. Measurements are performed for lean H<sub>2</sub>-air flames with equivalence ratio Φ = 0.56 and compared to CH<sub>4</sub>-air flames with Φ = 1. T<sub>wall</sub> signatures for H<sub>2</sub>-air Φ = 0.56 exhibit alternating high and low-temperature vertical streaks associated with finger-like flame structures, while CH<sub>4</sub>-air flames exhibit larger scale wrinkling with identifiable crest/cusp regions that exhibit higher/lower wall temperatures, respectively. The underlying differences in flame morphology and T<sub>wall</sub> distributions observed between the CH<sub>4</sub>-air and lean H<sub>2</sub>-air mixtures are attributed to the differences in their Lewis number (CH<sub>4</sub>-air Φ = 1: Le = 0.94; H<sub>2</sub>-air Φ = 0.56: Le = 0.39). Findings are presented at two different passage spacings to study the increased wall heat loss with larger surface-area-to-volume ratios. Additional experiments are conducted for H<sub>2</sub>-air mixtures with Φ = 0.45, where flame propagation was slower and was more suitable to resolve the wall heat signatures associated with thermodiffusive instabilities. These unstable flame features impose similar wall heat fluxes as flames with 2–3 times greater flame power. In this study, these flame instabilities occur within a small space/time domain, but demonstrate the capability to impose appreciable heat fluxes on surfaces.</p></div>","PeriodicalId":559,"journal":{"name":"Flow, Turbulence and Combustion","volume":"113 4","pages":"1161 - 1188"},"PeriodicalIF":2.0,"publicationDate":"2024-08-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10494-024-00571-1.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141937210","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Antony Premkumar, Francesca Loffredo, Heinz Pitsch, Markus Klein
{"title":"Towards Direct Numerical Simulations of Reactivity-Controlled Compression Ignition Engine Using n-Octanol/Ethanol Fuel Blends","authors":"Antony Premkumar, Francesca Loffredo, Heinz Pitsch, Markus Klein","doi":"10.1007/s10494-024-00570-2","DOIUrl":"https://doi.org/10.1007/s10494-024-00570-2","url":null,"abstract":"<p>The results of a two-dimensional direct numerical simulation of a lean n-octanol-ethanol fuel blend under Reactivity Controlled Compression Ignition (RCCI) conditions are presented in this paper. Stratified temperature and high reactivity fuel fields of Gaussian, bi-modal, and log-normal distributions are studied for uncorrelated and correlated scenarios. The pimple loop is made to run twice to achieve compression heating. A chemical mechanism with 171 species and 734 reactions was developed to capture autoignition characteristics and flame propagation reasonably well. Diagnosing techniques published in the literature are used to determine whether the flame fronts are spontaneously propagating or not. As reported previously for other fuel blends under RCCI conditions, both deflagration and spontaneous ignition flame fronts are observed to co-exist. Gaussian, bi-modal, and log-normal fields respectively move towards a spontaneously igniting mode. Correlating temperature and high reactivity fuel fields not only makes combustion more spontaneously igniting but also more premixed. The analysis reveals the sensitivity of the DNS results with respect to the initial conditions which accordingly should be chosen with care.</p>","PeriodicalId":559,"journal":{"name":"Flow, Turbulence and Combustion","volume":"3 1","pages":""},"PeriodicalIF":2.4,"publicationDate":"2024-07-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141869098","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Spatial Averaging Effects in Adverse Pressure Gradient Turbulent Boundary Layers","authors":"Fermin Mallor, Ramis Örlü, Philipp Schlatter","doi":"10.1007/s10494-024-00568-w","DOIUrl":"https://doi.org/10.1007/s10494-024-00568-w","url":null,"abstract":"<p>Thermal anemometry sensors for time-resolved velocity measurements average the measured signal over the length of their sensor, thereby attenuating fluctuations stemming from scales smaller than the wire length. Several compensation methods have emerged for wall turbulence, the most prominent ones relying on the small-scale universality in canonical flows or on the reconstruction based on two attenuated variance profiles obtained with sensors of different length. To extend these methods to non-canonical flows, the present work considers various adverse-pressure gradient (APG) turbulent boundary layer (TBL) flows in order to explore how the small-scale energy is affected in the inner and outer layer and how the two prominent correction methods perform as function of wall-distance, wire length and flow condition. Our findings show that the increased levels of small-scale energy in the inner, but also outer layer associated with APG TBLs reduces the applicability of empirical methods based on the universality of small-scale energy. On the other hand, a correction based on the relationship between the spanwise Taylor microscale and the two-point streamwise velocity correlation function, is able to correct the attenuated profiles of non-canonical cases. Combining the strength of both methods, a composite profile for the spanwise Taylor microscale is suggested, which then is used for the correction of probe-length attenuation effects across a multitude of flow conditions.</p>","PeriodicalId":559,"journal":{"name":"Flow, Turbulence and Combustion","volume":"45 1","pages":""},"PeriodicalIF":2.4,"publicationDate":"2024-07-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141741238","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Nikolaos Papafilippou, Francesco Pignatelli, Arman Ahamed Subash, Muhammad Aqib Chishty, Rikard Gebart
{"title":"LES of Biomass Syngas Combustion in a Swirl Stabilised Burner: Model Validation and Predictions","authors":"Nikolaos Papafilippou, Francesco Pignatelli, Arman Ahamed Subash, Muhammad Aqib Chishty, Rikard Gebart","doi":"10.1007/s10494-024-00558-y","DOIUrl":"10.1007/s10494-024-00558-y","url":null,"abstract":"<div><p>In this work, numerical investigations were performed using large eddy simulations and validated against detailed measurements in the CeCOST swirl stabilised burner. Both cold and reactive flow have been studied and the model has shown a good agreement with experiments. The verification of the model was done using the LES index of quality and a single grid estimator. The cold flow simulations predicted results closely to experiments setting baseline for the reactive simulations. Coherent structures like the vortex rope above the swirler and a precessing vortex core in the combustion chamber were identified. The reactive conditions were modelled with the Flamelet generated manifold and artificially thickened flame models. Simulations were performed for an experimental syngas composition from black liquor gasification at three different CO<sub>2</sub> dilution levels. Three different Reynolds numbers were investigated with the model matching closely to experimentally detected 2D flow field and OH for the most CO<sub>2</sub> diluted mixture. It was found that the opening angles of the flames differ by a maximum of 13% between experiments and simulations. The most diluted fuel investigated experienced a liftoff distance of 23.5 mm at the Re 25 k. This was also the highest liftoff distance experienced in this cohort of fuels. The same fuel also proved to have the thickest flame annulus at 78.5 mm. Overall, in cases with no experimental data available the predictions made by the model follow the same trends which hints its applicability to higher Re cases.</p></div>","PeriodicalId":559,"journal":{"name":"Flow, Turbulence and Combustion","volume":"113 4","pages":"1189 - 1214"},"PeriodicalIF":2.0,"publicationDate":"2024-07-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10494-024-00558-y.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141746132","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A Comparison of Evaluation Methodologies of the Fractal Dimension of Premixed Turbulent Flames in 2D and 3D Using Direct Numerical Simulation Data","authors":"Marco Herbert, Nilanjan Chakraborty, Markus Klein","doi":"10.1007/s10494-024-00560-4","DOIUrl":"10.1007/s10494-024-00560-4","url":null,"abstract":"<div><p>A Direct Numerical Simulation (DNS) database of statistically planar flames ranging from the wrinkled flamelets to the thin reaction zones regime and DNS data for a Bunsen premixed flame representing the wrinkled flamelets regime have been utilised to evaluate the fractal dimensions of flame surfaces using the filtering dimension method, the box-counting algorithm and the correlation dimension approach. The fractal dimension evaluated based on the fully resolved three-dimensional data has been found to be reasonably approximated by adding unity to the equivalent fractal dimension evaluated based on two-dimensional projections irrespective of the methodology of extracting fractal dimension. This indicates that the flame surface can be approximated as a self-similar fractal surface for the range of Karlovitz and Damköhler numbers considered here. While all methods, provide results identical to each other for benchmark problems, it has been found that the fractal dimension evaluation based on box-counting method provides almost identical results as that obtained using the filtering dimension method for both three and two dimensions, while the fractal dimensions based on the correlation dimension tend to be slightly smaller. The findings of the current analysis have the potential to be used to reliably estimate the actual fractal dimension in 3D based on experimentally obtained 2D binarised reaction progress variable field. The inner cut-off scales estimated based on all three methodologies yield comparable results in terms of order of magnitude with the box-counting method predicting a smaller value of inner cut-off scale in comparison to other methods. The execution times for fractal dimension extraction based on filtering dimension and box-counting methodologies are found to be comparable but the correlation dimension method is found to be considerably faster than the two alternative approaches and provides results consistent with theoretical bounds in all cases.</p></div>","PeriodicalId":559,"journal":{"name":"Flow, Turbulence and Combustion","volume":"113 4","pages":"1145 - 1160"},"PeriodicalIF":2.0,"publicationDate":"2024-07-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10494-024-00560-4.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141721988","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Instability Modes and Scaling Analysis During Electro-Hydro-Dynamic-Atomization: Theoretical and Experimental Study","authors":"Alok Kumar Ray","doi":"10.1007/s10494-024-00567-x","DOIUrl":"10.1007/s10494-024-00567-x","url":null,"abstract":"<div><p>The electro-hydro-dynamic-atomization (EHDA) is a well-established technology with numerous micro/nanoparticle fabrication applications. However, a consistent method for explaining the physics behind the process has yet to be established. The present study aims to report a comprehensive non-dimensional analysis to develop a correlation between different process parameters. The dimensionless numbers derived from Buckingham’s pi theorem match well with those derived from the Navier–Stokes equation, establishing the forces involved in EHDA. Flow instability modes during the EHDA process are experimentally visualized using the flow visualization technique and characterized using a microscope. The instability modes are described using derived non-dimension numbers, and results closely align with Ganan-Calvo’s findings. Derived scaling for the current is in good agreement with Ganan-Calvo (1997), which complies with the condition if δ<sub>μ</sub> × (Q/Qo)<sup>1/3</sup> > > 1, then I/Io = 11 × (Q/Qo)<sup>1/4</sup> -5. Moreover, the ratio of ln (Ehd)/ ln (Md) in cone jet mode is found to be ≈2, irrespective of fluids.</p></div>","PeriodicalId":559,"journal":{"name":"Flow, Turbulence and Combustion","volume":"113 4","pages":"947 - 974"},"PeriodicalIF":2.0,"publicationDate":"2024-07-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141647399","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Assessment of Wall Modeling With Adverse Pressure Gradient for High Reynolds Number Separated Flows","authors":"Sajad Mozaffari, Jérôme Jacob, Pierre Sagaut","doi":"10.1007/s10494-024-00562-2","DOIUrl":"10.1007/s10494-024-00562-2","url":null,"abstract":"<div><p>This paper applies a recently developed approach for modeling turbulence near wall regions within a lattice Boltzmann solver, in combination with a Hybrid RANS/LES turbulence model, to study turbulent separated flows at high Reynolds numbers. To simulate unsteady detached flows on a non-body-fitted Cartesian grid, wall models are employed to estimate the effects of unresolved near-wall turbulence on the overall flow. The article presents the extension of an equilibrium power law wall model to handle adverse pressure gradients and its application in simulating external aerodynamic flows. Hybrid RANS/LES simulations are conducted for two challenging test cases: a 3D NACA-4412 airfoil near stall and a complex Ahmed body configuration. Comparison with a reference simulation involving resolved boundary layers and experimental data demonstrates the strong performance of the wall model, when considering adverse pressure gradients, in simulating turbulent boundary layers under various conditions, ranging from fully attached to mild to high adverse pressure gradients.</p></div>","PeriodicalId":559,"journal":{"name":"Flow, Turbulence and Combustion","volume":"113 4","pages":"923 - 945"},"PeriodicalIF":2.0,"publicationDate":"2024-07-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141648305","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Maximilian Bambauer, Michael Pfitzner, Markus Klein
{"title":"LES of Premixed Turbulent Combustion Using Filtered Tabulated Chemistry","authors":"Maximilian Bambauer, Michael Pfitzner, Markus Klein","doi":"10.1007/s10494-024-00563-1","DOIUrl":"10.1007/s10494-024-00563-1","url":null,"abstract":"<div><p>The filtered tabulated chemistry (FTACLES) approach utilizes data from pre-tabulated explicitly filtered 1D flame profiles for closure of the LES-filtered transport terms. Different methodologies are discussed to obtain a suitable progress variable <i>c</i> from detailed chemistry calculations of a methane/air flame. In this context, special focus is placed on the analytical modeling of the reaction source term using series of parameterized Gaussians. For increasing effective filter sizes in LES (i.e. including the flame thickening) the precise shape of the reaction rate profile becomes less and less relevant. In particular, it is shown that for one-step chemistry, a single Gaussian is sufficient to derive an explicitly expressible 1D flame profile with a prescribed laminar flame speed and thermal flame thickness. The resulting artificial flame profile is shown to have similarities with profiles based on carbon chemistry and detailed reaction mechanisms. Next, the behavior of the filtered <i>c</i>-transport equation is analyzed and several possible closure methods are compared for a wide range of filter widths. It is shown that the unclosed contribution of the filtered diffusion term can be combined with the subgrid convection term, thus simplifying the FTACLES formulation. The model is implemented in OpenFOAM and validated in 1D for a variety of LES filter sizes in combination with artificial flame thickening. A power-law-based wrinkling model is modified for use with artificial flame thickening and combined with the FTACLES model to enable 3D simulations of a premixed turbulent Bunsen burner. The comparison of 3D Large Eddy Bunsen flame simulations at increasing levels of turbulence intensity shows a good match to experimental results for most investigated cases. In addition, the results are mostly insensitive to the variation of the mesh size.</p></div>","PeriodicalId":559,"journal":{"name":"Flow, Turbulence and Combustion","volume":"113 4","pages":"1111 - 1143"},"PeriodicalIF":2.0,"publicationDate":"2024-07-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10494-024-00563-1.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141570184","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}