Airflow Velocity Research Articles

Simulation and optimization on a novel air-cooled photovoltaic/thermal collector based on micro heat pipe arrays

To further enhance the performance of an air-cooled photovoltaic/thermal (PV/T) system, this study integrates micro heat pipe arrays (MHPAs) and serrated fins as the primary heat dissipation components. The lumped-parameter method was utilized to establish a mathematical model for the previously proposed PV/T system based on MHPAs. The model reliability was validated with the test data (errors within ± 10 %). Subsequently, the influence of solar irradiance, ambient temperature, airflow velocity inside the duct, and component structure on the air-cooled PV/T system performance was analyzed and optimized using the method of controlling variables. The impact of varying the condensing section length of the MHPA from 4 to 28 cm on the system performance was analyzed, and it was found that the condensing section length had a relatively small effect on the system thermoelectric performance. When the optimal comprehensive performance was achieved with fin heights ranging from 6 to 18 mm, heat collecting efficiency remained stable at around 19.5 %. When the cross-sectional dimension of the duct was 25 cm, the power generation efficiency reached its maximum value of 12.59 %. Such research findings provide data support for the promotion, development, and wider applications of air-cooled PV/T systems in energy utilization.

Numerical Simulation and Analysis of the Airflow Field in the Crushing Chamber of the Hammer Mill.

The airflow dynamics within hammer mills' crushing chambers significantly affect material crushing and screening. Understanding the crushing mechanism necessitates studying the airflow distribution. Using a self-built crushing test platform and computational fluid dynamics (CFD) simulations, we investigated the impact of screen aperture size, rotor speed, hammer-screen clearance, hammer quantity, and mass flow rate on airflow distribution within the rotor region, circulation layer, and screen apertures. Results indicated generally uniform axial static pressure distribution within the rotor region, with radial gradients. Increased rotor speed improved radial static pressure gradients, while higher mass flow rates reduced them. The highest airflow velocity within the circulation layer reached approximately 83.46% of the hammer tip's tangential velocity. Greater rotor speed and hammer quantity intensified circulation airflow, whereas increased mass flow rate decreased it. Eddies formed within screen apertures with higher rotor speeds and hammer quantities but diminished with larger apertures and higher mass flow rates. Static pressure differences across screen apertures increased with mass flow rate and rotor speed but decreased significantly with larger apertures. This systematic examination provides insights into airflow distribution within hammer mill crushing chambers, offering a theoretical foundation for improving and designing hammer mills.

Pneumatic separation and process simulation of waste printed circuit boards pyrolysis residue

In this study, the application of Z-shape pneumatic separator for separating metal and nonmetal particles in pyrolysis residue of waste printed circuit boards was investigated, and the flow field characteristics and particle motion behavior in the separator were revealed through flow field simulation. The pneumatic separation results indicated that pre-screening significantly enhanced the separation efficiency of metals and nonmetals in pyrolysis residue. The simulation results showed that the large gradient of airflow velocity distribution and the turbulent flow field within the Z-shape separation cavity are beneficial for dispersing particles and reducing the influence between particles of different properties, thus enhancing separation efficiency. Using Rocky DEM for modeling of different shapes and properties of particles, and then combined with ANSYS Fluent can better simulate the particle separation behavior in the pneumatic separator.

Application of CFD Simulation to the Design of an Innovative Drying Chamber

Drying and sanitising equipment has been very common in industrial plants since the pandemic. These are devices that consume significant amounts of energy. The best solution is to use a drying chamber equipped with a heat pump, which allows partial recovery of the energy. In the design of the drying chamber, the drying time is important, which depends both on the parameters of the heat pump itself, and the geometry and airflow of the drying chamber. The geometry and airflow supply should be arranged to ensure a uniform distribution of velocity throughout the drying area. For this purpose, the use of CFD simulations has been proposed. A model was developed in ANSYS/FLUENT where all model parameters, including the optimal mesh density, turbulence models, etc., were determined. The model was verified on the experimental results of the basic design of the chamber. Then an innovative design was proposed that was modelled and optimised in terms of the distribution of the inlet’s perforation. The final design was made, and, at the same time, the simulation’s results were verified by measuring the velocity of airflow in the new design. Together with the optimisation of the heat pump, this made it possible to reduce the drying time by 50%, with a simultaneous reduction in the energy consumed.

Stack-level analysis of the performance variation in air-cooled PEMFC with Z-type anode manifold

