Laboratory and Numerical Modelling Compared with Field Trials of Dilution at the Sydney Desalination Plant Outfall

Similar documents
AN INVESTIGATION OF LONGITUDINAL VENTILATION FOR SHORT ROAD TUNNELS WITH HIGH FIRE HRR

Experiment of a new style oscillating water column device of wave energy converter

INTERACTION BETWEEN WIND-DRIVEN AND BUOYANCY-DRIVEN NATURAL VENTILATION Bo Wang, Foster and Partners, London, UK

AIRFLOW GENERATION IN A TUNNEL USING A SACCARDO VENTILATION SYSTEM AGAINST THE BUOYANCY EFFECT PRODUCED BY A FIRE

TRIAXYS Acoustic Doppler Current Profiler Comparison Study

A Computational Assessment of Gas Jets in a Bubbly Co-Flow 1

Offshore Wind Energy Stringent quality assurance and quality control. Coastal and Freshwater Fast responding and flexible organisation

Experimental Determination of Temperature and Pressure Profile of Oil Film of Elliptical Journal Bearing

LABORATORY EXPERIMENTS ON WAVE OVERTOPPING OVER SMOOTH AND STEPPED GENTLE SLOPE SEAWALLS

Presented by. Mr.Danish.D.R. M.Tech Coastal Management Institute for Ocean Management Anna University, Chennai Tamil Nadu, India.

AIR-WATER FLOW STRUCTURES AT AN ABRUPT DROP WITH SUPERCRITICAL FLOW

Lagrangian Tracer Transport and Dispersion in Shallow Tidal Inlets & River Mouths

Columbia River Plume 2006 OSU Ocean Mixing Nash, Kilcher, Moum et al. CR05 Cruise Report (RISE Pt. Sur, May )

MODELLING OF FUME EXTRACTORS C. R.

Fine-Scale Survey of Right and Humpback Whale Prey Abundance and Distribution

High Frequency Acoustical Propagation and Scattering in Coastal Waters

Effect of Argon Gas Distribution on Fluid Flow in the Mold Using Time-Averaged k-ε Models

NUMERICAL INVESTIGATION OF THE FLOW BEHAVIOUR IN A MODERN TRAFFIC TUNNEL IN CASE OF FIRE INCIDENT

High Frequency Acoustical Propagation and Scattering in Coastal Waters

MAPCO2 Buoy Metadata Report Project Title:

LABORATORY STUDY ON TSUNAMI REDUCTION EFFECT OF TEIZAN CANAL

Temperature, salinity, density, and the oceanic pressure field

Student name: + is valid for C =. The vorticity

TRANSPORT OF NEARSHORE DREDGE MATERIAL BERMS

ANALYSIS OF THE POSITIVE FORCES EXHIBITING ON THE MOORING LINE OF COMPOSITE-TYPE SEA CAGE

7.4 Temperature, Salinity and Currents in Jamaica Bay

A Tracer Investigation of. Aerated Water Dispersion and Tidal Exchange in. the San Joaquin River and Stockton Ship Channel

HYDRAULICS. H89.8D - Hydraulic Bench

PHYS 101 Previous Exam Problems

Experimental Analysis on Vortex Tube Refrigerator Using Different Conical Valve Angles

Free Surface Flow Simulation with ACUSIM in the Water Industry

EXPERIMENTAL STUDY ON BEESWAX USING WATERJET DRILLING

Currents measurements in the coast of Montevideo, Uruguay

Multiple port vessels.

485 Annubar Primary Flow Element Installation Effects

Greenup Lock Filling and Emptying System Study

Air entrainment in Dip coating under vacuum

Flow transients in multiphase pipelines

Lab #4 Pipe Flow, Minor and Major Losses, and Walking in Osborne Reynolds Shoes CEE 331 Fall 2006

14/10/2013' Bathymetric Survey. egm502 seafloor mapping

The water supply for a hydroelectric plant is a reservoir with a large surface area. An outlet pipe takes the water to a turbine.

