Examination of flow birefringence induced by the shear components along the optical axis using a parallel-plate-type rheometer | Scientific Reports

Scientific Reports volume 14, Article number: 21931 (2024) Cite this article

742 Accesses

2 Altmetric

Metrics details

In this study, the concept of rheo-optics was applied to explore the flow birefringence induced by the stress components along the camera’s optical axis since it is often overlooked in the traditional theories of photoelastic flow measurement. A novel aspect of this research is that it involved conducting polarization measurements on simple shear flows, specifically from a perspective in which a shear-velocity gradient exists along the camera’s optical axis. A parallel-plate-type rheometer and a polarization camera are employed for these systematic measurements. The experimental findings for dilute aqueous cellulose nanocrystal suspensions demonstrate that the flow birefringence can be expressed as a power law based on the power of the second invariant of the deformation-rate tensor. This suggests that flow birefringence can be universally characterized by the coordinate-independent invariants and a pre-factor determined by the direction of polarization measurement. By adjusting the nonlinear term in the stress-optic law, its applicability could be expanded to include three-dimensional fluid stress fields in which the stress is distributed along the camera’s optical axis.

The non-invasive measurement of the internal stress states of fluids in flow is of great importance for flow engineering and biomechanics. For bulk fluid pressure measurements, the method of measuring velocity fields with particle tracking velocimetry1 or particle image velocimetry2 is commonly used. However, they require spatial derivatives of velocity fields to determine the stress, in which the noise amplification problem appears. Other methods have also been proposed to evaluate the information corresponding to the stress based on unique phenomena, such as light scattering and magnetic resonance3,4,5. These methods can convey information relating to the stress driving the fluid intact. Among these approaches, measurement of stress or velocity gradient distributions by the photoelastic method6 has been applied and is attractive due to its high measurement sensitivity and experimental simplicity.

The photoelastic method was developed as a stress-measurement technique for solids, and it has been further developed and widely studied for 50 years7. Among solid materials, various studies considering the measurement of stresses and residual stresses, especially in glass, have been undertaken6,8. Stress-loaded materials change their refractive index with regard to the direction of polarization vibration in response to strain. Therefore, when two orthogonally polarized light beams are transmitted through a strained material, retardation will appear between them. The incident and emitted polarized light is represented by a composite vector of the two linearly polarized light beams, and their trajectories will differ due to the retardation of the composite vector. Depending on the trajectory difference, the result can be classified as linearly, elliptically, or circularly polarized light. When circularly polarized light is incident onto a stress-loaded material, it will be emitted as elliptically polarized light with retardation \(\Delta\) and orientation angle \(\phi\). The values of \(\Delta\) and \(\phi\) (photoelastic parameters) correspond to the principal-stress difference and the principal-stress direction, respectively6. In the photoelastic method, stress can be estimated from \(\Delta\) based on the stress-optic law (SOL)8,9:

where: \({{{\text {d}}}}h\) is the infinitesimal thickness of the material; \(C_1\) is the stress-optic coefficient; and \(\sigma _1\) and \(\sigma _2\) are the maximum and minimum principal stresses, respectively. However, Eq. (1) holds only for two-dimensional (2D) stress fields with no stress along the optical axis, or, if there is stress along the optical axis, then it must be uniform. This means that \(\Delta\) can be obtained by only considering the secondary principal stress difference10,11. Here, the secondary principal stress difference (\(\sigma _1 - \sigma _2\)) is the principal stress difference projected onto a plane perpendicular to the camera’s optical axis. In the case of 3D stress fields, applying the SOL is more complicated than in the above equation. For stress distributions along the camera’s optical axis, it is necessary to introduce the concepts of the “optically equivalent model”12 and “integrated photoelasticity”8,13. This makes it possible to consider that the 3D stress fields consist of sufficiently thin linear polarizers that can be assumed to be 2D stress fields. Thus, the polarization state transmitting through the 3D stress field can be calculated by multiplying the Mueller matrices of each optically equivalent model13.

McAfee and Pih observed 3D flow in a channel and described the dependence of the “isoclinic” (optical-anisotropy in this study) on shear strain rates perpendicular and parallel (camera’s optical axis) to the light beam3. However, they did not provide a direct relationship between the SOL and the measurements. Some attempts have been made in the past to extend the SOL by including stress components along the camera’s optical axis as a second-order term14,15. However, this has not been experimentally validated and the impact of the stress components along the optical axis continues to be neglected in photoelastic methods. Despite its recognised importance, it has not been incorporated into many studies.

In recent years, research has been conducted on the application of photoelastic methods to fluids, demonstrating that the principle of Eq. (1) can be applied to quasi-2D flows16,17. Despite this, several studies using quasi-2D flow channels have highlighted discrepancies between the birefringence predicted by SOL in the 2D case and the experimental data, particularly in the centre of the channel18,19,20. These reports can be considered experimental evidence that as the flow field becomes more 3D (e.g., circular tubes21 or rectangular channels with an aspect ratio close to 120,22), the stress distribution along the optical axis becomes more critical. These experimental findings have also been corroborated by numerical calculations, which have established that regions, where the impact of 3D effects is dominant, are close to the “plane of symmetry”23.

By introducing the concept of “rheo-optics”, studies of optical-anisotropy responses to the fluid’s deformation have been conducted. This is the simultaneous measurement of rheological and optical-anisotropy and has contributed to the understanding of the internal structures of complex fluids. Starting from the apparatus proposed by Lodge24 and Philipoff25, the filament stretching rheometer26 and the capillary breakup extensional rheometry dripping-onto-substrate (CaBER-DoS) method27,28 were developed to perform birefringence measurements on liquid polymers under uniaxial extension flow. To establish the relationship between the shear rate (shear stress) and optical-anisotropy, which is the aim of the present study, the flow given to fluids should be simple and well-controlled. Consequently, the use of rheometers (concentric cylinder (CC)-type17,29,30,31, parallel plate (PP)-type32,33,34,35, and cone plate (CP)-type36) has been reported, as they are easy to incorporate into optical-anisotropy measurement systems. The optical-anisotropy induced by secondary principal stress difference (principal stress difference acting in a plane perpendicular to the optical axis) can be measured from the shear-vorticity direction, e.g., in uniaxial extension systems and CC-type rheometers. Conversely, measurements from the shear direction provide a direct visualization of the stress along the optical axis rather than showing the principal stress difference. This is due to the existence of a velocity gradient along the optical axis, which is of interest to the present study.

