Journal of Korea TAPPI. 30 April 2019. 26-39
https://doi.org/10.7584/JKTAPPI.2019.04.51.2.26

ABSTRACT


MAIN

1. Introduction

Refining is an energy intensive process in the pulp and paper industry. The spindle power of a disc refiner is transferred to the refining area through the bars of the plate and the morphology of the pulp fibers is changed by the mechanical friction between the bar and fibers or between the fibers themselves to achieve specific pulp properties. The refining intensity of the refiner is greatly influenced by the bar structure of the plate, which plays an important role in changing the pulp properties and fiber morphology.

Many studies have conducted quantitative analyses of the refining process to find an accurate intensity to measure refining based on energy. Brecht1) proposed specific refining energy (SRE) to represent the net energy consumption of a unit mass of absolutely dry pulp, which is widely used in mills; however, this measure is not the most effective means to measure the intensity of the refining process because it just represents the amount of energy consumed and does not reflect how it is consumed.

Wultsch et al.2) took into account the number of bars, their average length, and the rotation speed of the disc to propose the “spezifische Kantenbelastung”. Brecht et al.3) presented a new concept, namely, cut edge length (CEL), and made an improvement on Wultsch’s work with the specific edge load (SEL), which is widely used today. The SEL is defined as the energy consumed per unit bar cut edge length, and the CEL is considered an effective parameter when using the SEL to describe the refining process. Lumiainen,4,5) Meltzer et al.,6,7) Musselman et al.8) accounted for many bar parameters, such as the bar width and angle, which are not considered in the SEL, and defined the special surface load (SSL), modified specific edge load (MEL) and modified specific surface load (MSSL) based on the SEL. Roux et al.9) modified the SEL by accounting for additional geometric parameters and identified the reference specific load (RSL) by replacing the CEL with the reference length Lref. The four intensities mentioned above are all based on the SEL and generate a further extension of the SEL, known as SEL-based refining intensities by Liu et al.10)

Some researchers have proposed a new way to qualify the refining process by considering the energy consumed from a single impact to the fiber through the unit bar length. Lewis and Danforth11) proposed that the SRE of the refining process is the product of the impact number N and the impact intensity I of the plate on the fiber, and determined that the formulas for the dimensionless quantity N and the fiber intensity I which are inaccurate. Leider and Nissan12,13) deduced more accurate expressions for N and I. The variable I is estimated by measuring the energy required to stretch the fiber to its elastic limit, and N is calculated by estimating the residence time of the fiber in the refining area. Kerekes14) used the C factor to characterize the ability of a refiner to apply an impact to the fibers and derived formulas for N and I based on many refining process parameters, such as plate parameters and gap and fiber properties. This approach is the most rigorous and accurate of all the refining intensities. However, the calculation of C factor in steady state is normally complicated and made more complex by the dynamic nature of the actual refining process. Therefore, the C factor can only predict an average value and cannot be used for real-time monitoring of process parameters, which limits its practical application. Andritz proposed the MagnusTM theory, which describes the refining process using two parameters: the compression index and the number of treatments.15-17) According to this theory, the refining effects of the bar on the fiber are classified by cutting, squeezing and brushing. This theory shows that using the reasonable groove, bar and gap parameters can improve the quality and efficiency of refining.

The refining process can be classified as being high or low intensity. As shown in Fig. 1, low intensity refining can further increase the specific surface area of the fiber by promoting the internal fibrillation of the fiber and increase the strength of the paper. However, external fibrillation of the fibers is promoted and greater cutting occurs during high intensity refining, meaning that the strength and length of the fibers decrease as the printing performance of the paper improves.

https://cdn.apub.kr/journalsite/sites/ktappi/2019-051-02/N0460510203/images/JKTAPPI_2019_v51n2_26_f001.jpg
Fig. 1.

Effect of different refining intensities on fibers.

