Estimation of Entropy Generation in a SCR-DeNOx System with AdBlue Spray Dynamic Using Large Eddy Simulation

In this work, the entropy generation analysis is extended to the multi-phase fluid flow within a Large Eddy Simulation (LES) framework. The selected study case consists of a generic selective catalytic reduction (SCR) configuration in which the water/AdBlue is injected into a cross-flow of the internal combustion (IC) engine exhaust gas. The adopted numerical modules are first assessed by comparing with experimental data for film thickness in the case of AdBlue injection and then with H₂O mass fraction and temperature for water injection case. Subsequently, the impact of heat transfer, fluid flow, phase change, mixing and chemical reaction due to AdBlue injection on the entropy generation is assessed. Hence, the individual contributions of viscous and heat dissipation together with the species mixing, chemical reaction during the thermal decomposition of urea into NH₃ and dispersed phase are especially evaluated and analysed. In comparison to the shares of the viscous and mixing processes, the entropy generation is predominated by the heat, chemical and dispersed phase contributions. The influence of the operating parameters such as exhaust gas temperature, flow rate and AdBlue injection on entropy generation is discussed in details. Using a suitable measures, the irreversibility map and some necessary inferences are also provided

The Exergy Losses Analysis in Adiabatic Combustion Systems including the Exhaust Gas Exergy

The entropy generation analysis of adiabatic combustion systems was performed to quantify the exergy losses which are mainly the exergy destroyed during combustion inside the chamber and in the exhaust gases. The purpose of the present work was therefore: (a) to extend the exergy destruction analysis by including the exhaust gas exergy while applying the hybrid filtered Eulerian stochastic field (ESF) method coupled with the FGM chemistry tabulation strategy; (b) to introduce a novel method for evaluating the exergy content of exhaust gases; and (c) to highlight a link between exhaust gas exergy and combustion emissions. In this work, the adiabatic Sandia flames E and F were chosen as application combustion systems. First, the numerical results of the flow and scalar fields were validated by comparison with the experimental data. The under-utilization of eight stochastic fields (SFs), the flow field results and the associated scalar fields for the flame E show excellent agreement contrary to flame F. Then, the different exergy losses were calculated and analyzed. The heat transfer and chemical reaction are the main factors responsible for the exergy destruction during combustion. The chemical exergy of the exhaust gases shows a strong relation between the exergy losses and combustion emission as well as the gas exhaust temperature

Thermal Decomposition of a Single AdBlue® Droplet Including Wall–Film Formation in Turbulent Cross-Flow in an SCR System

Apart from electric vehicles, most internal combustion (IC) engines are powered while burning petroleum-based fossil or alternative fuels after mixing with inducted air. Thereby the operations of mixing and combustion evolve in a turbulent flow environment created during the intake phase and then intensified by the piston motion and influenced by the shape of combustion chamber. In particular, the swirl and turbulence levels existing immediately before and during combustion affect the evolution of these processes and determine engine performance, noise and pollutant emissions. Both the turbulence characteristics and the bulk flow pattern in the cylinder are strongly affected by the inlet port and valve design. In the present paper, large eddy simulation (LES) is appraised and applied to studying the turbulent fluid flow around the intake valve of a single cylinder IC-engine as represented by the so called magnetic resonance velocimetry (MRV) flow bench configuration with a relatively large Reynolds number of 45,000. To avoid an intense mesh refinement near the wall, various subgrid scale models for LES; namely the Smagorinsky, wall adapting local eddy (WALE) model, SIGMA, and dynamic one equation models, are employed in combination with an appropriate wall function. For comparison purposes, the standard RANS (Reynolds-averaged Navier–Stokes) k- ε model is also used. In terms of a global mean error index for the velocity results obtained from all the models, at first it turns out that all the subgrid models show similar predictive capability except the Smagorinsky model, while the standard k- ε model experiences a higher normalized mean absolute error (nMAE) of velocity once compared with MRV data. Secondly, based on the cost-accuracy criteria, the WALE model is used with a fine mesh of ≈39 millions control volumes, the averaged velocity results showed excellent agreement between LES and MRV measurements, revealing the high prediction capability of the suggested LES tool for valve flows. Thirdly, the turbulent flow across the valve curtain clearly featured a back flow resulting in a high speed intake jet in the middle. Comprehensive LES data are generated to carry out statistical analysis in terms of (1) evolution of the turbulent morphology across the valve passage relying on the flow anisotropy map, (2) integral turbulent scales along the intake-charge stream, (3) turbulent flow properties such as turbulent kinetic energy, turbulent velocity and its intensity within the most critical zone in intake-port and along the port length, it further transpires that the most turbulence are generated across the valve passage and these are responsible for the in-cylinder turbulence

