Transient elongational viscosities of aqueous polyacrylamide solutions measured with an optical rheometer

An anionic polyacrylamide solution was characterized in elongational flow by combining laser-Doppler velocimetry to determine the strain rate in the flow direction and the two-color flow-induced birefringence method to measure the first normal stress difference along the axial centerline of a hyperbolic die. The elongational rate was constant along the axial centerline of the planar hyperbolic die as long as vortices at the die entrance did not occur. The transient elongational viscosity μ ⁺ was determined as a function of the elongational rate. The parameters varied are the Hencky strain rate and the polymer concentration. μ ⁺ showed a pronounced increase over the linear viscoelastic behavior above critical Hencky strains of 1.2 to 1.5; that is, a significant strain hardening could be observed for polyacrylamide solutions. This strain hardening is stronger the higher the elongational rate. A slight enhancement of strain hardening was found by increasing the concentration from 0.5 to 1 g/l. The stress optical coefficient was determined as 1.8 × 10⁻⁷ Pa⁻¹ (0.5 g/l) and 1.2 × 10⁻⁷ Pa⁻¹ (1 g/l).

Flow-induced crystallization of high-density polyethylene: the effects of shear and uniaxial extension

The effects of shear, uniaxial extension and temperature on the flow-induced crystallization of two different types of high-density polyethylene (a metallocene and a ZN-HDPE) are examined using rheometry. Shear and uniaxial extension experiments were performed at temperatures below and well above the peak melting point of the polyethylenes in order to characterize their flow-induced crystallization behavior at rates relevant to processing (elongational rates up to 30 s^− 1 and shear rates 1 to 1,000 s^− 1 depending on the application). Generally, strain and strain rate found to enhance crystallization in both shear and elongation. In particular, extensional flow was found to be a much stronger stimulus for polymer crystallization compared to shear. At temperatures well above the melting peak point (up to 25°C), polymer crystallized under elongational flow, while there was no sign of crystallization under simple shear. A modified Kolmogorov crystallization model (Kolmogorov, Bull Akad Sci USSR, Class Sci, Math Nat 1:355–359, 1937) proposed by Tanner and Qi (Chem Eng Sci 64:4576–4579, 2009) was used to describe the crystallization kinetics under both shear and elongational flow at different temperatures.     

Large amplitude oscillatory elongation flow

A filament stretching rheometer (FSR) was used for measuring the elongation flow with a large amplitude oscillative elongation imposed upon the flow. The large amplitude oscillation imposed upon the elongational flow as a function of the time t was defined as where ε is the Hencky strain, is a constant elongational rate for the base elongational flow, Λ the strain amplitude (Λ ≥ 0), and Ω the strain frequency. A narrow molecular mass distribution linear polystyrene with a molecular weight of 145 kg/mol was subjected to the oscillative flow. The onset of the steady periodic regime is reached at the same Hencky strain as the onset of the steady elongational viscosity ( Λ = 0). The integral molecular stress function formulation within the ‘interchain pressure’ concept agrees qualitatively with the experiments.     

Reynolds Number Effects on the Coherent Dynamics of the Turbulent Horseshoe Vortex System

The adverse pressure gradient induced by a surface-mounted obstacle in a turbulent boundary layer causes the approaching flow to separate and form a dynamically rich horseshoe vortex system (HSV) in the junction of the obstacle with the wall. The Reynolds number of the flow (Re) is one of the important parameters that control the rich coherent dynamics of the vortex, which are known to give rise to low-frequency, bimodal fluctuations of the velocity field (Devenport and Simpson, J Fluid Mech 210:23–55, 1990; Paik et al., Phys Fluids 19:045107, 2007). We carry out detached eddy simulations (DES) of the flow past a circular cylinder mounted on a rectangular channel for Re = 2.0 × 10⁴ and 3.9 × 10⁴ (Dargahi, Exp Fluids 8:1–12, 1989) in order to systematically investigate the effect of the Reynolds number on the HSV dynamics. The computed results are compared with each other and with previous experimental and computational results for a related junction flow at a much higher Reynolds number (Re = 1.15 × 10⁵) (Devenport and Simpson, J Fluid Mech 210:23–55, 1990; Paik et al., Phys Fluids 19:045107, 2007). The computed results reveal significant variations with Re in terms of the mean-flow quantities, turbulence statistics, and the coherent dynamics of the turbulent HSV. For Re = 2.0 × 10⁴ the HSV system consists of a large number of necklace-type vortices that are shed periodically at higher frequencies than those observed in the Re = 3.9 × 10⁴ case. For this latter case the number of large-scale vortical structures that comprise the instantaneous HSV system is reduced significantly and the flow dynamics becomes quasi-periodic. For both cases, we show that the instantaneous flowfields are dominated by eruptions of wall-generated vorticity associated with the growth of hairpin vortices that wrap around and disorganize the primary HSV system. The intensity and frequency of these eruptions, however, appears to diminish rapidly with decreasing Re. In the high Re case the HSV system consists of a single, highly energetic, large-scale necklace vortex that is aperiodically disorganized by the growth of the hairpin mode. Regardless of the Re, we find pockets in the junction region within which the histograms of velocity fluctuations are bimodal as has also been observed in several previous experimental studies.     

