RNG Model Research Articles

Comparing and analyzing the accuracy and stability of calculating SF6 tracer gas mixing results in mine tunnels using RNG model and DES model through experiments

The tracer gas wind measurement method must measure the average concentration based on the known mixing distance of the tracer gas in order to accurately calculate the air volume. we established a T-shaped cross-ventilation duct with a square cross-sectional side of 0.3 m to simulate mine tunnels experimentally, and created the above duct model in Fluent software. Using the renormalization group (RNG k-ε) and detached eddy simulation (DES) models in Fluent to calculate the concentration and mixing distance of sulfur hexafluoride (SF6) tracer gas in air collection duct, and conducted SF6 blending experiments. In the first 7m stage of air collection duct, the SF6 mixing rate in the experiment was the fastest, while the RNG k-ε model was the slowest. However, after 7 meters, the mixing rate in the experiment rapidly slowed down, while there was no significant change in the RNG k-ε model. The rate of change of the DES model was between the above two results, but the concentration change was the most unstable. And the mixed distance results of the three were highly consistent. The results of this study provide a reference for selecting appropriate numerical simulation models for future calibration work of mine ventilation systems.

An Analysis of Air Conditioning System Operation Patterns on Droplet Particle Distribution in a Classroom-A CFD Approach

In the classroom context, the transmission of pathogens among students is a significant concern. Therefore, it is important to determine appropriate airflow patterns, the placement of supply and exhaust ventilation, and the optimization of classroom design to reduce the risk of pathogen transmission. This research aims to determine the performance of two air conditioning (AC) operating patterns—low speed and high speed—in six scenarios involving window and door configurations and identify the most effective strategies for minimizing virus exposure to occupants in classrooms. The method used in this research is numerical simulation with the 3D unsteady k-ε RNG model to simulate air flow and the Eulerian-Lagrange approach to capture the movement of SARS-CoV-2 aerosol droplets. The results of this research show that of the six scenarios determined by the researchers, the low-speed AC operating pattern with an incoming air speed of 3.5 m/s occurs in scenario 5, that is, all windows open and doors closed. This is based on the lowest number of students exposed to the virus, which is 22.22%. Meanwhile, the high-speed AC operating pattern with an incoming air speed of 6 m/s occurs in scenario 2, that is, all windows closed and doors open. This is based on the lowest number of students exposed to the virus, which is 22.22%, so it can be concluded that increasing the air flow speed originating from the AC will speed up the droplets to leave the room through the outlet. Meanwhile, increasing the outlet capacity will shorten the particle path, thereby shortening the time when the droplets are in the classroom.

Performance enhancement of solar air heater with two-sided curvilinear transverse rib: Experimental and numerical investigation

Solar energy as prime energy source draws attention of world due to its availability and eco-friendliness. Low performance of solar air heater create curiosity for researchers to increase its performance. To enhance performance of solar air heater active and passive techniques used by researcher. Surface modification is one of the most prominent passive techniques to improve the performance. Present experimental and numerical investigation explore thermal characteristics of two-sided curvilinear rib roughened solar air heater. An indoor experiment was performed on two side curvilinear rib roughened absorber plate having constant heat flux of 1000 W/m2. Parallelly 2-dimensional numerical simulation was also performed to explore flow behavior and to complete investigation within optimum cost. ANSYS Fluent 2021.R was used for performing this simulation with K-ε RNG model with enhanced wall treatment. Velocity inlet and pressure outlet was selected as boundary condition. Main roughness and flow parameters were relative roughness height ( e/ D h = 0.21–0.042), relative roughness pitch ( p/ e = 7.14–35.71), and Reynolds number ( Re = 3800–18,000). Maximum value of Nusselt number improvement ratio and friction factor improvement ratio was 3.17 and 3.26 at roughness and flow parameter of ( e/ D h = 0.042, p/ e = 15, Re = 15,000) and ( e/ D h = 0.042, p/ e = 7.14, Re = 3800) respectively. Thermohydraulic performance parameter achieved its maximum value of 2.77 at e/ D h = 0.042, p/ e = 17.85, and Re = 18,000.

Optimizing of heat transfer and flow characteristics within a roughened solar air heater duct with compound turbulators