Kerekes18) concluded that some definition of the refining intensity is the machine intensity and others provide fiber intensity. For machine intensity, such as SEL, SSL, etc., the result of the refining process depends on the strength of them. When considering the SEL, a value of 3-5 J/m is suitable for pulp with long fibers, and strong cutting effect on the fibers occurs during high intensity refining. An acceptable range for low intensity refining of short fibers is 0.5-2 J/m, which results in a better fibrillation effect.19) However, for fiber intensity, the SRE is the product of N and I when the fibers undergo continuous refining, meaning that different refining effects can be obtained for different N and I values that still yield the same SRE,20,21) as shown in Fig. 2.

https://cdn.apub.kr/journalsite/sites/ktappi/2019-051-02/N0460510203/images/JKTAPPI_2019_v51n2_26_f002.jpg
Fig. 2.

The refining results for different N and I values.

Previous research on refining intensity has mainly focused on the machine intensity, as represented by the SEL. However, some studies have examined the refining intensity from the fiber perspective, which requires further study. The focus of our research is to seek a reasonable refining intensity to predict or measure the changing properties of pulp fibers, such as the beating degree, average fiber length, etc. The characterization parameters for the refining intensities and BI-based refining intensities are proposed in this article, and the characterization performances of six refining intensities are obtained using data from a previous study.22,23)

2. Materials and Methods

2.1 Characterization parameters of refining intensity

The machine intensity is broadly defined in Eq. 1 as:

[1]
I=PnetnA

where Pnet represents the net power, n is the rotation rotor speed, and A is a relative parameter called the characterization parameter of refining intensity.

The intensity of the refining process is very complex and can be influenced by many factors, including the design variables of the disc refiner, control parameters and pulp properties. Most of the measures of the refining intensity use plate parameters, such as bar width, bar height, groove width, etc., as starting points. By analyzing the refining intensities, studies of characterization parameters can be carried out through two aspects: the characterization parameter based on the SEL and the characterization parameter based on the BI considering the status of the plate.

2.1.1 Characterization parameters based on CEL

Brecht et al.4) proposed the SEL theory to characterize the effect of the plate on the fibers, as described by Eq. 2:

[2]
SEL=PnetnZsZnL

where Pnet is the net power consumed during refining, n is the rotation speed, Zs and Zr are the number of rotor and stator bars, respectively, and L is the average length of the bars.

The variable ZsZrL is the cut edge length (CEL) and is a characterization parameter of the SEL via a comparison to the definition of A in Eq. 1. The calculation of the CEL is defined in Eq. 320) as:

[3]
CEL=1iZrrZsrΔr

where r is the radius of the plate.

Only the effect of the bar edge is taken into account in the SEL, and other plate parameters, such as bar width, groove width, disc radius, etc., are not taken into consideration. To make the SEL more complete, Lumiainen5) proposed the new SSL, MEL, and MSSL intensities based on plate parameters, such as the bar width and bar interaction angle γ. By similar methods, the characterization parameters of these intensities can be obtained, as shown in Table 1. The characterization parameters of the SSL, MEL and MSSL are extensions of the CEL and are still based on the CEL. The characterization parameters of the SSL and MSSL are the same for the same rotor and stator parameters.

Table 1.

Characterization parameters of the SEL and its extended intensities

IntensitySELSSLMELMSSL
Characterization
parameters
CELCELBs+Br2cosγ2CEL2tanγBB+GCELBcosγ2

Bs: the bar width of the stator, Br: the bar width of the rotor, γ: the interaction angle, B: the average bar width, and G: the average groove width.

2.1.2 Characterization parameters based on bar interaction

The bars and grooves of rotor and stator that intertwine with each other exert shearing and squeezing effects on the pulp or fibers during refining and result in changes to the fiber morphology. Therefore, it is particularly important to analyze the refining intensity in terms of the bar interactions. Elahimher et al.22,23) studied the bar interaction process for the first time at the mosaic point and deduced the parameters of the bar interactions.

Assuming that the total area of the disc is the working area AT and the bar interaction quadrilateral is defined as ASC, as shown in Fig. 3, then the interaction point N is calculated by Eq. 4:

https://cdn.apub.kr/journalsite/sites/ktappi/2019-051-02/N0460510203/images/JKTAPPI_2019_v51n2_26_f003.jpg
Fig. 3.

