Heat flux measurement inside internal combustion engine with gradient heat flux sensor.

This master’s thesis is devoted to study different heat flux measurement techniques such as differential temperature sensors, semi-infinite surface temperature methods, calorimetric sensors and gradient heat flux sensors. The possibility to use Gradient Heat Flux Sensors (GHFS) to measure heat flux in the combustion chamber of compression ignited reciprocating internal combustion engines was considered in more detail. A. Mityakov conducted an experiment, where Gradient Heat Flux Sensor was placed in four stroke diesel engine Indenor XL4D to measure heat flux in the combustion chamber. The results which were obtained from the experiment were compared with model’s numerical output. This model (a one – dimensional single zone model) was implemented with help of MathCAD and the result of this implementation is graph of heat flux in combustion chamber in relation to the crank angle. The values of heat flux throughout the cycle obtained with aid of heat flux sensor and theoretically were sufficiently similar, but not identical. Such deviation is rather common for this type of experiment.

Effect of different surface roughnesses on a turbulent boundary layer

The classical treatment of rough wall turbulent boundary layers consists in determining the effect the roughness has on the mean velocity profile. This effect is usually described in terms of the roughness function delta U+. The general implication is that different roughness geometries with the same delta U+ will have similar turbulence characteristics, at least at a sufficient distance from the roughness elements. Measurements over two different surface geometries (a mesh roughness and spanwise circular rods regularly spaced in the streamwise direction) with nominally the same delta U+ indicate significant differences in the Reynolds stresses, especially those involving the wall-normal velocity fluctuation, over the outer region. The differences are such that the Reynolds stress anisotropy is smaller over the mesh roughness than the rod roughness. The Reynolds stress anisotropy is largest for a smooth wall. The small-scale anisotropy and interniittency exhibit much smaller differences when the Taylor microscale Reynolds number and the Kolmogorov-normalized mean shear are nominally the same. There is nonetheless evidence that the small-scale structure over the three-dimensional mesh roughness conforms more closely with isotropy than that over the rod-roughened and smooth walls.

Modeling of turbulent flow through radial diffuser

The work considers the modeling of turbulent flow in radial diffuser with axial feeding. Due to its claimed capability to predict flow including features such as separation, curvature and adverse pressure gradient, the RNG k-epsilon model of Orzag et al. (1993) is applied in the present analysis. The governing equations are numerically solved using the finite volume methodology. Experiments were conducted to assess the turbulence model. Numerical results of pressure distribution on the front disk surface for different flow conditions when compared to the experimental data indicated that the RNG k-epsilon model is adequate to predict this class of flow.

Turbulent shallow-water model for orographic subgrid-scale perturbations

A parallel pseudo-spectral method for the simulation in distributed memory computers of the shallow-water equations in primitive form was developed and used on the study of turbulent shallow-waters LES models for orographic subgrid-scale perturbations. The main characteristics of the code are: momentum equations integrated in time using an accurate pseudo-spectral technique; Eulerian treatment of advective terms; and parallelization of the code based on a domain decomposition technique. The parallel pseudo-spectral code is efficient on various architectures. It gives high performance onvector computers and good speedup on distributed memory systems. The code is being used for the study of the interaction mechanisms in shallow-water ows with regular as well as random orography with a prescribed spectrum of elevations. Simulations show the evolution of small scale vortical motions from the interaction of the large scale flow and the small-scale orographic perturbations. These interactions transfer energy from the large-scale motions to the small (usually unresolved) scales. The possibility of including the parametrization of this effects in turbulent LES subgrid-stress models for the shallow-water equations is addressed.

An experimental/numerical study of the turbulent boundary layer development along a surface with a sudden change in roughness

A theory for the description of turbulent boundary layer flows over surfaces with a sudden change in roughness is considered. The theory resorts to the concept of displacement in origin to specify a wall function boundary condition for a kappa-epsilon model. An approximate algebraic expression for the displacement in origin is obtained from the experimental data by using the chart method of Perry and Joubert(J.F.M., vol. 17, pp. 193-122, 1963). This expression is subsequently included in the near wall logarithmic velocity profile, which is then adopted as a boundary condition for a kappa-epsilon modelling of the external flow. The results are compared with the lower atmospheric observations made by Bradley(Q. J. Roy. Meteo. Soc., vol. 94, pp. 361-379, 1968) as well as with velocity profiles extracted from a set of wind tunnel experiments carried out by Avelino et al.( 7th ENCIT, 1998). The measurements are found to be in good agreement with the theoretical computations. The skin-friction coefficient was calculated according to the chart method of Perry and Joubert(J.F.M., vol. 17, pp. 193-122, 1963) and to a balance of the integral momentum equation. In particular, the growth of the internal boundary layer thickness obtained from the numerical simulation is compared with predictions of the experimental data calculated by two methods, the "knee" point method and the "merge" point method.

