Numerical Investigation of SCR Mixer Design Optimization for Improved Performance

The continuous increase in the number of stringent exhaust emission legislations of marine Diesel engines had led to a decrease in NOx emissions at the required level. Selective catalyst reduction (SCR) is the most prominent and mature technology used to reduce NOx emissions. However, to obtain maximum NOx removal with minimum ammonia slip remains a challenge. Therefore, new mixers are designed in order to obtain the maximum SCR efficiency. This paper reports performance parameters such as uniformity of velocity, ammonia uniformity distribution, and temperature distribution. Also, a numerical model is developed to investigate the interaction of urea droplet with exhaust gas and its effects by using line (LM) and swirl (SM) type mixers alone and in combination (LSM). The urea droplet residence time and its interaction in straight pipe are also investigated. Model calculations proved the improvement in velocity uniformity, distribution of ammonia uniformity, and temperature distribution for LSM. Prominent enhancement in the evaporation rate was also achieved by using LSM, which may be due to the breaking of urea droplets into droplets of smaller diameter. Therefore, the SCR system accomplished higher urea conversion efficiency by using LSM. Lastly, the ISO 8178 standard engine test cycle E3 was used to verify the simulation results. It has been observed that the average weighted value of NOx emission obtained at SCR outlet using LSM was 2.44 g/kWh, which strongly meets International Maritime Organization (IMO) Tier III NOx (3.4 g/kWh) emission regulations.


Introduction
Environmental safety is one of the hottest research areas nowadays due to increased public awareness. To achieve this goal, efforts are continuously made to reduce pollution and develop green processes. Exhaust emissions from marine Diesel engine are responsible for producing severe environmental pollution, especially nitrogen oxide (NO x ) emissions [1,2]. The automobile Diesel engine produces exhaust emissions species such as carbon monoxide, carbon dioxide, and hydrocarbon in abundance compared to NO x . On the other hand, marine Diesel engine produces more detrimental NO x emissions [3]. Hence, in order to reduce exhaust emissions from ships, many national and international organizations have promulgated regulations and also enforced strict requirements on NO x emissions in Emission Control Areas [4]. In 2016, International Maritime organization (IMO) Tier III regulations on NO x emissions have already been enforced in North America emission control areas, including the East and West Coast of USA and Caribbean. As reported, it will also be enforced in North Sea and Baltic Sea in future [5]. Both high pressure fuel injection and exhaust gas recirculation systems have potential to reduce NO x emissions, but due to poor engine performance results and continuous increase in engine emission legislations, more primitive and improved processes are needed to overcome the said issues [6]. The most convenient and easy option is the treatment of exhaust gas. Selective Catalyst Reduction (SCR) is technically mature, and is the most prominent after-treatment technology used to meet the latest NO x emission regulations due to its high NO x removal efficiency, cost-effectiveness, and good fuel economy [7]. However, challenges related with SCR system include improper mixture of urea water solution (UWS) with exhaust gas and ammonia leakage [8]. UWS (32.5% urea) is injected into exhaust gas [9,10]. Urea decomposition occurs in three steps: firstly, water is evaporated from urea water solution; secondly, the urea pyrolysis reaction occurs, which results in the decomposition of urea into isocyanic acid (HNCO) and ammonia (NH 3 ); and lastly, hydrolysis of HNCO occurs which produces NH 3 and CO 2 [11,12]. Ammonia is used as a reducing agent and cannot be used directly due to poisonous nature, storing and handling difficulties [13]. It has been proved that NO x is mainly 90% composition of NO at the exhaust of the marine Diesel engine. The key reactions involved in SCR system are described as follows [14][15][16][17]. CO A static mixer is commonly adopted to generate uniform distribution of ammonia at the inlet of SCR catalyst [18]. CFD code (Fire 8. 3,2004) has been used to optimize the design of SCR system. The authors have done many studies about the decomposition and evaporation of urea water solution without using static mixer. Birkhold [19] studied the urea droplet's evaporation and decomposition at different temperatures of exhaust gas. Strom [20] studied the effects of turbulent velocity and different forces on distribution and movements of urea water droplet in the exhaust species. In addition, other authors have established various studies to evaluate the effect of static mixer on the SCR performance. Sivanandi Rajadurai [21] studied the distribution of ammonia by using wire mesh mixer. Zhang [22] investigated the uniformity index of ammonia by adopting delta wing mixer. Shazam Williams and Ming Chen studied the velocity and ammonia distribution in the straight pipe together with static mixer having two rows of tabs [23]. Azael and Ibarra investigated the interaction of fluid with a double vortex mixer to improve the performance of SCR. The evaporation effect and droplet crushing was significantly improved. The distribution of NH 3 was also studied, revealing a smooth and uniform distribution using a double vortex mixer; however, the velocity uniformity distribution, reaction temperature, and droplet distribution time was not studied [24]. One of the disadvantages of the SCR system is that it occupies the additional space on vessels. Mostly, authors are working to simplify the SCR system according to the space requirements. One of the authors used an 18 mm distance between the two mixers without considering the impact of velocity on SCR system [18]. Park et al. reported that if the continuous uneven impact of velocity occurs, which causes excessive aerodynamic velocity and temperature, it will lead to the thermal fatigue failure, which reduces the service life [25]. The first mixer can create turbulence intensity and sudden impact on velocity, which directly affects the second mixer. By considering both of the references, it can be perceived that the optimum condition which can satisfy both of the conclusions from different studies can lie within limits. It therefore looked suitable to trial for 0.2 m distance. One of the study used series of only swirl mixer (SM) near the injection point and in line [18]. However, use of a line mixer (LM) can result in increasing the mixing flow (exhaust gas plus urea droplets) in the blind corner near the pipe wall, while SM results in increasing the mixing flow in the center [25]. Some studies investigated LM with different blade angles [25,26], and others used SM with different blade angles [18,27]. In this study LM was used in front of SM. Initially, LM and SM were used separately and then analyzed for the combined effect of both (LSM). In this work, the blade angle was 45 • for both mixers because it has been proven that static mixers consisting of bitched blades with an angle of inclination angle of 45 • can generate higher turbulence intensity and a swirling flow in the pipe [28]. Furthermore, many authors have used a greater number of blades than this study. One author used 36 blades in only LM [25]. In this study, 18 (LM:12, SM:6) blades were used and distributed in two difference places in the pipe. In addition, if a single mixer has many blades at a certain location, it tends to decrease the wall temperature, which ultimately results in deposit formation [29]. Hence it is necessary to distribute number of blades in two mixers with different locations to increase the mixing performance and prevents deposit formation. Generally, SCR performance depends upon the velocity uniformity and ammonia uniformity near the inlet of the SCR catalyst. Therefore, a suitable design of mixers for increasing the efficiency of SCR is needed. Furthermore, a numerical model has been developed, which describes the impact of mixer on the ammonia and velocity uniformity, droplet residence time, and temperature distribution in the pipe. In addition, it is also necessary to consider reaction temperature and wall temperature, as it affects the catalyst performance and deposit formation [1]. In this paper, two different types of mixers-line (LM) and swirl (SM) type mixers alone and in combination (LSM)-were used to indicate the effects on the performance of SCR. The main objective of this investigation is to create the optimum SCR design, to achieve higher uniformity index for velocity and ammonia distribution, better evaporation rate, and droplet distribution by consideration of the reaction temperature and wall temperature distribution based on CFD. Furthermore, for the verification of the simulated results of ISO 8178 a standard marine Diesel engine test cycle E3 was used [30].