ITTC Recommended Procedures and Guidelines

BERNOULLI EFFECTS ON PRESSURE.ACTIVATED W ATER LEVEL GAUGES

Technical Brief - Wave Uprush Analysis Island Harbour Club, Gananoque, Ontario

Robin J. Beaman. School of Earth and Environmental Sciences, James Cook University, Cairns, Qld 4870, Australia.

An experimental study of internal wave generation through evanescent regions

Abstract. The main conclusions from analysis for the data collected are:

Numerical and Experimental Investigation of the Possibility of Forming the Wake Flow of Large Ships by Using the Vortex Generators

DEVELOPMENT OF HIGH ALTITUDE TEST FACILITY FOR COLD JET SIMULATION

Learn more at

THE EXTERNAL EXPLOSION CHARACTERISTICS OF VENTED DUST EXPLOSIONS.

Large-scale Field Test

R/V Point Sur Pre-Cruise Questionnaire

Digiquartz Water-Balanced Pressure Sensors for AUV, ROV, and other Moving Underwater Applications

1. All fluids are: A. gases B. liquids C. gases or liquids D. non-metallic E. transparent ans: C

AIR EJECTOR WITH A DIFFUSER THAT INCLUDES BOUNDARY LAYER SUCTION

Undertow - Zonation of Flow in Broken Wave Bores

Flow in a shock tube

Old-Exam.Questions-Ch-14 T072 T071

Development of High-speed Gas Dissolution Device

CHAPTER 68. RANDOM BREAKING WAVES HORIZONTAL SEABED 2 HANS PETER RIEDEl. & ANTHONY PAUL BYRNE

PRESSURE DISTRIBUTION OF SMALL WIND TURBINE BLADE WITH WINGLETS ON ROTATING CONDITION USING WIND TUNNEL

Long-Term Autonomous Measurement of Ocean Dissipation with EPS-MAPPER

Autonomous Environmental Profiling Moorings for Coastal Monitoring

Exploring Localized Mixing Dynamics During Wet Weather in a Tidal Fresh Water System

Deep Water Currents Lab

Ermenek Dam and HEPP: Spillway Test & 3D Numeric-Hydraulic Analysis of Jet Collision

Hamurana Stream water movement in Lake Rotorua

Chapter 2 Hydrostatics and Control

Helium Mass Spectrometric Leak Detection In Large Size Process Plants

Angela Lane, Lowe Environmental Impact / Karen Akuhata (WDC) The Wairoa wastewater treatment system requires a replacement consent by May 2019.

Aalborg Universitet. Published in: Proceedings of Offshore Wind 2007 Conference & Exhibition. Publication date: 2007

Lab # 03: Visualization of Shock Waves by using Schlieren Technique

. In an elevator accelerating upward (A) both the elevator accelerating upward (B) the first is equations are valid

Experimental Verification of Integrated Pressure Suppression Systems in Fusion Reactors at In-Vessel Loss-of -Coolant Events

Exercise 2-2. Second-Order Interacting Processes EXERCISE OBJECTIVE DISCUSSION OUTLINE. The actual setup DISCUSSION

Offshore engineering science

WAVE OVERTOPPING AND RUBBLE MOUND STABILITY UNDER COMBINED LOADING OF WAVES AND CURRENT

The GCOOS Mooring Plan Element Draft, 19 February 2011

THEORETICAL EVALUATION OF FLOW THROUGH CENTRIFUGAL COMPRESSOR STAGE

Axial Base (PN3A) Medium-Power (MP) J-box

NTHMP - Mapping & Modeling Benchmarking Workshop: Tsunami Currents

STUDY OF SLUG CONTROL TECHNIQUES IN PIPELINE SYSTEMS

INTRODUCTION TO COASTAL ENGINEERING

CORESTA RECOMMENDED METHOD N 6

Argon Injection Optimization in Continuous Slab Casting

Cover Page for Lab Report Group Portion. Drag on Spheres