Bar interaction quadrilateral.

[4]
N=ATASC

Using the definition of bar interaction quadrilateral, the sum of the perimeter and area of all bar interaction quadrilaterals are defined as the bar interaction length and area, which can be obtained Eq. 5 and Eq. 6, respectively:

[5]
A=B2sinγN
[6]
B=4BsinγN

Elahimher believes that the entire area of the disc is effectively refining area when refining, but another results obtained by reasonable analyses and simplification of the actual bar interaction process present in this paper, which is elaborated on below.

2.2 Modeling the intensity based on the bar interactions

2.2.1 BIL and BIA

During the actual refining process, the pulp fed into the gap between the rotor and stator is continuously squeezed and sheared under the action of the edge and surface of the bars, which changes the fiber morphology. Therefore, defining the bar surface area (BIA) as the refining area would be more accurate than using the total disc area.

The bar interaction parallelogram (BIP) formed by the interaction of a single bar on the rotor (or stator) and the adjacent groove on stator (or rotor) is defined as being composed of two interaction points: a bar interaction point and a bar and groove interaction point. The sum of the adjacent two sides of the bar interaction point in the BIP is defined as the single bar interaction length CBIL. The single bar interaction area AB is the area of a single bar interaction point. These variables can be obtained by Eq. 7 and Eq. 8, respectively:

[7]
CBIL=2Bsinγ
[8]
AB=B2sinγ

By calculating the effective refining area Ae, the total number of bar interaction points NT of the refining area can be calculated by Eq. 9:

[9]
NT=NSN

The BIL and BIA of the entire refining area can be expressed by Eq. 10 and Eq. 11, respectively:

[10]
BIL=2BNSHψ,βX1β,ψ,Ri,R02B+G2
[11]
BIA=B2NSHψ,βX1β,ψ,Ri,R02B+G2

where H and X1 are functions describing the bar profile of refining plate. The detailed derivation process associated with BIA and BIL is shown in Appendix A.

2.2.2 The refining intensity based on the bar interactions

The refining intensity characterizes the degree of the impact on the pulp or fibers during refining. The SEL defines the bar edge as the effective part of the plate, and the bar edge and width are concerned by the SSL, which makes a improvement to the SEL. Many SEL-based refining intensities have been proposed considering the patterns of bars and grooves to modify the SEL. However, confirming the bar interaction parameters is more important for refining through the analysis of previous studies and the refining process. The BIA and BIL that characterize the bar interactions can be obtained in this way, and the refining intensities are established considering the general formula for refining intensity.

The refining intensity based on the bar interaction area BIA can be expressed as follows:

[12]
I=PnetnBIA

and the intensity that considers the BIL is defined as:

[13]
I=PnetnBIL

2.3 Experimental

To analyze the characterization of the different refining intensities for the refining process, the experimental data mentioned in Appendix C of reference 23, which were obtained from a 16-inch single disc refiner with a variable speed drive that can be adjust to 800, 1000 and 1,200 rpm using 3% CTMP pulp (Quesnel River Pulp), were used in this paper, and the initial beating degree of the pulp was 32.8 SR°, the average fiber length was 1.9 mm. All the trials were conducted at the same conditions, namely, the stock temperature was held constant at 60℃, and the flow rate was 250 rpm. A low consistency refiner was used in these experiments and is shown in Fig. 4.

https://cdn.apub.kr/journalsite/sites/ktappi/2019-051-02/N0460510203/images/JKTAPPI_2019_v51n2_26_f004.jpg
Fig. 4.

Illustration of the low consistency refiner used for experiments.

The power was varied over practical plate gap ranges with a minimum of 15 samples collected per trial to provide repeatability and high resolution. The trials were run at three (3) different rotational speeds (as mentioned above) and for a wide range of plate geometries. The fiber lengths and the beating degrees of the pulp were measured by a Fiber Quality Analysis and CSF tester. The refiner is equipped with flow meters, power meters, and plate actuation and pressure sensors to easily collect the power, gap and the other required parameters.

