Premixed Combustion Model Research Articles

Numerical Evaluation of Sodium Droplet Swarm Combustion and Its Modeling Optimization

Because of its superior thermal-hydraulic qualities, liquid sodium has been applied to a variety of industries, including energy storage, solar energy, sodium-cooled fast reactors, and aerospace. However, fires brought on by sodium leaks at high pressure can have major thermodynamic repercussions and put employees and equipment in use at risk directly or indirectly. As a result, a realistic and accurate forecast of the combustion behavior of sodium droplet swarm can offer technical backing for the use of liquid sodium in engineering as well as a way of sodium fire prevention and control. Spray dynamics (droplet settling, droplet particle size distribution, etc.), combustion kinetics (premixed combustion, gas phase combustion, etc.), sodium aerosol diffusion, and other specialized phenomena all contribute to the complex process of sodium droplet swarm combustion. The NACOM code created by Brookhaven National Laboratory for sodium droplet swarm combustion is utilized in this paper as a framework. The code is first validated using the benchmark of the sodium droplet swarm combustion tests carried out by prestigious institutions. The validation results demonstrate that the code’s drag model, droplet combustion model, and heat transfer model are to blame for the significantly overestimated thermodynamic effects of sodium droplet swarm combustion. NACOM is subsequently developed twice for the authors’ previously developed vapor-liquid two-layer-structure drag model, chemical kinetic combustion model, and suitable heat transfer coefficient. It is then thoroughly assessed for the separate-effects tests and integral-effects test. The evaluation results demonstrate that the optimized drag model accelerates the settling of sodium droplets due to the consideration of the sodium-vapor drag reduction effect, reducing the thermodynamic effects of liquid sodium combustion; the optimized premixed combustion model can accurately predict the low-temperature sodium droplet swarm combustion conditions, resolving the issue of serious misvaluation of the original version of NACOM. The associated research findings can serve as valuable resources and tools for deeper comprehension of the combustion effects and mechanisms of sodium droplet swarm under various operating settings (such as leakage rate and oxygen concentration).

Modeling of non-homogeneous premixed hydrogen-air flame acceleration and deflagration to detonation transition in an obstructed channel

This paper uses Large Eddy Simulation (LES) technology to study the dynamic behavior of flame acceleration (FA) and deflagration-to-detonation transition (DDT) in a duct with different obstacle layouts in non-homogeneous concentration fields. A turbulent premixed combustion model with three additional source terms was established using the OpenFoam code library. It was found that the predicted results of the coupled model were in good agreement with the measured values of the experiment. Numerical results show that turbulent combustion and self-ignition effects are the main factors affecting flame acceleration. Still, self-ignition effects are a necessary and sufficient condition for inducing DDT phenomena. Detonation flames have around-flow and diffraction effects and do not evolve into jet and vortex flames, which differ from deflagration flames. Two initiation mechanisms were found under different obstacle layouts: hotspots were generated in the area swept by the Mach stem at the flame front, triggering detonation near the turbulent flame through the reactive gradient mechanism; focusing of sufficiently strong reflected waves deposited enough energy in a particular area near the flame brush and directly triggered detonation. In addition, the central obstacle layout is most conducive to FA and DDT, while the symmetrical layout is the least favorable.

Surrogate modeling of parameterized multi-dimensional premixed combustion with physics-informed neural networks for rapid exploration of design space

Parametric optimization is a critical component in designing and prototyping combustion systems. However, existing parametric optimization methods often suffer from either insufficient prediction accuracy or computationally prohibitive data generation costs. This study develops and validates an efficient and robust field-resolving surrogate modeling framework for combustion systems based on physics-informed neural networks (PINNs). This study comprises three progressive steps. (1) We propose a PINNs framework for premixed combustion and validate its accuracy with direct numerical simulation (DNS) for one-dimensional and two-dimensional canonical premixed combustion problems. Results show a robust predictive capability of PINNs with R2 generally above 0.99. (2) We continue to develop a field-resolving surrogate modeling framework for parameterized two-dimensional methane-air jet combustion system with five input parameters: inlet equivalence ratio, temperature, velocity, diameter, and coflow temperature. The PINNS surrogate model after unsupervised training is used to predict the combustion fields for 32 input parameter combinations, and the results are compared against corresponding numerical simulations. The surrogate model accurately predicts the temperature fields with R2 over 0.95, except in case of flashback, where the accuracy is slightly lower. Good consistency between PINNs and numerical simulations is also demonstrated for four combustion performance metrics: flame length, maximum temperature, outlet temperature, and outlet CO2 mass fraction. (3) Based on the validated field-resolving PINNs surrogate model, sensitivity analysis of the parameterized combustion system is performed, and its implications for engineering applications are discussed. Most notably, distinct from common data-driven and machine learning approaches, the proposed PINNs models are entirely physics-driven and do not use any training datasets, saving tremendous amount of time in the design of experiments (DOE) and case-by-case simulations. Overall, the results of this study demonstrate PINNs can be an efficient and robust surrogate modeling approach for parameterized combustion system visualization, design, optimization, and active control.