AbstractThermal systems for solar air heating have been widely used in both industrial and residential contexts, and are essential for converting and recovering solar energy. Thermal performance in solar air heaters (SAHs) can be improved through the repetitive application of artificial roughness to the surfaces. This research work includes a numerical evaluation of SAH performance with artificial rough surfaces made up of combined transverse trapezoidal ribs and chamfered grooves. The ANSYS Fluent software version 2023 R1 was used to simulate SAH with varying relative roughness pitch (), relative roughness heights (), and Reynolds number (). The RNG model was chosen to forecast an enhancement in Nusselt number (), friction factor (), and thermohydraulic performance factor (TPF) for the proposed roughness. Out of multiple roughness parameters analyzed, it was determined that the compound turbulator with and , were the most effective. The TPF for this scenario was determined to be 2.12 at . Finally, a numerical based empirical correlations for and in terms of Re, , and were developed.

CFD Turbulence Models Assessment for the Cavitation Phenomenon in a Rectangular Profile Venturi Tube

This study investigates cavitation in a rectangular-profile Venturi tube using numerical simulations and four turbulence models. The unsteady Reynolds-averaged Navier–Stokes technique is employed to simulate vapor cloud formation and compared against experimental data. κ-ε realizable, κ-ε RNG, κ-ω SST, and κ-ω GEKO models are evaluated. The simulation results are analyzed for pressure, turbulence, and vapor cloud formation. Discrepancies in cavitation cloud formation among turbulence models are attributed to turbulence and vapor cloud interactions. RNG and SST models exhibit closer alignment with the experimental data, with RNG showing a superior performance. Key findings include significant vapor cloud shape differences across turbulence models. The RNG model best predicts velocity at the throat exit with an error of 4.145%. Static pressure predictions include an error of 4.47%. The vapor cloud length predictions show variation among models, with the RNG model having a 0.386% error for the minimum length and 4.9845% for the maximum length. The SST model exhibits 4.907% and 13.33% errors for minimum and maximum lengths, respectively. Analysis of the cavitation number reveals agreement with the experimental data and sensitivity to cavitation onset. Different turbulence models yield diverse cloud shapes and detachment points. Weber number contours illustrate the variation in the cavitation cloud behavior under different turbulence models.

Thermal analysis of a turbulent wall jet over an adiabatic wavy surface

The impact of wavy wall amplitude on the turbulent characteristics of a wall jet has been investigated experimentally, where the wavy wall is given an adiabatic condition. In the experimental study, the mean and turbulent velocity are measured using a constant temperature anemometer manufactured by Dantec Dynamics. For temperature measurement over the wall and within the flow domain, the FLIR camera and k-type thermocouple are used, respectively. For the current study, uniformly heated at △T0 = (T0 − T∞) = 30 ± 1.2 0C with a Reynolds number of 15,000 is utilized to flow over a sinusoidal wavy wall (y = A*sin (ωx)), where “A = Amplitude ⁄ a” is the amplitude normalised by the nozzle height “a". The nozzle height is 20 mm and to maintain the Reynolds number 15,000 at the inlet, velocity is set at 10.95 ± 0.36 m ⁄ s. The amplitude of the wavy wall is varied between 0.2 and 0.6 and the results are compared to those of the plane wall jet. The current problem is also numerically solved using the k − ε RNG and k − ω SST models. Based on the results, it is concluded that the k − ε RNG model predicts the wavy wall results quite well. It is also discovered that the thermal self-similar profile deviates from the Gaussian curve at the crest of the wavy wall. The influence of the wavy surface on the thermal characteristics of the turbulent wall jet is more pronounced in the near flow field. It is observed that the bottom wall temperature, fluctuating temperature and turbulent heat flux decrease with the increase in amplitude. This indicates that for the higher amplitude, the turbulence increases, which increases the intermixing of jet fluid with the colder ambient fluid. The present results can be used as a benchmark for validating numerical models working in this area.

SRANS Simulation on a Helical-Segmented Finned Tube Bank

A numerical simulation on a six-row helical-segmented finned tube in a staggered layout is performed with the Reynolds Averaged Navier-Stokes (RANS) approach. The turbulence effect is represented with the k-ε RNG model, and a fixed temperature in each tube row considers the inside fluid temperature. The tube bank is represented by a small finned tube bank and a single interacted module. The small finned tube bank considers symmetry conditions on the sides of the computational domain. The single interacted module considers periodic boundary conditions for velocity, pressure, and temperature. Both simulations are compared to develop the basis of complex simulations such as repetitive finned tube modules. Average values on parallel planes to the streamwise direction for velocities, pressures, and temperatures are compared in both simulations. The comparative analyses show deviations lower than 11% for the velocity field, 12% for the pressure field, and 10% single for the temperature field. Therefore, results show an appropriate flow representation with periodic boundary conditions.

