## Journal of the Brazilian Chemical Society

*On-line version* ISSN 1678-4790

### J. Braz. Chem. Soc. vol.17 no.6 São Paulo Sept./Oct. 2006

#### http://dx.doi.org/10.1590/S0103-50532006000600007

**ARTICLE**

**Linear viscoelasticity of textured carbonaceous mesophases**

**Luiz Rogério P. de Andrade Lima ^{I,}^{*}; Alejandro D. Rey^{II} **

^{I}Departamento de Ciência e Tecnologia dos Materiais, Universidade Federal da Bahia CP 6974, 41810-971 Salvador-BA, Brazil

^{II}Department of Chemical Engineering McGill University 3610 University Street, Montreal, Canada, H3A 2B2

**ABSTRACT**

High performance carbon fibers can be produced from anisotropic petroleum pitch using the fiber spinning process. The mesophase pitch is a thermotropic discotic nematic liquid crystal possessing orientational order. Liquid crystals are textured, anisotropic, viscoleastic materials. An important phenomenon of liquid crystals is the interaction between orientational textures and rheology, where the texture is due to the presence of orientational defects. The fiber spinning process used in the carbon fiber manufacture from nematic precursors uses a cascade of processing flows that involves capillary, convergent, and extensional flows, which promotes the formation of a variety of textures that will give specific properties to the resultant carbon fiber. Hence a fundamental understanding of the impact of textures on rheology is needed for eventual process control and product optimization. In this paper we study the interaction of rheology and texture and develop a model of anisotropic linear viscoelasticity of textured mesophases. The model is able to predict all the features observed experimentally, including the frequency dependency of the loss and storage modulus and their temperature sensitivity, and provides a better understanding of the role of defects and anisotropy on pitch viscoelasticity.

**Keywords: **mesophase pitch, rheology, nematic discotic liquid crystal, carbon fiber, Leslie-Ericksen theory, texture, anisotropy

**RESUMO**

Fibras de carbono de alto desempenho podem ser produzidas a partir de piche de petróleo anisotrópico usando o processo de extrusão-fiação. O piche mesomórfico é um cristal líquido discótico nemático termotrópico, que possui ordem na orientação das moléculas. Cristais líquidos são materiais que possuem textura, anisotropia e visco-elasticidade. Um fenômeno importante que é encontrado em cristais líquidos é a interação entre a textura devido à orientação das moléculas e o comportamento reológico, no qual a textura ocorre devido à presença de defeitos na orientação das moléculas. O processo de extrusão-fiação usado na fabricação de fibras de carbono a partir de precursores nemáticos usa uma cascata de escoamentos que incluem fluxo capilar convergente e extensional, os quais promovem a formação de uma variedade de estruturas, incluindo radial, concêntrica, randômica e bipolar, que darão propriedades especificas para as fibras de carbono resultantes. Portanto, a compreensão fenomenológica da interação entre a textura e a reologia do precursor nemático é necessária para o controle do processo e para a otimização do produto. Neste artigo, a interação entre a reologia e a textura de cristais líquidos discóticos nemáticos será enfocada e um modelo anisotrópico será desenvolvido para descrever o regime linear visco-elástico de piches mesofásico com textura. O modelo é capaz de prever todas as características observadas experimentalmente, incluindo a dependência dos módulos de armazenamento e de perda com a freqüência de oscilação e a sua dependência com a temperatura, assim como promover uma melhor compreensão do papel dos defeitos e da anisotropia no comportamento visco-elástico do piche anisotrópico.

**Introduction**

High performance carbon fibers can be produced from liquid crystalline petroleum pitches, also known as carbonaceous mesophase, using the melt spinning process, which was introduced in the 1970s to produce carbon fibers with ultra-high modulus and very high thermal and electrical conductivities. The melt spinning process consists of a flow sequence that induces unique textural transformations in the mesophase.^{1-4} Improvements in the high performance carbon fiber quality require a better understanding of the interaction between the rheology and texture of carbonaceous mesophase precursors.