A novel diffusion flux modeling for laminar premixed flame simulation with OpenFOAM

In this study, we developed a mass diffusion coefficient model for premixed combustion to reduce the computational time, in mixture-averaged diffusivities to serve as a microcomputer resource. In the development, a preliminary theory was derived to construct the associated variables and corrected using a benchmark model as the mixture-averaged model based on Chapman–Enskog binary diffusion. In the validation, the calculation results obtained using the proposed model with a benchmark case, regardless of thermal diffusion and for several fuels at various pressures, temperatures, and initial mixture compositions, showed good agreement. Validation was also performed with the two-dimensional premixed and partially premixed combustion models by implementing the novel diffusion coefficient model in OpenFOAM software. The OpenFOAM solver results were in good agreement with the experimental results. This solver showed good simulation performance at both micro and macro scales. In addition, the processing time was reduced by approximately 70% in all the simulated premixed flame models. Based on these findings, we conclude that the developed model allows researchers to work with more complex premixed combustion modeling with less computational time.

LES Investigation of Kerosene Spray Flame Emission Characteristics in a Staged Combustor

ABSTRACT Staged combustion technology is widely applied in modern low-emission energy and power systems. This study employs large eddy simulations (LES) to investigate the flame dynamics and emission characteristics in a staged two-phase partially premixed model combustor. To improve the suitability of the thickened flame model (TFM) for turbulent flames in the thin-reaction-zone regime, a new model for efficiency function is formulated such that the effects of small-scale turbulence on flame propagation are taken into account. The enhanced LES-TFM approach is then integrated with an efficient NOx model, in which an additional NOx transport equation is solved with pre-tabulated NOx formation source terms. Five LES runs are performed to investigate the evolution of flame dynamics and emission characteristics with staging ratios at the scaled cruise operating conditions, in which the optimal staging ratio is critical for NOx reduction. Results show that the trend of flame morphology and NOx emission predictions with staging ratios are consistent with the experimental results. The flame zones in the staged model combustor shift from a pilot-stage-only flame to a pilot-primary-stage flame, then to a primary-stage-only flame with decreasing fuel staging ratio. Due to the redistribution of the equivalence ratio and the NOx emission characteristics trade-off between the two-stage flows, there is an optimal staging ratio for NOx reduction in the staged model combustor. The proposed LES-TFM-NOx approach, which has both the benefits of efficient flame dynamics modeling with the global reaction mechanism and the advantages of prediction accuracy with the detailed mechanism for NOx emission, demonstrates its advantages in simulating flame and emission characteristics in low-emission combustion devices.

Effects of fuel-hydrogen levels on combustion, operability, and emission parameters of CH4/H2/O2/CO2 stratified flames in a dual-swirl gas turbine burner

Towards the development of emission-free hydrogen gas turbines with wider operability limits, the combined effects of lean premixed combustion, dual combustion (combustion staging), hydrogen enrichment, and oxy-fuel combustion on the combustion and emission characteristics of CH4/H2/O2/CO2 flames are investigated in a dual annular counter-rotating swirl burner. The central stream is of a higher equivalence ratio (∅) of 0.9 for stable flame ignition, whereas the annular stream is of a lower ∅ of 0.65 to maintain the burner’s environmental performance. The stratified flame reactions are modelled using the partially premixed combustion model, which has been validated using the available experimental data. The inner recirculation zone becomes smaller as the hydrogen fraction (HF) becomes higher until it vanishes and there exists only an outer recirculation zone when the secondary hydrogen fraction approaches 100%. For 100% HF in both primary and secondary flow, the maximum combustion temperature within the combustor was observed to be as high as 2350 K. Increasing HF makes the flame more compact and anchored to the burner centre body with faster kinetics, less ignition delay time, and, accordingly, with increased potential for flashback. The product formation rate was found to increase in the inner shear layer at higher HF, implying a higher flame intensity and a more compact flame. The hydrogen fraction of the primary stream (maintained at higher equivalence ratio) is observed to be dominant in regulating CO mole fractions. Stable well-anchored flames are obtained over the considered fuel mixtures, indicating the role of stratified combustion in widening the operability of gas turbine combustors under stratified combustion conditions.

