ABSTRACT This work applies the concepts of the attainable region for process synthesis in comminution. The attainable region analysis has been successfully applied for process synthesis of reactor networks. The Attainable Region is defined as the set of all possible output states for a constrained or unconstrained system of fundamental processes (Horn, 1964). A basic procedure for constructing the attainable region for the fundamental processes of reaction and mixing has been postulated in reaction engineering (Glasser et al., 1987). This procedure has been followed in this work to construct the candidate attainable region for size reduction processes as found in a size reduction environment. A population balance model has been used to characterise the evolution of particle size distributions from a comminution event. Herbst and Fuerstenau (1973) postulated the dependency of grinding on the specific energy. A specific energy dependent population balance model was used for the theoretical simulations and for the fitting of experimental data. A new method of presenting particle size distributions as points in the Euclidian space was postulated in place of the traditional cumulative distribution. This allows successive product particle size distributions to be connected forming a trajectory over which the objective function can be evaluated. The curve connects products from successive batch grinding stages forming a pseudo-continuous process. Breakage, mixing and classification were identified as the fundamental processes of interest for comminution. Agglomeration was not considered in any of the examples. Mathematical models were used to describe each fundamental process, i.e. breakage, mixing and classification, and an The application of the attainable region analysis in comminution Abstract algorithm developed that could calculate the evolution of product particle size distributions. A convex candidate attainable region was found from which process synthesis and optimisation solutions could be drawn in two dimensional Euclidian space. As required from Attainable Region Theory, the interior of the bounded region is filled by trajectories of higher energy requirements or mixing between two boundary optimal points. Experimental validation of the proposed application of the attainable region analysis results in comminution was performed. Mono-sized feed particles were broken in a laboratory ball mill and the products were successfully fitted using a population balance model. It was shown that the breakage process trajectories were convex and they follow first order grinding kinetics at long grind times. The candidate attainable region was determined for an objective function to maximise the mass fraction in the median size class 2. It was proved that the same specific energy input produces identical products. The kinematic and loading conditions are supposed to be chosen as a subsequent event after the required specific energy is identified. Finally the fundamental process of classification was added to the system of breakage and mixing. The attainable regions analysis affords the opportunity to quantify exactly the reduction in energy consumption due to classification in a comminution circuit, thus giving optimal targets. Classification showed the potential to extend the candidate attainable region for a fixed specific energy input. The boundary of the attainable region is interpreted as pieces of equipment and optimum process conditions. This solves both the original process synthesis and successive optimisation problems.

The theoretical and experimental investigation of the novel in-line CMOS image sensor was performed. This sensor is aimed to investigate particle size distribution, particle concentration and shape in hydraulic liquid in order to implement the proactive maintenance of hydraulic equipment. The existing instruments such as automatic particle counters and techniques are not sufficiently enough to address this task because of their restricted sensitivity, limit of concentration to be measured and they cannot determine particle shape. Other instruments cannot be used as inline sensors because they are not resistant to the arduous conditions such as high pressure and vibration. The novel mathematical model was proposed as it is not possible to use previously developed techniques based on using optical system and complicated algorithms. This model gives the output signal of the image sensor depending on the particle size, its distance from the light source (LED) and image sensor. Additionally, the model takes into account the limited exposure time and particle track simulation. The results of simulation based on the model are also performed in thesis. On the basis of the mathematical model the image processing algorithms were suggested in order to determine particle size even when this size is lower than pixel size. There are different approaches depending on the relation between the size of the particle and the pixel size. The approach to the volume of liquid sample estimation was suggested in order to address the problem of low accuracy of concentration measurement by the conventional automatic particle counters based on the single photodiode. Proposed technique makes corrections on the basis of particle velocity estimation. Approach to the accuracy estimation of the sensor was proposed and simulation results are shown. Generally, the accuracy of particle size and concentration measurement was considered. Ultimately, the experimental setup was used in order to test suggested techniques. The mathematical model was tested and the results showed sufficient correlation with the experiment. The zinc dust was used as a reference object as there are the particles within the range from 1 to 25 microns which is appropriate to check the sensitivity. The results of experiments using reference instrument showed the improved sensitivity and accuracy of volume measured compared to the reference one.