Direct comparison of equibiaxial elongational viscosity measurements from lubricated squeezing flow and the MultiAxiales Dehnrheometer

We present the first direct comparisons of rheological data from the lubricated squeezing flow (LSF) technique and the MultiAxiales Dehnrheometer (MAD) instrument developed by Meissner and coworkers (J Rheol 47:989–1010, 2003). Comparisons of transient equibiaxial elongational viscosity are carried out at strain rates well into the nonlinear regime on low-density polyethylene and polystyrene melts. We find data obtained using LSF deviate from the MAD data when the Hencky strain reaches a value of approximately 1, which we interpret as a failure of the LSF technique. The strain at which the LSF technique fails is relatively insensitive to experimental parameters including strain rate. For Hencky strains larger than 1, LSF data display behavior that could easily be mistaken for the phenomenon of strain hardening.     

Shear and elongational flow behavior of acrylic thickener solutions

We investigate the effect of hydrophobic aggregation in alkali-swellable acrylic thickener solutions on shear and extensional flow properties at technically relevant polymer concentrations using the commercial thickener Sterocoll FD as model system. Apparent molecular weight of aggregates in water is M _w  ≈ 10⁸ g/mol and decreases by more than an order of magnitude in ethanol. Zero shear viscosity η ₀ is low and shear thinning is weak compared to the high molecular weight of the aggregates. Linear viscoelastic relaxation is described by the Zimm theory up to frequencies of 10⁴ rad/s, demonstrating that no entanglements are present in these solutions. This is further supported by the concentration dependence of η ₀ and is attributed to strong association within the aggregates. Extensional flow behavior is characterized using the capillary break-up extensional rheometry technique including high-speed imaging. Solutions with ϕ ≥ 1% undergo uniform deformation and show pronounced strain hardening up to large Hencky strains. Elongational relaxation times are more than one order of magnitude lower than the longest shear relaxation times, suggesting that aggregates cannot withstand strong flows and do not contribute to the elongational viscosity.

Modelling elongational and shear rheology of two LDPE melts

Experimental data of two low-density polyethylene (LDPE) melts at 200°C for both shear flow (transient and steady shear viscosity as well as transient and steady first normal stress coefficient) and elongational flow (transient and steady-state elongational viscosity) as published by Pivokonsky et al. (J Non-Newtonian Fluid Mech 135:58–67, 2006) were analysed using the molecular stress function model for broadly distributed, randomly branched molecular structures. For quantitative modelling of melt rheology in both types of flow and in a very wide range of deformation rates, only three nonlinear viscoelastic material parameters are needed: Whilst the rotational parameter, a ₂, and the structural parameter, β, are found to be equal for the two melts considered, the melts differ in the parameter describing maximum stretch of the polymer chains.     

Hydrodynamic diffusivity of spherical particles in polymer solution

In this research experiments were performed to examine the hydrodynamic diffusion of spherical particles in a highly filled suspension. The suspension consisted of nearly monodisperse polymethylmethacrylate spheres in a density matched polymer solution. The polymer solution was prepared by dissolving 0–700 ppm of polyacrylamide in a mixture of ethyleneglycol and glycerine. The polymer solution did not show appreciable shear thinning. The particle loading was varied from 30 to 55%. The hydrodynamic diffusivity was estimated by measuring the time-dependent viscosity when the suspension was subjected to a circular Couette flow with an air bubble trapped under the rotor of the Couette apparatus. The results show that the dimensionless diffusivity (D/γ˙a ²) of particles in polymer solution is not proportional to shear rate (γ˙), as in the case of a Newtonian fluid, but that it decreases with increasing shear rate. The diffusivity also decreases with increasing polymer concentration. It is suggested that the elongational thickening behaviour and the increased lubrication force due to the first normal stress difference may be responsible for the reduction of diffusivity in the polymer solution. Received: 18 January 2000 Accepted: 6 April 2000     