Carbonaceous mesophases are textured anisotropic viscoelastic materials formed by flat aromatic molecules, see Figure 1a for a representative molecule.^{1-4} This material exhibits orientational order and positional disorder and is a discotic nematic thermotropic liquid crystal.^{5,6} The description of complex fluids systems such as discotic nematic liquid crystals by molecular dynamic or other microscopic approach is a complex task due to difficulties of computational implementation and the correct definition of the interaction potentials. A convenient and widely used approach to describe anisotropic liquids employs the average description of the molecular orientation. The theory most widely used to describe liquid crystalline system is that of Ericksen-Leslie, in which mass and linear momentum balance are coupled to angular momentum balance.^{5-7} The latter is basically a torque balance in which viscous torques are balanced by elastic torques. The torque balance equation describes the spatio-temporal evolution of the average molecular orientation under the action of externally imposed flow fields. Figure 1b shows a schematic of the discotic nematic liquid crystalline phase, where the unit normals to the discs-like molecules (**u**) orient more or less parallel to the average molecular orientation, known as the director **n**, which is a unit vector: **n**.**n**=1.

Since discotic mesophases have orientational order, the average molecular orientation is affected by external flow, a mechanism known as flow-induced orientation.^{5-7} These materials are also elastic, such that spatial gradients of the average molecular orientation increases the energy. The three basic elastic storage modes for discotic nematics are splay, twist, and bend, and are shown in Figure 2.^{5-7} Therefore, any imposed processing flow creates through the balance between flow-induced orientation and elastic torques a unique texture that is defined by spatial characteristics of the director field such as orientation gradients, boundary layers, topological defects, and defect networks. The characterization and identification of liquid crystallinity is evaluated by polarized optical microscopy, since different textures correspond to different optical outputs.^{1-4}

Due to the practical importance of rheology in the fiber spinning process, and the eventual cost reduction, the experimental characterization of isotropic and mesophase pitches using steady and transient simple shear and Poiseuille flows has received some attention since the 1980s.^{8-17} The theoretical analyses of the experimental data, using constitutive equations appropriate to liquid crystalline systems, has received much less attention despite the fact that it is crucial to characterize the behavior of this complex material.^{18-23}