This study aims to understand the 3D effect of birefringence by revisiting the stress-optic law in its original form (including the stress components along the optical axis) and verifying it quantitatively. Our study is novel and important because it is the first quantitative verification of this phenomenon. To simplify the discussion, fluid with properties similar to a Newtonian fluid was used as a first step in the study. This is to minimise the contribution of birefringence due to elastic or normal stress during shear imposition. As a validation of the stress-optic law, the tendency of birefringence induced by the stress along the optical axis was investigated experimentally as well as quantitatively and compared with the results obtained in previous studies. Note that the optical axis in our experiments is along the velocity gradient direction perpendicular to the parallel plates in a PP-type rheometer.

This section outlines the basic theory of polarization measurements, mainly focusing on flow birefringence and the SOL9.

Birefringent fluids are composed of crystals or polymer chains with large aspect ratios. This is the key to changing the refractive indices \(n_{\perp }\) and \(n_{\parallel }\) perpendicular and parallel to the direction of vibration of the transmitted light, respectively. One of the optical-anisotropies, birefringence \(\delta _n\), is a physical quantity that indicates the magnitude of the anisotropy of the refractive index. Birefringence is defined as the absolute value of the difference between the major and minor diameters of the index ellipsoid:

As long as there is no stress loading, birefringent fluids show optically isotropic properties (\(n_{\perp } = n_{\parallel }\)) because the particles are randomly oriented by Brownian motion. However, when shear is applied, the crystals or polymers become aligned in the direction corresponding to the stress, resulting in optical-anisotropy (\(n_{\perp } \ne n_{\parallel }\)). When the particles are strongly oriented in a particular direction, the anisotropy of the refractive index becomes stronger and the value of the birefringence increases (see Fig. 1a). Conversely, when the applied shear is reduced, the orientation becomes random again and birefringence is no longer induced (see Fig. 1b). This phenomenon is known as flow birefringence37; as examples of birefringent fluids, aqueous cellulose nanocrystal (CNC) suspensions38, wormlike micelle solutions18, and xanthan gum solutions39 are well known for showing birefringence.

Schematic of nanocrystals in solution (a) when no stress is applied, (b) under stress loading. Although needle-like nanocrystals are shown, the mechanism of onset is the same with macromolecular chains.

The retardation \(\Delta\) obtained by the photoelastic method is the summation (integrated value) of the birefringence along the optical axis. The retardation \(\Delta\) caused by flow birefringence \(\delta _n\) is related to the strain rate \(\dot{e}_{ij}\) inside the fluid14:

Here, the optical axis is defined as the z axis of a Cartesian coordinate system, and \(\alpha _1\) and \(\alpha _2\) are functions of the physical properties of the fluid. For Newtonian fluids, the stress is proportional to the strain rate. Therefore, Eqs. (3) and (4) can be expressed using stress20:

In these equations, \(C_1 =\alpha _1/\eta\) and \(C_2 = \alpha _2/\eta ^2\), in which \(\eta\) is the shear viscosity of the fluid. Here, \(C_1\) represents the sensitivity of optical-anisotropy to secondary principal stress difference, whereas \(C_2\) represents the sensitivity to stress along the optical axis. Aben and Puro15 also discussed the optical relationship based on Eqs. (5) and (6) and assumed that the stress components along the optical axis, i.e., \(\sigma _{xz}\), \(\sigma _{zy}\), and \(\sigma _{yz}\), were negligible. In other words, they made the assumption that \(C_2 = 0\), which leads to the proposal of:

It should be emphasized again that Eq. (7) is a relation that holds only for 2D stress fields. Moreover, Eq. (7) is an often-used expression in the solid-state photoelastic method15,40.

In this section, the details of the experiments are presented. Unless otherwise indicated, all experiments were repeated at least three times per case at \(25^\circ\)C.

Figure 2a shows a schematic of the rheo-optical measurement system. The stress-controlled rheometer (MCR 302, Anton Paar Co., Ltd.) was equipped with a parallel plate (PP43/GL-HT, Anton Paar Co., Ltd., plate radius \(R_o = 21.5\) mm) made of quartz glass and a flat glass plate (PTD200/GL, Anton Paar Co., Ltd.). In this system, shear flow is induced by clockwise rotation of the plate, and the average shear stress and torque are logged. Left-handed circularly polarized light is generated by attaching a polarizer and 1/4-wave plate to an LED light source (SOLIS-525C, Thorlabs Co., Ltd., wavelength \(\lambda = 525\) nm). This polarized light is reflected from the top of the plate by a mirror and emitted as elliptically polarized light with retardation \(\Delta\) and orientation angle \(\phi\). The transmitted elliptically polarized light is reflected by the mirror again and enters the polarization camera (CRYSTA PI-1P, Photron Co., Ltd.), which is equipped with a 524-nm band-pass filter. The polarization camera has a spatial resolution of up to \(512 \times 512\) pixels and a temporal resolution of 1.55 Mfps. In these experiments, all measurements were made with a resolution of \(512\times 512\) pixels (44.5 \(\upmu\)m/pixel) at 1000 fps. Additionally, the gap height between the plates was always fixed at \(H = 100\ \pm 5\) \(\upmu\)m, and the shear rate was set to 1000–10,000 \(\hbox {s}^{-1}\). Note that the shear rate is specified at 2/3 of the plate radius. This is the range in which it is expected that no artefacts will appear in the rheological measurement results (for further details, please see Appendix 1).

