A Study on the Structure of Paper: The Links between Paper and Fibre Properties by Francisco J. Fernandez Madrigal B.Sc. Mechanical Engineering, Universidad Simon Bolivar, 2010 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF Master of Applied Science in THE FACULTY OF GRADUATE STUDIES (Mechanical Engineering) The University Of British Columbia (Vancouver) April 2013 c Francisco J. Fernandez Madrigal, 2013 Abstract The work shown in this thesis, focuses on the effects of fibre geometry (fibre length, width and coarseness) to the paper strength and bulk. It is believed that fibre networks, for a given furnish, can be optimized through fractionation to create a stronger and bulkier paper. A first set of experi- ments was performed to determine the correlations among fibre properties of a softwood kraft pulp furnish. The fibre properties were measured using a Fibre Quality Analyzer (FQA). Coarseness (w) and width (D) were found to increase linearly with length (L f ). These correlations are thought to be influenced by the tree species and the pulping process. A second set of experiments was aimed to determine empirical expressions of bulk and tensile index (TI) in terms of the fibre geom- etry distribution (w , L f and D) and the press-drying pressure (P). Bulk and TI was measured for handsheets made at different pressures from different size distributions. These distributions were created by a combination of fractionation on a Bauer McNett Classifier (BMC) and fibre cutting. The determined relations show agreement with experimental work from other researchers. In or- der to provide insight into the causes of such behavior, simulations of two-dimensional (2D) and three-dimensional (3D) random networks were performed. Matlab was used to write the simula- tion codes. Mechanical response of fibres (affected by the fibre geometry) to the forces induced by drying or pressure is not considered directly in the simulations. However, 2D simulation models are a good representation of high press-drying conditions and high fibre flexibility, whereas 3D models are better predicting air-dried (low pressure) paper made from fibres of high rigidity. Geo- metric statistics of random networks explain some of the experimental observations. The statistical outputs of the 2D simulations were; network coverage (c), network thickness (tN), number of fibre crossings (Nc) and the relative bonded area (RBA). In the case of 3D simulations, only tN was determined. The 2D and 3D outputs were measured for different fibre size distributions (L f , D and w). ii Table of Contents Abstract . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ii Table of Contents . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . iii List of Tables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . v List of Figures . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . vi Nomenclature . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . xi Acknowledgments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . xiii 1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 1.1 The Structure of Pulp Suspensions, Wet Webs and Paper . . . . . . . . . . . . . . 1 1.2 Fibre Morphology and Fibre Properties . . . . . . . . . . . . . . . . . . . . . . . 3 1.3 A Brief Description of the Paper-making Process . . . . . . . . . . . . . . . . . . 5 1.4 Fibre and Paper Properties . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8 1.4.1 Property Distributions . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8 1.4.2 Basis Weight . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9 1.4.3 Paper Strength . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9 1.4.4 Fibre Strength . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10 1.4.5 Paper Thickness and Bulk . . . . . . . . . . . . . . . . . . . . . . . . . . 10 1.4.6 Relative Bonded Area . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10 1.5 Brief Summary of Determined Models and Studies on Paper and Fibre Properties . 11 1.5.1 Experimental Predictions . . . . . . . . . . . . . . . . . . . . . . . . . . . 11 1.5.2 The Page Equation for Strength . . . . . . . . . . . . . . . . . . . . . . . 11 iii 1.5.3 Statistical Models and Computer Simulations . . . . . . . . . . . . . . . . 12 1.6 Research Purpose . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 2 Predicting Bulk and Strength of Paper . . . . . . . . . . . . . . . . . . . . . . . . . . 16 2.1 Experimental Procedure . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16 2.2 Assumptions and Limitations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17 2.3 The Fiber Distributions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18 2.3.1 Length-Weighted Average Values . . . . . . . . . . . . . . . . . . . . . . 19 2.3.2 Property Distributions . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19 2.4 The Natural Links in Fiber Properties . . . . . . . . . . . . . . . . . . . . . . . . 23 2.5 Multi-variable Regression . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25 2.6 Zero Span Tensile Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29 2.7 Combining the Natural Links with the Regression . . . . . . . . . . . . . . . . . . 30 3 Computer Simulations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33 3.1 Paper as a 2D Random Network . . . . . . . . . . . . . . . . . . . . . . . . . . . 33 3.1.1 Assumptions and Limitations of the 2D Model . . . . . . . . . . . . . . . 34 3.1.2 Paper Thickness and the Pixel Value distribution . . . . . . . . . . . . . . 36 3.1.3 The Degree of Bonding . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38 3.2 Two-Dimensional (2D) Deposition . . . . . . . . . . . . . . . . . . . . . . . . . . 42 3.3 The Three-Dimensional (3D) Structure of Paper . . . . . . . . . . . . . . . . . . . 46 3.3.1 Assumptions and Limitations . . . . . . . . . . . . . . . . . . . . . . . . 46 3.3.2 The Three-Dimensional (3D) Algorithm . . . . . . . . . . . . . . . . . . . 48 3.3.3 Length Effect . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49 3.3.4 Degree of Bonding of a Three-Dimensional (3D) Network . . . . . . . . . 49 4 Conclusions and Final Remarks . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51 4.1 Recommendations for Future Work . . . . . . . . . . . . . . . . . . . . . . . . . . 53 Bibliography . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55 A Experimental Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57 B Matlab Codes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62 iv List of Tables Table 2.1 Size andw measurements of different pulps obtained from FQA. Errors are from standard deviation between FQA measurements. The table shows consistent correlations among all properties (increasing D and w with L f ) except for CI. Fibre cutting or fractionation did not affectCI significantly. . . . . . . . . . . . 18 Table 2.2 Coefficients for proposed bulk and tensile index models (equations 2.4 and 2.5) determined through linear regression. . . . . . . . . . . . . . . . . . . . . . . . 26 Table 2.3 Standard deviation of length-weighted pulp distributions. Obtained using equa- tion 1.8 and the length-weighted distribution shown in equation 1.9 . . . . . . . 26 Table 2.4 Improved coefficients for bulk and tensile index models (2.4 and 2.5) as ob- tained from regression. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27 Table 3.1 Description of the pixel value distribution P(x;y). The bonding area factor is multiplied by the pixel area (Ap) to calculate the bonded area at the (x;y) location. 34 Table A.1 Experimental handsheet data: bulk, TI and ZI for the feed pulp and different cuts 58 Table A.2 Experimental handsheet data: bulk, TI and ZI for the M14 fraction and different cuts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59 Table A.3 Experimental handsheet data: bulk, TI and ZI for the M28 fraction and different cuts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60 Table A.4 Experimental handsheet data: bulk, TI and ZI for the M48 fraction and different cuts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61 v List of Figures Figure 1.1 Micrograph (obtained using light microscopy) of the top side of a standard kraft pulp handsheet. The ribbon-like structures are collapsed wood fibres, held together by hydrogen bonds. . . . . . . . . . . . . . . . . . . . . . . . . 1 Figure 1.2 Micrograph (obtained using light microscopy) of the cross-section of a stan- dard kraft pulp handsheet. Although most fibres are collapsed, lumen is visible in some of them. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 Figure 1.3 Micrograph (obtained through light microscopy) of a kraft pulp fibre. The fibre is in a water solution. At this state the fibre is swollen and has a cylindrical shape. The lumen is irregular in size and changing along the length of the fibre. 3 Figure 1.4 Fibre geometric models. Kraft pulp fibres in a suspension are cylindrical (left figure). Mechanical pulp fibres tend to remain rectangular (right figure). How- ever, after drying, both can be modeled as rectangular because of the high collapse index (f ) of kraft pulp fibres due to their low lignin content. . . . . . 4 Figure 1.5 Depiction of a curled fibre. The pulping process causes the fibre to curl affect- ing the paper properties. The amount of curl is characterized by the curl index, CI ( SC f 1). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 Figure 1.6 Flow diagram of the paper-making process (commercial and standard hand- sheet). After the fibres are separated from the wood chips to form a suspension, water is taken from the mixture by a combination of gravity, suction, pressure and heat. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6 vi Figure 2.1 Flow diagram of experimental procedure. The original feed pulp was separated into five fractions. Handsheets were made at different pressures from some of the fractions (Feed, M14,M28 and M48). Air-dried samples were cut and repulped to create new handsheets. Tensile index and bulk for all handsheets was recorded and fibre properties were measured for all fractions. . . . . . . . 17 Figure 2.2 Length distributions (length-weighted) of different fibre fractions obtained from a BMC fractionation. The mean of each distribution shifts as the mesh is finer. All appear to have Gaussian behavior, with a consistent reduction in standard deviation as the mesh fineness increases. . . . . . . . . . . . . . . . . . . . . . 20 Figure 2.3 Length distribution (length-weighted) of a BMC fraction and its respective cuts. Fibre cutting was effective by reducing LWL of each fraction, shown here as shifts in the distributions. Furthermore, they remain Gaussian without considerably increasing standard deviation. . . . . . . . . . . . . . . . . . . . 20 Figure 2.4 Length-weighted pulp distributions for different cuts (in sequential order: M14, M28, M48 and Feed Pulp). As the mesh is finer the cutting effect is lost and the resulting distribution is wider. . . . . . . . . . . . . . . . . . . . . . . . . 21 Figure 2.5 Scatter plot of width D and length from fibres for different distributions (Feed Pulp, M14, M28, M48, M100 and M200 in sequential order). Although there is a broad range of width for fibres at certain length, there is a higher concen- tration at a given value. This is an indication that the width and length are correlated. This is regardless of the fibre fraction. Thus, there is a relation among fibre properties, that may be unaffected by the fractionation process. . . 22 Figure 2.6 Relation between fibre coarseness and length. Coarseness was obtained from the total measured length and the dry weight of each distribution sample. Thus, it is a mean value since there is a length distribution. A linear fit provided the highest coefficient of determination (R2 = 0:986). The standard errors, shown as bars in the plot, were obtained from standard deviation among FQA measurements. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23 vii Figure 2.7 Relation between fibre length and width (D). Constructed from raw data shown in Figure 2.5. The scatter points represent arithmetic averages of width among fibres in 0:2 mm bins. Left figure (a) shows width values for each BMC frac- tion. Right figure (b) shows the combined data of all BMC fractions. Mean fibre widths reported by FQA for all fractions have also been included. Fit 1 (equation 2.1 proposed by Madani [15]). Fit 2 (linear fit). The best fit was found to be linear (equation 2.2). . . . . . . . . . . . . . . . . . . . . . . . . 24 Figure 2.8 Effect of reducing the number of variables in the bulk model. The measured values of bulk were obtained from calliper and grammage measurements. Left figure (a) considers length and s , right figure (b) shows the improved model by eliminating length. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28 Figure 2.9 Effect of incorporating the standard deviation (s ) to the tensile index model. The measured values of tensile index were obtained from a tensile tester. Left figure (a) does not consider s as variable, right figure (b) shows the improved model by adding s . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28 Figure 2.10 The predicted value obtained through a linear regression of the zero span ten- sile index (ZI) against the measured values. The result of the large randomness of the measured ZI is a fit of very poor quality. . . . . . . . . . . . . . . . . . 30 Figure 2.11 Surface plot showing bulk for studied pulp in terms of pressure and length. Ob- tained by combining the model from regression (equation 2.6) and determined correlations between fibre properties (equations 2.2 and 2.3). . . . . . . . . . 31 Figure 2.12 Surface plot showing tensile index for studied pulp in terms of pressure and length at two s values ((a) 0.1 mm and (b) 10%). Obtained by combining the model from regression (equation 2.6) and determined correlations between fibre properties (equations 2.2 and 2.3). . . . . . . . . . . . . . . . . . . . . . 32 Figure 3.1 Simulated fibre networks at different coverage values. Left figure (a) c= 0:24 and right figure (b) c = 12. For both cases L f = 1 mm, D = 30 mm and the sample area is a 5 mm X 5 mm square . . . . . . . . . . . . . . . . . . . . . . 34 Figure 3.2 Difference in fibre contact between rigid and flexible fibres. a) Flexible fibre. b) Rigid Fibre. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35 Figure 3.3 The Poisson distribution of the pixel value (P(x;y)) obtained from simulations for two different coverage values. L f = 1 mm, D= 25 mm, N f = 3000 (a) and 10000 (b). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36 viii Figure 3.4 The effect of coverage (c) on the maximum pixel value (PMAX ) for different lengths (L f ) and widths (D). Values obtained from 2D simulations. PMAX changes only with c and not with L f or D. . . . . . . . . . . . . . . . . . . . . 38 Figure 3.5 The response of coverage(c) on total number of fibre-to-fibre crossings (Nc). Nc increases as fibre width (D) is reduced and coverage (c) is increased. . . . . 39 Figure 3.6 Comparison of number of fibre crossings obtained from simulations with sta- tistical formula (1.24) by M.Deng and Dodson [16]. The scattered values were obtained from simulations and the solid lines represent equation 1.24 at differ- ent widths. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40 Figure 3.7 The relative bonded area (RBA) against coverage (c) of simulated 2D networks for different lengths (L f ) and widths (D). RBA is only a function of c and unaffected by L f or D. The area along the thickness of the fibre (2 t f L f ) was not considered for RBA calculation. The solid line represents the statistical approximation provided by M.Deng and Dodson [16] as equation 1.23 . . . . . 42 Figure 3.8 Simulated network constructed from a Gaussian width fibre distribution. c = 0:969, bD= 45 mm and sD = 40 mm. . . . . . . . . . . . . . . . . . . . . . . . 43 Figure 3.9 Effect of width standard deviation (sD) of a normal distribution on the number of crossings Nc. The values of sD are in mm. sD shows no effect on Nc. Filled data points are for an average width (bD) of 30 mm and data points with no fill are for bD= 45 mm. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44 Figure 3.10 Simulated two-dimensional (2D) deposition. Fibres are considered to be rigid beams with high friction factor and high bonding force. Total fibre area is 74103 mm2, L f = 1 mm and t f = 15 mm. Resolution is 25 mm2 per pixel . . 44 Figure 3.11 Mean height of a rigid fibre 2D deposition against total area of deposited fibres. Data obtained from simulations. Each point represents an average of three simulations repeated under the same conditions. The fibre thickness t f is 15 mm. 45 Figure 3.12 Maximum height of a rigid fibre 2D deposition against total area of deposited fibres. Data obtained from simulations. Each point represents an average of three simulations repeated under the same conditions. The fibre thickness t f is 15 mm. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46 Figure 3.13 Simulated random fibre network generated in Matlab. The model used consid- ers fibres as rigid beams of high friction factor and bonding force. The fibre dimensions are L f = 1mm, D= 25 mm and t f = 15 mm. The figure shows 300 fibres deposited in a 5 mm square (10 GSM). . . . . . . . . . . . . . . . . . . 47 ix Figure 3.14 Top view image of the simulated random fibre network shown in Figure 3.13. Color has been added to show the height distribution in pixels. Resolution is 5:8107 mm3 per voxel. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47 Figure 3.15 Mean height of a simulated 3D random network (ctN) against fibre volume (N f L f D t f ). Fibre thickness and width are t f = 15 mm and D = 15 mm. Each point represents the average of three simulations and error bars are from standard deviation. The figure shows an increase inctN with L f and fibre volume. 49 Figure 3.16 Maximum height of simulated 3D random network (tN;MAX ) against fibre vol- ume (N f L f D t f ). t f = 15 mm and D = 15 mm. Each point represents the average of three simulations and error bars are from standard deviation. Comparing with Figure 3.15, as more fibres are added the difference between tN;MAX and ctN is reduced, yet for low fibre volumes, tN;MAX is in some cases ten foldsctN . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50 x Nomenclature Ac Fibre cross-sectional area g Acceleration of gravity Ap Pixel area GSM Paper grammage b Bond strength gD D exponential coefficient (bulk fit) BW Basis Weight gL L f exponential coefficient (bulk fit) bD D exponential coefficient (TI fit) gw w exponential coefficient (bulk fit) bL L f exponential coefficient (TI fit) gP P exponential coefficient (bulk fit) bw w exponential coefficient (TI fit) LA Dry lumen area bP P exponential coefficient (TI fit) LAo Wet lumen area bs s exponential coefficient (TI fit) LWL Length weighted average length C Pulp consistency m Distribution mean value c Coverage N Number of pixels ĉ Average coverage Nc Number of fibre crossings CB Bulk fit coefficient N f Number of fibres CI Curl index w Fibre coarseness C f Fibre chord P Press-drying pressure CT TI fit coefficient Pf Fibre cross-sectional perimeter D Fibre width PMAX Maximum pixel valuebD Mean fibre width P(x j) Probability of bin j Do Wet fibre width P(x;y) Pixel value distribution d Paper strip width RBA Relative bonded area F Breaking force r f Fibre density f Dried fibre voids fraction S Fibre chord length f (x) Probability density function s Standard deviationef (x) Weighted probability density function sD Width standard deviation f Collapse index T Breaking length xi TI Tensile index tN;MAX Maximum network thickness t f Fibre thickness Z Zero span breaking length tN(x;y) Network thickness distribution ZI Zero span tensile indexctN Mean network thickness xii Acknowledgments Many people have contributed in the development and review of the research shown in this thesis. I would like to acknowledge the assistance provided by the following: My supervisors: Dr. Mark Martinez and Dr. James Olson, for their continuous technical input. My co-researchers at the Pulp and Paper Center: Masita Mohamed, Jaime Sossa, Dr. Xinmei Shi, George Soong, Ario Madani and Ali Elahimehr. The broad variety in their backgrounds and re- search areas provided support in many aspects; solving technical issues, designing the experiments and assisting in the writing. Dr. Junaid Qazi, for his extensive technical support in analyzing all the data as well as his careful review of the content in this thesis. The thesis examiners: Dr. Peter Englezos and Dr. Douglas Romilly. My Friends and Family: For their unconditional support. Francisco J. Fernandez M. xiii Chapter 1 Introduction Paper is made mostly from wood fibres. Wood fibres are obtained by the pulping process of wood chips. In the wood chips, lignin, a glue-like substance holds the fibres together. By the use of me- chanical forces or chemicals, the chips are converted into individual fibres. Through this process, water is added to form a fibre-water mixture that is dried to make paper. Figure 1.1 and Figure 1.2 show the top view and cross-section of a standard handsheet made from softwood pulp. Figure 1.1: Micrograph (obtained using light microscopy) of the top side of a standard kraft pulp handsheet. The ribbon-like structures are collapsed wood fibres, held together by hydrogen bonds. 1.1 The Structure of Pulp Suspensions, Wet Webs and Paper A pulp suspension is a mixture of water and wood fibres. Ideally the fibres are sufficiently dispersed so that a small volume of the suspension is representative of the entire mixture. Consistency (C) is 1 Figure 1.2: Micrograph (obtained using light microscopy) of the cross-section of a standard kraft pulp handsheet. Although most fibres are collapsed, lumen is visible in some of them. commonly used to describe pulp suspensions. It is defined as the mass ratio of dry wood fibres to the total mass of the suspension: C = Dry Fiber Mass Contained in Suspension Total Mass o f Suspension 100 (1.1) Even at low consistencies, fibres have a tendency to cluster or flocculate due to fibre-fibre attractive forces. At very low consistencies, fibres (or fibre flocs) do not form a single connected structure. The way to characterize the connectivity or the formation of a fibre network is through the crowding number described by Kerekes and Schell [13]. As consistency is increased (approximately 30%) the mixture is no longer described as a sus- pension but as a wet web. In a wet web all the fibres are connected. There may not be physical contacts between the fibre surfaces, but the distance is small enough that the surrounding fibres interact constraining the motion of the fibre. However, in wet webs, consistency is low enough so that air and water is still present in the voids between fibres. The forces providing the bonding of the structure are, in this case, from surface tension. Lyne and Gallay [14], found that as water is removed from the web it became stronger. They attributed the increase in strength to the increase of the air-water interface caused by the removal of moisture. Nevertheless, they found that a further decrease in moisture content (approximately 35%), can actually reduce the amount of air-water interface making the web weaker. At higher consistencies (>90%) the structure can be defined as paper. As was previously stated, an increase in consistency reduces the surface tension bonding forces. However, in the case 2 of paper, strength increases because the nature of the bonds changes. In paper, distance between fibre surfaces is below 0.2 mm. This allows hydrogen bonding, stronger in nature than surface tension, to act as the connecting mechanism between fibres. Hydrogen bonding is what provides the strength of paper and makes it unique from other nonwovens. For example, in fibre glass, fibres do not form hydrogen bonds and need the aid of a bonding substance or glue. 1.2 Fibre Morphology and Fibre Properties If a piece of paper, or a pulp suspension was observed through a microscope, two different things may be perceived. The observations will depend on the pulping process. There are basically two types of pulping processes; mechanical pulping and chemical (or kraft) pulping. The main difference between the two is that lignin, a substance present in the tree to hold the fibres together, is removed in the kraft process. Kraft pulp fibres tend to have a cylindrical shape with a hollow lumen through its centerline. A dark field micrograph of a kraft pulp fibre is shown in Figure 1.3. As the fibres are conformed into the paper, the hollow cylinders shape into flattened ribbons because of the high flexibility. On the other hand, mechanical pulp fibres present in a suspension may seem more rectangular than cylindrical, yet still having a lumen. This is because of the presence of lignin in mechanical pulp fibres. This compound fills the voids between fibres in the tree making the cross-section rectangular. Furthermore, mechanical fibres being more rigid because of the lignin, remain rectangular after drying. They will only collapse into ribbons if the flexibility is increased through refining or the cell wall is thin enough. Figure 1.4 shows the cylindrical and rectangular fibre models. Figure 1.3: Micrograph (obtained through light microscopy) of a kraft pulp fibre. The fibre is in a water solution. At this state the fibre is swollen and has a cylindrical shape. The lumen is irregular in size and changing along the length of the fibre. 3 Figure 1.4: Fibre geometric models. Kraft pulp fibres in a suspension are cylindrical (left figure). Mechanical pulp fibres tend to remain rectangular (right figure). However, after drying, both can be modeled as rectangular because of the high collapse index (f ) of kraft pulp fibres due to their low lignin content. The collapse index (f ) characterizes the degree in which fibres collapse after drying. It can be measured using confocal laser scanning microscopy (Jang et al. [7]). It is defined as: f = 1 LA LAo (1.2) where LA and LAo are the cross-sectional areas of the lumen, before and after drying respectively. When observing the fibres along their length, it is realized that fibres are not straight but are curled sometimes having kinks. A depiction of a curled fibre is shown in Figure 1.5. A way to characterize the amount of curl of a fibre is by the dimensionless curl index (CI) defined as: CI = S C f 1 (1.3) Where S is the curved length and C f is the chord. Another important fibre property is coarseness (w). It is defined as the ratio of the fibre mass to its length: w = M f L f [ mg m ] (1.4) The properties described are for a single fibre. Since in a pulp suspension there is a combination 4 Figure 1.5: Depiction of a curled fibre. The pulping process causes the fibre to curl affecting the paper properties. The amount of curl is characterized by the curl index,CI ( SC f 1). of fibres of different properties, distributions and average values are used. Most of these properties can be measured quickly using a Fibre Quality Analyzer (FQA). However, FQA measurements are for fibres in a suspension. These properties will change when fibres are dried. Accurate, yet time consuming methods exist to determine these properties in the final state of paper (Chinga-Carrasco [1]). 1.3 A Brief Description of the Paper-making Process A description of the paper-making process is shown in Figure 1.6. The pulp suspension is created by separating the fibres in the wood chips be a mechanical (mechanical pulp) or chemical (kraft pulp) process. After that, paper is made by draining water from the suspension. This is done through a screen or a wire. On a paper machine the suspension flows through a headbox aligning the fibres to some extent followed by a dispersion from a rectangular jet to a moving wire. The result is an isotropic structure due to the preferred fibre orientation (machine direction). In the case of standard handsheets, the procedure is different. Water is drained vertically in a batch process and the fiber orientation is random. However, in both cases, the position of the fibre centers on the wire are not randomly distributed. This is due to a phenomena known as flocculation. This is the clustering of fibers caused by electrostatic and possibly other weak fiber attractive forces. It causes an uneven distribution of fibers in the forming wire and is a source of paper weakness. As water is drained by suction and/or gravity the consistency increases to 20% and the force providing the connectivity of the web is mostly surface tension. 5 Figure 1.6: Flow diagram of the paper-making process (commercial and standard handsheet). After the fibres are separated from the wood chips to form a suspension, water is taken from the mixture by a combination of gravity, suction, pressure and heat. 6 The process to follow in paper-making is press-drying. In the case of the paper machine, this is done by sets of rolls that apply pressure as the web goes through the gap between the rolls. In handsheet-making, the process is done through a vertical press. Two circular discs remove water by compressing the paper. In the press-drying process, there is further dewatering. The fibers slide and pack to a denser structure. The consistency is increased up to 50% depending on the pressure and residence time. At this stage, it is believed that fibers remain similar to hollow cylinders and that water is still contained inside the lumen. The final step in the paper-making process is drying. This is done by heated rolls in the paper machine. In contrast, standard handsheets are left to air-dry at a controlled temperature and humid- ity. After drying, consistency reaches 93% as there is still water contained within the fiber cell wall. At this stage, if the structure is observed in a microscope, many of the fibers will appear similar to flattened ribbons instead of hollow cylinders. Fibres have the tendency to collapse by dewatering and by surface tension forces. The tendency to collapse is a function of the fibre thickness and fibre flexibility (Paavilainen [20]). A similar phenomena is observed when drinking liquid from a very thin wall straw. The low pressure, in this case caused by our lungs, collapses the straw. However, when we stop drinking, the straw recovers to its original shape. In the case of fibers, as the gap between the inner surfaces reduces, hydrogen bonding and other weak bonds (Van der Waals) pre- vent it from recovering. The outer surface tension acts in a similar way bringing different fibers together followed by hydrogen bonding. As it was stated in the first section of the chapter, the hydrogen bonding is present only when the distances among the surfaces is less than 0.2 mm. At consistencies lower than 90%, water surrounds the fibers and the lumen, preventing these types of bonds from forming. Similar to the analogy of the collapsing straw, some fibers may be too rigid to collapse and prevent hydrogen bonding. Thinner fibres collapse easier than thicker fibers. The composition of the fibers also affects collapsibility. Mechanical pulp fibers tend to collapse less because the lignin in the lamellas of the fibers increase rigidity. Kraft pulp, on the other hand, is highly collapsible because of the low lignin content. Another source of collapsibility increase is refining. In this process, fibres are delaminated and fibrillated by the action of rotating discs that compress the fibres as the suspension flows through the discs. The result of the forming process described is a structure that is mostly two-dimensional. This is because fibers have high aspect ratios and are more likely to lay flat on the mat. Microscopy measurements show that fibres are aligned between 1 to 3 fibres thick along the plane of the mat. Difference between the structure and the two-dimensional approximation is apparent as we increase the basis weight and as the shape of the fibers deviates from the ribbon-like structures of highly collapsible fibers. 7 1.4 Fibre and Paper Properties It seems reasonable to say that the properties of paper should depend on the properties of the fibres that constitutes the paper and on the process used to consolidate the structure. The following are common paper properties used in the pulp and paper industry: 1.4.1 Property Distributions A tree is composed of fibres of different shapes and sizes. The seasons affect the tree growth such that during the winter, fibres grow smaller in length and cross-section than they do during the spring. This is to accommodate the different rates of liquid transport. Furthermore, pulp mills will use wood chips from different species. Fibre properties are known to change among species and by the pulping process used by the mill (USDA [23]). Softwood trees, such as pines, tend to have longer fibres than hardwood trees like birch or maple. When all these fibres are mixed together in the pulping process, the result is a pulp suspension that has a distribution of size properties (length, width and fibre wall thickness). These properties are known to affect the paper properties. The natural question that arises is; what statistical value is a representative quantity to predict the paper behavior? An arithmetic mean of the property distribution does not correlate strongly with the properties. However, length-weighted average, show better correlation. This is because of the strong weight dependency and the fact that the fibre mass for the most part is linearly proportional to its length. We should now introduce the length-weighted average length (LWL): LWL= åL f 2i åL fi (1.5) where L fi is the length of the i th fibre in the distribution. Similarly for all other properties g , such as width and wall thickness: LWg = å gi L fi åL fi (1.6) where gi is the fibre property of the i th fibre in the distribution. In the case of Probability Density Functions (PDF), we can compute the mean and variance in the following way: m = Z ¥ ¥ f (x) x dx (1.7) 8 s2 = Z ¥ ¥ f (x) (xm)2 dx (1.8) LWL can be calculated by applying equation 1.7 to the length-weighted distribution ef (x) de- fined as: ef (x) = 1R ¥ ¥ f (x) x dx f (x) x (1.9) The unity-weighted distribution can be calculated from the inverse transformation of the length- weighted distribution: f (x) = 1R ¥ ¥ ef (x) x dx ef (x) x (1.10) The previous equations (1.7 to 1.10) can be used for discrete probability in bins by changing f (x)dx to P(x j) and R ¥ ¥ to å N j=1. P(x j) is the probability of bin j. 1.4.2 Basis Weight It is the ratio of the dry fibre mass to the flat area that the fibres cover and commonly reported in grams per square meter. BW = Total Dry Fibre Mass Area [ g m2 ] (1.11) 1.4.3 Paper Strength The paper strength is usually characterized by the tensile index (TI). Other measures of strength are the tensile energy absorbed (TEA) and the breaking length (T ). The tensile index is defined as follows: TI = F d BW [ N m g ] (1.12) where F is the breaking force, d is the strip width, and BW is the basis weight in [ gm2 ]. The testing is performed in a horizontal tester that holds the paper strip with pneumatic clamps. The tester records the force as a standard specimen is pulled apart. A standard specimen has a length of 10 cm and a width of 15 mm. TI is divided by BW , since increasing the total fibre mass will increase the strength. In most cases, this division allows comparison of strength of papers with different BW . However, the comparison is not always valid. IAnson et al. [6] has found that for a given furnish, TI changes 9 with BW , and has an optimal value where it reaches a maximum. 1.4.4 Fibre Strength The fibre strength is usually measured through the Zero Span Tensile Index (ZI). The testing apparatus is very similar the one used for measuring TI. However, for ZI the initial distance between the clamps is very small (< 1 mm). Thus, when the strip is pulled apart, the fibres break instead of the bonds. The value is obtained using the same equation for TI (equation 1.12) except that by reducing the span, the recorded breaking force is usually higher. ZI is also affected by the clamping force of the sample. If the pressure is too low, the strip will slide recording a low value. If the pressure is too high, damage will be done to the fibres such as to reduce the required breaking force. Thus, the value that should be reported is the maximum obtained by changing the clamping pressure. 1.4.5 Paper Thickness and Bulk Calliper or paper thickness is measured with a micrometer applying a very low pressure (< 5psi). The measurement is done at different points of the sample and an average value is recorded. The bulk of the paper refers to the apparent density of paper. It is calculated as the ratio of the calliper to the basis weight of the sample: Bulk = Paper Thickness BW [ cm3 g ] (1.13) 1.4.6 Relative Bonded Area As was previously mentioned, fibres are joined to others fibres in the network through hydrogen bonding. When the paper is forced in tension, usually the bonds will start breaking before the actual fibres. However, it depends on how the strength of the fibre compares to the strength of the bond. Additionally during breakage, the stress can redistribute in such a way that even though bonds started breaking, the fibres will break as well. A way to characterize the degree of bonding is through the relative bonded area (RBA) defined by equation 1.14. RBA= Total Bonded Sur f ace Total Available Sur f ace f or Bonding (1.14) 10 1.5 Brief Summary of Determined Models and Studies on Paper and Fibre Properties For many years, researchers have tried to relate the fibre properties to the paper properties. Un- fortunately, the final product is highly dependent on the complex forming process. Therefore, an applicable model that perfectly relates the paper to the fibres and the process has not been fully developed. Nevertheless, there have been great contributions. A solution to the problem seems closer due to the growing capabilities of computers and simulations. 1.5.1 Experimental Predictions O’neil et al. [19] studied hardwood fibres and experimentally found that tensile index was linearly proportional to length. They did not find any dependency with other dimensional fibre properties. The explanation of the increase in tensile index was attributed to the increase in probability of the fibre-to-fibre bonds along the length of the fibre. However, the length variation in their study, was less than 200 mm. Clark [2] studied mixtures made from sulphite pulp and viscose fibres and determined through regression that the tensile index and bulk behaved as follows: TI w0:6 L0:5f (1.15) Bulk w0:15 L0:05f (1.16) where w is the fibre coarseness and L f is the length-weighted fibre length (LWL). 1.5.2 The Page Equation for Strength Page [21] developed a semi-empirical formula to describe the tensile strength of paper. It is known as the Page equation and is shown as equation 1.17. 1 T = 9 8 Z + 12 Ac r f g b Pf L f (RBA) (1.17) where: T is the breaking length Z is the zero span breaking length Ac is the fibre cross-sectional area 11 g is the acceleration of gravity r f is fibre density Pf is the perimeter of the fibre L f is fibre length b is bond strength RBA is the relative bonded area The Page equation can be expressed in terms of the tensile index by multiplying the breaking length by the acceleration of gravity g. This would lead to equation 1.18 TI = 8 3 b L f RBA 3 b L f RBA+ 4 Ac r f ZI Pf ZI (1.18) where TI and ZI are the tensile index and zero span tensile index respectively. Unfortunately, the use of this equation is limited due to the fact that some of the variables are difficult to measure. Additionally, the variables are interdependent in a complex way. For constant bond force (b), and constant zero span tensile index (ZI), this is an equation of the following form: TI = a x b x+ c y ZI (1.19) where x= (L f RBA) and y= AcPf . Under the assumption that zero span tensile index (ZI) is independent of length, the page equa- tion shows that paper strength increases linearly for small values of length and at a slower rate at higher lengths. It grows with length from zero reaching asymptotically a fraction of the zero span tensile index(ab = 8 9 ). This holds true under the additional assumption that RBA is independent of length. A prediction of the behavior of tensile index with coarseness and width is difficult to obtain from equation 1.18 because RBA is itself a function of the fibre coarseness, width, and collapsibility. 1.5.3 Statistical Models and Computer Simulations If several assumptions and simplifications are made, then the fibre network complexity can be reduced to the point where we can explain the behavior and run simulations to construct a network. More assumptions translates into more limitations in the model. A network can be considered to be composed of a determined distribution of fibres. They 12 can be collapsible, flexible and hollow tubes with cross-section varying along the chord of the fibre. A three-dimensional (3D) simulation can then be performed. However, This would take a considerable amount of effort, not only in the code development but also in computer processing time. The simplest of the models would be to consider fibres, rectangles of a constant shape and where all the fibres would have the same shape. The justification of the simplification in shape is the tendency of fibres to collapse into a ribbon-like structures. A two-dimensional (2D) model may seem unrealistic. However, it can predict valuable information on the paper structure. Additional to the complications arising from the shape and distributions of the fibres, the way in which the network is constructed should also be considered. Deposition of the fibres strongly depends on the flow characteristics and the paper-making process. In a paper-making machine, there is a strong fibre alignment in the machine direction caused by the headbox jet, whereas in a standard handsheet-maker the fibre direction is random because of the vertical drainage. Never- theless, in both cases there is flocculation effects in which the fibre center positions are correlated. If the suspension is very dilute and the formation process is fast enough, the flocculation effects can be minimized. Under all these simplifications, a fibre network can be considered to be a two- dimensional random network. Because of the paper’s random nature, it is not surprising that there is considerable amount of work done by mathematicians and statisticians. Kallmes, Benier and Perez have published several technical papers relating fiber geometry and coverage of two-dimensional (2D) networks to statistical values that may be linked to the paper mechanical properties (Kallmes and Corte [10], Kallmes and Bernier [9], Kallmes et al. [11]). Their theory became known as the KBP theory. Related work on the two-dimensional (2D) paper structure can also be found in the book of M.Deng and Dodson [16]. In order to understand the modeling of paper as a 2D random fibre network, we must introduce some concepts from the literature that define the random network. Coverage Coverage (c), as shown in equation 1.20, is defined as the ratio of area covered by the fibres pro- jected into the sample sheet, to the area of the sample sheet. It is important to notice that coverage includes the total area of fibres, not only the projected area. Hence, a coverage value greater than unity is possible. c= Fibre Area Sample Sheet Area (1.20) 13 The basis weight, BW , and the coverage, c are easily related: for a sample made from fibres of the same thickness, t f : BW = Fibre Area t f (1 f ) r f Sample Sheet Area (1.21) where f represents the voids fraction of the fibre and r f the fibre density. A common practice is to assume a fibre density of 1200 kgm3 instead of 1500 kg m3 , and a value of zero for the voids fraction. Equation 1.21 can be written in terms of the coverage using equation 1.20 as follows: BW = c t f (1 f ) r f c= BW t f (1 f ) r f (1.22) It should be noted that f and t f are functions of the fibre cross-section in the wet state (w and Do) and the fibre flexibility. Decreasing coarseness, w , will result in reducing fibre thickness, t f , (if fibres are collapsible) or increasing the voids fraction, f (if fibres do not collapse). Thus, regardless of fibre collapsibility, reducing the fibre coarseness will increase the coverage. One important relation shown by the KBP theory and also by M.Deng and Dodson [16] is that RBA of a random 2D network depends only on the average coverage (ĉ). They showed that it can be approximated by the following equation: RBA= 1 (1 e ĉ) ĉ (1.23) Another important parameter of a 2D network is the number of fibre-to-fibre contacts or crossings. The expected density of crossings, m , for low coverage is given by M.Deng and Dodson [16]: m = eĉ p ĉ2 4 D2 ( 1 2 + ĉ 3 ) (1.24) where ĉ is the average coverage and D is the fibre width. In three dimensions (3D), Gates and Westcott [4] using a orthogonal fibre model, estimated the 14 bonded area in terms of the fibre coverage and flexibility. One of the most common 3D models mentioned in literature, was developed by Nilsen et al. [17] and coding of a similar model is explained by Conceio et al. [3]. It consider fibres to have a flexible parameter. Fibres are deposited vertically and as the fibre encounters the network, they bend following a flexing rule. 1.6 Research Purpose It is believed that, if it is possible to determine the effect of fibre geometry to the paper properties, the fibre network can be optimized. Paper properties can then be enhanced by fractionation. The objective of the work presented in this thesis is to determine the effects of fibre geometry on the paper strength and bulk. The methods used to determine this relation will be experimental and by computer generated simulations. Predicting the properties of paper experimentally is proven to be a very difficult task. The difficulties arise from the numerous fibre properties (dimension and shape, flexibility, collapsibility, bond strength, pulp freeness) and process variables (flocculation, pressure, jet and wire speeds, etc.) that affect the final paper product. Other studies have determined fibre geometry effects by using different pulps. However different types of pulp differ in other properties besides shape, such as flexibility and chemistry, that can affect bond strength and collapsibility and therefore affect the results. This may be the reason why the literature is not consistent as to the effects of the fibre geometry to the paper properties. Furthermore, the closest empirical formula (Page equation) is an expression that requires knowledge of properties of the final paper. Experimentally we will study a single pulp furnish so that the differences in behavior can only be attributed to the difference in geometry and process variables. Experimental results are supported by computer simulations. There have been simulations done by other researchers, using a Lagrangian approach to solve the motion of fibres constructing the fibre mat. However, fibre-to-fibre effects are usually neglected because of the complications in the differential equations. A complete simulation considering the fibre to fibre effects has great com- putational costs. In the simulation studies presented here, the goal was to determine the statistics of fibre networks without involving the mechanics of the process. 15 Chapter 2 Predicting Bulk and Strength of Paper In this chapter, the relations of bulk and strength to fibre geometry (length, width and coarseness) and press-drying pressure are determined. This is done by producing handsheets from different fibre geometry at different pressures. The pulp used was softwood kraft pulp. This type of pulp provided a broader range of lengths compared to hardwood pulps. Kraft pulp was preferred because of the lower presence of fines. 2.1 Experimental Procedure Three hundred (300) grams of softwood kraft pulp was fractionated using a Bauer Mcnett Classifier (BMC) into five fractions (M14, M28, M48, M100 and M200). The coarseness, length, width and curl index of all fractions and feed pulp was measured using a Fibre Quality Analyzer (FQA). Only the M14, M28 and M48 fractions together with the feed were used for making handsheets as the mass of the collected samples of the finer screens (M100 and M200) was not enough to make handsheets. The handsheets were made at different press-drying pressures while others were air- dried. Some of the air-dried handsheets were cut using a guillotine into small squares and repulped. The cutting was done in a 1 mm and 3 mm grid pattern to reduce the fibre length. By doing this different values of fibre length were obtained at a constant coarseness. Length, width and curl index of the cuts was measured. From the repulped cuts, handsheets were made at different press- drying pressures. The residence time in the drying press was five (5) minutes for all samples. After handsheets were air-dried at controlled humidity and temperature, they were tested for bulk and tensile index. Figure 2.1 shows a diagram of how the pulp was fractionated and cut to prepare the handsheets. 16 Figure 2.1: Flow diagram of experimental procedure. The original feed pulp was separated into five fractions. Handsheets were made at different pressures from some of the frac- tions (Feed, M14,M28 andM48). Air-dried samples were cut and repulped to create new handsheets. Tensile index and bulk for all handsheets was recorded and fibre properties were measured for all fractions. 2.2 Assumptions and Limitations The results obtained from the experimental procedure described here are valid under certain as- sumptions. By air-drying the handsheets there should be no drying effects (Jentzen [8]) since the drying is done under minimum load. Therefore, after cutting, fibres have the same properties as before being cut except for the decrease in length. Fibre size distributions are assumed to be sufficiently narrow so that the average values can be representative of the distribution. The pulp used is a mixture of softwood species. It is assumed that during the fractionation 17 process, the species are not separated as to affect the data by a difference in flexibility, or any other property that may differ from one species to the other. Thus, the observed differences account exclusively for changes in the geometry. Aminimum of three (3) handsheets were prepared for the evaluation of handsheet properties. The errors for bulk and strength were calculated, assuming T-Distribution and propagation of error. Table 2.1: Size and w measurements of different pulps obtained from FQA. Errors are from standard deviation between FQA measurements. The table shows consistent correlations among all properties (increasing D and w with L f ) except for CI. Fibre cutting or frac- tionation did not affectCI significantly. Mesh Size Length, LWL (mm) Fibre Width (mm) Curl Index Average Coarseness(mgm ) Original Pulp 2:40:2 281 0:170:01 0:200:01 -(coarse cut) 1:500:08 - 0:1050:007 - -(Fine cut) 1:20:1 - 0:100:01 - 14 Mesh 3:660:03 31:50:3 0:1860:006 0:2570:004 -(coarse cut) 2:330:07 - 0:1240:004 - -(Fine cut) 1:510:02 - 0:0920:004 - 28 Mesh 2:540:04 29:20:7 0:170:02 0:2030:003 (coarse cut) 2:010:03 - 0:1570:004 - (Fine cut) 1:330:03 - 0:1170:004 - 48 Mesh 1:650:03 26:70:8 0:1390:009 0:1840:004 (coarse cut) 1:40:1 - 0:140:07 - (Fine cut) 1:20:3 - 0:120:06 - 100 Mesh 0:9410:008 24:60:4 0:1080:001 0:1600:007 200 Mesh 0:570:02 22:80:3 0:1270:001 0:1360:002 2.3 The Fiber Distributions In this section, the measurements of fibre size (L f , D andCI) and w are described. 18 2.3.1 Length-Weighted Average Values Using BMC and FQA, size distributions and coarseness of the different fractions and their respec- tive cuts were measured. Table 2.1 shows length-weighted averages of these measurements. The errors were obtained from analysis of different FQA measurements of the same sample. This error is not to be confused with the deviation or narrowness of each distribution. 2.3.2 Property Distributions The data from FQA measurements yields the length-weighted length distributions presented in Figure 2.2. The recorded distributions appear to be Gaussian. The literature (M.Deng and Dodson [16]) shows that fibre length distributions from the tree tend to be log-normal. However, perhaps from the pulping process, length-weighted distributions shift to Gaussian. Furthermore, Fig 2.2 shows that as the mesh turns finer, the distribution is narrower. This is commonly observed in multi-stage fractionation. It is due to the effect of the probability screening. Even though the grid openings are always larger than the fiber width, it is easier for shorter fibers to pass through the screen and long fibers to be rejected. Gooding and Olson [5] analyzed the fractionation process and proposed a plug flow model. They found that for a given type of screen aperture geometry (slots, holes, etc), the passage of long fibres cannot be reduced without reducing the passage of short fibres as well. Thus, by making the downstream distribution narrower the upstream distribution grows wider by the short fibres retained. The effect is consistent through each stage of the screening process. Another important observation of Fig 2.2 is, that even though Mesh 28 has a length- weighted average length (LWL) very close to that of the feed pulp, the distributions are different. The feed pulp distribution, being a mixture of all the mesh sizes, is very wide compared to the Mesh 28 distribution. Figure 2.3 displays the length distribution for Mesh 28 and the two cuts. The figure shows that for the coarse cut, LWL reduces but the standard deviation increases when compared with M28. In the fine cut, the distribution is narrower. This shows that for the coarse cut there are two effects occurring simultaneously. In the cutting, there is a reduction in length of some of the fibres but the rest of the fibres remain unchanged. Thus, it is as if short fibres were added to the uncut distribution. This is why the coarse cut is wider. In contrast, in the cutting process for the fine cut, length for all fibres is reduced. This creates a narrow distribution. Since the size of the cuts was the same for all fractions, the cutting effect is lost as we decrease the size of the mesh opening. This is shown in Figure 2.4. 19 Figure 2.2: Length distributions (length-weighted) of different fibre fractions obtained from a BMC fractionation. The mean of each distribution shifts as the mesh is finer. All appear to have Gaussian behavior, with a consistent reduction in standard deviation as the mesh fineness increases. Figure 2.3: Length distribution (length-weighted) of a BMC fraction and its respective cuts. Fibre cutting was effective by reducing LWL of each fraction, shown here as shifts in the distributions. Furthermore, they remain Gaussian without considerably increasing standard deviation. 20 Figure 2.4: Length-weighted pulp distributions for different cuts (in sequential order: M14, M28, M48 and Feed Pulp). As the mesh is finer the cutting effect is lost and the resulting distribution is wider. Figure 2.5 shows the scatter plot of length and width of each fiber obtained from FQA data. It can be seen that, even though the range is broad, there is a relation among width and length for each fibre. Madani [15], separating the scatter values into bins and averaging the values in each bin, proposed a relationship (implying length being linear to the fibre aspect ratio) of the form: D= L f ao L f +a1 (2.1) 21 Figure 2.5: Scatter plot of width D and length from fibres for different distributions (Feed Pulp, M14, M28, M48, M100 and M200 in sequential order). Although there is a broad range of width for fibres at certain length, there is a higher concentration at a given value. This is an indication that the width and length are correlated. This is regardless of the fibre fraction. Thus, there is a relation among fibre properties, that may be unaffected by the fractionation process. 22 2.4 The Natural Links in Fiber Properties It has been found that fiber properties are linked (Madani [15],Karenlampi and Suur-Hamari [12]). Length, coarseness and width are usually correlated. Longer fibers tend to have a higher coarseness and a larger width. Figure 2.6 shows the relation found between coarseness and length. FQA does not measure coarseness of each fibre but total length of a distribution. By knowing the dry weight of the FQA sample, a mean coarseness can be calculated as the ratio of sample dry weight to total length. The trend for coarseness was found to be linear (2.2) for the range of values in the study with a coefficient of determination R2 = 0:986. This is in agreement with the findings of Karenlampi and Suur-Hamari [12]. The reason of the behavior is unknown, yet it must be related to the tree biology, genetics, and climate. Another source for the correlation, may be due to the fact that coarser fibres are stronger and have a greater resistance to fibre cutting in the pulping process. Figure 2.6: Relation between fibre coarseness and length. Coarseness was obtained from the total measured length and the dry weight of each distribution sample. Thus, it is a mean value since there is a length distribution. A linear fit provided the highest coefficient of determination (R2 = 0:986). The standard errors, shown as bars in the plot, were obtained from standard deviation among FQA measurements. w = 0:036 mg m mm L f +0:12 mg m (2.2) 23 At a 95% confidence level, the coefficients of the linear fit (equation 2.2) are 0:036 0:006 and 0:120:01. Figure 2.5 shows that for all fractions there is a wide range of width for a given length. How- ever, there is a higher concentration of fibres of a certain width. The width value with higher concentration is increasing as the fibres get longer. Following the procedure used by Madani [15], the raw data of all fractions (Figure 2.5) was separated into 0:2 mm length bins. An average width for each bin was calculated. This calculation was done for all BMC fractions, generating curves of length against width for each. The curves are shown in Figure 2.7a. Regardless of the frac- tion, the generated curves are statistically the same. Thus, this correlation is independent of the fractionation process. Furthermore, the same observation proves that the BMC does not fractionate according to width which is consistent with the findings of Gooding and Olson [5]. In an attempt to reproduce the results of Madani [15], a fit of the form of equation 2.1 was proposed. Nevertheless, for the range of values in our studies and the pulp used, a linear fit (equation 2.3) better predicted the behavior as shown by the coefficients of determination. Figure 2.7: Relation between fibre length and width (D). Constructed from raw data shown in Figure 2.5. The scatter points represent arithmetic averages of width among fibres in 0:2mm bins. Left figure (a) shows width values for each BMC fraction. Right figure (b) shows the combined data of all BMC fractions. Mean fibre widths reported by FQA for all fractions have also been included. Fit 1 (equation 2.1 proposed by Madani [15]). Fit 2 (linear fit). The best fit was found to be linear (equation 2.2). 24 D= 2:2 mm mm L f +24:0mm (2.3) At a 95% confidence level, the coefficients of the linear fit (equation 2.3) are 2:20:3 and 241. As is the case of the determined relationship for fibre coarseness, it is difficult to establish the cause of correlation between length and width without involving tree genetics and biology. A study of this type is beyond the scope of this research, but a review of genetic and geographical studies on the fibre correlations is provided by Via et al. [24]. 2.5 Multi-variable Regression From the different BMC fractions and cuts, handsheets were tested for bulk and strength. The results are shown in Tables A.1 to A.4 of the Appendix. The most simple models for the tensile index and bulk, in terms of the fibre geometry and press-drying pressure, are the following: Bulk =CB L gL f D gD wgw Pgp (2.4) TI =CT L bL f D bD wbw Pbp (2.5) Applying logarithm to both sides of the previous expressions, the values of the unknown coef- ficients are obtained through a simple multi-variate linear regression. The regression yields the values shown in Table 2.2. The values of CB and CT depend on the units used for the other vari- ables of the models (length, width, coarseness and pressure). They are a function of the elasticity of the fibers and possibly other variables like curl, kappa number and fibre fibrillation. Table 2.2 shows considerable error in the determined coefficients. In the case of bulk, gL is negligible. Thus, the model shown in equation 2.4 can be improved by removing length as a variable. In the case of tensile index, the data shown in Appendix (Tables A.3 and A.1), reflects great difference in results for Mesh 28 and the feed pulp. Both pulps have similar LWL, yet the strength values are very different. The observed differences are due to the difference in the shape of the distribution. It is known that fines affect the strength of the sheet yet have negligible effect on LWL (Seth [22]). Thus, LWL is not a sufficient measure of length distribution. This means that the model shown in equation 2.5 can be improved by adding the distribution standard deviation (s ) shown in Table 2.3. The coefficients eliminating length as a variable for bulk and adding s as a variable for TI, 25 Table 2.2: Coefficients for proposed bulk and tensile index models (equations 2.4 and 2.5) determined through linear regression. Unknown Value (95% confidence level) CB 0:10:7 gL 0:030:09 gD 22 gw 0:10:8 gp 0:200:03 CT 8e104e13 bL 0:700:09 bD 72 bw 0:40:8 bp 0:310:02 Table 2.3: Standard deviation of length-weighted pulp distributions. Obtained using equation 1.8 and the length-weighted distribution shown in equation 1.9 Pulp Sigma(s ) (mm) Feed Pulp 1:30 Feed Pulp Coarse Cut 0:95 Feed Pulp Fine Cut 0:74 14 Mesh 1:03 14 Mesh Coarse Cut 1:12 14 Mesh Fine Cut 0:75 28 Mesh 0:84 28 Mesh Coarse Cut 0:90 28 Mesh Fine Cut 0:67 48 Mesh 0:56 48 Mesh Coarse Cut 0:72 48 Mesh Fine Cut 0:61 26 Table 2.4: Improved coefficients for bulk and tensile index models (2.4 and 2.5) as obtained from regression. Unknown Value (95% confidence level) CB 0:080:4 gD 1:30:6 gw 0:10:3 gp 0:200:01 CT 5e78e9 bL 0:40:1 bs 0:60:1 bD 52 bw 0:60:7 bp 0:310:02 are those shown in Table 2.4. By comparison of the errors in the coefficients shown in Table 2.2 and 2.4, a clear improvement of the model can be identified. To assess the improvement quality of the fit of the models to the data, the predicted values of the dependent variables with the measured values are plotted. Figure 2.8 shows the quality of the bulk model fit, considering length and s and without considering length at all. Clearly, length has negligible effect on bulk. Figure 2.9 shows the quality of the tensile index model fit under the same mentioned conditions. In this case, the tensile index model is improved by adding length and s as variables. Comparing the figures, not only is there an improvement in the deviation of the tensile index fit, but now the data from feed pulp distributes evenly around the diagonal line. 27 Figure 2.8: Effect of reducing the number of variables in the bulk model. The measured values of bulk were obtained from calliper and grammage measurements. Left figure (a) considers length and s , right figure (b) shows the improved model by eliminating length. Figure 2.9: Effect of incorporating the standard deviation (s ) to the tensile index model. The measured values of tensile index were obtained from a tensile tester. Left figure (a) does not consider s as variable, right figure (b) shows the improved model by adding s . 28 Thus, bulk and tensile index models (equations 2.4 and 2.5) can be rewritten as equations 2.6 and 2.7. Bulk = 0:1 D1:3 w0:1 P0:2 (2.6) T:I:= 5e7 L0:4f s 0:6 D5 w0:6 P0:3 (2.7) These correlations predicted the bulk and tensile index within 20% and 25% respectively of the actual values. From the determined relations, one of the findings is that the bulk and TI are not intrinsically related. Even for the same furnish, the pressure coefficient of the bulk is not the inverse of the pressure coefficient for tensile index. It is true that generally a reduction of one of the two leads to an increase of the other, since the power coefficients of equations 2.6 and 2.7 have opposite signs. Nevertheless, this is not always the case. A change in length may considerably affect strength and have little effect on bulk of the sheet. Furthermore, equations 2.6 and 2.7 show that bulk is much less sensitive than tensile index to being affected by fibre properties. 2.6 Zero Span Tensile Index Zero span tensile index (ZI) was measured for each sample. This property is dependent on the clamping pressure of the tester. ZI was measured at two clamping pressures (60 psi and 70 psi). The recorded value was the maximum of the two. The values are shown in Tables A.1 to A.4 of the Appendix. A regression performed on the zero span tensile index data yields a very poor correlation as shown in Figure 2.10. Due to the large randomness in recorded values, little is to be said about how the fiber geometry affects the zero span strength. The zero span tensile index is more sensitive to the clamping pressure and grammage than to fibre geometry (L f , D or w). In regards of the press-drying pressure, for most of the pulps, it was found that increasing pressure beyond a certain value (60-100 psi) caused a reduction of ZI (data shown in the Appendix). This was not the case for the tensile index. It is believed that increasing press-drying pressure, while it reduces the voids and increases the degree of bonding, still damages the fibres when it is excessive. 29 Figure 2.10: The predicted value obtained through a linear regression of the zero span tensile index (ZI) against the measured values. The result of the large randomness of the measured ZI is a fit of very poor quality. 2.7 Combining the Natural Links with the Regression Given the relations found in the regression analysis (equations 2.6 and 2.7) and the relation of coarseness and width to length (equations 2.2 and 2.3) for the studied pulp, both bulk and TI can be related to two (2) single variables (L f and P): Bulk = 0:1 (2:2L f +24:0)1:3 (0:036L f +0:12)0:1 P0:2 (2.8) TI = 5e7 L0:4f s 0:6 (2:2L f +24:0)5 (0:036L f +0:12)0:6 P0:3 (2.9) This can be done provided there is no mechanical modification of fibers (e.g. cutting or refining) since modified fibres will not follow correlations given in equations 2.2 and 2.3. Equations 2.8 and 2.9 can be represented as surface plots. Figure 2.11 shows the behavior of bulk (equation 2.8 ), while Figure 2.12 corresponds to tensile index (equation 2.9) for two (2) different s values. From the surface plots, in the case of bulk, we see monotonic behavior with respect to length and pressure. It increases with length and decreases with pressure. In the case of tensile index, it 30 Figure 2.11: Surface plot showing bulk for studied pulp in terms of pressure and length. Ob- tained by combining the model from regression (equation 2.6) and determined corre- lations between fibre properties (equations 2.2 and 2.3). increases with pressure but shows an optimal length. The optimal value for studied pulp furnish depends on the standard deviation. It is around 2.6 mm for a s of 10%. The two values of s are presented to show that the location of the optimal fraction will also depend on the shape of the distribution (or the fines content). The linear regression shows no optimal value when the variables are considered to be independent. The regression predicts a stronger paper, the longer the fibre. However, longer fibres are coarser and have a larger width. This acts against the strength. The result of the opposing effects is the existence of an optimal length value. 31 Figure 2.12: Surface plot showing tensile index for studied pulp in terms of pressure and length at two s values ((a) 0.1 mm and (b) 10%). Obtained by combining the model from regression (equation 2.6) and determined correlations between fibre properties (equations 2.2 and 2.3). 32 Chapter 3 Computer Simulations In this chapter, two-dimensional (2D) and three-dimensional (3D) simulations performed using Matlab are shown. The developed codes are shown in the Appendix B. 3.1 Paper as a 2D Random Network In this section, the paper structure is modeled as a 2D constant size (length and width) random fibre network. The random networks were created using Matlab imaging toolbox. Each fiber is placed as a 2.25 Megapixel (1500p X 1500p) binary image. Resolution of the image was set to a 5 mm square. Thus, the fibre dimensions were converted to pixels accordingly. An average fibre width of 30 mm is represented as a rectangle of 9 pixels in width. A pixel value of 1 represents fibre area and 0 represents a void. By adding the images we can form the random 2D structure. To avoid edge effects, once the network is constructed, the image was cropped to 80%. The resulting image will not be binary. Thus, there is a pixel value distribution along the plane of the sheet which will be referred as P(x;y). Figure 3.1 shows random networks created using this method for two different coverage values. The colors of the pixels are assigned according to P(x;y). Thus, we are not only looking at the structure of the simulated paper, but also the P(x;y) distribution. P(x;y) gives the total number of fibres stacked on top of each other at a given location. This means that information on the number of fibre-to-fibre crossings and the resulting bonded area can be obtained for that location. This is shown in Table 3.1. 33 Figure 3.1: Simulated fibre networks at different coverage values. Left figure (a) c = 0:24 and right figure (b) c= 12. For both cases L f = 1 mm, D= 30 mm and the sample area is a 5 mm X 5 mm square Table 3.1: Description of the pixel value distribution P(x;y). The bonding area factor is mul- tiplied by the pixel area (Ap) to calculate the bonded area at the (x;y) location. Pixel Value (P(x;y)) Region Description Bonding Area Factor 0 Void Regions 0 1 Fibres 0 2 Single Fibre-Fibre Crossing Region 2 3 Double Fibre-Fibre Crossing Region 4 N (N-1) Fibre-Fibre Crossing Region 2(N1) 3.1.1 Assumptions and Limitations of the 2D Model Usually, fibres are horizontally aligned within a thickness scale of 2-3 fibre thickness (M.Deng and Dodson [16]). Consequently, it is reasonable to analyze the paper structure as a 2D formation. In fibre deposition on the paper machine, fibre centers distribution and fibre orientation are influenced by several phenomena. These include; flocculation of the suspension before the wire, local basis weight unevenness caused by the flow through the forming fabric and preferred orienta- tion caused by the headbox jet. There are simulation studies that take these effects into account. An example is the work done by N.Provatas et al. [18]. The study predicts the statistics of disordered networks and devises a simulation parameter to characterize flocculation. For standard handsheet making there is no preferred orientation, thus the distributions tend to be more random. The simu- 34 lations shown here focus on standard handsheets, thus each fibre center position and orientation is assumed to be random. As fibres are stacked, a problem arises when more than two fibres are overlapping. If we consider fibres to be rigid, then as the fibres are stacked, they cannot be in contact with fibres in different layers other than the immediately adjacent layer. However, flexible fibres can bend so that they come into contact with fibres in different layers. Figure 3.2 shows the differences between the two cases. This can lead to issues when counting the fibre crossings within a randomly generated 2D network. For this study, the counting of the crossings is done in the following manner; If there is a double fibre crossing, it counts as two crossings, a triple crossing as three crossings and so on. This is equivalent to considering fibres as infinitely flexible. Thus, in this analysis, fibres are considered to be infinitely flexible and stacked in such, that there are no internal voids. The voids can only be on one of the sides of the sheet and not between the fibres of different layers. In reality fibres are not infinitely flexible. The fibres ability to bend will depend on the loads (induced by pressing and drying) and on the second moment of area of the fibres cross-section calculated from the horizontal axis. However, the model described here is the limiting case. Thus, we expect to overestimate the number of fibre-to-fibre crossings (Nc) and relative bonded area (RBA) and to obtain a lower paper thickness than that of the real paper. Figure 3.2: Difference in fibre contact between rigid and flexible fibres. a) Flexible fibre. b) Rigid Fibre. Finally, this initial study, only considers constant fibre property distribution. In order to obtain a representative statistical value, five simulations were conducted in the same conditions and the values were averaged. The error bars were obtained from standard deviation among simulations. 35 3.1.2 Paper Thickness and the Pixel Value distribution The pixel value P(x;y) represents how many fibres are stacked on top of the other at that location. Thus it is also the local coverage at that location. M.Deng and Dodson [16] have shown that local coverage behaves as a Poisson distribution. Consequently, the pixel value distribution is also a Poisson distribution as shown in Figure 3.3. Figure 3.3: The Poisson distribution of the pixel value (P(x;y)) obtained from simulations for two different coverage values. L f = 1 mm, D= 25 mm, N f = 3000 (a) and 10000 (b). Coverage (c) can be calculated from the pixel value distribution to show that is exactly equal to the mean pixel value\P(x;y): c= å N i=1Pi Ap SampleArea = åNi=1Pi N =\P(x;y) (3.1) Furthermore, if P(x;y) is multiplied by fibre thickness (t f ), the result is the network thickness dis- tribution tN(x;y). The mean network thicknessctN can be obtained by integrating over the sample area in the following way: ctN = R R tN(x;y)dxdySampleArea (3.2) in a discrete form: ctN = åNi=1 tNi ApSampleArea (3.3) 36 where N is the number of pixels and Ap is the pixel area. Equation 3.3 can be rewritten as follows: ctN = t f åNi=1Pi ApSampleArea = t f c (3.4) Thus, ctN is simply the product of fibre thickness t f and coverage c. Finally using equation 1.22,ctN can be expressed as: ctN = BW (1 f ) r f (3.5) where BW is the basis weight, f is the voids fraction in the dried state and r f is the fibre density. If our assumptions are valid, this means that bulk for the network is equal to fibre bulk in the dried state ( 1(1 f ) r f ). Thus, in order to make a bulkier paper, uncollapsible bulkier fibres must be used. It should also be noted that the voids fraction in the dried state is a function of the fibre collapsibility (f ) and the cross-section geometry (w and Do) in the wet state. Consequently, when the mean thickness is used to calculate the bulk of the network, the bulk only depends on the fibre cross-section (w and Do) and not its length. When measuring paper thickness experimentally, the measuring cylinder that drops down to the paper has a large diameter compared to the dimensions of the fibres. Thus, a more representa- tive paper thickness comparable to that measured experimentally would be the maximum network thickness , tN;MAX . There is no analytical expression for the maximum pixel value (PMAX = tN;MAX t f ). However, the behavior obtained from simulations is shown in Figure 3.4. Since the sample area is a 5 mm square, this would be the thickness measured by a device with an indenter of the same size. Figure 3.4 shows that similar to\P(x;y), PMAX depends only on coverage and not the fibre dimensions (L f and D). By performing a linear regression, PMAX was shown to behave as: PMAX 10 c0:6 (3.6) Similar to the way equation 3.5 was obtained, expression 3.6 can be rewritten to show the following: Bulk 10 c0:4 (1 f ) r f (3.7) 37 Figure 3.4: The effect of coverage (c) on the maximum pixel value (PMAX ) for different lengths (L f ) and widths (D). Values obtained from 2D simulations. PMAX changes only with c and not with L f or D. or: Bulk 10 t f 0:4 BW 0:4 (1 f )0:6 r f 0:6 (3.8) This shows, that when tN;MAX is used, bulk decreases as BW increases. Furthermore, f and t f are functions of the fibre cross-section (w and Do) in the wet state. Thus, regardless of the thickness used to calculate bulk (mean or maximum), the network bulk only depends on the fibre cross- section (w and Do) and not length. 3.1.3 The Degree of Bonding The strength of the network can be related to the number of fibre-to-fibre crossings and the relative bonded area. 38 Number of Crossings When a paper strip is loaded, the tension is distributed through the fibres and the bonds between the fibres. Usually prior to fracture, if each fibre is isolated, the transmitted load is considerably smaller than that required for fibre breakage. Therefore, bonds are usually the first to break. Consequently, it is important to maximize the number of fibre crossings (Nc) so that the load distributes through more bonds even if the total bonded area is the same. By increasing Nc, the probability of finding a crossing (along a fracture line) in which the fibres and the load are aligned is higher. Figure 3.5 shows the effect on Nc by the addition of fibres to the structure for different lengths and widths. The figure shows that Nc is only a function of width (D) and coverage (c). Nc increases with increasing c and decreases with increasing D. Thus, paper made from fibres of smaller widths have more fibre crossings at a given basis weight. Length (L f ) has no effect on Nc. From the figure, it can be seen that even for high coverage values, the number of crossings is only a function of coverage and fibre width. The dependence on fibre length is weak because fibres, regardless of length, have high aspect ratios. Figure 3.5: The response of coverage(c) on total number of fibre-to-fibre crossings (Nc). Nc increases as fibre width (D) is reduced and coverage (c) is increased. 39 At a constant coverage, the number of fibres (N f ) decreases as length increases. Since Nc is independent of length, at a constant coverage, increasing length will increase the number of cross- ings per fibre (NcN f ). Thus, networks made from longer fibres will be stronger because of the better alignment of the loads in the direction of the fibres and the fact that if a particular bond breaks, the fibre is still held by other bonds that can activate. Furthermore, reducing fibre coarseness will result in a higher coverage. Thus, the total number of crossings will increase. By increasing the total number of crossings, the load is transmitted through more bearing points, and the strength of the paper network increases. The effect of width, at first, is to decrease the number of crossings. However, reducing width (at the same coarseness) may increase or decrease the coverage depend- ing on the flexibility of fibres. Thus we cannot precisely determine the effect of width without considering the network flexibility. Deng and Dodson (1.24) give an expression (1.24) for the density of the number of crossings for small coverage values. This formula shows that Nc is independent of L f . Figure 3.6 shows a comparison of the formula with values from simulation. As shown, the formula is a better predic- tion as fibre width increases. Furthermore, the formula underpredicts Nc as c increases. Figure 3.6: Comparison of number of fibre crossings obtained from simulations with statis- tical formula (1.24) by M.Deng and Dodson [16]. The scattered values were obtained from simulations and the solid lines represent equation 1.24 at different widths. 40 The Relative Bonded Area (RBA) The relative bonded area (RBA) can be obtained by the pixel values. The total fibre area within the sample can be calculated through the following: Total Fibre Area= N å i=1 2PixelValuei PixelArea (3.9) Equation 3.10 is used to calculate the total unbonded area. Total Unbonded Area= 2Number o f pixels(PixelValue> 0)PixelArea (3.10) Therefore, RBA can be obtained by expression 3.11. RBA= Total Fibre AreaTotal Unbonded Area Total Fibre Area RBA= 1 Total Unbonded Area Total Fibre Area RBA= 1 Number o f pixels(PixelValue> 0) åNi=1PixelValuei (3.11) Figure 3.7 shows the RBA vs. fibre coverage. From the figure, it is observed that RBA is not a function of length or width of the dried state, but only a function of basis weight or coverage. The blue line in Figure 3.7 represents the theoretical value shown in Chapter 1, in equation 1.23. This shows that length does not affect the total degree of bonding, it only changes the efficiency in the way that load is transmitted (by affecting NcN f ). On the other hand, decreasing coarseness not only increases the total number of crossings, but also increases the amount of bonded area. This is because decreasing coarseness will increase coverage as shown in equation 1.22. Thus, reducing the coarseness should increase the strength of the paper network. Size Distributions Since width affects the number of crossings Nc, the effect of the width distribution was studied. This was done by creating random networks made from fibres with a normal width distribution. The simulated structure of a normal width distribution is shown in Figure 3.8. Several simulations 41 Figure 3.7: The relative bonded area (RBA) against coverage (c) of simulated 2D networks for different lengths (L f ) and widths (D). RBA is only a function of c and unaffected by L f or D. The area along the thickness of the fibre (2 t f L f ) was not considered for RBA calculation. The solid line represents the statistical approximation provided by M.Deng and Dodson [16] as equation 1.23 at different width standard deviation (sD) and mean width (bD) were performed. The number of crossings was recorded at different coverage values. The results of these simulations are shown in Figure 3.9. The figure shows that only the mean width (bD) affects the number of fibre crossings. The width standard deviation (sD) has no effect on the number of crossings. 3.2 Two-Dimensional (2D) Deposition From the 2D analysis, under the validity of assumptions (low coverage or high press-drying pres- sure), it is observed the length has no considerable effect on; the bulk, the number of Crossings (Nc) or RBA. Fibre width (D), shows an effect on Nc, increasing as D reduces. However, experimental 42 Figure 3.8: Simulated network constructed from a Gaussian width fibre distribution. c = 0:969, bD= 45 mm and sD = 40 mm. evidence shows length effects on both bulk and tensile strength. In the initial assumptions, the packing density changes due to geometry, were not considered. Thus length and width may have effects on the packing density, specially at low press-drying pressures, such as to change thickness and number of crossings. A 2D random deposition brings light into the answer of the previous question. In the previous sections of this chapter, random network simulations were performed to study the plane of the sheet (x,y). In this section, the focus is the vertical plane (x,z) along the thickness of the sheet. Matlab was used to develop a code to simulate the fibre deposition in this vertical plane. The deposition model considers fibres to be rigid beams of large friction factor and bonding force being deposited by the action of gravity. Figure 3.10 shows a deposition generated using 43 Figure 3.9: Effect of width standard deviation (sD) of a normal distribution on the number of crossings Nc. The values of sD are in mm. sD shows no effect on Nc. Filled data points are for an average width (bD) of 30 mm and data points with no fill are for bD= 45 mm. this code. The total fibre area is 74 103 mm2 or around 100 fibres, the fibre length is 1 mm and thickness is 15 mm. Figure 3.10: Simulated two-dimensional (2D) deposition. Fibres are considered to be rigid beams with high friction factor and high bonding force. Total fibre area is 74 103 mm2, L f = 1 mm and t f = 15 mm. Resolution is 25 mm2 per pixel Each simulated beam was deposited vertically from a random x position and at a horizontal orientation. Once the beam encountered the formed mat, if equilibrium was not reached, it rotated until finding another contact point. At this point another fibre is deposited a so on. It must be said 44 that in the deposition model, the fibres are considered rigid with a large friction factor and bonding strength. By performing several simulations, a statistical average can be obtained to measure the effect of length on the mean and maximum height of the deposition. This is shown in Figures 3.11 and 3.12. The figures show that for both heights (mean and maximum), there is a slight increase with length. This is an indication that packing density reduces as fibre length is increased. In the next section, this question will be answered by performing a 3D simulation. Figure 3.11: Mean height of a rigid fibre 2D deposition against total area of deposited fibres. Data obtained from simulations. Each point represents an average of three simulations repeated under the same conditions. The fibre thickness t f is 15 mm. 45 Figure 3.12: Maximum height of a rigid fibre 2D deposition against total area of deposited fibres. Data obtained from simulations. Each point represents an average of three simulations repeated under the same conditions. The fibre thickness t f is 15 mm. 3.3 The Three-Dimensional (3D) Structure of Paper In the previous sections we have studied the in plane and vertical statistics of randomly deposited fibres by 2D simulations. In this section the the focus is on the statistics of the 3D random network. The developed 3D Matlab code is an extension of the 2D codes. A 3D network can be simulated as a 3D image. In this case, a pixel value of unity in the 3D array represents fibre material, whereas a pixel value of zero is empty space. In contrast with the 2D random network code, each fibre is a 3D array representing a volume. The network is created by translation and rotation of the volume. As the fibre is being deposited, intersection is detected when pixel values greater than unity appear. Once the intersection is detected, the fibre will rotate until coming into contact with another point of the mat. Both of these operations can be broken down into translations and rotations of 2D images. Figure 3.13 show a 3D deposition of 300 fibres in a 5 mm square using the developed code (L f = 1 mm and GSM 10). Figure 3.14 shows a top view of the simulated network, where the height (in pixels) of the network is represented by a color scale. 3.3.1 Assumptions and Limitations The starting point for deposition is picked randomly at a random orientation within the horizontal plane. The intersection, as the fibre is being deposited, is obtained by adding the pixel value and 46 Figure 3.13: Simulated random fibre network generated in Matlab. The model used considers fibres as rigid beams of high friction factor and bonding force. The fibre dimensions are L f = 1 mm, D= 25 mm and t f = 15 mm. The figure shows 300 fibres deposited in a 5 mm square (10 GSM). 10 20 30 40 50 60 Figure 3.14: Top view image of the simulated random fibre network shown in Figure 3.13. Color has been added to show the height distribution in pixels. Resolution is 5:8 107 mm3 per voxel. detecting a value greater than unity. Once the intersecting height is determined, the fibre is rotated up to point where it will come into contact with another point of the existing network. The sign of the rotation angle is chosen by comparing the free lengths on each side of the contact pivot point. After that, the fibre becomes permanently attached to the network. Thus, our assumptions are the following: 47 Fibres are considered to be rigid beams of a known width, fibre thickness and length. The size and shape for all fibres are the same. Property distributions are not considered. The friction factor of the fibres is large enough so that there is no sliding. The bond force is large enough such that there is no sudden restructuring due to the fibres own weight. That is, the deposition of each fibre does not change the previously formed mat. The deposition of each fibre is done at a random location and orientation. Thus flocculation or other forming effects are not considered. 3.3.2 The Three-Dimensional (3D) Algorithm A fibre center position and orientation angle is selected randomly. The fibre is deposited vertically until there is contact with the network (Initially it is a flat surface). After contact, a 2D height map similar to Figure 3.14, is cropped to the size of the fibre. The cropped image is used to predict the possible fibre rotation. The free lengths of the fibre are compared and side of rotation is determined. The rotation angle is calculated from the 2D height map and the pivot point. The fibre is rotated and the first step is repeated until a certain number of fibres is reached. For every iteration, the total volume of the fibre network is obtained by summing the total pixels. The height is calculated by creating a 2D array of the maximum height of the 3D structure. From the 2D array, the maximum and mean height are estimated. The developed code, differs from those explained in the introduction (Conceio et al. [3]). In the code used here, rotation has been incorporated and fibres are considered rigid. It is believed, that this is a more realistic simulation for high rigidity and low pressure cases. The other code would overpredict the bulk by preventing rotation. 48 3.3.3 Length Effect By adding fibres to the mat, it is obvious that the mean height increases. This is because adding fibres is equivalent to increasing the basis weight of the mat. If the mean height of the mat is plotted against the fibre total volume, the effect of length at a constant fibre total volume can be observed. Figure 3.15 shows the effect of length at a given fibre total volume. It can be inferred from the figure, that length has the effect of increasing the mean height of the mat after 1 mm3 of fibre volume. This volume is equivalent to a basis weight of approximately 50 GSM. 3.3.4 Degree of Bonding of a Three-Dimensional (3D) Network Unfortunately, obtaining the number of crossings is a difficult task in 3D. Nevertheless, some information can be inferred. The 3D simulations (voids are considered), show the mean height increasing with length. In contrast, 2D simulations (voids are not considered), show a constant mean height. Thus, the amount of voids, or in other words packing density, is dependent on length and pressure. At low pressure, length has the effect of producing a higher bulk. However, at higher pressure, length has no effect. Figure 3.15: Mean height of a simulated 3D random network (ctN) against fibre volume (N f L f D t f ). Fibre thickness and width are t f = 15 mm and D = 15 mm. Each point represents the average of three simulations and error bars are from standard deviation. The figure shows an increase inctN with L f and fibre volume. 49 Figure 3.16: Maximum height of simulated 3D random network (tN;MAX ) against fibre vol- ume (N f L f D t f ). t f = 15 mm and D= 15 mm. Each point represents the average of three simulations and error bars are from standard deviation. Comparing with Figure 3.15, as more fibres are added the difference between tN;MAX andctN is reduced, yet for low fibre volumes, tN;MAX is in some cases ten foldsctN . 50 Chapter 4 Conclusions and Final Remarks The combination of fibre fractionation and guillotine cutting, proved to be a successful method for the study of the paper structure and the effects of fibre geometric properties. Furthermore, the 2D imaging method used here to simulate fibres, shows to be relatively fast and successful. It has the benefit compared to other methods, that fibre-to-fibre contacts can be easily computed. Both simulation and experiments, show that some of the effects can be attributed to the statistics while others are due to mechanical responsive properties. By relating the fibre properties for the studied pulp, fibre coarseness (w) and width (D) were found to increase linearly with fibre length (L f ). Regardless of the fraction, the correlations were the same. This shows that the Bauer McNett Classifier (BMC), separates the fibres according to their length (L f ) independent of their width (D) or coarseness (w). The results from the performed linear regression, are consistent with the findings in the liter- ature (Clark [2]). It was found that length, when considered as an independent property, has no significant effect on the bulk of the sheet. However, the same is not true for the strength. Length and the standard deviation of the length distribution (s f ), were found to considerably increase ten- sile index (TI µ L f 0:4s f 0:6). Width and coarseness were found to increase bulk (Bulk µ w0:1D1:3) and to decrease strength (TI µ w0:6D5). It should be said that the proposed equation models (equations 2.4 and 2.5) are a first approach for studying the bulk and tensile index behavior. This type of empirical model allows a comparison with the work developed by Clark [2] and O’neil et al. [19]. This type of fit implies that the effect of one of the variables is independent of the other and forces monotonic behavior. Thus, the model needs to be improved. An approach to improve the model is mentioned in the recommendations for future work. When the correlations between fibre geometric properties (L f , D and w) are combined with 51 the empirical relations determined by linear regression, the result is an optimal length fraction for a stronger paper. The linear regression shows no optimal value when the variables are considered to be independent. The regression predicts a stronger paper, the longer the fibre. However, longer fibres are coarser and have larger widths. This acts against the strength. The result of the opposing effects is the existence of an optimal length. This explains why the original feed, being shorter than the longest mesh (M14) yet longer than the shortest mesh (M48), made the strongest paper. Thus, the tensile index can be increased through the use of fibre fractionation. However, the determined fibre property correlations are not universal among all fibre species and types of pulp furnish. The existence of an optimal value on other types of pulps will depend on how the fibre geometric properties are correlated. The performed two-dimensional (2D) simulations provide insight into the observed strength behavior. Increasing fibre length, by increasing the number of fibre crossings per fibre (NcN f ), in- creases strength. This is because by increasing the number of crossings per fibre, since the fibre is mostly straight, will increase load alignment. Furthermore, during elongation if a bond between fi- bres breaks, each fibre has other bonds that can activate to support the load. Increasing coarseness, reduces the network coverage (equation 1.22). This in turn, decreases the degree of bonding (RBA and Nc), making a weaker paper. The effect of width is difficult to obtain from 2D simulations without involving the fibres ability to collapse and bend. Width in the wet state (cylindrical fibre) is different than width in the dried state (rectangular fibre). Increasing width in the wet state (at the same coarseness) will increase or decrease coverage at the dried state depending on collapsibility. Furthermore, 2D simulations provide an overestimate of Nc and RBA. The initial assumptions of the 2D model is that fibres are infinitely flexible. In reality this is not the case. Reducing width will reduce the area moment of inertia allowing the fibres to bend and better approximate the 2D structure. The effect of fines is also difficult to relate with 2D simulations. Fine coarseness is very low, thus the coverage is increased. However, they are also very short. So, why is it that fines enhance strength? We believe that this can be explained by the filling effect of fines rather than by the change in average length or coarseness. We have said that in reality the fibres are not infinitely flexible. Thus, at a particular fibre crossing the fibre does not bend at a 90o angle to come into contact with fibres of another level (see Figure 3.2). Since fines are small they fill in these voids such as bring the structure closer to the 2D approximation. This means that adding fines after a certain threshold has little benefit in terms of enhancing strength. In fact, if fibres were infinitely flexible it is predicted that adding fines would have little to no effect. However, it is recognized that this explanation is simplistic and requires further studying to prove. Simulations also explain the bulk behavior. We studied the two limiting cases of infinitely 52 flexible and infinitely rigid fibre. Two-dimensional simulations show that when fibres are very flexible, length had no effect on bulk. On the other hand, three-dimensional (3D) simulations show that when fibres are very rigid, increasing length increases bulk. This attributes to the notion that paper is for the most part, a 2D structure since the 2D model better predicts the behavior. The stresses induced by drying and fibre flexibility, are sufficient to cause the fibres to assemble into a dense 2D structure. However, the 2D model is still a limiting case and the approximation of paper as a 2D structure is itself dependent on the fibre geometry. 4.1 Recommendations for Future Work The experimental work to follow is to study the application of the determined empirical relations to other types of furnish and to account for the effect of fibre flexibility. This can be done by examining data from other researchers or by repeating the experimental methodology described here. Flexibility can be changed by beating or by chemical treatments. However, fibre cutting with guillotine proved to be a tedious task. Thus for future work, perhaps it is necessary to develop an apparatus to cut the fibres in a controlled manner. Furthermore, the simple empirical models (equations 2.4 and 2.5) can be improved by replacing them with equations that can predict bulk and strength on a broader range of values. The form of equations 2.4 and 2.5 allow a comparison with the work done by other researchers. However, we know that in order to be consistent with the physics, the form of equations 2.4 and 2.5 should be different. These equations assume that the effect of each variable is independent and monotonic. In reality, we know for example that changing the coarseness at a given width changes the fibre collapsibility. A change in coarseness will have a smaller effect on wall thickness the larger the width. Furthermore, it is expected that as pressure is increased the effect of the rest of the variables should be smaller. Thus, a better approach is to preform an analysis of variance and use a non-linear model that accounts for interaction between variables. We predict strong interactions between pressure (P) and the rest of the variables (L f , w and D) as well as between D and w . However, this will also require increasing the data points and varying the properties in ranges that may not be possible. Finally, it has been seen that statistics of random networks explain qualitatively the effects of fibre geometry to bulk and tensile index. However, there is a need to link them directly in a quantitative way. This will allow the development of a wholesome theory to predict the behavior of any pulp. Nevertheless, it will require an analysis of the mechanical response of fibres to; the flow field during deposition, the applied pressure, the forces induced by the drying process, and the 53 tensile force during elongation. Even in a two-dimensional random network predicting breakage is difficult. The stress at a particular bond is not constant but a distribution that depends on the tension of the bonded fibres, the fibre dimensions, and the orientation between fibres. All of which are changing from bond to bond. Thus, computer solving through FEA (Finite Element Analysis) may be useful relating network statistics (Nc, RBA, tN) to the strength. 54 Bibliography [1] G. Chinga-Carrasco. Microscopy and computerized image analysis of wood pulp fibres multiscale structures, volume 3 of 4. Formatex Research Center, 2010. ISBN 9788461461912. ! pages 5 [2] J. D. Clark. Effects of fiber coarseness and length. i. bulk, burst, tear, fold, and tensile tests. TAPPI, 45(8):628, 1962. ! pages 11, 51 [3] E. L. T. Conceio, J. M. R. Curto, R. M. S. Simes, and A. A. T. G. Portugal. Coding a Simulation Model of the 3D Structure of Paper, volume 6026 of Lecture Notes in Computer Science. Second International Symposium, CompIMAGE 2010, Buffalo, NY, USA, May 5-7, 2010. Proceedings, 2010. ISBN 9783642127113. ! pages 15, 48 [4] D. Gates and M. Westcott. Predicting fiber contact in a three dimensional model of paper. Journal of Statistical Physics, 94(1):31, 1999. ! pages 14 [5] R. Gooding and J. Olson. Fractionation in a bauer-mcnett classifier. JPPS, 27(12):423, 2001. ! pages 19, 24 [6] S. IAnson, W. Sampson, and S. Savani. Density dependent influence of grammage on tensile properties of handsheets. JPPS, 34(3):182, 2008. ! pages 9 [7] H. F. Jang, R. Amiri, R. Seth, and A. Karnis. Fiber characterization using confocal microscopycollapse behavior of mechanical pulp fibers. TAPPI, 79(4):203, 1996. ! pages 4 [8] C. Jentzen. The effect of stress applied during drying on some of the properties of individual pulp fibres. TAPPI, 47(7):412, 1964. ! pages 17 [9] O. Kallmes and G. Bernier. The structure of paper:the absolute, relative, and maximum bonded areas of random fiber network. TAPPI, 45(11):867, 1962. ! pages 13 [10] O. Kallmes and H. Corte. The structure of paper:the statistical geometry of an ideal two dimensional fiber network. TAPPI, 43(9):737, 1960. ! pages 13 [11] O. Kallmes, G. Bernier, and M. Perez. The theory of the tensile failure of paper. Paper Technology and Industry, 18(11):222, 1977. ! pages 13 55 [12] P. Karenlampi and H. Suur-Hamari. Paperi Puu, 79(6):404, 1997. ! pages 23 [13] R. Kerekes and C. Schell. Characterization of fibre floculation regimes by a crowding factor. JPPS, 18(1):J32, 1992. ! pages 2 [14] L. Lyne and W. Gallay. Fiber properties and fiber-water relationships in relation to the strength and rheology of wet webs. TAPPI, 37(12):581, 1954. ! pages 2 [15] A. Madani. Fractionation of Particle Suspensions in a Viscoplastic Fluid: Towards a Novel Process. UBC Thesis, 2011. ! pages viii, 21, 23, 24 [16] M.Deng and C. Dodson. Paper An Engineered Stochastic Structure. TAPPI PRESS, 1994. ISBN 0898522838. ! pages ix, 13, 14, 19, 34, 36, 40, 42 [17] N. Nilsen, K. Niskanen, and M. Zabihian. Kcl-pakka: A simulation tool for paper properties. TAPPI Proceedings:1997 Engineering and Papermakers Conference, page 1311, 1997. ! pages 15 [18] N.Provatas, M.Haataja, E.Seppala, S.Majaniemi, M.Alava, and T.Ala-Nissila. Structural properties of disordered fibre networks. Physica A, 239:304, 1997. ! pages 34 [19] P. O’neil, S. Mohiddin, and A. J. Michell. Exploring data for relationships between wood, fibre and paper properties. Appita, 52(5):358, 1999. ! pages 11, 51 [20] L. Paavilainen. Conformability-flexibility and collapsibility-of sulphate pulp fibres. Paperi Ja Puu, 75(9):689, 1993. ! pages 7 [21] D. Page. A theory for the tensile strength of paper. TAPPI, 52(4):674, 1969. ! pages 11 [22] R. Seth. The measurement and significance of fines. TAPPI, 104(2):T47, 2003. ! pages 25 [23] USDA. Density, fibre length, and yields of pulp for various species of wood. United States Department of Agriculture Technical Note, No. 191, 1953. ! pages 8 [24] B. Via, M. Stine, T. Shupe, C.-L. So, and L. Groom. Genetic improvement of fibre length and coarseness based on paper product performance and material variability-review. IAWA Journal, 25(4):401, 2004. ! pages 25 56 Appendix A Experimental Data 57 Ta bl e A .1 :E xp er im en ta lh an ds he et da ta :b ul k, T I an d Z I fo rt he fe ed pu lp an d di ff er en tc ut s M es h Si ze L en gt h (m m ) W id th ( m m ) C ur lI nd ex C oa rs en es s (m g m ) Pr es su re (p si ) B ul k (c m 3 g ) T. I. (N m g ) Z er o Sp an (N m g ) Fe ed Pu lp 2: 4 0: 2 28 1 0: 17 0: 01 0: 20 0: 01 A ir dr ie d (' 5) 3: 6 0: 1 11 :9 0: 5 16 6 - - - - - 20 2: 62 0: 05 17 1 18 4 - - - - - 40 2: 50 0: 02 17 :3 0: 8 18 4 - - - - - 60 2: 44 0: 08 19 1 18 3 - - - - - 65 2: 31 0: 02 19 1 19 8 - - - - - 80 2: 32 0: 05 20 1 15 5 -( co ar se cu t) 1: 50 0: 08 - 0: 10 5 0: 00 7 - 20 2: 47 0: 06 13 :0 0: 9 12 2 - - - - - 40 2: 49 0: 09 12 :9 0: 6 13 2 - - - - - 60 2: 32 0: 07 13 :7 0: 7 14 3 - - - - - 80 2: 3 0: 1 14 1 16 4 -( Fi ne cu t) 1: 2 0: 1 - 0: 10 0: 01 - A ir dr ie d (' 5) 3: 7 0: 1 6: 3 0: 3 18 1 - - - - - 20 2: 72 0: 03 9: 1 0: 4 11 3 - - - - - 40 2: 48 0: 05 10 :9 0: 7 16 6 - - - - - 60 2: 36 0: 08 11 :9 0: 6 14 6 - - - - - 80 2: 29 0: 08 12 1 18 4 58 Ta bl e A .2 :E xp er im en ta lh an ds he et da ta :b ul k, T I an d Z I fo rt he M 14 fr ac tio n an d di ff er en tc ut s M es h Si ze L en gt h (m m ) W id th (m m ) C ur lI nd ex C oa rs en es s (m g m ) Pr es su re (p si ) B ul k (c m 3 g ) T. I. (N m g ) Z er o Sp an (N m g ) 14 M es h 3: 66 0: 03 31 :5 0: 3 0: 18 6 0: 00 6 0: 25 7 0: 00 4 A ir dr ie d (' 5) 4: 9 0: 3 3: 4 0: 2 21 4 - - - - - 20 3: 62 0: 09 6: 3 0: 8 19 5 - - - - - 40 3: 3 0: 2 7 1 27 6 - - - - - 60 3: 0 0: 2 9 1 26 6 - - - - - 80 2: 8 0: 2 10 2 25 4 -( co ar se cu t) 2: 33 0: 07 - 0: 12 4 0: 00 4 - A ir dr ie d (' 5) 4: 7 0: 2 4: 4 0: 4 19 2 - - - - - 20 3: 6 0: 2 7 1 17 6 - - - - - 40 3: 2 0: 2 8 1 12 4 - - - - - 60 3: 0 0: 2 10 1 25 5 - - - - - 80 2: 7 0: 2 11 1 22 5 -( Fi ne cu t) 1: 51 0: 02 - 0: 09 2 0: 00 4 - A ir dr ie d (' 5) 4: 54 0: 02 2: 6 0: 3 15 1 - - - - - 20 3: 49 0: 06 4: 2 0: 5 13 2 - - - - - 40 3: 22 0: 09 5: 2 0: 6 14 4 - - - - - 60 2: 9 0: 2 6 1 19 1 - - - - - 80 2: 7 0: 2 7 1 19 3 59 T a bl e A .3 :E xp er im en ta lh an ds he et da ta :b ul k, T I an d Z I fo rt he M 28 fr ac tio n an d di ff er en tc ut s M es h Si ze L en gt h (m m ) W id th (m m ) C ur lI nd ex C oa rs en es s (m g m ) Pr es su re (p si ) B ul k (c m 3 g ) T. I. (N m g ) Z er o Sp an (N m g ) 28 M es h 2: 54 0: 04 29 :2 0: 7 0: 17 0: 02 0: 20 3 0: 00 3 A ir dr ie d (' 5) 4: 82 0: 09 5: 5 0: 3 9 3 - - - - - 20 3: 2 0: 1 10 1 17 3 - - - - - 40 2: 8 0: 2 13 1 19 5 - - - - - 60 2: 6 0: 1 14 2 20 4 - - - - - 80 2: 37 0: 08 16 1 16 3 - - - - - 10 0 2: 21 0: 08 19 1 14 2 (c oa rs e cu t) 2: 01 0: 03 - 0: 15 7 0: 00 4 - A ir dr ie d (' 5) 4: 35 0: 02 5: 9 0: 5 15 5 - - - - - 20 2: 93 0: 09 10 :1 0: 7 16 3 - - - - - 40 2: 66 0: 09 11 :3 0: 8 19 2 - - - - - 60 2: 40 0: 05 13 :0 0: 7 21 3 - - - - - 80 2: 32 0: 05 13 :3 0: 9 20 1 - - - - - 10 0 2: 21 0: 01 15 1 17 3 (F in e cu t) 1: 33 0: 03 - 0: 11 7 0: 00 4 - A ir dr ie d (' 5) 4: 25 0: 04 3: 5 0: 2 17 4 - - - - - 20 2: 9 0: 1 5: 9 0: 8 17 2 - - - - - 40 2: 6 0: 1 7: 6 0: 8 20 1 - - - - - 60 2: 5 0: 1 8: 0 0: 6 15 2 - - - - - 80 2: 4 0: 2 8: 8 0: 6 16 3 - - - - - 10 0 2: 25 0: 02 10 :1 0: 4 13 3 60 T a bl e A .4 :E xp er im en ta lh an ds he et da ta :b ul k, T I an d Z I fo rt he M 48 fr ac tio n an d di ff er en tc ut s M es h Si ze L en gt h (m m ) W id th (m m ) C ur lI nd ex C oa rs en es s (m g m ) Pr es su re (p si ) B ul k (c m 3 g ) T. I. (N m g ) Z er o Sp an (N m g ) 48 M es h 1: 65 0: 03 26 :7 0: 8 0: 13 9 0: 00 9 0: 18 4 0: 00 4 A ir dr ie d (' 5) 4: 4 0: 1 7: 0 0: 4 9 2 - - - - - 20 2: 79 0: 05 13 :5 0: 6 9 3 - - - - - 40 2: 5 0: 1 15 :1 0: 6 10 4 - - - - - 60 2: 4 0: 1 16 :0 0: 7 12 5 - - - - - 80 2: 2 0: 1 19 :4 0: 4 20 4 (c oa rs e cu t) 1: 4 0: 1 - 0: 14 0: 07 - A ir dr ie d (' 5) 3: 8 0: 1 4: 8 0: 5 16 1 - - - - - 20 2: 71 0: 06 8: 3 0: 5 14 8 - - - - - 40 2: 43 0: 04 10 :1 0: 4 13 2 - - - - - 60 2: 31 0: 05 10 :9 0: 8 12 2 - - - - - 80 2: 17 0: 03 12 1 13 1 - - - - - 10 0 2: 15 0: 01 12 :2 0: 3 17 4 (F in e cu t) 1: 2 0: 3 - 0: 12 0: 06 - A ir dr ie d (' 5) 3: 83 0: 08 7: 4 0: 6 13 2 - - - - - 20 2: 66 0: 05 12 1 18 :0 0: 7 - - - - - 40 2: 3 0: 1 14 1 18 1 - - - - - 60 2: 26 0: 07 14 :6 0: 7 19 2 - - - - - 80 2: 15 0: 06 16 1 17 :8 0: 5 - - - - - 10 0 2: 04 0: 01 16 :5 0: 4 19 4 61 Appendix B Matlab Codes 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- A study on the structure of paper : the links between...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
A study on the structure of paper : the links between paper and fibre properties Fernandez Madrigal, Francisco J. 2013
pdf
Page Metadata
Item Metadata
Title | A study on the structure of paper : the links between paper and fibre properties |
Creator |
Fernandez Madrigal, Francisco J. |
Publisher | University of British Columbia |
Date Issued | 2013 |
Description | The work shown in this thesis, focuses on the effects of fibre geometry (fibre length, width and coarseness) to the paper strength and bulk. It is believed that fibre networks, for a given furnish, can be optimized through fractionation to create a stronger and bulkier paper. A first set of experiments was performed to determine the correlations among fibre properties of a softwood kraft pulp furnish. The fibre properties were measured using a Fibre Quality Analyzer (FQA). Coarseness (ω) and width (D) were found to increase linearly with length (Lf ). These correlations are thought to be influenced by the tree species and the pulping process. A second set of experiments was aimed to determine empirical expressions of bulk and tensile index (TI) in terms of the fibre geometry distribution (ω, Lf and D) and the press-drying pressure (P). Bulk and TI was measured for handsheets made at different pressures from different size distributions. These distributions were created by a combination of fractionation on a Bauer McNett Classifier (BMC) and fibre cutting. The determined relations show agreement with experimental work from other researchers. In order to provide insight into the causes of such behavior, simulations of two-dimensional (2D) and three-dimensional (3D) random networks were performed. Matlab was used to write the simulation codes. Mechanical response of fibres (affected by the fibre geometry) to the forces induced by drying or pressure is not considered directly in the simulations. However, 2D simulation models are a good representation of high press-drying conditions and high fibre flexibility, whereas 3D models are better predicting air-dried (low pressure) paper made from fibres of high rigidity. Geometric statistics of random networks explain some of the experimental observations. The statistical outputs of the 2D simulations were; network coverage (c), network thickness (τN), number of fibre crossings (NC) and the relative bonded area (RBA). In the case of 3D simulations, only τN was determined. The 2D and 3D outputs were measured for different fibre size distributions (Lf, D and ω). |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Date Available | 2013-04-17 |
Provider | Vancouver : University of British Columbia Library |
Rights | Attribution-NonCommercial-NoDerivatives 4.0 International |
DOI | 10.14288/1.0073696 |
URI | http://hdl.handle.net/2429/44227 |
Degree |
Master of Applied Science - MASc |
Program |
Mechanical Engineering |
Affiliation |
Applied Science, Faculty of Mechanical Engineering, Department of |
Degree Grantor | University of British Columbia |
GraduationDate | 2013-05 |
Campus |
UBCV |
Scholarly Level | Graduate |
Rights URI | http://creativecommons.org/licenses/by-nc-nd/4.0/ |
AggregatedSourceRepository | DSpace |
Download
- Media
- 24-ubc_2013_spring_fernandez_francisco.pdf [ 2.85MB ]
- Metadata
- JSON: 24-1.0073696.json
- JSON-LD: 24-1.0073696-ld.json
- RDF/XML (Pretty): 24-1.0073696-rdf.xml
- RDF/JSON: 24-1.0073696-rdf.json
- Turtle: 24-1.0073696-turtle.txt
- N-Triples: 24-1.0073696-rdf-ntriples.txt
- Original Record: 24-1.0073696-source.json
- Full Text
- 24-1.0073696-fulltext.txt
- Citation
- 24-1.0073696.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
https://iiif.library.ubc.ca/presentation/dsp.24.1-0073696/manifest