Com a crescente diminuição de recursos minerais e o alto custo envolvido na construção da estrutura de uma mina, a recuperação máxima possível de uma jazida vem se tornando fundamental. Para isso além da escolha do método de lavra ter a necessidade de ser feito cautelosamente, é possível lançar mão de métodos adicionais de recuperação, como por exemplo, a recuperação de pilares. Essa pesquisa foi baseada na determinação da caracterização do enchimento (backfill cimentado) utilizado nas câmaras vazias que possibilita a posterior recuperação dos pilares. A caracterização do enchimento é composta da determinação da resistência simples do backfill necessária para que o enchimento cumpra com seu objetivo, desenvolvimento da classificação granulométrica ótima para os agregados e dosagem de cimento e água para alcançar a resistência proposta. A metodologia desenvolvida para obter a nova caracterização é composta de várias etapas que incluem pesquisas em campo e trabalhos em laboratório. Primeiramente, foram obtidos através de análise em campo os parâmetros de dosagem de cimento e classificação granulométrica dos agregados já utilizados na planta de fabricação do enchimento, bem como sua resistência correspondente. Em seguida definições teóricas da dosagem de cimento ideal e classificação granulométrica ótima foram realizadas com base na resistência à compressão simples que foi identificada como necessária para cumprir com as solicitações geomecânicas do maciço rochoso, então posteriormente, a nova caracterização definida teoricamente foi posta à prova através da confecção de corpos de prova de backfill, seguido de execuções de ensaios de compressão. Durante a primeira etapa da metodologia, já se pôde identificar que os agregados possuíam um alto índice de partículas tamanho argila que estavam afetando os resultados de resistência obtidos com a caracterização empregada inicialmente. A partir disso se optou por construir a curva granulométrica ótima sem essa fração. A resistência à compressão simples calculada de 2,69 MPa, foi obtida com base no planejamento de longo prazo que prevê a total recuperação dos pilares existentes na mina. Dessa maneira toda a área que será minerada foi considerada como um único bloco. Finalmente, foi identificada a dosagem de cimento sendo de 4% em peso, que juntamente com a granulometria ótima é capaz de alcançar os valores esperados de resistência. Para que o planejamento da produção da mina durante os próximos anos de vida útil seja efetivamente cumprido, o enchimento deverá prover à mina estabilidade geomecânica local a nível de câmaras abertas com paredes verticais de backfill estáveis e também estabilidade global a nível de contato entre níveis e galerias de acesso. Isso somente será alcançado se a nova caracterização for corretamente aplicada. / As a consequence of the ongoing reduction of mineral resources and the high cost involved in the construction of a mine, the maximum recovery of a mineral deposit becomes a fundamental issue. Therefore, besides the need of caution on the choice of the mining method, it is possible to make use of additional recovery methods, such as the recovery of pillars. This research was based on the determination of the characterization of the fill (cemented backfill) used in avoid stopes that allows the subsequent recovery of adjacent pillars. The characterization of the fill consists of determining the uniaxial compressive strength of the backfill required for an efficient filling, developing an optimal particle-size distribution for the aggregates and finding the cement-water ratio necessary to reach the desired resistance. The methodology developed to obtain the new characterization is comprised of several steps which include field work and laboratory tests. First, cement dosing parameters and particle size of the aggregates (already used at the filling manufacturing plant), as well as their corresponding strength, were obtained through analyses in the field work. Then, theoretical definitions of the ideal cement dosing and optimal particle-size analysis were carried out based on the uniaxial compressive strength that has been identified as necessary to comply with the geomechanical requests from the rock mass, and then later, the new theoretical characterization was tested by making backfill samples, followed by execution of compression tests. During the first stage of this methodology, it has been identified a high proportion of clay particle size for the aggregates, that have affected the strength results obtained from the characterization used initially. From this point, we decided to build the optimal particle-size curve without this fraction. Uniaxial compressive strength, calculated as 2.69 MPa, was obtained from the long-term planning that determines the full recovery of the existing pillars in the mine. In this way, the entire area to be mined was considered as a single block. Finally, the cement dosing has been identified as 4% by weight, which together with the optimal particle size, is able to achieve the expected strength values. In order to effectively fulfill the mine production planning over the next years of lifespan, the filling should provide the mine local geomechanical stability at open stopes level, with vertical walls of stable backfill, and also global stability at the contacts between levels and access galleries. This will only be achieved if the new characterization is correctly applied.