Turbulent heat transfer and pressure drop in pinned annular regions

Experiments were performed to determine average heat transfer coefficients and friction factors for turbulent flow through annular ducts with pin fins. The measurements were carried out by means of a double-pipe heat exchanger. The total number of pins attached to the inner wall of the annular region was 560. The working fluids were air, flowing in the annular channel, and water through the inner circular tube. The average heat transfer coefficients of the pinned air-side were obtained from the experimental determination of the overall heat transfer coefficients of the heat exchanger and from the knowledge of the average heat transfer coefficients of the circular pipe (water-side), which could be found in the pertinent literature. To attain fully developed conditions, the heat exchanger was built with additional lengths before and after the test section. The inner circular duct of the heat exchanger and the pin fins were made of brass. Due to the high thermal conductivity of the brass, the small tube thickness and water temperature variation, the surface of the internal tube was practically isothermal. The external tube was made of an industrial plastic which was insulated from the environment by means of a glass wool batt. In this manner, the outer surface of the annular channel can be considered adiabatic. The results are presented in dimensionless forms, in terms of average Nusselt numbers and friction factors as functions of the flow Reynolds number, ranging from 13,000 to 80,000. The pin fin efficiency, which depends on the heat transfer coefficient, is also determined as a function of dimensionless parameters. A comparison of the present results with those for smooth sections (without pins) is also presented. The purpose of such a comparison is to study the influence of the presence of the pins on the pressure drop and heat transfer rate.

Is there correlation between the turbulent eddies size and mechanical hemolysis ?

Hemolytic profile of an artificial device chronically implanted in the cardiovascular system may represent the difference between the success and failure in its long-term performance. Last decades have witnessed efforts on the development of methods capable of predicting red blood cell damage in artificial organs. However, all of them have had limited success to predict hemolysis. The primary cause of this problem is that such models do not take into consideration structures of turbulent flow. The present paper demonstrates that microscopic measurable occurrences of the turbulent flow may be linked to red blood cell trauma. This study suggests that if the smallest turbulent eddies dimension is under 10 m m hemolysis is not dependent on the exposure time and the red blood cells damage depends only on the dissipation of the turbulent energy in the erythrocyte membrane. The analysis reported here opens the possibility of mapping the flow field in artificial assist devices based on the smallest eddy length scales. This is a promising new trend and should be considered in the designing requirements of the next generations of artificial organs.

Fluid dynamic modelling and simulation of circulating fluidized bed reactors: importance of the interface turbulence transfer

The main objective of this work is to analyze the importance of the gas-solid interface transfer of the kinetic energy of the turbulent motion on the accuracy of prediction of the fluid dynamic of Circulating Fluidized Bed (CFB) reactors. CFB reactors are used in a variety of industrial applications related to combustion, incineration and catalytic cracking. In this work a two-dimensional fluid dynamic model for gas-particle flow has been used to compute the porosity, the pressure, and the velocity fields of both phases in 2-D axisymmetrical cylindrical co-ordinates. The fluid dynamic model is based on the two fluid model approach in which both phases are considered to be continuous and fully interpenetrating. CFB processes are essentially turbulent. The model of effective stress on each phase is that of a Newtonian fluid, where the effective gas viscosity was calculated from the standard k-epsilon turbulence model and the transport coefficients of the particulate phase were calculated from the kinetic theory of granular flow (KTGF). This work shows that the turbulence transfer between the phases is very important for a better representation of the fluid dynamics of CFB reactors, especially for systems with internal recirculation and high gradients of particle concentration. Two systems with different characteristics were analyzed. The results were compared with experimental data available in the literature. The results were obtained by using a computer code developed by the authors. The finite volume method with collocated grid, the hybrid interpolation scheme, the false time step strategy and SIMPLEC (Semi-Implicit Method for Pressure Linked Equations - Consistent) algorithm were used to obtain the numerical solution.