Geometric Model
The complications in the design and arrangements of SCR system leads to difficulties in the uniformity of flow field and ammonia in the exhaust flow by using only urea water solution (UWS) and species diffusion. It is therefore necessary to use static mixer to optimize the mixing of exhaust gas flow with UWS. It has been proved by previous studies that uniformity of flow velocity and ammonia distribution can be increased with a suitable mixer, which directly enhances the performance of SCR catalyst. The optimized mixers have ability of complete rapid mixing of UWS droplet with exhaust gas which accelerates the pyrolysis reaction. Therefore, the catalyst conversion efficiency was improved [31].
In this study, two different types of mixers are used-line type mixer (LM) and swirl type (SM). In the first part of study, the SCR performance was analyzed by LM and SM separately. The second part of study focused on the combined effects of LM and SM. The distance between two mixers was kept 0.2 m. The pipe diameter was 100 mm. Injection point was located 0.1 m away from the LM and 0.3 m away from the SM, when both mixers were used together. For the single mixer, the position was taken as 0.1 m away from the injection point. The complete systematic design of SCR system is represented in Figure 1.
The design of LM includes 4 × 4 rows of blades and the angle of blades with horizontal plane was 45 • . This design represents the improved heat transfer efficiency and abolishes the blind corners. Four blades were used at the center rows while two at the corners. One column of blades shows a downward direction and the other shows an upward direction. The design construction of SM was adopted in order to enhance the breakup effect of droplets and increase the turbulence intensity which ultimately results rapid mixing and pyrolysis of urea. The SM mixer possess six number of blades tilted at an angle of 45 • with the horizontal. By analyzing the result of individual mixers on SCR performance, a system was designed using both LM and SM in combination to achieve an improvement in performance. Considering the complexity of SCR system, the mesh grid is generated by tetrahedral method for adequate adaptation as shown in Figure 2. Away from the solid wall boundary, the grid was stretched which was the good comprise between the mesh points and computational costs. Furthermore, mesh independence test has been performed to calculate the ammonia uniformity without using any mixer, as shown in Figure 3. After mesh independence analysis, the total number of cells in the mesh system was 1.8 million; because a further increased in mesh cells has no significant effect on the uniformity index.
Processes 2019, 7, x FOR PEER REVIEW 4 of 22 ammonia uniformity without using any mixer, as shown in Figure 3. After mesh independence analysis, the total number of cells in the mesh system was 1.8 million; because a further increased in mesh cells has no significant effect on the uniformity index.   ammonia uniformity without using any mixer, as shown in Figure 3. After mesh independence analysis, the total number of cells in the mesh system was 1.8 million; because a further increased in mesh cells has no significant effect on the uniformity index.

Numerical Procedure
Numerical studies have been done for the optimum design of SCR system in the form of flow uniformity and the reduction of NOX emissions [22,32]. In this study, the flow field distribution of velocity was taken uniform and thermal expansion in the catalyst was negligible. Exhaust gases from the engine were adopted as the ideal gas mixture in the CFD. Incompressible steady flow was used. The flow was simulated by using the standard K-epsilon model (K-ε) in combination with wall treatment function to evaluate the turbulent flow velocity in the exhaust system because the Reynolds number is usually very high in an exhaust pipe. It is a two equation model that gives a general description of turbulence by means of two transport partial differential equations. The first transported variable is the turbulence kinetic energy (K) and the second transported variable is the rate of dissipation of turbulence energy (ε) [25,33]. The injection and decomposition process of UWS were simulated by transport model. Evaporation and decomposition performances were described by using a mixture model of multicomponents. The concentration of 32.5 wt% UWS was injected at the temperature window of 313 K. Pressure swirl atomizer was used for injection of UWS. KH-RT model was used for breakup of droplet. Firstly, water was evaporated from urea droplets during decomposition. The rate of water evaporation was predicted by Raouli's law of multi component mixture. Raoult's law is based on the assumptions that the vapor phase behaves as an ideal gas and the liquid phase is an ideal solution in mixing of the flow. The rate of water evaporation can be obtained by Equation (1).
, water sat water water c water P If the fraction of water in UWS is going beyond the threshold value (0.01), then the UWS droplets are considered as the particles of urea and it can be easily calculated by the decomposition reactions of urea particles. As the urea particles react with exhaust gas species, it results in NH3 and HNCO. The mass change rate can be obtained by Equation (2).
The species transport and continuity equation are used to define the variation of gases produced from decomposition reactions such as NH3, CO2, and HNCO (Equations (3) and (4)).

Numerical Procedure
Numerical studies have been done for the optimum design of SCR system in the form of flow uniformity and the reduction of NO X emissions [22,32]. In this study, the flow field distribution of velocity was taken uniform and thermal expansion in the catalyst was negligible. Exhaust gases from the engine were adopted as the ideal gas mixture in the CFD. Incompressible steady flow was used. The flow was simulated by using the standard K-epsilon model (K-ε) in combination with wall treatment function to evaluate the turbulent flow velocity in the exhaust system because the Reynolds number is usually very high in an exhaust pipe. It is a two equation model that gives a general description of turbulence by means of two transport partial differential equations. The first transported variable is the turbulence kinetic energy (K) and the second transported variable is the rate of dissipation of turbulence energy (ε) [25,33]. The injection and decomposition process of UWS were simulated by transport model. Evaporation and decomposition performances were described by using a mixture model of multicomponents. The concentration of 32.5 wt% UWS was injected at the temperature window of 313 K. Pressure swirl atomizer was used for injection of UWS. KH-RT model was used for breakup of droplet. Firstly, water was evaporated from urea droplets during decomposition. The rate of water evaporation was predicted by Raouli's law of multi component mixture. Raoult's law is based on the assumptions that the vapor phase behaves as an ideal gas and the liquid phase is an ideal solution in mixing of the flow. The rate of water evaporation can be obtained by Equation (1). dm dt = MW water k c x water P sat,water RT P − C water (1) If the fraction of water in UWS is going beyond the threshold value (0.01), then the UWS droplets are considered as the particles of urea and it can be easily calculated by the decomposition reactions of urea particles. As the urea particles react with exhaust gas species, it results in NH 3 and HNCO. The mass change rate can be obtained by Equation (2).
The species transport and continuity equation are used to define the variation of gases produced from decomposition reactions such as NH 3 , CO 2 , and HNCO (Equations (3) and (4)).
∂ρ ∂t Processes 2019, 7, 168 6 of 22 By concentration gradients, the mass diffusion in laminar flows used and it was modeled by Fick's law, as shown in Equation (5).
The turbulence due to turbulent diffusion is directly responsible for the quick mixing and species transport in turbulent flows therefore the mass diffusion can be obtained by Equation (6).
The diffusive process is modeled by using turbulent diffusivity (D T ) in the K-ε model. The D T can be obtained from Equation (7) [34].
Schmidt turbulence number (Sc t ) is usually taken as 0.7. In the turbulence model, the turbulent diffusivity D T directly depends upon the use of turbulent viscosity V T . Turbulent viscosity is obtained from given formula in the K-ε model Equation (8).
v T = C u k 2 /ε (8) In the K-ε model, closure coefficient (C u ) is taken as 0.09. Energy dissipation rate and turbulent kinetic energy are obtained from a standard two-layer K-ε turbulent model as showed in Equations (9) and (10), respectively.

