Effect of multihole baffle-induced lobe flow structures on a high efficiency micro-thermophotovoltaic system

For a combustion-based thermophotovoltaic generator (TPV), the baffle-induced reacting flows are numerically investigated in the micro combustor with a multihole baffle representing a mixed feature of bluff body and multiple jets. The lobed flow feature is realized by the momentum difference between fuel and air streams and the baffle wall, not the geometrical conditions like a lobed nozzle. Compared with the H2–air combustor, the lobed flame and the votical structure are enhanced for the CH4–air combustor. To analyze the lobed structure, the lobe intensity and the equivalent perimeter are defined for the stoichiometric mixture line. The maximum lobe intensity is observed at the beginning zone of flame, whereas the equivalent perimeter is increased near the active reaction zone. As the baffle thickness decreases, the lobe structure is enhanced due to the increased vortices. Also, a close relationship between the maximum vorticities at the baffle exit and the maximized lobe structure is observed. The intensified lobe structure contributes to a shorter flame length and more efficient combustion, resulting in increased radiation available for a micro-TPV system. When the lobe structure is enhanced, the high emitter efficiency for an energy conversion device is obtained as 0.319–0.326 for the H2–air combustor.


Introduction
Recently, with the increasing demand of micro or small devices, various applications of micro combustors are being explored. The combustion-based micropower systems (milliwatts to watts) have been proposed in various forms, including micro-gas turbine, micro-rotary engines, micro-thermophotovoltaic (TPV) and microthermoeclectric generators (TEG). In the systems of micro-gas turbine (Waitz et al., 1998) and micro-rotary engine (Fernandez-pello, 2002), the power is generated directly by using the heat and kinetic energy of combustion gas. Whereas TPV and TEG systems do not use the kinetic energy of combustion gas and convert the combustion heat into the electrical energy through the wall heat transfer. Also, the combined system of microgas turbine and TPV or TEG devices is possible because of being the difference of the power generation form (Kim & Park, 2022). In these systems, the role or design of micro combustor can be changed depending on the form of the micropower system, but the efficient micro combustors are commonly required to improve the system performance whether the system is the form of an engine or the form of a combined system. It is because the micro combustor in micropower systems is a core component to convert the chemical energy of fluids into heat energy (Chou et al., 2011). However, the micro combustor faces difficulties owing to the reduced space for the mixing of fuel and oxidizer, increased heat loss, insufficient residence time, weakened turbulent fluctuation, and inefficient combustion. In general, these problems can be resolved by creating recirculation flows in the combustion chamber. The recirculating flows improve the mixing and the residence time of fuel and oxidizer. It is possible to develop these features inside the combustor by changing the combustion chamber based on bluff bodies (Bagheri et al., 2014;Fan et al., 2013;Zhang et al., 2015), cavities (Wan et al., 2015), baffle plates (Choi et al., 2008;Choi & Park, 2009;Kim & Park, 2018, and backward-facing steps (Akhtar et al., 2019). The main idea of these studies is to change the main streams using the secondary flows depending on the combustor shape. Among them, the baffle plate with multiple holes has a mixed feature of bluff body and multiple jets, because the baffle wall wake is considered as the flow near a bluff body and the flow past the baffle hole is regarded as a jet. Therefore, if the flow passes through the baffle plate, threedimensional recirculations can naturally occur due to the streamwise momentum difference between jet flows and the baffle wall. According to Park (2018, 2019), two large vortical regions called wall and center recirculation are formed by the interaction between center jet and annularly placed jets for multihole baffled micro combustors. And, the lobed flow structure of the streamwise velocity is commonly observed. In the case of the combustor having a high combustion efficiency, such lobed flow becomes more and more clear. This feature of the micro combustor with the baffle plate can be useful to control various potential possibilities for an energy conversion device. Accordingly, the generation of recirculating zones and the development of the lobed flow can be an important feature of the baffled micro combustor, but generation conditions and secondary flow characteristics related to the lobed flow structure have not been closely studied in a micro combustor.
From a literature survey, the lobed flow structure is dominated by the secondary flows due to the lobed shape like lobed mixer and lobed nozzle. It has been widely used to get efficient combustion, reduced pollutants, and reduced jet noise. Early studies are concentrated on lobed types comprised of corrugated surfaces like sinusoidal and scalloping shapes to understand the geometrical feature for mixing enhancement. Smith et al. (1997) conducted the experimental studies of the nonreactive mixing processes in a lobed fuel injector, which consists of two corrugated plates. The lobed injector showed a 30% increase in mixing performance compared to a non-lobed injector. Hu et al. (2000Hu et al. ( , 2001Hu et al. ( , 2002 studied exhaust ejector/mixer systems based on lobed nozzles with rectangular mixing tubes. The lobed nozzles enhance the mixing process between the hightemperature and high-speed gas plume and ambient cold air. Such lobed nozzles generate large-scale streamwise vortices and highly strained flowfields for enhanced mixing. Also, Cai et al. (2010) studied the generation and evolution of vortical flows in the jets from notched circular nozzles. The notched nozzle jets show the mixing enhancement due to the development of secondary flows. Their study found that the variation of streamwise vortices and secondary flows is similar to the lobed flow structure. Based on these studies, we deduce that the distorted azimuthal vortices in a jet flow induce additional streamwise vortices. Many methods can generate additional vortices to the mainstream. In the case of the baffled micro combustor, it is interesting that the lobed flow structure can be implemented by the momentum difference of fuel and air jets through the baffle plate.
When the lobed flow develops, the secondary motions are strengthened than the non-lobed flow. The central fluid moves outward by such secondary flows, and then the ambient fluid flows into the central region. Streamwise vortices are also generated as a result of these features, and mixing is improved. In the baffled micro combustor, the special characteristics of the lobed flow are mainly attributed to the baffle hole configuration, whereas the general lobed flow is dominated by the noncircularity degree of corrugated surfaces. However, the relationship between the hole condition and the generation of the lobed flow has received little attention. The related results may provide useful information for expanding the production method of the lobed flow structure using a flat plate with multiple holes. Furthermore, many studies on lobed flow have concentrated on the flow mixing of hot and cold fluids, but few studies have focused on reacting flow coupled to lobed flow. Since the combustion reaction is the important phenomenon of various applications, understanding reacting flows based on the lobed flow structure would be very meaningful to get an efficient combustor.
In the present study, the influence of the lobed flow structure on turbulent combustion characteristics is investigated when the baffle plate with multiple holes is introduced in a micro combustor. The previous study (Kim & Park, 2018) shows that large recirculation regions are mainly changed by the baffle hole conditions of diameter, shape, and position. As a result, the lobed flow structure can be differently obtained according to the feature of baffle configuration. Meanwhile, the typical flow structure in the lobed jet flows is determined by the inflow conditions such as number, surface shape, and penetration angle of lobes (Paul et al., 2020;Sheng et al., 2015). Based on these findings, the geometrical condition of the baffle plate has a great influence on the axial flow patterns because the baffle exit can be considered as the inlet of the combustion chamber. Therefore, the inflow condition of the baffled micro combustor is controlled by changing baffle shapes. The main shape variable is chosen to be the baffle thickness. Changing the baffle thickness results in a variety of jet flows past the baffle hole. The jet flow has orifice characteristics for the thin baffle and fully developed pipe characteristics for the thick baffle. As a result, it is worthwhile to investigate the effect of baffle thickness on reaction flows. For combustion with air, the common fuels hydrogen and methane are chosen. In lobed jets, the lobe structure significantly affects the flow and thermal field because of adding streamwise vortices. Therefore, it is necessary to examine the coupling of the lobed flow structure and two large vortical regions of wall and center recirculation. From the results, the variables related to the lobed flow structure are discussed for different fuels. The lobe structure of the stoichiometric mixture line is investigated for reacting flows to explain combustion characteristics for different baffle thicknesses. The lobe intensity and the equivalent perimeter of such lobed stoichiometric line are introduced. Various parameters of the flame length, combustion efficiency, and wall heat flux are discussed with the axial variations of the lobe intensity.
Based on such lobed flames, the present micro combustor is combined with micro-TPV as shown in Figure 1. A typical TPV system consists of three parts: a heat source, an emitter, and a photovoltaic cell with the low bandgap energy (0.4-0.7 eV). When various types of heat energy such as solar heat, combustion heat, nuclear heat, waste heat etc. are provided, the radiant heat is emitted from the emitter. The bounded electron in the PV cell is excited by higher thermal radiation than the bandgap energy. Figure 1(b) shows the combination of present micro combustor and micro-TPV. According to Wien's law, as the radiation heat increases, the bandgap energy is reduced. Therefore, as the emitter temperature increases, the performance of the TPV system can be improved.
When a combustor is adopted as the heat source, the emitter temperature is closely related to the combustion efficiency. The present combustor can get various heat transfer rates depending on changes of the lobed flame structure. Based on such flame and heat transfer characteristics, the feature for an energy conversion device is explored.

Governing equations
When the Reynolds number in the micro combustor with a bluff body is greater than 500, a proper turbulence model is required to predict flow mixing and combustion characteristics (Kuo & Ronney, 2007). As previously shown in the study (Choi et al., 2008;Kim & Park, 2018, the multihole baffled micro combustor has the interaction of multiple jets and various vortical structure near the baffle plate of x < 5D. The increased strain rates become the significant source of turbulent kinetic energy. Therefore, the turbulent simulation is required even if the inflow has laminar flow characteristics. The Reynolds stress model (RSM) is adopted for better prediction of three-dimensional turbulent recirculations. From Kim and Park (2022), the baffled micro combustor has unsteady features and it is dominated by a monotone frequency of three-dimensional vortices. They showed that the steady field is almost similar to the time-averaged field. Therefore, the study on various shapes of micro combustor is based on the time-averaged field, because the investigation of dynamic vortices is not purpose. For steady-state and incompressible turbulent reacting flows, the governing equations of the continuity, momentum, mass fraction of species i, enthalpy, and Reynolds stress are as follows where ρ, U i , u i u j , τ ij , P, Y i , D i,m , R i , h, T, D ij , P ij , φ ij , and ε ij are density, velocity, Reynolds stress, stress tensor, pressure, mass fraction of species i, diffusion coefficient, chemical reaction source, enthalpy, temperature, diffusion of Reynolds stress, production of Reynolds stress, pressure-strain correlation, and anisotropic dissipation rate, respectively. The anisotropic dissipation rate is obtained by ε ij = 2ρεδ ij /3, where δ ij is the kronecker delta. The turbulent viscosity μ t is calculated as μ t = ρC μ k 2 /ε (C μ = 0.09), where the dissipation rate of turbulent kinetic energy (ε) is computed from the transport equation. The turbulent Schmidt and Prandtl numbers are Sc t = 0.7 and Pr t = 0.85, respectively. For simplicity, details of turbulence model functions and constants are omitted. To examine the flame characteristics preferentially, the radiating heat from the combustion gases in the flame to the combustor wall is neglected (Kim & Park, 2018;Zhang et al., 2015). The governing equations are solved by ANSYS 13.0 (2006). The SIMPLEC algorithm for the pressure-velocity coupling and the second-order upwind scheme for the convection terms are adopted. The convergence criterion is set as 10 −6 .

Computation domain and boundary conditions
The computational domain of the micro combustor with a multihole baffle plate is depicted in Figure 2. The baffle has a fuel hole in the center and annular six air holes around it. According to Mi and Nathan (2010) and Park (2020, 2021), the inflow types of a smooth nozzle, a sharp-edged orifice, and a long pipe significantly affect jet development and axis-switching flow. Each inlet has its own set of secondary flow components, and these secondary flows at the inlet dominate the development of the jet flow. In this study, the flows past the baffle plate are considered as the inflow of the combustion chamber.
The nature of such inflow can be changed by the baffle thickness (b t ). The fuel and air jets passing through the thin baffle have orifice inflow characteristics, whereas the thick baffle jets have long pipe inflow characteristics. Based on this, the variation of the baffle thickness is tried to get different fuel and air jets. The b t change induces different typical vortices because the inlet of the micro combustor is differently assigned by such change. Also, the lobed flow structure is related to the secondary flow development in the combustor cross-section. Accordingly, we comprehend that the b t change is meaningful to  see the lobed structure and its effect on the flame structure. In terms of baffle design, various parameters such as the hole diameter, the number of air holes, their position, etc., should be considered for efficient combustion. However, for reacting flows, it is difficult to distinguish the characteristics of one variable while considering many variables at the same time. Therefore, a baffle configuration of seven holes is selected based on the previous studies of Park (2018, 2019). In the preliminary process, the combustor length is changed for L = 50D f , 100D f 150D f , and 200D f . As a result, the length is determined to be more than the maximum flame length. The geometrical condition is summarized in Table 1. The baffle thickness is changed for b t /D f = 0.01, 0.05, 0.1, 0.3, 0.5, 0.75, and 1.5. The constant mass flow rate is applied for the inlet boundary. At the inlet, the mass flow rates of fuel and air are an important factor in determining the flame pattern. The global equivalence ratio based on the inlet mass flow rate is set to 1. It is calculated by φ G = AF st /AF where AF st denotes the stoichiometric air-fuel mass flow ratio and AF is the inlet air-fuel mass flow ratio m air /m fuel . Then, the inlet Reynolds number based on the combustor diameter (D = 2.8 mm) is Re D = ρ m U m D/μ m = 500, where ρ m , U m, and μ m are the mean density, streamwise velocity, and dynamic viscosity considering the inflow of fuel and air. For φ G = 1, AF st = AF = ρ air A air U air /ρ fuel A fuel U fuel . Here, ρ, A, and U are density, cross-sectional area, and streamwise velocity at the inlet, respectively. The air velocity is expressed by The mixture density satisfies the relation of ρ air m air + ρ fuel m fuel = ρ m (m air + m fuel ) based on the fluid masses. After replacing the fluid mass with the mass flow rate, the mean density ρ m = (ρ air AF + ρ fuel )/(1 + AF) is calculated by using AF st = 17.11 for the CH 4 -air combustion and AF st = 34.11 for the H 2 -air combustion. The mean viscosity μ m is similarly calculated. First, U m is obtained by Re D = 500. From a combination of U air = AF st (ρ fuel /ρ air )(A fuel /A air )U fuel and U m = (U air AF + U fuel )/(1 + AF), the velocities of air and fuel are determined. In Table 2, velocities, momentum fluxes, and Reynolds numbers at the inlet are presented. To get the spatial distributions of velocity and turbulence Notes: Re fuel = ρ fuel D f U fuel /μ fuel , Re air = ρ air D a U air /μ air , M fuel = ρ fuel πD 2 f U 2 fuel /4, M air = ρ air πD 2 a U 2 air /4, ρ CH4 = 0.6517kg/m 3 , μ CH4 = 1.09 × 10 −5 kg/m · s, ρ air = 1.1709kg/m 3 , μ air = 1.7219 × 10 −5 kg/m · s.
variables, the inflow generator of Kim and Park (2021) is adopted as Figure 2. The inlet temperatures of fuel and air are 300 K. At the outlet, a constant pressure condition of P = P out is applied and the others are used by the Neumann boundary condition. For wall boundaries, the enhanced wall treatment for turbulence modeling and zero fluxes of species are adopted (Khalil et al., 2019;Sui et al., 2021). The reacting flow and heat transfer are analyzed by the conjugate heat transfer including the baffle plate and the combustor wall. The temperature field of the solid zone is calculated by solving λ w ∂ 2 T/∂x 2 i = 0, where λ w is the wall thermal conductivity. The balance of heat fluxes at the solid-fluid interface is treated by where λ fluid is the fluid thermal conductivity. The resulting temperature is commonly assigned at the boundaries of solid and fluid zone. The boundary conditions and the wall property are summarized in Table 3. To estimate the heat loss from the combustor wall, the combustor wall is treated by the convective and radiative heat transfer. The wall heat flux of the combustor wall is an important factor for the development of an efficient micro combustor. The growing heat loss means a reduction in the stability and efficiency of the combustion chamber. However, to configure energy conversion devices, the increased wall heat transfer is a better condition. Therefore, such conflicting characteristics need to be reviewed for the change of baffle shape.

Combustion model
To get reacting flows, the reaction mechanism is selected as 325 reactions of GRI-Mech 3.0 (1999) for CH 4 -air combustion and 19 reactions of Giovangigli and Smooke (1987) for H 2 -air combustion. The turbulent flame structures of non-premixed combustions are investigated by the eddy dissipation concept (EDC). The mixture fraction Z is the local mass ratio of material originating from the fuel stream to the mixture of fuel and air. In general, the flame structure is estimated by the stoichiometric mixture fraction Z st , because the chemically complete combustion is made in the vicinity of the stoichiometric condition of Z = Z st . Based on this feature, the flame length L f is obtained from the position of the stoichiometric condition Z st along the centerline (Kim & Park, 2018). The stoichiometric conditions of hydrogen and methane reactions are obtained by two global reactions The mixture fraction is locally related to the mass fractions of the fuel and oxygen.
Here, ν is the stoichiometric mass ratio ν = ν s The stoichiometric coefficients are ν s O 2 = 2 and ν s fuel = 1 for CH 4 -air combustion and ν s O 2 = 0.5, ν s fuel = 1 for H 2 -air combustion. Y fuel,1 denotes the mass fraction of fuel in the inlet fuel flow and Y O 2 ,2 is the mass fraction of oxygen in the inlet oxidizer flow (e.g. Y O 2 ,2 = 0.232 for air). For a stoichiometric mixture, the fuel and oxidizer are completely consumed and νY fuel = Y O 2 and Z st = Y O 2 ,2 /(νY fuel,1 + Y O 2 ,2 ) are obtained. From this definition, the stoichiometric mixture fraction is calculated by Z st = 0.232/(4 + 0.232) = 0.0548 for CH 4 -air combustion and Z st = 0.232/(8 + 0.232) = 0.0282 for H 2 -air combustion. The basic characteristics of CH 4 -air and H 2 -air reactions are summarized in Table 4. The gas density is calculated by the ideal equation of state at 101,325 Pa.

Model validation
Before proceeding further, the grid dependency is performed for five cases. The number of control volumes (CV) is changed by 200,000 (i = 1), 400,000 (i = 2), 600,000 (i = 3), 800,000 (i = 4), 1,000,000 (i = 5). The near wall modeling of turbulence equations is based on the enhanced wall treatment (Espinoza-Jara et al., 2022). To solve the viscous sublayer, the first grid points near the wall are located at y + ∼ 1, where y + is the maximum of non-dimensional normal lengths at the first grid from the wall. If y + ∼ 1 is not satisfied when the converged solution is obtained, the grid spacings are readjusted. The minimum axial spacing is almost maintained at x/D = 0.005 and the maximum axial spacing is changed for x/D = 0.042 (200,000 CVs), 0.036 (400,000 CVs), 0.028 (600,000 CVs), 0.023 (800,000 CVs), and 0.016 (1,000,000 CVs). In Figure 3, the relative error of flame length is presented for different grids. Here, the error is defined as As the grid resolution increases, the relative error is gradually decreased and the flame length approaches to a constant value with a specific flame pattern. For the grid larger than 800,000 CVs, the difference between flame lengths is less than 1%. Therefore, the final grid resolution is maintained at more than 800,000CVs.
Since there was no data for the flow field of the multihole baffled micro combustor, the flow and temperature fields of jet flame, bluff-body burner, micro cylindrical combustor, micro can combustor were chosen for comparison (Kim & Park, 2018. Similarly, to check the adopted turbulence and combustion models, several calculations are performed for the experiment of Veríssimo et al. (2015) because their combustor is similar to the present baffled micro combustor. As shown in Figure 4(a), the combustor consists of a central air orifice and around 16 methane orifices. The streamwise velocity variation from the multiport inlet is similar to that of the multiple air jets past the baffle plate. The thermal field of a combustor is nonlinearly coupled to the flow and turbulent field. The combination of combustion and turbulent models needs to be examined for the change of the temperature field. In Figure 4(b), temperature distributions of the standard k-ε model (SKE), the shear stress transport k-ω model (SST), and the RSM model are compared. Here, D a is the diameter of the air orifice. For the RSM model, the turbulent kinetic energy k is defined by The mixing of fuel and air is enhanced near the injection surface because the fuel and air through the multihole plane are provided to the combustion chamber. The resulting temperature rise is observed for x/D a < 14. The RSM model well predicts it. For the CH 4 -air combustor, the combustion model is tested by the 2-step EDC (Westbrook & Dryer, 1981), the 6-step EDC (Jones & Lindstedt, 1988), the 41-step EDC (Yang & Pope, 1998), and the detailed mechanism of GRI-Mech 3.0 (Smith et al., 1999). As can be seen in Figure 4(c), the GRI-Mech 3.0 mechanism gives higher accuracy than other mechanisms. Based on these results, it is considered that the combination of RSM and the GRI-Mech 3.0 mechanism can be effectively used to capture the multihole effect on reacting flows in the baffled micro combustor.

Flow characteristics
As explained in Park (2018, 2019), the multihole baffled micro combustor generally has two large recirculation zones derived by the baffle plate. In Figure 5(a), the recirculation zone characteristics of x-y and x-z planes are explained by velocity vectors, streamlines, and contours of U/U c . As the fuel and air jets progress downstream, the shear layers grow radially. Such development of jets entrains the ambient fluids. As a result, the wall  recirculations like backward-facing flows are generated near the region between fuel and air jets and the combustor wall (see blue lines). In the central region, the fuel flow after the baffle plate shows the general jet flow. When the axial momentum of the fuel jet is weakened, the fuel flow is entrained by the air flow. As a result, the center recirculation is formed at the combustor center region in the front of the fuel jet (see red lines). Such flow characteristics in the baffled micro combustor can be identified in the streamlines and velocity fields. Figure 5(b) shows contours and axial distributions along the centerline of streamwise velocity. The vortical regions are additionally represented by streamlines. Two combustions of CH 4 -air and H 2 -air are selected, and U c is the center velocity at the exit of the baffle fuel hole. The general feature for the axial development of U/U c is observed.
The negative region of U/U c is formed for x/D f = 5-9. For the CH 4 -air combustor, the center recirculation is developing a little larger in spite of the same shape. The stoichiometric air-fuel ratio is 17.11 for the CH 4 -air combustor and 34.33 for the H 2 -air combustor. Considering the CH 4 density greater than H 2 under the constant mass flow rate, the CH 4 velocity diminishes to 1/8 of the H 2 velocity. This momentum change of the fuel jet is considered to be the main factor in the difference of the center recirculation zone. Therefore, the geometry effect for the flow field needs to be considered with the fuel density and the corresponding momentum. According to Kim and Park (2018), as the negative region of U/U c is enlarged, better combustion is obtained. Based on their research, the above result shows that the heavy fuel gives a more effective combustion condition under the same geometrical condition and the inlet of the fixed mass flow rate.
The developmental patterns of reacting flow can be examined by recirculating flows and temperature field. Figure 6 shows streamwise velocity contours, velocity vectors, and temperature contours on the x-y and x-z planes. In Figure 6(a), the x-y plane passes through the air hole at z = 0, and the x-z plane is the side between the air holes at y = 0. Two large recirculation zones are observed. As can be seen, the wall recirculation zone is related to the air jets, and the center recirculation zone is formed in the fuel jet front. The flame zone starts at the beginning region of the center recirculation zone, and the maximum temperature region is formed a little downstream. The wall recirculation zone produces reversed air flows, and the fuel flow is bifurcated by the center recirculation zone toward the air flows. As the residence time and the reactant mixing is improved by this flow structure, the reacting region becomes an elongated candle flame. Therefore, the efficient mixing and combustion in the baffled micro combustor are significantly affected by wall and center recirculating regions behind the baffle plate. In particular, the center recirculation zone enhances the fuel-air mixing by promoting the movement of the fuel toward the air. In Figure 6(b), the stoichiometric mixture line of Z = Z st is plotted with temperature contours to see the flame structure. For various geometry and boundary conditions, the stoichiometric mixture fraction is generally adopted to interpret the flame structure in the Z space. Four axial positions of x/D f = 1, 4, 10, 20 are selected to check the flame development. The velocity vectors of secondary flows are classified into two parts. The red vectors represent the fuelrich mixture of Z > Z st, which is originated from the fuel stream. And the blue vectors are the fuel-lean mixture of Z < Z st originated from the air stream. At x/D f = 1, the lobed Z = Z st line develops weakly. Since the axial position is very near the baffle, the secondary flows and the fuel-rich stream are weak. But the lobed features of the Z = Z st line and the fuel-rich stream are observed at x/D f = 4. This is considered that the lobed flame shape is actively developed. In the downstream x/D f = 10, the lobed flame is very weak, and after that, it is completely disappeared. These changes in the lobed Z = Z st line and the fuel-rich zone are attributed to the development of the wall and center recirculations. The lobed flame shape is roughly maximized at the streamwise position x/D f ≈ 4 corresponding to the center of the wall recirculation or the beginning of the center recirculation. When the wall recirculation zones develop, the spreading of the fuel-rich stream is enhanced than the case of no wall recirculations because the central stream tries to compensate for the mass reduction by the reversed flow of the wall recirculation. Also, as the center recirculation develops, the radially outward flow of the fuel-rich stream is enhanced. As a result, it can be comprehended that the lobed flame shape is developed. The lobed structure enhances the mixing of fuel and air due to the secondary flows of moving the air to the center of the combustion chamber and the fuel to the combustion chamber wall. The magnitude of such secondary flows is reduced while passing through the flame region, because the streamwise flow is accelerated across the flame region. Then, the lobed structure is weakened, and the cross-sectional flame shape is formed into a circular shape. Based on these characteristics, the baffle-induced lobed flame significantly influences the initial development of the flame zone. As a result, the developing pattern of the lobed flame is critical in altering the efficient combustion and heat transfer characteristics of the baffled micro combustor.
The secondary flows associated with the lobed structure in the cross-section of the combustion chamber can be examined by the streamwise and azimuthal vortices. In lobed jets, the flow mixing is improved by the increased such vortices (Hu et al., 2000(Hu et al., , 2001(Hu et al., , 2002. From a similar perspective, it is meaningful to examine the lobed structure and the related vortices. For reacting flows of CH 4 -air and H 2 -air, Figure 7 shows the distribution of streamwise (ω x ) and azimuthal vorticity (ω o = (ω y 2 + ω z 2 ) 1/2 ) on the y-z plane at x/D f = 1.0, 2.0, 4.0, and 8.0. In Figure 7(a), streamwise vortices are axially developed in the form of six pairs. At x/D f = 1.0, the high vorticity is concentrated near air holes. As the flow develops downstream, the region of the high vorticity is expanded toward the fuel hole and combustion chamber wall. In Figure 7(b), the high azimuthal vorticity at x/D f = 1.0 is mainly observed near the fuel and air holes. At x/D f = 2.0, 4.0, the azimuthal vorticity is more and more widely distributed around the holes. As a result of such changes, counter-rotating streamwise vortices are generated at x/D f = 2.0 and 4.0. These vortical motions are related to the wavy development of the stoichiometric mixture line. Regardless of the reaction type, the development pattern of streamwise and azimuthal vortices is  almost similar, but the vortex intensity of CH 4 -air combustor is higher than that of H 2 -air combustor. Therefore, it can be considered that the vortical motions due to the multihole baffle increase with high fuel density.
From the literature on the lobe jet (Cai et al., 2010;Hu et al., 2000Hu et al., , 2001Hu et al., , 2002Paul et al., 2020), it is well known that the streamwise vortices in the lobe crest or trough have an important role in enhancing the flow mixing. The baffle-induced secondary flows are very similar to the flow structure of the lobed jet. However, the present lobed structure is generated by the momentum difference between fuel and air jets from the baffle hole and the baffle wall, unlike the vortices of the lobed jet due to the specific shape like the corrugated surface. As a result, secondary flows in the baffled micro combustor develop more slowly than secondary flows in the lobed jet. However, the lobed jet has the feature that the secondary flow immediately develops, because it is based on the specific nozzle. Considering this difference, the axial variation of the streamwise vorticity magnitude is compared to the experimental results of Hu et al. (2000Hu et al. ( , 2001Hu et al. ( , 2002. The streamwise vorticity magnitude is an important feature representing the lobed flow structure, but it is difficult to directly compare the vorticity magnitude for different geometries. In Figure 8, non-dimensionalized vorticities are compared. It is based on that the axial development of the lobed flow is similar even if the lobed flow can be obtained from various methods. Although the production process of the lobed flow structure is different, the two results show that they have very similar developmental process. From the results, we understand that the baffle-induced lobed flow has a secondary flow structure that is very similar to the lobed jet.

Development of lobed flames
To examine the development of lobe flame structure, in Figure 9, contours of temperature and OH mass fraction at several streamwise locations are plotted with isosurfaces of Z = Z st . The black line of Z = Z st is inserted in the contour plots. In general, the high-temperature zone and the large OH (hydroxyl radical) mass fraction are formed near the stoichiometric condition of Z = Z st . The large OH mass fraction represents the front flame region with the peak temperature. Therefore, the change of the flame shape can be analyzed by the pattern of the Z = Z st line and the distribution of the OH mass fraction. In the figure, the flame length and the flame surface are understood by the isosurface of Z = Z st . The flame length of CH 4 -air is decreased to 60% in the flame length of H 2 -air. As shown in Figure 5, the CH 4 -air case has a larger center recirculation than the H 2 -air case. The center recirculation induces the movement of the fuel toward the air. Eventually, the fuel-air mixing gets faster and the flame length decreases due to the increased center recirculation. The change of the center recirculation zone is strongly dependent on the baffle hole conditions of diameter, shape, and position (Kim & Park, 2018). For example, when the air hole position (r a ) changes, the flame length of the 2r a /D = 0.643 case is reduced to 58% of the 2r a /D = 0.429 case. And, the triangular hole baffle induces a 40% reduced flame length compared to the circular hole baffle. In the present result, two fuel types of CH 4 and H 2 are compared under the same baffle plate. Interestingly, the flame length change due to the fuel type is comparable to the flame length reduction according to the geometrical change of the baffle plate. In Figure 9(b), it is confirmed that the flame shape of the Z = Z st line becomes a lobed structure. In the lobed Z = Z st line, as the distance difference between crest and trough increases, it can be determined that the lobed flame is enhanced. For the CH 4 -air flame, such difference is more severe than that of the H 2 -air flame at the same streamwise position. Also, the reaction zone is checked from the radical distributions of OH in Figure 9(c). In the reactive zone, radicals are created to maintain the chemical reactions, and the reaction heat is emitted with the recombination reactions using the radicals. The Z = Z st line represents the boundary of the high-temperature zone, and the large mass fraction of OH is observed in the high-temperature region. As can be seen, the radical mass fractions increase near the crest of the Z = Z st line, and the temperature increase. Based on these features, the baffle-induced lobed flame can be identified. However, we conclude that the analysis of the lobed flame is appropriate to be based on the Z = Z st line because it is difficult to determine the flame boundary by the distributions of temperature and radical.
For reacting flows, the change of the flow structure mixes the fuel and air differently, and the different reactions are eventually obtained. Such variations of reacting flows are reflected in the distributions of temperature and mass fractions of reactants and products. Figure 10 shows axial distributions of temperature and mass fractions of reactants and products. To help the understanding of axial temperature variations, streamlines in the x-y and x-z planes are added with temperature contours. Based on the temperature rise, the combustion reaction is slowly initiated at x/D f = 2-4. This location is consistent with the position where the streamwise and azimuthal vorticity distributions begin to become the lobed pattern (See Figures 6 and 7). After passing through the location, the temperature drastically increases because of the increased reaction rate. As can be seen, the CH 4 consumption is axially faster than the H 2 consumption. For the H 2 -air case, the radial flow development and the reactant mixing become weak because of the reduced center recirculation zone. As a result, the H 2 -air case is less activated, and the reaction zone is positioned further downstream than that of the CH 4 -air. The steep increase of temperature appears at x = 4-6 D f for the CH 4 -air case and x = 4-10 D f for the H 2 -air case. These regions approximately correspond to the axial position of the center recirculation. For different fuel types, the wall recirculation size is weakly changed, but the center recirculation is more severely changed. The change of the center recirculation zone depending on the fuel property is comparable to the variation due to the geometrical change of the baffle plate (Kim & Park, 2018. Therefore, it is considered that the center recirculation has a crucial role in the change of combustion characteristics. In this sense, the lobed flame structure is strongly influenced by the development of the center recirculation. In a baffled combustor, the production of recirculating flows depends on the geometrical conditions of the baffle plate, and different recirculations are closely related to changes in the lobed flow structure. Among the geometrical variables, when the air hole position (r a ) and the number of air holes change, the distributions of temperature and OH mass fraction and streamlines of secondary flows are compared in Figure 11. To see the development of the lobed flame shape, the fuel-lean regions of Z < Z st and the fuel-rich regions of Z > Z st are plotted by blue and red vectors, respectively. When r a changes, the lobed flame structure is enhanced for the 2r a /D = 0.643 case than the 2r a /D = 0.471 case. For the 2r a /D = 0.471 case, the fuel-lean region increases and the flame length is reduced due to the degraded mixing efficiency. Meanwhile, the number of air holes is changed. The 5-hole case is developing the lobed flame structure just like the 6-hole case. Compared to the maximized lobe shape, the wavy shape of the lobe is sharp for the 6-hole case, whereas it is blunt for the 5-hole case. However, the region satisfying the Z = Z st condition for the 5-hole case is more than the 6-hole case at the same streamwise location. This means that the well-mixed state is quickly obtained for the 5-hole case. It is confirmed by the flame length comparison. Considering the well-mixed state, getting the Z = Z st condition in as many places as possible is better. Therefore, the lobed Z = Z st line represents a wellmixed mixture for efficient combustion than the circular Z = Z st line. Based on flow characteristics, the different lobes of the lobed Z = Z st line can be obtained when the number of air holes is changed. But many lobes do not necessarily indicate a well-mixed state because the difference between crest and trough of the lobed Z = Z st line may be small. Therefore, an appropriate parameter representing such a situation needs to be defined.
From a closer inspection of Figures 6-7 and 11, the mixing structure of fuel and air is systematically explained in Figure 12. In the present combustor, the fuel hole is located at the center of the baffle plate, and six air holes surround it. The fuel-air mixing and the corresponding mixture fraction are determined by the geometrical feature of the baffle plate. The lobed shape of the Z = Z st line, i.e. lobed stoichiometric line or lobed flame, is observed due to the interaction of radially outward fuel flow and inward airflow. As can be seen, the lobed Z = Z st line has six crests and troughs because the baffle has six air holes. As discussed in Figure 11, if the number of air holes changes, the number of crests and troughs will be changed by it. The secondary flows are divided into a fuel-rich stream of red vectors and a fuel-lean stream of blue vectors, as shown in Figures 6  and 11. First, between the air jets and the combustion chamber wall, recirculating wall flows are created. Due to the constraint of the baffle plate and combustion chamber walls, reversed flows produce radially inward flows. The fuel jet exhibits azimuthally differing spreading rates in the interval region of air jets and non-air jets. In the azimuthal interspace of air jets, the radial spreading of the fuel jet is more free because the influence of the air jets is small. Therefore, as the fuel stream is pushed to the combustion chamber wall, the radial outflows are developed between air holes, and the crest pattern of the lobed structure is obtained. Because of this flow, the inward flows near the air jets become stronger than those of the azimuthal region between air jets. As a result, the trough pattern of the lobed structure is obtained. Based on these flow structures, the flows pushed out of the crests are directed to the air holes, and the counter-rotating vortex pairs are generated near both sides of crests. Since the production of the lobed flame structure is mainly attributed to the radially outward motion of the fuel stream, the existence of the center recirculation zone is highly critical to the development of lobed flame. The above pattern of secondary flows and six pairs of counter-rotating streamwise vortices are confirmed from the velocity vectors and contours of streamwise vorticity. Comparing the CH 4 -air case to the H 2 -air case, the lobed structure is significantly developed, and the streamwise vorticity is enhanced. However, their overall feature is commonly observed regardless of the fuel type. Therefore, the lobed flame can be considered as a typical feature of the multihole baffled micro combustor.

Shape parameter of lobed flame
The lobed flow and flame development are closely related to the variation of recirculating flows in the combustion chamber. For b t /D f = 0.1, Figure 13 shows the relation between the axial change of the lobed flame and wall friction coefficient on the combustion chamber. To quantify the lobed flame structure, the lobe intensity I lob is newly defined for the lobed stoichiometric line. Here, I lob = a/b, where a and b are radii to the crest and the trough of the lobed Z = Z st line, respectively. Also, we can determine the size of the wall recirculation zone from the negative region of the wall friction coefficient. From the I lob definition, the increased I lob means an enhanced lobed shape, while I lob = 1 represents a non-lobed case of circular shape. I lob increases to x/D f = 4.0-5.0, then decreases. Most changes of I lob are occurring in the region of x/D f < 8.0 where the wall recirculation is developing. Based on the characteristics of I lob change, the lobed position X lob is defined, and the maximum lobe intensity I lob,max is obtained. X lob is estimated from the axial position of the maximum lobe intensity. The maximum lobe intensity is quickly obtained for the CH 4 -air case with a larger center recirculation than the H 2 -air case. That is, X lob ≈ 4.0 D f for CH 4 -air and X lob ≈ 5.4 D f for H 2 -air. For the H 2 -air case, X lob is shifted downstream compared to that of the CH 4 -air case, but there is little change in the sizes of the wall recirculation. Comparing X lob to the recirculation zones in Figures 5 and 6, the position is roughly consistent with the middle of wall recirculation and the beginning region of center recirculation. In terms of reacting flows, X lob is important because the active reaction can be started at the position. In Figure 10, it is found that X lob is approximately matched to the positions of x ≈ 3.3 D f for CH 4 -air and x ≈ 4.5 D f for H 2 -air where temperatures start to rise sharply. From this result, I lob and X lob are considered as meaningful factors representing combustion characteristics in the multihole baffled micro combustor.
For efficient combustion, it is necessary to obtain a lot of Z = Z st regions inside the combustion chamber. From this sense, the lobed stoichiometric line shows that the mixture is close to a well-mixed state than the circular stoichiometric line in the cross-section of a combustor. As discussed previously, the shape of the lobed flame depends on the number of air holes. The maximized I lob is close to a well-mixed mixture, but it may not represent the best condition for Z = Z st . Therefore, to quantify the Z = Z st region, the equivalent perimeter of the lobed Z = Z st line is introduced as P eq = 2π r eq . The r eq value is obtained by the area integral.
Here R Z st means the radius of the Z = Z st line. After the Z > Z st area surrounded by the Z = Z st line is numerically integrated, the equivalent radius r eq and perimeter P eq are calculated. The high P eq represents a better mixing of reactants regardless of the Z = Z st shape. Figure 14 shows the evolution of the equivalent perimeter and the secondary flows based on the lobe intensity. As the secondary flow develops, the lobed structure is initiated near the baffle plate. The secondary flow becomes stronger as it goes downstream. The resulting Z = Z st line is radially expanded and the lobe intensity is enhanced. For the CH 4 -air case, the lobe intensity is maximized at x/D f = 3.8, i.e. X lob /D f = 3.8, but P eq is maximized at x/D f = 7. The two positions are differently observed. From the baffle shape feature, the central jet dominates the radial expansion of the lobed Z = Z st line. Accordingly, as the lobe intensity or the radial expansion of the Z = Z st line increases, P eq increases. If there is no radial expansion of the Z = Z st line, P eq is maximized at X lob because the lobe intensity is a quantitative expression of the wavy Z = Z st line. However, the peak P eq is obtained after the lobe intensity is the maximum value. Therefore, it is considered that P eq is more largely influenced by the radially outward fuel flow. The maximum of P eq represents a well-mixed state of reactants, and its axial position is the place where the reaction is maximized. As can be seen in Figure 10, the position coincides with the location having the maximum temperature. Also, the location is almost identical to the end of the center recirculation region. After passing through X lob , the flow is further developed in an axial direction, and the lobe intensity is reduced. Finally, the Z = Z st line becomes a circle. From a detailed comparison of Figures 10 and 14, X lob is close to the beginning of the flame zone, and the position of the peak P eq is related to the location after the temperature rises sharply. Therefore, the positions are useful to find the region of high wall heat flux.
In the current baffled micro combustor, the lobed structure can be obtained by changing the baffle shape. One of the possible methods is to change the baffle thickness while fixing other shape conditions. As discussed earlier, the lobed flow structure is related to the variation of radial flows depending on the center and wall recirculations development. From Kim and Park (2019), the vortical structures are differently obtained for the inflow type of a sharp-edged orifice and a long pipe. For the baffled micro combustor, this situation is realized by variable baffle thickness. The jet flow through the baffle hole represents characteristics of an orifice type as the baffle thickness lowers. With increasing baffle thickness, however, it becomes the feature of fully developed pipe flows. By this background, the vorticity magnitude and the lobe intensity are examined for seven different baffle thicknesses. For the CH 4 -air case, Figure 15 shows the magnitudes of maximum vorticities (ω x,max , ω o,max ) at the baffle exit and I lob,max for different baffle thicknesses. The streamwise velocity contours for b t /D f = 0.01 and 1.5 are added in the figure. For b t /D f = 0.01, the velocity profile has off-centered peaks called the saddle-shaped profile. It is because the flow passage is abruptly reduced, and its length is short. At the same time, the velocity profile for b t /D f = 1.5 shows a typical feature of pipe flow. As b t /D f decreases, I lob,max and vorticity magnitudes increase. For b t /D f > 0.5, the values converge to constants. Compared to the thin baffle of b t /D f = 0.01, the vorticity values of b t /D f = 1.5 are 60%-70%. Similarly, I lob,max of b t /D f = 1.5 decreases to 60% for the thin baffle of b t /D f = 0.01. Therefore, it can be understood that the reduction of baffle thickness produces the vortical flows promoting the lobed structure generation, and the realization of the orifice flow from the baffle plate is possible for b t /D f < 0.5. Also, in Figure 15(b), ω x,max and ω o,max have the linear relation for I lob,max . This means that that there is a close relationship between the maximum vorticities (ω x,max , ω o,max ) at the baffle exit and the maximized lobe structure. Generally, when the vortical structures are enhanced in the combustion chamber, the mixing of fuel and oxidizer is improved. From the above similarity, it can be seen that the mixing structure for reacting flows is examined in terms of I lob,max .
To see the mixing performances for different baffle thicknesses, Figure 16(a) shows probability density function (PDF) of the local equivalence ratio φ local for the range of 0 ≤ φ local ≤ 2.0, when φ G = 0.7. φ local is related to the local mixture fraction.
As can be seen, the local equivalence ratio of 0.7 ≤ φ local ≤ 0.75 is frequently observed for b t /D f = 0.01, 0.1, 0.5, and 1.5. It depends on the inlet condition of φ G = 0.7. As b t increases, the PDF value decreases for 0.7 ≤ φ local ≤ 0.75 and more enhanced for φ local = 1.0. For a well-mixed condition, the probability density function becomes the single mode of φ local = φ G . However, since the present baffled combustor is a nonpremixed type, the PDF distributions of φ local depending on the fuel-air mixing. To examine the mixing characteristics, the probability distributions for different baffle thicknesses are divided into three parts: φ local < 0.7, 0.7 ≤ φ local ≤ 0.75, and φ local > 0.75. Figure 16(b) represents the summated probability densities for different I lob,max . As I lob,max decreases, the center parts decrease and the off-center parts are increased. From these results, the high I lob,max case has a good mixing performance, whereas the performance is reduced with decreasing I lob,max and the mixture becomes fuel-rich states.

Energy conversion characteristics
For micropower devices, the energy conversion media of a thermal emitter and n-and p-type semiconductors are installed to use the wall heat transfer. Most such systems require high wall temperature to increase the system performance. In particular, the TPV system uses a conversion process from heat to electricity via photons. If the high wall temperature maintains, the photon emission through the radiative heat transfer of the thermal emitter becomes active. Temperature contours of the baffle plate and the combustion chamber wall, the azimuthally areaaveraged wall temperature, and temperature contours of the x-z plane are plotted in Figure 17 to investigate such characteristics. In the temperature contour plot of the x-z plane, the stoichiometric condition of Z = Z st is represented by the black line.
where A w and T w are the surface area and temperature on the outer combustor wall. The azimuthal A w and T w at the axially same grid cell are used for T avg . In Figure 17(a), the wall temperatures are azimuthally distributed in a wavy form. The high temperature regions are consistent with the number of crests in the lobe structure. As b t increases, the size of the high temperature region and the baffle temperature are reduced by the axially elongated flame. Since the fuel and air of 300 K flow into the combustion chamber, the preheating effect of the baffle plate is confirmed. The preheating effect is more enhanced for the thin baffles. This feature is reflected in the averaged wall temperature of Figure 17(b). That is, as b t increases, the averaged wall temperature is decreased. And, it is large in the region where the lobe structure is dominant. Regardless of the fuel type, the axial variation of the area-averaged wall temperature is almost similar for decreasing b t . But, the temperature values of the H 2 -air case are larger than those of the CH 4 -air case. It is because the H 2 -air case has a longer flame zone. For example, in the case of b t /D f = 0.01, the flame length of the H 2 -air case is increased to 180% of the CH 4 -air case. For a micro combustor, the short flame length is better to get efficient combustion. However, the proper long flame is required to improve energy conversion even if the combustion efficiency is reduced. Further analysis is performed to examine the relationship between the lobe structure and variables representing the reacting flows for the flame length, heat loss, and energy conversion rate. In Figure 18, X lob and the heat loss are plotted for different I lob,max . Here, the heat loss q is estimated by the sum of convective (q convection ) and radiative (q radiation ) heat transfer through the outer combustor wall from inside the combustion chamber to the ambient.
where h w , T ∞ , ε w , and σ are the convection heat transfer coefficient of 20 Wm −2 K −1 , the ambient temperature of 300 K, the wall emissivity of 0.85 (Yang et al., 2014), and the Stephan-Boltzman constant of 5.67 × 10 −8 Wm −2 K −4 , respectively. In the present baffled micro combustor, the increased I lob indicates a wellmixed state of the mixture. Therefore, the flame characteristics and the heat loss caused by convection and radiation are dependent on the change of I lob . Under the given geometry and inlet conditions, X lob and the peak P eq are related to the best mixing state of reactants.
As I lob,max increases, X lob is decreased, but the heat loss increases. This means that the lobed flame structure is quickly obtained for increasing the lobe intensity. When I lob increases, the flame zone is enlarged to the combustor wall, and the wall heat transfer increases due to temperature gradients. At the outer combustor wall, the radiative heat transfer is 5 times larger than the convective heat transfer. For example, the case of b t /D f = 0.01 has the radiation heat rate of 80.1 W for CH 4 -air and 71.8 W for H 2 -air and it has the convection heat rate of 15.3 W for CH 4 -air and 13.7 W for H 2 -air. Compared with the H 2 -air case, the CH 4 -air case has lower X lob and higher heat loss. As discussed in Figure 14, the flame zone is almost initiated at X lob . Therefore, when X lob decreases, the flame length decreases. From the temperature distributions depending on the flame shape, it is comprehended that the shorter flame length gives the wall temperatures of the combustor smaller than the longer flame length. For I lob,max = 1.1, X lob of H 2 -air is 1.5 times larger than that of CH 4 -air. This is consistent with the situation in which the longer flame is obtained for the H 2 -air combustor. To cool the combustor wall, the wall material with the high specific heat capacity may be selected or the heat exchanger like a regenerative cooling may be introduced. However, the increased heat loss can be used for a micro-device such as TPV and TEG. Also, as b t /D f increases, the radiation heat rate on the outer combustor wall is reduced. It represents that the thin baffle is proper for a TPV system. Considering the heat losses and the radiation heat rates for two cases of H 2 -air and CH 4 -air, the CH 4 -air case has the relative advantage for energy conversion.
For the H 2 -air and CH 4 -air cases, Figure 19 shows the variation of the flame length, the combustion efficiency (η), and the emitter efficiency (η e ) for different b t /D f .
Here, two efficiencies are defined as where Y f ,in and Y f ,out are the mass fraction of inlet and outlet (Akhtar et al., 2019;Kim & Park, 2019) and LHV is fuel's lower heating value. For a TPV system, the combustion heat is absorbed into the PV cells in the radiative heat transfer because the photons are promoted by the thermal radiation. On the basis of this, the radiative heat transfer q radiation is used for Equation (14). As I lob,max increases, the lobed shape is enhanced, the flame length decreases, and the combustion efficiency increases. The increased I lob,max can give the better mixing for efficient combustion. Therefore, when I lob,max gets larger, the flame length decreases, and higher η is obtained. As shown in Figure 15, the lobe intensity is enhanced for lower b t cases. From the result, more efficient combustion is obtained for lower b t cases. The flame zone with active reaction rate is described by the distribution of the OH radical. The flame patterns are compared for two CH 4 -air combustors with different baffle thicknesses. The reaction zone of b t /D f = 0.01 is formed than the b t /D f = 1.5 case in the front region. It is related that the thin baffle promotes the vortical flows and the lobe intensity. In Figure 19(b), as the combustion efficiency get betters, the emitter efficiency is improved. The emitter efficiency of the present baffled micro combustor is similar to or better than the micro combustor of Akhtar et al. (2019). For the H 2 -air combustor, when b t /D f < 0.1, the emitter efficiency is η e = 0.319-0.326. When baffle thickness changes, the emitter efficiency of the CH 4 -air combustor is 14.3-25.9% higher than that of the H 2 -air combustor. These results confirm that the thin baffle is favorable for efficient combustions in the baffled micro combustor.

Conclusions
For a multihole baffled micro combustor with a mixed feature of bluff body and multiple jets, the influence of lobe structure on reacting flows is numerically investigated by the detailed chemistry kinetics for CH 4 -air and H 2 -air. From the model tests, it is considered that the combination of RSM and the GRI-Mech 3.0 mechanism is suitable for capturing the multihole effect on reacting flows in the baffled micro combustor. Flow and combustion characteristics for various baffle thicknesses and fuel properties are investigated by the numerical procedure.
Regardless of the fuel type, the recirculation zones and the lobed flame shapes are variously formed for the baffle thicknesses. The lobed structure of several crests and troughs related to air holes is strengthened at the streamwise position near the beginning of center recirculation. The counter-rotating vortex pairs are generated near both sides of crests. They enhance the flows of moving the air to the center of the combustion chamber and the fuel to the combustion chamber wall. Such secondary flows are similar to the flow structure of the lobed jet or nozzle. After the lobe structure is weakened and the reactant mixing is terminated, the maximum temperature is obtained.
To characterize the lobed structure of reacting flows, the lobe intensity I lob and the equivalent perimeter P eq of the lobed Z = Z st line are newly defined. The maximum I lob is observed at the beginning position of the flame zone, its change is related to the wall heat transfer rate. A close relationship between the maximum vorticities at the baffle exit and the maximized lobe structure is also confirmed. When P eq becomes maximum, the temperature rises sharply. Such positions are useful to find the region of high wall heat flux. The enhanced lobe structure is not desirable due to the increased radiative heat transfer, but it is a good feature for a micro-device such as TPV and TEG system.
As the lobe intensity increases, the flame length decreases, and the combustion efficiency is increased. Although the short flame length is better for efficient combustions, it is observed that the longer flame is a good feature for energy conversion. For a micro-TPV system, the radiative heat transfer through the combustor wall is 5 times larger than the convective heat transfer. For b t /D f < 0.1 showing the strong lobe structure, the emitter efficiency of energy conversion is obtained as η e = 0.319-0.326 for the H 2 -air combustor. The efficiency is increased by 26% for the CH 4 -air combustor. The results confirm that the multihole baffled micro combustor can be efficiently utilized for an energy conversion device through the baffle shape change.

Disclosure statement
No potential conflict of interest was reported by the author(s).