Modelling a wool scour bowl

Similar documents
Natural convection adjacent to a sidewall with three fins in a differentially heated cavity

Robust fuzzy control of an active magnetic bearing subject to voltage saturation

An investigation into the effects of heat transfer on the motion of a spherical bubble

Modelling the induced force of attraction in electrorheological nanofluids

The lens of freshwater in a tropical island: the two interface case

A new moving boundary shallow water wave equation numerical model

Science 14 Unit A: Investigating Properties of Matter Chapter 3 Mixtures and Their Uses pp WORKBOOK Name:

An experimental and simulated investigation of particle flow through a conveyor rock box

Semi-analytical solutions for cubic autocatalytic reaction-diffusion equations; the effect of a precursor chemical

Controlling Downstream Particle Segregation with a Twisted Standpipe

Lecture 25: Manufacture of Maleic Anhydride and DDT

Experimental Procedure. Lab 406

Problem Set 4 Issued: Wednesday, March 18, 2015 Due: Wednesday, April 8, 2015

Factors That Affect Solubility

Two dimensional particle solution of the extended Hamilton Jacobi equation

Analyzing Mass and Heat Transfer Equipment

AquaFlux A New Instrument for Water Vapour Flux Density Measurement

2. Conservation of Mass

Influence of forced material in roller compactor parameters I.

Chemistry 1B Experiment 17 89

Practical reliability approach to urban slope stability

Process Design Decisions and Project Economics Prof. Dr. V. S. Moholkar Department of Chemical Engineering Indian Institute of Technology, Guwahati

7 Mathematical Methods 7.6 Insulation (10 units)


B-1. Attachment B-1. Evaluation of AdH Model Simplifications in Conowingo Reservoir Sediment Transport Modeling

Principles of Convection

LAMMPS Simulation of a Microgravity Shear Cell 299r Progress Report Taiyo Wilson. Units/Parameters:

Recovering the absolute phase maps of two fringe patterns with selected frequencies

University of Bristol - Explore Bristol Research. Link to publication record in Explore Bristol Research PDF-document.

An investigation on the impacts of density currents on the sedimentation in dam reservoirs using TCM model; case study: Maroon dam

Darcy s Law, Richards Equation, and Green-Ampt Equation

Module 15 : Grit Chamber. Lecture 19 : Grit Chamber

Exercise sheet 5 (Pipe flow)

15. GRIT CHAMBER 15.1 Horizontal Velocity in Flow Though Grit Chamber

FOLLOW THE ENERGY! EARTH S DYNAMIC CLIMATE SYSTEM

GENERALISATION OF THE TWO-SCALE MOMENTUM THEORY FOR COUPLED WIND TURBINE/FARM OPTIMISATION

Plumes and jets with time-dependent sources in stratified and unstratified environments

Chromatography Extraction and purification of Chlorophyll CHM 220

Flow analysis in centrifugal compressor vaneless diffusers

MERGING OF SHEET PLUMES IN TURBULENT CONVECTION

The scientific side of cleaning-in-place

MAE 598 Project #1 Jeremiah Dwight

A polyhedral model for carbon nanotubes

Centrifugation. Tubular Bowl Centrifuge. Disc Bowl Centrifuge

Circular dispersive shock waves in colloidal media

Ohio s State Tests ITEM RELEASE SPRING 2018 GRADE 5 SCIENCE

Deformability Modulus of Jointed Rocks, Limitation of Empirical Methods and Introducing a New Analytical Approach

Closed duct flows are full of fluid, have no free surface within, and are driven by a pressure gradient along the duct axis.

Part I.

A model of an aircraft towing a cable-body system

Boundary element methods in the prediction of the acoustic damping of ship whipping vibrations

Production planning of fish processed product under uncertainty

Vertical Mantle Heat Exchangers for Solar Water Heaters

Recovery of hydrochloric acid from metal pickling solutions by membrane distillation

Chemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore

arxiv: v1 [physics.flu-dyn] 11 Oct 2012

Model completion and validation using inversion of grey box models

Comments on "Generating and Counting Binary Bent Sequences"

This section develops numerically and analytically the geometric optimisation of

A First Course on Kinetics and Reaction Engineering Example 23.1

A note on Devaney's definition of chaos

Experiment (4): Flow measurement

INDBOND 3000 Dry Strength Resin for Paper

Thermo-Hydro-Mechanical-Chemical (THMC) Modelling of the Bentonite Barrier in Final Disposal of High Level Nuclear Waste

CADE/CAODC DRILLING CONFERENCE October 20 & 22, 2003 Calgary, Alberta, Canada

Infiltration from irrigation channels into soil with impermeable inclusions

EXPERIMENT No.1 FLOW MEASUREMENT BY ORIFICEMETER

An-Najah National University Civil Engineering Department. Fluid Mechanics. Chapter 1. General Introduction

Comparing Color and Leader Line Approaches for Highlighting in Geovisualization

Numerical Study on the Condensation Length of Binary Zeotropic Mixtures

Investigation into Sand Deposition and Transportation in Multiphase Pipelines Phase 2

Flux Enhancement By Shear Free Surfaces In A Turbulent Convection

Is the Basis of the Stock Index Futures Markets Nonlinear?

For example an empty bucket weighs 2.0kg. After 7 seconds of collecting water the bucket weighs 8.0kg, then:

One student complains that this is impossible because the car would not be able to stop before hitting the cab of the lorry.

Unit 4, Lesson 04, Reaction Mechanisms O 2 (g) 7 CO 2 (g) + 8 H 2 O (g)

Lake Pollution Model

The distortion observed in the bottom channel of Figure 1 can be predicted from the full transport equation, C t + u C. y D C. z, (1) x D C.

Hydraulic Fundamentals Hydraulics Definition Advent of Oil Hydraulics Pascal s Law Pressure

Some long-period random number generators using shifts and xors

On integer matrices obeying certain matrix equations

Axial profiles of heat transfer coefficients in a liquid film evaporator

Numerical Investigation of Conjugate Natural Convection Heat Transfer from Discrete Heat Sources in Rectangular Enclosure

The Solution Process. Section 2. Increasing the Surface Area of the Solute

Qualitative Analysis I - Cations

1. Introduction. Technical presentation of the Bernoulli Filter

IGCSE(A*-G) Edexcel - Chemistry

Chapter 6: Solid-Liquid Separation in WWTPs. Raúl Muñoz Pedro García Encina

Abstract. Introduction

An Interruption in the Highway: New Approach to Modeling Car Traffic

Adsorption (Ch 12) - mass transfer to an interface

Module 15 : Grit Chamber. Lecture 19 : Grit Chamber

On the INFOMAX algorithm for blind signal separation

Best Practices 8.0. Understanding and treating an industrial static problem

AN INVESTIGATION INTO THE EFFECTS OF HEAT TRANSFER ON THE MOTION OF A SPHERICAL BUBBLE

Answers. Mixtures. Year 7 Science Chapter 4. p75. p77. p79

ORGANIZATION OF THE SCHEME FOR QUALITATIVE ANALYSIS

THE MOST IMPORTANT BIT

On small defining sets for some SBIBD(4t-1, 2t-1, t-1)

Diffraction of ocean waves around a hollow cylindrical shell structure

Transcription:

University of Wollongong Research Online Faculty of Informatics - Papers (Archive) Faculty of Engineering and Information Sciences 26 Modelling a wool scour bowl James Caunce University of New South Wales Geoff Mercer University of New South Wales S. I. Barry University of New South Wales Timothy R. Marchant University of Wollongong, tim@uow.edu.au Publication Details Caunce, J. F., Mercer, G., Barry, S. & Marchant, T. R. (26). Modelling a wool scour bowl. ANZIAM Journal, 47 (EMAC25), C9-C33. Research Online is the open access institutional repository for the University of Wollongong. For further information contact the UOW Library: research-pubs@uow.edu.au

Modelling a wool scour bowl Abstract Wool scouring is the process of washing dirty wool after shearing. Our model simulates, using the advectiondiffusion equation, the movement of contaminants within a scour bowl. The effects of varying the important parameters are investigated. Interesting, but simple, relationships are found which give insight into the dynamics of a scour bowl. Disciplines Physical Sciences and Mathematics Publication Details Caunce, J. F., Mercer, G., Barry, S. & Marchant, T. R. (26). Modelling a wool scour bowl. ANZIAM Journal, 47 (EMAC25), C9-C33. This journal article is available at Research Online: http://ro.uow.edu.au/infopapers/2636

ANZIAM J. 47 (EMAC25) pp.c9 C33, 26 C9 Modelling a wool scour bowl J. F. Caunce S. I. Barry G. N. Mercer T. R. Marchant (Received 2 October 25; revised 3 June 26) Abstract Wool scouring is the process of washing dirty wool after shearing. Our model simulates, using the advection-diffusion equation, the movement of contaminants within a scour bowl. The effects of varying the important parameters are investigated. Interesting, but simple, relationships are found which give insight into the dynamics of a scour bowl. Contents Introduction C2 School of Physical, Environmental and Mathematical Sciences, University of New South Wales at adfa, Canberra, Australia. mailto:s.barry@adfa.edu.au School of Mathematics and Applied Statistics, University of Wollongong, Wollongong, Australia. See http://anziamj.austms.org.au/v47emac25/caunce3 for this article, c Austral. Mathematical Soc. 26. Published June 22, 26. ISSN 446-8735

Contents C2 2 Discretisation model C24 3 Results C27 4 Conclusion C32 References C32 Introduction Wool scouring is the process of washing shorn wool to remove unwanted contaminants. These contaminants fall into three categories: suint, grease and dirt. Suint is water soluble and removed in the first bowl. Grease is removed with the aid of detergents. Insoluble dirt (dust, stones and vegetable matter) is removed from the wool through agitation. The conventional wool scouring machine usually consists of two to four scour bowls and three rinse bowls. Figure shows a schematic of a scouring machine with three scour and three rinse bowls. Wool enters at the feed hopper at the left of the figure, travels through the scour and rinse bowls and leaves through the dryer. Water is transported from bowls 6 5 for efficiency. However, within each bowl, water flows in the same direction as the wool, hence the flow is not a classic separation process. For example, in Figure the water in bowl 3 proceeds from left to right within the bowl before being pumped to the left of bowl 2. Scour bowls are filled with a liquor primarily consisting of water, but also soaps, detergents and the contaminants that have been removed from the wool. The wool is driven as a saturated mat along the top of the bowl, by sets of harrows. Both the wool and water enter from the same side of the bowl. The agitation frees the dirt which settles down the bowl to a set of three to four settling tanks as shown in Figure 2. Periodically each tank

Introduction C2 Wool Flow Fresh Water Flowback Dryer Feed Hopper Scour Bowls Rinse Bowls HST Centrifuge Wool Grease Small Solids and Grease Heavy Solids Non Grease Wastes Waste Water Figure : The conventional scouring system Wool region Cross-flow area Settling tanks x y Squeeze rollers Figure 2: The conventional scour bowl. With dots showing the discretisation points of the model and arrows showing the liquor flow.

Introduction C22 is drained to purge the collecting contaminants. The liquor is cleaned and recycled to the bowls. The purges occur according to a drainage pattern, that specifies a purge cycle for each tank, shown in Figure 3, where the time of purges are indicated by the peaks. Grease is swollen by the detergents in the bowls, but most of the grease is removed at the squeeze rollers by the rush of water through the wool caused by the rollers. The grease rich liquor at the top of the bowls and the squeeze rollers is sent to a centrifuge where the grease is separated to be sold as lanolin. Previous work has modelled the scouring system to find the average conditions of the whole system at steady state. These models treat each scour bowl or even the whole scouring machine as a single entity [3, 4]. The limitations of these models are that: they cannot find the state of the system at any time before it reaches steady state or explore the effect of the periodic drainage; and it is not possible to see how the contaminants move within each scour bowl. A study group report, [], briefly considered a time dependent model of a scour bowl compartmentalised into eight regions, but a more complete model is needed. The aim of our model is to evaluate the flow of contaminants accurately enough to reflect the important features of the system, but be efficient enough to allow a computer simulation to run quickly. The scour bowl is not an exact system, many of the flows and parameters associated with the scour bowl are very variable, with some of the flows turbulent. Thus, fully three dimensional fully turbulent flow models, for example, would provide few useful results compared to simple models whilst taking far longer to compute. We believe that the simpler model developed in this paper allows a better understanding how contamination in wool is effected by operating parameters. For example, Figures 5 7 in the results section, show physical output of contamination in the wool versus the key operating parameters of wool and settling velocities. Our model looks at only one scour bowl, rather than the entire scour

Introduction C23 tank tank 2 tank 3 tank 4 t = time Figure 3: The purge cycles; the rate of drainage from each tank against time. machine as in [3, 4]. Each bowl consists of four tanks separated into wool, cross-flow area and settling tank regions (Figure 2). The model allows the variation of many of the important parameters, such as; wool velocity w wool ; the initial contamination concentration of the wool c w in; the input flow rate of liquor into the bowl q in ; the concentration of input liquor c in ; the particle settling velocity v set ; and the drainage pattern, which specifies the purge cycles. This paper discusses first the mathematical model, then the numerical discretisation method, followed by results and conclusions.

Introduction C24 2 Discretisation model The concentration of contaminants c(x, y, t) throughout the bowl with time t is modelled by the advection-diffusion equation ( ) c 2 t = D c x + 2 c 2 y 2 x (vc) (wc), () y where v is vertical velocity and w is the horizontal velocity of the contaminants and D is a small constant diffusion term that reflects not only pure diffusion, but also some realistic mixing within the fluid due to the harrows. Numerical finite differences are used to discretise the scour bowl into finite points, shown in Figure 2. Each point represents the concentration of contamination in a region around it for the width of the tank. The x discretisation width x is decided by balancing the conflicting aspects of accuracy and calculation time. As there is one column of points in each tank, the y discretisation width y is rather crude. It is assumed the diffusion term D 2 c y 2 in equation () is negligible, given the length scales y, as the effect of horizontal advection w c is far greater than diffusion. For simplicity, y c k j,i = c(j x, i y, k t) (2) is used to denote the contamination at each discretised point, and similarly for other variables. Before the advection-diffusion equation is solved, the velocities v and w need to be found for all times. The vertical velocity v takes into account both the flow rate of the liquor and the settling velocity of the contamination v set. The vertical velocity at a particular height is calculated by dividing the drainage flow rate d k i by the cross-sectional area A j of the tank at that height, to find the flow s velocity, and adding the settling velocity: v k j,i = dk i A j + v set. (3)

2 Discretisation model C25 The horizontal flow rates through the cross flow area, q i,i+ from tank i to i+, must be calculated given a known input flow rate q in and the known rates of drainage d k i. The flow rates are found by the simultaneous solving of a set of conservation of mass equations, one for each tank. For a bowl of four tanks the four equations are q in q 2 d =, q 2 q 23 d 2 =, (4) q 23 q 34 d 3 =, q 34 q out d 4 =, (5) where q out is the unknown output flow rate from tank 4. Once the flow rates are known the horizontal velocities are found by dividing by the horizontal cross-sectional area of the cross-flow area B w k j,i = q i,i+ B. (6) To solve this system satisfactorily, the bowl is divided into three regions, the wool region, cross-flow area and the settling tanks. The regions in each tank are solved with a slightly different numerical method. In the settling tank region there is no flow in the y direction, so, for a rectangular settling tank, the advection-diffusion equation simplifies to c t = D 2 c x (vc). (7) 2 x However, in a more realistic pyramid shaped tank, the concentration of contamination increases as the the cross-sectional area decreases down the tank. An area-weighted advection-diffusion equation, which takes into account the shape of the tank is c t = D 2 c x 2 + A A D c x A A vc v c x, (8) where A is the cross-sectional area of the tank at height x and A is its derivative.

2 Discretisation model C26 The boundary condition at the top of the region is simply the continuous contamination flow down from the cross-flow area. The boundary condition at the bottom depends on the drainage regime in the tank. At times when there is no drainage occurring from the tank, the no flux boundary condition is D c vc =. (9) x At other times, when drainage does occur, the boundary condition c = x allows contamination to advect unhindered through the bottom. Tests were run on the accuracy of second to fifth order finite difference schemes for the settling tanks, comparing results with simple exact solutions. A fourth order scheme was chosen because it gives a satisfactory degree of accuracy at the bottom boundary and is not too computationally intensive. The cross-flow area is the region of the scour bowl between the wool and the settling bowls. In this region the scour liquor flows both horizontally from tank to tank and vertically down from the wool and into the settling tanks as they drain, so the full two dimensional advection-diffusion equation is used, as in equation (). The equation is easiest to solve in this region because it has no complicated boundary conditions, so it can be numerically approximated with first order finite difference methods for time and y-space and a second order method for x-space. The boundary conditions on the top and bottom are obtained by continuity with the settling tank region and the wool region. The boundary condition at y = is c(x,, t) = c in () where c in is the concentration of contamination of the input flow into the cross-flow area. As there is only a first order derivative in y, c/ y, the first order finite difference method used in y-space means no boundary condition is required at the end of the bowl. Unlike the rest of the bowl which is discretised into points x apart, the wool is represented by a single point in each tank. This is reasonable

2 Discretisation model C27 because the wool region is well mixed by the action of the harrows. This also makes the no-flux boundary condition on top of the wool easier to apply; with the same discretisation as the rest of the bowl, this region would have required a more accurate and computationally intensive method than has been used. The governing equation for the wool region is a simple mass balance of incoming and outgoing contamination, c t = F (, y) H where H is the height of the wool, and w wool (c(x, i y) c(x, )), () y F (x, y) = D c v(x, y)c(x, y) x represents the flux of contamination between the wool and the scour liquor. In the finite difference scheme, equation () becomes c k+ n,i = w ( wool t c k n,i + w ) wool t c k n,i t y y H F n,i k. (2) The only boundary condition is at the beginning of the wool to set the contamination of the incoming wool. Various values of x, t are used to generate results, with these terms modified to maximise convergence and minimise iteration time depending on the parameter values. 3 Results Figure 4 illustrates the typical results the model generates, showing concentration of contamination c in the first tank versus time t. These results are, as indicated by the schematic beside the graph, of the concentration at the bottom of a rectangular shaped tank. The three lines represent simulations

rge 3 Results C28 2 4 6 8.5 c = concentration (kg/m 3 ).5 long purge typical purge short purge.5.5.2.4.6.8 -.5 2 4 6 8 t = time (mins) Purge Cycle Figure 4: The build up of contamination at the bottom of the first tank with three different purge cycles.

nts 3 Results C29.5.2.4.6.8 Concentration of contamination (kg/m 3 ).4.3.2. Wool velocity (top) Fitted function Data points w 5 5 2 25 3 Velocity (m/min) Figure 5: The effect of wool velocity on the contamination of the wool in the first tank at steady state. of three different purge cycles. In these simulations the tank begins clean of contamination at time t =, but with increasing time contamination settles down the bowl and concentration builds up, eventually reaching a quasi steady state. This is particularly clear in the simulation with the longest purge. The purge cycle for each simulation is shown beneath the graph and the effect of each purge can be seen as sharp decreases in contamination, as contaminants are flushed out. In the simulations, the more of the purge cycle that is spent purging, the lower the accumulated contamination in the tank. Often, only the steady state concentration in a simulation is important, not what happens beforehand or the effect of the drainage pattern. To find the steady state concentration, a simulation is run until steady state