∂k ∂t
∂ε ∂t where, For describing the conservation equation of gas-liquid phase and spray phenomenon, the discrete phase model was used. By the discrete phase model under the Lagrangian coordinate system, it could be constructed a governing equation to analyze the motion law of droplet and motion of liquid droplet which can be given by Equation (11).
Gradient diffusion is responsible for the rate of vaporization at low vaporization rates [35]. Transportation of urea droplet vapor was made at the bulk of gas and droplet surface because of the concentration gradient. At high vaporization rates, the fluctuation of urea water droplets into exhaust gases is directly associated with the effects of evaporating species under convective flow [36]. UWS droplets include the mixtures of two parts; therefore, the evaporation of UWS is separated into two different periods. During the first stage of decomposition, contents of water evaporate quickly from UWS droplets. In the later stage, molten urea evaporates [37,38]. In this study, UWS droplets were taken as multicomponent, therefore the total vaporization rate of UWS is the sum of vaporization rates for all UWS components. The evaporation of urea solution is the heat and mass transfers between the urea droplet and exhaust gases. Diffusion is controlled by concentration gradients described by Fick's law which was used as diffusion also occurs in multispecies systems that are experiencing the Stefan flow also and the model included the internal recirculation and Stefan flow effects [39].Thus the evaporation rate can be described by Equation (12).
Radiation properties are very complex to be described by the flow and have less importance. It was therefore neglected in our study. Moreover, some studies also have neglected this effect due to the complexity of system [25,26].
A honeycomb-type catalyst is often used in the SCR system, and its single channel is very small to the outline size of the catalyst [32]. In this study, the catalytic reaction section was taken as a multiple holes section, and the porous medium model was adopted for evaluation. In the porous medium model, the momentum source term was included to the momentum equation and showed in Equation (13).
If each position of catalyst has the same nature in gas convection and diffusion in the channel are ignored, the above equation can be simplified as Equation (14).

Boundary Conditions
CFD boundary conditions were selected according to the ISO standard 8178 marine Diesel engine. The initial temperature and exhaust flow rate before the catalyst reactor weretaken as 327 • C and 1.6 kg/s, respectively. The quantity of NO x emission was 256 ppm. The inlet and outlet diameter of catalyst reactor were 100 mm. The boundary conditions at the inlet of catalytic reactor was set to the mass flow rate, turbulent kinetic energy was set as 6% of the average velocity, and the characteristic length was 10% of the inlet diameter. The boundary condition at the outlet of catalytic reactor was defined as outlet pressure and the outlet pressure was same as atmospheric pressure. The boundary conditions at the wall of catalytic reactor were set as nonslip velocity and frictionless. The material is defined as iron.
The working parameters of the SCR system are shown in Tables 1 and 2. The working process of the SCR system was studied and simulated in this paper. Flow rate, temperature at the diesel exhaust, exhaust pressure, and urea injection rate increased with respect to engine load. Furthermore, urea water solution was injected by using air assisted injector at the upstream of mixer with the ten numbers of streams. Pressure atomizer injector was used with the injection hole diameter of 0.0007 m. The spray angle was kept 60 • with spray pressure of 6 bar. The rate of multitude hole was 0.89 in the porous medium model. The thermal conductivity of solid phase coefficient was taken as 1.7 W/(m·k) and the specific heat capacity was 1016 J/(kg·k).

Uniformity
Catalyst efficiency depends upon the two important parameters: flow velocity uniformity and uniform distribution of ammonia. Long-time irregular scattering and non uniformity of flow velocity will cause the excessive aerodynamic velocity and temperature, which directly impacts the structure and performance of catalyst. As the result, thermal stress can be developed, which leads to the fatigue failure and effects the service life of catalyst [25]. It is therefore recommended that the velocity of UWS must be well distributed as much as possible, at the front of the SCR catalyst reactor inlet. Equation (15) can be used to calculate the uniformity index of flow velocity and ammonia distribution.