Uniaxial elongational behavior of poly(vinyl chloride) physical gel

A poly(vinyl chloride) (PVC,  M_w = 102×10³)(\mbox{PVC,}\;{\rm M}_{\rm w} =102\times 10^3) di-octyl phthalate (DOP) gel with PVC content of 20 wt.% was prepared by a solvent evaporation method. The dynamic viscoelsticity and elongational viscosity of the PVC/DOP gel were measured at various temperatures. The gel exhibited a typical sol–gel transition behavior with elevating temperature. The critical gel temperature (T_gel) characterized with a power–law relationship between the storage and loss moduli, G^′ and G^″, and frequency ω, G￠=G^￠￠/tan  ( np/2 ) μ wⁿ{G}^\prime={G}^{\prime\prime}{\rm /tan}\;\left( {{n}\pi {\rm /2}} \right)\propto \omega ^{n}, was observed to be 152°C. The elongational viscosity of the gel was measured below the T_gel. The gel exhibited strong strain hardening. Elongational viscosity against strain plot was independent of strain rate. This finding is different from the elongational viscosity behavior of linear polymer solutions and melts. The stress–strain relations were expressed by the neo-Hookean model at high temperature (135°C) near the T_gel. However, the stress–strain curves were deviated from the neo-Hookean model at smaller strain with decreasing temperature. These results indicated that this physical gel behaves as the neo-Hookean model at low cross-linking point, and is deviated from the neo-Hookean model with increasing of the PVC crystallites worked as the cross-linking junctions.     

Three-dimensional Turbulent Boundary Layer in a Shrouded Rotating System

A thre-dimensional direct numerical simulation is combined with a laboratory study to describe the turbulent flow in an enclosed annular rotor-stator cavity characterized by a large aspect ratio G = (b − a)/h = 18.32 and a small radius ratio a/b = 0.152, where a and b are the inner and outer radii of the rotating disk and h is the interdisk spacing. The rotation rate Ω considered is equivalent to the rotational Reynolds number Re = Ωb ²/ν= 9 .5 × 10⁴ (ν the kinematic viscosity of water). This corresponds to a value at which experiment has revealed that the stator boundary layer is turbulent, whereas the rotor boundary layer is still laminar. Comparisons of the computed solution with velocity measurements have given good agreement for the mean and turbulent fields. The results enhance evidence of weak turbulence by comparing the turbulence properties with available data in the literature (Lygren and Andersson, J Fluid Mech 426:297–326, 2001). An approximately self-similar boundary layer behavior is observed along the stator. The wall-normal variations of the structural parameter and of characteristic angles confirm that this boundary layer is three-dimensional. A quadrant analysis (Kang et al., Phys Fluids 10:2315–2322, 1998) of conditionally averaged velocities shows that the asymmetries obtained are dominated by Reynolds stress-producing events in the stator boundary layer. Moreover, Case 1 vortices (with a positive wall induced velocity) are found to be the major source of generation of special strong events, in agreement with the conclusions of Lygren and Andersson (J Fluid Mech 426:297–326, 2001).     

Small and large strain rheology of wheat gluten

 The rheological properties of wheat gluten were studied under both small and large deformation and compared with those of the parent flours. The limiting strain of linear viscoelastic behaviour of gluten doughs, 3 × 10⁻², was an order of magnitude larger than that of the flour doughs, 10⁻³. The role of starch in the lower limiting strain of flour doughs was indicated by the exponential decrease in the limiting strain of gluten-starch mixtures with greater quantities of starch. Large strain measurements showed gluten doughs possessed greater shear and elongational viscosities than flour doughs and these differences were greatest at lower shear and elongation rates (0.01 and 0.1 s⁻¹). The larger viscosities of flour and gluten doughs at the low strain rates help to stabilise and prevent the collapse of gas bubbles during bread fermentation and baking. Increasing starch levels in gluten-starch mixtures, at either constant or optimal water levels, lowered the elongational viscosity. Dynamic measurements were, however, more sensitive to the level of water added to the gluten-starch mixtures. The storage modulus decreased with increasing starch levels when constant water levels were used to prepare the mixtures, but when optimal water levels were used the storage modulus increased. Gluten and starch are major contributors to the large and small strain rheological properties of flour doughs; however, gluten-starch mixtures were unable to duplicate exactly the rheological properties of flour doughs, indicating that other flour components such as pentosans, lipids and water soluble proteins also influence dough rheology. Received: 20 March 2001 Accepted: 11 July 2001     