3 Results C3.5..2.4.6.8 Concentration of contamination (kg/m 3 ).2.8.6.4.2 Tank Tank 2 Tank 3 Tank 4 Fitted function.2.4.6.8 Settling velocity (m/min) Figure 6: The effect of settling velocity of the contamination at the bottom of the four tanks at steady state. v is reached and then averaged over the last purge cycle. This is especially useful when assessing the effect of varying the parameters. Figure 5 shows the steady state concentrations in the wool region as the velocity of the wool w wool is varied, with data from eleven simulations. At w wool =, the contamination is zero as the wool is stationary and will eventually be completely cleaned. At the other extreme, as w wool, the wool is moving so fast that none of the contamination is cleaned and the contamination in the wool remains at the boundary condition, in this case c =.4. Figure 6 illustrates the effect of varying particle settling velocity, showing the concentration at the bottom of each of the four tanks against settling velocity v set. The fitted functions are best fit quadratics. Slower settling

3 Results C3.5.4.2.4.6.8 Concentration of contamination (kg/m 3 ).2.8.6.3.45.4.6.2.75.9 2 4 6 8 Distance along the bowl (m) v Figure 7: The contamination at the bottom of the bowl at steady state for different settling velocities. particles are carried along by the cross-flow and settle in later tanks, while faster settling particles move quickly into the earlier tanks. So, each settling tank collects particles that settle at a particular speed more than others. For example, in tank four the optimal velocity is about.7 m/min. Figure 7 shows the same data as Figure 6, but arranged as concentration against distance along the bowl for different settling velocities. For all the velocities shown, the particles preferentially settle to one distance along the bowl and their concentration decays on either side. Particle settling velocity is linked to particle size, so Figure 6 could be seen as representing the distribution of particle sizes collected in each tank, assuming there is an equal amount of each particle size in the wool.

3 Results C32 Some of the results from this model can be compared with another simple model of the wool in the scour bowl discussed in [2]. 4 Conclusion The movement of contamination, especially dirt, in a wool scour bowl can be modelled using discretised advection-diffusion equations, so as to give time dependent results, revealing, for the first time, how contamination builds up in the tanks and the effect of the drainage pattern. The steady state results show the relationships between the variable parameters and contamination in each tank and at different heights in the bowl. In particular, the relationship between particle settling velocity and contamination at the bottom of the tanks reveals the distribution of particle sizes that collect in each tank, which could be used to set appropriate operating parameters for wools containing different sized contaminants. The relatively complex model has revealed some simple relationships between the operating parameters of the scour machine and contamination throughout the wool and scour bowl. References [] S. I. Barry, T. R. Marchant and G. N. Mercer. Grease recovery and dirt removal in wool scouring. Mathematics in Industry Study Group, Proceedings, ed: J. Hewitt and K. White, 22. C22 [2] J. F. Caunce, S. I. Barry, and G. N. Mercer. A simple mathematical model of wool scouring. In Andrew Stacey, Bill Blyth, John Shepherd, and A. J. Roberts, editors, Proceedings of the 7th Biennial Engineering Mathematics and Applications Conference, EMAC-25, volume 47 of ANZIAM J., pages C34 C47, June 26.

References C33 http://anziamj.austms.org.au/v47emac25/caunce2 [June 22, 26]. C32 [3] D. B. Early. A theoretical analysis of equilibrium conditions in the continuous scouring of raw wool. J. Textile Inst., 68(2/3), 6 67, 978. C22, C23 [4] G. F. Wood. Aqueous jet scouring of raw wool. Part IV. Continuous scouring The approach to equilibrium. J. Textile Inst., 58, 49 53, 967. C22, C23