Numerical study of wedge supported oblique shock wave-oblique detonation wave transitions

The results of a numerical study of premixed Hydrogen-air flows ignition by an oblique shock wave (OSW) stabilized by a wedge are presented, in situations when initial and boundary conditions are such that transition between the initial OSW and an oblique detonation wave (ODW) is observed. More precisely, the objectives of the paper are: (i) to identify the different possible structures of the transition region that exist between the initial OSW and the resulting ODW and (ii) to evidence the effect on the ODW of an abrupt decrease of the wedge angle in such a way that the final part of the wedge surface becomes parallel to the initial flow. For such a geometrical configuration and for the initial and boundary conditions considered, the overdriven detonation supported by the initial wedge angle is found to relax towards a Chapman-Jouguet detonation in the region where the wedge surface is parallel to the initial flow. Computations are performed using an adaptive, unstructured grid, finite volume computer code previously developed for the sake of the computations of high speed, compressible flows of reactive gas mixtures. Physico-chemical properties are functions of the local mixture composition, temperature and pressure, and they are computed using the CHEMKIN-II subroutines.

Characterisation of waste for combustion : with special reference to the role of zinc

Waste combustion has gone from being a volume reducing discarding-method to an energy recovery process for unwanted material that cannot be reused or recycled. Different fractions of waste are used as fuel today, such as; municipal solid waste, refuse derived fuel, and solid recovered fuel. Furthermore, industrial waste, normally a mixture between commercial waste and building and demolition waste, is common, either as separate fuels or mixed with, for example, municipal solid waste. Compared to fossil or biomass fuels, waste mixtures are extremely heterogeneous, making it a complicated fuel. Differences in calorific values, ash content, moisture content, and changing levels of elements, such as Cl and alkali metals, are common in waste fuel. Moreover, waste contains much higher levels of troublesome trace elements, such as Zn, which is thought to accelerate a corrosion process. Varying fuel quality can be strenuous on the boiler system and may cause fouling and corrosion of heat exchanger surfaces. This thesis examines waste fuels and waste combustion from different angles, with the objective of giving a better understanding of waste as an important fuel in today’s fuel economy. Several chemical characterisation campaigns of waste fuels over longer time periods (10-12 months) was used to determine the fossil content of Swedish waste fuels, to investigate possible seasonal variations, and to study the presence of Zn in waste. Data from the characterisation campaigns were used for thermodynamic equilibrium calculations to follow trends and determine the effect of changing concentrations of various elements. The thesis also includes a study of the thermal behaviour of Zn and a full—scale study of how the bed temperature affects the volatilisation of alkali metals and Zn from the fuel. As mixed waste fuel contains considerable amounts of fresh biomass, such as wood, food waste, paper etc. it would be wrong to classify it as a fossil fuel. When Sweden introduced waste combustion as a part of the European Union emission trading system in the beginning of 2013 there was a need for combustion plants to find a usable and reliable method to determine the fossil content. Four different methods were studied in full-scale of seven combustion plants; 14Canalysis of solid waste, 14C-analysis of flue gas, sorting analysis followed by calculations, and a patented balance method that is using a software program to calculate the fossil content based on parameters from the plant. The study showed that approximately one third of the coal in Swedish waste mixtures has fossil origins and presented the plants with information about the four different methods and their advantages and disadvantages. Characterisation campaigns also showed that industrial waste contain higher levels of trace elements, such as Zn. The content of Zn in Swedish waste fuels was determined to be approximately 800 mg kg-1 on average, based on 42 samples of solid waste from seven different plants with varying mixtures between municipal solid waste and industrial waste. A review study of the occurrence of Zn in fuels confirmed that the highest amounts of Zn are present in waste fuels rather than in fossil or biomass fuels. In tires, Zn is used as a vulcanizing agent and can reach concentration values of 9600-16800 mg kg-1. Waste Electrical and Electronic Equipment is the second Zn-richest fuel and even though on average Zn content is around 4000 mg kg-1, the values of over 19000 mg kg-1 were also reported. The increased amounts of Zn, 3000-4000 mg kg-1, are also found in municipal solid waste, sludge with over 2000 mg kg-1 on average (some exceptions up to 49000 mg kg-1), and other waste derived fuels (over 1000 mg kg-1). Zn is also found in fossil fuels. In coal, the average level of Zn is 100 mg kg-1, the higher amount of Zn was only reported for oil shale with values between 20-2680 mg kg-1. The content of Zn in biomass is basically determined by its natural occurrence and it is typically 10-100 mg kg-1. The thermal behaviour of Zn is of importance to understand the possible reactions taking place in the boiler. By using thermal analysis three common Zn-compounds were studied (ZnCl2, ZnSO4, and ZnO) and compared to phase diagrams produced with thermodynamic equilibrium calculations. The results of the study suggest that ZnCl2(s/l) cannot exist readily in the boiler due to its volatility at high temperatures and its conversion to ZnO in oxidising conditions. Also, ZnSO4 decomposes around 680°C, while ZnO is relatively stable in the temperature range prevailing in the boiler. Furthermore, by exposing ZnO to HCl in a hot environment (240-330°C) it was shown that chlorination of ZnO with HCl gas is possible. Waste fuel containing high levels of elements known to be corrosive, for example, Na and K in combination with Cl, and also significant amounts of trace elements, such as Zn, are demanding on the whole boiler system. A full-scale study of how the volatilisation of Na, K, and Zn is affected by the bed temperature in a fluidised bed boiler was performed parallel with a lab-scale study with the same conditions. The study showed that the fouling rate on deposit probes were decreased by 20 % when the bed temperature was decreased from 870°C to below 720°C. In addition, the lab-scale experiments clearly indicated that the amount of alkali metals and Zn volatilised depends on the reactor temperature.