NSF's Ocean Observatories Initiative: Building Research Infrastructure for the Pacific Northwest and the Broader Community

BILLY BISHOP TORONTO CITY AIRPORT PRELIMINARY RUNWAY DESIGN COASTAL ENGINEERING STUDY

Pressure distribution of rotating small wind turbine blades with winglet using wind tunnel

Chapter 3 EXPERIMENTAL: EQUIPMENT, SET-UP, AND PROCEDURE

Experiment 8: Minor Losses

Development of an Oil Delivery System to Create Neutrally Buoyant Oil Plumes

Wade Reynolds 1 Frank Young 1,2 Peter Gibbings 1,2. University of Southern Queensland Toowoomba 4350 AUSTRALIA

DETRMINATION OF A PLUNGER TYPE WAVE MAKER CHARACTERISTICE IN A TOWING TANK

South Bay Coastal Ocean Observing System California Clean Beaches Initiative

Gravity wave effects on the calibration uncertainty of hydrometric current meters

Mixing and Transport in the Surfzone

Waves, Bubbles, Noise, and Underwater Communications

Undertow - Zonation of Flow in Broken Wave Bores

Transcription:

International Symposium on Outfall Systems, 2016 Organized by: IAHR-IWA Joint Committee on Marine Outfall Systems May 10 th -13 th, 2016, Ottawa, Canada Laboratory and Numerical Modelling Compared with Field Trials of Dilution at the Sydney Desalination Plant Outfall B.M. Miller, G.P. Smith and J.E. Ruprecht. Water Research Laboratory, School of Civil and Environmental Engineering, UNSW Australia ABSTRACT Physical and numerical modelling was undertaken as part of the diffuser design of the Sydney Desalination Plant brine multiport outfall. Subsequent to the construction, field investigations to validate the diffuser performance were undertaken with the desalination plant operating at full capacity. This paper presents the laboratory and field methodologies and provides a comparison between the model predictions and field performance. Laboratory testing of many diffuser configurations were undertaken in a 4.5m x 4.5m x 0.6m tank. Two full days of field dilution testing were undertaken where a known concentration of Rhodamine WT fluorescent dye was injected into the seawater concentrate discharge and the concentrations in the receiving waters were measured using both diver-held and boat-cast fluorimeters. Measurements of the plume depth and geometry were also derived by analyzing the backscatter from ADCP current meter transects. The dilutions in the receiving waters were effected by varying currents. These field observations show how actual conditions can vary from laboratory conditions, but also provide reassurance that the laboratory and numerical predictions are slightly conservative. Keywords: Dense Jets, Field Trials, Physical Modelling, Plumes, Rhodamine, Outfall, Desalination 1 THE SYDNEY DESALINATION PLANT OUTFALL The Sydney Desalination Plant is located at Kurnell towards the southern end of the Sydney coastline. The diffuser has two risers each 25m apart and located approximately 500m offshore in approximately 25m of water (Figure 1). The risers are located on a gradually sloping rocky seabed. Each riser cap is 7m diameter, with four 370mm diameter nozzles angled upwards at 60 degrees from the horizontal. For a plant production of 250 ML/day, the discharge is 4.0 m 3.s -1 and the exit velocity is approximately 4.8 m.s -1. The discharge volume and salinity can vary with production, but the brine salinity is typically between 58 ppt and 62 ppt whereas the receiving water salinity is 35 ppt.