The basic parameters of the four plates used in the experiments are shown in Table 2. All the plates were treated with an abrasive material in a stock suspension before the experiments were conducted to ensure parallelism.

Table 2.

The basic parameters of 4 plates used in the experiments

NumberRi (mm)R0 (mm)B (mm)G (mm)α (°)φ (°)β (°)CEL (km/rev)
1111.125203.23.24.815157.50.89
2111.125203.223.615157.51.80
3111.125203.21.63.215157.52.45
4111.125203.212.415157.55.00

3. Results and Discussion

3.1 Relationship between the CEL and the BI-based characterization parameters

According to the formula of the SEL, SSL, MEL and MSSL and the characterization parameters of them in Table 1, it can be concluded that the last three values are all SEL-based refining intensities, and the characterization parameters for these are all CEL-based. Both the MEL and SEL represent the energy consumption for the unit bar length, resulting in their characterizations and dimensions being similar and their values being proportional for the same patterns of bars and grooves. However, the SSL represents the energy consumed per unit area considering both the bar width and bar edge, meaning that its dimension is different than the SEL and MEL. The BIA and BIL based on the bar interaction point N and bar interaction parallelogram BIP are different from the SEL-based intensities due to the different analysis methods. However, the BI-based intensities are also machine intensities.

The values for N, BIL and BIA are calculated using the data in references 22 and 23, which showed that their values were slightly higher than the results obtained by Elahimehr. From an analysis of the data shown in Fig. 5, we see that the relationship between N2 and N1 calculated by Elahimehr can be expressed by: N2=1.42237N1+0.01888.

https://cdn.apub.kr/journalsite/sites/ktappi/2019-051-02/N0460510203/images/JKTAPPI_2019_v51n2_26_f005.jpg
Fig. 5.

The relationship among the different CEL, N, BIL and BIA values for the same conditions.

Meanwhile, the bar interaction points N and CEL are positively correlated, indicating that the CEL can better reflect the bar interaction point N. Additionally, the value of BIL increases with CEL, and BIA decreases as CEL increases; however, no obvious proportional relationship is found between these parameters and CEL, meaning that these values are complex and influenced by various parameters, such as the bar width, groove width and bar interaction angle. Therefore, we can assume that the characterization performance of the BI-based refining intensities will be different during the refining process.

3.2 The characterization performance of different refining intensities

Refining is a process that reforms the fiber properties and morphology to achieve specific fibrillation and cutting. The beating degree and average fiber length are employed to analyze changes in the fiber properties, and predictions and characterization performance of the beating degree and the average fiber length are studied considering different intensities, such as the SRE, SEL, SEL-based intensities and BI-based intensities.

Fig. 6 shows that the SRE can predict the beating degree of the pulp well without considering the patterns of bars and grooves used as part of refining; however, a deviation existed when predicting it using SEL-based and BI-based intensities. The trend in the changes in the beating degree predicted by the SEL-based intensities is consistent with the actual measured value but there is a large deviation in the prediction of accurate values. However, the SSL has a better forecasting ability when predicting the beating degree compared with other SEL-based intensities because it considers bar parameters more comprehensive to produce a higher relative accuracy. Considering the bar interactions, the refining intensity based on BIA is consistent with the predictions of the changing beating degree when P/(nBIA)<0.05 kW/m2, and when the measurement and prediction performance of the beating degree for P/(nBIL) and SSL are determined to be similar through a comparison of Fig. 6(b) and Fig. 6(e).

https://cdn.apub.kr/journalsite/sites/ktappi/2019-051-02/N0460510203/images/JKTAPPI_2019_v51n2_26_f006.jpg
Fig. 6.

The relationship between the refining intensities based on the different characterization parameters and changes in the beating degree of the pulp under constant conditions.

In a comprehensive analysis, the SRE can better measure and predict the change in the beating degree and Pt/(nBIA) can better predict the degree of beating within a certain range. Compared with other SEL-based intensities, the prediction and measurement of P/(nBIL) and SSL is slightly better.