Model based analysis of the post-combustion calcium looping process for carbon dioxide capture

This thesis presents a one-dimensional, semi-empirical dynamic model for the simulation and analysis of a calcium looping process for post-combustion CO2 capture. Reduction of greenhouse emissions from fossil fuel power production requires rapid actions including the development of efficient carbon capture and sequestration technologies. The development of new carbon capture technologies can be expedited by using modelling tools. Techno-economical evaluation of new capture processes can be done quickly and cost-effectively with computational models before building expensive pilot plants. Post-combustion calcium looping is a developing carbon capture process which utilizes fluidized bed technology with lime as a sorbent. The main objective of this work was to analyse the technological feasibility of the calcium looping process at different scales with a computational model. A one-dimensional dynamic model was applied to the calcium looping process, simulating the behaviour of the interconnected circulating fluidized bed reactors. The model incorporates fundamental mass and energy balance solvers to semi-empirical models describing solid behaviour in a circulating fluidized bed and chemical reactions occurring in the calcium loop. In addition, fluidized bed combustion, heat transfer and core-wall layer effects were modelled. The calcium looping model framework was successfully applied to a 30 kWth laboratory scale and a pilot scale unit 1.7 MWth and used to design a conceptual 250 MWth industrial scale unit. Valuable information was gathered from the behaviour of a small scale laboratory device. In addition, the interconnected behaviour of pilot plant reactors and the effect of solid fluidization on the thermal and carbon dioxide balances of the system were analysed. The scale-up study provided practical information on the thermal design of an industrial sized unit, selection of particle size and operability in different load scenarios.

Behavior of black liquor nitrogen in combustion : formation of cyanate