Grafting of polyamide 6 on a styrene–acrylonitrile maleic anhydride terpolymer: melt rheology at the critical gel state

In this work, we studied the melt rheology of multigraft copolymers with a styrene–acrylonitrile maleic anhydride (SANMA) terpolymer backbone and randomly grafted polyamide 6 (PA 6) chains. The multi-grafted chains were formed by interfacial reactions between the maleic anhydride groups of SANMA and the amino end groups of PA 6 during melt blending. Because of the phase separation of SANMA and PA 6, the grafted SANMA backbones formed nearly circular domains which were embedded in the PA 6 melt with a diameter in the order of 20 to 40 nm. The linear viscoelastic behaviour of PA 6/SANMA blends at a sufficiently large SANMA concentration displayed the characteristics of the critical gel state, i.e. the power relations G′ ∝ G′′ ∝ ω ^0.5. In elongation, the PA 6/SANMA blend at the critical gel state showed a non-linear strain hardening behaviour already at a very small Hencky strain. In contrast to neat PA 6, the elasticity of the PA 6/SANMA blends was strongly pronounced, which was demonstrated by recovery experiments. Rheotens tests agreed with the linear viscoelastic shear oscillations and the measurements using the elongational rheometer RME. Increasing the SANMA concentration led to a larger melt strength and a reduced drawability. The occurrence of the critical gel state can be interpreted by the cooperative motion of molecules which develops between the grafted PA 6 chains of neighbouring micelle-like SANMA domains.     

The MSF model: relation of nonlinear parameters to molecular structure of long-chain branched polymer melts

The elongational viscosity data of model PS combs (Hepperle J, Einfluss der Molekularen Struktur auf Rheologische Eigenschaften von Polystyrol- und Polycarbonatschmelzen. Doctoral Thesis, University Erlangen-Nürnberg, 2003) are reconsidered by including the interchain pressure term of Marrucci and Ianniruberto [Macromolecules 37:3934–3942, 2004] in the Molecular Stress Function model [Wagner et al., J Rheol 47(3):779–793, 2003, Wagner et al., J Rheol 49:1317–1327, 2005d]. Two nonlinear model parameters are needed to describe elongational flow, β and . The parameterβ determines the slope of the elongational viscosity after the inception of strain hardening. It is directly related to the molecular structure of the polymer and represents the ratio of the molar mass of the (branched) polymer to the molar mass of the backbone alone. β follows from the hypothesis of Wagner et al. [J Rheol 47(3):779–793, 2003] that side chains are compressed onto the backbone. We consider also the case that side chains are oriented by deformation, but not stretched, and found little difference in the model predictions. The parameter represents the maximum strain energy stored in the polymeric system and determines the steady-state value of the viscosity in extensional flows. The relation of this energy parameter to the molecular structure is discussed. Good correlations between the energy parameter and different coil contraction ratios, as determined either experimentally or calculated theoretically by considering the topology of the macromolecule, are found. The smaller the relative size of the polymer coil, the larger is the energy parameter and the more strain energy can be stored in the polymeric system. Presented at the 3rd Annual European Rheology Conference, AERC2006, Crete, Greece.     

Evaluation of drag correction factor for spheres settling in associative polymers

Drag correction factors are calculated for the creeping motion of spheres descending in various associative polymers of different concentration with various sphere-container ratios and Weissenberg numbers. The simple-shear rheology and linear viscoelasticity of these polymeric fluids have been previously presented and modeled with the BMP (Bautista–Manero–Puig) equation of state (Mendoza-Fuentes et al., Phys Fluids 21:033104, 2009). The drag on the sphere is initially kept nearly constant for small Weissenberg numbers, We < 0.1. As the Weissenberg number increases, We < 0.1, a reduction in drag is found. Experimental results show the presence of a critical Weissenberg number at which a drag reduction occurs. The reduction in the drag correction factor is associated to the onset of extension-thinning, which coincides with the formation of a negative wake. No increase in the drag correction factor was observed, due to the simultaneous opposing effects of extension-thickening and shear-thinning viscosity. The shape of the drag correction factor curve may be predicted considering the extensional properties of the solutions, as suggested elsewhere (Chen and Rothstein, J Non-Newton Fluid Mech 116:205–215, 2004).     