Impact toughness analysis of offshore wind power structures under the influence of long-period waves

Abstract In this paper, the performance of offshore wind turbine structures under long-period wave impacts is investigated, and a numerical model of long-period waves is developed to simulate the wave motion and fluid seepage in the pore medium by using the VARANS equation with the OlaFlow solver, and various turbulence models such as the model, RNG model, and the VOF method is applied to capture free surfaces, which can accurately simulate wave generation, propagation, reflection, breaking, and fluid seepage in the pore medium. These methods can accurately simulate the wave generation, propagation, reflection, breaking, and fluid seepage in the pore medium, and the accuracy of the numerical simulation is verified by comparing the results with those of the physical experiment. The results show that the wind farm exhibits good impact toughness under the influence of long period waves, and its overturning stability and slip stability are better than the safety coefficient required by the specification.

Systematic assessment of computational models for the prediction of flow fields in bubble columns

Systematic changes in interfacial forces, turbulence models for the continuous and dispersed phases, lift force coefficient, and bubble size distribution were used to evaluate, in terms of complexity, the performance of nine computational models in the simulation of the time-averaged axial liquid velocity and holdup profiles in a bubble column. The oscillation period of the gas plume was also included in this evaluation. The mesh used for the simulations was chosen after a convergence analysis carried out with four meshes. The finest mesh displayed the highest divergence from the experimental data. Such behavior indicates that the proper mesh for a simulation of a multiphase flow must consider other factors in addition to element size and the quality indicators. The standard, RNG, and realizable κ-ε models were used for the liquid phase, while an algebraic and the standard κ-ε were the models for the dispersed phase. The isotropic RNG model for the liquid phase achieved accurate predictions in the anisotropic region near the gas distributor. The algebraic approach for the gas phase produced data similar to those calculated with the standard model. The simultaneous presence of the lift and virtual mass forces led to reliable predictions, and a lift coefficient of 0.2 displayed better performance than 0.14. This coefficient produced an excessive dampening of the lift force and caused the simulations to lose accuracy. The population balances did not improve the predictions.

Enhancement of heat transfer by flowing turbulent jet on a linearly decaying (LD) wavy wall

The influence of a partial linearly decaying (LD) wavy wall on different fluid flow and thermal parameters has been investigated for the turbulent jet. For the present study, a segment of wall from the leading edge is made with a constant amplitude 0.8a (a = nozzle height), followed by a segment of LD wavy wall. In the region of LD wavy wall, the amp (amplitude) of wavy wall is linearly reduced from 0.8a amplitude to 0a in the downstream direction. The LD wavy wall segment varies from 100% to 0%. For 100% LD wavy wall, the whole wall is made of linearly decaying amplitude, whereas for, 0% LD wavy wall, the whole wall is made wavy with a constant amplitude. The Reynolds number at the inlet is same for all the cases, i.e. 15,000 to maintain the fully turbulent condition at the exit of nozzle. The problem is solved numerically using k−ε RNG model. The results obtained from this study show that for 100% LD wavy wall the velocity decay is minimum and as the % of LD wavy wall decreases, the velocity decay increases. Whereas the thermal decay of jet is maximum in the case of 100% LD wavy wall. Similarly, Ymax, YTmax, and the velocity and thermal jet spread are found to be minimum for the 100% LD wavy wall and all these parameters increase as the % of LD wavy wall decreases. The maximum increment in the heat transfer rate is 27.9% for 100% LD wavy wall in regards to the baseline case (plane wall). However, the maximum increment in thermal hydraulic performance (THP) is noticed for the wavy wall with 50% LD wavy wall, which is equal to 6.7%. In a previous study (ArchanaKumari, 2021 [1]), the increment of 19.0% and 5.3% is observed for the heat transfer and THP respectively. In the present study by introducing LD wavy wall, the increment in the heat transfer and THP are quite significant.

Experimental and numerical investigation of the thermal and dynamic behavior of a heated vortex multijet system impacting a flat plate