Figure 1. Location of the Outfall Risers Relative to the Coastline. The performance criteria for the diffuser are that the seawater concentrate meets water quality criteria for salinity and relevant treatment chemicals by the edge of the near-field and the visual amenity of the sea surface is maintained. The edge of the near field was specified as the impact point where the brine returned to the seabed. At the end of the nearfield: salinity must be within 1ppt of background; and the dilution must be at least 30 times (as determined by eco-toxicity testing of other chemicals in the brine stream). The adopted visual amenity criteria was that the maximum plume rise height should be at least 2m below the water surface. Due to a change in design during construction, the discharge through each riser nozzle increased and there was no opportunity to increase the spacing between risers. As such the inner jets, being those oriented towards the other riser, would interact at the impact point. 2 LABORATORY TESTING Laboratory testing (Figure 2) was undertaken in a 4.5m x 4.5m x 0.6m tank with similitude of the densimetric Froude number. A 10 cm 10 cm grid placed on the marine plywood floor with string lines aided in observing and undertaking measurements. The receiving waters were filled and drained by an orifice in the tank floor with the tank cleaned and refilled with mains supply water between each test. High definition digital video cameras were mounted to visually record each experiment. Cameras were mounted overhead the tank providing a plan view of the plumes and beside the tank to provide a profile view of the plumes. Brine was mixed into a head tank and discharged into the receiving water tank through valves and rotameters. Twelve monitoring points consisting of a single microelectrodes mi-900 series conductivity electrode were available for the model tests. The electrodes consisted of two small probes and a cylindrical body of 30 mm length and 4 mm diameter, hence minimising the potential for disruption to flow. The electrodes measured the conductance across the two probes using the bipolar current technique (Fan and Brown, 2006) and have a nominal detection limit of 0.002 volts and WRL in-house testing has confirmed that repeat measurements are within this detection limit. The electrodes were mounted in two array frames of six electrodes. Each array was mounted on a brass base, which was heavy enough to ensure the arrays did not move during the experiments. Individual electrodes were then mounted on aluminium tabs with the electrode sensors offset from the brass base. The electrodes were set in line and horizontally spaced at 50 mm, which is equivalent to 3 m at prototype

scale. Sensors were set to be 8 mm above the model floor, which is equivalent to 0.5 m above the prototype sea bed. The arrays were designed to allow measurement of EC both longitudinally and laterally through the plume at the point of impact. Electrical conductivity was related to density through chemistry laboratory tests using the specific electodes and pool salt used in the main experiments. Dilutions were calculated on the basis of density. The modified gravity (g0 ) in the model was kept the same as the modified gravity in the prototype. As such, the density difference between the seawater concentrate and the receiving waters in the ocean were the same as in the physical modelling experiments (i.e. the term Δr = 1). Mixing processes in jets are similar and fully turbulent providing that the Reynolds number (below) has a value greater than 2000 (Fischer et al. 1979). However, it has been reported that this may be as high as 4000, so the Reynolds number in the models tests was maintained above 4000. Figure 2. Laboratory Testing Facility. Many tests were undertaken and are reported in Miller et al (2009). Table 1 presents the various rates, sizes and parameters from the configuration that was finally constructed (prototype) and the scaled model. Individual jets were found to have impact dilution and height of rise as predicted by Roberts (1999). However, the inner jets, being those from each riser that interacted (Figure 3), achieved a lower dilution. The combined dilution results of three repeated tests are presented in Figure 4.

Table 1. Prototype and Model Parameters Prototype Units Model Units Scale 1 59.68 Number of risers tested 2 2 Number of nozzles per riser 4 4 Nozzle diameter 370 mm 6.2 mm Nozzle angle 60 degrees 60 degrees Volume 16688 m3/h 10.1 L/min Water Depth above Seabed 25 m 419 mm Salinity of Brine 62 ppt Density Brine 1046.4 kg/m3 1018.7 kg/m3 Nozzle exit velocity 5.4 m/s 0.70 m/s Density 1026.2 kg/m3 999 kg/m3 Density Difference (brine sea water) 20.2 kg/m3 19.7 kg/m3 Modified Gravity 0.2 0.2 Port Froude Numer 20.1 20.1 Reynolds Number 2.00E+06 4.33E+03 Figure 3. Overhead photograph of laboratory testing of the two risers, showing the four outer jets being separate and the four inner jets interacting.

