A Foul-Weather Fountain 251 a Foul-Weather fountain lyan K Card R Ernie e esser Jeffrey H Giansiracusa University of Washin Seattle Wa Advisor: James Allen Morrow Introduction We devise a fountain control algorithm to monitor wind conditions and ensure that a fountain at the center of a plaza fires water high enough to be dazzling while not drenching the pedestrian areas surrounding the fountain We construct a model of a fountain based on the physics of falling water under various wind conditions through computer simulation. Using complex analytic techniques, we model the wind flow through the plaza and estimate how anemometer readings from a nearby rooftop relate to plaza conditions We construct four algorithms-two intelligent algorithms, a conservatiy approach, and an enthusiastic system-to control the fountain We devise a measure of unacceptable spray levels outside the fountain and use this criterion to compare performance. First, we examine the behavior of these algorithms under general abstract wind conditions. Then we construct a wind signal generator that simulates the conditions of several major cities from meteorological database data, and we compare the performance of our control systems in each city Simulations show that the Conservative and enthusiastic algorithms both perform unacceptably in realistic conditions. The Weighted Average Algorithm works best in gusty cities such as Chicago, but the Averaging Algorithm is superior in calmer cities such as Los Angeles and Seattle The control algorithm cannot possibly respond to changes in conditions at anything below the 10 s scale, since wind is highly variable and the response of the anemometer is somewhat slow [Industrial Weather Products 2002]. The goal is therefore to design the algorithm to operate on a time scale of 10 s up to a couple of hours and adapt the height of the fountain to a maximum safe level The UMAP Journal 23(3)(2002)251-266. Copyright 2002 by COMAP, Inc. Allrights reserved Permission to make digital or hard copies of part or all of this work for personal or classroom use is granted without fee provided that copies are not made or distributed for profit or commercial dvantage and that copies bear this notice. Abstracting with credit is permitted, but copyrights for components of this work owned by others than COMAP must be honored. To copy otherwise, to republish, to post on servers, or to redistribute to lists requires prior permission from COMAP
A Foul-Weather Fountain 251 A Foul-Weather Fountain Ryan K. Card Ernie E. Esser Jeffrey H. Giansiracusa University of Washington Seattle, WA Advisor: James Allen Morrow Introduction We devise a fountain control algorithm to monitor wind conditions and ensure that a fountain at the center of a plaza fires water high enough to be dazzling while not drenching the pedestrian areas surrounding the fountain. We construct a model of a fountain based on the physics of falling water droplets considered as a particle system. We examine the behavior of a fountain under various wind conditions through computer simulation. Using complex analytic techniques, we model the wind flow through the plaza and estimate how anemometer readings from a nearby rooftop relate to plaza conditions. We construct four algorithms—two intelligent algorithms, a conservative approach, and an enthusiastic system—to control the fountain. We devise a measure of unacceptable spray levels outside the fountain and use this criterion to compare performance. First, we examine the behavior of these algorithms under general abstract wind conditions. Then we construct a wind signal generator that simulates the conditions of several major cities from meteorological database data, and we compare the performance of our control systems in each city. Simulations show that the Conservative and Enthusiastic algorithms both perform unacceptably in realistic conditions. The Weighted Average Algorithm works best in gusty cities such as Chicago, but the Averaging Algorithm is superior in calmer cities such as Los Angeles and Seattle. The control algorithm cannot possibly respond to changes in conditions at anything below the 10 s scale, since wind is highly variable and the response of the anemometer is somewhat slow [Industrial Weather Products 2002]. The goal is therefore to design the algorithm to operate on a time scale of 10 s up to a couple of hours and adapt the height of the fountain to a maximum safe level. The UMAP Journal 23 (3) (2002) 251–266. c Copyright 2002 by COMAP, Inc. All rights reserved. Permission to make digital or hard copies of part or all of this work for personal or classroom use is granted without fee provided that copies are not made or distributed for profit or commercial advantage and that copies bear this notice. Abstracting with credit is permitted, but copyrights for components of this work owned by others than COMAP must be honored. To copy otherwise, to republish, to post on servers, or to redistribute to lists requires prior permission from COMAP
252 The UMAP Journal 23.3 (2002) Model of the water jet Ne model the spray from the fountain as a particle system. As water droplets spew forth from the nozzle, they are subjected to forces(gravity, air drag, turbulence, etc. ) We formulate a simplified differential equation governing the motion and then numerically integrate to find the trajectory for each droplet This equation is based on a physically realistic model of small droplets(around 1 mm radius)and we scale it up to an effective model for larger clumps of water(up to 10 cm across) because the physics of turbulence and viscosity at he larger scale cannot be computed accurately We need the following assumptions The drag force is proportional to the square of the speed and to the square f the radius [ Nasa 2002 Droplets break into smaller droplets when subjected to wind. Breakup rate is proportional to relative wind speed and surface area [INobauer 1999 When a droplet breaks turbulence causes the new droplet fragments to move slightly away from their initial trajectory. Modeling a Single Droplet We formulate the motion of a water droplet a m dt=-mgi+mlwul-ir where vis the velocity w is the wind velocity relative to the motion of the droplet (wind vector minus velocity vector), m and r are the droplet's mass and radius, and n is a constant of proportionality. According to the Virtual Science Center Project Team[2002 a raindrop with radius 1 mm falls at a terminal velocity of 7 m/s; so we determine that n=0.855 kg/m. Large drops fall quickly; very tiny drops fall very slowly, mimicking a fine mist that hangs in the air for a long time We assume droplet breakup is a modified Poisson process, with rate λoulr If the breakup rate did not depend on variable parameters w and r2, the process would be a standard Poisson process. We l determine Ao by fitting the water ream of our fountain to the streams of two real fountains: the jet D'Eau of Geneva, Switzerland, and the Five Rivers Fountain of Lights in Miami, Florida When a breakup occurs, we split the droplet into two new droplets and divide the mass randomly, using a uniform distribution. Air turbulence tends to impart to the two new droplets a small velocity component perpendicular to the relative wind direction w. This effect causes a tight stream of water to spread
252 The UMAP Journal 23.3 (2002) Model of the Water Jet We model the spray from the fountain as a particle system. As water droplets spew forth from the nozzle, they are subjected to forces (gravity, air drag, turbulence, etc.). We formulate a simplified differential equation governing the motion and then numerically integrate to find the trajectory for each droplet. This equation is based on a physically realistic model of small droplets (around 1 mm radius) and we scale it up to an effective model for larger clumps of water (up to 10 cm across) because the physics of turbulence and viscosity at the larger scale cannot be computed accurately. We need the following assumptions: • The drag force is proportional to the square of the speed and to the square of the radius [Nasa 2002]. • Droplets break into smaller droplets when subjected to wind. Breakup rate is proportional to relative wind speed and surface area [Nobauer 1999]. • When a droplet breaks, turbulence causes the new droplet fragments to move slightly away from their initial trajectory. Modeling a Single Droplet We formulate the motion of a water droplet as m dv dt = −mgzˆ + η|w| 2wrˆ 2, wherev is the velocity, w is the wind velocity relative to the motion of the droplet (wind vector minus velocity vector), m and r are the droplet’s mass and radius, and η is a constant of proportionality. According to the Virtual Science Center Project Team [2002], a raindrop with radius 1 mm falls at a terminal velocity of 7 m/s; so we determine that η = 0.855 kg/m3. Large drops fall quickly; very tiny drops fall very slowly, mimicking a fine mist that hangs in the air for a long time. We assume droplet breakup is a modified Poisson process, with rate λbreakup = λ0|w|r2. If the breakup rate did not depend on variable parameters |w| andr2, the process would be a standard Poisson process. We l determine λ0 by fitting the water stream of our fountain to the streams of two real fountains: the Jet D’Eau of Geneva, Switzerland, and the Five Rivers Fountain of Lights in Miami, Florida. When a breakup occurs, we split the droplet into two new droplets and divide the mass randomly, using a uniform distribution. Air turbulence tends to impart to the two new droplets a small velocity component perpendicular to the relative wind direction w. This effect causes a tight stream of water to spread
A Foul-Weather Fountain 253 out as it travels, even under zero-wind conditions. We let this velocity nudge have magnitude 2% of the particle' s speed relative to the air and a random perpendicular direction. We give the two drops equal and opposite nudges Putting water drops Together to make a Fountain We define the waterjet as a stream of large water drops. Their size is roughly the size of the nozzle and they leave with an initial velocity equal to the nozzle output velocity( Figure 1) oO A real four Representati。n。 fa fountain Figure 1. A continuous water jet is approximated by a discrete stream of water blobs The water blobs leave at a rate such that the flux of water is equal to the flux given by a nozzle-sized cylindrical stream moving at the same speed To model the turbulence in the jet as the water leaves the nozzle, we give each water blob a normal distribution of radius and initial speed . The standard deviation of blob radii is 10% of the nozzle size The standard deviation of initial speeds is 5% of the initial speed The blobs leave with an angular spread of 3, consistent with industrial high-pressure nozzles [Spray Nozzles 2002] Wind drag in particle streams is significantly reduced for particles follow ing one another closely (nASCar drivers and racing cyclists are intima familiar with this phenomenon). These effects are already incorporated into the dynamics of large water blobs(which can be thought of as representing many small drops moving together
A Foul-Weather Fountain 253 out as it travels, even under zero-wind conditions. We let this velocity nudge have magnitude 2% of the particle’s speed relative to the air and a random perpendicular direction. We give the two drops equal and opposite nudges. Putting Water Drops Together to Make a Fountain We define the water jet as a stream of large water drops. Their size is roughly the size of the nozzle, and they leave with an initial velocity equal to the nozzle’s output velocity (Figure 1). Figure 1. A continuous water jet is approximated by a discrete stream of water blobs. The water blobs leave at a rate such that the flux of water is equal to the flux given by a nozzle-sized cylindrical stream moving at the same speed. To model the turbulence in the jet as the water leaves the nozzle, we give each water blob a normal distribution of radius and initial speed: • The standard deviation of blob radii is 10% of the nozzle size. • The standard deviation of initial speeds is 5% of the initial speed. • The blobs leave with an angular spread of 3◦, consistent with industrial high-pressure nozzles [Spray Nozzles 2002]. Wind drag in particle streams is significantly reduced for particles following one another closely (NASCAR drivers and racing cyclists are intimately familiar with this phenomenon). These effects are already incorporated into the dynamics of large water blobs (which can be thought of as representing many small drops moving together)
254 The UMAP Journal 23.3(2002) Fitting the Fountain The Five Rivers Fountain of Lights in Daytona, Florida, is one of the largest fountains in the world. It consists of several water jets, and on low-wind days each propels a water stream 60 m high and 120 m out. The jet D Eau in geneva Switzerland, anoth pressive fountain, shoots a 30 cm-diameter stream of water at 60 m/s straight up. The water reaches a height of 140 m and on an average breezy day (wind speed 5 m/s)returns to earth approximately 35 m downwind from the nozzle [micheloud cie 2002(Figure 2 45d-西 30 m Figure 2. The jet D Eau and the Fountain of lights To determine Xo, we first match our geometry to the Five rivers fountain of Lights. We fix Ao so that with an initial velocity such that the stream reaches height of 60 m, it returns to the ground at a distance of just over 100 m. Too large a Ao results in the water breaking up too quickly into tiny droplets, which have a much lower terminal velocity and thus fail to reach the desired distance; if the value is too small, then an unrealistically small amount of spray is produced and the water blob travels too far. The results are summarized in Table 1 We set Ao 5000. The results are highly insensitive to this parameter varying o by a factor of 2 cause only a 15% changes in the distances. Therefore, even though our method for determining this parameter is fairly rough, the important behavior is much more strongly affected by other parameters Table 1 Comparison between real fountains and our model Jet DEa eal model model Height(m) 140 121 Distance() 3 120 100 We conclude from this comparison that our model reproduces the spray patterns of extreme fountains accurate to within about 15%%. We expect that for
254 The UMAP Journal 23.3 (2002) Fitting the Fountain The Five Rivers Fountain of Lights in Daytona, Florida, is one of the largest fountains in the world. It consists of several water jets, and on low-wind days each propels a water stream 60 m high and 120 m out. The Jet D’Eau in Geneva, Switzerland, another impressive fountain, shoots a 30 cm-diameter stream of water at 60 m/s straight up. The water reaches a height of 140 m and on an average breezy day (wind speed 5 m/s) returns to earth approximately 35 m downwind from the nozzle [Micheloud & Cie 2002] (Figure 2). Figure 2. The Jet D’Eau and the Fountain of Lights. To determine λ0, we first match our geometry to the Five Rivers Fountain of Lights. We fix λ0 so that with an initial velocity such that the stream reaches a height of 60 m, it returns to the ground at a distance of just over 100 m. Too large a λ0 results in the water breaking up too quickly into tiny droplets, which have a much lower terminal velocity and thus fail to reach the desired distance; if the value is too small, then an unrealistically small amount of spray is produced and the water blob travels too far. The results are summarized in Table 1. We set λ0 = 5000. The results are highly insensitive to this parameter; varying λ0 by a factor of 2 cause only a 15% changes in the distances. Therefore, even though our method for determining this parameter is fairly rough, the important behavior is much more strongly affected by other parameters. Table 1. Comparison between real fountains and our model. Jet D’Eau Five Rivers Fountain real model real model Height (m) 140 121 60 62 Distance (m) 35 30 120 100 We conclude from this comparison that our model reproduces the spray patterns of extreme fountains accurate to within about 15%. We expect that for
A Foul-Weather Fountain 255 a plaza-sized fountain, our model will be more accurate, since our formulas for breakup and drag force are derived under less extreme conditions Wind Flow Through the Plaza Buildings and other structures in an urban environment can cause signif- icant disturbances to wind flow patterns; rooftop and street-level conditions can often be quite different. To model the plaza wind, we assume There are no significant structures between the buildings beside of the plaza The plaza is large, so effects caused when wind flow leaves the plaza negligible at the plaza center; the significant effects are entirely caused at inward boundary passage The air flow is smooth enough so that turbulent vortices are negligible Formulation We approximate the geometry of the plaza as in Figure 3 and use complex analytic flow techniques[ Fisher 1990, 225 Plaza Figure 3. Schematic representation of the relevant features of the plaza With a Schwarz-Cristoffel mapping of a smooth horizontal flow from the upper half of the complex plane onto the region above the plaza, we obtain a flow function for the wind as it enters the plaza area r0={(+(2-1y12+kg(++(+i92-1y/)} where t parametrizes a streamline for each value of c. These streamlines are plotted in Figure 4, where the acceleration of the wind as it passes over the building edge and the decreased velocity in the plaza are both clearly visible The flow velocity v is inversely proportional to the streamline spacing, so the horizontal component of it is O dc
A Foul-Weather Fountain 255 a plaza-sized fountain, our model will be more accurate, since our formulas for breakup and drag force are derived under less extreme conditions. Wind Flow Through the Plaza Buildings and other structures in an urban environment can cause significant disturbances to wind flow patterns; rooftop and street-level conditions can often be quite different. To model the plaza wind, we assume: • There are no significant structures between the buildings beside of the plaza. • The plaza is large, so effects caused when wind flow leaves the plaza are negligible at the plaza center; the significant effects are entirely caused at the inward boundary passage. • The air flow is smooth enough so that turbulent vortices are negligible. Formulation We approximate the geometry of the plaza as in Figure 3 and use complex analytic flow techniques [Fisher 1990, 225]. Figure 3. Schematic representation of the relevant features of the plaza. With a Schwarz-Cristoffel mapping of a smooth horizontal flow from the upper half of the complex plane onto the region above the plaza, we obtain a flow function for the wind as it enters the plaza area: Γc(t) = h0 π (t + ic) 2 − 1 1/2 + log t + ic + (t + ic) 2 − 1 1/2 , where t parametrizes a streamline for each value of c. These streamlines are plotted in Figure 4, where the acceleration of the wind as it passes over the building edge and the decreased velocity in the plaza are both clearly visible. The flow velocity v is inversely proportional to the streamline spacing, so the horizontal component of it is vx = Im ∂Γc ∂c .