A curva de retenção da água no solo é um dos principais instrumentos para avaliar a qualidade física dos solos e possibilitar seu manejo adequado. Por meio da Teoria da Capilaridade vários equipamentos foram desenvolvidos para determinar a intensidade com que a água está retida ao solo, porém, pouco se tem dado atenção para verificar se os pressupostos para o real funcionamento da teoria estão sendo atendidos. Um aspecto refere-se ao tamanho da amostra utilizada para determinar a curva de retenção, de modo que haja continuidade dos feixes capilares na amostra e placa porosa. Desta forma, este trabalho propõe avaliar diferentes tamanhos de amostra indeformada para a determinação da curva de retenção. Para isso, coletaram-se amostras em anéis volumétricos cilíndricos de três tamanhos (altura) diferentes (T1 - 0,075 m; T2 - 0,05 m; T3 - 0,025 m;) e mesmo diâmetro interno (0,07 m), dos horizontes diagnósticos de um Latossolo e um Nitossolo em áreas experimentais da Escola Superior de Agricultura Luiz de Queiroz (Esalq/USP), Piracicaba - SP. Realizou-se a caracterização física destes solos, por meio da análise granulométrica, densidade do solo, densidade de partículas, porosidade total e teor de carbono orgânico. As curvas foram determinadas para cada tamanho de amostra, utilizando-se o Funil de Haines, para as tensões 0,5, 1, 4, 6 e 10 kPa, e a Câmara de Pressão de Richards para 33, 100 e 500 kPa. As curvas de retenção foram ajustadas pelo modelo utilizado por van Genuchten. Estimadas as curvas, avaliou-se a distribuição de poros do solo das amostras, determinando-se a curva de frequência acumulada de poros em função do logaritmo do raio e, depois pela diferenciação das equações de ajuste das curvas de retenção, a curva diferencial de frequência acumulada de poros. Os resultados mostram que o Latossolo, por ter textura arenosa no horizonte estudado, não apresentou diferença significativa nas curvas de retenção para os tamanhos das amostras estudadas. Verificou-se pouca modificação na distribuição dos poros deste solo, que possui teor elevado das frações areia fina e muito fina, e desenvolveram papel importante para a retenção de água. O Nitossolo, por sua vez, apresentou diferença significativa da curva obtida pela amostra de menor tamanho (T3), havendo maior retenção de água com a diminuição do tamanho da amostra. Devido a sua textura muito argilosa, o arranjo estrutural deste solo foi diferenciado ao se utilizar as amostras maiores, com provável interrupção e descontinuidade dos feixes capilares. Consequentemente, houve também alteração na distribuição dos poros, com redução dos mesoporos e aumento dos microporos. Desta forma, pode-se concluir que o tamanho da amostra influenciou a curva de retenção da água devido à complexidade estrutural do solo, que provavelmente é diferente nas amostras maiores por causa da continuidade dos feixes capilares, principalmente no Nitossolo. Em outras palavras, quanto menor o tamanho da amostra há menor diferenciação no arranjo de poros, ou seja, maior proximidade da real condição do solo e, assim, uma interpretação da retenção de água \"mais correta\" por meio da Teoria da Capilaridade. / The soil water retention curve is one of main tools to assess the physical quality of the soil and to make possible its adequate management. By means of the Capillary Theory, many instruments have been developed to determine the water retention forces in soil, but, little attention has been given to check whether the assumptions for the application of the theory are being attended. One aspect relates to the sample size used to determine the retention curve, so that there is capillary continuity of the sample and porous plate. Thus, this study aimed to evaluate different sizes of undisturbed cylindrical samples for determination of the retention curve. The samples were collected from diagnostic horizons of Latosol and a Nitosol, in experimental areas of the Escola Superior de Agricultura Luiz de Queiroz (Esalq/USP), Piracicaba - SP. Three volumetric rings with three different heights (T1 - 0,075 m; T2 - 0,05 m; T3 - 0,025 m;), were used the diameter of the rings were the same for the three sizes (0,07 m). The physical characterization soils were made by the granulometric analysis, bulk density, particle density, porosity and organic carbon. The curves were determined for each sample size, using the Haine\'s funnel, for tensions 0.5, 1, 4, 6 and 10 kPa, and Richard\'s pressure chamber for 33, 100 and 500 kPa. The pore size distribution of the soil was evoluated first by determining the cumulative frequency curve of pore radius, and then, by differentiating fitling equation of fluis curve to obtain a differential curve of pore cumulative frequency. The retention curves were fitted by the model used by Van Genuchtem. The results showed that in the Latosol there was no significant difference in retention curves for the sizes of the samples studied. There was little change in pore distribution of the fluis soil with high content of fractions fine and very fine sand, that developes significant role in soil water retention. The Nitosol exhibited significant difference of the curve obtained by the smallest size sample (T3), with greater water retention with the decreasing of the sample size. Because of its very clayey texture, the structural arrangement of this soil was different when larger samples were used, with probable interruption and discontinuity of capillaries. Consequently there was also a change in the pore distribution, with reduction of mesopores and an increase of micropores. Thus, itean be conclude that the sample size influenced the soil water retention curve due to the structural complexity of the soil that probably is different in the bigger sample because of the continuity of the capillary budles fluit was more affected in the Nitosol. In other words, the smaller the sample size, the smaller the difference in the arrangement of pores, that is, more closeness of the actual condition of the soil and so a \"more correct\" interpretation of soil water retention by capillary theory.