(a) Schematic of the rheo-optical setup, in which \(R_i\) and \(R_o\) are the inner and outer radii of the transparent part of the plate, respectively. (b) Representative intensity image for CNC suspension 1.0 wt% at shear rate \(\dot{\gamma }= 1000\) \(\hbox {s}^{-1}\). As shown in inset 1 of panel (b), retardation is obtained from four neighbouring polarizers. (c) Temporal evolution of retardation at no flow.

A polarization camera (CRYSTA PI-1P, Photron Co., Ltd.) was used to detect the retardation \(\Delta\), which is the integrated value of birefringence along the optical axis of the light transmitted through the apparatus. Using the phase-shifting method41, this was obtained from the radiance through linear polarizers oriented in four different directions \((0^\circ\), \(45^\circ\), \(90^\circ\), and \(135^\circ\)) in an area of \(2 \times 2\) pixels (as shown in inset 1 of Fig. 2b). Defining the light intensities detected at each of these pixels as \(I_1\), \(I_2\), \(I_3\), and \(I_4\), respectively, the retardation can then be given by:

where \(\lambda\) [m] is the wavelength of the light source. As the retardation is obtained from four linear polarizers, the spatial resolution was \(512 \times 512\) pixels, which is 1/4 of the \(1024 \times 1024\)-pixel light-intensity image. As there is no stress distribution along the optical axis, the retardation \(\Delta\) can be calculated as the product of the birefringence \(\delta _n\) and the gap height H. Therefore, in this study, we calculated the birefringence by dividing the retardation measurements by H.

An arc-shaped region of interest (ROI-A) located at a distance \(L\approx 3\) mm from the edge of the plate (outlined in yellow in Fig. 2b) was chosen, and the background-subtracted retardation field during flow was provided by the CRYSTA Stress Viewer software package (Photron Co., Ltd.). This background-subtraction process is based on the presence of a certain degree of unevenness on the plate surface or birefringence in the plate or stage itself, which can cause a spatial distribution of retardation even when there is no flow. The position of ROI-A was chosen to avoid possible interference at the edges of the plate by the liquid–air interface or the liquid meniscus42. ROI-A was also chosen to keep the uniformity of the lighting system, as the polarized image sensors use the light intensity itself to calculate the birefringence, as explained in the previous section. Another arc-shaped region, ROI-B (shown in red in inset 2 of Fig. 2b), with a width of 10 pixels (440 \(\upmu\)m), was also chosen for quantitative characterization of the flow birefringence. This was selected to ensure that the shear-rate variation along the radial direction was kept to a minimum, with \(\approx\)0.8% difference. Therefore, within ROI-B, the shear rate applied to the microscopic fluid is \(\dot{\gamma }_{{\text {rep}}} = 0.75R_0{\Omega }/H\), and this was considered as the representative shear rate in the polarization measurements. Here, \(\Omega\) [rad/s] is the angular velocity of the rotating plate.

The parallelism of the plate to the stage was not perfect. Therefore, due to changes in the integral of the birefringence, the measured retardation also changed with time. To exclude this effect, measurements were taken after a sufficient time had passed from the start of the plate’s rotation and averaged over 1000 frames. As the polarization cameras measure the light intensity distribution as an image, minute noise in the intensity values is included in the calculation results as a constant retardation value, even if the fluid is not flowing. Figure 2c shows the time variation of the spatial average retardation in ROI-B measured with no flow. The retardation at a certain level was measured to be 1.46 nm; therefore, the birefringence was calculated on the assumption that this value (\(\delta _{n, {\text {error}}}\approx 1.5\times 10^{-5}\)) was included as a measurement error or uncertainty.

Suspensions of CNCs (Alberta Pacific Co., Ltd.) of two different concentrations (weight percentage, wt%) were studied: 0.5 wt% and 1.0 wt%. These suspensions were prepared by mixing CNCs with ultrapure water using a hot stirrer at \(40^\circ\)C and 650 rpm for more than 24 h. The shear viscosities of these prepared CNC suspensions were then measured using a rheometer equipped with a 50-mm-diameter cone plate (CP50-0.5, Anton Paar Co., Ltd.) to obtain \(\eta\) values. The measurement results are shown in Fig. 3. Both CNC suspensions showed a significant shear-thinning nature at low shear rates (\(\dot{\gamma }< 10~{\text {s}}^{-1}\)) whereas they exhibited relatively weak shear-thinning at high shear rates (\(10^2~{\text {s}}^{-1}<{\dot{\gamma }}\)). Shear thinning of CNC suspensions is commonly observed as typical of lyotropic liquid crystals. The shear-thinning nature exhibited in the low shear rate region is due to the alignment of the chiral nematic liquid crystal domains. In contrast, at high shear rates, it is due to the disruption of the liquid crystal domains and the orientation of the individual CNC rods along the shear flow direction43. Here, the crowding number N is introduced to note the particle-particle interaction of the CNCs. N is expressed as \(N=2c_v(l/d)^2/3\) [–] using the particle volume concentration \(c_v\) (can be estimated from the CNC’s density), particle length l and particle width d44. For \(N<1\), particles are relatively free to move; conversely, collisions between particles occur for \(N>1\). The aspect ratio of the CNCs used in the present study is comparable to those used in previous studies (\(l/d\approx 16\))17,31,43. Referring to morphology and the CNC density of 1500 kg/m338, the N for 0.5 wt% and 1.0 wt% suspension was 0.56 and 1.13, respectively. \(N>1\) may be one of the reasons for the non-Newtonian nature of the 1.0wt% CNC suspension as shown in Fig. 3. The normal stress (Weissenberg effect) was also measured using the rheometer and confirmed to be small enough to be regarded as a measurement error.

Steady shear viscosity \(\eta\) for CNC suspensions at different concentrations. The inset shows an enlarged region at \(\dot{\gamma }= 10^2\)–\(10^4\) \(\hbox {s}^{-1}\).