This study concerns the experimental and numerical study of the thermal and dynamic behavior of a configuration of a system of vortex jets impacting a flat plate, The objective of this study is to study the behavior of the thermal and dynamic field of vortex blowing of hot air from a multi-jet system impacting a flat plate. The experimental test bench comprising a support of three diffusers of diameter D, impacting the perpendicular plate. A uniform inlet temperature (T, T, T) is imposed such that the impact height H = 4D. The vortex is obtained by a vortex generator made up of 12 fins arranged at 60° from the vertical, placed just at the outlet of the diffuser. A thermo-anemometer device, to measure the blowing temperature at the point in question. The system was numerically simulated by the fluent code using a k-ε RNG turbulence model. It should be noted that the multi-jet system first appears as a free jet: going from the injection orifice to the impact zone, the axial velocity weakens, the jet undergoes considerable deflection, this is the deflection zone the velocities become mainly radial and the thickness of the boundary layer increases radially: this is the parietal flow zone, the structure of the velocity field has two zones of intense deflection with a wall jet on both sides other, favoring a good development of the resulting jet. The results show that this configuration (T, T, T) gave a better optimized distribution of temperature and velocity on the surface of the plate. This homogenization of the temperatures results from a better thermal transfer of the plate.The k-ε RNG model gave acceptable results, which coincide with those of the experimental results.

CFD SIMULATION AND VALIDATION FOR MIXING VENTILATION SCALED-DOWN EMPTY AIRCRAFT CABIN USING OPENFOAM

An investigation into the spread of the COVID-19 virus within a confined space including an aircraft cabin is essential in order to find out the mechanism. However, this is time-consuming and limited in scope, so a computational fluid dynamics (CFD) simulation is used instead. Therefore, a prior study and an appropriate choice of turbulence model are required before the simulation. The main objective of this study is to validate and evaluate the results predicted by the Open Source Field Operation and Manipulation (OpenFOAM) software through comparison with the experimental data from the literature which was conducted using particle image velocimetry (PIV) measurement. Three different Reynolds-averaged Navier-Stokes turbulence models were selected; Re-normalisation Group k - ɛ (RNG), Realizable k - ɛ (RLZ) and Low Reynold Number (LRN) to simulate a mixing ventilation system of a scaled-down model of empty aircraft cabin. In the RNG and LRN model cases, a fairly large circulation flows were observed on the right and left sides of the model representing the passenger area. The results were also evaluated quantitatively using the factor of two of observations (FAC2) for the velocity components and turbulent kinetic energy (TKE) with root mean square error (RMSE) for the former and normalised mean square errors (NMSE) for the latter. The simulation results showed that RNG and LRN are capable of predicting the flow field well. However, for TKE prediction LRN performed better than RNG which concluded that LRN is the suitable turbulence model in simulating flow fields in investigated case.

Double-skin façade simulation with computational fluid dynamics: A review of simulation trends, validation methods and research gaps

Dynamic simulation of a double-skin façade (DSF) with computational fluid dynamics (CFD) can be challenging due to the lack of validated models and benchmarking datasets. Furthermore, there is a lack of consensus in the scientific community on what constitutes a successfully validated DSF model. The present review study identifies simulation trends and research gaps for DSFs simulated with CFD. Additionally, this article presents a series of CFD simulations in which key aspects of the DSF modelling are varied: 2D or 3D modelling approaches, turbulence viscosity models (TVMs), radiation models, and wall function. These simulation results are compared to the empirical data (both temperature and velocity fields) of a benchmark test with laboratory-controlled boundary conditions. This analysis shows that using the k-ε RNG model with enhanced wall treatment and surface-to-surface (S2S) radiation model yields the best results for the 2D case of natural convection flow. Moreover, it is shown that accounting for the velocity field in the validation process is essential to ensure the suitability of a model. Finally, the authors advocate for the use of selected dimensionless numbers to improve the comparability of the different DSF scientific studies. This would also help to identify relevant experimental datasets for validation and suitable CFD simulation settings for specific DSF cases.

Thermal Hydraulic Performance of a Turbulent Wall Jet Flowing Over a Wavy Wall