Small amplitude oscillatory shear flow is a main rheological tool used to characterize the viscoelaticity in terms of storage modulus (G'), loss (G") modulus, and loss angle (d) as a function of frequency (w). The characterization of isotropic and mesophase pitches using this technique has recently received some attention.^{24-30} The main features of the experimental measurements in textured carbonaceous mesophase are summarized in Figure 3, which shows schematics of the master curves for the frequency dependence of the storage G' and loss G" linear viscoelastic moduli using the Williams-Landel-Ferry (WLF) superposition.^{11,25-27} Figure 4 shows an schematic of a 2D view of a textured carbonaceous mesophases displaying the polydomain state;^{7} H is the average domain size within which a more or less uniform **n** is found. The key observations that can be extracted from experimental data (for the reduced frequency less than 100(rad/sec)) are: (*i*) The loss modulus scales with G"~w and is independent of temperature; (*ii*) the storage modulus has a temperature-sensitive terminal zone (small w) exhibiting solid-like behavior followed by a temperature insensitive viscoelastic, where G'~w^{2}.

The experimental results can be summarized as follows:^{24-27 }G" = G" (w, á**n**ñ), G' = G' (w, T, á**n**ñ, H),

where á**n**ñ is an average orientation that reflects splay, twist, and bend deformation modes.

Previous theoretical work^{31} on linear viscoelasticty was able to predict the frequency scalings (G'~w^{2}, G"~w) using a single mode (*i.e.* pure splay) anisotropic model. In this paper we present a generalized model that takes into account splay and bend distortions, and include the role of defects and textures, thus explaining the solid-like behavior in the terminal zone, as shown in the experiments.^{24-27} Hence in this paper we present a mesoscopic model that includes the required information of director orientation á**n**ñ and texture length scale H on the linear viscoelastic moduli. It should be mentioned that twist distortions have contribution to G' and their contribution to G" has the same frequency dependency as splay and bend.^{31} Hence there is no need to include the twist mode.

This paper is organized as follows: Section 2 presents the governing equations used to describe liquid crystalline system. The section 3 presents the results of the linear viscoelasticity for anisotropic textured carbonaceous mesophase. Section 4 presents the conclusions.

**Methods**

*Governing equations*

The continuum theory of elasticity of liquid crystals takes into account external forces and torques that introduce deformations in the relative molecular orientations and can distort the equilibrium configurations of the molecules. The elastic free energy density, F_{d,} a nematic uniaxial liquid crystal material is given by:

where the three basic modes of elastic storage are the splay (K_{11}), twist (K_{22}) and bend (K_{33}) modes, as mentioned above (see Figure 2).

The continuum theory of uniaxial nematic liquids consists of the linear and angular momentum balances, and constitutive equations for the stresses, viscous and elastic torques that takes into account external forces distort the spatially uniform equilibrium configurations of liquid crystals molecules. For incompressible isothermal conditions the general conservation of linear and angular momentum are given by the follows equations:^{5-7}

where r is the density,** v** is the velocity vector, **f** is the body force per unit volume vector, **s** is the total stress, r_{1} is the moment of inertia per unit volume, **G **is the external director body force vector, **g** is the intrinsic director body force vector, and** p** is the director stress tensor. The following constitutive equations for the stress tensor and the director body force that describes anisotropic liquids was found using transversely isotropic tensor coefficients, which reflect the material symmetry:

where,

p is the pressure, **I** is the unit tensor, {a_{i} : i=1, 2, 3, 4, 5, and 6}, are the Leslie viscosity coefficients that describes an anisotropic liquid, **A** is the rate of deformation tensor, **N** is the corotational derivative of the director vector, **b** is an arbitrary vector, "a" is an arbitrary scalar used to constrain the director (**n**) to be a unit vector, F_{d} is the elastic free energy density known as Frank elasticity, g_{1 }is the rotational viscosity, g_{2} is the irrotational torque coefficient, **W** is the vorticity tensor, and l is the reactive parameter.

Experiments show that the average molecular orientation of mesophase pitch on solid surfaces can be either planar (edge-on) or homeotropic (face-on) depending of the nature of the substrate.^{32,33} The edge-on anchoring (see Figure 5) is the more frequently found, and is the result of strong non-covalent interactions involving aromatic p clouds.^{32,33} In the analytical equations derived in our previous studies to describe the orientation field in capillary flow the average molecular orientation at the boundaries are assumed to be planar, and this also will be the case in this study.^{20,31,34-36}

A detailed discussion about the viscoelastic properties involved in the small amplitude oscillatory shear of discotic nematics can be found in our previous papers,^{20,31,34-36} here we briefly present the essential features of the material properties. The main properties include: (*i*) the Frank elastic moduli with the three distortion modes, splay (K_{11}), twist (K_{22}), and bend (K_{33}) (see equation 1); (*ii*) the Miesowicz shear viscosities that characterize viscous anisotropy and are measured in a steady simple shear flow between parallel plates with fixed director orientations along three characteristic orthogonal directions (h_{1} : when the director is parallel to the velocity direction, h_{2} : when it is parallel to the velocity gradient, and h_{3} : when it is parallel to the vorticity axis); (*iii*) the reactive parameter (see equation 8c) and the flow-alignment angle (q_{s }= 1/2 cos^{-1} (1/l)) that obtains when l<-1; (*iv*) the director re-orientation viscosities associated with splay, twist, and bend deformations (h_{tiwst}= g_{1}, h_{splay }= g_{1 }- a_{3}^{2} / h_{1}, and h_{bend}= g_{1 }- a_{2}^{2} / h_{2}), which are given by the rotational viscosity (g_{1}) decreased by a factor introduced by the backflow effect.^{7} Backflow is re-orientation driven flow and is essentially the reverse effect to flow-induced orientation.^{7} Note that twist is the only mode that creates no backflow.^{7}

**Results and Discussion**

*Linear viscoelasticity of monodomain and polydomain carboanceous mesophase*

The solution of equations 1 to 8c for simple shear flow and capillary Poiseuille flow, assuming linear behavior that happens at small deformation regime as in the case of small amplitude oscillatory flows, allows the derivation of viscoelastic functions, namely the storage modulus (G'), the loss modulus (G"), and loss tangent ( tan d = ):^{31,34-36}

where U_{1} and U_{2} are frequency-dependent functions, and P_{1} and P_{2} are the material properties ratios that assume distinct values for the distinct distortion modes (*i.e.* splay, bend, and twist) as shown in Table 1. As mentioned above, the twist mode is not included because this mode does not contribute to elastic storage. Since the material parameters {P_{1}(T,**n**), P_{2}(T,**n**)} are functions of (T,**n**) (see Table 1),^{31,36} the linear viscoelastic functions (equations 9 to 11) are also functions of (T,**n**). In particular, when **n** is along the flow direction the distortion mode bend and when is along the shear rate direction, the distortion mode is pure splay.

In this paper the calculations are performed using a set of estimated viscoelastic material parameters listed in Table 2.^{15,18} The characteristic size (L) is the gap between plates or the capillary radius in the rheometer. The essential features of the solutions and the main conclusions of this paper will remain unaffected if other parameters are used but the reactive parameter (l) is less than minus one and if the Miesowicz inequalities hold.^{20-23}

Figures 6a and 6b shows the dynamic moduli (G" and G') and the phase lag as a function of the frequency (w) for a single distortion mode (splay and bend, respectively) of a monodomain (defect free). In both splay and mode cases, the low frequency (terminal) regime is classic of a viscous fluid, and the loss modulus is always greater than the storage modulus. The characteristic slopes in both cases are:

These single mode results capture features of the experimental data shown in Figure 3 but miss the temperature dependence of G', as well as the solid-like terminal frequency plateau of G'. On the other hand, the single mode model reproduces the experimentally observed G" response. Since the slope of G" is insensitive to the deformation mode, there is no need to developed a multimode equation. Hence we conclude from the present theory (Figures 6a and 6b) and from experiments (Figure 3) that G" is insensitive to the distortion mode and to textural processes. Next we focus attention to G' and develop a model based on liquid crystal physics that reproduces the features of G' shown in Figure 3.

In this study the storage modulus G' of textured carbonaceous mesophase, with splay, and bend deformation domains is estimated using the aggregate model, which is widely used in stress analysis of anisotropic materials.^{27} In the present paper the polydomain material behavior is assumed to be given by weighted sum of both modes and a defect contribution:

where G'_{texture} is the storage modulus of the textured sample, G'_{bend} is the storage modulus for the monodomain in bend distortion, G'_{splay} is the storage modulus for the monodomain in splay distortion, G'_{defect} is the elasticity stored by the orientation defects (disclinations) that create the texture, and j is the volume fraction of each distortion mode. Equation 14 is derived using the uniform stress principle of textured materials, and it is reasonably convenient to represent the mesophase pitch in the multi-distortion mode. In this study we assume that the mesophase pitch has an equal amount of domains that will be distorted in bend and splay modes, therefore j_{splay} = j_{bend} = 0.5.

The elastic storage of orientational defects (disclination) of nematics liquid crystals is a well known phenomena.^{5,6,38,42-44} The amount of elastic storage depends on the domain size (H), the defect charge (S), and the temperature, as follows:^{5,6,38}

where c is a constant, K(T) is the characteristic Frank elastic constant or average of splay, bend, and twist. In the present study we take into account the defect elastic storage for three temperatures, T_{1}<T_{2}<T_{3}. As well known, K(T) is a decreasing function of increasing temperature, equation 15 predicts that G'_{defect} decreases with increasing temperature. No exact temperature is indicated in this paper because no comprehensive data exist for the temperature dependency of K for carbonaceous mesophases. The selected K values are equal to 10^{-7 }Pa (T=T_{3}), 10^{-6} Pa (T=T_{2}), and 10^{-5 }Pa (T=T_{1}), respectively, which imply that domain size is in the micron range, as observed experimentally.^{44}

Figure 7 shows the computed storage modulus (G'_{mix}) as a function of frequency (w (rad/sec)) of a textured carbonaceous mesophase in the small frequency regime at three temperatures (T_{1}<T_{2}<T_{3}), computed using equation 14. The figure shows a plateau in the terminal zone (w < 0.1), whose amplitude decreases with increasing temperature, in perfect agreement with Figure 3. The solid like behavior is due to the domain texture, whose elastic solid-like nature is due to the presence of orientational defects.^{44}

**Conclusions**

This paper extends the theoretical methods used to describe the linear viscoelasticity of anisotropic textured carbonaceous mesophases. Previous models focused on single deformation modes in the absence of defects. It is well known that carbonaceous mesophases are textured materials, that when sheared store elasticity due to orientation gradients and orientational defects. A linear viscoelasticity model based on the constant stress principle used to described anisotropic textured materials is developed based on liquid crystal physics. The model predict that at low frequencies the loss modulus is insensitive to temperature, textures and deformation modes, since dissipative mechanisms are insensitive to the presence of defects. On the other hand, the model predicts that the storage modulus has a strong sensitivity to temperature, and texture. The solid like behavior observed in the terminal zone is attributed to topological defects. The model predictions are in excellent agreements with experiments. This work provides a new way to incorporate textural information into linear viscoelasticity and hence provides new tools to understand and use experimental linear viscoelastic data. As shown by Equation 15 the solid like behavior displayed by G' can easily be converted into a texture scale (H) or into a Frank elastic constant.

**Acknowledgments**

This research was supported by a grant from Engineering Research Center program of National Science Foundation (USA) under award number EEC 9731680.

**References**

1. Singer, L. S.; *Fuel* **1981**, *60*, 839. [ Links ]

2. Marsh, H.; Latham, C. S.; *ACS Symposium Series* **1986**, *303*, 1. [ Links ]

3. Edie, D. D.; Diefendorf, R. J. In *Carbon-Carbon Materials and Composites*; Buckley, J. D; Edie, D. D., eds.; Noyes Publications, 1993, p. 19. [ Links ]

4. Edie, D. D. In *Design and Control of Structure of Advanced Carbon Materials for Enhanced Performance*; Rand, B.; Appleyard, S. P.; Yardim, M. F., eds.; Kluwer Academic Publishers: Netherlands, 2001, p. 163. [ Links ]

5. Chandrasekhar, F. R. S.; *Liquid Crystals*, 2^{nd} ed., University Press: Cambridge, 1992. [ Links ]

6. de Gennes, P. G.; Prost, J.; *The Physics of Liquid Crystals*, 2^{nd} ed., Oxford University Press: London, 1993. [ Links ]

7. Rey, A. D.; Denn, M. M.; *Annu. Rev. Fluid Mech.* **2002**, *34*, 233. [ Links ]

8. Nazem, F. F.; *Fuel*** 1980**, *59*, 851. [ Links ]

9. Sakai, M.; Inagaki, M.; *Carbon* **1981**, *19*, 37. [ Links ]

10. Nazem, F. F.; *Carbon* **1982**, *20*, 345. [ Links ]

11. Nazem, F. F.; Lewis I. C.; *Mol. Cryst. Liq. Cryst.* **1986**, *139*, 195. [ Links ]

12. Oh, K. H; Ryu, S. K.; Rhee, B. S.; *Korean J. Chem. Eng.* **1986**, *3*, 93. [ Links ]

13. Rand, B.; *Fuel* **1987**, *66*, 1491. [ Links ]

14. Fathollahi, B.; White, J. L.; *J. Rheol. ***1994**, *38*, 1591. [ Links ]

15. Fleurot O.; Edie, D. D.; *J. Rheol.* **1998**, *42*, 781. [ Links ]

16. Cato, A. D.; Edie, D. D.; *Carbon* **2003**, *41*, 1411. [ Links ]

17. Cato, A. D.; Edie, D. D.; Harrison, G. M.; *J. Rheol.*** 2005**, *49*, 161. [ Links ]

18. Grecov, D.; Rey, A. D.; *Mol. Cryst. Liq. Cryst.* **2003**, *391*, 57. [ Links ]

19. Grecov, D.; Rey, A. D.; *J. Rheol. ***2005**, *49*, 175. [ Links ]

20. de Andrade Lima, L. R. P.; Rey, A. D.; *J. Non-Newton. Fluid Mech.* **2003**, *110*, 103. [ Links ]

21. de Andrade Lima, L. R. P.; Rey, A. D.; *J. Rheol.* **2003**, *47*, 1261. [ Links ]

22. de Andrade Lima, L. R. P.; Rey, A. D.; *Mater. Res.* **2003**, *6*, 285. [ Links ]

23. de Andrade Lima, L. R. P.; Rey, A. D.; *Carbon* **2004**, *42*, 1263. [ Links ]

24. Sakai, M.; Sato, Y.; Inagaki, M.; *Carbon'90: Extended Abstracts and Program – International Carbon Conference,* Paris, 1990, p. 298. [ Links ]

25. Hara, R.; *Carbon'90: Extended Abstracts and Program – International Carbon Conference*, Paris, 1990, p. 290. [ Links ]

26. Hara, R.; *Mitsubishi Kasei R&D Review* **1992**, *6*, 19. [ Links ]

27. Hara, R.; *Nihon Reoroji Gakkaishi* **1993**, *21*, 74. [ Links ]

28. Turpin, M.; Cheung, T.; Rand, B.; *Carbon* **1994**, *32*, 225. [ Links ]

29. Cheung, T.; Turpin, M.; Rand, B.; *Carbon* **1996**, *34*, 265. [ Links ]

30. Wu, X. X.; Wang, C. Y.; Guo, C. T.; *Carbon* **1998**, 36, 1291. [ Links ]

31. de Andrade Lima, L. R. P.; Rey, A. D.; *Chem. Eng. Sci.* **2004**, 59, 3891. [ Links ]

32. Hurt, R.; Krammer, G.; Crawford, G.; Jian, K. Q.; Rulison, C.; *Chem. Mater.* **2002**, *14*, 4558. [ Links ]

33. Jian, K.; Shim, H-S.; Tuhs-Dubrow, D.; Bernstein, S.; Woodward, C.; Pfeffer, M.; Steingart, D.; Gournay, T.; Sachsmann, S.; Crawford, G. P.; Hurt, R. H.; *Carbon* **2003**, *41*, 2073. [ Links ]

34. de Andrade Lima, L. R. P.; Rey, A.D.; *Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys.* **2004**, *70*, 011701-1. [ Links ]

35. de Andrade Lima, L. R. P.; Rey, A.D.; *J. Rheol.* **2004**, *48*, 1067. [ Links ]

36. de Andrade Lima, L. R. P.; Rey, A.D.; *Rheol. Acta* **2006**, in press. [ Links ]

37. Hyer, M. W.; *Stress Analysis of Fiber-Reinforced Composite Materials*, WCB/McGraw-Hill: New York 1998. [ Links ]

38. Chandrasekhar, S.; Ranganath, G. S.; *Adv. Phys.* **1986**, *35*, 507. [ Links ]

39. Rey, A. D.; *Rheol. Acta* **1995**, *34*, 461. [ Links ]

40. Rey, A. D.; *Mol. Cryst. Liq. Cryst.* **1996**, *281*, 155. [ Links ]

41. Rey, A. D.; *Liq. Cryst. ***1996**, *20*, 147. [ Links ]

42. Currie, P. K; Leslie, F. M.; *Arch. Ration. Mech. Anal.* **1984**, *87*, 35. [ Links ]

43. Yada, M.; Yamamoto, J.; Yokoyama, H.; *Langmuir* **2003**, *19*, 3650. [ Links ]

44. Zimmer, J. E.; White, J. L.; *Adv. Liq. Cryst.* **1982**, *5*, 157. [ Links ]

Received: January 22, 2006

Published on the web: June 29, 2006

* e-mail: lelo@ufba.br