Visualized birefringence fields at different shear rates are shown in Fig. 4a,b. As can be seen, the birefringence increased significantly as the shear rate was increased. This is because the CNCs attain more uniform orientations with increasingly strong shear leading to a stronger optical-anisotropy. The magnitude of birefringence obtained was \(\delta _n \sim O(10^{-5})\), comparable to the values found in previous studies29,31,45,46. In the experiments, we observed an uneven distribution of birefringence outside ROI-A, which was seen to vary unsteadily. This was considered to be an effect of interference at the liquid–air interface due to the plate’s rotation.

We now consider an analytical model to help understand the measurement results. A simple analytical model of shear flow between the plates can be expressed:

where H is the gap height. The strain tensor \(\dot{e}\) is:

This means that the birefringence measured in the present study was induced by strain (stress) components along the optical axis. Note that the birefringence measured by the present experimental system cannot be explained as far as the SOL used previously (which does not include the stress component along the optical axis, Eq. (1) or (7)) is applied. Substituting each component of the above tensor into Eqs. (3) and (4) yields:

Thus, the birefringence is particularly dependent on the square of the radius r and the square of the plate’s angular velocity \(\Omega\), as long as \(C_2\) and \(\eta\) are constant.

Visualized birefringence fields within ROI-A for (a) 0.5 wt% and (b) 1.0 wt% CNC suspensions under steady-state conditions. Radial birefringence distributions of the plate in (c) 0.5 wt% and (d) 1.0 wt% CNC suspensions, in which the grey shaded areas correspond to ROI-A. The error bars show the standard deviations and all data were time and circumferential averaged.

Line profiles of birefringence distributed across the radial direction of ROI-A at varying shear rates are shown in Fig. 4c,d. These show that the shear rate increases outwardly from the centre of the plate, which leads to an increase in the birefringence. When these results are compared with the analytical model, the quadratic dependence of the birefringence on r is evident in the experimental results; however, the birefringence is clearly not proportional to the square of \(\Omega\). This deviation indicates that \(C_2\) is not a constant value, and in particular that it varies with the shear rate. Details are given later in this section.

To further investigate the variation of the stress-optic coefficient \(C_2\), its magnitude was estimated and compared. From Eq. (11), the magnitudes of each parameter are \(\delta _n \sim O(10^{-5})\), \(\eta \sim O(1)\) mPa\(\cdot\)s, \(r\sim O(10)\) mm, \({\Omega } \sim O(10^1\)-\(10^2)\) rad/s, and \(H\sim O(100)~\upmu\)m. By magnitude comparison, the magnitude of \(C_2\) was estimated to be \(C_2\sim O(10^{-7}\text{--}10^{-6})~\rm Pa^{-2}\). To the best of our knowledge, this is the first systematic measurement report to identify the \(C_2\) value of a birefringent fluid. The stress-optic coefficient \(C_1\), generally known as the photoelastic modulus, has been widely investigated in solid polymers, and it is known to have a magnitude of \(O(10^{-12})~\rm Pa^{-1}\)47. In contrast, there have been a few reports on investigations of the \(C_1\) values of fluids. For reference, the \(C_1\) values obtained for fluids are: worm-like micelles \(O(10^{-7})~\rm Pa^{-1}\)18,45,48, aqueous xanthan gum solutions \(3.3\times 10^{-8}~\rm Pa^{-1}\)39, and a 0.5 wt% CNC suspension \(O(10^{-5})~\rm Pa^{-1}\)20. Although simple comparisons are difficult because of the different units, i.e., \(\hbox {Pa}^{-1}\) and \(\hbox {Pa}^{-2}\), the magnitude of \(C_2\) was found to be smaller than that of \(C_1\). Nevertheless, it is inappropriate to assume that \(C_2\) can be ignored in the presence of stress distribution along the optical axis as discussed in the Introduction of this paper.

Next, the trend of birefringence to the shear rate was investigated. For nanofiber suspensions, it has been reported that the no-slip condition in the velocity profile near the rotating plate is no longer applied49. Although the present experiment used a rigid-rod suspension, it could be worthwhile to consider the wall slip. To correct the shear rate applied to the suspension, wall slip was quantified by following the analysis described in the literature50, as only one of the correction methods. Hereafter, the discussion will be based on this corrected shear rate. The detail of the present analysis is described in Appendix 2. In Fig. 5, the vertical axis shows the spatiotemporally averaged birefringence \(\delta _{n,{\text {ave}}}\) in ROI-B, while the horizontal axis shows the representative shear rate \(\dot{\gamma }_{{\text {rep}}}\). It can be seen that the birefringence is increasing exponentially. When modelling the relationship between the flow birefringence and shear rate, the following empirical nonlinear model was proposed by Lane et al.31:

where w [–] is the mass fraction of CNC defined as \(w = ({\text {weight of CNC}})/({\text {total weight of suspension}})\), and A [s], n [–], and m [–] are fitting parameters. The experimental results were fitted using Eq. (12), and the black dash-dotted lines in Fig. 5 show the results. The fitting parameters were \(A=4.80\times 10^{-7}\) s, \(n=0.574\) and \(m = 1.28\). It should be emphasized that this model is based on the results of polarization measurements conducted from the vorticity direction to the shear using a CC-type rheometer, which is different from that used in the present study.

As described in the previous section, flow birefringence is induced by the aligned orientation of crystals or polymer chains dispersed in a fluid. For generalization, and to provide a better prospect for the physical interpretation of the present experimental results, we attempt to disentangle the birefringence \(\delta _n\) from the suspension concentration c and the direction of the polarization measurement from the effect of CNC alignment due to the flow. We used tensor invariants that are independent of the coordinate system. Therefore, the discussion henceforth will be based on the second invariant of the deformation-rate tensor, \(\Pi ~[\rm s^{-2}]\), i.e.,

Here, \(B~\hbox {[s}^{2}\)], \(\alpha\) [–], and \(\beta\) [–] are defined as new parameters. In the present system, \(\Pi\) can be calculated as:

where \(\textbf{E}\) is the deformation-rate tensor:

Since \(\Pi\) is proportional to the square of the shear rate \(\dot{\gamma }\) from Eq. (14), \(n = 2\alpha\) holds here. The \(\delta _n\) measurements were fitted using Eq. (13), and the results are shown by the black dashed lines in Fig. 6. The fitting parameters were \(B = 2.10\times 10^{-13}~\rm s^2\), \(\alpha = 0.287\) (\(n = 0.574\)), and \(\beta = 1.28\). Also represented in Fig. 5, it has been well reported that birefringence generally increases nonlinearly as the shear rate is increased29,31,51,52. This trend remains unchanged when organized by invariant, and the increase in birefringence per invariant decreases with decreasing gradient at higher invariants.

We focus our discussion on the exponent which characterizes the tendency of birefringence with the deformation of the fluid (solvent). Calabrese et al.38 studied the birefringence of 0.1 wt% CNC suspensions, and they proposed a proportionality of \(\delta _n\propto {\dot{\gamma }}^{0.9}\) from the results of their experiments. They reported CNC geometries with an average length of 260 nm and an average width of 4.8 nm. Note that the 0.1 wt% CNC suspension was stated to be a dilute region with no particle interactions. Lane et al. also found the relationship \(\delta _n\propto {\dot{\gamma }}^{0.537}\) in part of an investigation into whether CNC suspensions (0.7–1.3 wt%) could be used in studies of flow birefringence31. These suspensions were at concentrations above the dilute region, where particle interactions need to be considered. The exponent of 0.9 differs significantly from the present results, while 0.574 is relatively close. We assume that the difference in n results from different CNC particle-interaction behaviours. In the existing literature, the relationship between CNC particle shear alignment and rotational diffusion is described by the Péclet number \(Pe=|{{\text {E}}}|/D_r\), where \(|{{\text {E}}}|\) is the characteristic deformation rate53. When \(Pe\ge 1\), convective forces are strong enough to align CNCs to the flow direction, eventually inducing birefringence. Here, \(D_r\) is the rotational diffusion coefficient which plays an important role in flow birefringence. The \(D_r\) is determined by the CNC rod length l and the suspension concentration c, giving \(D_r\propto c^{-2}l^{-9}\)53,54. This indicates that the optical properties of CNC suspensions are due to different rod lengths and concentrations, resulting in different exponents. The present experimental results and those of Lane et al.31 indicate that regardless of the direction of polarization measurement with respect to shear, there seems to be a common physical background that leads to the flow birefringence following a power law.

Spatiotemporally averaged birefringence \(\delta _{n,{\text {ave}}}\) and corresponding fit using Eq. (12) with \(A=4.80\times 10^{-7}\), \(n=0.574\) and \(m = 1.28\), in which the error bars represent the standard deviation.

Spatiotemporally averaged birefringence \(\delta _{n,{\text {ave}}}\) and corresponding fit using Eq. (13) with \(A=2.10\times 10^{-13}\), \(\alpha =0.287\) and \(\beta = 1.28\), in which the error bars represent the standard deviation.

The consistency between the SOL (Eq. (11)) and the empirical relation (Eq. (13)) is now discussed. Since both equations are equally related by \(\delta _n\),

It can be seen from Fig. 3, the change the shear viscosity \(\eta\) was only O(0.1) mPa\(\cdot\)s. for changes in shear rate \(\dot{\gamma }\) between \(10^3\)-\(10^4~[\rm s^{-1}]\). To simplify the discussions, if we assume that the CNC suspension in the present study behaves as a Newtonian fluid, i.e., the change in the viscosity coefficient \(\eta\) with the shear rate is sufficiently small,

From Eq. (17), the coordinate-independent invariant \({\Pi }\) and the pre-factor B, which reflects the direction of polarization measurement, along with the mass fraction w, disentangle the magnitude of the stress-optic coefficient \(C_2\), which depends on the coordinate. In other words, birefringence (optical-anisotropy) can be universally described by an invariant and certain kinds of biases. The specific values of \(C_2\) in the present study are given in Fig. 7 using the infinite-shear viscosity at \(\dot{\gamma }= 10^4~\rm s^{-1}\) with 1.3 and 1.7 mPa\(\cdot\)s for 0.5 wt% and 1.0 wt% CNC suspensions, respectively.

Determined values of \(C_2\) for each CNC suspension in base-10 logarithmic scale calculated from Eq. (17).

These values are consistent with the range of the earlier magnitude-approximation result. Interestingly, as can be seen from Eq. (16) and Fig. 7, \(C_2\) decreased with the invariant. This agrees with the result in Fig. 4c,d, where it was experimentally shown that \(C_2\) is not a constant value. As described previously, \(C_2\) represents the sensitivity of stresses along the optical axis. In the present experiment, the sensitivity depended on the strain rate. Previous studies have reported that \(C_1\) is determined by the internal structure inside the respective fluid, suggesting that the value changes with structural changes. Ito et al. conducted rheo-optical experiments on worm-like micelles and reported that \(C_1\) value changed with shear-induced structure (SIS)45. Muto et al. also performed experiments on worm-like micelles and suggested that \(C_1\) depended only on the total amount of worm-like micelles28. Sridhar et al. reported that the \(C_1\) value of polymer decreased with increasing extension strain55. Therefore, the decrease in \(C_2\) values may reflect the internal structure of the CNC suspensions. Additionally, \(C_2\) was found to increase with increasing mass fraction of CNC suspension. This means that the sensitivity to stress (degree of optical-anisotropy) increased with increasing concentration, which is a well-known trend.

In summary, as in the PP-type rheometer used in the present experiments, \(C_2\) can be calibrated using a flow field with a shear-velocity distribution along the optical axis. The polarization measurement direction-dependent pre-factor B is then incorporated into the SOL as \(C_2\). We suggest that \(C_2\) as obtained in the present method can be replaced by \(C_2=f(B,~\Pi ,~ {w})\) in the SOL to give a good expression for the optical-anisotropy due to the stress distribution along the optical axis. A detailed study of the physical quantities dominating the pre-factor B, which determines \(C_2\), is a subject for future work. We believe that a systematic investigation of B (\(C_2\)) also has the potential to contribute to our understanding of the flow dynamics of complex fluids.