The Kraft pulping process is the dominant chemical pulping process in the world. Roughly 195 million metric tons of black liquor are produced annually as a by-product from the Kraft pulping process. Black liquor consists of spent cooking chemicals and dissolved organics from the wood and can contain up to 0.15 wt% nitrogen on dry solids basis. The cooking chemicals from black liquor are recovered in a chemical recovery cycle. Water is evaporated in the first stage of the chemical recovery cycle, so the black liquor has a dry solids content of 65-85% prior to combustion. During combustion of black liquor, a portion of the black liquor nitrogen is volatilized, finally forming N2 or NO. The rest of the nitrogen remains in the char as char nitrogen. During char conversion, fixed carbon is burned off leaving the pulping chemicals as smelt, and the char nitrogen forms mostly smelt nitrogen (cyanate, OCN-). Smelt exits the recovery boiler and is dissolved in water. The cyanate from smelt decomposes in the presence of water, forming NH3, which causes nitrogen emissions from the rest of the chemical recovery cycle. This thesis had two focuses: firstly, to determine how the nitrogen chemistry in the recovery boiler is affected by modification of black liquor; and secondly, to find out what causes cyanate formation during thermal conversion, and which parameters affect cyanate formation and decomposition during thermal conversion of black liquor. The fate of added biosludge nitrogen in chemical recovery was determined in Paper I. The added biosludge increased the nitrogen content of black liquor. At the pulp mill, the added biosludge did not increase the NO formation in the recovery boiler, but instead increased the amount of cyanate in green liquor. The increased cyanate caused more NH3 formation, which increased the NCG boiler’s NO emissions. Laboratory-scale experiments showed an increase in both NO and cyanate formation after biosludge addition. Black liquor can be modified, for example by addition of a solid biomass to increase the energy density of black liquor, or by separation of lignin from black liquor by precipitation. The precipitated lignin can be utilized in the production of green chemicals or as a fuel. In Papers II and III, laboratory-scale experiments were conducted to determine the impact of black liquor modification on NO and cyanate formation. Removal of lignin from black liquor reduced the nitrogen content of the black liquor. In most cases NO and cyanate formation decreased with increasing lignin removal; the exception was NO formation from lignin lean soda liquors. The addition of biomass to black liquor resulted in a higher nitrogen content fuel mixture, due to the higher nitrogen content of biomass compared to black liquor. More NO and cyanate were formed from the fuel mixtures than from pure black liquor. The increased amount of formed cyanate led to the hypothesis that black liquor is catalytically active and converts a portion of the nitrogen in the mixed fuel to cyanate. The mechanism behind cyanate formation during thermal conversion of black liquor was not clear before this thesis. Paper IV studies the cyanate formation of alkali metal loaded fuels during gasification in a CO2 atmosphere. The salts K2CO3, Na2CO3, and K2SO4 all promoted char nitrogen to cyanate conversion during gasification, while KCl and CaCO3 did not. It is now assumed that cyanate is formed when alkali metal carbonate or an active intermediate of alkali metal carbonate (e.g. -CO2K) reacts with the char nitrogen forming cyanate. By testing different fuels (bark, peat, and coal), each of which had a different form of organic nitrogen, it was concluded that the form of organic nitrogen in char also has an impact on cyanate formation. Cyanate can be formed during pyrolysis of black liquor, but at temperatures 900°C or above, the formed cyanate will decompose. Cyanate formation in gasifying conditions with different levels of CO2 in the atmosphere was also studied. Most of the char nitrogen was converted to cyanate during gasification at 800-900°C in 13-50% CO2 in N2, and only 5% of the initial fuel nitrogen was converted to NO during char conversion. The formed smelt cyanate was stable at 800°C 13% CO2, while it decomposed at 900°C 13% CO2. The cyanate decomposition was faster at higher temperatures and in oxygen-containing atmospheres than in an inert atmosphere. The presence of CO2 in oxygencontaining atmospheres slowed down the decomposition of cyanate. This work will provide new information on how modification of black liquor affects the nitrogen chemistry during thermal conversion of black liquor and what causes cyanate formation during thermal conversion of black liquor. The formation and decomposition of cyanate was studied in order to provide new data, which would be useful in modeling of nitrogen chemistry in the recovery boiler.

Heterogeneous mass transfer in fluidized beds by computational fluid dynamics

The main objective of this research is to estimate and characterize heterogeneous mass transfer coefficients in bench- and pilot-scale fluidized bed processes by the means of computational fluid dynamics (CFD). A further objective is to benchmark the heterogeneous mass transfer coefficients predicted by fine-grid Eulerian CFD simulations against empirical data presented in the scientific literature. First, a fine-grid two-dimensional Eulerian CFD model with a solid and gas phase has been designed. The model is applied for transient two-dimensional simulations of char combustion in small-scale bubbling and turbulent fluidized beds. The same approach is used to simulate a novel fluidized bed energy conversion process developed for the carbon capture, chemical looping combustion operated with a gaseous fuel. In order to analyze the results of the CFD simulations, two one-dimensional fluidized bed models have been formulated. The single-phase and bubble-emulsion models were applied to derive the average gas-bed and interphase mass transfer coefficients, respectively. In the analysis, the effects of various fluidized bed operation parameters, such as fluidization, velocity, particle and bubble diameter, reactor size, and chemical kinetics, on the heterogeneous mass transfer coefficients in the lower fluidized bed are evaluated extensively. The analysis shows that the fine-grid Eulerian CFD model can predict the heterogeneous mass transfer coefficients quantitatively with acceptable accuracy. Qualitatively, the CFD-based research of fluidized bed process revealed several new scientific results, such as parametrical relationships. The huge variance of seven orders of magnitude within the bed Sherwood numbers presented in the literature could be explained by the change of controlling mechanisms in the overall heterogeneous mass transfer process with the varied process conditions. The research opens new process-specific insights into the reactive fluidized bed processes, such as a strong mass transfer control over heterogeneous reaction rate, a dominance of interphase mass transfer in the fine-particle fluidized beds and a strong chemical kinetic dependence of the average gas-bed mass transfer. The obtained mass transfer coefficients can be applied in fluidized bed models used for various engineering design, reactor scale-up and process research tasks, and they consequently provide an enhanced prediction accuracy of the performance of fluidized bed processes.

