Temperature Dependence of Nylon and PTFE Triboelectrification
Isaac A. Harris, Melody X. Lim, Heinrich M. Jaeger

TL;DR
This study investigates how temperature affects tribocharging between nylon and PTFE, revealing that higher temperatures reduce charge transfer, likely due to changes in surface water adsorption, with implications for industrial applications.
Contribution
It demonstrates the temperature dependence of triboelectrification in nylon and PTFE and highlights the role of surface water in charge transfer mechanisms.
Findings
Increasing temperature decreases charge transfer magnitude.
Surface water adsorption is crucial for charge exchange.
Results are supported by Monte-Carlo simulations and TGA measurements.
Abstract
Contact electrification, or tribocharging, is pertinent to a broad range of industrial and natural processes involving dielectric materials. However, the basic mechanism by which charge is transferred between insulators is still unclear. Here, we use a simple apparatus that brings two macroscopic surfaces into repeated contact and measures the charge on the surfaces after each contact. We vary the temperature of the surfaces, and find that increasing temperature leads to a decrease in the magnitude of charge transfer. When paired with a Monte-Carlo simulation and TGA measurements, our results support a mechanism where adsorbed surface water is crucial for charge exchange. Our setup is easily extendable to a variety of industrially relevant materials.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Temperature Dependence of Nylon and PTFE Triboelectrification
Isaac A. Harris, Melody X. Lim, Heinrich M. Jaeger
Department of Physics and James Franck Institute, The University of Chicago, 5720 S. Ellis Ave, Chicago, Illinois 60637, USA
Abstract
Contact electrification, or tribocharging, is pertinent to a broad range of industrial and natural processes involving dielectric materials. However, the basic mechanism by which charge is transferred between insulators is still unclear. Here, we use a simple apparatus that brings two macroscopic surfaces into repeated contact and measures the charge on the surfaces after each contact. We vary the temperature of the surfaces, and find that increasing temperature leads to a decrease in the magnitude of charge transfer. When paired with a Monte-Carlo simulation and TGA measurements, our results support a mechanism where adsorbed surface water is crucial for charge exchange. Our setup is easily extendable to a variety of industrially relevant materials.
I Introduction
Contact charging forms the basis for many industrial applications, including powder coating bailey1998science , cleaning surfaces kawamoto2015electrostatic , and air filtration jaworek2007modern . At the same time, tribocharging also presents a challenge to powder handling in industrial processes, as the separation of charge may cause the formation of clusters ciborowski1962electrostatic ; boland1969static ; liang1996electrostatic ; cocco2010particle , or create sparks that damage electronic components voldman1999state . In addition to its widespread industrial relevance, the buildup of charge from frictional collisions is thought to be responsible for natural phenomena ranging from lightning inside volcanic ash clouds brook1974lightning ; pahtz2010particle to potentially also the very early stages of the formation of planetesimals from interstellar dust blum2008growth .
Despite the importance of phenomena related to triboelectric charging, the underlying mechanism for charge exchange is not well understood. In particular, it is unclear how insulators, which have very low charge mobility, can transfer charge. Several mechanisms for tribocharging have been proposed, including the transfer of electrons in trapped states lowell1986triboelectrification ; liu2010electrostatic ; lacks2008nonequilibrium , ion transfer mccarty2008electrostatic ; lee2018collisional , and the transfer of nanoscale pieces of charged material baytekin2011mosaic .
Tribocharging has been shown to be sensitive to a wide variety of conditions, including the ambient environment conditions greason1999investigation ; matsuyama1995charge ; schella2017influence , impact velocity xie2016effect , material strain sow2012strain , the contact force liu2013contact , and the microscopic details of the surface contact williams2012triboelectric ; horn1992contact . Previous work has reduced these effects by making highly controlled, precision measurements of the impact charging of a small particle matsuyama2003impact ; watanabe2006measurement ; lee2018collisional ; haeberle2018double , or by using single-crystal materials that are flat to nanometer precision stern1988deposition ; collins2018simultaneous . However, even with well-controlled surfaces, the inherent stochasticity of tribocharging can lead to a large variance in the results.
We present an experiment that takes a different approach by measuring the average charge transfer between macroscopic materials. In doing so, we effectively average over many microscopic charge transfer processes. Despite the fact that we do not control all microscopic aspects of the surface contact, we show that we are nevertheless able to extract robust statistical trends from our data. Here, we measure the charge transfer between materials on opposite ends of the triboelectric series mccarty2008electrostatic : nylon and polytetrafluouroethylene (PTFE). As we increase the surface temperature of the macroscopic samples, we observe a decrease in the overall magnitude of charge transfer. Our data, in conjunction with simulations and thermogravimetric analysis (TGA) measurements, suggest that adsorbed surface water plays a key role in the charge transfer between dielectric surfaces.
II Experimental methods and results
Our setup is illustrated in Fig. 1a. A macroscopic piece of PTFE ( 30303.2 mm3) is housed in a Faraday cup, which is shielded by a Faraday cage. At the same time, a matching piece of nylon ( 25253.2 mm3) is attached to the end of a pneumatic actuator, which lowers the nylon into the Faraday cup, making contact with the PTFE. Both the nylon and PTFE are cleaned before an experiment (see Appendix A for details). Our experiments take place in a humidity-controlled (42-44% RH) environment. An electrometer (Keithley 6514) measures the total charge on nylon and PTFE during contact (10 seconds), then the charge on the PTFE after contact (10 seconds). These two measurements are then subtracted, such that the final measurement is the charge on the nylon after each collision. The apparatus thus allows for the automated control of a sequence of collisions, and the measurement of the charge as a function of collision number.
In addition to automating repeated contact between a pair of dielectrics, our apparatus also allows for the independent control of the surface temperature of both the nylon and the PTFE during an experiment. A PID-controlled heating element and thermistor were placed on the back side of each dielectric (see Appendix B for calibration curves). During an experiment, the nylon and PTFE are both held at a specified temperature for 15 minutes, then undergo a sequence of 30 collisions over a time span of about 10 minutes. We vary the surface temperature of the dielectrics from room temperature (25∘C) to 80∘C.
Representative examples of charging curves are shown in Fig. 1b. When the PTFE is removed from the Faraday cup (so that the nylon does not touch anything when the linear actuator moves downward) the charge on the nylon does not systematically change at any temperature (grey data), suggesting that over the timescales of the experiment (30 collisions), the dissipation of charge from the nylon is negligible. In contrast, when the nylon and PTFE are brought into repeated contact, the nylon acquires a positive charge that increases with collision number (green and orange data). This is consistent with previous measurements of the triboelectric series clint2001acid ; diaz2004semi ; mccarty2008electrostatic ; lacks2011contact ; galembeck2014friction . Qualitatively, the asymptotic accumulation of charge on the nylon at 78∘C is smaller than at 40∘C. These results are consistent with the trends reported in Ref. greason1999investigation .
These trends are reproducible over many experiments, despite the relatively uncontrolled nature of the contact between the nylon and PTFE. Fig. 2 shows examples of the charge time-series obtained for several different experiments, each with freshly prepared samples, at two different temperatures (34 ∘C in part a, and 77 ∘C for part b). There is variation between each time-series, due to microscopic differences between the surfaces. At the same time, the overall shape of the curves is highly consistent for each temperature. In light of the consistency in the data, we focus on average measures for each temperature, thus providing insight to the underlying statistical processes during contact charging.
Our data, as shown in Figs. 1b and 2 exhibit a rapid increase of charge as a function of collision number, a behavior which has been reported in other contact electrification experiments as well collins2018simultaneous ; liu2013contact ; greason1999investigation ; apodaca2010contact ; matsusaka2010triboelectric . In addition to the initial increase, we typically observe a much more slowly increasing background, which similar experiments by others did not report. We make these observations quantitative by fitting our experimental data to a variation of the condenser model proposed in Refs. matsusaka2000electrification ; matsusaka2007control ; matsusaka2010triboelectric . This phenomenological model provides a mechanism for an exponential buildup of charge towards some saturation value . Although the model also includes a second exponential term to represent charge dissipation, our results (Fig. 1b) suggest that this term is not relevant within the timescales of our experiment. In addition, we include a linear term to account for the slower background that we observe (as we will discuss below, we link this background to slow drift in the contacting portions of the surfaces being pushed together by the actuator). We thus fit the charge as a function of collision number to the following function:
[TABLE]
where , , and are fitting parameters.
Figure 3 shows the variation of with surface temperature. Our surfaces are macroscopic and thus may vary widely in their microscopic properties, leading to previously observed stochastic charging processes medley1950frictional ; horn1992contact ; matsusaka2010triboelectric ; apodaca2010contact ; haeberle2018double . Despite this statistical nature of contact charging, we nevertheless observe an average decrease of with increasing temperature. Given that corresponds to the asymptotic saturation of charge on the surface in the absence of any other effects, Fig. 3 suggests that increasing the surface temperature reduces the capacity of the surfaces to exchange charge. This result suggests that contact charging in our experiment is mediated by a mechanism that is responsive to this small temperature change (room temperature to 80∘C).
In contrast, neither nor show systematic variation with surface temperature (Fig. 4a and b respectively): fluctuates around 2 collisions, and fluctuates around 2 pC/collision. Returning to the phenomenological model in Eq. (1), we interpret as a characteristic timescale for charge exchange, and as a slow background increase in the total charge. The fact that does not vary with temperature suggests that this linear term originates in a systematic feature of the experiment. Similarly, the fact that does not show systematic variation with temperature suggests that although the temperature affects the total capacity to exchange charge (Fig. 3), it does not affect the timescale (and, by implication, the mechanism) by which this charge is exchanged.
III Model for temperature dependence of charge exchange
In order to explain our experimental results, we develop a simple computational model for contact charging. We base our model on the charging mechanisms proposed in Refs. wiles2004effects ; thomas2008patterns ; mccarty2008electrostatic ; lee17 ; lee2018collisional , where the transfer of charge is associated with the collision of nanoscale dry and wet patches on the surfaces. Dry (wet) here refers to the absence (presence) of molecularly thin surface layers of water, which are ubiquitous under ambient conditions. During such a collision, adsorbed OH- ions can be transferred from a wet area on one surface to a dry area on another surface yatsuzuka1996electrification ; zimmermann2001electrokinetic ; kudin2008water . Depending on the number of dry and wet patches, which depends on the hydrophobicity of the surface ewing2006ambient ; kudin2008water ; xu2010graphene , there may be a net charge transfer as one surface transfers its wet patches to another surface yatsuzuka1994electrification ; yatsuzuka1996electrification ; ravelo2011demonstration ; burgo2016water .
We conduct a Monte Carlo simulation of the two insulator surfaces in order to model this charge transfer process. These surfaces are partitioned into wet and dry patches. As an initial simplification, we treat the hydrophobic PTFE surface as if it did not have any wet patches (see Appendix D for supporting data). By contrast, the nylon surface is initially covered with a random distribution of wet patches (see Appendix E for details). Here, we use the parameter (which increases the proportion of wet to dry patches) as the control variable for the amount of surface water. We assume that the amount of negative charge (OH- ions) available for transfer is directly proportional to the amount of water on the insulator surface. An example of the simulated initial distribution of surface water (and thus negative surface charge) on the nylon and PTFE is shown in Fig. 5.
Given the initial distribution of charges on the nylon and PTFE pieces, we then implement collisions by transferring a fraction of the charge on a wet patch to a corresponding dry patch on the PTFE. In addition, we stochastically translate the area of contact between the nylon and PTFE surfaces for each collision. This mimics our experimental conditions, where the areas of contact can vary from collision to collision because of inadvertent, small amounts of drift in the linear actuator. To implement this in the simulations, we used the following scheme. We defined a square with side length , called the contact area, and using a random walk, we displaced the contact area a distance in a random direction each collision. During each collision, only the elements within the contact area were allowed to exchange charge. In the simulation, this produces a slowly increasing background that can be approximated by a linear term in the charge accumulation, similar to what we see in our experiments. Examples of the evolution of the surface charge distribution over several collisions are shown in Fig. 5.
In order to calibrate the simulations to the effect of temperature on the nylon and PTFE surfaces, we measured the mass of water adsorbed by a piece of nylon at different temperatures using thermogravimetric analysis (TA Instruments Q600 SDT). We assume that the mass of surface adsorbed water decreases proportionally to the bulk water content, since the mass of surface adsorbed water is too small to measure directly. We heat a piece of nylon ( mm3, mass = 6 mg) at a rate of 20 ∘C/minute, until it reaches a specified temperature between 30 and 80 ∘C. The temperature is then held constant, and the mass of the sample is measured as a function of time. Data was taken at atmospheric pressure, under inert nitrogen (flow rate 100mL/minute). Prior work has shown that the onset decomposition temperature of nylon is over 400∘C carrizales2008thermal , suggesting that mass loss over the range of temperatures that we probe can be attributed to the evaporation of adsorbed water. We confirm that the mass loss we observe is due only to the evaporation of water by comparing the mass loss curves from a piece of nylon soaked in water overnight, and a piece of nylon dried in a vacuum oven (see Appendix C).
We first measure the total mass percentage of water adsorbed by a piece of nylon stored in ambient conditions by heating the sample to 100∘C and waiting for the mass of the sample to equilibrate (see Appendix C for data). Denoting the total initial mass of the nylon sample to be , including a mass of adsorbed water, we measure the initial fractional water content of our nylon samples to be . We then heat samples to the temperatures we use in the experiment, measuring the fractional mass of water as a function of time at each temperature (inset to Fig. 6). We mirror the conditions of the charging experiment by measuring the percentage of remaining adsorbed water as a fraction of the original mass of water, after 15 minutes at the target temperature. Data are shown in the main panel of Fig. 6. We find that the mass of water remaining in the sample decreases with temperature, and that heating to just 80 ∘C for 15 minutes evaporates 40% of the water content of a nylon sample, while there is no measurable effect for the PTFE (see Fig. 10 and Appendix D for data on the water content of a PTFE sample).
We use these results to proportionally decrease the proportion of wet patches on the simulated nylon by setting , (see Appendix E). This reproduces the effect of temperature on the transfer of charge. Fig. 7a shows simulated data for the total (dimensionless) charge on the nylon as a function of collision number , for two different temperatures (green, orange data points). The initial distribution of wet patches on the surface of the nylon is shown for each temperature (inset). Qualitatively, we reproduce the shape of the experimental results shown in Fig. 1b, and so we fit the simulation data also to Eq. (1) (black lines). The variation of the fit parameters with temperature is shown in Figs. 7b-d. In particular, we find that decreases linearly with temperature, in proportion to the experimental results in Fig. 3 (a decrease by a factor of two over the same range of temperatures), suggesting that the simulations have captured the effect of temperature in the experiments. In addition, the simulated value of fluctuates just below 2 collisions, again in good agreement with the experiments.
We note that in the simulation there is a decreasing trend in with temperature, reflecting the decreasing proportions of water on the surfaces. We do not see this trend in the experimental data in Fig. 4b; however, it is possible that larger temperatures increase the drift, which could compensate at least for part of the loss of surface water.
Our model suggests a means to reconcile previously contradictory results. While some experiments have observed an exponential increase of charge during contact charging experiments collins2018simultaneous ; liu2013contact ; greason1999investigation ; apodaca2010contact ; matsusaka2010triboelectric , other measurements on single grains have observed a linear, rather than exponential, charging when the collision area is randomised lee2018collisional . In our model, most charge is exchanged after a few collisions, with a slow, stochastic background motion of the contact area. After 2-3 collisions, the two surfaces have exchanged most of their wet patches within the contact area, leading to an exponential saturation. This local saturation is offset on further collisions by the gradual motion of the contact area, making new dry and wet patches accessible. Since the new wet patches are revealed slowly, the result is that the charge increases much more slowly than on initial contact. We propose that this separation of timescales between collisions and the motion of the contact area leads to the observed exponential charging with linear term that we observe. If the contact area motion is much faster, i.e., a new set of surface patches comes into contact every collision, then the surfaces will not locally equilibrate and thus only the linear term will be observed. In contrast, we recover purely exponential charging if the timescale for contact area motion becomes very large. Within this framework, then, the interplay of local saturation and stochastic contact motion is able to explain the variety of observed charging trends.
IV Conclusions
We have presented an experiment and simulation to study collisional charging between a pair of macroscopic dielectrics. Although the overall nature of contact charging is stochastic, and depends on the microscopic details of the two surfaces that come into contact, our results nevertheless show that statistically significant trends can be extracted from macroscopic experiments. In particular, we demonstrate that the magnitude of charge transfer between a piece of nylon and a piece of PTFE decreases as their surface temperature is increased, while the charge transfer trends remain constant. The agreement between our experiments and simulations, combined with TGA measurements, suggest a contact charging mechanism for nonionic insulators that relies on the presence of surface water. Charge is carried by OH- ions that are transferred at the sites where a wet patch on one of the surfaces collides with a dry patch on the other surface. Evaporating the surface water via a modest increase in temperature thus leads to a decrease in the overall magnitude of charge transfer, while not affecting the initial rate of exponential charge saturation.
Our experimental method can be used straightforwardly for measurements on other materials, in a wide range of different environments. We anticipate that similar experiments will continue to inform an understanding of the basic mechanism of collisional charge transfer. In addition, we envision that such macroscopic, statistical measurements will shed light on processes in tribocharging that are especially relevant to industrial applications.
Acknowledgments
We thank Grayson Jackson for his assistance with the TGA measurements, and Victor Lee for his work on the early stages of the experiment. We thank Kieran Murphy, Leah Roth, Abhinendra Singh, and Nicole James for insightful discussions. This research was supported by the National Science Foundation through grant DMR-1810390. I.A.H. acknowledges primary support from the Chicago MRSEC, funded by the NSF through grant DMR-1420709.
Appendix A: Sample preparation
Samples of Nylon 6,6 (McMaster Carr) were cut to size, then rinsed with toluene, ethanol, and then DI water. Samples of PTFE (McMaster Carr) were cut to size, then rinsed with toluene, isopropanol, ethanol, and DI water. For both samples, excess water was dried with high-purity N2 gas and then placed in a vacuum chamber with an oil-free pump for several hours. The plastics were then transferred into the apparatus and were left to discharge under ambient conditions for 5 or more hours. Using heating pads attached to the mounts for the nylon and PTFE, we then heated the plastics to the desired temperature with PID controllers, attached to thermistors on the underside of the plastics.
Appendix B: Temperature calibration
In order to calibrate our temperature sensors, we ran tests outside of the apparatus with the PID controllers set to heat to a specified temperature, while we also monitored the surface temperature of the plastics with a separate thermistor placed on the sample surface. Our results for PTFE are shown in Fig. 8, on the basis of which we wait 15 minutes for the temperature of the samples to equilibrate (and thus take readings after 15 minutes at constant temperature in the TGA measurements). The results are identical for the nylon samples. We use the actual surface temperature for the data shown in the main text.
Appendix C: TGA measurements for dry and wet nylon
In order to verify that the mass loss we measure is due to the evaporation of bulk water, we performed control measurements with three nylon samples ( mm3, mass = 6 mg). Our results are plotted in Fig. 9 . One of the samples was soaked in DI water for 12 hours (yellow curve). The other nylon sample was dried in a vacuum oven for 2 hours at 100-140∘C (uncertainty due to the vacuum oven), then held at 10 hours under vacuum (15 psi) in the same oven (purple curve). These were also compared to a sample stored at ambient conditions (same as the other samples we test in the main text, magenta curve). Samples were heated at a rate of 10 ∘C/minute, until a temperature of 100∘C was reached. The temperature was then held constant, and the mass of the sample measured as a function of time. Data was taken at atmospheric pressure, under inert nitrogen (flow rate 100mL/minute).
We find that the sample that was dried shows very little mass loss compared to the sample stored at ambient conditions, and that the mass loss is greater (by a factor of 2) from the soaked nylon. We thus conclude that the temperature-dependent mass loss we measure can be attributed to the loss of bulk water, within the accuracy of our experiment.
Appendix D: TGA measurements for PTFE
In order to justify our approximation in the simulations that the PTFE begins with no wet patches, we measured the mass loss of a PTFE sample ( mm, mass = 18 mg) using TGA. Samples (stored in ambient conditions) were heated at a rate of 2 ∘C/minute, until a temperature of 100∘C was reached. The temperature was then held constant, and the mass of the sample measured as a function of time. Data were taken at atmospheric pressure, under inert nitrogen (flow rate 100mL/minute). Data are shown in Fig. 10. The mass loss from a nylon sample stored in the same conditions is also plotted for comparison. Over the course of 5 hours, the mass of the PTFE equilibrated to 99.93% of its original mass (in contrast to 97% for the nylon sample). Our results suggest that the mass of water adsorbed by the PTFE is 50 times smaller than the mass of water adsorbed by the nylon.
Appendix E: Simulation protocol
We stochastically generate the wet patches on the nylon. We initialise a 100100 matrix , representing the nylon surface, whose entries are either 1 or 0, with probability of generating a 1 in any given entry. The nonzero entries of the matrices serve as wet patches. Typical values of in our simulations were between 0.3 and 0.6. We choose the proportion of wet patches on room temperature nylon as . Since we are only interested in the proportionality in across different temperatures as measured by the TGA, the value of is not important. As described in the main text, we simulate the effect of increasing temperature (and thus decreasing surface water) by decreasing in proportion to our TGA measurements (Fig. 6): .
We stochastically translate the area of contact between two surfaces in order to simulate the drift of the actuator, or the general change in charging surfaces over repeated contacts. For any given collision, we model the contact area as a square of side length . Surfaces outside of this contact area do not exchange charge. Over the course of the simulated experiment, this contact area is displaced over the nylon matrix as a random walk with a displacement each collision, with the direction of the displacement also chosen randomly from the four cardinal directions (for instance means that the nylon contact area shifts by 6 matrix elements in the or direction for each collision). For the results shown in the main text, and . For a discussion of how the results change for different parameter choices, see Appendix F.
After the area of contact has been stochastically translated in this way, we implement our charge transfer rules: we compare with only within the designated contact area. If the local charge on the PTFE is smaller than the corresponding local charge on the nylon, they exchange a fraction of the charge difference between them. Formally, within the contact area, the local charge on the PTFE after collision is given by
[TABLE]
To conserve charge, the corresponding local charge on the nylon after collision is
[TABLE]
The total charge on the PTFE (as plotted in Fig. 7a) is then . By conservation of charge, the total charge on the nylon .
As shown in Fig. 11(c), controls the rate of the exponential charge saturation in our simulation: as decreases, increases. In order to match our simulation (Fig. 7) to our experiments (Fig. 4), we set to 0.27. We note that in the experiment, several different physical mechanisms could control . In the simulation, we summarise the combined effect of these mechanisms with the single variable .
Appendix F: Robustness of model results
In order to test the robustness of the results that we present in Sec. III of the main text, we reproduce our simulation model results across a wide range of parameters. In Fig. 11, we vary the displacement per collision , the proportion of wet patches on the nylon , and the fraction of charge transferred per collision . In order isolate the effect of each parameter, the (random) motion of the contact area was generated for one charge sequence, then used for all other parameter choices.
Fig. 11(a) shows normalised charge time-series for varying values of , the ratio of the width of the contact area and the displacement of the contact area per collision . For small , the displacement of the contact area and the width of the contact area are of comparable size, leading to small overlaps in the areas that exchange charge from collision to collision. This lack of self-overlap leads to a lack of saturation in the underlying surface, and thus to a linear increase in charge with collision number. As is increased, the overlap of the contact area from collision to collision increases. Correspondingly, the charge time-sequences become increasingly exponential.
Fig. 11(b) shows charge time-series for different values of , the proportion of wet patches on the nylon matrix. We note that increasing changes only the magnitude of the final charge ( ) and not the timescale of saturation or the functional form of the charge increase. In contrast, Fig. 11(c) shows that the timescale of the exponential charge saturation is controlled by , the fraction of charge transferred each collision.
Appendix G: Surface Characterization
Images of the microscopic surface structure of the samples (see Fig. 12) were obtained using an Olympus LEXT OLS5000 surface scanner.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1(1) Adrian G Bailey. The science and technology of electrostatic powder spraying, transport and coating 1. Journal of electrostatics , 45(2):85–120, 1998.
- 2(2) Hiroyuki Kawamoto and Takuya Shibata. Electrostatic cleaning system for removal of sand from solar panels. Journal of Electrostatics , 73:65–70, 2015.
- 3(3) Anatol Jaworek, Andrzej Krupa, and Tadeusz Czech. Modern electrostatic devices and methods for exhaust gas cleaning: A brief review. Journal of electrostatics , 65(3):133–155, 2007.
- 4(4) J Ciborowski and A_ Wlodarski. On electrostatic effects in fluidized beds. Chemical Engineering Science , 17(1):23–32, 1962.
- 5(5) D Boland, QAW Al-Salim, and D Geldart. Static electrification in fluidised beds. Chemical Engineering Science , 24(8):1389–1390, 1969.
- 6(6) Shu-Chien Liang, Jian-Ping Zhang, and Liang-Shih Fan. Electrostatic characteristics of hydrated lime powder during transport. Industrial & engineering chemistry research , 35(8):2748–2755, 1996.
- 7(7) Ray Cocco, Frank Shaffer, Roy Hays, SB Reddy Karri, and Ted Knowlton. Particle clusters in and above fluidized beds. Powder Technology , 203(1):3–11, 2010.
- 8(8) Steven H Voldman. The state of the art of electrostatic discharge protection: Physics, technology, circuits, design, simulation, and scaling. IEEE Journal of Solid-State Circuits , 34(9):1272–1282, 1999.