Dilution 50 45 40 35 30 25 Individual Plume Merging Plume Impact Point 20 15 10 15 20 25 30 35 Distance (m) from Nozzle Figure 4. Results of laboratory testing of the constructed configuration showing the decrease in dilution from merging plumes. With a brine of 62 ppt and receiving waters of 35ppt, a dilution of 26 is required to achieve 1ppt above background. This was predicted by the physical modelling to be achieved at the impact point for the outer plumes and within 26m for the inner plumes. The 30 times dilution target was predicted to occur within 20m for the outer plumes and 33m for the inner plumes. 3 FIELD TESTING 3.1 Methodology The technique to measure outfall plume dilution performance was to inject Rhodamine WT (a fluorescent tracer dye) into the ocean outfall chamber at the Sydney Desalination Plant at a known rate and concentration and then to measure the dye tracer concentration in receiving waters. Our past experience indicated that detailed measurement of the plume in the near vicinity of the outfall risers and seabed at 25m depth would be too difficult using instruments cast from a surface vessel. To enable detailed measurement of the plume close to the outfall risers and the seabed, divers were used to swim the instruments directly into and around the outfall diffuser jets. This allowed for direct and targeted measurement of the plume dilution at the plume impact point and from the plume impact point through the intermediate field into the far-field provides a direct assessment of the increase in plume dilution as the plume transitions from the near-field into the far-field. Rhodamine WT was adopted as it is not naturally occurring and detection cannot be confused with a background concentration of any other tracer that might be considered. Experimental dye concentrations of 0 5 ppb are safely within the adopted instrument detection accuracy (a sensitivity of 0.01 ppb).

The first vessel was a dive vessel with a qualified commercial dive team on board. This vessel was moored between the two risers to enable a stable point to conduct the diving. The second vessel was a roving vessel equipped with a surface suspended flourimeter, an Acoustic Doppler Current Profiler (ADCP) and a Seabird 19Plus for taking salinity and temperature profiles. The divers were cabled back to the surface so that engineers could see measurements in real time. In addition to the helmet mounted video equipment, a WetLabs ECO fluorometer was attached to the divers Buoyancy Control Device (BCD) and was used to sample the dye concentration in the plume. The diver was directed by engineers on board the dive vessel to make direct and targeted measurements of the plume rise height and plume impact point by swimming into and around the plume jets. The diver s location at each measurement point was noted using a distance tagged grid line attached to the outfall diffuser along with a diver recorded compass bearing and a digital depth gauge. Plume measurements by the diver were focussed on the near-field around the plume impact point. The diver also obtained grab samples as a secondary means to confirm dilution. Figure 5 shows photographs from the field work. 3.2 Results Field work was undertaken on the 29 th and 30 th March 2011. Although there were some differences in the findings, the results were similar on both days. Due to space limitations in this paper, only results from Day 1 are presented here, but the full results can be found in Smith et al (2011). The desalination plant production rate was set at 250 ML/day. The discharge to the outfall remained constant (within +/- 1.5%) at 4.00 m 3 /s, 24.1 o C and 57.9ppt. Rhodamine WT was injected into the outfall flow at a rate of 7.1 L/h, providing an active Rhodamine WT concentration in the seawater concentrate flow of 98 ppb. Rhodamine WT dye injection at the plant began at 9:16 a.m. and ceased at 12:39 p.m. (3 hours 23 minutes). Calm offshore winds prevailed and wave conditions were slight (<0.5m). The tide fell from +0.9m to +0.5m above mean sea level throughout the experiment. Ocean currents at the site were measured using a boat-mounted ADCP as being between 0.16 m/s and 0.30 m/s in a south-westerly direction. There was negligible stratification of temperature or salinity throughout the 25m of water column. Water temperature was 22.7 o C and salinity was 35.0 ppt. On Day 1, all measurements were taken on the two southern offshore jets (SE and SW) from Riser 3 due to the currents observed onsite and the bottom time limitation of the divers at the outfall depth. Only the two outer risers were constructed for the diffusers, while the other two locations were for potential future plant expansion. Figure 6 shows a schematic of the locations. Concentration readings close to the riser were measured by the diver who was able to safely approach the riser structure with the handheld Wet Labs ECO fluorometer. The diver was also able to visually sight individual riser jets issuing from the risers and swim into and around these jets. Concentration measurements further afield from the riser were obtained using the Seabird mounted Chelsea Minitracka fluorometer cast from the roving vessel. Profile plots showing mean concentrations for readings from the diver-held fluorometer from the SE offshore merging jet and SW onshore jet are provided as Figure 7 and Figure 8 respectively. These figures show mean measured concentration values plotted at the measured distance and depth relative to the riser nozzle, which is represented on the bottom left of the figure. All depths were obtained from the calibrated CTD probe attached to the divers BCD unit and cross-checked with WRL staff field notes. Figure 4.1 and Figure 4.3 also present an envelope for the 1 ppt within background salinity contour (22

times dilution) and a contour showing 30 times dilution. These envelopes have been manually interpolated based on the recorded values and the plume behaviour observed during the experiment. Solid lines on the contour envelope represent sections of the contour that were directly interpreted from the measured data. Dashed lines on sections of the contour envelope represent sections of the envelope that were inferred on the basis of nearby concentration measurements and the observations of the divers and the engineers during the experiment. The observed lack of additional dilution between 0 and 40 m for the south west plume can be associated with the likely interaction of plumes issuing from the northern side of the riser and folding back with the currents. Table 2 presents the mean concentrations and dilutions with distance from the outfall for both plumes. Figure 5. (Clockwise from top left) A diver with camera on head and fluorimeter on shoulder; Dye coloured brine jet; the Rhodamine WT injection at the Desalination Plant; and the Roving vessel with shoreline in the background.