Analysis of Local Exergy Losses in Combustion Systems Using a Hybrid Filtered Eulerian Stochastic Field Coupled with Detailed Chemistry Tabulation: Cases of Flames D and E

A second law analysis in combustion systems is performed along with an exergy loss study by quantifying the entropy generation sources using, for the first time, three different approaches: a classical-thermodynamics-based approach, a novel turbulence-based method and a look-up-table-based approach, respectively. The numerical computation is based on a hybrid filtered Eulerian stochastic field (ESF) method coupled with tabulated detailed chemistry according to a Famelet-Generated Manifold (FGM)-based combustion model. In this work, the capability of the three approaches to capture the effect of the Re number on local exergy losses is especially appraised. For this purpose, Sandia flames D and E are selected as application cases. First, the validation of the computed flow and scalar fields is achieved by comparison to available experimental data. For both flames, the flow field results for eight stochastic fields and the associated scalar fields show an excellent agreement. The ESF method reproduces all major features of the flames at a lower numerical cost. Next, the second law analysis carried out with the different approaches for the entropy generation computation provides comparable quantitative results. Using flame D as a reference, for which some results with the thermodynamic-based approach exist in the literature, it turns out that, among the sources of exergy loss, the heat transfer and the chemical reaction emerge notably as the main culprits for entropy production, causing 50% and 35% of it, respectively. This fact-finding increases in Sandia flame E, which features a high Re number compared to Sandia flame D. The computational cost is less once the entropy generation analysis is carried out by using the Large Eddy Simulation (LES) hybrid ESF/FGM approach together with the look-up-table-based or turbulence-based approach

Prediction of Heat Transfer and Fluid Flow Effects on Entropy Generation in a Monolithic Catalytic Converter Using Large-Eddy Simulation

In the present work, heat transfer and fluid flow and their effects on entropy generation in a realistic catalytic converter of a Lada Niva 21214 vehicle are studied using large eddy simulation. At first, the pressure drop over the catalytic converter is measured for dry air at constant temperature (T=298 K), different volumetric flow rates, and extrapolated to large volumetric flow rates for dry air (T=298 K) and for the exhaust gas under realistic engine conditions (T=900 K) using the Darcy–Forchheimer relation. Then, coupled heat and fluid flow phenomena inside the catalytic converter are analyzed for nonreacting isothermal conditions and nonreacting conditions with conjugate heat transfer by using the large-eddy simulation. The predicted pressure drop agrees well with the measured and extrapolated data. Based on the obtained numerical results, the characteristic flow features are identified, namely: the impinging flow with stagnation, recirculation, flow separation and laminarization within the fine ducts of the monolith, which depends on the heat transfer through temperature-dependent thermophysical properties of exhaust gas. Moreover, due to high-velocity gradients at the wall of the narrow ducts in the monolith, entropy production by viscous dissipation is observed predominantly in the monolith region. In contrast, entropy production due to heat transport is relatively small in the monolith region, while it overwhelms viscous dissipation effects in the pipe regions

Database of Near-Wall Turbulent Flow Properties of a Jet Impinging on a Solid Surface under Different Inclination Angles