A reasonable characterization of the average fiber length during the refining process is essential for controlling the pulp and its properties. By analyzing the relationship between the refining intensities based on the different characterization parameters and the average fiber length, it can be concluded that the SSL, which was established using a wide variety of parameters and was compared with other SEL-based intensities and the intensity based on the BIL(P/(nBIL)), can better indicate changes in the average fiber length during refining, as shown in Fig. 7. However, P/(nBIA) has a better prediction performance when its value is less than 0.045 kW/m2. The characterization performance of the SRE on fiber average length is similar to that of P/(nBIA), which can predict the average fiber length within a certain range, and a larger deviation will occur if the SRE is beyond a critical value.

https://cdn.apub.kr/journalsite/sites/ktappi/2019-051-02/N0460510203/images/JKTAPPI_2019_v51n2_26_f007.jpg
Fig. 7.

The relationship between the refining intensities based on the different characterization parameters and the change in the average fiber length for the same conditions.

It is not difficult to comprehensively conclude that the change in P/(nBIL) and the SSL correlate well with the average fiber length. Meanwhile, P/(nBIA) and the SRE can better characterize the change in the average fiber length when its value is less than the critical value.

The above analysis includes the relationship between various intensities based on the different characterization parameters and the properties of the pulp and fiber and the relationship between the CEL and the characterization parameters based on the bar interactions. However, there are a variety of parameters describing the fiber and pulp properties that need to be characterized in addition to the beating degree and average fiber length. The intensities mentioned above cannot be used to characterize all the properties of the fiber and pulp, meaning further studies are needed to better understand the refining intensity.

4. Conclusions

In this paper, the characterization parameters associated with refining intensities were studied by considering the CEL and bar interactions. Based on previous studies, the BSA was used as the effective refining area instead of the entire disc area. Two characterization parameters based on bar interactions, BIL and BIA, were calculated by defining a bar interaction parallelogram. The intensities based on BIA and BIL were established and computed using experimental data from previous studies. The following conclusions can be drawn from the data analysis:

  • (1) The SEL-based refining intensities and BI-based intensities, both of which consider the patterns of bars and grooves and their individual natures are slightly different. The value of the CEL can directly reflect the number of bar interaction points N, but do not describe the changes in BIA and BIL, where the CEL is only related to the bar length and bar angle. Additionally, the two characterization parameters based on bar interactions are complicated by the influence of many other parameters, such as bar width, groove width and bar angle.

  • (2) The SRE is the best indicator when measuring changes in beating degree or freeness during refining, and the SSL can predict these properties slightly better than other SEL-based intensities, such as the SEL and MEL. Meanwhile, the characterization performance of P/(nBIA) is better than that of P/(nBIL), but it can only partially reflect the changes in the beating degree.

  • (3) The SSL and P/(nBIL) are superior for predicting the average fiber length and can accurately measure the changes in the fiber length before and after refining. However, the SRE and P/(nBIA) can only predict the fiber length within a certain range, which is better than the SEL and MEL, which cannot measure fiber length very well.

  • (4) The SSL, SRE, P/(nBIA), and P/(nBIL) have their advantages and disadvantages and are not precise and comprehensive refining intensities. Meanwhile, CEL(Bs+Br)/[2cos(r/2)], BIA and BIL are also not optimal characterization parameters for describing the refining process.

Acknowledgements

This work was supported by the National Natural Science Foundation of China (Grant No. 50745048).

Appendix A

In this project, a plate with a central angle of β is assumed and the main parameters are for the analysis shown in Fig. 8.

https://cdn.apub.kr/journalsite/sites/ktappi/2019-051-02/N0460510203/images/JKTAPPI_2019_v51n2_26_f008.jpg
Fig. 8.

Micro-element dr of the bar surface area.

The calculation of the BSA is similar to the previous method for calculating the CEL, and the segment with the specific patterns of bars and grooves can be radially divided into numerous tiny annular areas. Assuming that the number of bars for a micro-element is n(r), the area of a single bar interaction micro element quadrilateral dS can be estimated by Eq. 14:

[14]
dS=drcosαB

Usually, α is 10-20°, and the formula for calculating the area of the bar surface area dBSAs of the ring-shaped micro-element for the specific segment is:

[15]
dBSAS=nrBcosαdr

There is an angle Ψ between the bar and the edge of the segment due to its design and manufacture. Therefore, in the actual bar interaction process, only part of the bar surface area is effective, and it is necessary to eliminate the invalid area when calculating the effective BIA. The area of the ineffective refining area of the ring-shaped micro-element is dAs, as shown in Fig. 9.

https://cdn.apub.kr/journalsite/sites/ktappi/2019-051-02/N0460510203/images/JKTAPPI_2019_v51n2_26_f009.jpg
Fig. 9.

Invalid bar interaction area.

The variable dAs is a function of the radius of the elementary annulus and the angle Ψ, which can be determined by solving the system of equations in Eq. 16 and Eq. 17:

[16]
dAs=rϕdr
[17]
ϕ=arcsinR0sinϕr-ψ

The invalid bar surface area dAsb can be approximated by calculating the dAs of the ineffective refining area considering the bar width B and the groove width G:

[18]
dAsb=BB+GdAs

Then, the expression of the effective refining area micro-element dAe can be replaced by Eq. 19:

[19]
dAe=dBSAs-dAsb

To calculate the number of bar interaction points N, it is necessary to assume that the bar interaction parallelogram BIP is formed by the interaction of a single bar on the rotor (or stator), and the adjacent groove on the stator (or rotor) is composed of two interaction points: a bar interaction point and a bar and groove interaction point. The BIP is shown in Fig. 10.

https://cdn.apub.kr/journalsite/sites/ktappi/2019-051-02/N0460510203/images/JKTAPPI_2019_v51n2_26_f010.jpg
Fig. 10.

BIP and its area calculation.

The BIP area is the sum of the single bar interaction area Ab and the area of the bar and groove

interaction point Ag. As depicted in Fig. 10, the interaction angle between the rotor bar and stator bar is γ, and the ABIP is obtained based on the calculation of the area of the parallelogram:

[20]
ABIP=BB+Gsinγ

And the range of γ is 0-π.

The number of bar interaction points in a single segment is dN, which can be solved by using Eq. 4, Eq. 19, and Eq. 20:

[21]
dN=dAeBB+Gsinγ

In fact, there is a link between the variation in the bar interaction angle γ and the existing plate patterns of the bars and grooves during the bar interaction process. Fig. 11 shows that γ is divided into the bar interaction angle of the rotor γr and bar interaction angle of the stator γs according to the actual bar interactions. When the stator and rotor plates are the same, the relationship between γ, Ψ and β is easy to from Fig. 8, which can be expressed as Eq. 22 and Eq. 23:

https://cdn.apub.kr/journalsite/sites/ktappi/2019-051-02/N0460510203/images/JKTAPPI_2019_v51n2_26_f011.jpg
Fig. 11.

Bar interaction angle.

[22]
ψ<rs<ψ+β
[23]
ψ<rr<ψ+β

Therefore, the average value for total number of bar interaction points on the elementary annulus dN can be obtained by the integration of γ from 2Ψ to 2Ψ+2β:

[24]
dN=dAeBB+G2ψ2ψ+2βsinγdr2β

To simplify the calculation, the arc length composed of a bar and a groove can be estimated by the sum of the bar width B and the groove G. Then, the number of bars in a segment can be expressed as Eq. 25:

[25]
nr=βrB+G

Then, by substituting Eq. 19 and Eq. 25 for dAe in Eq. 24, the bar interaction points N for the segment with a central angle β may be calculated by integration between Ri and Ro:

[26]
N=Hψ,βX1β,ψ,Ri,R02B+G2

where and are defined in Eq. 27 and Eq. 28:

[27]
Hψ,β=cos2ψ-cos2ψ+2β2β
[28]
X1β,ψ,Ri,R0=4βR02-Ri2cosα-R02ψ+Ri2arcsinR0sinψRi+R0cosψ-Ri2-R02sin2ψ

Taking into account the number of segments, the total bar interaction points for refining can be calculated by Eq. 29:

[29]
NT=NSN

where the NS is the number of segments.