Analysis and modelling of chemical looping combustion process with and without oxygen uncoupling

Effective control and limiting of carbon dioxide (CO₂) emissions in energy production are major challenges of science today. Current research activities include the development of new low-cost carbon capture technologies, and among the proposed concepts, chemical combustion (CLC) and chemical looping with oxygen uncoupling (CLOU) have attracted significant attention allowing intrinsic separation of pure CO₂ from a hydrocarbon fuel combustion process with a comparatively small energy penalty. Both CLC and CLOU utilize the well-established fluidized bed technology, but several technical challenges need to be overcome in order to commercialize the processes. Therefore, development of proper modelling and simulation tools is essential for the design, optimization, and scale-up of chemical looping-based combustion systems. The main objective of this work was to analyze the technological feasibility of CLC and CLOU processes at different scales using a computational modelling approach. A onedimensional fluidized bed model frame was constructed and applied for simulations of CLC and CLOU systems consisting of interconnected fluidized bed reactors. The model is based on the conservation of mass and energy, and semi-empirical correlations are used to describe the hydrodynamics, chemical reactions, and transfer of heat in the reactors. Another objective was to evaluate the viability of chemical looping-based energy production, and a flow sheet model representing a CLC-integrated steam power plant was developed. The 1D model frame was succesfully validated based on the operation of a 150 kWth laboratory-sized CLC unit fed by methane. By following certain scale-up criteria, a conceptual design for a CLC reactor system at a pre-commercial scale of 100 MWth was created, after which the validated model was used to predict the performance of the system. As a result, further understanding of the parameters affecting the operation of a large-scale CLC process was acquired, which will be useful for the practical design work in the future. The integration of the reactor system and steam turbine cycle for power production was studied resulting in a suggested plant layout including a CLC boiler system, a simple heat recovery setup, and an integrated steam cycle with a three pressure level steam turbine. Possible operational regions of a CLOU reactor system fed by bituminous coal were determined via mass, energy, and exergy balance analysis. Finally, the 1D fluidized bed model was modified suitable for CLOU, and the performance of a hypothetical 500 MWth CLOU fuel reactor was evaluated by extensive case simulations.

Relationship between pulmonary function and indoor air pollution from coal combustion among adult residents in an inner-city area of southwest China

Few studies evaluate the amount of particulate matter less than 2.5 mm in diameter (PM2.5) in relation to a change in lung function among adults in a population. The aim of this study was to assess the association of coal as a domestic energy source to pulmonary function in an adult population in inner-city areas of Zunyi city in China where coal use is common. In a cross-sectional study of 104 households, pulmonary function measurements were assessed and compared in 110 coal users and 121 non-coal users (≥18 years old) who were all nonsmokers. Several sociodemographic factors were assessed by questionnaire, and ventilatory function measurements including forced vital capacity (FVC), forced expiratory volume in 1 s (FEV1), the FEV1/FVC ratio, and peak expiratory flow rate (PEFR) were compared between the 2 groups. The amount of PM2.5 was also measured in all residences. There was a significant increase in the relative concentration of PM2.5 in the indoor kitchens and living rooms of the coal-exposed group compared to the non-coal-exposed group. In multivariate analysis, current exposure to coal smoke was associated with a 31.7% decrease in FVC, a 42.0% decrease in FEV1, a 7.46% decrease in the FEV1/FVC ratio, and a 23.1% decrease in PEFR in adult residents. The slope of lung function decrease for Chinese adults is approximately a 2-L decrease in FVC, a 3-L decrease in FEV1, and an 8 L/s decrease in PEFR per count per minute of PM2.5 exposure. These results demonstrate the harmful effects of indoor air pollution from coal smoke on the lung function of adult residents and emphasize the need for public health efforts to decrease exposure to coal smoke.