This paper provides the results of an investigation into the effect of the stress distribution along the camera’s optical axis on optical-anisotropy. The original form of the stress-optics law including the stress components along the camera’s optical axis, which has not been quantitatively verified, was revisited. In the experiments, rheo-optical measurements were performed using a PP-type rheometer and a polarization camera. For the birefringent fluid, a dilute aqueous CNC suspension was used, showing properties close to those of a Newtonian fluid. The birefringence was found to be induced by the shear-stress distribution along the camera’s optical axis, rather than by a principal stress difference as considered in conventional photoelastic theory in a 2D stress field. Note that in the present system, the optical axis is along the direction of the velocity gradient perpendicular to the parallel plate of the PP-type rheometer. The birefringence was found to increase nonlinearly and monotonically with the shear rate, following a power function of the second invariant of the deformation-rate tensor. In addition, from the well-known SOL and an empirical equation proposed in a previous study, the contribution of stress to optical-anisotropy can be determined by a pre-factor representing the direction of polarization measurement and the invariant. The degree of optical-anisotropy (birefringence) in the present study corresponds to the coefficient (stress-optic coefficient \(C_2\)) of the nonlinear term of the SOL. Moreover, it was found that the \(C_2\) is not constant to the second invariant (strain rate). This significant finding demonstrates the need for the SOL that accounts for these components. Based on the present experimental results, it can be said that the assumption that \(C_2=0\) is not appropriate, especially when the SOL is applied to 3D fluid stress fields in which the stress is distributed along the camera’s optical axis. By defining \(C_2\) as a non-zero invariant function, the contribution of the stress component along the camera’s optical axis in each experimental system can be quantitatively described. In the future, SOL with \(C_2\) term may provide a deeper understanding of complex fluids showing non-Newtonian properties such as particle-particle interactions and polymer extension and contraction.

Data sets generated during the current study are available from the corresponding author Y.T. upon reasonable request.

Adamczyk, A. A. & Rimai, L. 2-dimensional particle tracking velocimetry (PTV): Technique and image processing algorithms. Exp. Fluids6, 373–380. https://doi.org/10.1007/BF00196482 (1988).

Article CAS Google Scholar

van Oudheusden, B. W. PIV-based pressure measurement. Meas. Sci. Technol. 24, 032001. https://doi.org/10.1088/0957-0233/24/3/032001 (2013).

Article ADS CAS Google Scholar

McAfee, W. J. & Pih, H. Scattered-light flow-optic relations adaptable to three-dimensional flow birefringence. Exp. Mech. 14, 385–391. https://doi.org/10.1007/BF02324941 (1974).

Article Google Scholar

Martins, A. F., Esnault, P. & Volino, F. Measurement of the viscoelastic coefficients of main-chain nematic polymers by an NMR technique. Phys. Rev. Lett. 57, 1745–1748. https://doi.org/10.1103/PhysRevLett.57.1745 (1986).

Article ADS CAS PubMed Google Scholar

Odagiri, K. et al. Non-invasive evaluation of pulmonary arterial blood flow and wall shear stress in pulmonary arterial hypertension with 3D phase contrast magnetic resonance imaging. Springerplus 5, 1071. https://doi.org/10.1186/s40064-016-2755-7 (2016).

Article PubMed PubMed Central Google Scholar

Aben, H. & Guillemet, C. Photoelasticity of Glass (Springer Science & Business Media, 2012).

Lautre, N. K., Sharma, A. K., Kumar, P. & Das, S. A photoelasticity approach for characterization of defects in microwave drilling of soda lime glass. J. Mater. Process. Technol. 225, 151–161. https://doi.org/10.1016/j.jmatprotec.2015.05.026 (2015).

Article CAS Google Scholar

Ramesh, K. Developments in Photoelasticity: A Renaissance (IOP Publishing, 2021).

Prabhakaran, R. On the stress-optic law for orthotropic-model materials in biaxial-stress fields. Exp. Mech. 15, 29–34. https://doi.org/10.1007/BF02318522 (1975).

Article Google Scholar

Doyle, J. F. & Danyluk, H. T. Integrated photoelasticity for axisymmetric problems. Exp. Mech. 18, 215–220. https://doi.org/10.1007/BF02328416 (1978).

Article Google Scholar

Sampson, R. C. A stress-optic law for photoelastic analysis of orthotropic composites. Exp. Mech. 10, 210–215. https://doi.org/10.1007/BF02324034 (1970).

Article Google Scholar

Srinath, L. S. & Sarma, A. V. S. S. S. R. Determination of the optically equivalent model in three-dimensional photoelasticity. Exp. Mech. 14, 118–122. https://doi.org/10.1007/BF02324775 (1974).

Article Google Scholar

Yokoyama, Y. et al. Integrated photoelasticity in a soft material: Phase retardation, azimuthal angle, and stress-optic coefficient. Opt. Lasers Eng. 161, 107335. https://doi.org/10.1016/j.optlaseng.2022.107335 (2023).

Article Google Scholar

Doyle, J. F. On a nonlinearity in flow birefringence. Exp. Mech. 22, 37–38. https://doi.org/10.1007/BF02325702 (1982).

Article Google Scholar

Aben, H. & Puro, A. Photoelastic tomography for three-dimensional flow birefringence studies. Inverse Prob. 13, 215–221. https://doi.org/10.1088/0266-5611/13/2/002 (1997).

Article ADS MathSciNet Google Scholar

Noto, D., Tasaka, Y., Hitomi, J. & Murai, Y. Applicability evaluation of the stress-optic law in Newtonian fluids toward stress field measurements. Phys. Rev. Res. 2, 043111. https://doi.org/10.1103/PhysRevResearch.2.043111 (2020).

Article CAS Google Scholar

Lane, C., Baumann, F., Rode, D. & Rösgen, T. Two-dimensional strain rate imaging study using a polarization camera and birefringent aqueous cellulose nanocrystal suspensions. Exp. Fluids 65, 8. https://doi.org/10.1007/s00348-023-03730-8 (2023).