Approximately 1500 km (43%) of the strategic Scottish trunk road network drainage asset takes the form of Highway Filter Drains (HFDs). However, despite their popularity, they are prone to clogging, therefore they have an estimated operational life-cycle of ten-years. This research was undertaken to investigate the complex inter-relationship between catchment characteristics, road dynamics and the physical characteristics of the road to establish the key factors that govern the generation and spatial variability of Road-Deposited-Sediment (RDS). The aim being to establish the impact RDS Particle Size Distribution (PSD) has on clogging and the operational life-cycle of HFDs. The research adopted an integrated approach, incorporating: (i) a Field Study to investigate RDS PSD grading envelopes across the Scottish trunk road network, (ii) a HFD Field Survey (HFD-FS) to investigate the current condition of HFDs across a range of catchments, (iii) a HFD Field Study to establish whether, or not, graded stone PSD envelopes comply with specification requirements when first placed in the trench, (iv) 3 HFD Field Studies to assess the level of risk of system failure through evaluating the causes and quantifying the individual, cumulative and influencing factors which contribute to the evolution of clogging in HFDs, and (v) a Field Study utilising Ground Penetrating Radar data to explore why stratified (clogged) layers evolve within some HFDs. Based on the results of this research, design and maintenance procedures were then identified to improve the performance of HFDs. Results established that Scottish trunk roads operate under variable catchment characteristics and road dynamics, with the condition and specification of the road surface, volume of traffic, road geometry, number of running lanes and adjacent land use governing the generation and spatial variability of RDS. Five of the nine trunk road RDS PSD profiles shared a peak particle diameter of 425 μm, two had a peak of 600 μm, one had a peak of 1180 μm and one resulted in a peak of 2120 μm. Particles >1000 μm were mostly mineral or asphalt and it was shown that there is a direct link between the factors that govern the generation and spatial variability of RDS and those that govern the evolution of clogging and actual operational life-cycle of HFDs. The HFD-FS revealed that 69% were assigned Filter Drain Condition Index ratings of 3 or below, which identifies these as having exceeded the ten-year estimated operational life-cycle. 94% of those deemed to have reached the end of their operational life-cycle were over-the-edge (OTE) HFDs, which supports the assertion that pre-treatment would increase their operational life-cycle. Results also established that introducing a kerb-line and gully-pots or grass-strip between the road and the HFD significantly reduced the indices of particle size composition d50 and d90 and percentage of RDS retained at depths spanning 0 - 400 mm, compared to OTE HFDs with comparable catchment characteristics and road dynamics. Based on these results, the operational life-cycle of HFDs with a kerb-line and gully-pots and OTE HFDs with a grass-strip can be expected to exceed twenty-years, if catchment characteristics and road dynamics are representative of those in this study. This research also identified that compacting Type B graded stone with ‘heavy vibrating machinery’ during construction could potentially contribute to clogging. It was also established that HFD harrowing may exacerbate clogging because the process of disintegrating the cake-layer mobilises an otherwise rigid and compacted RDS mass and this is more likely to penetrate deeper into the HFD and inundate the HFD during a storm event. These findings indicate that current HFD construction and maintenance practice could have a detrimental impact on the effective operational life-cycle of HFDs. Overall, this research study has demonstrated that there are considerable uncertainties related to PSD grading envelopes and percentage of RDS migrating from roads to HFDs. It is clear therefore that one of the most notable findings of this research is that given the scale of strategic trunk road networks, assuming a single HFD operational life-cycle profile, for a trunk road or trunk road network, is highly unlikely to be representative of a HFD at the local level. It follows then that the widely accepted estimated ten-year operational life-cycle for HFDs, does not reflect the actual operational life-cycle of HFDs.