Velocity Uniformity
The distribution of velocity field of LM, SM and the combination of both (LSM) is represented in Figure 4. The exhaust gases flow with LM was initially very uneven; downstream of mixer at the center of pipe. Furthermore, a baffle gap was also generated at the downstream but the velocity distribution inside the baffle seems smooth. As the distance was increased to the reactor inlet, uniformity was also increased, as shown in Figure 4a. SM effectively created the rotational air flow generating swirling mixing flow at the pipe wall. Swirl flow shows well mixing of exhaust gases and urea solution. Moreover, a high baffle gap was produced at the downstream of mixer, but the velocity distribution inside the baffle was highly uneven. However, flow velocity of mixture at the middle seems comparatively low as in Figure 4b. Figure 4c showed the combined effects of LM and SM. In comparison with pure rotating blade and line type blade structures of SM and LM, respectively, the velocity uniformity was enhanced when LSM was employed. The arrangement of LSM not only generated a strong swirl flow beside the mixer wall, but also increased the velocity near the wall of pipe.
Catalyst effectiveness and utilization rate directly depends on the velocity uniformity index of UWS and exhaust gas mixture along the axial direction at the catalyst reactor inlet. The relation of the velocity uniformity index from the mixer downstream to catalyst reactor inlet is shown in Figure 5. Uniformity index of flow velocity of LM, SM, and LSM mixers were calculated as 0.93, 0.86, and 0.95 at the upstream of reactor inlet (0.8 m), respectively. Uniformity index of LSM at the downstream of mixers (0.1 m) is relatively low because of continuous and sudden impacts of two mixers on the flow velocity. However, with the increase of distance the uniformity of velocity for LSM was greatly improved due to the turbulence outcome produced from the sudden impact effect of two mixers. Therefore, the overall uniformity index of flow velocity of LSM before the reactor inlet was relatively high as compared to LM and SM.   Figure 4c showed the combined effects of LM and SM. In comparison with pure rotating blade and line type blade structures of SM and LM, respectively, the velocity uniformity was enhanced when LSM was employed. The arrangement of LSM not only generated a strong swirl flow beside the mixer wall, but also increased the velocity near the wall of pipe.
Catalyst effectiveness and utilization rate directly depends on the velocity uniformity index of UWS and exhaust gas mixture along the axial direction at the catalyst reactor inlet. The relation of the velocity uniformity index from the mixer downstream to catalyst reactor inlet is shown in Figure  5. Uniformity index of flow velocity of LM, SM, and LSM mixers were calculated as 0.93, 0.86, and 0.95 at the upstream of reactor inlet (0.8 m), respectively. Uniformity index of LSM at the downstream of mixers (0.1m) is relatively low because of continuous and sudden impacts of two mixers on the flow velocity. However, with the increase of distance the uniformity of velocity for LSM was greatly improved due to the turbulence outcome produced from the sudden impact effect of two mixers. Therefore, the overall uniformity index of flow velocity of LSM before the reactor inlet was relatively high as compared to LM and SM.    Figure 4c showed the combined effects of LM and SM. In comparison with pure rotating blade and line type blade structures of SM and LM, respectively, the velocity uniformity was enhanced when LSM was employed. The arrangement of LSM not only generated a strong swirl flow beside the mixer wall, but also increased the velocity near the wall of pipe.
Catalyst effectiveness and utilization rate directly depends on the velocity uniformity index of UWS and exhaust gas mixture along the axial direction at the catalyst reactor inlet. The relation of the velocity uniformity index from the mixer downstream to catalyst reactor inlet is shown in Figure  5. Uniformity index of flow velocity of LM, SM, and LSM mixers were calculated as 0.93, 0.86, and 0.95 at the upstream of reactor inlet (0.8 m), respectively. Uniformity index of LSM at the downstream of mixers (0.1m) is relatively low because of continuous and sudden impacts of two mixers on the flow velocity. However, with the increase of distance the uniformity of velocity for LSM was greatly improved due to the turbulence outcome produced from the sudden impact effect of two mixers. Therefore, the overall uniformity index of flow velocity of LSM before the reactor inlet was relatively high as compared to LM and SM.  To illuminate the SCR reactor internal stability and uniformity index of flow velocity, the distribution at the reactor entrance was extracted and is shown in Figure 6. The main objective of displaying the velocity uniformity at the reactor inlet was to remove the turbulence effect produced by mixers and to show actual velocity distribution more clearly. Figure 6a,c shows that the exhaust gases gained high turbulent kinetic energy, resulting in the better velocity uniformity. However, SM (Figure 6b) shows that the radial velocity gradient was high and resulted in poor uniformity, but the velocity flow was slightly decreased with the increase of distance. distribution at the reactor entrance was extracted and is shown in Figure 6. The main objective of displaying the velocity uniformity at the reactor inlet was to remove the turbulence effect produced by mixers and to show actual velocity distribution more clearly. Figure6a,c shows that the exhaust gases gained high turbulent kinetic energy, resulting in the better velocity uniformity. However, SM (Figure6b) shows that the radial velocity gradient was high and resulted in poor uniformity, but the velocity flow was slightly decreased with the increase of distance.

Ammonia Uniformity
The irregular and insufficient ammonia distribution will result in NOX escape [40,41]. The NH3 uniformity distribution is an important index of the SCR system for estimating performance [42][43][44]. The water from the urea droplets evaporates entirely and converts urea in NH3gas when it reaches to reactor inlet. Reductant uniformity index is the most significant factor used to evaluate, whether the catalyst will achieve the highest denitrification rate and lowest ammonia leakage. The plane was created in CFD at a certain distance/area and then divided the plane into the number of points and finds the mass fraction of ammonia on each single point. After that, the standard deviation and average value of all the points were calculated. Finally, the formula for the uniformity index was used to find ammonia uniformity at certain distances. The ammonia uniformity index from the mixer downstream to catalyst reactor inlet in the axial direction was calculated and is shown in Figure 7. The NH3uniformity distribution of LM, SM, and LSM were calculated as 0.87, 0.94, and 0.96 at the upstream of reactor inlet, respectively.
Ammonia uniformity index of LSM at the downstream of mixers is relatively low due to uninterruptedly and sudden impacts of two mixers continuously. However, the uniformity index of ammonia of two mixers was greatly enhanced with the increase of distance due to high turbulence outcome produced from the impact effect of two mixers. Hence, generally, the ammonia uniformity index is the combination of two (line and swirl type) mixers before the catalyst reactor inlet is relatively high, as compared with the separately use of line and swirl type mixers.
Uniformity index of ammonia distribution at the outlet of reactor inlet was extracted to clarify the result as much as possible and showed in Figure 8. At first, the plane was created in CFD at the catalyst inlet, plane was divided into number of cells and value of molar concentration of ammonia at each single point was calculated. Next, the standard deviation and average value of ammonia concentration at the catalyst inlet were calculated. Finally, the uniformity index was used to calculate the ammonia uniformity at catalyst inlet. In LM, a high concentration of ammonia