Air-cooled proton exchange membrane fuel cell (PEMFC) stack holds significant promise due to its potential to reduce system complexity. The anode manifold type leads to gas maldistribution in an air-cooled stack, and affects the efficiency and life span of the stack. In this study, 3D multi-physics model of stacks ranging from 10 to 40 cells are established to investigate the gas distribution characteristic of Z-type anode manifold. Considering the heat dissipation aspects only, it is advisable to keep a thin endplate while meeting the required stress tolerances. The core temperature of the stack with 20 mm endplates measures 43.8 °C, while the stack with 6 mm endplates reaches a core temperature of 45.5 °C. Under a cathode air flow velocity condition of 4 m‧s−1, the maximum net oxygen consumption reaches 1.69 mol‧s−1 in Cell-4. Furthermore, the electrochemical reaction intensity decreases sequentially from Cell-4 and Cell-5 in the middle to Cell-1 and Cell-10 in the edges of the stack. Additionally, cells near the anode inlet and outlet are prone to have the insufficient gas supply in the Z-type anode manifold stack.

Numerical study on the gas-solid diffusion characteristics of airflow and dust particles in metal mine blasting excavation under three hybrid ventilation modes

Dust pollution has long been a challenge in the blasting excavation of underground metal mines, and uncontrolled diffusion of dust particles seriously endangers the health and safety of workers as well as the environment in underground space. Most existing research on dust diffusion focuses on fully mechanized coal mining face or borehole drilling, where the dust source emission characteristics differ significantly from metal mine blasting excavation. Furthermore, it is still unclear how the ventilation modes of Far-Pressing-Far-Absorption (FPFA), Near-Pressing-Far-Absorption (NPFA), and Far-Pressing-Near-Absorption (FPNA) affect the mechanism of dust diffusion in the breathing zone. In this work, gas-solid flow characteristics in a typical metal mine blasting tunnel are numerically studied based on Euler-Lagrange method. The interphase forces between airflow and dust particles are comprehensively modeled, and the particle diffusion effect caused by fluid turbulence is described by the discrete random walk model. Five typical regions for airflow disorder phenomena in the working tunnel are predicted, including the jet development region, the return airflow core region, the airflow reversal region, the airflow recirculation region, and the secondary flow region. The spatial distribution of dust exhibits nonuniform characteristics under turbulent airflow transport. The time results of the monitored dust concentration reveal that the junction of the working tunnel and the return airway is a key area of focus for the prevention and control of dust pollution. Among the three ventilation modes, the Far-Pressing-Near-Absorption (FPNA) ventilation mode has the best dust removal efficiency. The supply vent further from the blasting face allows the jet to fully develop, while the exhaust vent closer to the blasting face allows dust particles to be drawn in and discharged earlier. Under the condition of fixed total power consumption for the supply and exhaust fans, a larger exhaust airflow velocity helps remove dust from the breathing zone.

Motion behaviors of droplets containing Au nanoparticles on a superhydrophobic laser-induced graphene surface

The movement of nanoparticle-containing droplets on solid surfaces significantly affects the distribution of the nanoparticles and is of great interest in the fields of two-phase separation, biosensing detection, inkjet printing, and microarrays. There has been little research on the initiation and motion behaviors of colloidal droplets containing nanoparticles on superhydrophobic surfaces. Here, we prepare superhydrophobic laser-induced graphene (LIG) surfaces with excellent depinning effects using an extremely simple method and explore the formation mechanism of the depinning-LIG surfaces. The reduction of nano-graphene fibers and the increased hydroxyl group ratio after alcohol modification further enhance the hydrophobic properties of depinning-LIG, reducing its surface adhesion. The initial and continuous motion of droplets containing Au nanoparticles (AuNPs) on these superhydrophobic surfaces under airflow is studied using high-speed microscopy. The coupling effects of the droplet size, surface properties, airflow velocity, and nanoparticles on the droplet motion behaviors are analyzed. The dimensionless parameter G is incorporated to obtain the partition diagram of AuNP droplet motion behaviors on depinning-LIG surfaces, which delineate the critical conditions for droplet “oscillation,” “initiate sliding,” and “continuous rolling” as a function of system parameters. For AuNP droplets, the viscous force Fγ,p exerted by the nanoparticles on the contact line significantly affects the droplet movement behaviors. In addition, a mathematical model about the competition of dynamic forces and resistance is established to describe the motion of AuNP droplets, and the critical conditions for different motion behaviors of the droplet are clarified to guide practical applications.