Analysis on roles of thermal radiation to evaporation and combustion of fuel droplets

Thermal radiation is a significant mode of energy transmission besides conduction and convection in most combustion systems. Accurate treatment of thermal radiation is a crucial element for predicting the flame structure and pollutant formation. However, thermal radiation to fuel droplet evaporation is often neglected in the numerical simulations of two-phase spray combustion systems. In this paper, computational analysis on the thermal radiation to evaporation and combustion of fuel droplets in two-phase spray combustion is conducted. The influence of flame radiation on droplet evaporation and combustion is studied by adding radiation source term to the energy equation describing droplet evaporation process. Radiation transfer equation is solved, in which the radiation of CO2, H2O and soot is imposed to consider radiation effect of the flame. The partially premixed combustion model is employed, in which radiation heat transfer between droplet and ambient gases is taken into account. A reduced chemical reaction mechanism for n-heptane consisting of 76 species and 349 reactions is used. Four different numerical setups have been used and compared, and the simulation results coincide with the experiment data in terms of spray statistics and flame shapes. The peak flame temperature increases by 28 K when the radiative heat source term in the evaporation process is considered. The droplet evaporation is significantly enhanced by thermal radiation, leading to a promotion of 34.8% in C7H16 concentration in the center region of combustor. The concentration of combustion products, such as H2O and CO2, increases by 13.2% in the reaction region.

An improved TRF mechanism for a new turbulent premixed combustion model with application to engine combustion CFD

In the present work, an existing TRF (toluene reference fuel) chemical kinetic mechanism has been improved. In the improvement, the existing TRF mechanism was found to have unphysically higher laminar flame speeds at temperatures greater than 750 K which are typical in engine compression-combustion operating conditions. Eight dominant reactions which mainly control the laminar flame speeds under higher temperatures ( T > 750 K) have been identified. Based on a recently published power-law laminar flame speed correlation which is able to provide good predictions under a wide range of engine conditions, correct reaction rate constants for the eight high-temperature dominant reactions have been determined for physical laminar flame speeds. The identification and improvement of the eight high-temperature dominant reactions are the first novelty of this work. In order to simulate ethanol and gasoline blends, a latest ethanol sub-mechanism has been implemented into the TRF mechanism. The improved TRF mechanism was validated using available experimental data in the literature. Based on the improved TRF mechanism, a new mechanism library has been generated, which can be used for an improved mechanism-dynamic-selection turbulent premixed combustion model in which turbulent diffusivity is constructed as a function of local turbulence/thermodynamics conditions. The dynamic turbulent diffusivity sub-model is the second novelty of this work. The improved TRF mechanism and its combination with the improved mechanism-dynamic-selection turbulent premixed combustion model were successfully applied to the prediction of combustion and emissions of GTDI (gasoline turbocharged direct injection) engines under both warm-up and transient cold start operating conditions, indicating that the improved models in this work are beneficial for predictive engine combustion CFD (computational fluid dynamics).

Gas Temperature Distribution and Fluctuation in a Lab-Scale Fire Whirl

Fire whirls cause an increase in fire damage. This study clarified the unsteady behavior of fire whirls, considering that instantaneous changes in the temperature and flame shape of fire whirls can affect the damage to the surrounding area. Numerical simulations of a lab-scale flame that simulates a fire whirl were performed to investigate the changes in gas temperature and velocity fields under various fuel inflow velocities. The flow field was obtained by solving a continuity equation and a three-dimensional Navier-Stokes equation, and the turbulence was resolved using a large eddy simulation. A chemical equilibrium partially premixed combustion model was used, and radiation effects were considered. The time-averaged gas temperature distribution along the burner central axis revealed that the gas temperature decreased monotonically from upstream to downstream. The time-averaged velocity distribution along the burner central axis showed that the velocity decreased as one moved downstream, but the decrease was uneven. The time variation of the gas temperature demonstrated that the higher the fuel inflow velocity, especially near the burner, the greater the gas temperature flutter. Furthermore, the larger the fuel inflow velocity, the larger the flame swell and wobble. The results showed that the fuel inflow velocity affected temperature fluctuation and flame undulating movement.