Article Google Scholar

Ober, T. J., Soulages, J. & McKinley, G. H. Spatially resolved quantitative rheo-optics of complex fluids in a microfluidic device. J. Rheol. 55, 1127–1159. https://doi.org/10.1122/1.3606593 (2011).

Article ADS CAS Google Scholar

Kim, S. G., Ok, C. M. & Lee, H. S. Steady-state extensional viscosity of a linear polymer solution using a differential pressure extensional rheometer on a chip. J. Rheol. 62, 1261–1270. https://doi.org/10.1122/1.5033499 (2018).

Article ADS CAS Google Scholar

Nakamine, K., Yokoyama, Y., Worby, W. K. A., Muto, M. & Tagawa, Y. Flow birefringence of cellulose nanocrystal suspensions in three-dimensional flow fields: Revisiting the stress-optic law. Cellulose 31,7405–7420. https://doi.org/10.1007/s10570-024-06045-x (2024).

Article CAS Google Scholar

Alizadehgiashi, M. et al. Shear-induced alignment of anisotropic nanoparticles in a single-droplet oscillatory microfluidic platform. ACS Publ. 34, 322–330. https://doi.org/10.1021/acs.langmuir.7b03648 (2018).

Article CAS Google Scholar

Kim, J. et al. Monitoring the orientation of rare-earth-doped nanorods for flow shear tomography. Nat. Nanotechnol. 12, 914–919. https://doi.org/10.1038/nnano.2017.111 (2017).

Article ADS CAS PubMed Google Scholar

Clemeur, N., Rutgers, R. P. G. & Debbaut, B. Numerical evaluation of three dimensional effects in planar flow birefringence. J. Nonnewton. Fluid Mech. 123, 105–120. https://doi.org/10.1016/j.jnnfm.2004.07.002 (2004).

Article CAS Google Scholar

Lodge, A. S. A network theory of flow birefringence and stress in concentrated polymer solutions. Trans. Faraday Soc. 52, 120–130. https://doi.org/10.1039/TF9565200120 (1956).

Article CAS Google Scholar

Philippoff, W. Stress-optical analysis of fluids. Rheol. Acta 1, 371–375. https://doi.org/10.1007/BF01989069 (1961).

Article Google Scholar

Rothstein, J. P. & McKinley, G. H. A comparison of the stress and birefringence growth of dilute, semi-dilute and concentrated polymer solutions in uniaxial extensional flows. J. Nonnewton. Fluid Mech. 108, 275–290. https://doi.org/10.1016/S0377-0257(02)00134-9 (2002).

Article CAS Google Scholar

Muto, M. & Tagawa, Y. Unsteady rheo-optical measurements of uniaxially extending liquid polymers. arXiv preprint (2022). https://doi.org/10.48550/arXiv.2204.13450.

Muto, M., Yoshino, T. & Tamano, S. Simultaneous measurement of extensional stress and flow birefringence field for uniaxially extending worm-like micellar solutions, (2024) https://doi.org/10.48550/arXiv.2404.17643.

Decruppe, J. P., Cressely, R., Makhloufi, R. & Cappelaere, E. Flow birefringence experiments showing a shear-banding structure in a CTAB solution. Colloid Polym. Sci. 273, 346–351. https://doi.org/10.1007/BF00652348 (1995).

Article CAS Google Scholar

Tanaka, R., Li, S., Kashiwagi, Y. & Inoue, T. A self-build apparatus for oscillatory flow birefringence measurements in a co-cylindrical geometry. Nihon Reoroji Gakkaishi 46, 221–226. https://doi.org/10.1678/rheology.46.221 (2018).

Article CAS Google Scholar

Lane, C., Rode, D. & Rösgen, T. Birefringent properties of aqueous cellulose nanocrystal suspensions. Cellulose 29, 6093–6107. https://doi.org/10.1007/s10570-022-04646-y (2022).

Article CAS Google Scholar

Hausmann, M. K. et al. Dynamics of cellulose nanocrystal alignment during 3D printing. ACS Nano 12, 6926–6937. https://doi.org/10.1021/acsnano.8b02366 (2018).

Article CAS PubMed Google Scholar

Kádár, R., Fazilati, M. & Nypelö, T. Unexpected microphase transitions in flow towards nematic order of cellulose nanocrystals. Cellulose 27, 2003–2014. https://doi.org/10.1007/s10570-019-02888-x (2020).

Article CAS Google Scholar

Detert, M., Santos, T. P., Shen, A. Q. & Calabrese, V. Alignment-rheology relationship of biosourced rod-like colloids and polymers under flow. Biomacromol 24, 3304–3312. https://doi.org/10.1021/acs.biomac.3c00347 (2023).

Article CAS Google Scholar

Sato, T. et al. Two-dimensional rheo-optical measurement system to study dynamics and structure of complex fluids. Appl. Rheol. 34, 20240006 (2024). https://doi.org/10.1515/arh-2024-0006.

Oba, N. & Inoue, T. An apparatus for birefringence and extinction angle distributions measurements in cone and plate geometry by polarization imaging method. Rheol. Acta 55, 699–708. https://doi.org/10.1007/s00397-016-0952-5 (2016).

Article CAS Google Scholar

Maxwell, J. C. On double refraction in a viscous fluid in motion. Proc. R. Soc. Lond. 22, 46–47. https://doi.org/10.1098/rspl.1873.0011 (1874).

Article Google Scholar

Calabrese, V., Haward, S. J. & Shen, A. Q. Effects of shearing and extensional flows on the alignment of colloidal rods. Macromolecules 54, 4176–4185. https://doi.org/10.1021/acs.macromol.0c02155 (2021).

Article ADS CAS Google Scholar

Yevlampieva, N. P., Pavlov, G. M. & Rjumtsev, E. I. Flow birefringence of xanthan and other polysaccharide solutions. Int. J. Biol. Macromol. 26, 295–301. https://doi.org/10.1016/S0141-8130(99)00096-3 (1999).

Article CAS PubMed Google Scholar