Ammonia Uniformity
The irregular and insufficient ammonia distribution will result in NO X escape [40,41]. The NH 3 uniformity distribution is an important index of the SCR system for estimating performance [42][43][44]. The water from the urea droplets evaporates entirely and converts urea in NH 3 gas when it reaches to reactor inlet. Reductant uniformity index is the most significant factor used to evaluate, whether the catalyst will achieve the highest denitrification rate and lowest ammonia leakage. The plane was created in CFD at a certain distance/area and then divided the plane into the number of points and finds the mass fraction of ammonia on each single point. After that, the standard deviation and average value of all the points were calculated. Finally, the formula for the uniformity index was used to find ammonia uniformity at certain distances. The ammonia uniformity index from the mixer downstream to catalyst reactor inlet in the axial direction was calculated and is shown in Figure 7. The NH 3 uniformity distribution of LM, SM, and LSM were calculated as 0.87, 0.94, and 0.96 at the upstream of reactor inlet, respectively.
Ammonia uniformity index of LSM at the downstream of mixers is relatively low due to uninterruptedly and sudden impacts of two mixers continuously. However, the uniformity index of ammonia of two mixers was greatly enhanced with the increase of distance due to high turbulence outcome produced from the impact effect of two mixers. Hence, generally, the ammonia uniformity index is the combination of two (line and swirl type) mixers before the catalyst reactor inlet is relatively high, as compared with the separately use of line and swirl type mixers.
Uniformity index of ammonia distribution at the outlet of reactor inlet was extracted to clarify the result as much as possible and showed in Figure 8. At first, the plane was created in CFD at the catalyst inlet, plane was divided into number of cells and value of molar concentration of ammonia at each single point was calculated. Next, the standard deviation and average value of ammonia concentration at the catalyst inlet were calculated. Finally, the uniformity index was used to calculate the ammonia uniformity at catalyst inlet. In LM, a high concentration of ammonia distribution is located at the upper side of the pipe wall, and some places at the mid of pipe also show that ammonia distribution was not smooth (Figure 8a), resulting in poor catalyst performance. In SM, ammonia seems well distributed as a whole but the small area at the upper side of pipe shows the high concentration (Figure 8b). Distribution of ammonia was greatly improved when using LSM and much less ammonia was deposited at sides and center (Figure 8c). distribution is located at the upper side of the pipe wall, and some places at the mid of pipe also show that ammonia distribution was not smooth (Figure 8a), resulting in poor catalyst performance. In SM, ammonia seems well distributed as a whole but the small area at the upper side of pipe shows the high concentration (Figure8b). Distribution of ammonia was greatly improved when using LSM and much less ammonia was deposited at sides and center (Figure8c). The values of indexes of the velocity and ammonia uniformity of LM, SM, and LSM from the location of mixer downstream to catalyst reactor inlet are tabulated in Table 3. It is very clear that LM reveals a good uniformity index for velocity but, on the contrary, it exhibits a poor uniformity index for ammonia. SM behaves opposite to the LM: it possesses good ammonia uniformity but shows poor velocity distribution.
It is difficult to obtain the good uniformity of the both important indicators at the same time while using only single mixer. In comparison, LSM can be categorized with the good uniformity for both, the velocity and ammonia distribution due to high turbulent intensity. The values calculated distribution is located at the upper side of the pipe wall, and some places at the mid of pipe also show that ammonia distribution was not smooth (Figure 8a), resulting in poor catalyst performance. In SM, ammonia seems well distributed as a whole but the small area at the upper side of pipe shows the high concentration (Figure8b). Distribution of ammonia was greatly improved when using LSM and much less ammonia was deposited at sides and center (Figure8c). The values of indexes of the velocity and ammonia uniformity of LM, SM, and LSM from the location of mixer downstream to catalyst reactor inlet are tabulated in Table 3. It is very clear that LM reveals a good uniformity index for velocity but, on the contrary, it exhibits a poor uniformity index for ammonia. SM behaves opposite to the LM: it possesses good ammonia uniformity but shows poor velocity distribution.
It is difficult to obtain the good uniformity of the both important indicators at the same time while using only single mixer. In comparison, LSM can be categorized with the good uniformity for both, the velocity and ammonia distribution due to high turbulent intensity. The values calculated The values of indexes of the velocity and ammonia uniformity of LM, SM, and LSM from the location of mixer downstream to catalyst reactor inlet are tabulated in Table 3. It is very clear that LM reveals a good uniformity index for velocity but, on the contrary, it exhibits a poor uniformity index for ammonia. SM behaves opposite to the LM: it possesses good ammonia uniformity but shows poor velocity distribution. It is difficult to obtain the good uniformity of the both important indicators at the same time while using only single mixer. In comparison, LSM can be categorized with the good uniformity for both, the velocity and ammonia distribution due to high turbulent intensity. The values calculated for uniformity were 0.95 and 0.96 for the velocity and ammonia, respectively. For the purpose of mass transfer and homogenization, when the mixture species passes through the second mixer of LSM, it creates swirl and endorse the mass transfer effect near the pipe wall. The shear and swirl flows with the pipe wall produce molecular diffusion and eddy diffusion under forced convection. Hence, the mixing of ammonia with exhaust gases is comparatively sufficient which results the effective improvement in the uniformity of velocity and ammonia, simultaneously.

Droplet Residence Time
UWS injection in the exhaust gases lowers the temperature of the gas phase due to heat transfer phenomena. Figure 9 represents the UWS droplets residence time for diameter ranges from 0.007 mm to 0.07 mm. The plane was located at the downstream of injector. By counting the frequency of droplets with different diameters to pass through the plane for calculation purpose, the plane was divided into number of sections with each section having a fixed width. Each section was analyzed for the droplet, which was assumed to be spherical. After that, each droplet diameter was calculated from droplet volume obtained by taking into account the volume fraction occupying the computational meshes. Also, the cumulative probability distribution function was used to describe the probability of finding the diameter of the secondary droplets in a sample of splashing drops. The mass of splashing droplets from the wall depends upon the splashing energy of the droplets [45]. In the system of LM, SM, and LSM, while the UWS droplets were injected into the gas phase, the UWS droplets residence time ranged from 0.004 s to 0.121 s, 0.004 s to 0.091 s, and 0.004 s to 0.064 s, respectively. From this analysis, the LSM-installed system showed obvious differences of 47% and 29% decreases of residence time compared to LM and SM, respectively. for uniformity were 0.95 and 0.96 for the velocity and ammonia, respectively. For the purpose of mass transfer and homogenization, when the mixture species passes through the second mixer of LSM, it creates swirl and endorse the mass transfer effect near the pipe wall. The shear and swirl flows with the pipe wall produce molecular diffusion and eddy diffusion under forced convection. Hence, the mixing of ammonia with exhaust gases is comparatively sufficient which results the effective improvement in the uniformity of velocity and ammonia, simultaneously.

Droplet Residence Time
UWS injection in the exhaust gases lowers the temperature of the gas phase due to heat transfer phenomena. Figure 9 represents the UWS droplets residence time for diameter ranges from 0.007 mm to 0.07 mm. The plane was located at the downstream of injector. By counting the frequency of droplets with different diameters to pass through the plane for calculation purpose, the plane was divided into number of sections with each section having a fixed width. Each section was analyzed for the droplet, which was assumed to be spherical. After that, each droplet diameter was calculated from droplet volume obtained by taking into account the volume fraction occupying the computational meshes. Also, the cumulative probability distribution function was used to describe the probability of finding the diameter of the secondary droplets in a sample of splashing drops. The mass of splashing droplets from the wall depends upon the splashing energy of the droplets [45].In the system of LM, SM, and LSM, while the UWS droplets were injected into the gas phase, the UWS droplets residence time ranged from 0.004 s to 0.121 s, 0.004s to 0.091s, and 0.004s to 0.064s,respectively. From this analysis, the LSM-installed system showed obvious differences of 47% and 29% decreases of residence time compared to LM and SM, respectively. The angle of injection of UWS was 30, and the injection was done at the center of pipe; therefore, it is difficult for droplets to move far in radial directions but as an alternative of it, droplets will try to moves toward center. The droplets with bigger diameter will not follow the flow but start to move in The angle of injection of UWS was 30, and the injection was done at the center of pipe; therefore, it is difficult for droplets to move far in radial directions but as an alternative of it, droplets will try to moves toward center. The droplets with bigger diameter will not follow the flow but start to move in radial directions. Small droplets evaporate faster than the bigger diameter and can easily be blown to the center of pipe. However, with the use of LSM, the mixing ability of UWS droplets and exhaust gas enhance sufficiently, which resist the gathering of larger UWS droplets and avert the temperature drop at the center of pipe. Moreover, the combination of two mixers results to increase the turbulence intensity which produces high evaporation rate. Overall, in the system of LSM, the residence time of UWS droplets decreases not only due to the improvement in the distribution of droplets and but also due of high turbulence impact of two mixers.