Abstract The thermal and fluid flow characteristics of a two-dimensional turbulent wall jet have been studied numerically for the partial wavy wall. The partial wall is created by giving a segment of wall a wavy pattern from the leading edge followed by the plane wall; the wavy wall segment varies from 10% to 100% of the wall. The amplitude of wavy surface has also been varied from 0.2a to 0.8a with the interval of 0.2a, where “a” is the nozzle height. The results of the present problem have been compared with the results of a fully plane wall jet. The Reynolds number at the nozzle exit is constant, i.e., 15000 for all the cases to achieve fully turbulent jet. To solve this problem, low Reynolds number RNG model has been used. The results obtained from the present study show that the heat transfer rate remains almost the same for 10% to 100% wavy wall for 0.2a amplitude. In the case of amplitude 0.8a, the heat transfer rate is maximum for 30% wavy wall case; the heat transfer rate reduces further for higher wavy wall %. There is a 26.27% increment in heat transfer for the 30% wavy wall with 0.8a amplitude relative to the fully plane wall jet. The maximum increment in the thermal hydraulic performance (THP) of 5.3% is achieved for 70% wavy wall portion for 0.8a amplitude and it remains the same for further increase in the % of wavy wall. It must be noted here that earlier a 19.08% increase in average heat transfer is achieved for the wavy wall with an amplitude of 0.7 in the previous study. However, in this paper, an increase of 26.27% is noticed with the new design of wavy wall. This study will help in designing a highly efficient cooling or heating system in the engineering application by modifying just only the leading surface of the wall.

Numerical study on pollutant emissions characteristics and chemical and physical exergy analysis in Mild combustion

This analytical study appraises the emission characteristics of different pollutants such as CO and NOX as well as conducting chemical and physical exergy analyses on conventional and Mild combustion regimes. The RANS turbulence approach and the k-ε RNG model were utilized in this study to simulate turbulence, whereas the eddy dissipation concept (EDC) model and the GRI–Mech 2.11 chemical mechanism were adopted to treat the turbulence-chemistry interaction. This study evaluated local and global pollutant emissions and exergy efficiency along the reacting jet and explored the effects of temperature and chemical species on the physical and chemical exergies in MILD combustion in comparison to conventional combustion. It is found that the overall exergy efficiency of MILD combustion rises by 1, 3, and 4% for oxygen mass fractions of 3, 6, and 9% in hot oxidizer coflow, with the NOx emissions are 81, 72, and 55% lower than conventional combustion, respectively. Chemical exergy is greater than physical exergy near the fuel nozzle due to high CH4 and H2 concentrations. Chemical exergy reduces along the burner due to chemical reactions, while physical exergy increases due to the increased temperature. However, owing to the trade-off between these two exergies along the flame jet, the total exergy approaches a constant value at large distances from the burner. The maximum total exergy occurs at a distance of 0.12 m from the fuel nozzle, with chemical exergy accounting for 83% of the total exergy. In the outlet (0.5 m from the fuel nozzle), chemical exergy accounts only for 30% of the total exergy. It can also be concluded that the exergy is stabilized at a distance of 0.3 m from the nozzle (60% of the domain length). Despite a higher local CO emission in conventional combustion than in MILD combustion, the CO emission substantially decreases at a smaller distance from the nozzle in conventional combustion due to the high O2 concentration.

The effect of airflow guiding components on effective ventilation rates in single-sided ventilation applications

Wind-driven single-sided ventilation (SSV) is present in many existing buildings across Europe, and with new near-zero energy building (NZEB) regulations for the refurbishment of the existing building stock, its attractiveness as a noninvasive, low-energy solution is set to continue. As a strategy, however, in addition to its air change rate capacity, the distribution of fresh air is an important evaluation criterion for its performance. Airflow guiding components located in the external opening that enhance the effectiveness of the wind-driven flow in ventilating the occupied zone could improve the quality of indoor environments. To our knowledge, the literature is sparse on the practical implications for ventilation when adopting guiding components such as louvres, an increasingly popular approach. In the present study, the performance of wind-dominant SSV was simulated using RNG and RSM CFD models, with and without louvres at three building orientations, for example, windward, parallel and leeward. The purpose of this study was to investigate whether louvres installed in the opening would improve both the effective ventilation rate and the penetration depth of the flow into the indoor space. The performance of SSV was evaluated using the age of air and interpreting the secondary air circulation inside the room affected by louvres. As a result of these investigations, a newly configured airflow guiding component was designed and compared to the other cases. Results show louvres can play a crucial role in controlling the secondary air circulation inside the room, and they could either improve or worsen the performance of SSV in terms of air-exchange efficiency. It was shown that in most cases, if louvres were the cause of incremental changes in turbulent intensity within the indoor space, then they are effective as an air-exchange efficiency improvement strategy.