Multiple Mapping Conditioning Mixing Time Scales for Turbulent Premixed Flames

A novel multiple mapping conditioning (MMC) mixing time scale model for turbulent premixed combustion has been developed. It combines time scales for the flamelet and distributed flame regimes with the aid of a blending function. The blending function serves two purposes. Firstly, it helps to identify zones where the premixed flame resides and where the time scale associated with the premixed flame shall be used. Secondly, it uses the Karlovitz number to identify the turbulent premixed combustion regime and to reduce the weighting of the premixed flame time scale if Karlovitz numbers are high and deviations from the flamelet regime are expected. A series of three-dimensional direct numerical simulations (DNS) of statistically one dimensional, freely propagating turbulent methane-air flames provides a wide range of turbulent combustion regimes for the mixing model validation. The new mixing time scale provides correct predictions of the flame speed of freely propagating turbulent flames which could not be matched by most recognized mixing models. The turbulent flame structure predicted by the new model is in good agreement with DNS for all combustion regimes from flamelet to the thickened reaction zone.

Stratified Flames in Dual Annular Counter-Rotating Swirl Burner for Wider Operability Gas Turbines

Abstract This study investigates numerically the effects of equivalence ratio (Φ) on flow/flame interactions and emissions of stratified oxy-methane (CH4/O2/CO2) flames in a dual annular counter-rotating swirl (DACRS) burner for wider operability and environmental-friendly gas turbines. The flow mixture entering the combustor is split into two coaxial streams of different equivalence ratios. The central stream is characterized by higher Φ to continuously ignite the flame for enhanced flame stability, whereas the annular stream is a highly lean mixture to sustain the environmental performance of the combustor. The partially premixed combustion model is adopted in the ansys-fluent 2021-r1 software to model the reaction kinetics of the generated stratified flames in the two-dimensional axisymmetric domain. Nine cases of the same inlet velocity ratio (primary stream to secondary stream) of 3.0 are examined at a fixed oxygen fraction (OF: volumetric percentage of oxygen in the O2/CO2 mixture) of both streams of 30%. Flame stratification is achieved by varying the equivalence ratios of the primary (Φp = 0.9, 0.8, and 0.7) and secondary (Φs = 0.7, 0.55, and 0.4) streams. The results indicate effective flame/flow interactions, complete combustion, and reduced emissions for the DACRS stratified flames.

Evaluation of flamelet-based partially premixed combustion models for simulating the gas phase combustion of a grate firing biomass furnace

Depending on the mass flow rates of fuel, the primary air (injected usually from underneath of the fuel bed) as well as the composition of volatile gases extracted from the bed, it is believed that the volatile gases taking part in the gas-phase combustion undergo a partially premixed process close to the bed section. Therefore, the aim of this study is to evaluate the prediction capability of flamelet-based partially premixed combustion models in simulating the gas-phase combustion process of a grate firing biomass furnace. Additionally, the effects of the adopted premixed models (i.e., C Equation based model and extended coherent flame model) and non-premixed models (i.e., steady diffusion flamelet (SFM) and unsteady diffusion flamelet (UFM)) on the overall prediction of partially premixed model are assessed. The predicted temperature field and species concentrations are compared with published experimental measurements and also with numerical simulations which use other combustion models (i.e., EDC/Flamelet hybrid model, SFM and UFM). The results of this study revealed that except for slow forming and chemically dominated species, partially premixed combustion models (both extended coherent flame model/SFM and C-equation/SFM based partially premixed model) are capable of reproducing the experimental temperature and major species with reasonable accuracy. For instance, the predicted temperature deviates from the experimental counterpart by a maximum of 9.01%, and the major species by 15.44% (CO2) and 21.60% (O2). In addition, C-equation/UFM-based partially premixed model is found to be the most optimum combination of the examined partially premixed models for overcoming the deficiency faced while predicting the slow-forming species.

Numerical simulation of methane-hydrogen-air premixed combustion in turbulence

