Sheldon Green


Research Interests

building ventilation

Relevant Degree Programs



Complete these steps before you reach out to a faculty member!

Check requirements
  • Familiarize yourself with program requirements. You want to learn as much as possible from the information available to you before you reach out to a faculty member. Be sure to visit the graduate degree program listing and program-specific websites.
  • Check whether the program requires you to seek commitment from a supervisor prior to submitting an application. For some programs this is an essential step while others match successful applicants with faculty members within the first year of study. This is either indicated in the program profile under "Requirements" or on the program website.
Focus your search
  • Identify specific faculty members who are conducting research in your specific area of interest.
  • Establish that your research interests align with the faculty member’s research interests.
    • Read up on the faculty members in the program and the research being conducted in the department.
    • Familiarize yourself with their work, read their recent publications and past theses/dissertations that they supervised. Be certain that their research is indeed what you are hoping to study.
Make a good impression
  • Compose an error-free and grammatically correct email addressed to your specifically targeted faculty member, and remember to use their correct titles.
    • Do not send non-specific, mass emails to everyone in the department hoping for a match.
    • Address the faculty members by name. Your contact should be genuine rather than generic.
  • Include a brief outline of your academic background, why you are interested in working with the faculty member, and what experience you could bring to the department. The supervision enquiry form guides you with targeted questions. Ensure to craft compelling answers to these questions.
  • Highlight your achievements and why you are a top student. Faculty members receive dozens of requests from prospective students and you may have less than 30 seconds to peek someone’s interest.
  • Demonstrate that you are familiar with their research:
    • Convey the specific ways you are a good fit for the program.
    • Convey the specific ways the program/lab/faculty member is a good fit for the research you are interested in/already conducting.
  • Be enthusiastic, but don’t overdo it.
Attend an information session

G+PS regularly provides virtual sessions that focus on admission requirements and procedures and tips how to improve your application.


Master's students
Doctoral students
Postdoctoral Fellows
I support experiential learning experiences, such as internships and work placements, for my graduate students and Postdocs.
I am open to hosting Visiting International Research Students (non-degree, up to 12 months).
I am interested in hiring Co-op students for research placements.

Graduate Student Supervision

Doctoral Student Supervision (2008-2018)
Paper drying : experimental studies on the influence of dryer fabric (2018)

This research is an experimental investigation on paper drying that primarily focuses on the effects of the dryer fabric on the drying process of paper. A novel method for moisture content measurement is presented. The working principle of this method is the strong correlation between the optical transparency of paper and its moisture content due to the refractive index matching role of water in wet paper. Spectrographic and microscopic measurement techniques were employed to characterize the relation of moisture content and relative transparency of paper. As optical access to the paper is restricted by the dryer fabric, the optical transparency of paper should be measured only with one-sided optical access. To achieve this goal, a novel technique of transmittance measurement is developed that is able to determine the transparency of thin film objects (i.e. paper) with only one-sided optical access. Employing a fluorescence imaging method, this optical configuration eliminates the spurious effect of reflection of the incident light by filtering the excitation wavelength before reaching the optical detector.To study the paper drying process in a multi-cylinder dryer, an experimental setup is designed to simulate realistic conditions of a typical paper dryer while providing optical access for the measurement system. Ten commercially available fabric types manufactured by weaving synthetic filaments are used in the investigations. It is shown that the fabric structure affects the drying progression and the drying time significantly. The contact area and three-dimensional arrangement of the filaments have the greatest impact on the drying process. To study through air drying (TAD), another experimental apparatus is designed to perform drying under controlled conditions of air temperature and mass-flowrate. Four commercially available TAD fabrics with different structural designs and characteristics are used in the investigations. It is shown that the geometry of the contact spots of the fabrics has a significant impact on the drying time at high drying intensities. Comparing the spatial maps of moisture content with the paper grammage distribution reveals that there is a correlation between the local grammage and the local moisture in a paper sheet during the drying process.

A novel SPH method for investigating the role of saliva in swallowing using 4D CT images (2017)

The thesis presents novel computer methods towards simulation of oropha-ryngeal swallowing. The anatomy and motion of the human upper airwaywas extracted from dynamic Computed Tomography (CT) data using a noveltool and workflow. A state-of-the-art SPH method is extended to accommo-date non-Newtonian materials in the extracted geometries. A preliminarynumerical experiment of six human oropharyngeal swallows using SmoothedParticle Hydrodynamics (SPH) demonstrates that the methods are robustand useful for simulation of oropharyngeal swallowing.The presence of saliva is well known to be important for mastication,swallowing, and overall oral health. However, clinical studies of patientswith hyposalivation are unable to isolate the effect of saliva from other con-founding factors. The simulation presented in this thesis examines fluidboluses under lubricated and non-lubricated boundary conditions. Upon comparison with medical image data, the experiments suggest that salivadoes not provide a significant lubricative effect on the bolus transit times,but it may serve to reduce residue and therefore improve overall swallowingefficacy. Our findings, while preliminary, corroborate with existing clinicalresearch that finds that groups with hyposalivation do not have significantlydifferent transit times with control groups, but that residue may be increased in the hyposalivation group.Previous studies using computer simulation of fluid flow in the orophar-ynx typically make use of simplified geometries. Our work uses dynamic320-row Area Detector Computed Tomography (ADCT) images as the ba-sis for the simulations, and therefore does not require simplifying geometricassumptions. Since the data are dynamic, motion trajectories are all sup-plied by the ADCT data, and extrapolation from 2D sources such as bi-planevideofluoroscopy is not required. Processing the image data required the de-velopment of a novel workflow based on a new tool, which we call BlendSeg.We utilize and extend Unified Semi-Analytic Wall (USAW) SPH methodsso that orophrayngeal swallowing simulations may be performed. Theseextensions include the simulation of non-Newtonian boluses, and moving3D boundaries. Partial validation of the extended USAW SPH method isperformed using canonical flows.

Air filtration : predicting and improving indoor air quality and energy performance (2015)

Air filtration is used to reduce particle concentrations in the indoor environment to provide improved occupant health due to reduced exposure. Increased focus on occupant health in emerging design standards is leading to the installation of higher efficiency filtration systems. These systems generally have higher resistance to flow and therefore impose a greater energy penalty. Previous air filter models have used simplified assumptions with regards to the dynamics of filter operation, which have limited the potential to determine energy efficiency or optimization approaches to system design and operation. This dissertation focuses on developing an improved air filter model to investigate the potential for system modifications to reduce energy consumption and improve indoor air quality (IAQ) within commercial buildings.A new air filter performance model was developed using generalizable results from ASHRAE Standard 52.2-2012 and validated against laboratory and real-world experiments. The results showed better agreement with laboratory tests than with real operation. The filter model was combined with existing indoor particle dynamics and epidemiological models to determine the impacts of changes to system operation through monetization of operation costs and health benefits. Laboratory experiments were performed to evaluate the role that particle properties and relative humidity play in determining the filter performance changes with the aim of better understanding the reasons for discrepancies in operation between laboratory and field filter tests.Operation can now be optimized by accounting for dynamic characteristics of filter performance. Benefits of improved filtration efficiency were found to outweigh added costs. Adopting specific indoor particle concentration limits is recommended to replace existing specifications relying on filter efficiency. System designs can then be optimized to account for local particle concentration and energy costs. A number of system design changes have been highlighted that allow for simultaneous reduction in operation cost and indoor particle concentrations. Relative humidity has been identified as a critical parameter in filter performance and standardized tests should be modified to account for variability in relative humidity and particle characteristics typical of real operation to allow for improvements to future model predictions. Supplementary materials:

Modeling the fluid-structure interaction of the upper airway : towards simulation of obstructive sleep apnea (2014)

Obstructive Sleep Apnea (OSA) is a syndrome in which the human Upper Airway (UA) collapses during sleep leading to frequent sleep disruption and inadequate air supply to the lungs. OSA involves Fluid-Structure Interaction (FSI) between a complex airflow regime and intricate mechanics of soft and hard tissue, causing large deformation of the complicated UA geometry. Numerical simulations provide a means for understanding this complex system, therefore, we develop a validated FSI simulation, composed of a 1D fluid model coupled with a 3D FEM solid solver (Artisynth), that is applied to a parameterized airway model providing a fast and versatile system for researching FSI in the UA.The 1D fluid model implements the limited pressure recovery model of Cancelli and Pedley [28] using a dynamic pressure recovery term, area function corrections allowing complete closure and reopening of fluid geometries, and discretization schemes providing robust behavior in highly-uneven geometries. The fluid model is validated against 3D fluid simulations in static geometries and simple dynamic geometries, and proves reliable for predicting bulk flow pressure. Validation of simulation methods in Artisynth is demonstrated by simulating the buckling, complete collapse, and reopening of elastic tubes under static pressure which compare well with experimental results.The FSI simulation is validated against experiments performed for a collapsible channel (a "2D" Starling resistor) designed to have geometry and characteristics similar to the UA. The observed FSI behaviors are described and compared for both experiment and simulation, providing a quantitative validation of the FSI simulation. The simulations and experiments agree quite well, exhibiting the same major FSI behaviors, similar progression from one behavior to another, and similar dynamic range.A parameterized UA model is designed for fast and consistent creation of geometries. Uniform pressure and dynamic flow FSI simulations are performed with this model for numerous parameters associated with OSA. Uniform pressure simulations compare well to clinical data. Dynamic flow results demonstrate airflow limitation and snoring oscillations. The simulations are fast, simulating 1 s of FSI in 30 minutes. This model is a powerful tool for understanding the complex mechanics of OSA.

Dispersion of expiratory droplets in a model single patient hospital recovery room with stratified ventilation (2013)

No abstract available.

Investigation of fibre interaction with a forming fabric (2013)

The forming stage of the papermaking process plays an important role in shaping the quality of final paper sheets. This thesis focuses on studies of fibre motion in the forming section. Wire mark was investigated both numerically and experimentally. Initial sheet forming was simulated with hundreds of fibres of random initial distribution placed into the flow above the fabric and advected onto the fabric. The surface roughness of the resulting fibre mat was calculated. The results show that during initial formation, topographic wire mark is caused in part by fibre bending and in part by the geometry of the fabric. For the specific fibres and sinusoidal forming fabric examined, more than 50% of topographic wire mark was caused by fabric geometry, with the remainder caused by fibre bending. In the experiments, the surface roughness of paper sheets made from different fibre properties was studied using an X-ray tomography device. Light-weight fibre mats were made in a handsheet former machine. A surface map of the wire side of the paper was produced via image analysis. The results reveal that increasing fibre coarseness decreases the surface smoothness of paper. As fibre length increases, surface roughness decreases slightly. Both fine and coarse forming fabrics were used in the sheet forming section. The surfaces of fibre mats made from finer forming fabrics were found to be smoother.The fibre orientation distribution of final paper sheets is closely dependent on the physical properties of the sheets. Fibre orientation in the forming section was studied numerically. In the simulations, one end of each fibre was held by the wire/fibre mat, with the other end carried in the flow. In the uniform flow, analyzed solution from the analysis was obtained. The fibre angle after deposition was only found to be a function of flow direction and initial fibre position. In the shear flow, a dimensionless group of γL/u_z was defined. As the value of γL/u_z increases, fibre mats increase in anisotropy. Fibre properties such as flexural rigidity and aspect ratio were found to have a insignificant effect on fibre orientation.

Three-dimensional flow through forming fabrics and the motion of flexible fibres in flow (2010)

This dissertation addresses the 3D flow through forming fabrics. In the first part of this thesis, the single phase 3D flow through certain specific fabrics was modeled. In practice, the Reynolds number of the fabric flow, based on paper-side filament diameter, is around 100. Consequently, it is not too surprising that the permeability of the fabrics was found to vary approximately as Re ⁻⁰˙⁴, which is intermediate between the expected low Re (Re⁻¹) and high Re (Re⁰) limits. The resistance of a multilayer fabric was found to be nearly equal to the sum of the resistance of each layer considered in isolation. The effect of filament-scale and weave-scale flow non-uniformity on the fiber distribution in the finished paper was considered. For one specific fabric, there was 3 times more chance for short fibers to accumulate initially over openings than blockages of the fabric. Jet-to-wire speed ratio was found to have an insignificant effect on permeability results, but a marked effect on the Machine Direction shear stress in the vicinity of the paper-side filaments. In an attempt to model sheet formation, numerical simulations of the motion of a single fiber in the flow field of a cylinder was carried out as a canonical test case of the fiber/filament interaction system. Seven dimensionless groups govern the problem. A range of dimensionless groups were found for which the fiber hung up on the cylinder, whereas for other values of the dimensionless groups the fiber slid over the cylinder. In general, longer and more flexible fibers had a greater likelihood to be caught by the filament. Yawed finite aspect ratio cylinders at moderate Reynolds numbers are approximate representations of fibers in the flow field of forming fabrics. No analytical solution or any experimental data are reported in the literature to predict the drag and lift force on such particles in such flows. Computational Fluid Dynamics (CFD) simulations were conducted to find the drag and lift coefficients of inclined finite circular cylinders at Reynolds numbers in the range 1-40. The simulations showed that the Independence Principle was highly inaccurate for low inclination angles.

Master's Student Supervision (2010-2017)
Understanding the wheel/rail transfer mechanism in liquid friction modifier carry-down (2017)

In the rail road industry liquid friction modifiers (LFM’s) are used on thetop of rail (TOR) between the wheel/rail interface to reduce curve noise,lateral forces, rail wear and fuel consumption. The friction modifier may beapplied to the rail via a track side applicator and is carried down by thetrain into curved sections of the track where the greatest benefit is seen.A custom laboratory scale machine was designed and built for the purpose of conducting experiments to study the behaviour of LFM carry-downover a large number of wheel/rail interactions. The machine was also designed so that the film transfer at the wheel/rail interaction location couldbe studied.The use of a fluorescent agent to enhance the ability to visualize LFMcarry-down showed promising results, enabling small amounts of carry-downthat couldn’t otherwise be seen under ambient light conditions to now beseen under fluorescence.Qualitative experiments using the machine were performed showing thatan increase in the wheel speed results in an increase in the amount of frictionmodifier transferred from the rail to the wheel at the initial pickup location,thus increasing the carry-down. Increasing the applied load had the oppositeeffect and reduced the amount of friction modifier initially transferred fromthe rail to the wheel, and thus reducing the carry-down. The profile of thewheel was observed to effect the initial transfer amount and the ensuingcarry-down due to high/low pressure zones along the wheel/rail interface.

Full scale train underbody aerodynamic evaluation for top of rail friction modifier application (2015)

Liquid jet impingement is employed in the rail industry to apply friction modifier to the rail surface. Use of the friction modifier is known to reduce wear and improve fuel efficiency. L.B. Foster ® deploys friction modifier using a nozzle located downwind of the wheels on freight trains. Understanding the aerodynamic environment of the nozzle is important for researching how to maximize the deposition of the liquid friction modifier from the nozzle to the tracks. The air pressure and velocity at the location of the nozzle was evaluated experimentally at full scale in field trials. The pressure at a fixed ground location was measured by transducers as the train passed. The air velocity in the reference frame of the moving vehicle was measured using a fiber-film anemometer at the location of the liquid-friction-modifier spray nozzle, 0.4 wheel diameters downwind of the wheel center.The measured air speeds scales linearly with the train speed, and the measured pressure scales linearly with the dynamic pressure, implying that Reynolds number effects are negligible. The pressure distribution showed an initial pressure increase just downwind of the leading edge of the vehicle followed by a spike in suction. The pressure distribution was found to depend on the orientation of the vehicle. With a rail car leading the vehicle, the spike in suction produced was about 50% larger than the suction spike produced when a locomotive, lower to the ground, was leading the vehicle. The mean air speed was measured to be approximately 29% of the train speed. The mean air speed the same distance upwind of the wheel was measured to be approximately 38% of the train speed. Turbulence intensity levels were measured to be about 0.15. Cross wind effects became much less significant when the train speed was equal to or greater than the cross wind speed.The train undercarriage airflow was modeled numerically using Autodesk Simulation CFD™ software. The CFD simulations were in approximate agreement (typically, within 2%) with experimental measurements and confirmed that the presence of the support bracket for the anemometer had limited impact on the measured wind speed.

Newtonian and viscoelastic liquid jet impingement on a moving surface (2015)

Motivated by the need to improve transfer efficiencies of liquid coatings from jet impingement in railroad friction control applications, an experimental investigation into Newtonian and viscoelastic liquid jet impingement on moving surfaces is presented.Seven PEO-glycerine-water solutions and three commercial liquid friction modifiers were tested with a variety of jet speeds, jet diameters, surface speeds and surface roughnesses. The effects of these test conditions on jet impingement splash behaviours as well as jet and lamella geometries were studied. High-speed imaging was employed to visualize the interaction between the impinging jet and the moving surface.Experiments on the effect of modest surface roughness revealed that, while jet and surface speed were both important factors, splash was more likely to occur on surfaces with lower roughness levels. By analyzing experimental results for Newtonian liquids, a relation between lamella geometry and test conditions were found, which can be used to predict lamella dimensions. Three types of non-Newtonian behaviours were observed at high surface speed and low jet speed: jet necking, jet bending and jet stretching.

Transient jet impaction on a high speed moving surface (2015)

In the railroad industry liquid friction modifiers are applied on the rail track in the form of a liquid jet in order to reduce the friction and fuel consumption. In this application, the transfer efficiency of the liquid on the rail track is very important. To maximize the transfer efficiency, Newtonian and non-Newtonian transient liquid jet impingement on a dry moving surface was studied. Five different water glycerin solutions with widely varying shear viscosities were used as Newtonian test liquids to isolate the effect of shear viscosity from other fluid properties. Furthermore, the effect of surface roughness on the impaction was investigated using four different roughness heights. The effects of jet velocity and surface speed were also studied. High speed imaging was performed to visualize the interaction between the jet and the moving surface. For surface roughness heights between 0.02 μm and 0.64 μm, it was found that as the roughness increases, the jet becomes more prone to splash. It was also shown that increased jet and surface speeds trigger the splash. The transient jet characteristics were also investigated for Newtonian liquids at different nozzle back pressures. It was found that at higher Reynolds and Weber numbers the transient jet breaks up downstream of the nozzle; However, it was shown that the Weber number has the dominant role in jet break-up compared to the Reynolds number. A numerical study was also undertaken to determine the drag force exerted on the plunger of the solenoid valve in the nozzle. The simulation results were in reasonable (16% on average) agreement with experiment.

Micro-PIV measurement of flow through forming fabrics (2014)

This thesis describes an experimental investigation of the flow field upstream of forming fabrics that are typically used in the paper making process. Micro Particle Image Velocimetry was used to measure the velocity distribution upstream of a forming fabric. The velocity upstream of two different types of forming fabrics, namely Monoflex D60TM and IntegraTM, was studied. As expected, the experiments show the existence of a highly variable drainage velocity field upstream of both fabrics. The drainage velocity over the holes can be several times greater than the drainage velocity above the fabric filaments. Since fines and filler tend to follow fluid streamlines, one would therefore expect substantially higher fines and filler concentrations in the holes between the filaments as compared over fabric knuckles. The decay in drainage velocity variations can be represented by the equation Aexp(-Bz/D)+C, where A and B are constant and C is the uncertainty in the experimental setup. D and z represent the fabric’s filament diameter and distance above the fabric surface.It is expected that the response of pulp fibers to the velocity variations caused by the fabric’s weave structure is strongly correlated to their length. The fibers with a length greater than 1.5 mm experience a weighted-average velocity field along their length that is approximately uniform. The deposition of short fibers with length

Design of an optical uroflowmeter and assessing bladder pressure through video analysis of the male urine stream (2013)

Voiding dysfunction, such as benign prostatic hyperplasia and impaired detrusor contractility, affects more than half of men over the age of 50. Uroflowmetry provides quantitative information of flow dysfunction by measuring the flow rate and total volume of urine expelled by the body. The majority of existing clinical uroflowmeters determine flow rate using a scale to measure increasing mass of urine expelled with time; however, they are expensive and typically found only in specialists’ offices, making it difficult for patients to receive testing. An opportunity therefore exists to develop a much more affordable device which would allow flow rate testing to become a part of routine care and to be conducted in a wider variety of environments such as General Practitioners’ offices and home monitoring. The high demand for digital cameras, particularly due to their extensive use in mobile devices, has resulted in their accelerated advancement and cost reduction. Therefore, a device based on this technology is investigated. In addition to the development of this device, a study was conducted to investigate if information regarding bladder pressure may be obtained by analyzing digital images of the urine stream. Voiding dysfunction may result in abnormally high bladder pressure, caused by urinary obstruction, or low bladder pressure which may be caused by reduced detrusor contractility. The clinical implications and treatment for these two cases are very different; however, they present with similar low flow rate voiding patterns and cannot currently be distinguished non-invasively. It was postulated that increased inertia and turbulence may exist in high pressure flows, and may be identifiable in the digital images of the urine stream.

Hydrodynamic interaction between cylinders at moderate Reynolds numbers (2013)

The hydrodynamic interaction between two cylinders perpendicular to the freestream, in a tandem arrangement was studied for moderate Reynolds numbers (1≤Re≤40). The influence of multiple geometric variables was considered: separation distances between the cylinders, ellipticity of the cylinders, the cylinder aspect ratio, and the angular inclination between the cylinders. In the first part of this study, a numerical investigation of the two-dimensional steady flow past cylinders was carried out. The characteristic length, D, in all simulations was taken to be twice the major axis of the cylinder cross-section, (i.e. equal to the diameter for cylinders of circular cross-section). The two-dimensional flow was studied for separations up to 50D. Four different ellipticities were studied. The drags experienced by front and rear cylinders were compared with that experienced by a single cylinder of the same cross-section. The second part of the study consisted of the steady three-dimensional flow analysis for parallel cylinders in tandem for separations ranging from 2D to 20D and cylinder lengths up to 20D. In the third part of this thesis, a steady flow analysis was done for two circular cylinders in tandem with lengths equal to 5D but with the cylinder axes in different orientations relative to the plane normal to the flow. This angular separation between the cylinders produces a hydrodynamic moment, which is dependent on the geometry and the flow Reynolds number.The fourth and final part of this work is the study of the unsteady three-dimensional flow that would result from the hydrodynamic moment discussed in relation to the third part of the thesis. The thesis closes with some remarks on the implications of these findings to papermaking and recommendations for future work.

An experimental study of liquid jet impingement on a moving surface : the effects of surrounding air pressure and fluid properties (2012)

The impingement of a high-speed liquid jet on a moving surface and the resulting deposition or splash is important in a variety of technical and industrial processes. Of particular interest is the coating of the top-of-rail surface, in the rail road industry, with a thin film of viscoelastic liquid friction modifier, by liquid jet impingement, to control friction and reduce wear at the wheel-rail interface, thereby reducing fuel consumption and maintenance costs. For effective operation it is required that the fluid deposited by the jet adhere to the surface after impingement. An experimental investigation into the effect of surrounding air pressure and fluid properties on liquid jet impingement on a moving surface was performed. The study was carried out with Newtonian liquids impacting smooth, dry surfaces. A variety of ambient air pressures, jet speeds, surface speeds, surface tensions, and liquid viscosities were studied. The interaction between the impinging jet and the moving surface was analysed through high-speed imaging. It was observed that, as is the case for Newtonian droplet impact, the surrounding air pressure plays a crucial role in the splashing behaviour of jet impingement. There exists a threshold pressure below which splash does not occur. It is proposed that for certain impingement conditions lamella detachment from the surface occurs due to aerodynamic forces acting on the leading edge of the lamella, which destabilizes the balance between surface tension and fluid pressure forces. It was observed that both the Reynolds number and Weber number were salient to the occurrence of lamella detachment, with lamella detachment having a non-linear dependence on the Reynolds number. Lamella detachment was prone to occur for intermediate Reynolds numbers as the Weber number was increased, bounded by regions of deposition at higher and lower Reynolds numbers.

An experimental study of the surface distribution of filler material in paper (2012)

A series of experiments were conducted to observe the effect of different pulp suspension and formation characteristics on the variation in filler concentration on the surface of paper. Hand sheet samples were formed in a laboratory apparatus. The surface distribution of two types of filler material was investigated: Precipitated Calcium Carbonate (CaCO₃) and Kaolin Clay (Al₂Si₂O₅(OH)₄). The effect of retention aids, dewatering rate, and forming fabric geometry on filler distribution was tested.The analysis was focused on the variation in surface filler concentrations on the scale of individual strands of the forming fabric. The procedure involved locating areas of interest in the samples, particularly the area of paper formed over knuckles, threads, and openings in the fabric, at which point the sample was analysed with Energy Dispersive X-Ray Spectroscopy and image analysis techniques to determine relative filler concentrations.In samples formed by gravity and vacuum drainage, Kaolin displayed a significantly greater variation in local surface filler concentration than PCC, though the difference was reduced under vacuum drainage conditions. This effect is attributed to the electrostatic attraction between PCC and pulp fibres in contrast with the repulsion felt by Kaolin filler. The attraction resists the distributing forces of the flow at low drainage velocities but the bonds are broken by large shear forces. The distribution on the top side of the paper was comparable between the filler types, due to the more uniform flow field at a distance from the forming fabric. Vacuum drainage increased the spatial variation of both fillers by a similar amount. It was found that under vacuum drainage, retention aids did not improve filler uniformity on the wire side. However, on the top side of the paper, a moderate reduction in spatial variation was observed. Additionally, on the wire side of samples made with gravity drainage, it was found that the addition of retention aids produced a significant improvement in the uniformity of the filler material. Finally, it was found that a finer forming fabric improved the uniformity of filler distribution.

An experimental study on jet impingement on a very high speed moving surface (2012)

Motivated by the need to improve transfer efficiencies of liquid coatings from jet impingement, an experimental investigation into jet impingement on very high speed moving surfaces is presented. Seven different Newtonian liquids with widely varying shear viscosities were made to impinge on a surface which could be made to move at speeds up to 350km/hr. Tests for the Newtonian liquids were done with several modified surfaces to study the effects of roughness and surface inconsistencies. Nozzle sizes and impingement angles were varied to interrogate their effects on the interaction of the impacting jet and moving surface while high speed photography was employed to capture these interactions. Spread radii and spread widths were measured for viscous fluids which deposited.While it was observed that stable jets of fluids with sufficiently high viscosities almost always deposited, tests with water indicate that the effects of the impingement angle as well as jet diameter significantly alter the locations of boundaries between deposition, spatter and lamella lift-off. Impingement angles that result in jet velocities with large components of velocity parallel to the surface velocity are prone to deposit. Jets of smaller diameters are also prone to deposit. It was observed that both the jet velocity and surface velocity are important determining factors in the likelihood of deposition.The deposition of viscous fluids demonstrated that it is possible to observe transitions from deposition to lift-off and vice versa through mechanisms that trigger random fluctuations in the lamella. The track distance covered before a transition from lift-off to deposition occurs is shown to be a Poisson Process.

Newtonian and elastic liquid jet interaction with a moving surface (2011)