Hydrodynamic Analysis Inside a Circulating Fluidized Bed Boiler Based on Time Change Using CFD

The power generation industry has progressively improved the design of power plants to tackle global climate change. Coal-fired power plants will be cleaner and more efficient with the CFB (Circulating Fluidized Bed) technology in the boiler. CFB boilers have the advantage where the resulting carbon dioxide, SOx, and NOx emissions are significantly reduced. This makes research on CFB boilers widely carried out, especially those using CFD (Computational Fluid Dynamic). This study aims to determine the most appropriate turbulence model using CFD simulations on CFB boilers to analyze the hydrodynamic elements in the boiler against time changes, including the distribution of solid volume fractions, pressure, velocity, and wall shear stress. Three turbulence models were tested as part of the simulation: standard , RNG , and RSM, which were then compared to earlier research. The RNG model produces the closest result to the experimental data, with an error value of 6.24%. Solid volume fraction area is increasingly widespread with values ranging from 0.05 to 0.1 from the furnace’s bottom to a height of 15 m. The high-pressure area further expands at the bottom of the furnace to a height of 15 m, with values ​​ranging from 3 kPa to 10 kPa. The outlet line velocity has increased from 43 m/s to 50 m/s. The wall shear stress increases from 0.75 Pa to 1.1 Pa at the outlet wall.

Analysis of wave added resistance and motion response of fishing boat sailing against waves

In order to accurately predict the on-wave resistance and responses to hull motions of ships in actual sea conditions, the k-ε method of the RNG model is adopted on the basis of the unsteady RANS method. The twoformula turbulence model deals with the viscous flow, the VOF method captures the free surface, the velocity boundary method makes waves, the artificial damping method is used to eliminate waves, and the nested grid technology is used to deal with the motion response of ships on waves. Combined with the 6-DOF motion formula, a three-dimensional numerical wave cell for regular waves is established. For one example, taking a KCS Container ship and fishing boat sailing at a mid-high-speed, the increase of wave resistance and motion response at different wavelengths are analyzed, and the simulation results are compared with the experimental value, the content of strip theory in potential flow theory and the panel method to prove the reliability of CFD method in predicting ship motion.

Numerical simulation of the heat transfer characteristics of the second and third superheaters in a coal-fired thermal power plant radiant boiler

Superheaters are essential components in a coal-fired power plant because they keep the highest tube wall temperature points in the radiant boiler. The pressurized steam flowing in the superheater tubes is heated via the thermal radiation and convection from the combustion gas. To prevent the bursting of superheater tubes caused by the thermal fatigue for stable and safe plant operation, precisely predicting complex heat transfer characteristics, and estimating the local temperature and heat flux of the superheater are necessary. In this study, a computational fluid dynamics model of the boiler and the second and third superheaters in a coal-fired power plant is developed using radiation model (DO model) and turbulence models (k-ε RNG model). In the second superheater installed near the burner, the influence of radiation from the combustion gas is strong, and the temperature and heat flux of the outermost tube are highest. In the third superheater, the tube wall temperature and heat flux becomes large in the upper and middle parts of the superheater, where the combustion gas flow velocity around the superheater is high, whereas those at the bottom of the superheater are significantly small due to the dead water region of the gas. Therefore, to consider the combustion gas flow is important for prediction of the heat transfer characteristics of the third superheater.

CFD analysis of karanja oil methyl ester (KOME) blend B20 and hydrogen using ANSYS-Fluent

Combustion and emission characteristics are plays an important role in diesel engines, which have higher combustion efficiency meets the stringent emission norms. Combustion is a very complex phenomenon, which cannot be analyzed easily with analytical method or different visualization technique. Computational Fluid Dynamics(CFD) has becoming useful tool in understanding the fluid dynamics of IC Engines for design purposes. In the current study, CFD code is used to perform 3D simulation of mixture formation, combustion and emissions of K20 and K20 with hydrogen fuels simulated by applying the ANSYS-FLUENT non-premixed k-ε RNG model with combustion chamber moving mesh. The parameters like Pressure, Temperature, Parcel diameter, CO, CO2 and NOx generated inside the combustion chamber is compared with K0. The variation in the combustion behavior is find when the fuel is used likeK0, K20 and K20 with hydrogen. Finally, the emission results of simulation are validated with experimental results. The variations of CFD and experimental results ofemissions are find acceptable.