This study compares four different methods to measure hydraulic conductivity (K) at two sites on the East Tennessee State University Valleybrook Campus. It compares the K values to each other, to the different K values between the two sites, and to United States Department of Agriculture (USDA) K values. Two field methods, Well Bail Test and Auger Hole Test, and two lab methods, Constant Head Permeameter Test and Grain Size Distribution Test (GSD), were performed on the clay rich Ultisol soils on an East Tennessee hillslope in the Valley and Ridge Physiographic Province. One site was located close to a monitoring well and the other on the floodplain of an existing stream. The Hazen, Alyamani & Sen, and Slichter methods were used to compute K from the GSD Test. The Alyamani & Sen, Slichter, and permeameter methods produced similar K values ranging from 9.52 x 10-6 to 1.25 x 10-3 cm/sec. These are similar to the USDA K values ranging from 9.17 x 10-4 to 2.82 x 10-4 cm/sec. The Hazen method overestimated K and ranged from 8.10 x 10-3 to 1.09 x 10-1 cm/sec. The Well Bail Test yielded a lower K value (ranging from 8.16 x 10-9 to 1.19 x 10-8 cm/sec) than the USDA values as expected for water flow in deeper soil horizons at a depth of 8.50 meters. Comparing these values helped to better understand the difference between various methods to compute the hydraulic conductivity.

The interception and redirection of rainfall by vegetation has implications for many fields such as remote sensing of soil moisture, satellite observation of rainfall, and the modeling of runoff, climate, and soil erosion. Although the modeling of rainfall partitioning by forests has received attention in the past, partitioning caused by crops has been overlooked. The present work proposes a two front experimental and computational methodology to comprehensively study rainfall interception and partitioning by the maize canopy. In the experimental stage, we deployed two compact weather stations, two optical disdrometers, and five tipping bucket rain gauges. Two of the tipping bucket rain gauges were modified to measure throughfall while two were adapted to measure stemflow. The first optical disdrometer allowed for inspection of the unmodified drop-size and velocity distributions, whereas the second disdrometer measured the corresponding distributions under the canopy. This indicates that the outcome of the interaction between the hydrometeors and the canopy depends on the drop diameter. In the computational stage, we created a model that uses drop-size and velocity distributions as well as a three-dimensional digital canopy to simulate the movement of raindrops on the surfaces of leaves. Our model considers interception, redirection, retention, coalescence, breakup, and re-interception of drops to calculate the stemflow, throughfall, and equivalent height of precipitation stored on plants for a given storm. Moreover, the throughfall results are presented as two-dimensional matrices, where each term corresponds to the accumulated volume of drops that dripped at a given location. This allows insight into the spatial distribution of throughfall beneath the foliage. Finally, we examine the way in which the maize canopy modifies the drop-size distribution by recalculating the drop velocity based on the raindrop's size and detachment height and by storing the counts of drops in diameter-velocity classes that are consistent with the classes used by disdrometers in the experimental study.

