Study on Methane Distribution in the Face Zone of the Fully Mechanized Roadway with Overlap Auxiliary Ventilation System

: An overlap auxiliary ventilation system is very often used for driving roadways in methane-rich coal seams. An overlap zone between the outlets of the forcing duct ends with a whirl ﬂow air-duct (WFAD) and the exhaust duct ends with a dust scrubber that is created by applying the overlap system. This study examines the distribution of methane concentrations at various distances in the overlap zone. Maintaining a long overlap zone could increase the advance of the face. Therefore, the impact of overlap zone length on the methane concentration distribution, in and beyond the overlap zone, is investigated. The evaluation of methane concentrations is performed utilizing a well-established computational ﬂuid dynamics (CFD) approach. The mathematical model of methane emissions into the roadway is adopted. Moreover, the CFD model is validated. A vortex of the return air, caused by the free airstream ﬂowing out of the dust scrubber, is found. This air vortex is responsible for higher methane concentrations at the end of the overlap zone. Therefore, the conclusion can be drawn that maintaining the length of the overlap zone at 5 m to 10 m should be done to control permissible methane concentrations.


Introduction
Appropriate working conditions in mining excavations are ensured using throughflow and auxiliary ventilation systems. There is the risk of dangerous methane concentrations during the advance of development headings in a methane-rich coal seam. Therefore, an effective auxiliary ventilation system is vital. The purpose of using an auxiliary ventilation system during roadway drivage in a methane-rich coal seam is to supply the desired amount of air to the face so that a permissible methane concentration is not exceeded along the entire length of the roadway [1][2][3].
There are three auxiliary ventilation systems in the driven roadway: force, exhaust and overlap systems. There are two technical solutions among overlap auxiliary ventilation systems. One of them is called a primary forcing with an exhaust overlap system, and the other is called a primary exhaust with a forcing overlap system [1,[4][5][6][7]. These ventilation systems are selected depending on the technology of roadway drivage and environmental threats, such as methane, coal dust, and thermal conditions, which can result in explosions and/or health issues for miners.
A primary forcing with an exhaust overlap ventilation system is usually used for methane emissions and airborne dust generation from the fully mechanized face of the roadway. Essential components of the system include the primary forcing duct and a short exhaust duct. The advantage of the system is the possibility of dedusting air, taken from the face to the exhaust duct, by including a dust collector. The disadvantages of the system are the zones with various velocities in the cross-section of the roadway and the necessity of constantly moving the exhaust duct, which slows down the advance of the mining roadway. In practice, the overlap zone of the ducts may be lengthened, which may lead to methane layering and the appearance of dangerous concentrations under the roof of the roadway.
For over three decades, numerical modelling has been used to improve safety in driving roadways in methane-rich coal seams. In addition to ensuring mine safety, it is also essential to reduce energy costs for ventilation purposes. Practical solutions for face zone ventilation are mainly sought to remove the dust generated during mechanical mining of rocks in the face of advanced roadways [8,9], methane emitted from rock masses into the roadway [10], or blasting gases while using explosives [11]. One part of the published research concerns developments that are not equipped with machines and other devices. However, there are also research results that take into account the equipment of a roadway face. For example, the influence of a scrubber in the face zone on the distribution of methane concentrations in advancing faces was investigated by, among others, Wala et al. [12,13] and Petrov et al. [14].
Many studies on the auxiliary ventilation of driven roadways in coal mines show dead zones where higher methane concentrations may form [15][16][17]. Cao and Liu [10] pointed out that the mean velocity across a tunnel should be kept at 0.5 m/s to maintain a fully developed tube flow to prevent methane build-up under the tunnel roof.
Kurnia et al. [18] presented some studies on methane dispersion in an underground mine tunnel, assuming a constant methane emission rate and constant methane source velocity. The influence of the number and location of emission sources on the distribution of methane concentrations in the tunnel's cross-section with forced ventilation was demonstrated. However, the impact of the machinery and equipment of the tunnel was not taken into account.
Torano et al. [1] published results of a computational fluid dynamics (CFD) simulation of the methane distribution in a roadway with a forcing ventilation system. The k-epsilon turbulence model was selected and verified using field measurements. They indicated three points of air velocity inversion in the tested excavation and showed the variability of methane concentration along a 60-m section from the face of a roadway. In addition, they demonstrated the possibility of identifying zones with an increased methane concentration. Finally, they suggested the option of using CFD simulations to identify zones with increased methane concentration for reinforced ventilation, such as a forcing system with an exhaust overlap.
Sasmito et al. [19] studied the behavior of methane flow at the blind faces of a roomand-pillar underground coal mine to effectively remove methane. They indicated that the Spallart-Almaras turbulence model is helpful for this type of mine excavation due to the low computational time. They analyzed both exhausting and forcing ventilation systems with the use of brattices.
Hasheminasaba et al. [20], using CFD methods, investigated the effect of cooperation of a ventilation brattice with an exhausting duct system on the distribution of methane concentration in a 12-m section of a dead-end, which did not have machines or operating devices. Rahimi et al. [3] presented methane distribution based on the results of CFD analyses from the same analyzed mine excavation. They used a Monte Carlo simulation to determine gas emissions based on the uncertainty of the gas content. As a result, the optimum ventilation air characteristics were found to minimize high costs and provide a safe environment at the working face.
Numerical simulations of more than one gas distribution in roadways or tunnels are rarely considered [2].
The most popular solution for auxiliary ventilation is a primary forcing system with an exhaust overlap, including a whirl flow-air duct (WFAD) installed at the end of the forcing duct. The WFAD ensures a swirling airflow both towards the face and in the overlap zone of air ducts. To the best of the authors' knowledge, the influence of the length of the overlap zone on the distribution of methane concentration released during working at the face has not been investigated. Airflow simulations in a fully mechanized tunnelling face with Energies 2021, 14, 6379 3 of 23 similar devices were carried out by Hua et al. [6], using the so-called "multi-radial swirling flow generator" and by Hua et al. [7], employing a multi-radial vortex airflow generator. In general, the impact of the location of generators on the effectiveness of dedusting was examined in these works.
The main objective of auxiliary ventilation with an overlap system is to maintain permissible methane concentrations in both the face and overlap zone. From the point of view of work organization at the face, it would be advisable to minimize scrubber shifting frequency and to simultaneously extend the exhaust duct. As mentioned above, such organizational activities increase the length of the ducts' overlap zones and the appearance of potential methane concentrations. Therefore, the question arises, "what length of the zone should be kept to carry out work safely". A method based on CFD modelling was chosen to solve this issue. The model is developed based on the driven roadway in a coal mine. The results of the measurements obtained from that roadway allow to validate the CFD model.

Research Background
The primary forcing with an exhaust overlap system consisted of a forcing duct ending in a WFAD. The secondary exhaust duct was connected to the scrubber and was used to dilute methane and remove dust from the face zone. The inlet to the exhaust duct was located near the face. Air ducts were located partially in parallel, which formed the overlap zone. The WFAD was not located in the overlap zone, and was designed to create swirling air movement towards both the face and the overlap zone.
A duct cassette or storage duct can be placed in front of the WFAD, ending the primary forcing duct. Using the WFAD aimed to remove methane accumulations from under the roof and direct the dust cloud towards the inlet, to the short exhaust duct containing the scrubber. An electrically powered rotor was mounted at the inlet of the whirl-flow air duct. The swirling made the air flow through a 2-cm-wide slot located along the whole whirl-flow air duct. The air velocity was constant along the whole length of the slot. The outlet of the WFAD was closed with a cut-off damper. The opened damper resulted in air outflow towards the face, not through the slot. Such a ventilation option is possible in the case of a non-working face. The WFAD can work on the left and right sides of the roadway without needing to change its structure. The rotor of the WFAD was not a typical fan rotor, and it did not create resistance to the airflow along the forcing duct. When the rotor was switched off, the appropriate structure of the rotor caused it to stop after about 20 s while maintaining the same airflow rate from the forcing duct. Figure 1 shows the principle of operation of such an auxiliary ventilation in a roadway face as well as the WFAD. More detailed information on the WFAD and the dedusting equipment can be obtained from the authors.
In Polish mines, the face zone's required air quantity should be no lower than the mean velocity of 0.3 m/s in the roadway, except for the overlap zone. However, at least 20% of the airflow in the exhaust duct of the scrubber must flow into that zone.
Currently, it is believed the overlap zone should not be longer than 10 m. These recommendations are based on the belief that a more extended zone may increase methane concentrations under the roof due to lower air velocities. However, this concept has not been confirmed on site.
Maintaining a predetermined distance of the overlap zone requires, in turn, a frequency shift of the scrubber together with the advancing face. From the perspective of work organization at the face, it would be advisable to minimize the frequency of shifting and, at the same time, extend the forcing duct ending in the WFAD. If the scrubber is immobilized for a longer time, the distance between the face and the exhaust duct will increase together as the face advances. This can be remedied by adding another exhaust duct; thus maintaining a constant distance between the duct inlet and the face. However, the overlap zone of the systems will then lengthen. The heuristic approaches of researchers have determined the effectiveness of the overlap system with a WFAD so far; however, no one has thoroughly investigated the gas concentration fields of such a ventilation system for an advancing face of a roadway. The question arises as to how the increase in the length of an overlap zone affects the distribution of methane concentrations within it. In response to this, a numerical model was built to simulate the field of both velocity and methane concentrations, considering the dimensions and equipment of a particular roadway. In Polish mines, the face zone's required air quantity should be no lower than the mean velocity of 0.3 m/s in the roadway, except for the overlap zone. However, at least 20% of the airflow in the exhaust duct of the scrubber must flow into that zone.
Currently, it is believed the overlap zone should not be longer than 10 m. These recommendations are based on the belief that a more extended zone may increase methane concentrations under the roof due to lower air velocities. However, this concept has not been confirmed on site.
Maintaining a predetermined distance of the overlap zone requires, in turn, a frequency shift of the scrubber together with the advancing face. From the perspective of work organization at the face, it would be advisable to minimize the frequency of shifting and, at the same time, extend the forcing duct ending in the WFAD. If the scrubber is immobilized for a longer time, the distance between the face and the exhaust duct will increase together as the face advances. This can be remedied by adding another exhaust duct; thus maintaining a constant distance between the duct inlet and the face. However, the overlap zone of the systems will then lengthen.
The heuristic approaches of researchers have determined the effectiveness of the overlap system with a WFAD so far; however, no one has thoroughly investigated the gas concentration fields of such a ventilation system for an advancing face of a roadway. The question arises as to how the increase in the length of an overlap zone affects the distribution of methane concentrations within it. In response to this, a numerical model was built to simulate the field of both velocity and methane concentrations, considering the dimensions and equipment of a particular roadway.

Geometric Model
The analyzed Belt Roadway was driven in one of coal mines, at a depth of 950 m, in the Upper Silesia Coal Basin (USCB), Poland. The roadway was run, among other things, in a 1.0-m thick coal seam 405/1, which was predicted to release 5.0 m 3 /min of methane. The total target length of the roadway driven was 1400 m. The roadway was driven with a 5 per mile slope. The roadway, with an arched section of 17.6 m 2 , was supported with arch support V36/12/4 yielding steel sets. The spacing between the yielding sets was 0.8 m. The cross and longitudinal sections of the 100-m distance to the face are presented in Figure 2. The average advance rate of the roadway driven by an Alpine AM-75 reached 8 m per day. in a 1.0-m thick coal seam 405/1, which was predicted to release 5.0 m 3 /min of methane. The total target length of the roadway driven was 1400 m. The roadway was driven with a 5 per mile slope. The roadway, with an arched section of 17.6 m 2 , was supported with arch support V36/12/4 yielding steel sets. The spacing between the yielding sets was 0.8 m. The cross and longitudinal sections of the 100-m distance to the face are presented in Figure 2. The average advance rate of the roadway driven by an Alpine AM-75 reached 8 m per day. Analysis of methane distribution in the roadway required a precise determination of the devices used; namely, the ventilation devices, mining machines, and conveyors. The following devices were used in this case:  road header machine, AM75;  whirl-flow air duct, WIR-1000 (WFAD);  dust collector, UO 800-1000;  belt conveyor, GWAREK-1000;  chain conveyor, SKAT 80/KJ.
The equipment was located centrally in the cross-section of the roadway. Furthermore, ventilation devices were located at the roadway walls, following the rules of the overlap ventilation system. The scrubber was suspended 6 m from the face, and its axis was 2 m above the floor. The WFAD was located on the floor at a distance of 12 m from the face.

Analyzed Cases
In practice, it is easier to lengthen both the primary forcing duct and the exhaust duct during roadway driving than to move the scrubber behind the advancing face of the roadway. However, this results in lengthening the overlap zone, as such, there is a possibility Analysis of methane distribution in the roadway required a precise determination of the devices used; namely, the ventilation devices, mining machines, and conveyors. The following devices were used in this case: • road header machine, AM75; • whirl-flow air duct, WIR-1000 (WFAD); • dust collector, UO 800-1000; • belt conveyor, GWAREK-1000; • chain conveyor, SKAT 80/KJ.
The equipment was located centrally in the cross-section of the roadway. Furthermore, ventilation devices were located at the roadway walls, following the rules of the overlap ventilation system. The scrubber was suspended 6 m from the face, and its axis was 2 m above the floor. The WFAD was located on the floor at a distance of 12 m from the face.

Analyzed Cases
In practice, it is easier to lengthen both the primary forcing duct and the exhaust duct during roadway driving than to move the scrubber behind the advancing face of the roadway. However, this results in lengthening the overlap zone, as such, there is a possibility of increased methane concentrations along the overlap zone. Therefore, four cases of different lengths of overlap zones were considered: 5, 10, 15, and 20 m, respectively. The assumptions for the CFD calculations are presented in Table 1.  The model takes into consideration well-established sources of steady and homogeneous methane emissions with the assumption of the continuous roadheader operation.

Model Development
Preparation of the computational model required making a fluid domain based on the physical model of the roadway. In this model, considering one of the flow conditions, two fluid domains were made without voids, but sharing the same topology. This allowed to improve the quality of the computational mesh, limit its elements, and thus reduce the computation time. The first fluid object was extracted by cutting out a fragment of the devices in each roadway cross-section. Then, the second fluid object was inserted as a separate object into the place of the scrubber cut-out.
The numerical model, based on computational fluid mechanics, was performed using ANSYS Workbench version 19.3. The geometric model of the roadway was developed in SolidWorks 2018 and was adapted for simulation in Ansys Space Claim. The actual dimensions of the part of the roadway were applied to build the spatial geometry. The analysis consists of four numerical models varying in exhaust duct length, ending in the scrubber. Figure 3 presents the location of the equipment in the face zone of the model. of increased methane concentrations along the overlap zone. Therefore, four cases of different lengths of overlap zones were considered: 5, 10, 15, and 20 m, respectively. The assumptions for the CFD calculations are presented in Table 1.  The model takes into consideration well-established sources of steady and homogeneous methane emissions with the assumption of the continuous roadheader operation.

Model Development
Preparation of the computational model required making a fluid domain based on the physical model of the roadway. In this model, considering one of the flow conditions, two fluid domains were made without voids, but sharing the same topology. This allowed to improve the quality of the computational mesh, limit its elements, and thus reduce the computation time. The first fluid object was extracted by cutting out a fragment of the devices in each roadway cross-section. Then, the second fluid object was inserted as a separate object into the place of the scrubber cut-out.
The numerical model, based on computational fluid mechanics, was performed using ANSYS Workbench version 19.3. The geometric model of the roadway was developed in SolidWorks 2018 and was adapted for simulation in Ansys Space Claim. The actual dimensions of the part of the roadway were applied to build the spatial geometry. The analysis consists of four numerical models varying in exhaust duct length, ending in the scrubber. Figure 3 presents the location of the equipment in the face zone of the model.

Discretization of the Geometrical Model
The same meshing software was used in simulations for all cases with the same settings. Fluent Mesher was applied to create the computational mesh instead of the Ansys Mesher Module; it was used to get a better-quality mesh.
The convergence test for the computing mesh was conducted for a few cases. Total computing time (from 22 h to 112 h), net balance of the flow (from −0.0000314 kg·s −1 to −0.0000076 kg·s −1 ), number of elements (from 6.5 million to 29 million), single iteration computing time, and the quality of the mesh (from 0.42 to 0.43), were checked in the test, respectively.
Three different meshes were made in this analysis for one case (that of the 5 m overlapping zone). It was assumed that the mesh size increased twice between tested cases. Three models with meshes of 6.5 million elements, 14 million elements, and 29 million elements were obtained. The net flow balance in the computing models, between the inlets and outlets, was conducted for the airflow without methane.
The Watertight scenario discretized the prepared geometry. All four computational models were digitized using the same mesher, with the same settings for each case. The size of the mesh for objects in the geometry was determined in the first stage, as follows: When the local sizing mesh was completed, the mesh size was varied from 0.0035 m to 0.35 m. The surface mesh was made in the whole numerical domain using settings presented in Table 2. However, the generated surface mesh was of low quality, which influenced the quality of the volumetric mesh; therefore, its quality was improved using algorithms upgrading the surface mesh, according to the following settings: • quality improvement of max angle: 0.85; • quality improvement of iterations: 500; • quality improvement of collapse skewness limit: 0.85.
As a result, surface meshes with a maximum skewness of 0.62 ± 0.2 for all cases were obtained. In the following step, the geometries were determined, i.e., their physical features, boundary conditions, external interfaces, and fluid regions. Then, volume discretization of the computational domain was performed usig the parameters presented in Table 3. The obtained computational mesh was still of low quality; therefore, the re-algorithms were used to improve the quality of the computational mesh, this time for the volumetric mesh, with the following settings: The view of the mesh of an exemplary cross-section is shown in Figure 4. Table 4 shows the obtained results for the parameters of the volumetric computational grids.

Sources of Methane Emission
With the advancing faces of mine roadways, methane is emitted from many sources, including exposed planes of the coal seam and mined coal. The total volumetric emission of methane, Q t , to the roadway was the sum of methane emissions from coal mined, Q c , from the coal surface in roadway walls, Q w , and from the exposed coal face, Q f .

Methane Release from Freshly Mined Coal Transported by a Conveyor
Determined in the laboratory, the methane content of coal, V i in the developed coal seam 405/1 was 7.7 m 3 /t daf . The initial methane content, V i is the sum of desorbable V d and residual V r gas contents. Desorbable gas content was 6.2 m 3 /t daf . The change in the desorbable gas content of coal with time, τ, can be expressed by a differential equation: where V(τ) is the methane content of desorbed methane remaining in one ton of mined coal after time τ (instantaneous methane content of coal), m 3 /t, V d is the initial value of desorbable gas content. The constant characterizing the rate of methane desorption from the mined coal was determined by c d , 1/min. The constant c d depends on the properties of the carbon-methane phase system, i.e., mainly the fragmentation of coal and the size of the developed area, temperature, and pressure of the adsorbed gas, τ is the time counted from the beginning of methane separation from the mined coal, min. The solution to Equation (1) was given in [21,22]. Assuming the initial condition in the form: τ = 0 and V(τ) = V d , the solution of Equation (1) is obtained as follows: Assuming the continuous operation of both the roadheader and the conveyor, amount of methane Q c , m 3 /min, which is released from the mined coal for time τ can be presented as follows: where v t is the transport speed of the output, m/min, x is the coordinate of the length measured along the axis of the roadway from the coal loading site, m and A is the continuous coal output from the face, t/min. Constant c d ranges from 1.4 to 1.6 1/min for the Upper Silesian Coal Basin [23].

Methane Release from Coal Surfaces along Walls of Roadway
The rate of gas emission change from the exposed coal surface of roadway walls, q w m 3 /(m 2 min), can be determined as follows [24]: where Q o is the initial methane emission at the time of exposure of the coal surface, m 3 /(m 2 min), k is the constant characterizing the rate of methane emission from exposed coal surface, 1/min, which depends on the rate of methane desorption from coal but also the permeability of the fracture zone around the driven roadway [25], τ is the time that has passed since the exposure of the coal surface, min. Constant k assumes values in the range of 0.0001-0.0002 1/min [21]. By assuming a constant face advance, the following expression is used to determine the volumetric flow rate of methane from two surfaces of coal, located on the left and right of the roadway [24]: where m is the thickness of the exposed coal seam, m, v d is the advance of the roadway face, m/min. Methane emission at the moment of exposing the coal surface, Q o , is the function of initial methane content, V i [26], and is given as follows:

Methane Release from Exposed Coal Surface in Roadway Face
Methane emission from the surface of the roadway face, Q f at the time of its exposure can be calculated according to formula: where b is the width of the roadway, m. Methane emission from the face is a relatively small value. Table 5 shows the data results for determining methane emission sources that were included in the numerical model. The methane emission model was implemented in the CFD simulation model. Several conditions for maintaining the flow and dispersion of methane for four cases of the length of the overlap zones were investigated to obtain an answer to the influence of increasing the overlap distance on the distribution of methane concentrations and safety of mining works.

Boundary Conditions and Solver
A Fluent solver was used to perform numerical calculations. The simulation was performed as a stationary one, taking gravity into account. The following data were selected in the settings of the computational model: • k-e Realizable flow model with enhancement wall treatment function; • model species for the methane-air mixture (components: methane, air) with the density model set to volume-weighted mixing law; • energy equation was excluded due to the change in the method of calculating the density of gas mixture as the isothermal condition was assumed.
The model was prepared in the form of two volumes sharing the same topology. This assumption permitted for using the following conditions for cell zones: • For the first zone-the roadway-two sources were added. The first one relates to methane emission to the roadway, while the second one results from the emission source, i.e., the source of mass. The model of methane emission presented in Section 3.3 was implemented using the User Defined Function, prepared especially for this case. Its task is to ensure the optimal emission along the whole length of the excavation from selected surfaces. The UDF overlays the first discretization layer because there is only methane in each iteration (boundary condition). This function is used for methane emission from the walls, the face, and the upper surfaces of the belt and scraper conveyors. The total value for the selected areas is expressed by: where x is the distance of the unit surface from the face of the roadway and depends on the distance of a given mesh from the face and is proportional to the size of the mesh itself.
• For the second zone-the dust scrubber-the source of momentum towards axis Z of value 9 N/m 3 is related to the flow rate in the scrubber of value 6.64 m 3 /s ± 0.03 m 3 /s.
The boundary conditions were defined as a mass-flow-inlet for the WFAD of 8 kg/s in the normal direction to the slot of the WFAD. The pressure-outlet condition for this slot was left as the default value.

Governing Equations
The k-ε model of turbulence for the mixing law of the mixture of methane-air without additional factors was used for considering the roadway model. Many factors affect the motion of methane, and airflow's carrying effect is the most significant. The simulation results of the airflow field determine the accuracy of the methane motion calculation. The mathematical model of the airflow movement in a fully-mechanized roadway face was accurately established using ANSYS Fluent software [27]. The equation for the conservation of mass, or continuity equation, can be expressed as: Equation (9) is the most common form of the mass conservation equation, and it is applicable to compressible and incompressible fluids. Source S m is the mass added to the continuous phase from other dispersed phases or any user-defined source (UDS).
The momentum equation of airflow can be expressed as follows: where p is the static pressure, τ is the stress tensor, ρ Stress tensor τ can be expressed as follows: where µ is the molecular viscosity, l is the unit tensor and the second term on the right-hand side of the formula is volume dilatation. The turbulence kinetic energy and rate of dissipation are obtained as follows: where G k represents the generation of turbulence kinetic energy due to mean velocity gradients, G b is the generation of turbulence kinetic energy due to buoyancy, Y M represents the contribution of the fluctuating dilatation incompressible turbulence to the overall dissipation rate, C 2 and C 3ε are constants, σ k and σ ε are the turbulent Prandtl numbers for k and ε, respectively, S k and S ε , are user-defined source terms.

Validation of Turbulence Model
The filed measurements were made during the Belt Roadway drilling. The length of the overlap zone was kept at 10 m. The measurements were carried out during the first working shift from 6:00 a.m. to 2:00 p.m. They included manual methane concentration measurements using a Drager X-AM 7000 methane meter with 2.5% accuracy of the measured value. Measurements were taken at selected points of the roadway cross-section at the height of 1.0 m and 2.0 m above the roadway floor. Regardless of the manual measurements, the methane concentration data from sensors of the real-time monitoring system in the roadway were obtained. The accuracy of methane monitor sensors was also 0.1% vol. Three typical cross-sections of the roadway with built-in methane sensors were selected to validate the model. The methane sensors were located at a distance of 0.1 m from the roof of the roadway. Selected sensors in the cross-sections are marked in Figure 2. The measurement points in the cross-sections are shown in Figures 5-7, which present the distribution of methane concentrations obtained from the CFD simulation. Manual measurements were performed three times during roadheader working. Results presented in Figures 5-7 have average values obtained from the measurements mentioned above. Manual measurements were performed three times during roadheader working. Results presented in Figures 5-7 have average values obtained from the measurements mentioned above.   Manual measurements were performed three times during roadheader working. Results presented in Figures 5-7 have average values obtained from the measurements mentioned above.   The airflow rate at the exit of the roadway was also measured. The absolute methane emission into the roadway was calculated based on the determined airflow rate and indications of methane concentration of the M6 sensor. The methane emission ranged from 4.85 m 3 /min to 5.15 m 3 /min with an average value of 5.00 m 3 /min.
The comparison of methane concentrations is presented in Figures 5-7. The differ- The airflow rate at the exit of the roadway was also measured. The absolute methane emission into the roadway was calculated based on the determined airflow rate and indications of methane concentration of the M6 sensor. The methane emission ranged from 4.85 m 3 /min to 5.15 m 3 /min with an average value of 5.00 m 3 /min.
The comparison of methane concentrations is presented in Figures 5-7. The difference between the results of measurements and calculations does not generally exceed 10%; therefore, it is assumed that the validation of the adopted turbulence model is approved.
As the behavior of overall airflow and methane concentration for design purposes is investigated, model k-ε realizable with functional enhancement wall treatment was considered to be sufficient.

Results and Discussion
The boundary conditions are the same for all four analyzed cases. A total of 80% of the fresh air flowing out of the WFAD moved toward the face, ensuring methane dilution and proper transport of return air to the scrubber. A total of 20% of the fresh air flowed through the overlap zone in a swirling motion.
The only difference that affects the flow of fluids in the roadway is the overlap zone's length. The simulation results indicate that some of the air leaving the scrubber returns along the opposite scrubber's sidewall to the overlap zone. This effect of the vortex crating air backflow has not been detected so far, although previously Li et al. [28] described a similar backflow, but in the face zone when using the forcing ventilation system.
The air vortex is visible at the roadway's left sidewall in Figure 8a On the other hand, for = 15 m and = 20 m cases, a much shorter range of air backflow to the overlap zone is visible than for the cases = 5 m and = 10 m. This is because the collision of the backflow vortex with fresh air occurs at the end of the overlap zone.
The excavation equipment also influences the range of air backflow effects, which is confirmed by the results of other studies [29,30]. Therefore, it is essential to consider equipment in the roadway during numerical modelling, especially with such a complex The excavation equipment also influences the range of air backflow effects, which is confirmed by the results of other studies [29,30]. Therefore, it is essential to consider equipment in the roadway during numerical modelling, especially with such a complex ventilation as in overlap auxiliary ventilation systems.
The observed phenomena of airflow in such a ventilation system significantly affect the distribution of methane concentration. Figure 9a-d shows that methane concentrations are highest in the face zone, which results from the location and intensity of methane emission sources. A decrease in gas concentration is noticeable in the overlap zone compared to the one in the face zone. In the case of L o = 5 m, the vortex of air flowing along the left sidewall achieved the WFAD disturb the swirl profile of air flowing out (Figures 8a and 9a). The alignment of the methane concentration concerning the behind-face zone is visible at the end of the overlap zone. The length of the overlap zone influences the profile of the backflow vortex when mixed with the swirling flow of fresh air. As the overlap zone's length increases, the air backflow has more space to spread (Figure 10b-d). The alignment of the methane concentration concerning the behind-face zone is visible at the end of the overlap zone. The length of the overlap zone influences the profile of the backflow vortex when mixed with the swirling flow of fresh air. As the overlap zone's length increases, the air backflow has more space to spread (Figure 10b-d). The alignment of the methane concentration concerning the behind-face zone is visible at the end of the overlap zone. The length of the overlap zone influences the profile of the backflow vortex when mixed with the swirling flow of fresh air. As the overlap zone's length increases, the air backflow has more space to spread (Figure 10b-d).  The alignment of the methane concentration concerning the behind-face zone is visible at the end of the overlap zone. The length of the overlap zone influences the profile of the backflow vortex when mixed with the swirling flow of fresh air. As the overlap zone's length increases, the air backflow has more space to spread (Figure 10b-d). For the case of L o = 5 m, the vortex flows, even up to 24 m from the face, i.e., 2 m from the beginning of the overlap zone (Figure 10a). The uniform methane concentration in that cross-section is visible (Figure 11a), with a value close to the average methane concentration in the other part of the roadway. This results from the inflow of vortex of return air with an increased methane content. Figure 11b shows the velocity profile at the distance of 10 m from the scrubber outlet. The increase in velocity is visible on the left side of the roadway, delivering air into the overlap zone. For the case of = 5 m, the vortex flows, even up to 24 m from the face, i.e., 2 m from the beginning of the overlap zone (Figure 10a). The uniform methane concentration in that cross-section is visible (Figure 11a), with a value close to the average methane concentration in the other part of the roadway. This results from the inflow of vortex of return air with an increased methane content. Figure 11b shows the velocity profile at the distance of 10 m from the scrubber outlet. The increase in velocity is visible on the left side of the roadway, delivering air into the overlap zone.   For the case of = 10 m, the effect of mixing the vortex with the swirling flow of fresh air appears at the distance of 29 m to the face, i.e., precisely 5 m further than in the case of = 5 m. Moreover, in this case, the methane concentration was leveled in the cross-section of the roadway (Figure 12a). The effect of turbulence on the airflow velocity is visible at the distance of 38 m to the face, i.e., 11 m from the beginning of the overlap zone (Figure 12b).      To sum up, it should be stated that, until now, a belief has been held that methane concentrations in an overlap zone result from only a small amount of air flowing through this zone with a swirling motion. The results of numerical simulations have shown that the air-methane mixture flowing out of the scrubber returns to the overlap zone. Therefore, there might be an additional reason for an increase in the methane concentration in that zone. If the length of the overlap zone is short (in the case of = 5 m), the velocity profile of the swirling air, leaving the WFAD, may be disturbed. It may disturb the swirl of fresh air and increase methane concentration under the roof of the overlap zone. In turn, in cases of = 15 m and = 20 m, the backflow zone in the behind-face zone is lengthened, where methane concentration may increase compared to the case of = 10 m. However, in the behind-face zone, methane concentration is the lowest in all the cases. In the case of increasing the overlap zone's length, the distance in which the concentration To sum up, it should be stated that, until now, a belief has been held that methane concentrations in an overlap zone result from only a small amount of air flowing through this zone with a swirling motion. The results of numerical simulations have shown that the air-methane mixture flowing out of the scrubber returns to the overlap zone. Therefore, there might be an additional reason for an increase in the methane concentration in that zone. If the length of the overlap zone is short (in the case of L o = 5 m), the velocity profile of the swirling air, leaving the WFAD, may be disturbed. It may disturb the swirl of fresh air and increase methane concentration under the roof of the overlap zone. In turn, in cases of L o = 15 m and L o = 20 m, the backflow zone in the behind-face zone is lengthened, where methane concentration may increase compared to the case of L o = 10 m. However, in the behind-face zone, methane concentration is the lowest in all the cases. In the case of increasing the overlap zone's length, the distance in which the concentration equalization occurs moves further from the face (Figure 15). The advance rate of the heading face is practically in proportion to the increase in the length of the overlap zone. The mean methane concentrations in the analyzed cross-sections indicate that extending the overlap zone would be beneficial because it would increase the region with the reduced methane concentration. However, methane concentrations in the behind-face zone are higher.  Therefore, the methane concentration was analyzed at 11 characteristic points of the cross-section within the length of the roadway (Figure 16). The locations of the analyzed places in the cross-section are shown in Figure 16. Results show the lowest concentration of methane for case = 10 m at most points of the cross-section; therefore, it should be  Therefore, the methane concentration was analyzed at 11 characteristic points of the cross-section within the length of the roadway (Figure 16). The locations of the analyzed places in the cross-section are shown in Figure 16. Results show the lowest concentration of methane for case L o = 10 m at most points of the cross-section; therefore, it should be assumed that the best case is L o = 10 m. Figure 17a,k shows methane concentration at measurement points within the roadway.   The authors realize that they obtained the results for specific dimensions of the excavation and its equipment. Nevertheless, the detected backflows in the behind-face zone will always limit the swirling flow of intake air in the overlap zone. The qualitative effect of this phenomenon may be different in other roadways. However, the computational fluid dynamics (CFD) approach can be used to study this phenomenon. However, there is no doubt that a minor impact of return air backflows will occur if the length of the ducts' overlap zones range from 5 to 10 m. The authors realize that they obtained the results for specific dimensions of the excavation and its equipment. Nevertheless, the detected backflows in the behind-face zone will always limit the swirling flow of intake air in the overlap zone. The qualitative effect of this phenomenon may be different in other roadways. However, the computational fluid dynamics (CFD) approach can be used to study this phenomenon. However, there is no doubt that a minor impact of return air backflows will occur if the length of the ducts' overlap zones range from 5 to 10 m.

Conclusions
The selection of an auxiliary ventilation system for development workings predicted for the conditions of methane emission should consider the methane distribution along the advancing roadway equipped with ventilation and technological devices. The release of methane and its distribution with the air in a run-of-mine roadway was simulated using the computational fluid dynamics (CFD) approach for primary forcing with an exhaust overlap system. Investigation into the length of the overlap zone between the force and exhaust ducts under the conditions of methane emission were shown: 1. When considering the distribution of methane concentrations in roadways with auxiliary ventilation, it is essential to properly reflect the geometry and equipment of the roadway and the number, location, and efficiency of methane emission sources. The roadway ventilation model, validated by measurements, confirmed the validity of the assumed functions of changes in methane emission from three sources: coal surface in the face, coal surfaces in the sidewalls of the roadway, and mined coal transported by a conveyor.
2. The whirl-flow air duct (WFAD) at the end of the force duct ensures swirling airflow both to the face and in the overlap zone. The supply of 20% more air through the forcing duct than the air discharged through a short exhaust duct ensures effective mixing of gases with a swirling motion and the removal of methane concentration from under the roof of the roadway.
3. The length of the overlap zone is essential to the distribution of methane concentration in the roadway. The simulation cases carried out for all four lengths of 5 m, 10 m, 15 m and 20 m of the overlap zone showed that methane concentrations in that zone are effectively lowered. However, the results obtained have also shown an air backflow along the sidewall opposite the exhaust overlap duct. It is due to the formation of a negative pressure zone around the free stream, leaving the exhaust overlap duct, in which the air backflow phenomenon occurs.
4. If the length of the overlap zone is only 5 m, the return air backflow may disturb the stabilized velocity profile of the intake air, leaving the WFAD. If the overlap zone is 15 m or 20 m, the range of the air backflow is greater and results in higher methane concentrations in the region from the scrubber outlet to the point of free stream die down.
5. The test results have shown that an overlap distance ranging from 5 to 10 m should be maintained between the force and the exhaust ducts for methane control in a face zone. The roadway driving organization should consider the necessity of simultaneous reconstruction of a forcing duct ended with the WFAD and an exhaust duct equipped with a scrubber.