Urea Conversion
Better conversion or evaporation of UWS results the enough production of ammonia for SCR reaction. With the use of mixers, primary UWS droplets can be distributed and broken up into secondary droplets resulting in quick evaporation of UWS. Conversion of urea from the mixer downstream to the catalyst inlet (Figure 1 blue arrow) for LM, SM, and LSM are shown in Figure 10. It has been observed that there is deficient conversion of urea with the use of single mixer. As the decomposition distance increases, the conversion of urea also increases gradually. For the system of LM and SM, the conversion of urea is 76.1% and 83.2%, respectively. However, for the LSM system, the conversion increases up to 95.4%. Urea conversion is directly related with the mixing of flow distribution and droplet residence time in the SCR system. radial directions. Small droplets evaporate faster than the bigger diameter and can easily be blown to the center of pipe. However, with the use of LSM, the mixing ability of UWS droplets and exhaust gas enhance sufficiently, which resist the gathering of larger UWS droplets and avert the temperature drop at the center of pipe. Moreover, the combination of two mixers results to increase the turbulence intensity which produces high evaporation rate. Overall, in the system of LSM, the residence time of UWS droplets decreases not only due to the improvement in the distribution of droplets and but also due of high turbulence impact of two mixers.

Urea Conversion
Better conversion or evaporation of UWS results the enough production of ammonia for SCR reaction. With the use of mixers, primary UWS droplets can be distributed and broken up into secondary droplets resulting in quick evaporation of UWS. Conversion of urea from the mixer downstream to the catalyst inlet (Figure 1 blue arrow) for LM, SM, and LSM are shown in Figure 10. It has been observed that there is deficient conversion of urea with the use of single mixer. As the decomposition distance increases, the conversion of urea also increases gradually. For the system of LM and SM, the conversion of urea is 76.1% and 83.2%, respectively. However, for the LSM system, the conversion increases up to 95.4%. Urea conversion is directly related with the mixing of flow distribution and droplet residence time in the SCR system. Conversion of urea takes place in two stages, one is water vaporization stage, while the other urea vaporization stage as shown in Figure 11. Initially, the evaporation of pure water occurs from UWS droplets; therefore it becomes unstable in the exhaust gas in the later stage. As a result, the decomposition of urea produces NH3 and HCNO. The droplets with bigger diameter can only undergo water vaporization stage and hence does not produces sufficient ammonia in long distances. With the use of LSM, bigger droplets are break up into small droplets and need short time to evaporate which ultimately results in increased evaporation rate. Conversion of urea takes place in two stages, one is water vaporization stage, while the other urea vaporization stage as shown in Figure 11. Initially, the evaporation of pure water occurs from UWS droplets; therefore it becomes unstable in the exhaust gas in the later stage. As a result, the decomposition of urea produces NH 3 and HCNO. The droplets with bigger diameter can only undergo water vaporization stage and hence does not produces sufficient ammonia in long distances. With the use of LSM, bigger droplets are break up into small droplets and need short time to evaporate which ultimately results in increased evaporation rate.

Relation of Temperature with SCR performance
Two important factors relating to the SCR system temperature include deposit formation and catalyst reaction efficiency. SCR is a means of converting NOx emissions into N2 and H2O by using the catalyst and O2 in the temperature window of 280 to 420 °C [11]. The rate of reaction will be slow and unwanted reactions will occur, if the temperature goes below 280 °C, resulting in poor SCR catalyst performance. Ammonia will start to burn without reacting with the NOx emissions if the temperature goes above 420°C. It is therefore recommended to control the reaction temperature of SCR system [46].
The wall temperature is the prominent source for finding the deposit formation. Urea starts to decompose rapidly, associated with the secondary reactions, if the temperature is more than 163 °C. If the temperature is in between the range of 133 to 163 °C, pyrolysis of urea occurs slowly. Urea crystals are formed certainly as temperature is reduced below 133 °C. Once the temperature goes down, urea crystals are produced and the exhaust pipe will be blocked, which is responsible for decreasing the mixing performance and increasing the back pressure [18]. Concurrently, the overall catalytic reactor efficiency is decreased considerably. Once the urea injected into the system, a number of small droplets produce collisions with the pipe wall and mixer. The increase in collision strength of droplets with pipe wall decreased the temperature. Under the Leidensforst temperature, droplets of urea are separated into four boiling phenomenon as maintained by three values of temperature: 140, 190, and 300 °C. Heat transfer to the liquid film from pipe wall increased at ~180 °C [9].If the wall temperature is lower than the boiling point temperature (Tb) of the UWS droplet, the droplets will stick on the wall and, if the temperature is higher than Tb, the droplets will rebound after striking with the wall [45]. Therefore, it is important to obtain maximum liquid film temperature for improved performance of urea decomposition and decrease the deposit formation.

Temperature Distribution along the Mixer Downstream
The temperature distribution along the radial direction at the mixer downstream (Figure 1 red arrow) as represented in Figure 12. The pipe diameter is 10 cm; therefore, the initial position (0 cm) Figure 11. Urea water solution (UWS) droplet at two vaporization stages.

