Sensors: Compact Viscometer Prototype For Remote in Situ Analysis of Sludge
Sensors: Compact Viscometer Prototype For Remote in Situ Analysis of Sludge
Article
Compact Viscometer Prototype for Remote In Situ
Analysis of Sludge
Tomas Fried 1,2, * , David Cheneler 1 , Stephen D. Monk 1 , C. James Taylor 1 and
Jonathan M. Dodds 2
1 Engineering Department, Lancaster University, Lancaster LA1 4YW, UK
2 National Nuclear Laboratory, Workington CA14 3YQ, UK
* Correspondence: [email protected]
Received: 6 June 2019; Accepted: 24 July 2019; Published: 26 July 2019
Abstract: On the Sellafield site there are several legacy storage tanks and silos containing sludge of
uncertain properties. While there are efforts to determine the chemical and radiological properties of
the sludge, to clean out and decommission these vessels, the physical properties need to be ascertained
as well. Shear behaviour, density and temperature are the key parameters to be understood before
decommissioning activities commence. However, limited access, the congested nature of the tanks
and presence of radioactive, hazardous substances severely limit sampling and usage of sophisticated
characterisation devices within these tanks and therefore, these properties remain uncertain. This
paper describes the development of a cheap, compact, and robust device to analyse the rheological
properties of sludge, without the need to extract materials from the site in order to be analysed.
Analysis of a sludge test material has been performed to create a suitable benchmark material for the
rheological measurements with the prototype. Development of the device is being undertaken with
commercial off the shelf (COTS) components and modern rapid prototyping techniques. Using these
techniques, an initial prototype for measuring shear parameters of sludge has been developed, using
a micro-controller for remote control and data gathering. The device is also compact enough to fit
through a 75 mm opening, maximising deployment capabilities.
1. Introduction
The landscape of the nuclear industry is changing nuclear power plants, and as reprocessing
facilities approach the end of their planned lifespan, decommissioning has become a new focal point
for the nuclear industry. This is reflected in academic research as well, with research moving from
reprocessing operations to decommissioning [1]. An example of a decommissioning facility in the
United Kingdom is the Sellafield site. Once the home of the first commercial nuclear reactor producing
electricity on an industrial scale, the nuclear fuel reprocessing capability is approaching the end of its
life cycle with the Thermal Oxide Reprocessing Plant (Thorp) facility stopping shearing in 2018 and
Magnox reprocessing will end in 2020 [2,3].
Decommissioning of nuclear facilities requires a post-operational cleanout of active and hazardous
materials. Using the Sellafield site as an example, there are a number of legacy tanks, silos and other
containers with nuclear waste, often in the form of suspension, that will need to be cleaned out as
part of the decommissioning plans. However, to accomplish this, analysis of remnant material must
be performed to determine and minimise the risks, waste packages and ascertain the appropriate
processes. Due to the hazardous nature of nuclear materials, this often has to be done remotely. Remote,
in situ characterisation is, therefore, one of the key aspects of nuclear decommissioning [4].
The current strategy for analysis of remnant materials mainly focuses on visual inspection,
radiological measurements and chemical composition analysis. Radiological measurements vary
from using scintillators in storage tanks to determine the activity of remnant materials, pipe crawling
robots assessing the contamination of pipework or autonomous platforms for analysis of floor
contamination [5–8]. Current academic research focuses on imaging systems and the deployment of
these in challenging environments. Advances in chemical composition analysis in nuclear facilities
have introduced laser-induced breakdown spectroscopy (LIBS) as one potential method of remote
analysis [9,10]. Visual inspection aids with assessing structures, waste locations and remnant material
quantities in locations inaccessible to human workers. Most recently, advances in 3D scanning
technologies (such as Light Detection and Ranging, or LiDAR) have been the focal point of research
for visualizing remote areas. Other areas of research interest focus on the deployment platforms for
both traditional camera systems and LiDAR in remote areas, whether through limited openings or
underwater [11].
Currently used methods to ascertain rheological parameters of materials are predominantly
rotational and oscillation rheometry using various geometries, either with samples or on-line [12–14].
Squeeze flow and capillary rheometry is frequently used in industries other than nuclear and academic
research has been exploring these areas, too [15–17]. However, all of these methods utilise large,
expensive equipment and development focus on increasing the precision and measurement range
of devices.
In situ rheology in hazardous environments is, however, an unexplored research area. Research
has been limited to test materials, or it relies on collecting samples from storage containers [18–20].
Novel methods for rheological measurements that have been developed are focused on miniaturisation
and small sample methods, such as quartz crystal microbalance, but these are still in early research
stages [21].
This paper presents work undertaken to develop a compact, cheap and robust device for shear
behaviour analysis and its initial validation in a non-active laboratory environment. The aims are to
propose a novel prototype device; to manufacture this device using rapid prototyping methods and
commercial off the shelf (COTS) components; to calibrate the new device using industrial standard
materials and to validate its performance at a laboratory bench scale. Section 2 provides a summary of
materials and conventional instruments that have been used for analysis and manufacturing. Section 3
provides a description of the designed device, the tests performed, and the results of these tests. Finally,
the discussion and conclusions are presented in Section 4.
The second material was zirconium-molybdate (ZM) suspended in 0.5 M nitric acid. The sample
was synthesised by Johnson Matthey and supplied as a suspension. Some of the supernate was
decanted, and a sample was dried to ascertain the volume fraction of the suspension. ZM is an example
of a higher viscosity test material. It is the most novel simulant being used in the nuclear industry [22].
A Malvern Mastersizer 3000 (Malvern Panalytical Ltd, Malvern, UK) particle characterisation
system with a wet dispersion unit was used to ascertain the particle size distribution (PSD) of the
two test materials, as shown in Figure 1 Both materials were dispersed in de-ionized water (refractive
index 1.33). TiO2 was stirred at 1800 rpm and measured at approximately 3.7% laser obscuration,
using a refractive index 2.493 and an absorption index 1. ZM was stirred at 2400 rpm and measured at
approximately 30.7% laser obscuration, using a refractive index 1.19 and an absorption index 0.
Figure 1. Particle size distribution (PSD) of the titanium dioxide and zirconium molybdate samples
used in this research.
A majority of the structural components and other parts were manufactured using a fused
deposition modelling 3D printer (Ultimaker 3 extended, Ultimaker, Utrecht, The Netherlands).
Two different materials were utilised as follows. Linkages and outer casings not coming into contact
with chemical substances were made using Ultimaker PLA. PLA is a biodegradable, cheap and
straightforward to use material offering high stiffness. It is the most appropriate material for rapidly
prototyping initial development prototypes, as required for the present work. Measurement geometries
and their corresponding outer cups were manufactured with XSTRANDTM GF30-PP—glass fibre
reinforced polypropylene (Owens Corning, Toledo, OH, USA). This material offers much higher
chemical resistance and thermal stability compared to PLA, making it suited for use with chemically
aggressive substances, such as ZM in nitric acid. Table 1 compares the mechanical and other properties
of these materials.
Table 1. Comparison of the parameters of Ultimaker PLA and XSTRANDTM GF30-PP 3D printing
filaments [23,24].
Figure 3. Cutaway render of all the geometries used with the prototype.
Two of the geometries and their respective cups used with the device are designed based on the
dimensions proposed by the ISO 3219 standard, widely used in the design of conventional benchtop
rheometers [25]. They are manufactured with the same dimensions as the V25 vane and C25 bob used
on the Bohlin CVO100 rheometer. The device, however, diverges from the ISO standard in two main
regards. One, since it is deployed by inserting it in a sample, the cup is open on the bottom. This allows
the material to enter the cup and surround the measurement geometry as the device is inserted into the
material. Secondly, the immersion height is not controlled by sampling but by deploying the device.
In the experiments described below, the material height is controlled in the test material vessel.
Sensors 2019, 19, 3299 5 of 14
In potential deployment scenarios, the material height may not be sufficient to fully immerse long
geometries. For this reason, three other geometries have been designed and used, diverging from
the ISO standard. Another consideration is preventing the slip of the geometry when it rotates in the
material. Vanes are often used with suspensions and higher viscosity samples, as they slip less than
bobs when shearing these samples. Finally, inserting the geometry itself affects the structure of the
sample. Vanes penetrate samples much easier than bobs. Based on all of these considerations, all three
non-standard geometries (i.e., Vane 1, 2 and 3 in Figure 3) take the form of four-bladed vanes.
As described in the ISO standard, the representative shear rate in the middle of the gap between
.
the geometry and the cup γrep is determined as shown in Equation (1).
2
. 1 + ( rre )
γrep = Ω· 2
i
, (1)
( rrei ) − 1
where Ω is the angular velocity of the geometry, re is the radius of the outer cylinder (cup) and ri is the
radius of the inner cylinder (geometry). It is apparent that the shear rate is dependent on the ratio of
the dimensions of the geometry and the cup, and thus, the ratio has been maintained constant with all
the geometries. To accommodate potential low material heights, the non-standard vanes have been
designed to be shorter and thus require less material to be fully immersed. To ensure that the torque is
not below the detection threshold when performing the measurement, the radius of the geometry has
been increased in comparison with the ISO geometries.
A summary of the geometry dimensions that have been manufactured and used with the prototype
is outlined in Table 2, while the geometry and cup configuration is illustrated in Figure 3.
The outer case and internal linkages were 3D printed using a fused deposition modelling
printer with Ultimaker PLA. The bob and vane geometries and the outer cups were 3D printed with
XSTRANDTM GF30-PP. Overall, the largest configuration (ISO vane or bob) of the device is 219 mm long
(i.e., edge of the cable gland to end of the geometry cup) and 75 mm in diameter at the widest point.
The prototype was designed to remove as many sensitive electronic components from the
deployment environment. Hence, the electronic circuitry consists of two main assemblies (Figure 4).
The Control Unit houses all the data acquisition and control electronics necessary to drive the
mechanism and read the data from an on-board sensor. All of the components used are COTS products.
A Cortex-M4F (180 Mhz rated core speed) based microcontroller was selected as it supports the C++
based programming language, has integrated digital-to-analogue conversion (DAC) functionality and
16-bit resolution analogue-to-digital conversion (ADC) functionality. It controls the stepper motor
through a Pololu DRV8834 stepper motor driver (Pololu Corporation, Las Vegas, NV, USA) and
reads the output of the Bourns 6630S0D-C28-A103 continuous turn potentiometer (Bourns, Riverside,
CA, USA) in the viscometer device. Furthermore, it reads two buttons used to start and stop the
pre-tensioning and measurement procedures and controls the 74HC595 shift register used to display
the current state of the prototype through 4 LED indicators on the top panel of the Control Unit. All
of the Control Unit components are in a PLA housing. The overall dimensions of the housing are
100 mm × 80 mm × 53 mm. All of the units can be seen in Figure 5.
Figure 5. Experimental setup with the prototype. Devices from left to right: laptop, external power
supply, control unit, viscometer prototype device in a beaker immersed in a sample of TiO2 .
The prototype is programmed with two procedures. One is to zero the mechanism to establish
the origin of the measurement. This procedure runs the geometry in the opposite direction of the
measurement rotation in a ramp-up of rotational velocity, which deflects the mechanism. When the
mechanism deflects beyond a pre-programmed value, the program switches to a ramp down to allow
the mechanism to settle, in what will be used as a zero position. This procedure serves to eliminate all
potential deflection and tensioning of the mechanism during deployment.
The second procedure is the measurement itself. The prototype is programmed to supply a
ramp-up of rotational velocities in integer values. The stepper motor supplies this velocity for a total
of 12 s, measuring the deflection of the mechanism in the last six seconds to allow the mechanism to
fully deflect before the measurement. In those six seconds, the microcontroller reads and saves the
value from the potentiometer every 400 ms, logs the time since it started the measurement procedure
and saves the rotational velocity value it is currently supplying to the geometry.
Figure 6. Shear stress as a function of the rotational velocity of silicone viscosity standard oil on a
commercial rheometer with: (a) C25 bob and cup configuration; (b) V25 vane and cup configuration.
The same calibration oil was then used with the prototype and all of its geometries. The prototype
supplied rotational velocities in the range of 1 to 100 rpm with the bob and 2 to 200 rpm with the vanes
(corresponding to roughly 1 to 100 s−1 on the commercial rheometer). The spacing is approximately
logarithmic, limited by the fact that the prototype can supply only integer values of revolutions per
minute. Five runs with the standard oil were performed for each geometry and the raw data were
fitted with linear functions, as illustrated by Figure 7.
Figure 7. Potentiometer bit output as a function of revolutions per minute of silicone viscosity standard
oil with the prototype with (a) C25 bob and cup configuration; (b) Selection of vanes.
The raw data collected with the potentiometer in the prototype indicates the torque applied to the
vane. The raw data were fitted with a linear function and, using data collected with the conventional
rheometer, the relationship between the raw data and expected shear stress based on the conventional
rheometer measurement was determined. Using this approach, the conversion constants between raw
data and shear stress were obtained for each geometry. Summary of calibration data can be seen in
Table 3.
Figure 8 summarises the calibration results for all five geometries. Vane 3 offers the best low shear
stress detection rate. However, the linearity is limited in the higher shear stress regions. Smaller vanes,
namely the ISO and Vane 1 geometries exhibit better repeatability and linearity in higher rotational
velocities, but as the torques are lower, low shear stress detection is relatively poor. The bob geometry
and Vane 2 exhibit the poorest results overall. With the bob geometry, this can be attributed to the
potentially highest sensitivity to axial misalignment of the bob and the cup, with the results being
affected by the uneven gap between the bob and the cup.
Figure 8. Shear stress results obtained with all five prototype geometries during the calibration
procedure, in comparison to a conventional rheometer.
All the geometries exhibit repeatable, systematic deviations from the expected linear response to
the silicone oil. Therefore, a compensation procedure was created to remove the nonlinearities. For all
the geometries, the relative deviation from the expected value of shear stress was plotted against
the shear stress measured with the prototype. Polynomial functions were fitted to the values of the
deviation as a function of the measured shear stress, as illustrated in Figure 9. To ensure a good fit
using low order functions, only data in the range necessary to be compensated were fitted. These
functions are subsequently used to determine the compensation value for the measured results. Actual
reported shear stress values are, therefore, defined as the measured shear stress minus the value of the
polynomial function at that measured shear stress.
It is important to note that the results obtained with this procedure must be used carefully and
that the compensation should only be used within the range of the fitted data, see Figure 9. The results
discussed below illustrate instances in which the procedure appears to be appropriate or not.
One aspect of using vanes in rotational rheology that needs to be considered when interpreting
the results are secondary flows and vortices appearing in the sample. Research on quantifying the
onset of secondary flows using vanes is limited and only applied to one phase materials. Onset of
Taylor vortices can be described using the Taylor number. For concentric cylinders, with the outer
cylinder stationary and the inner cylinder rotating, in the middle of the gap between the geometries,
the Taylor number Ta is often calculated as shown in Equation (2) [26].
Sensors 2019, 19, 3299 9 of 14
ri (re − ri )3 Ω2i
Ta = , (2)
ν2
where ν is the kinematic viscosity of the sample.
Figure 9. Non-linearity of the prototype with Vanes 1, 2 (first five values are off the scale and not used in
calculations) and Vane 3 using a silicone viscosity standard oil and fitting functions for compensation.
The onset of turbulent flow in the same setup as considered in the calculation of the Taylor number
is referred to as the rotational Reynolds number and can be calculated as illustrated in Equation (3).
ri (re − ri )Ωi
Re = . (3)
ν
It is important to note that these parameters have been validated and used with concentric
cylinders and their validity for vanes with particulate suspensions has not been investigated in the
literature to date. Rotational rheology using vanes assumes the material shears along the blades of the
vanes, and thus vanes are often described using the same principles as concentric cylinders. Current
research and rheology laboratory equipment suppliers suggest that vanes are more susceptible to
both turbulent flows and Taylor vortices. However, it is apparent from Equations (2) and (3) that it
is assumed here that both vanes and bob geometries exhibit the same behaviour. Both Taylor and
Reynolds numbers have been calculated using the viscosity and angular velocity determined on
the rheometer with the calibration silicone oil and the dimensions of the geometries used with the
prototype. The results are illustrated in Figure 10.
Sensors 2019, 19, 3299 10 of 14
Figure 10. Potential indicators of secondary effects when measuring standard silicone oil: (a) Taylor
number for all geometries; (b) Reynolds number for all geometries.
Figure 11. Shear stress results obtained with a rheometer, the prototype (using vane 3) and compensated
prototype values with titanium dioxide suspension.
The Taylor and Reynolds numbers calculated using the rheometer collected data, but with the
dimensions of Vane 3, are illustrated by Figure 12. Secondary flow effects appear to start at around
80 rpm, corresponding to Ta = 6 and Re = 7.5. These values are low, compared to reported thresholds
indicating the onset of secondary effects in current literature.
Sensors 2019, 19, 3299 11 of 14
Figure 12. Potential indicators of secondary effects when measuring titanium dioxide suspension with
Vane 3: (a) Taylor number; (b) rotational Reynolds number.
3.4. ZM Suspension
ZM suspension has a much higher viscosity in the low rotational velocity regions. It is also
demonstrably shear thinning. It is apparent from Figure 13 that the values obtained with Vane 2 and
Vane 3 are higher than expected. With Vane 2, this can be attributed to the relatively poor calibration
results. With Vane 3, some over-reporting is expected due to the non-linearity of the response of the
prototype. Vane 1 appears to offer the best approximation of the results without compensating the
results. All three geometries report shear thinning behaviour.
Figure 13. Shear stress results obtained with the rheometer and the prototype with Vanes 1, 2 and 3:
(a) Measured values; (b) Compensated values.
The compensation procedure does not seem appropriate for Vane 1 in this case, as it changes the
slope of the curve. Vane 2 is affected similarly. However, the values obtained with the prototype are
still higher than expected across the measurement range. The procedure works well for Vane 3, as the
slope of the data is still the same, but results are much closer to the expected values.
The potential onset of secondary flow seems to be delayed as opposed to the lower viscosity
sample. This is consistent with the calculations of the Taylor and Reynolds numbers, as increased
viscosity is expected to delay the onset of secondary flows. The secondary flow effects appear only with
Vane 3, from approximately 150 rpm, corresponding to Ta = 3.5 and Re = 6, as illustrated in Figure 14.
These values are similar to those obtained with titanium dioxide suspension.
Sensors 2019, 19, 3299 12 of 14
Figure 14. Potential indicators of secondary effects when measuring zirconium molybdate suspension
with Vane 1, 2 and 3: (a) Taylor number; (b) rotational Reynolds number.
Multiple underwater remotely operated vehicles (ROVs) have been designed and deployed on various
sites and are capable of carrying payloads such as the proposed device [27]. The prototype could be
carried and detached from the carrying platform when it reaches the required position. Lastly, for
remote tanks, pipe crawler and tracked robots have been used to carry sensors and other payloads and
these platforms could deploy the device in harder to access areas [28,29]. In all instances, the tether
would serve not only a data transferring function but also as a means of recovery of the device from
these areas. The device is light (0.5 kg in the laboratory scale configuration), and the casing is robust
enough to support a recovery procedure using a tether. The simple assembly would make it easy to
reuse the internal components if necessary and decontaminate the outer casing, cups and geometries
after use. If the device cannot be decontaminated, the low price supports the case to use the device
only once before disposal.
Potential future development of the device could include lowering friction in the mechanism, which
would improve low shear stress detection and repeatability of the measurements. A more complex
mathematical model for the non-linearity compensation would be beneficial for the interpretation of
the results. Lastly, understanding the onset of secondary flows in various particulate suspensions
would be necessary to increase the operating range of the instrument.
Overall, the precision and performance of the proposed mechanism makes it suitable for
preliminary analysis of sludge in situ. The capability of the current stage of the prototype would make
it useful for determining the rough properties of the materials, to inform further decisions as necessary
in the upcoming steps of Post Operational Clean Out (POCO).
Author Contributions: T.F. designed and developed the device. T.F. developed the methodology and interpreted
the results with the assistance of J.M.D., D.C. and S.D.M. T.F. wrote the paper with assistance of J.M.D., D.C.,
S.D.M. and C.J.T.
Funding: This research was jointly funded by the CINDe consortium and by the UK Engineering and Physical
Sciences Research Council (EPSRC), via grant EP/R02572X/1.
Acknowledgments: This work has been supported by the Centre for Innovative Nuclear Decommissioning
(CINDe), which is led by the National Nuclear Laboratory, in partnership with Sellafield Ltd. and a network of
Universities that includes the University of Manchester, Lancaster University, the University of Liverpool and the
University of Cumbria.
Conflicts of Interest: The authors declare no conflict of interest.
References
1. NIRAB. The UK Civil Nuclear R&D Landscape Survey (February 2017). Available online: http://www.nirab.
org.uk/media/10671/nirab-123-4.pdf (accessed on 19 March 2019).
2. GOV.UK End of Reprocessing at Thorp Signals New Era for Sellafield. Available online: https://www.gov.uk/
government/news/end-of-reprocessing-at-thorp-signals-new-era-for-sellafield (accessed on 19 March 2019).
3. Sellafield Ltd. Corporate Plan 2016/17. Available online: https://www.gov.uk/government/uploads/system/
uploads/attachment_data/file/627566/SEL11098_corporate-plan_web.pdf (accessed on 19 March 2019).
4. OECD. R&D and Innovation Needs for Decommissioning Nuclear Facilities; Radioactive Waste Management,
OECD Publishing: Paris, France, 2014; ISBN 9789264222199.
5. Keele, B.D.; Addleman, R.S.; Blewett, G.R.; Mcclellan, C.S.; Troyer, G.L. Non-destructive in situ measurement
of radiological distributions in Hanford Site waste tanks. IEEE Trans. Nucl. Sci. 1996, 43, 1821–1826.
[CrossRef]
6. Maley, S.; Jones, H.; Whittaker, W. Novel Radiometry for In-Pipe Robotic Inspection of Holdup Deposits in
Gaseous Diffusion Piping. In Proceedings of the WM2019 Conference, Phoenix, AZ, USA, 3–7 March 2019.
7. Bird, B.; Griffiths, A.; Martin, H.; Codres, E.; Jones, J.; Stancu, A.; Lennox, B.; Watson, S.; Poteau, X.
Radiological Monitoring of Nuclear Facilities: Using the Continuous Autonomous Radiation Monitoring
Assistance Robot. IEEE Robot. Autom. Mag. 2018. [CrossRef]
8. Tsitsimpelis, I.; Taylor, C.J.; Lennox, B.; Joyce, M.J. A review of ground-based robotic systems for the
characterization of nuclear environments. Prog. Nucl. Energy 2019, 111, 109–124. [CrossRef]
Sensors 2019, 19, 3299 14 of 14
9. Baudelet, M. 17 Applications of laser spectroscopy in nuclear research and industry. In Laser Spectroscopy for
Sensing; Woodhead Publishing: Sawston, UK, 2014; pp. 522–543.
10. Martin, M.Z.; Allman, S.; Brice, D.J.; Martin, R.C.; Andre, N.O. Exploring laser-induced breakdown
spectroscopy for nuclear materials analysis and in-situ applications. Spectrochim. Acta Part B At. Spectrosc.
2012, 74, 177–183. [CrossRef]
11. Griffiths, A.; Dikarev, A.; Green, P.R.; Lennox, B.; Poteau, X.; Watson, S. AVEXIS-Aqua vehicle explorer for
in-situ sensing. IEEE Robot. Autom. Lett. 2016, 1, 282–287. [CrossRef]
12. Mould, S.; Barbas, J.; Machado, A.V.; Nóbrega, J.M.; Covas, J.A. Measuring the rheological properties of
polymer melts with on-line rotational rheometry. Polym. Test. 2011, 30, 602–610. [CrossRef]
13. Klein, C.O.; Spiess, H.W.; Calin, A.; Balan, C.; Wilhelm, M. Separation of the nonlinear oscillatory response
into a superposition of linear, strain hardening, strain softening, and wall slip response. Macromolecules 2007,
40, 4250–4259. [CrossRef]
14. Poulos, A.; Renou, F.; Jacob, A.; Koumakis, N.; Petekidis, G. Large amplitude oscillatory shear (LAOS) in
model colloidal suspensions and glasses: Frequency dependence. Rheol. Acta 2015, 54, 715–724. [CrossRef]
15. Schermeyer, M.; Sigloch, H.; Bauer, K.C.; Oelschlaeger, C.; Hubbuch, J. Squeeze flow rheometry as a novel
tool for the characterization of highly concentrated protein solutions. Biotechnol. Bioeng. 2016, 113, 576–587.
[CrossRef] [PubMed]
16. Kádár, R.; Naue, I.F.C.; Wilhelm, M. First normal stress difference and in-situ spectral dynamics in a high
sensitivity extrusion die for capillary rheometry via the ‘hole effect’. Polymer 2016, 104, 193–203. [CrossRef]
17. Morita, A.T.; Toma, M.S.; Paoli, M.-A. De Low cost capillary rheometer, transfer molding and die-drawing
module. Polym. Test. 2006, 25, 197–202. [CrossRef]
18. Paul, N.; Biggs, S.; Edmondson, M.; Hunter, T.N.; Hammond, R.B. Characterising highly active nuclear waste
simulants. Chem. Eng. Res. Des. 2013, 91, 742–751. [CrossRef]
19. Pease, L.F.; Daniel, R.C.; Burns, C.A. Slurry rheology of Hanford sludge. Chem. Eng. Sci. 2019, 199, 628–634.
[CrossRef]
20. Paul, N. Characterisation of Highly Active Nuclear Waste Simulants; University of Leeds: Leeds, UK, 2014.
21. Botha, J.A.; Ding, W.; Hunter, T.N.; Biggs, S.; Mackay, G.A.; Cowley, R.; Woodbury, S.E.; Harbottle, D.
Quartz crystal microbalance as a device to measure the yield stress of colloidal suspensions. Colloids Surf. A
Physicochem. Eng. Asp. 2018, 546, 179–185. [CrossRef]
22. Paul, N.; Hammond, R.B.; Hunter, T.N.; Edmondson, M.; Maxwell, L.; Biggs, S. Synthesis of nuclear waste
simulants by reaction precipitation: Formation of caesium phosphomolybdate, zirconium molybdate and
morphology modification with citratomolybdate complex. Polyhedron 2015, 89, 129–141. [CrossRef]
23. PLA|Professional 3D Printing Made Accessible|Ultimaker. Available online: https://ultimaker.com/en/
resources/49911-pla#technical (accessed on 17 April 2019).
24. XSTRANDTM GF30-PP-Owens Corning Composites. Available online: https://www.owenscorning.com/
composites/product/xstrand-gf30-pp (accessed on 17 April 2019).
25. British Standards Institution. Plastics—Polymers/Resins in the Liquid State or as Emulsions or
Dispersions—Determination of Viscosity Using a Rotational Viscometer with Defined Shear Rate. EU Patent
EN ISO 3219:1995, 15 March 1994.
26. White, F.M. Fluid Mechanics, 7th ed.; McGraw-Hill: New York, NY, USA, 2011; ISBN 978-0-07-352934-9.
27. Nancekievill, M.; Jones, A.R.; Joyce, M.J.; Lennox, B.; Watson, S.; Katakura, J.; Okumura, K.; Kamada, S.;
Katoh, M.; Nishimura, K. Development of a Radiological Characterization Submersible ROV for Use at
Fukushima Daiichi. IEEE Trans. Nucl. Sci. 2018, 65, 2565–2572. [CrossRef]
28. Tibrea, S.; Nance, T.; Kriikk, E. Robotics in hazardous environments—Real deployments by the savannah
river national laboratory. J. South Carol. Acad. Sci. 2011, 9, 5.
29. Rood, M.; Kelly, D. Hanford Double Shell Tank Integrity Inspection–Integration of Customized Robotics.
In Proceedings of the Waste Management Symposia, Phoenix, AZ, USA, 18–22 March 2018.
© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access
article distributed under the terms and conditions of the Creative Commons Attribution
(CC BY) license (http://creativecommons.org/licenses/by/4.0/).