In the present paper, direct numerical simulation (DNS) and particle image velocimetry (PIV) have been applied complementarily in order to generate a database of near-wall turbulence properties of a highly turbulent jet impinging on a solid surface under different inclination angles. Thereby, the main focus is placed on an impingement angle of 45°, since it represents a good generic benchmark test case for a wide range of technical fluid flow applications. This specific configuration features very complex flow properties including the presence of a stagnation point, development of the shear boundary layer and strong streamline curvature. In particular, this database includes near-wall turbulence statistics along with mean and rms velocities, budget terms in the turbulent kinetic energy equation, anisotropy invariant maps, turbulent length/time scales and near-wall shear stresses. These properties are useful for the validation of near-wall modeling approaches in the context of Reynolds-averaged Navier–Stokes (RANS) and large-eddy simulations (LES). From this study, in which further impingement angles (0°, 90°) have been considered in the experiments only, it turns out that (1) the production of turbulent kinetic energy appears negative at the stagnation point for an impingement angle other than 0° and is balanced predominantly by pressure-related diffusion, (2) quasi-coherent thin streaks with large characteristic time scales appear at the stagnation region, while the organization of the flow is predominantly toroidal further downstream, and (3) near-wall shear stresses are low at the stagnation region and intense in regions where the direction of the flow changes suddenly

Numerical Investigation of Flow through a Valve during Charge Intake in a DISI -Engine Using Large Eddy Simulation

Apart from electric vehicles, most internal combustion (IC) engines are powered while burning petroleum-based fossil or alternative fuels after mixing with inducted air. Thereby the operations of mixing and combustion evolve in a turbulent flow environment created during the intake phase and then intensified by the piston motion and influenced by the shape of combustion chamber. In particular, the swirl and turbulence levels existing immediately before and during combustion affect the evolution of these processes and determine engine performance, noise and pollutant emissions. Both the turbulence characteristics and the bulk flow pattern in the cylinder are strongly affected by the inlet port and valve design. In the present paper, large eddy simulation (LES) is appraised and applied to studying the turbulent fluid flow around the intake valve of a single cylinder IC-engine as represented by the so called magnetic resonance velocimetry (MRV) flow bench configuration with a relatively large Reynolds number of 45,000. To avoid an intense mesh refinement near the wall, various subgrid scale models for LES; namely the Smagorinsky, wall adapting local eddy (WALE) model, SIGMA, and dynamic one equation models, are employed in combination with an appropriate wall function. For comparison purposes, the standard RANS (Reynolds-averaged Navier–Stokes) k- ε model is also used. In terms of a global mean error index for the velocity results obtained from all the models, at first it turns out that all the subgrid models show similar predictive capability except the Smagorinsky model, while the standard k- ε model experiences a higher normalized mean absolute error (nMAE) of velocity once compared with MRV data. Secondly, based on the cost-accuracy criteria, the WALE model is used with a fine mesh of ≈39 millions control volumes, the averaged velocity results showed excellent agreement between LES and MRV measurements, revealing the high prediction capability of the suggested LES tool for valve flows. Thirdly, the turbulent flow across the valve curtain clearly featured a back flow resulting in a high speed intake jet in the middle. Comprehensive LES data are generated to carry out statistical analysis in terms of (1) evolution of the turbulent morphology across the valve passage relying on the flow anisotropy map, (2) integral turbulent scales along the intake-charge stream, (3) turbulent flow properties such as turbulent kinetic energy, turbulent velocity and its intensity within the most critical zone in intake-port and along the port length, it further transpires that the most turbulence are generated across the valve passage and these are responsible for the in-cylinder turbulence

Numerical Investigation of AdBlue Droplet Evaporation and Thermal Decomposition in the Context of NOx-SCR Using a Multi-Component Evaporation Model