Relation of Temperature with SCR Performance
Two important factors relating to the SCR system temperature include deposit formation and catalyst reaction efficiency. SCR is a means of converting NO x emissions into N 2 and H 2 O by using the catalyst and O 2 in the temperature window of 280 to 420 • C [11]. The rate of reaction will be slow and unwanted reactions will occur, if the temperature goes below 280 • C, resulting in poor SCR catalyst performance. Ammonia will start to burn without reacting with the NO x emissions if the temperature goes above 420 • C. It is therefore recommended to control the reaction temperature of SCR system [46].
The wall temperature is the prominent source for finding the deposit formation. Urea starts to decompose rapidly, associated with the secondary reactions, if the temperature is more than 163 • C. If the temperature is in between the range of 133 to 163 • C, pyrolysis of urea occurs slowly. Urea crystals are formed certainly as temperature is reduced below 133 • C. Once the temperature goes down, urea crystals are produced and the exhaust pipe will be blocked, which is responsible for decreasing the mixing performance and increasing the back pressure [18]. Concurrently, the overall catalytic reactor efficiency is decreased considerably. Once the urea injected into the system, a number of small droplets produce collisions with the pipe wall and mixer. The increase in collision strength of droplets with pipe wall decreased the temperature. Under the Leidensforst temperature, droplets of urea are separated into four boiling phenomenon as maintained by three values of temperature: 140, 190, and 300 • C. Heat transfer to the liquid film from pipe wall increased at~180 • C [9]. If the wall temperature is lower than the boiling point temperature (T b ) of the UWS droplet, the droplets will stick on the wall and, if the temperature is higher than T b , the droplets will rebound after striking with the wall [45]. Therefore, it is important to obtain maximum liquid film temperature for improved performance of urea decomposition and decrease the deposit formation.

Temperature Distribution along the Mixer Downstream
The temperature distribution along the radial direction at the mixer downstream (Figure 1 red arrow) as represented in Figure 12. The pipe diameter is 10 cm; therefore, the initial position (0 cm) stands for the upper edge of the pipe and 10 cm shows the position of lower edge. A plane was created in CFD at a certain distance and then divided into number of points. The value of temperature on each single point was found. Next, the standard deviation and average value of all the points were calculated. The minimum temperature of LM at the top edge is 274 • C, which is 32 • C less than SM and LSM.As it has been proved that after the urea injection process a large number of droplets collide with the mixer and the pipe wall: the wall temperature decreases with the increase in collision intensity [29]. In all types of mixers arrangement, the temperature difference at the center of the pipe is not obvious. However, the temperatures of LM and SM at the bottom edge of pipe were 276 • C and 281 • C, respectively. The temperatures at top edge and bottom edge of pipe are 304 • C and 301 • C respectively. Generally, the combination of two mixers (LSM) have good temperature distribution for both upper edge as well as lower edge, which is beneficial for catalyst reaction performance and also very suitable to prevent the deposit formation.
The axial wall temperature distribution from the mixer downstream to the catalyst inlet (blue arrow in Figure 1) as represented in Figure 13. Wall temperature distributions of LM and SM at the catalyst inlet were289 • C and 294 • C, respectively, but the combination of two mixers (LSM) was 300 • C due to better heat transfer effect. A continuous decrease in the temperature was observed from 0.3 to 0.5 m distance; with minimal temperature of 270 • C in SM. Low temperature region produces direct effect on the wall of pipe and reaction performance, without creates disturbance for the mixer.
Generally, temperature plays an important role in the reaction performance and deposit formation. With the use of LSM, the temperature remains above 280 • C in axial and radial directions, which is very helpful for preventing the unwanted reactions and deposit formations. stands for the upper edge of the pipe and 10 cm shows the position of lower edge. A plane was created in CFD at a certain distance and then divided into number of points. The value of temperature on each single point was found. Next, the standard deviation and average value of all the points were calculated. The minimum temperature of LM at the top edge is 274 °C, which is 32 °C less than SM and LSM.As it has been proved that after the urea injection process a large number of droplets collide with the mixer and the pipe wall: the wall temperature decreases with the increase in collision intensity [29]. In all types of mixers arrangement, the temperature difference at the center of the pipe is not obvious. However, the temperatures of LM and SM at the bottom edge of pipe were 276 °C and 281 °C, respectively. The temperatures at top edge and bottom edge of pipe are 304 °C and 301 °C respectively. Generally, the combination of two mixers (LSM) have good temperature distribution for both upper edge as well as lower edge, which is beneficial for catalyst reaction performance and also very suitable to prevent the deposit formation. The axial wall temperature distribution from the mixer downstream to the catalyst inlet (blue arrow in Figure 1) as represented in Figure 13. Wall temperature distributions of LM and SM at the catalyst inlet were289 °C and 294 °C, respectively, but the combination of two mixers (LSM) was 300 °C due to better heat transfer effect. A continuous decrease in the temperature was observed from 0.3 to 0.5 m distance; with minimal temperature of 270 °C in SM. Low temperature region produces direct effect on the wall of pipe and reaction performance, without creates disturbance for the mixer.
Generally, temperature plays an important role in the reaction performance and deposit formation. With the use of LSM, the temperature remains above 280 °C in axial and radial directions, which is very helpful for preventing the unwanted reactions and deposit formations.

Working Performance of LSM-basedSCR System
The standard ISO 8178 is the international marine Diesel engine test cycle used to measure the exhaust emissions from ships. By following the control requirements of exhaust emission as per IMO Tier III, standard marine Diesel engine test cycle has been divided into two parts one is ISO 8178 D2 test cycle for marine Diesel engine operated with constant speed and other is ISO 8178 E3 test cycle with propelling character for marine Diesel engine. In this study marine Diesel engine with propelling character was studied. Hence, the ISO 8178 E3 test cycle was used to measure the exhaust emissions. The weightage average value of NOx exhaust emission was calculated by ISO 8178 marine Diesel engine test cycle E3 as shown in Table 4 [30].
The main purpose of SCR system is to decrease the NOx emissions and to prevent the ammonia leakages responsible for air pollution. Figure 14 represents the NOx removal efficiency under different loading conditions. It was observed that, there is a small difference in between theoretical value (standard) and calculated value (simulation). At low load, the NOx removal efficiency is low, Figure 13. Axial distribution of temperature from mixer downstream to catalyst inlet in pipe wall.