3D modeling is based on an experimental platform named turbulence constant volume combustion chamber (CVCC). The turbulence model applies the realizable k-ε model. The combustion model applies the premixed combustion model. Simulation of methane-hydrogen-air turbulent premixed combustion is realized using the model. The initial condition is T0 °C = 298 K, P0 = 1 bar and the gas equivalence ratio is φ = 1. The characteristics of flame are analyzed by modifying the turbulence intensity (uRMS = 0 m/s, 0.7 m/s, 1.3 m/s, 2.0 m/s)and hydrogen fraction (R (H2) = 0%, 10%, 30%, 50%). It is calculated that in the process of spherical flame combustion and propagation, the pressure increase rate increases first and then decreases. Turbulence intensity and hydrogen fraction have similar effect on the flame propagation process. If the turbulence intensity is used as an independent variable, the flame propagation rate increases with its increase. When the flame radius is 30 mm, the flame propagation velocity of 2.0 m/s turbulence intensity is about 2.7 times that of laminar flow. Due to which both the maximum combustion pressure and the maximum pressure rise rate are monotonically increasing. As it increases, the combustion pressure rises and falls the fastest. As the R (H2) increases, the flame propagation rate can significantly increase. When R (H2) is 50%, the maximum velocity is 4.82 m/s, which is about 1.5 times that of the maximum of pure methane. When R (H2) is high, the greater the maximum pressure is, the shorter the time for the maximum pressure to occur will be. When R (H2) is 50%, the maximum pressure is about 9 bar, which is 2.25 times of the force of pure methane combustion.

Simulation of Combustion Flow of Methane Gas in a Premixed Low-Swirl Burner using a Partially Premixed Combustion Model

Simulation of Combustion Flow of Methane Gas in a Premixed Low-Swirl Burner using a Partially Premixed Combustion Model

A multidimensional combustion model for oblique, wrinkled premixed flames

A new premixed turbulent combustion model is proposed, which is based on one-dimensional (1D) filtering of density times progress variable and reaction source term of laminar premixed flame profiles using a filter kernel that reflects the variation of the area of planar flame fronts moving through filter volumes. It is shown that these genuinely multidimensional effects qualitatively change the relation between the filtered reaction source term and the Favre-filtered reaction progress variable compared to strictly 1D models, particularly at large filter widths. Analytical results for the filtered quantities are obtained by approximating the laminar flame profiles by suitable ansatz functions. Filtered data from Direct Numerical Simulations (DNS) of statistically planar turbulent premixed flames at a range of turbulence levels of the unburnt mixture and for two different heat release parameters is used to optimize and validate the model. For filter widths smaller than the laminar flame thickness, minor effects from subgrid flame folding are observed. For larger filters, the effect of subgrid flame folding can be represented well through a wrinkling factor, which reduces the effective filter width used in the 1D filtering operation. Two wrinkling factor models are proposed, using fractal dimensions patterned on literature models. The modelled reaction source term written as function of Favre averaged progress variable and filter width shows excellent agreement with filtered DNS data conditioned on the Favre averaged progress variable for all investigated turbulence levels, filter widths and both heat release parameters.

A Flamelet Based Partially Premixed Combustion Model for Simulating the Gas Phase Combustion of a Grate Firing Biomass Furnace

A Flamelet Based Partially Premixed Combustion Model for Simulating the Gas Phase Combustion of a Grate Firing Biomass Furnace

A Numerical Prediction of Stabilized Turbulent Partially Premixed Flames Using Ammonia/Hydrogen Mixture

The objective of this work is to computationally assess the performance of a carbon free ammonia-hydrogen mixture when burnt in a gas turbine like combustor. Recently, utilizing ammonia as an alternative carbon-free fuel for future power, industry applications and achieving clean energy attracted enormous interest. Pure ammonia oxidation is facing many challenges such as high NOx emissions, high ignition energy, slow reactivity and lower laminar flame speeds. Therefore, the use of ammonia/hydrogen mixture provides flame stability and increasing flame speed. In this manuscript a numerical study for a new swirl stabilized combustor for oxidizing ammonia/hydrogen mixture. Numerical two dimensional model simulations of a turbulent flame on Reynolds Averaged Navier Stokes (RANS) including a realizable k-e turbulent scheme with the aid of chemistry mechanism were performed under various conditions. Partially premixed combustion model with flame-let concept was selected and radiation effects are also considered. Validation for the predicted results showed a reasonable agreement when validated with the experimental data. The results discuss the influence of changing inlet pressure and equivalence ratio on the stability and the characteristics of unburnt NH3 and NO emissions. Results show that for constant operating conditions such as constant equivalence ratio of 0.8 that increasing hydrogen content resulted in increasing NO emission. Also, for constant ammonia/hydrogen concentrations, NO emissions increases with equivalence ratio then reduced at rich conditions and NH3 emissions are generally low. Equivalence ratio lower than 1.2 will be preferable to reduce the amount of unburnt NH3 formation.