Experimental investigation of interactions between a water droplet and an airflow boundary layer

The deformation and movement of droplets is widely relevant in many fields of research. The present work experimentally investigates the evolution of a single droplet interacting with an air boundary layer. A series of experiments are carried out using a high-speed photography technique to determine the effects of the airflow velocity, drop height, and droplet size. The morphological characteristics can be classified into three types according to the experiments. The outcomes indicate that both the drop height and the airflow velocity significantly influence the maximum streamwise spreading length, but only the drop height has an impact on the maximum lateral spreading width. The maximum streamwise spreading factor follows a power function relationship with WeRe−0.5. In addition, the crater maximum streamwise and lateral spreading diameters are mainly influenced by the drop height. An energy conversion model is established by considering the effects of the aerodynamic drag force, surface tension, and viscous force. This study provides experimental reference data for the scenario of a droplet interacting with an air boundary layer.

Heat transfer simulation for re-design of tray dryer to reduce the energy consumption in the cocoa bean drying process

Abstract The tray dryer for the cocoa bean drying process is a solution to overcome conventional drying methods that rely on solar heat. The tray dryer is expected to facilitate the cocoa farmers to dry cocoa beans because it is faster and does not depend on weather conditions. However, even though the current condition of the tray dryer can reduce the energy consumption of drying cocoa beans (from Gunung Kidul, Yogyakarta, Indonesia) compared to the solar heat method, it remains unknown whether the cocoa beans drying quality is good enough. Therefore, this research analyzed the heat transfer phenomena of the cocoa drying process using the computational simulation method to re-design the current tray dryer design so that the new design can fit the drying characteristic of the cocoa bean and reduce energy consumption. As a result, for the air velocity of 25.91 m/s and the ½ opening of the gas valve in 10 minutes, which is already at steady-state conditions, the temperature distribution in the existing tray dryer is not suitable with the characteristics of the cocoa beans drying process. The temperature can be increased by reducing the airflow velocity or increasing the heating value by opening the valve. However, with a new design, the temperature can be increased by up to 37.00% at a velocity of 25.91 m/s with a ½ gas valve opening. Furthermore, the new design can achieve the same conditions as the existing design by only using an opening for a 1/8 gas valve. Thus, using the new design, the tray dryer can achieve the drying characteristics of cocoa beans and reduce energy consumption.

Microscopic Mechanism and Percolation Model of Dynamic Deposition of Elemental Sulfur Particles in Acidic Gas Reservoirs.

The dissolution of elemental sulfur in acidic gas leads to its precipitation as gas pressure decreases, thereby causing potential damage to the formation due to the deposition of sulfur particles. Previous sulfur deposition prediction models often relied on the solubility of sulfur in acidic gas and the stress state of sulfur particles to determine the occurrence of deposition, thus establishing predictive models. However, in the presence of complex geological conditions, the multiphase flow through porous media and the adsorption of particles on pore throat walls can also influence sulfur particle deposition to some degree. It is well known that sulfur particle deposition during gas reservoir development exhibits instability, with multiple factors influencing the deposited sulfur particles. Particularly noteworthy is the influence of airflow velocity, which can resuspend sulfur particles that are physically adsorbed on pore throat surfaces, thereby reintegrating them into the gas phase. Additionally, the dynamic deposition of larger sulfur particles involves a dynamic process. This study elucidates the dynamic process of sulfur deposition by considering the diverse transport dynamics of sulfur particles. Physical adsorption and desorption behaviors of sulfur particles are determined based on variations in reservoir conditions. The desorption status of sulfur particles with different particle sizes within the formation is established by evaluating the equilibrium between the force exerted on the pore throat wall and the suspension force generated by gas flow. The critical conditions for sulfur deposition in Yuanba gas reservoir were obtained by substituting on-site parameters into calculations. Moreover, a mathematical model is proposed to describe the dynamic deposition and migration of sulfur particles, adopting principles from continuous porous media porous flow theory, fluid flow mass conservation, as well as sulfur particle desorption and migration. The formulated model is solved, and its resulting solution process and outcomes hold significant implications for numerical simulation and predictive assessment of the development impact on gas reservoirs, particularly in later stages.

Study on the influence of side-blown airflow velocities on plasma and combustion wave generated from fused silica induced by combined pulse laser