Working Performance of LSM-Based SCR System
The standard ISO 8178 is the international marine Diesel engine test cycle used to measure the exhaust emissions from ships. By following the control requirements of exhaust emission as per IMO Tier III, standard marine Diesel engine test cycle has been divided into two parts one is ISO 8178 D2 test cycle for marine Diesel engine operated with constant speed and other is ISO 8178 E3 test cycle with propelling character for marine Diesel engine. In this study marine Diesel engine with propelling character was studied. Hence, the ISO 8178 E3 test cycle was used to measure the exhaust emissions. The weightage average value of NO x exhaust emission was calculated by ISO 8178 marine Diesel engine test cycle E3 as shown in Table 4 [30]. The overall weighted NO x exhaust emission level in g/kWh can be calculated using Equation (16) [47].
The main purpose of SCR system is to decrease the NO x emissions and to prevent the ammonia leakages responsible for air pollution. Figure 14 represents the NO x removal efficiency under different loading conditions. It was observed that, there is a small difference in between theoretical value (standard) and calculated value (simulation). At low load, the NO x removal efficiency is low, because under lower loading conditions, the exhaust temperature of Diesel engine is comparatively low, which directly affects the catalyst performance; as a result, the catalyst efficiency decreases. Figure 15 shows the NH 3 escaping rate at different engine loading conditions. The rate of NH 3 escaping decreases with the increase of engine load. NH 3 escaping rate is relatively high at low load due to lower temperature of the exhaust gas, resulting in incomplete catalyst reaction. The reaction of NH 3 oxidization is neglected for design of SCR model parameters in this study; therefore, more NH 3 slipping occurs at low loading conditions. However, at higher loading conditions (75% and 100%), the NH 3 escaping rate is less than 10 ppm, which meets the design requirement of the SCR system. low, which directly affects the catalyst performance; as a result, the catalyst efficiency decreases. Figure 15 shows the NH3 escaping rate at different engine loading conditions. The rate of NH3escaping decreases with the increase of engine load. NH3 escaping rate is relatively high at low load due to lower temperature of the exhaust gas, resulting in incomplete catalyst reaction. The reaction of NH3 oxidization is neglected for design of SCR model parameters in this study; therefore, more NH3 slipping occurs at low loading conditions. However, at higher loading conditions (75% and 100%), the NH3 escaping rate is less than 10 ppm, which meets the design requirement of the SCR system.
Verification of simulated results was confirmed by using ISO 8178 standard marine Diesel engine test cycle E3.According to emission regulations of IMO Tier III, the value of NOx emission should be less than 3.4 g/kWh under all loading conditions. Figure 16 shows a NOx emission value of the LSM-based SCR system at different loading conditions. It was observed that NOx emission decreases as the load increases. The average weighted value of NOx emission was 2.44 g/kWh for four different loads at the downstream of SCR catalyst. Hence, the system based on using LSM strongly meets the standard of IMO Tier III NOx emission regulations effectively.   low, which directly affects the catalyst performance; as a result, the catalyst efficiency decreases. Figure 15 shows the NH3 escaping rate at different engine loading conditions. The rate of NH3escaping decreases with the increase of engine load. NH3 escaping rate is relatively high at low load due to lower temperature of the exhaust gas, resulting in incomplete catalyst reaction. The reaction of NH3 oxidization is neglected for design of SCR model parameters in this study; therefore, more NH3 slipping occurs at low loading conditions. However, at higher loading conditions (75% and 100%), the NH3 escaping rate is less than 10 ppm, which meets the design requirement of the SCR system.
Verification of simulated results was confirmed by using ISO 8178 standard marine Diesel engine test cycle E3.According to emission regulations of IMO Tier III, the value of NOx emission should be less than 3.4 g/kWh under all loading conditions. Figure 16 shows a NOx emission value of the LSM-based SCR system at different loading conditions. It was observed that NOx emission decreases as the load increases. The average weighted value of NOx emission was 2.44 g/kWh for four different loads at the downstream of SCR catalyst. Hence, the system based on using LSM strongly meets the standard of IMO Tier III NOx emission regulations effectively.   Verification of simulated results was confirmed by using ISO 8178 standard marine Diesel engine test cycle E3. According to emission regulations of IMO Tier III, the value of NO x emission should be less than 3.4 g/kWh under all loading conditions. Figure 16 shows a NO x emission value of the LSM-based SCR system at different loading conditions. It was observed that NO x emission decreases as the load increases. The average weighted value of NO x emission was 2.44 g/kWh for four different loads at the downstream of SCR catalyst. Hence, the system based on using LSM strongly meets the standard of IMO Tier III NO x emission regulations effectively.

Conclusions
In this work, numerical methods were used to analyze the mixing performance, effects of mixers on the evaporation rate of urea, residence time of urea droplet in the pipe, and temperature distribution for the catalyst reaction and deposit formation. The prime results are shown below.

•
For the in-line type mixer (LM), the uniformity index of velocity was good (0.93) but the uniformity of ammonia was poor (0.87). In contrary to LM, the swirl type mixer (SM) has good ammonia uniformity (0.94) but poor uniformity index of velocity (0.86).However, better values were observed by using combination of two mixers (LSM). The uniformity index of velocity and ammonia uniformity achieved the values of 0.95 and 0.96, respectively, for LSM-based SCR system.

•
The residence time UWS was studied. The results show that the residence time of urea droplets for LSM-based SCR system was 0.064 s, which represents47% and 29% decreases compared to LM and SM, respectively. Furthermore, the conversion of urea into ammonia is highly related with the residence time of urea droplets in the pipe. Hence, urea conversion achieves the value of 95.4% by using LSM, which is 19.3% and 12.2% higher than the value of LM and SM, respectively.

•
It was also observed that the combination of two mixers (LSM) have good temperature distribution than the LM and SM for radial and axial directions, the temperature at catalyst inlet in axial direction was 300 °C for LSM-based SCR system which is suitable for the catalyst reaction performance and prevents the deposit formation. • Finally, the simulated results of the model parameters were compared and verified by using ISO 8178 standard marine Diesel engine test cycle E3. The average weighted value of NOx emission was calculated as 2.44 g/kWh for four different loads. Hence, it is concluded that the system based on using LSM strongly meets the standard of IMO Tier III NOx emission regulations effectively.

Conflicts of Interest:
The authors declare no conflicts of interest.

Conclusions
In this work, numerical methods were used to analyze the mixing performance, effects of mixers on the evaporation rate of urea, residence time of urea droplet in the pipe, and temperature distribution for the catalyst reaction and deposit formation. The prime results are shown below.

•
For the in-line type mixer (LM), the uniformity index of velocity was good (0.93) but the uniformity of ammonia was poor (0.87). In contrary to LM, the swirl type mixer (SM) has good ammonia uniformity (0.94) but poor uniformity index of velocity (0.86). However, better values were observed by using combination of two mixers (LSM). The uniformity index of velocity and ammonia uniformity achieved the values of 0.95 and 0.96, respectively, for LSM-based SCR system.

•
The residence time UWS was studied. The results show that the residence time of urea droplets for LSM-based SCR system was 0.064 s, which represents47% and 29% decreases compared to LM and SM, respectively. Furthermore, the conversion of urea into ammonia is highly related with the residence time of urea droplets in the pipe. Hence, urea conversion achieves the value of 95.4% by using LSM, which is 19.3% and 12.2% higher than the value of LM and SM, respectively. • It was also observed that the combination of two mixers (LSM) have good temperature distribution than the LM and SM for radial and axial directions, the temperature at catalyst inlet in axial direction was 300 • C for LSM-based SCR system which is suitable for the catalyst reaction performance and prevents the deposit formation. • Finally, the simulated results of the model parameters were compared and verified by using ISO 8178 standard marine Diesel engine test cycle E3. The average weighted value of NO x emission was calculated as 2.44 g/kWh for four different loads. Hence, it is concluded that the system based on using LSM strongly meets the standard of IMO Tier III NO x emission regulations effectively.