To cope with the progressive tightening of the emission regulations, gasoline and diesel engines will continuously require highly improved exhaust after-treatment systems. In the case of diesel engines, the selective catalytic reduction (SCR) appears as one of the widely adopted technologies to reduce NOx (nitrogen oxides) emissions. Thereby, with the help of available heat from exhaust gas, the injected urea–water solution (UWS) turns inside the exhaust port immediately into gaseous ammonia (NH3) by evaporation of mixture and thermal decomposition of urea. The reaction and conversion efficiency mostly depend upon the evaporation and subsequent mixing of the NH3 into the exhaust gas, which in turn depends upon the engine loading conditions. Up to now, the aggregation of urea after evaporation of water and during the thermal decomposition of urea is not clearly understood. Hence, various scenarios for the urea depletion in the gaseous phase that can be envisaged have to be appraised under SCR operating conditions relying on an appropriate evaporation description. The objective of the present paper is therefore fourfold. First, a reliable multi-component evaporation model that includes a proper binary diffusion coefficient is developed for the first time in the Euler–Lagrangian CFD (computational fluid dynamics) framework to account properly for the distinct evaporation regimes of adBlue droplets under various operating conditions. Second, this model is extended for thermal decomposition of urea in the gaseous phase, where, depending on how the heat of thermal decomposition of urea is provided, different scenarios are considered. Third, since the evaporation model at and around the droplet surface is based on a gas film approach, how the material properties are evaluated in the film influences the process results is reported, also for the first time. Finally, the impact of various ambient temperatures on the adBlue droplet depletion characteristics as well as the effect of gravity is pointed out. The prediction capability of the model variants is assessed by comparing the achieved results to each other and with experimental data. It turns out that satisfactory agreement between experiment and numerical predictions is achieved for a wide range of operating temperatures by using correlations by “Wilke and Lee” for urea and by “Fuller et al.” for water. The results are essentially sensitive to gravity. From subsequent comparisons of different ways to account for the thermal decomposition in the gaseous urea, a significant difference is observed. Finally, the 1/3 film rule widely used for evaluating the material properties in the film shows accurate prediction of both evaporation and thermal decomposition regimes of urea

Modeling and unsteady simulation of turbulent multi-phase flow including fuel injection in IC-engines

In internal combustion engine (ICE), researchers have to face with stringent environment regulations concerning pollutants while improving engine thermal efficiency, making the engine design a complex task. To meet these requirements, an understanding of the salient features of all the engine processes are very important. Being the primitive process of engine operations, fuel injection influences whole engine cycle via fuel-air mixture preparation, thereby the combustion behavior and subsequently the emission performance. The inhospitable environment inside a combustion chamber makes the experimental investigations more complex and expensive. In contrast, a CFD based investigation can provide comprehensive insight about in-cylinder flow field, spray injection phenomena as encountered in IC-engine. In the present study, a CFD tool that enables to investigate the real unsteady behavior of realistic engine configuration is developed by coupling Large Eddy Simulation (LES) together with a spray module using the KIVA4-mpi Code. It is based on an Eulerian-Lagrangian framework to describe the spray evolution including primary and secondary atomization. A linear instability sheet atomization (LISA) based sub-model is integrated to represent the primary atomization. The secondary atomization is modeled by an available Taylor analogy break-up (TAB) model. In dense spray region, the droplet-droplet interaction considerably influences the overall spray dynamics. The first novelty of the proposed methodology is to include droplet-droplet interaction processes via an appropriate collision sub-model that is independent of mesh size and type. Thereby, taking account of different regimes, such as bouncing, separation, stretching separation, reflective separation and coalescence. The formation of wall film on hot cylinder surface is a critical process in an IC-engine, since it largely influences the engine performance and emission characteristics. The second novelty of this spray module is the implementation of an improved wall film model that includes the combined effects of droplet kinetic energy and wall temperature into KIVA4-mpi code. To perform an IC-engine simulation, a good quality mesh generation in ICEM-CFD for an engine geometry is challenging task. The KIVA4-mpi is compatible only with block structured mesh without any use of O-grid. Due to this reason, only certain degree of mesh refinement is possible. This makes it difficult to achieve a good quality fine mesh required for LES simulation. In the present study, a new meshing strategy is proposed to generate suitable mesh for real IC-engine configurations. The new method clearly demonstrates the improvement in resolving the in-cylinder flow structures. First, the simulated results for motored case (no fuel injection and no combustion) are compared with the experimental data for a transparent combustion chamber (TCC) engine configuration from Engine Combustion Network (ECN). Second, to demonstrate the importance of fuel injection sub-models, further simulations are carried out including the evolution of evaporating fuel spray with wall impingement. Third, using the new meshing strategy, simulations are also performed for a real complex canted 4-valve engine configuration. Simulated results are compared well with available experimental data