The sum of the adjacent two sides of the bar interaction point in the BIP, which is shown in Fig. 6, is defined as the single bar interaction length CBIL. The single bar interaction area Ab is defined above. These values can be obtained by Eq. 30 and Eq. 31.

[30]
CBIL=2Bsinγ
[31]
AB=B2sinγ

The calculation of CBIL and Ab both involve the bar interaction angle γ, and both of these values can provide relatively accurate values using the bar interaction points N. However, the bar width in the actual refining process is small and the effects of γ can be ignored.

The BIL and BIA of the entire refining area can be expressed by Eq. 32 and Eq. 33.

[32]
BIL=2BNSHψ,βX1β,ψ,Ri,R02B+G2
[33]
BIA=B2NSHψ,βX1β,ψ,Ri,R02B+G2

Literature Cited

1
W. Brecht, TAPPI Journal, A method for comparative evaluation of barequipped beating devices, 50(8); 40-44 (1967)

Brecht, W., A method for comparative evaluation of barequipped beating devices, TAPPI Journal 50(8):40-44 (1967).

2
F. Wultsch and W. Flucher, Papier, Der escher-wyss-kleinrefiner als standard-prüfgerät für moderne stoffaufbereitungsanlagen, 12(13); 13-18 (1958)

Wultsch, F. and Flucher, W., Der escher-wyss-kleinrefiner als standard-prüfgerät für moderne stoffaufbereitungsanlagen, Papier, 12(13): 13-18 (1958).

3
W. Brecht and W. Siewert, Papier, Zur Theoretisch-technischen beurteilung des mahlprozesses moderner mahlmaschinen, 20(1); 4-14 (1966)

Brecht, W. and Siewert, W., Zur Theoretisch-technischen beurteilung des mahlprozesses moderner mahlmaschinen, Papier 20(1):4-14 (1966).

4
J. Lumiainen, Pulp and Paper International, New theory can improve practice, 32(7); 46-47 (1990)

Lumiainen, J., New theory can improve practice, Pulp and Paper International 32(7):46-47 (1990).

5
J. Lumiainen, Specific surface load theoryProceeding of 3rd PIRA International Refining Conference, Atlanta; 56-60 (1995)

Lumiainen, J., Specific surface load theory, Proceeding of 3rd PIRA International Refining Conference, Atlanta, pp. 56-60 (1995).

6
F. P. Meltzer, Technologie der zellstoffmahlung, Aachen, Germany. RWTH Aachen University. (1994)

Meltzer, F. P., Technologie der zellstoffmahlung, Ph.D. Thesis, RWTH Aachen University, Aachen, Germany (1994).

7
F. P. Meltzer and P. W. Sepke, New ways to forecast the technological results of refiningProceeding of 3rd PIRA International Refining Conference, Atlanta; 20-22 (1995)

Meltzer, F. P. and Sepke, P. W., New ways to forecast the technological results of refining, Proceeding of 3rd PIRA International Refining Conference, Atlanta, pp. 20-22 (1995).

8
R. Musselman, D. Letarte and R. Simard, Third stage low consistency refining of TMP for energy savings and quality enhancementProceeding of 4th PIRA International Refining Conference, Fiuggi, Italy; 139-148 (1997)

Musselman, R., Letarte, D. and Simard, R., Third stage low consistency refining of TMP for energy savings and quality enhancement, Proceeding of 4th PIRA International Refining Conference, Fiuggi, Italy, pp. 139-148 (1997).

9
J. C. Roux and G. Joris, Angular parameters beyond specific edge load, TAPPSA Journal. (2005)

Roux, J. C. and Joris, G., Angular parameters beyond specific edge load, TAPPSA Journal (2005).

10
H. Liu, J. Dong, X. Guo, L. Qiao and H. Jing, China Pulp & Paper, Quantitative analysis of pulp refining and its research progress, 37(8); 66-71 (2018)

Liu, H., Dong, J., Guo, X., Qiao, L., and Jing, H., Quantitative analysis of pulp refining and its research progress, China Pulp & Paper 37(8):66-71 (2018).