Figure 6. Schematic showing the SE measurement line (plume merges with the next riser) and the SW measurement line (independent plume) along with indicative plume footprint and plume lines. Table 2 - Mean Concentration and Dilutions with Distance for SE and SW Plumes SE Line SW Line Distance from Riser 3 (m) Mean Concentration (ppb) Mean Dilution Distance from Riser 3 (m) Mean Concentration (ppb) Mean Dilution 0 2.3 44 0 2.4 43 10 3.4 30 10 2.9 35 30 3.2 32 25 3.3 31 50 2.8 36 45 4.0 26

Figure 7. South East Plume Concentrations of Rhodamine WT. Dilutions can be calculated by dividing by 98. Figure 8. South West Plume Concentrations of Rhodamine WT. Dilutions can be calculated by dividing by 98.

4 CONCLUSIONS Physical modelling and field observations can be compared for the same diffuser configuration, flow and densities. The physical modelling was under quiescent conditions whereas the field experiment was undertaken with currents approximately 0.23 m.s -1. The physical modelling predicted the impact zone to be 16m from the diffuser with a dilution of 27 times. The field observations indicated that the impact point was approximately 30m from the diffuser but with dilutions of approximately 30 times. Indeed, all observed field dilutions including those before and after the impact point were greater than the physical model predicted impact dilution. Clearly the current has a strong influence on the position of the plume, but it is uncertain as to whether it has a significant impact on the dilution. The plume dilutions in the ocean were always slightly greater that the physical modelling indicating that the model was conservative. ACKNOWLEDGEMENT The authors would like to acknowledge Sydney Water, the Bluewater Consortium and Veolia Water. REFERENCES Fischer, H, List, J, Koh, K, Imberger, J and Brooks, N (1979), Mixing in Inland and Coastal Waters. Academic Press, London. Miller, B.M. Smith, G.P. and Tarrade, L. (2009), Physical Modelling of Sydney s Desalination Outfall at 250ML/day Operation with Two Risers, WRL Technical Report No 2009/15. Smith G.P., Miller B.M., Ruprecht J.E. and Flocard F. (2011), Dye Tracer Testing of The Sydney Desalination Plant Outfall: 29th and 30th March 2011, WRL Technical Report No 2009/15. Roberts, J W, Ferrier, A and Daviero, G (1997), Mixing in Inclined Dense Jets. Journal of Hydraulic Engineering. Vol. 123, No. 8 August.