The interchain pressure effect in shear rheology

Recently, the tube diameter relaxation time in the evolution equation of the molecular stress function (MSF) model (Wagner et al., J Rheol 49: 1317–1327, 2005) with the interchain pressure effect (Marrucci and Ianniruberto, Macromolecules 37:3934–3942, 2004) included was shown to be equal to three times the Rouse time in the limit of small chain stretch. From this result, an advanced version of the MSF model was proposed, allowing modeling of the transient and steady-state elongational viscosity data of monodisperse polystyrene melts without using any nonlinear parameter, i.e., solely based on the linear viscoelastic characterization of the melts (Wagner and Rolón-Garrido 2009a, b). In this work, the same approach is extended to model experimental data in shear flow. The shear viscosity of two polybutadiene solutions (Ravindranath and Wang, J Rheol 52(3):681–695, 2008), of four styrene-butadiene random copolymer melts (Boukany et al., J Rheol 53(3):617–629, 2009), and of four polyisoprene melts (Auhl et al., J Rheol 52(3):801–835, 2008) as well as the shear viscosity and the first and second normal stress differences of a polystyrene melt (Schweizer et al., J Rheol 48(6):1345–1363, 2004), are analyzed. The capability of the MSF model with the interchain pressure effect included in the evolution equation of the chain stretch to model shear rheology on the basis of linear viscoelastic data alone is confirmed.     

Effect of melt viscosity of polypropylene on fibrillation of thermotropic liquid crystalline polymer in in situ composite film

 Various grades of polypropylene were melt blended with a thermotropic liquid crystalline polymer, a block copolymer of p-hydroxy benzoic acid and ethylene terephthalate (60/40 mole ratio). The blends were extruded as cast films at different values of draw ratio (slit width/film thickness). Fibrillation of TLCP dispersed phase with high fiber aspect ratio (length/width) was obtained with the matrix of low melt flow rate, i.e., high viscosity and with increasing film drawing. Melt viscosities of pure components and blends measured using capillary rheometer were found to decrease with increasing shear rate and temperature. Viscosity ratios (dispersed phase to matrix phase) of the systems being investigated at 255 °C at the shear rate ranged from 10² to 10⁴ s⁻¹, were found to lie between 0.04 and 0.15. The addition of a few percent of elastomeric compatibilizers; a tri-block copolymer SEBS, EPDM rubber and maleated-EPDM, was found to affect the melt viscosity of the blend and hence the morphology. Among these three compatibilizers, SEBS was found to provide the best fibrillation. Received: 10 January 2000/Accepted: 24 January 2000     

The effect of 1,3:2,4-<Emphasis Type="Italic">bis</Emphasis>-<Emphasis Type="Italic">O</Emphasis>-(<Emphasis Type="Italic">p</Emphasis>-methylbenzylidene)-<Emphasis Type="SmallCaps">d</Emphasis>-sorbitol (PDTS) on uniaxial elongational viscosity of polypropylene

We investigated the dynamic viscoelasticity and elongational viscosity of polypropylene (PP) containing 0.5 wt% of 1,3:2,4-bis-O-(p-methylbenzylidene)-d-sorbitol (PDTS). The PP/PDTS system exhibited a sol–gel transition (T _gel) at 193 °C. The critical exponent n was nearly equal to 2/3, in agreement with the value predicted by a percolation theory. This critical gel is due to a three-dimensional network structure of PDTS crystals. The elongational viscosity behavior of neat PP followed the linear viscosity growth function 3η ⁺ (t), where η ⁺ (t) is the shear stress growth function in the linear viscoelastic region. The elongational viscosity of the PP/PDTS system also followed the 3η ⁺ (t) above T _gel but did not follow the 3η ⁺ (t) and exhibited strong strain-softening behavior below T _gel. This strain softening can be attributed to breakage of the network structure of PDTS with a critical stress (σ _c) of about 10⁴ Pa.     

On modes and criteria of ABS melt failure in extension

