1. Introduction
Rheology is the study of the flow and deformation of liquids and soft solids and has many applications in the materials and processing industries. A large group of materials can be described as viscoelastic, that is, they exhibit behaviour that can be described by a combination of elastic and viscous stress when undergoing strain. To capture the stress/strain relationships for these materials, discrete models consisting of a combination of elastic springs (E) and viscous dampers (η) connected in different ways can be used. Typical examples would be Maxwell and Kelvin Voigt systems, shown in
Figure 1.
The mechanical impedance in the Laplace domain relating stress to strain in each case is given by:
To investigate material behaviour, dynamic material analysis is often employed where a time-varying stress or strain is applied to the material and the corresponding strain or stress is recorded from which the impedance is derived. Decomposing this into the relationship at individual frequencies, it is then possible to represent the material behaviour as a frequency response function (FRF). The FRF for both a typical Maxwell and a typical Kelvin Voigt system is shown in
Figure 2 in terms of magnitude and phase, and is otherwise known as a Bode plot. It is worth mentioning here that the phase in a bode plot is the frequency-equivalent of the time shift between the output signal when compared to the input signal as the signal passes through a system. This is not to be confused with the term phase used to describe the physical form of a material, for example, liquid or solid.
These FRFs can be understood by replacing the Laplace operator s in Equations (1) and (2) with the Fourier operator jω. Of note is the fact that the phase at low and high frequencies will tend to either zero or radians. At intermediate frequencies, the phase gradually shifts from the low-frequency level to the high-frequency level. This is standard behaviour for many viscoelastic systems and is characteristic of first-order systems.
Standard time-domain tests such as a creep test or a relaxation test include the application of a step input, which results in an exponential decay in the response from the system. However, there are materials that show a power-law time decay of the form
where α varies between 0 and 1, which implies a transfer function with fractional powers of
s corresponding to fractional derivative operators [
1,
2]. For these materials, the phase characteristic does not follow patterns such as those in
Figure 2 but instead shows a constant phase, intermediate between zero and
radians over a wide frequency range.
One of the first recorded observations of power-law decay was when monitoring the behaviour of cheese as it ripens during the maturation process [
3]. Since then, there have been many applications where fractional behaviour of materials has been studied [
4,
5,
6]. In these cases, modelling has involved the incorporation of a “springpot”, an element with an impedance described by
As the name suggests, this element has a characteristic somewhere between an elastic spring and a viscous dashpot where
k is a constant and
α describes the fractional order. A zero value for
α indicates pure elasticity while a value of one implies purely viscous behaviour. The springpot is typically added to a conventional model such as a standard linear solid (SLS) [
7,
8], replacing either a spring or a dashpot.
It has been shown mathematically that fractional order systems can be related to a fractal structure [
9,
10]. For materials, this is relevant because chemical bonds forming elastic connections in a viscous background can create self-similar structures within the material and, hence, give rise to fractional behaviour.
In this paper, a fractal structure is developed to represent the dynamic behaviour of phase changing materials. The fractal topology is then related to a technique [
11] used for modelling fractional systems. Finally, some examples of physical transformations are considered where the aforementioned technique is used to model material phase transformation. Liquid/solid transformations are considered, which raises the issue of gelation, which will be considered in
Section 5, when discussing results.
2. Fractal and Fractional Modelling
For simple elastic materials, a Hookean elastic model is often used where the stress is solely dependent on a constant, the elasticity, multiplied by the strain. For purely viscous materials, the stress is dependent on the viscosity multiplied by the strain rate. For both these models, the material is considered to be uniform and without structure. These models can be combined to produce simple structural models, termed viscoelastic models such as those in
Figure 1. These models are based on long length scales and result in integer order descriptions of the material dynamic response. If one includes more detailed structural topology when using springs and dashpots to model materials at small length scales, non-integer order behaviour can be included. One technique to do this is to use fractal patterns [
9,
10].
Figure 3 shows a simple fractal structure consisting of springs and dashpots and how it can be condensed to an equivalent impedance
X.
Looking at the structures in
Figure 3, one can write
so that
This gives rise to a fixed fractional order of 0.5. For phase changing materials, the fractional power changes as the transformation progresses, requiring a more sophisticated model. A recursive fractal model has been developed [
12] where, depending on the fractal structure, any fractional power value can be achieved. However, it was recognised that visualisation of the physical structure is obscured due to the recursive approach. Here, an iterative/recursive technique is used to establish fractal structures that are somewhat easier to understand and that can be related to the physical phenomenon seen in practice.
To begin this analysis, the impedance derived in Equation (5) will be termed
; the use of 4 as a subscript will become apparent later. Now, consider that all the elastic elements in
Figure 3 are replaced by an element of type
. Using the same approach to derive Equation (5), we can then find a new element
such that
Iterating once again by replacing the element
in the network described by Equation (6) by the element
, we can define a new fractal pattern
such that
On the other hand, one could take the fractal network in
Figure 3 and replace the viscous element with an element of type
to find a new fractal network of type
Now, iterating this procedure again by replacing the element
in the network described by Equation (9) by the element
, we can define a new fractal pattern
such that
This can be repeated to give
The process of phase change from liquid to solid could be continued to give ever smaller fractional increments, but if one assumes a start point with a purely viscous system
and a finish point with a purely elastic system
then, the intervening stages can be described by
to
.
If one regards the elastic elements as representing chemical bonds, the material transitions from a weakly connected system dominated by viscous behaviour to one that is highly connected having many chemical bonds and is dominated by an elastic response. From to , fractal patterns consolidate to form elastic elements, while, from to , further fractal components account for crosslinking between molecular chains.
The transformation process can be defined by taking the end state
and dividing it by the start state
to define a transition function
HT [
11] so that
Then, to determine the state of the system at any stage, the transition function raised to a fractional power
is multiplied by the start condition
Plotting the powers of
s as the transition occurs through stages 0 to 8, one obtains the points shown in
Figure 4.
In practice, it would not be expected that the material phase changes all occur in perfect synchrony. To account for this and still allow for a fractal pattern, it is more realistic to consider a random fractal pattern resulting from a distribution of mechanical impedances which shift as the transformation progresses. The continuous gradual change can then be represented by the red trend line in
Figure 4.
This trend is often seen in the cure characteristics of adhesives. An example of this can be seen in
Figure 5 where the strength of adhesive bonds is presented as the cure develops for a methacrylate adhesive at different bond gaps.
It seems reasonable to assume the strength would be related to the degree of phase change.
Figure 6 shows the fractional power used to model experimental results for a cyanoacrylate adhesive as it cures [
11].
Although the start and end values differ from the idealised case in
Figure 4, the shape is similar and can be described by a sigmoid function as follows:
where
is the fractional value at the start of the process,
is the final value,
is the stage of the process at the median value of
and
is a constant. The independent variable can be time or some other parameter such as temperature depending on the physical process involved.
The analysis so far results in a springpot of order
, which, as mentioned in the Introduction, is the element used in conjunction with springs and dashpots to model systems showing fractional behaviour. The fractal model presented in
Section 2 assumes a perfectly viscous material transitioning to a perfectly elastic model during the phase change process, resulting in a dashpot with varying fractional power. In practice, materials will often show various viscoelastic properties before and after phase change. To account for such behaviour, a more general approach may be used [
11], whereby the start and end state used to formulate the transition function in Equation (14) can be any appropriate integer order viscoelastic function such as a Maxwell or Kelvin Voigt element rather than just a purely viscous or purely elastic one. Thus, by inspecting the bode response of the material before and after transformation, the viscous and elastic elements can be substituted by the appropriate viscoelastic models to generate a transition function appropriate to the behaviour shown in the experimental data. This allows the modeller to create the overall description of the transformation process
based on some understanding of the physical characteristics of the material.
3. Test Methods and Materials
The dynamic mechanical analysis of liquids and soft solids is carried out using the Micro Fourier Rheometer (GBC Scientific). It is an oscillatory squeeze film rheometer that uses an axisymmetric geometry where flat circular plates of radius R move relative to one another along their axes (z-axis) with amplitudes up to 20 μm. A schematic of the test geometry is shown in
Figure 7.
The top plate moves with displacement zp, squeezing/stretching the material between it and the fixed lower plate. The displacement in the upper plate is measured to calculate the instantaneous height h, and the force induced in the lower plate is recorded from which the complex modulus of the sample under test is determined.
One of the issues when testing phase changing materials is that they will eventually bond the rheometer together, rendering it unusable. To preserve the integrity of the rheometer, a removable system using neodymium magnets is employed to hold the top and bottom plates in place. A representation of the top plate is shown in
Figure 8 and a similar setup is used for the bottom plate. The forces incurred during testing are well below the magnetic holding force on the removable plate, but once testing is finished, it is then possible to slide the bonded removable plates out from the rheometer and separate them.
The force and displacement data are Fourier-transformed and analysed using a non-linear spectral technique [
14] so that the material characteristics can be extracted in terms of complex modulus and phase. It is important to acquire the data quickly so that the transient nature of the material properties can be treated as being quasi-static. For this reason, a band-limited random signal between 0 and 20 Hz is used, which allowing for a 40 Hz sample rate, and 400 points requires 10 s to acquire. This is comparable to the high-speed techniques described in [
15] when using the swept sine technique. The spectral data are calculated using just one record rather than the usual ensemble average technique due to the transient nature of the material. Data at low frequencies are somewhat compromised by the piezoelectric force transducer, which is evident in the results presented.
6. Conclusions
A fractal topology has been used to describe the rheological phase transition in materials. This results in a sigmoidal-shaped evolution of material properties that is often seen in practise. The fractal structure can be related to a fractional derivative modelling approach that has been used previously when characterising adhesive cure and which is used here to model the whipping process of cream, the solidification of gelatine and the melting of EVA. The three materials were chosen to represent three different mechanisms of phase transition, rheopexy for whipping cream, chemical bonding for the solidification of gelatine and temperature-induced melting for EVA. Using the fractional derivative modelling approach, a good fit was obtained when modelling most of the material’s dynamic response during phase transition.
For solid-like behaviour, the modelling approach adopted here requires start and end points that behave as solids. Once a liquid-type model is incorporated at either or both of the terminal points, the overall system will not be able to sustain a yield condition and will, therefore, have liquid-like behaviour. When dealing with liquid/solid transitions, a non-zero magnitude of the frequency response function at 0 rad/s was used to indicate the onset of gelation. In the case of the melting process for EVA, it was not possible to accurately capture both the magnitude and the phase with the fractional model. This arose because the material exhibits a yield stress in terms of magnitude and, yet, does not have zero phase at 0 rad/s.