This study examines the impact of variations in side-blowing airflow velocity on plasma generation, combustion wave propagation mechanisms, and surface damage in fused silica induced by a combined millisecond-nanosecond pulsed laser. The airflow rate and pulse delay are the main experimental variables. The evolution of plasma motion was recorded using ultrafast time-resolved optical shadowing. The experimental results demonstrate that the expansion velocities of the plasma and combustion wave are influenced differently by the side-blowing airflow at different airflow rates (0.2 Ma, 0.4 Ma, and 0.6 Ma). As the flow rate of the side-blow air stream increases, the initial expansion velocities of the plasma and combustion wave gradually decrease, and the side-blow air stream increasingly suppresses the plasma. It is important to note that the target vapor is always formed and ionized into plasma during the combined pulse laser action. Therefore, the side-blown airflow alone cannot completely clear the plasma. Depending on the delay conditions, the pressure of the side-blowing airflow, the influence of inverse Bremsstrahlung radiation absorption and target surface absorption mechanisms can lead to a phenomenon known as the double combustion waves when using a nanosecond pulse laser. Both simulation and experimental results are consistent, indicating the potential for further exploration of fused silica targets in the laser field.

EXPERIMENTAL STUDIES OF THE PERFORMANCE EFFICIENCY OF A WIND TURBINE WITH COMBINED BLADES

This article studies the aerodynamic characteristics of a wind turbine at different streamline parameters. For this purpose, an experimental sample of the plant with combined power elements in rotating cylinders form with fixed blades was fabricated. The airflow velocity was varied, ranging from 3 to 12 m/s. The results of the experiment to study the variation of the angle α of the fixed blade position relative to the cylinder from the air flow velocity were analyzed. The changes in aerodynamic forces from the air flow velocity are graphically presented. It was found that at the optimal angle α = 0 ° of the fixed blade relative to the cylinder, the maximum values of aerodynamic forces were obtained. Graphs of the dependence of aerodynamic coefficients on the Reynolds number are also constructed, in which it is established that, at α= 0º, the minimum value of the lift coefficient is 0.012 N and the maximum value of the drag coefficient is 10.07 N at Re = 1∙105. The presented results show the effectiveness of the combined use of the and the fixed blade. The experimental results obtained on the aerodynamic parameters of a wind turbine can be used in the development of prototypes of installations designed for low wind speeds.

Research of dust removal performance and power output characteristics on photovoltaic panels by longitudinal high-speed airflow

Photovoltaic (PV) panels' photoelectric conversion efficiency will decrease as dust deposition on their surface. An approach to dust removal on the PV panel's surface by longitudinal high-speed airflow was investigated to increase the output power. In this paper, commercial CFD software was used to numerically simulate the characteristics of dust removal by longitudinal high-speed airflow. The PV panels' output characteristic model under the action of dust removal by high-speed airflow is established by Simulink software, and the influence of high-speed airflow on the output characteristic is studied. Firstly, the dust removal mechanism of the PV panel under high-speed airflow is studied. Secondly, the impact of different tilt angles of the PV panel, dust particle size, airflow velocity, and blowing time on the dust removal effect of the PV panel surface were studied. The optimal airflow velocity and blowing time were obtained according to the airflow velocity and blowing time and their influence on the dust removal rate. Finally, the effect of high-speed airflow dust removal on PV power generation's efficiency and output characteristics is studied. The findings demonstrate that as the tilt angle increases, so does the dust removal rate increases continuously. In addition, the rate at which dust is removed increases with airflow velocity and blowing time. When the tilt angle is 30°, the best airflow velocity of dust removal is 5 m/s, and the best blowing time is 8 s. When the tilt angle is 45°, the best airflow velocity of dust removal is 10 m/s, and the best blowing time is 5 s. With the increase in airflow velocity and blowing time, the output power of PV panels continues to increase. When the airflow velocity is 10 m/s, and the blowing time is 5 s, 7.5 s, 12.5 s, and 15 s, the maximum output power after dust removal is increased by 15.44 %, 23.05 %, 23.38 %, and 23.39 %. This paper's research results can guide the design of the practical engineering application of longitudinal blowing high-speed airflow in the dust removal of PV panels.

Modelling the Smoke Flow Characteristics of a Comprehensive Pipe Gallery Fire with Rectangular Section

In this study, a numerical model of the cable cabin of a comprehensive pipe gallery was established to study the smoke flow diffusion behaviour of a comprehensive pipe gallery fire under a rectangular cross-section. The effects of fire source power (Q = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0 MW) and fire source location (D = 10, 20, 40, 50, 60, 80, 100 m) on the smoke flow characteristics—such as smoke layer height and thickness, longitudinal airflow velocity, and ceiling temperature distribution—were analysed, and the corresponding prediction model was fitted. The results show the following: (1) The height of the smoke layer decreases with increasing fire power, and the predictive model of the smoke layer thickness obtained from the fitting is proportional to the smoke mass flow rate and inversely proportional to the aspect ratio of the pipe gallery. (2) Longitudinal air velocity prediction models of D < 50 m and D ≥ 50 m are fitted, and the average error between them and the numerical simulation values is 9.611%. (3) The temperature decay gradient of the smoke decreases gradually with increasing distance from the fire source, while there is a significant temperature difference between the two sides of the fire source. The average relative errors of the dimensionless temperature rise models fitted upstream and downstream of the fire source in the form of ΔTT0=AeBD−XH+C exponentials with respect to the numerical simulations were 11.688% and 7.296%, respectively. The results of the study can provide a reference for smoke flow and fire prevention and control in comprehensive pipe galleries.

Design and Evaluation of Face Mask Filtration: Mechanisms, Formulas, and Fluid Dynamics Simulations

The global adoption of face masks as a preventive measure against the spread of the SARS-CoV-2 virus (COVID-19) has spurred extensive research into their filtration efficacy. This study begins by elucidating various mechanisms of particle penetration and comparing filtration efficiency formulas with experimental data from prior studies. This is compared to the filtration efficiency experimental measurement developed in our previous study. Moreover, it delves into fluid dynamics simulations to examine different turbulent airflow models. Specifically, it contrasts the airflow velocity distribution of the k-ω and k-ε turbulent flow models with that of a quadrant-based average velocity model developed within this research. Furthermore, the study conducts fluid dynamic simulations to assess airflow profiles for six distinct medical and non-medical face masks. The results underscore substantial disparities among the simulations, emphasising the criticality of employing accurate fluid dynamics models for evaluating airflow patterns during diverse respiratory activities such as breathing, coughing, or sneezing, thereby enhancing environmental health in the realm of infectious disease prevention.

Vertical transmission of aerosols between building flats through drainage system: A review

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been demonstrate to be transmitted through fecal aerosols in building drainage systems (BDS), emphasizing the urgent need to understand its specific transmission pathways and mechanisms. This paper concludes the two periods of aerosols transmission in BDS: airflows caused by pressure fluctuations during toilet flushing and airflows induced by the chimney effect. It also compares the airflow velocity, air pressure, and transmission duration during these two periods. Subsequently, research methods (e.g., aerosol experiments, tracer gas experiments, and numerical simulations), influencing factors, and related conclusions regarding aerosols transmission in BDS were thoroughly concluded and analyzed. Finally, the advantages and limitations of the research methods are compared, offering valuable insights on essential future research directions and drainage system design. Although there is a qualitative understanding of the vertical transmission mechanism in BDS, which summarize the pathway of aerosols transmission as “respiratory/fecal aerosols-drainage pipes-faulty floor drains-bathrooms”, it is imperative to establish a quantitative understanding of the relative importance and duration of aerosols transmission during these two periods. Further research urgently needs to assess the risk levels of aerosols transmission during both periods and establish pressure-coupling curves for these two periods. Overall, both designers and scholars can utilize this review as a reference to mitigate the potential risk of vertical aerosols transmission in BDS, so that BDS can be appropriately designed to minimize pressure fluctuations and the water seal in floor drains can be regularly checked.

ИССЛЕДОВАНИЕ ТЕПЛООБМЕНА АГРЕГАТОВ ТРАНСМИССИИ И ДВИГАТЕЛЯ ГРУЗОВОГО АВТОМОБИЛЯ