Melt failure of a commercial ABS polymer in uniaxial extension over ranges of elongation rate ([(e)\dot] = 0.01 - 1.0 s ^{- 1}\dot \varepsilon = 0.01 - 1.0\,{\rm s}^{ - 1} ) and temperature (140-200 °C) was investigated. Four methods of experimental and numerical calculation for determination of modes and criteria of melt failure in uniaxial extension were investigated: 1) visual observation of necking; 2) visual observation of non-uniform flow during stress relaxation after cessation of steady elongation; 3) calculation of the Considère criterion from the measured elongational stress-strain curve; 4) numerical calculation of inflection point (‘C²/‘)²=0) from the tensile stress-strain curve. In addition, under higher Deborah number conditions the critical Hencky strains at Considère criterion were calculated using PSM model parameters (! and #) and were compared with those obtained from the measured elongational stress-strain curve. The relationship between these failure modes is discussed in terms of rheological properties of the polymer, putting emphasis on the relationship with the thermoforming process. The Considère criterion appears to be the most effective indicator of the non-uniform deformation of ABS melt in uniaxial extension under conditions where cohesive fracture does not occur. The rheological properties such as elongational viscosity, strain hardening and/or strain softening, and their temperature dependence play an important role in determining the growth and transition of melt failure of ABS polymer in uniaxial extension.     

Simultaneous determination of shear and extensional viscosities and wall slip by on-line rheometry with parameter fitting: an application to EPDM rubber

Shear and extensional viscosities and wall slip are determined simultaneously under extrusion processing conditions using an on-line rheometer. Because it is not possible to independently control flow rate and temperature, classical methods for interpretation of capillary data cannot be used with on-line rheometry. This limitation is overcome using computational optimization to fit parameters in a flow model. This consists of three parts, representing shear viscosity, extensional viscosity, and wall slip. Three-parameter, power law forms, based on local instantaneous deformation rates and including temperature dependence, are used for each, and analytic solutions applied for entry flow and flow in the capillary. For entry flow, the Cogswell–Binding approach is used, and for developed flow in the capillary a solution incorporating wall slip is derived. The rheometer, with interchangeable capillaries, is mounted in place of the die on a rubber profile extrusion line. Pressure drops and temperatures for extrusion of an EPDM rubber through 2 mm diameter capillaries of length 0, 2, 3, 4, and 5 mm are logged and flow rates determined for a range of extruder speeds (5 to 20 rpm). Pressures ranged from 60 to 75 bar and temperatures from 86 to 116 °C. Mean flow velocity in the capillaries was between 5 × 10⁻³ and 5 × 10⁻¹ m s⁻¹. The nine material parameters are optimized for best fit of the analytic pressure drops to experimental data, using about 100 data points, with the Levenberg–Marquardt method. It is concluded that flow is dominated by extension and wall slip. Shear flow appears to play little part. The slip model indicates that slip velocity increases much more rapidly than the wall shear stress (in the range 0.5–1 MPa) and decreases with temperature for a given stress level. Results for the (uniaxial) extensional viscosity represent an engineering approximation to this complex phenomenon at the high strains (approximately 200) and high extension rates (up to 800 s⁻¹) applying in the extrusion. Results indicate a slight extension hardening and a decrease with temperature. Results are put into the context of the available studies in the literature, which, particularly with regard to wall-slip and extensional flow, consider conditions far removed from those applying in industrial extrusion. The present methods provide a powerful means for flow characterization under processing conditions, providing data suitable for use in computer simulations of extrusion and optimization of die design.     

Rheology of SiO2/(acrylic polymer/epoxy) suspensions. I. Linear viscoelasticity

Linear viscoelastic properties of SiO₂/(AP/EP) suspension with various SiO₂ volume fractions (ϕ) in a blend of acrylic polymer (AP) and epoxy (EP) were investigated at various temperatures (T). The AP/EP contained 70 vol.% of EP. The SiO₂ particles were treated with epoxy silane coupling agent. The effects of the SiO₂ particles are more pronounced in the terminal zone: a transition from viscoelastic liquid (ϕ ≤ 30 vol.%) to viscoelastic solid (ϕ ≥ 40 vol.%) was observed which can be interpreted as a critical gelation occurring at a critical particle content and critical gel temperature. The SiO₂/(AP/EP) systems exhibited a critical gel behavior at ϕ ≅ 35 vol.% and T ≅ 100°C characterized with a power–law relationship between the storage and loss moduli (G ^′ and G ^″) and frequency (ω); G ^′ = G ^″/tan(nπ/2) ∝ ω ⁿ . The critical gel exponent (n) was estimated to be about 0.45. The gelation occurred with increasing T.