FlameFoam: An open source CFD solver for turbulent premixed combustion

Severe accidents in nuclear power plants may cause significant and long-lasting harm to the environment in the case of containment integrity loss. Numerical simulation of phenomena threatening containment integrity, e.g., hydrogen combustion and explosion, is essential for the effective severe accident management.International research projects and benchmarks are organized in order to promote improvement of simulation codes adequacy and accuracy in the severe accident relevant conditions. Results of the most recent benchmarks show that current state-of-the-art simulations were able to mostly reproduce the observable flame speed at the considered conditions, but there are still inadequacies in estimating velocity and its maximum value.Addressing a need for an accessible open source turbulent premixed combustion simulation tool, a new solver is being developed. An OpenFOAM version 7 based combustion solver called flameFoam has been built by implementing a progress variable approach and turbulent flame closure model for premixed combustion simulation. Turbulent flame speed is defined as a function of laminar flame speed and turbulence parameters, and is evaluated using one of three different correlations. Laminar flame speed can be set as a constant or calculated using a laminar flame speed correlation. The basis for the solver are standard compressible OpenFOAM solvers rhoPimpleFoam, buoyantPimpleFoam and chtMultiRegionFoam.The article presents initial public flameFoam version 0.8 and provides examples of its validation. Initial validation and verification of general solver features are represented by the simulations of shock tube and backward facing step cases. Solver capabilities in relation to the current turbulent premixed combustion modeling state-of-the-art are demonstrated with the simulations of ETSON-MITHYGENE benchmark experiments.It is shown that flameFoam version 0.8 is capable of adequate simulation of general flow features, turbulent compressible flow phenomena and turbulent premixed combustion of lean hydrogen-air mixtures at the considered conditions. Solver code is publicly hosted and being developed further.

Study on the flow characteristics and oscillating mechanism of a swirl flame

The larger eddy simulation (LES) method was used to simulate the combustion oscillation using a partially premixed combustion model and coupling with a user-defined function. The simulated results were post-processed by proper orthogonal decomposition (POD) and compared with the experimental results. The results showed that the asymmetry of the upper and lower vortex structures in the inner recirculation zone represented the destruction of the vortex structures and reflected the occurrence of a processing vortex core (PVC). The PVC structure produced by the swirler was more complex than a cyclone. The inner shear layer (ISL) and the outer shear layer (OSL) appeared in the inner and outer recirculation zones respectively, which represented the region with the strongest velocity pulsation; the velocity pulsation of axial centerline could be used to characterize the pulsation energy when the combustion oscillation occurred. First, necking occurred in the inner recirculation zone, and then it adhered to the center of the swirler outlet. The ISL and OSL had strong shearing actions that generating pulsating energy, then they gradually developed downstream and merged into the inner recirculation zone. The different rotation directions of the inner and outer recirculation zones increased the instability of the boundary of the inner recirculation zone. It could also increase the breakdown of the asymmetric vortex structures and exacerbated the occurrence of the PVC. In addition, the phase difference between the heat release rate and pressure pulsation in different regions was different, which leaded to inconsistent driving characteristics of the flame oscillation.

Energy, Exergy and Exergo-Economic Characteristics of Hydrogen Enriched Hydrocarbon-Based Fuels in a Premixed Burner

ABSTRACT The enrichment of conventional fuels by hydrogen is quite popular owing to positive environmental aspects compared to the main fuel. In most of the studies, the entropy analysis was not examined in detail. In this research, exergy (second law) analysis of hydrogen-enriched methane, natural gas, propane, LPG, and biogas were analyzed using the numerical model in a premixed burner due to their wide usage in the industry as well as in the household appliances. In the numerical model, the discrete ordinates radiation and Realizable k- ε turbulence models were coupled with the premixed combustion model. The rate of increase for hydrogen of each fuel is considered as 10%. The premixed burner exhibited better performance in terms of energy and exergy with the rise of hydrogen enrichment. The most improvement was observed at biogas tests. The results showed that the chemical composition of the base fuel is a significant parameter that affects the exergy and energy capability of the burner due to chain reactions of radicals in the base fuel with hydrogen. It can be concluded that the hydrogen enrichment of biogas in premixed burners is promising in comparison with other tested fuels from both energy and exergy perspective.