A significant share of goods in agriculture are transported by road at low ambient temperatures. Changing the properties of oils in these conditions significantly reduces the efficiency of the machines and requires additional measures for the thermal preparation of the units. (Research purpose) The research purpose is identifying heat flows between the engine and the transmission units of a truck and their effect on the thermal mode of the middle and rear axles. (Materials and methods) We performed experimental studies on a KamAZ 65115 car. Convective heat flows and heat transfer through thermal conductivity in the "engine-transmission units" system, as well as their effect on the thermal mode of transmission gearboxes, were studied. (Results and discussion) The effect of the engine on the thermal mode of the gearbox was established: the oil temperature rise was five degrees; the oil temperature stabilization time was 70-74 minutes. The effect of the engine on the thermal mode of the middle and rear axles was not revealed. It was shown that the main heat flows causing the gearbox to heat up are the air flows washing the block and crankcase of the engine. An increase in the air flow velocity to 25 meters per second leads to a decrease in temperature near the engine block by 2.5 percent to 278 kelvin, and near the crankcase by 3 percent to 253 kelvin. (Conclusions) Engine temperature has a significant effect only on the temperature of the gearbox. At air speeds above 10 meters per second, stabilization of transmission temperature is observed. Heat flows that cause additional heating of the gearbox are formed on the side surfaces of the cylinder block and the engine oil sump.

Removal efficiency of restroom ventilation revisited for short-term evaluation

Ventilation efficiency or contaminant removal efficiency is often evaluated using the ratio between the concentrations in the exhaust air and the room air. This ratio does not truly represent the expectation of ventilation in restrooms, where dynamic airflow fields and sources are more typical. This study focuses on a short-term (10 min) pollutant removal percentage in a residential restroom featuring a dynamic airflow field, particularly with the onset of window-induced stack ventilation during toilet uses. Thirteen ventilation scenarios of a residential restroom were studied using the numerical method that was validated by a mock-up experiment. The scenarios differed in the operation of the exhaust fan and window. Results show that the 10-min pollutant removal percentage of a typical exhaust ventilation system at 10 h-1 air change rate (ACH) is only 68.5%. Under exhaust ventilation, opening the window can introduce both adverse short circuit and favorable stack ventilation depending on the difference between the indoor and outdoor temperatures. As the temperature difference increases from 0 to 12.5 °C, the removal percentage increases from below 50%, a drop due to short circuit, to above 98% thanks to a tripled ventilation rate. The human thermal plume has notable effect on the removal percentage, but its effect can be neglected with the presence of stack ventilation. The hybrid ventilation strategy has impact on perceived air quality and thermal comfort. When the outdoor air is colder, opening the window under exhaust ventilation may increase the current sitting user’s exposure to the self-produced pollutants but can reduce the exposure of the next immediate standing user. In addition, opening the window in cold days will make the toilet user thermally uncomfortable with reduced local temperatures and increased airflow velocities. The study highlights the importance of using the short-term removal percentage to evaluate the performance of restroom ventilation.

Nasal Airflow Dynamics following LeFort I Advancement in Cleft Nasal Deformities: A Retrospective Preliminary Study.

Unilateral cleft lip and palate (UCLP) nasal deformity impacts airflow patterns and pressure distribution, leading to nasal breathing difficulties. This study aims to create an integrated approach using computer-aided design (CAD) and computational fluid dynamics (CFD) to simulate airway function and assess outcomes in nasal deformities associated with unilateral cleft lip and palate (UCLP) after LeFort I osteotomy advancement. Significant alterations were observed in nasal geometry, airflow velocity, pressure dynamics, volumetric flow rate, and nasal resistance postoperatively, indicating improved nasal airflow. The cross-sectional area increased by 26.6%, airflow rate by 6.53%, and nasal resistance decreased by 6.23%. The study offers quantitative insights into the functional impacts of such surgical interventions, contributing to a deeper understanding of UCLP nasal deformity treatment and providing objective metrics for assessing surgical outcome.

CFD MODELLING OF THE DUST AND AIR VELOCITY BEHAVIOUR AT AN UNDERGROUND COAL MINE ROADWAY

Dust pollution for ventilation systems is a big problem for mine workers in underground coal mines. The respirable dust concentration is important to the physical and spiritual health of miners. It is hazardous up to 2 mg/m3in the air in underground coal galleries. In the present work, the airflow velocity and dust concentration values were measured at Çay 1-Kartiye roadway in Kozlu, Zonguldak. The goal of this research is modelling to the dust concentration and airflow velocity values by the computational fluid dynamics (CFD) method. An airflow and dust fluid and solid model was done using the Eularian-Granular multiphase method with the k-epsilon two-equation model. The results showed that dust concentration values were too low and fixed. The error percentage of airflow velocity values changed between 5% and 17%. Thus, this roadway is safe for coal dust explorations.