Riera, J. D. & Mark, R. The optical-rotation effect in photoelastic shell analysis. Exp. Mech. 9, 9–16. https://doi.org/10.1007/BF02327872 (1969).

Article Google Scholar

Onuma, T. & Otani, Y. A development of two-dimensional birefringence distribution measurement system with a sampling rate of 1.3MHz. Optics Commun. 315, 69–73. https://doi.org/10.1016/j.optcom.2013.10.086 (2014).

Article ADS CAS Google Scholar

Ewoldt, R. H., Johnston, M. T. & Caretta, L. M. Experimental challenges of shear rheology: How to avoid bad data. In Spagnolie, S. E. (ed.) Complex Fluids in Biological Systems: Experiment, Theory, and Computation, Biological and Medical Physics, Biomedical Engineering, 207–241, https://doi.org/10.1007/978-1-4939-2065-5-6 (Springer, New York, NY, 2015).

Shafiei-Sabet, S., Hamad, W. Y. & Hatzikiriakos, S. G. Rheology of Nanocrystalline Cellulose Aqueous Suspensions. Langmuir 28, 17124–17133. https://doi.org/10.1021/la303380v (2012).

Article CAS PubMed Google Scholar

Kerekes, R. & Schell, C. Characterization of Fibre Floccula tion Regimes by a Crowding Factor. J. Pulp Pap. Sci. 18, 32–38 (1992).

CAS Google Scholar

Ito, M., Yoshitake, Y. & Takahashi, T. Shear-induced structure change in shear-banding of a wormlike micellar solution in concentric cylinder flow. J. Rheol. 60, 1019–1029. https://doi.org/10.1122/1.4961034 (2016).

Article ADS CAS Google Scholar

Lane, C., Rode, D. & Rösgen, T. Two-dimensional birefringence measurement technique using a polarization camera. Appl. Opt. 60, 8435–8444. https://doi.org/10.1364/AO.433066 (2021).

Article ADS PubMed Google Scholar

Koike, Y. & Tagaya, A. Photonics polymer (Kyoritsu Publishing, 2004).

Shikata, T., Dahman, S. J. & Pearson, D. S. Rheo-optical behavior of wormlike micelles. Langmuir 10, 3470–3476. https://doi.org/10.1021/la00022a019 (1994).

Article CAS Google Scholar

Lauri, J., Koponen, A., Haavisto, S., Czajkowski, J. & Fabritius, T. Analysis of rheology and wall depletion of microfibrillated cellulose suspension using optical coherence tomography. Cellulose 24, 4715–4728. https://doi.org/10.1007/s10570-017-1493-5 (2017).

Article CAS Google Scholar

Yoshimura, A. & Prud’homme, R. K. Wall Slip Corrections for Couette and Parallel Disk Viscometers. J. Rheol. 32, 53–67. https://doi.org/10.1122/1.549963 (1988).

Article ADS CAS Google Scholar

Santos, T. P., Calabrese, V., Boehm, M. W., Baier, S. K. & Shen, A. Q. Flow-induced alignment of protein nanofibril dispersions. J. Colloid Interface Sci. 638, 487–497. https://doi.org/10.1016/j.jcis.2023.01.105 (2023).

Article ADS CAS PubMed Google Scholar

Peebles, F. N., Prados, J. W. & Honeycutt, E. H. Jr. Birefringent and rheologic properties of milling yellow suspensions. J. Polymer Sci. Part C Polym. Symp. 5, 37–53. https://doi.org/10.1002/polc.5070050105 (1964).

Article Google Scholar

Doi, M., Edwards, S. F. & Edwards, S. F. The theory of polymer dynamics (oxford university press, 1988).

Maguire, J. F., McTague, J. P. & Rondelez, F. Rotational diffusion of sterically interacting rodlike macromolecules. Phys. Rev. Lett. 45, 1891–1894. https://doi.org/10.1103/PhysRevLett.45.1891 (1980).

Article ADS CAS Google Scholar

Sridhar, T., Nguyen, D. A. & Fuller, G. G. Birefringence and stress growth in uniaxial extension of polymer solutions. J. Nonnewton. Fluid Mech. 90, 299–315. https://doi.org/10.1016/S0377-0257(99)00080-4 (2000).

Article CAS Google Scholar

Download references

This work was funded by JSPS KAKENHI Grant Nos. JP20H00222 and JP20H00223, Presto Grant No. JPMJPR21O5, and Japan Agency for Medical Research and Development (Grant No. JP22he0422016). We thank Prof. T. Rösgen, Dr. C. Lane, and Prof. T. Shikata for their valuable discussions and suggestions. We also thank Dr. J. Yee for providing language help and writing assistance.

Department of Mechanical Systems Engineering, Tokyo University of Agriculture and Technology, Koganei, Tokyo, 184-8588, Japan

William Kai Alexander Worby, Kento Nakamine & Yoshiyuki Tagawa

Micro/Bio/Nanofluidics Unit, Okinawa Institute of Science and Technology Graduate University, Kunigami-gun, Okinawa, 904-0495, Japan

Yuto Yokoyama

Department of Electrical and Mechanical Engineering, Nagoya Institute of Technology, Nagoya, Aichi, 466-8555, Japan

Masakazu Muto

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

Y.T. and M.M. made substantial contributions to the design of the work. W.W. conducted the experiments, contributed to the data analysis and interpretation, and drafted the original manuscript. K.N. and Y.Y. made substantial contributions to the interpretation of the data and to the revision of the manuscript drafts. All authors agree to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

Correspondence to Yoshiyuki Tagawa.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

Worby, W.K.A., Nakamine, K., Yokoyama, Y. et al. Examination of flow birefringence induced by the shear components along the optical axis using a parallel-plate-type rheometer. Sci Rep 14, 21931 (2024). https://doi.org/10.1038/s41598-024-73087-3

Download citation

Received: 07 May 2024

Accepted: 13 September 2024

Published: 20 September 2024

DOI: https://doi.org/10.1038/s41598-024-73087-3

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative