Wind and Wate 235 Wind and waterspray Tate jarrow Colin landon Mike powell U.S. Military Academy West point n Advisor: David Sanders Introduction Given anemometer readings from a nearby building, the task is to devise an algorithm that controls the height of a fountain in an open square. Our mission is to keep passersby dry and yet have the fountain look as impressive as possible. With ever-changing winds, we must devise a scheme to regulate the flow of water through the fountain to ensure that the bulk of the water shot into the air falls back to the ground within the fountain basin boundary. Our model considers many factors and is divided into five basic parts The conversion of wind speed on top of the building to wind speed at ground level based on height and the force of drag The determination of initial velocity, maximum height, and time of flight from fountain nozzle characteristics, using Bernoulli's equation and the rate of flow equation of continui The assessment of the displacement effects of the wind on the water's ascent. The assessment of the displacement effects of the wind on the waters de- scent The calculation of the optimal flow rate by comparing the water's total hor izontal displacement to the radius of the fountain basir After creating this model in a Math CAD worksheet, we solved every func- tion involved in this model as a function of the water flow rate. This worksheet takes the input from several variables such as the nozzle radius, the maximum flow rate the fountain can handle, the dimensions of the building on which the The LIMAP Journal23(3)(2002)235-250. 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 for components of this work owned by others than COMAP must be ho permitted, but copyrights advantage and that copies bear this notice. Abstracting with credit is nored. To copy otherwise, to republish, to post on servers, or to redistribute to lists requires prior permission from COMAP
Wind and Waterspray 235 Wind and Waterspray Tate Jarrow Colin Landon Mike Powell U.S. Military Academy West Point, NY Advisor: David Sanders Introduction Given anemometer readings from a nearby building, the task is to devise an algorithm that controls the height of a fountain in an open square. Our mission is to keep passersby dry and yet have the fountain look as impressive as possible. With ever-changing winds, we must devise a scheme to regulate the flow of water through the fountain to ensure that the bulk of the water shot into the air falls back to the ground within the fountain basin boundary. Our model considers many factors and is divided into five basic parts: • The conversion of wind speed on top of the building to wind speed at ground level based on height and the force of drag. • The determination of initial velocity, maximum height, and time of flight from fountain nozzle characteristics, using Bernoulli’s equation and the rate of flow equation of continuity. • The assessment of the displacement effects of the wind on the water’s ascent. • The assessment of the displacement effects of the wind on the water’s descent. • The calculation of the optimal flow rate by comparing the water’s total horizontal displacement to the radius of the fountain basin. After creating this model in a MathCAD worksheet, we solved every function involved in this model as a function of the water flow rate. This worksheet takes the input from several variables such as the nozzle radius, the maximum flow rate the fountain can handle, the dimensions of the building on which the The UMAP Journal 23 (3) (2002) 235–250. 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
236 The UMAP Journal 23.3 (2002) anemometer is placed, and the dimensions of the fountain From theinputs, the model finds the maximum flow rate that keeps the water in the fountain basin As wind speed and direction vary, the model reacts to produce the optimal flow rate Testing the model shows that while the results are reasonable the main source of error results from our drag calculations due to the interaction between wind and the buildings. To solve this error, measurements should be taken at both the building roof and the fountain itself. Although future work would resolve this issue and improve the model, our current model still provides realistic results We provide in Table 1 a list of symbols used Problem Approach We break the overall problem down into several smaller pieces, solve the pieces separately, and put the pieces together to find the overall solution e How the wind is affected as it flows around the buildin How the wind varies with height off the ground How the buildings slow the wind How the wind affects the water from the fountain How the wind affects the water on the way up How the wind affects the water on the way down How to contain that total displacement within the basin Assumptions Overall assumptions The plaza has a fountain in the center with four surrounding buildings Other arrangements can be handled with slight modifications. The buildings are rectangular and have the same dimensions. Most build- ings are rectangular; for same-size buildings, we can use a single constant drag coefficient The distances from each building to the fountain are the same, so each building has the same effect on the fountain water The acceptable splash area is the radius of the fountain basin. A basin surrounds the water jet, and people walking outside the fountain do not want to get we
236 The UMAP Journal 23.3 (2002) anemometer is placed, and the dimensions of the fountain. From the inputs, the model finds the maximum flow rate that keeps the water in the fountain basin. As wind speed and direction vary, the model reacts to produce the optimal flow rate. Testing the model shows that while the results are reasonable, the main source of error results from our drag calculations due to the interaction between wind and the buildings. To solve this error, measurements should be taken at both the building roof and the fountain itself. Although future work would resolve this issue and improve the model, our current model still provides realistic results. We provide in Table 1 a list of symbols used. Problem Approach We break the overall problem down into several smaller pieces, solve the pieces separately, and put the pieces together to find the overall solution. • How the wind is affected as it flows around the buildings. – How the wind varies with height off the ground. – How the buildings slow the wind. • How the wind affects the water from the fountain. – How the wind affects the water on the way up. – How the wind affects the water on the way down. – How to contain that total displacement within the basin. Assumptions Overall Assumptions • The plaza has a fountain in the center with four surrounding buildings. Other arrangements can be handled with slight modifications. • The buildings are rectangular and have the same dimensions. Most buildings are rectangular; for same-size buildings, we can use a single constant drag coefficient. • The distances from each building to the fountain are the same, so each building has the same effect on the fountain water. • The acceptable splash area is the radius of the fountain basin. A basin surrounds the water jet, and people walking outside the fountain do not want to get wet
Wind and Waterspray Table 1 Table of symbols Symbol Meaning(units rate of flow of the fountain(m/s) Reynolds number flow speed(m/s) kinematic viscosity of the fluid force of drag(N) of the wind(kg/m) Wbh speed of the wind before the building at height h(m/s) Cd drag coefficient A surface area interacting with the wind(m2) vz wind speed measured by the anemometer at the height z(m) h height above ground(variable)(m) z height of the building(m) terrain constant number=0.105 hmax maximum height that the water reaches, a function of R(m) kinetic energy of the wind-building system before the wind hits the building o) Kf kinetic energy of the wind-building system after the wind passes the building 0) work done by nonconservative forces, drag of the building times the leng over which it is applied o) d distance over which drag acts, length and width of the building (m) width or half the length of one of the buildings(m) wh speed of the wind after it passes the building at a height h(m/s) mass of the air that interacts with the building in 1 s if the speed Ubh was constant over the face of the building(kg) angle at which the wind strikes the building(o) Ap cross-sectional area of the pipe at the nozzle tip(m2) f speed of the water as it leaves the nozzle(m/s) radius of the pipe at the nozzle tip(m) 9 acceleration due to gravity, 9.803 m/s rc radius of the column of water at a time t after leaving the nozzle with a rate of flow R(m) pressure on the water caused by the wind(N/s Ac surface area of the column of ascending water(m2) density of air(kg/m) mt total mass of the water in the air at a flow rate R(kg) TTotal total time that the water spends in the air with a flow rate R(s) Water density of water(kg/m horizontal acceleration of the water in the column with a flow rate of r and a wind of speed Uh(m/s2) Fc force on the column of water from the wind of speed vh(N) horizontal displacement of the ascending column of water with a flow rate R and vind speed Uh at a time t(m) ressure on a drop of water from wind of speed uh(N/m2) Fa force on the drop from wind of speed uh(N) area of a drop(m2) md mass of a drop of water(kg) horizontal acceleration of the drop of water as a function of rate of flow R and time aavg average horizontal acceleration of a drop during its descent at a rate of flow R and wind of speed vh(m/s2)
Wind and Waterspray 237 Table 1. Table of symbols. Symbol Meaning (units) R rate of flow of the fountain (m3/s) Re Reynolds number υ flow speed (m/s) d a relevant dimension (m) ν kinematic viscosity of the fluid FD force of drag (N) ρ density of the wind (kg/m3) vbh speed of the wind before the building at height h (m/s) Cd drag coefficient A surface area interacting with the wind (m2) vz wind speed measured by the anemometer at the height z (m) h height above ground (variable) (m) z height of the building (m) α terrain constant number = 0.105 hmax maximum height that the water reaches, a function of R (m) Ki kinetic energy of the wind-building system before the wind hits the building (J) Kf kinetic energy of the wind-building system after the wind passes the building (J) WNC work done by nonconservative forces, drag of the building times the length over which it is applied (J) d distance over which drag acts, length and width of the building (m) b width or half the length of one of the buildings (m) vh speed of the wind after it passes the building at a height h (m/s) m mass of the air that interacts with the building in 1 s if the speed vbh was constant over the face of the building (kg) θ angle at which the wind strikes the building (◦) Ap cross-sectional area of the pipe at the nozzle tip (m2) vf speed of the water as it leaves the nozzle (m/s) rp radius of the pipe at the nozzle tip (m) g acceleration due to gravity, 9.803 m/s2 rc radius of the column of water at a time t after leaving the nozzle with a rate of flow R (m) P pressure on the water caused by the wind (N/s2) Ac surface area of the column of ascending water(m2) ρ density of air (kg/m3) mT total mass of the water in the air at a flow rate R (kg) TTotal total time that the water spends in the air with a flow rate R (s) ρwater density of water (kg/m3) ac horizontal acceleration of the water in the column with a flow rate of R and a wind of speed vh (m/s2) Fc force on the column of water from the wind of speed vh (N) xc horizontal displacement of the ascending column of water with a flow rate R and wind speed vh at a time t (m) PD pressure on a drop of water from wind of speed vh (N/m2) Fd force on the drop from wind of speed vh (N)) Ad area of a drop (m2) md mass of a drop of water (kg) ad horizontal acceleration of the drop of water as a function of rate of flow R and time in air t (m/s2) aavg average horizontal acceleration of a drop during its descent at a rate of flow R and wind of speed vh (m/s2)
238 The UMAP Journal 23.3(2002) Figure 1. The fountain in the center of four buildings The fountain does not squirt water higher than the buildings, although shooting water over the roofs would indeed be spectacular The fountain shoots water straight into the air. This is important for our model so that we can predict how the water will flow up, how it will fall, and where it will fall The fountain nozzle creates a single sustained stream of water. This as- sumption enables us to neglect drag as the water reaches its peak height Furthermore, most fountains have a continuous flow of water Wind The pertinent wind flow is around the sides of the buildings, not over them. Since the fountain does not exceed the height of the buildings, it does not interact with wind that passes over the tops of the buildings. This assumption is important in calculating the drag caused by the building The flow of the wind continues in the same direction across the entire plaza. The wind flows through the plaza in a constant direction, goes around obstacles, and resumes the same direction of motion. The wind does not get stuck in the plaza nor react to cars, people, doors, or windows in the plaza Wakes caused by buildings are not factors. The wake that results when wind hits a building and goes around it does not change the velocity after the wake, so the wake force does not influence the wind s speed or direction The fountain is not in the wake of the buildings. With this assumption, there is no need to worry about wake in our model. This is important because wake is too complex to be modeled The change in wind velocity is due solely to drag. The reason that the wind decreases before and after hitting the building is because of drag. This
238 The UMAP Journal 23.3 (2002) Figure 1. The fountain in the center of four buildings. • The fountain does not squirt water higher than the buildings, although shooting water over the roofs would indeed be spectacular. • The fountain shoots water straight into the air. This is important for our model so that we can predict how the water will flow up, how it will fall, and where it will fall. • The fountain nozzle creates a single sustained stream of water. This assumption enables us to neglect drag as the water reaches its peak height. Furthermore, most fountains have a continuous flow of water. Wind • The pertinent wind flow is around the sides of the buildings, not over them. Since the fountain does not exceed the height of the buildings, it does not interact with wind that passes over the tops of the buildings. This assumption is important in calculating the drag caused by the buildings. • The flow of the wind continues in the same direction across the entire plaza. The wind flows through the plaza in a constant direction, goes around obstacles, and resumes the same direction of motion. The wind does not get stuck in the plaza nor react to cars, people, doors, or windows in the plaza. • Wakes caused by buildings are not factors. The wake that results when wind hits a building and goes around it does not change the velocity after the wake, so the wake force does not influence the wind’s speed or direction. • The fountain is not in the wake of the buildings. With this assumption, there is no need to worry about wake in our model. This is important because wake is too complex to be modeled. • The change in wind velocity is due solely to drag. The reason that the wind decreases before and after hitting the building is because of drag. This
Wind and Waterspray 239 assumption allows us to use the law of conservation of energy to predict the change in velocity The anemometer measures wind speed and direction at the top building before any effects of drag. The anemometer must be at the the building on the windward side, elevated above the height of the building so as not to measure any of the effects of the building. To simplify, we assume that it is at the height of the building The wind pattern is the same across the entire plaza as measured at th anemometer. If the pattern changed, the anemometer reading would be invalid The fountain is in a city or urban area. This assumption allows us to determine the effect of the ground on wind speed a given height The drag applied to wind at a certain height is equal to the average effect of drag, that is, to the total drag caused by the building at the velocity at that height divided by the height of the building. This is slightly inaccurate but still produces a reasonable model Water heig ht Water has laminar flow. Water has a constant velocity at any fixed point, regardless of the time. A fluid may actually have various internal flows that complicate the model, but we consider the flow as the jet of water ascends as constant so that we can model it as an ideal fluid Water has nonviscous flow. The water experiences no viscous drag force in the pipe or in the air. The outer edge of the column of water actually interacts with the air and loses some energy to due to the viscosity of both fluids; but since air and water both have a low viscosity, this loss is negligible Water is incompressible. The density of water is constant and does not change as the water moves up into the air and back down again Water Movement sideways The water jet upward flows as a cylinder. since the surface tension of the water holds it together unless it is acted upon by a force, the water should somewhat retain the dimensions of the nozzle from which it emerges The pressure of the wind is a force per area on the water column and on water drops. Wind and water are both fluids, so the interaction between them is a complex relationship of their viscosities; but we also know that wind creates a pressure difference that we can model. We model the force on the water as the pressure caused by a certain velocity of wind multiplied by the surface area of the body of water
Wind and Waterspray 239 assumption allows us to use the law of conservation of energy to predict the change in velocity. • The anemometer measures wind speed and direction at the top of the building before any effects of drag. The anemometer must be at the top of the building on the windward side, elevated above the height of the building so as not to measure any of the effects of the building. To simplify, we assume that it is at the height of the building. • The wind pattern is the same across the entire plaza as measured at the anemometer. If the pattern changed, the anemometer reading would be invalid • The fountain is in a city or urban area. This assumption allows us to determine the effect of the ground on wind speed a given height. • The drag applied to wind at a certain height is equal to the average effect of drag, that is, to the total drag caused by the building at the velocity at that height divided by the height of the building. This is slightly inaccurate but still produces a reasonable model. Water Height • Water has laminar flow. Water has a constant velocity at any fixed point, regardless of the time. A fluid may actually have various internal flows that complicate the model, but we consider the flow as the jet of water ascends as constant so that we can model it as an ideal fluid. • Water has nonviscous flow. The water experiences no viscous drag force in the pipe or in the air. The outer edge of the column of water actually interacts with the air and loses some energy to due to the viscosity of both fluids; but since air and water both have a low viscosity, this loss is negligible. • Water is incompressible. The density of water is constant and does not change as the water moves up into the air and back down again. Water Movement Sideways • The water jet upward flows as a cylinder. Since the surface tension of the water holds it together unless it is acted upon by a force, the water should somewhat retain the dimensions of the nozzle from which it emerges. • The pressure of the wind is a force per area on the water column and on water drops. Wind and water are both fluids, so the interaction between them is a complex relationship of their viscosities; but we also know that wind creates a pressure difference that we can model. We model the force on the water as the pressure caused by a certain velocity of wind multiplied by the surface area of the body of water