In the railroad industry a friction modifying agent may be applied to the rail or to the wheel in the form of a liquid jet. In this mode of application the interaction between the high speed liquid jet and a fast moving surface is important. Seven different Newtonian liquids with widely varying shear viscosities along with twelve different solutions of polyethylenoxide (PEO) and water with varying relaxation times were tested to isolate the effect of viscosity and elasticity from other fluid properties. Tests for the Newtonian liquids were done with five surfaces having different roughness heights to investigate the effects of surface roughness. High speed video imaging was employed to scrutinize the interaction between the impacting jet and the moving surface. For both Newtonian and Elastic liquids and all surfaces, decreasing the Reynolds number reduced the incidence of splash and consequently enhanced the transfer efficiency. At the elevated Weber numbers of the testing, the Weber number had a much smaller impact on splash than did the Reynolds number. The ratio of the surface velocity to the jet velocity has only a small effect on the splash, whereas increasing the roughness-height-to-jet-diameter ratio substantially decreased the splash threshold. Moreover, the Deborah number was also salient to the splash of elastic liquids.

PIV measurements of flow through forming fabrics (2011)

Three-dimensional velocity fields in the single phase approach flow to a multiple layer woven forming fabric were measured using Particle Image Velocimetry (PIV). The measurements were conducted on a scale model of a forming fabric in a water/glycerin flow loop. Each strand on the paper side of the model forming fabric had a filament diameter around 15.4mm, and the loop test section was 310mm squared, permitting the measurement of detailed velocity distributions over multiple strands of the fabric. The flow speed in the loop test section were varied to achieve screen Reynolds numbers (Res), calculated based on paper side filament diameter (d), between 15 and 65. PIV measurements showed that the normalized ZD velocity deviation decreases from 19.7% at a plane 0.25d upstream from the forming fabric to 4.2% at a plane 1.5d upstream. The normalized CMD velocity deviation decreases from 15.3% at a plane 0.25d upstream from the forming fabric to 1.9% at a plane 1.5d upstream. The normalized MD velocity deviation decreases from 14.5% at a plane 0.25d upstream from the forming fabric to 2.3% at a plane 1.5d upstream. The highest ZD velocity is about 3.3 times higher than the lowest ZD velocity at a plane 0.25d above the fabric. This ratio decreases to 1.2 at a plane 1.5d above the fabric. These findings show that the flow non-uniformity caused by the fabric weave is restrained to a short distance above the fabric. However, even this non-uniformity is not particularly felt by fibers, whose length scale results in an averaging of the local velocity field. CFD simulations of the same flow were consistent with the PIV measurements.

The role of suction boxes on forming section retention and filler migration (2010)

An apparatus was constructed to observe fibre mat formation under applied vacuum pressure comparable to that experienced in a suction box. The main focus of study was to determine the effects of system parameters on overall retention, filler distribution, and filler migration by analysing samples through weighing, ashing and EDX analysis. It was found that increasing suction pressure slightly decreased the overall retention, although the effect was greater at higher filler loadings. The forming fabric selection was statistically significant in determining the overall retention, and the significance of forming fabric selection grows with increasing filler loading. The most important factor in determining overall retention was the presence of retention aids. Measurement of the instantaneous permeability of the forming fabric and fibre mat during the forming process showed differences in the permeability curves between forming fabrics, particularly during the initial forming process. The permeability variations are correlated with differences in retention. Tests performed under gravity drainage, to provide a correlation between fabric properties and overall retention results, showed that air permeability had the highest correlation.Sheet splitting and ashing hand sheets formed using cationic PCC filler showed that applying a vacuum during the drainage process decreased the average filler content. The filler was preferentially removed from the wire side of the sheet. Similar tests performed with anionic filler showed that filler charge affects the distribution in a hand sheet. Anionic filler also showed little change in distribution when vacuum was applied during forming. Cationic and anionic filler have different retention mechanisms; attachment and filtration, respectively. Applied vacuum was found to affect the filler distribution of particles retained by attachment with the highest changes at the wire side. Applied vacuum was found to have only a small affect on particles retained by filtration. The reason for both is the compaction of the fibre mat, which occurs to the greatest extent at the wire side. The compaction reduces the size of openings in the mat which will further trap large particles retained by filtration but increase the drag on small particles, possibly exceeding the attractive forces responsible for attachment to the pulp fibres.



Membership Status

Member of G+PS

Program Affiliations


If this is your researcher profile you can log in to the Faculty & Staff portal to update your details and provide recruitment preferences.