The overall goal of this doctoral dissertation is to develop a prototype instrument, a Portable Aerosol Collector and Spectrometer (PACS), that can continuously measure aerosol size distributions by number, surface area and mass concentrations over a wide size range (from 10 nm to 10 µm) while also collecting particles with impactor and diffusion stages for post-sampling chemical analyses. To achieve the goal, in the first study, we designed, built and tested the PACS hardware. The PACS consists of a six-stage particle size selector, a valve system, a water condensation particle counter to measure number concentrations and a photometer to measure mass concentrations. The valve system diverts airflow to pass sequentially through upstream stages of the selector to the detectors. The stages of the selector include three impactor and two diffusion stages, which resolve particles by size and collect particles for chemical analysis. Particle penetration by size was measured through each stage to determine actual performance and account for particle losses. The measured d50 of each stage (aerodynamic diameter for impactor stages and geometric diameter for diffusion stages) was similar to the design. The pressure drop of each stage was sufficiently low to permit its operation with portable air pumps. In the second study, we developed a multi-modal log-normal (MMLN) fitting algorithm to leverage the multi-metric, low-resolution data from one sequence of PACS measurements to estimate aerosol size distributions of number, surface area, and mass concentration in near-real-time. The algorithm uses a grid-search process and a constrained linear least-square (CLLS) solver to find a tri-mode (ultrafine, fine, and coarse), log-normal distribution that best fits the input data. We refined the algorithm to obtain accurate and precise size distributions for four aerosols typical of diverse environments: clean background, urban and freeway, coal power plant, and marine surface. Sensitivity studies were conducted to explore the influence of unknown particle density and shape factor on algorithm output. An adaptive process that refined the ranges and step sizes of the grid-search reduced the computation time to fit a single size distribution in near-real-time. Assuming standard density spheres, the aerosol size distributions fit well with the normalized mean bias (NMB) of -4.9% to 3.5%, normalized mean error (NME) of 3.3% to 27.6%, and R2 values of 0.90 to 1.00. The fitted number and mass concentration biases were within ± 10% regardless of uncertainties in density and shape. With this algorithm, the PACS is able to estimate aerosol size distributions by number, surface area, and mass concentrations from 10 nm to 10 µm in near-real-time. In the third study, we developed a new algorithm–the mass distribution by composition and size (MDCS) algorithm–to estimate the mass size distribution of various particle compositions. Then we compared the PACS for measuring multi-mode aerosols to three reference instruments, including a scanning mobility particle sizer (SMPS), an aerodynamic particle sizer (APS) and a nano micro-orifice uniform deposit impactor (nanoMOUDI). We used inductively coupled plasma mass spectrometry to measure the mass of collected particles on PACS and nanoMOUDI stages by element. For the three-mode aerosol, the aerosol size distributions in three metrics measured with the PACS agreed well with those measured with the SMPS/APS: number concentration, bias = 9.4% and R2 = 0.96; surface area, bias = 17.8%, R2 = 0.77; mass, bias = -2.2%, R2 = 0.94. Agreement was considerably poorer for the two-mode aerosol, especially for surface area and mass concentrations. Comparing to the nanoMOUDI, for the three-mode aerosol, the PACS estimated the mass median diameters (MMDs) of the coarse mode well, but overestimated the MMDs for ultrafine and fine modes. The PACS overestimated the mass concentrations of ultrafine and fine mode, but underestimated the coarse mode. This work provides insight into a novel way to simultaneously assess airborne aerosol size, composition, and concentration by number, surface area and mass using cost-effective handheld technologies.

The objective of this work is to optimize high-power low-frequency sonication (US) pretreatment of sludge, and especially to investigate for the first time possible improvements by higher pressure and audible frequency. After a preliminary examination of regular process conditions (sludge conditioning, sludge type, prior alkalization, temperature control, etc), effects of US parameters (power -PUS, intensity -IUS, specific energy input -ES, frequency -FS, etc.) and of hydrostatic pressure (Ph) were specifically looked into, separately and in combination, first under cooling at constant temperature (28°C), then under the progressive temperature rise provoked by sonication. First, it was confirmed that specific energy input (ES) plays a key role in sludge US disintegration (i.e. solubilisation of organic matter) and that temperature rise during adiabatic-like sonication is beneficial through additional effects of thermal hydrolysis and cavitation. At a given ES value, low FS (12 kHz vs. 20 kHz) and high PUS enhance soluble chemical oxygen demand (SCOD) due to more violent cavitation, while hydrostatic pressure gives rise to an optimum value due to its opposite effects on cavitation threshold and intensity. One major result is that optimal pressure depends on IUS (P¬US) as well as temperature profile, but not on ES, FS, nor sludge type. Setting the other parameters at the most favorable conditions expected, i.e. 12 kHz, 360 W , 28 gTS/L, and adiabatic conditions, final optimization was achieved by searching for this pressure optimum and examining sequential procedure to avoid too high temperature dampening cavitation intensity and damaging the transducer. Such conditions with sequential mode and Ph of 3.25 bar being selected succeeded in achieving very high SCOD, but only marginally improved subsequent methanization yield.

In this study, we developed and standardized a sanding method to evaluate the emission of airborne particles from products that contain carbon nanotubes (CNTs) under different conditions, including three types of sandpaper and three sanding disc speed. We also characterized the emission of the airborne particles from one neat epoxy test sample, four CNTs-incorporating test samples with different CNTs loading, and two commercial products. The total number concentration, respirable mass concentration, and particle size number/mass distribution of the emitted particles were calculated and compared, followed by an electron microscopy (EM) analysis. These data suggest that the sanding process can produce substantial quantities of airborne particles. Also, the emission of airborne particles was associated with different test conditions. EM analysis of the airborne particle samples showed embedded CNTs protruding from the outer surface, which was different from CNTs-incorporating bulk material. Our study suggests a potential generation of particles during the life cycle event of sanding. Further studies should be carried out to investigate the potential human health hazard in other life cycle events.