11
J. Lewis and D. W. Danforth, TAPPI Journal, Stock preparation analysis, 45(3); 185-188 (1962)

Lewis, J. and Danforth, D. W., Stock preparation analysis, TAPPI Journal 45(3):185-188 (1962).

12
P. J. Leider and A. H. Nissan, TAPPI Journal, Understanding a disk refiner, 60(10); 85-89 (1977)

Leider, P. J. and Nissan, A. H., Understanding a disk refiner, TAPPI Journal 60(10):85-89 (1977).

13
P. J. Leider and A. H. Nissan, TAPPI Journal, The average number of impacts experienced by a fiber in a disk refiner: Response, 61(8); 91-92 (1978)

Leider, P. J. and Nissan, A. H., The average number of impacts experienced by a fiber in a disk refiner: Response, TAPPI Journal 61(8):91-92 (1978).

14
R. J. Kerekes, Nordic Pulp Paper Research Journal, Characterization of pulp refiners by a C-factor, 16(1); 3-8 (1990)

Kerekes, R. J., Characterization of pulp refiners by a C-factor, Nordic Pulp Paper Research Journal 16(1):3-8 (1990).

10.3183/npprj-1990-05-01-p003-008
15
M. Pfaffli, Rethinking the Art of Refining: Improving the Efficiency and Quality of Refining, Beijing, China. CIPTE. (2015)

Pfaffli, M., Rethinking the Art of Refining: Improving the Efficiency and Quality of Refining, CIPTE, Beijing, China (2015).

16
T. Berger, Papermaking success stories-Optimized refiner plate designsPaper Conference, Covington, USA (2011)

Berger, T., Papermaking success stories-Optimized refiner plate designs, Paper Conference, Covington, USA (2011).

17
H. Liu and J. Dong, China Pulp & Paper, Research progress of pulp flow in disc refiner, 36(9); 66-73 (2017)

Liu, H. and Dong, J., Research progress of pulp flow in disc refiner, China Pulp & Paper 36(9):66-73 (2017).

18
R. J. Kerekes, Journal of Pulp and Paper Science, Energy and forces in refining, 36(1-2); 10-15 (2010)

Kerekes, R. J., Energy and forces in refining, Journal of Pulp and Paper Science 36(1-2):10-15 (2010).

19
M. A. Meng and Z. Jin-qiang, China Pulp and Paper Industry, JC-04 conical refiner Excellent beating equipment in modern paper board mill, 23(10); 22-25 (2002)

Meng, M. A. and Jin-qiang, Z., JC-04 conical refiner Excellent beating equipment in modern paper board mill, China Pulp and Paper Industry 23(10):22-25 (2002).

20
R. J. Kerekes, Paper and Biomaterials, Quantifying the pulp refining process, 1(1); 38-41 (2016)

Kerekes, R. J., Quantifying the pulp refining process, Paper and Biomaterials 1(1):38-41 (2016).

21
H. Paulapuro and P. Hannu (ed.), 2nd Edition. Stock Preparation and Wet End; 42, Helsinki, Finland. Fapet Oy. (2000)

Paulapuro, H., Stock Preparation and Wet End, Hannu, P. (ed.), 2nd Edition, Fapet Oy, Helsinki, Finland, p. 42 (2000).

22
A. Elahimehr, J. A. Olson and D. M. Martinez, Nordic Pulp and Paper Research Journal, Estimating the area and number of bar crossings between refiner plates, 27(5); 836-843 (2012)

Elahimehr, A., Olson, J. A., and Martinez, D. M., Estimating the area and number of bar crossings between refiner plates, Nordic Pulp and Paper Research Journal 27(5):836-843 (2012).

23
A. Elahimehr, Low consistency refining of mechanical pulp: The relationship between plate pattern, operational variables and pulp properties; 95-101, Vancouver, Canada. The University of British Columbia. (2014)

Elahimehr A., Low consistency refining of mechanical pulp: The relationship between plate pattern, operational variables and pulp properties, Ph.D. Thesis, The University of British Columbia, Vancouver, Canada, pp. 95-101 (2014).

페이지 상단으로 이동하기