cancel
Showing results for 
Search instead for 
Did you mean: 
cancel
Showing results for 
Search instead for 
Did you mean: 

Community Tip - Visit the PTCooler (the community lounge) to get to know your fellow community members and check out some of Dale's Friday Humor posts! X

Using Mathcad Prime for Snow Farming

DanMarotta
2-Guest

Using Mathcad Prime for Snow Farming

Snow farming might only be a Boston term. It refers to the [re]moving of fallen snow in preparation for an upcoming snow storm.

See how Chris Hartmann used Mathcad Prime to calculate the magnitude of the shoveling challenge posed by a 14-inch snow storm.

http://bit.ly/goP4Jb

How are you using Mathcad to solve your home engineering problems?

26 REPLIES 26
MikeArmstrong
5-Regular Member
(To:DanMarotta)

It seems Chris has too much time on his hands - has he been SNOWED in recently?

Mike

http://communities.ptc.com/servlet/JiveServlet/showImage/2-149416-2532/bar_anim.gif

MikeArmstrong
5-Regular Member
(To:DanMarotta)

Very simple example attached. I was planning on looking at the friction losses and other parameters which effect the flow rate but ran out of time.

Mike

http://communities.ptc.com/servlet/JiveServlet/showImage/2-149416-2532/bar_anim.gif

Interesting pond you have there Mike. Not quite as interesting as your hose though!

MikeArmstrong
5-Regular Member
(To:RichardJ)

Interesting pond you have there Mike. Not quite as interesting as your hose though!

The text in the *.pdf is class, nearly fell off my chair laughing .

Ok very poor example - I was doing my university work at the time, maybe shouldn't have bothered. The pond was already in the garden when I moved in and it keeps several large Koi Carp. The depth is closer to a meter thinking about it and the hose is more like 20mm .

You could have looked at friction losses due to imperfections and length of pipe........ Come on!!!!!!!

Oh, by the way......We only get about 2 days of sunshine a year here, so I wouldn't be filling the pond on that rare occurrence, bur I would have a beer in my hand.

Mike

http://communities.ptc.com/servlet/JiveServlet/showImage/2-149416-2532/bar_anim.gif

You could have looked at friction losses due to imperfections and length of pipe........ Come on!!!!!

That formula does take that into account, but only as long the flow is laminar. Then the friction factor is just 64/Re, where Re is the Reynolds number (i.e. you don't actually need to worry about the size of the imperfections).

MikeArmstrong
5-Regular Member
(To:RichardJ)

That formula does take that into account, but only as long the flow is laminar. Then the friction factor is just 64/Re, where Re is the Reynolds number (i.e. you don't actually need to worry about the size of the imperfections).

Fair enough. Have a look below - an extract from one of our calculations.

Clipboard01.jpg

We take the roughness of the pipe/hose material into consideration in our calculations, which I agree, for a garden hose could be ignored.

Mike

http://communities.ptc.com/servlet/JiveServlet/showImage/2-149416-2532/bar_anim.gif

Well, this is something you know much more about than I do. I got the equations from here:

http://www.efunda.com/formulae/fluids/calc_pipe_friction.cfm

which says that as long as the flow is laminar the friction factor is simply 64/Re. If you simplify their equations using that assumption you get the equations given here:

http://en.wikipedia.org/wiki/Poiseuille%27s_law

So that's what I used.

There is one thing that I don't understand though. For my pond change the water pressure to zero and the elevation change to 1.5m. This represents a real case: when I drain the pond by siphoning over a retaining wall. It says it will take almost 90 hours to drain it, but I know that is not correct because a number of years ago I did it. It took a long time, but not that long! I am going to have to do it again this spring too, because the sludge has built up in the bottom to an unacceptable level

Just a short note.

There usuall isn't laminar floe; it transitions very quickly to turbulent, and Poiseuille flow assumptions fail. That may be the source of your 90 hour estimate. In the attached sheet I get 550 cubic meters drained in 90 hours, a pond one meter deep and 26 meters across.

RichardJ
19-Tanzanite
(To:f.kohlhepp)

Yes, I saw that equation, but I discarded it. The reason is that it's derived from Bernoulli's equation, which does not take into account the friction in the pipe. You have a 1.5m head of water giving a linear velocity of 5 m/s. That is a real jet of water! there is no way the water comes out of the hosepipe that fast. Your siphon would empty my pond in under 15 minutes. It didn't take 90 hours to empty, but it took a lot of hours, not 15 minutes!

Actually, since the viscosity and hose dimensions do not affect your calculations, if I filled my pond with treacle and used a 500ft hose to empty it, I would have a jet of treacle coming out of the end of the hose at 5 m/s

Richard Jackson wrote:

Yes, I saw that equation, but I discarded it. The reason is that it's derived from Bernoulli's equation, which does not take into account the friction in the pipe. You have a 1.5m head of water giving a linear velocity of 5 m/s. That is a real jet of water! there is no way the water comes out of the hosepipe that fast. Your siphon would empty my pond in under 15 minutes. It didn't take 90 hours to empty, but it took a lot of hours, not 15 minutes!

Actually, since the viscosity and hose dimensions do not affect your calculations, if I filled my pond with treacle and used a 500ft hose to empty it, I would have a jet of treacle coming out of the end of the hose at 5 m/s

You're right. (Sorry it took so long to get back; work keeps getting in the way!)

But this is the right equation, we just need to add in the friction term. And the friction term is where length of the hose, diameter of the hose, and viscosity of the treacle get jammed in. The attached file is different (with the same name.)

RichardJ
19-Tanzanite
(To:f.kohlhepp)

You're right. (Sorry it took so long to get back; work keeps getting in the way!)

I know that feeling, although in my case this has been bugging me enough that it has been getting in the way of work!

But this is the right equation, we just need to add in the friction term. And the friction term is where length of the hose, diameter of the hose, and viscosity of the treacle get jammed in. The attached file is different (with the same name.)

Yes, I see. Thanks. But with the friction factor taken into account, If everything is correct we should get the same answers, because all the equations are derived from the Bernoulli equation. The answers are different because both our worksheets have problems! I think I have it figured out now though (famous last words!).

You calculate the velocity, and then calculate the friction factor. Then the you calculate a new velocity. The problem is, the friction factor you calculated is then not correct. You could iterate everything I guess, but it's not necessary (see the attached worksheet). As an aside, why do you have the base length set equal to 1 inch? This effectively sets 1in=1m, and really screws up the units handling. If I evaluate your surface roughness I get 0.85mm, which is a pretty rough hose!

Your worksheet did make me realize what was wrong in mine. I had the viscosity of water too high by a factor of 1000. Oops! Not quite treacle, but not water either! When I fix that the flow is not laminar, so then I needed to fix everything else too, using the correct friction factor. I think the attached worksheet is now correct. The numbers all look more reasonable now. It takes about an hour to drain the pond. It actually took longer than that, but this does not take into account restrictions in the pipe, or the fact that hose inlet kept clogging up! I'm a bit surprised I can fill my pond so quickly, but if I remember then I'm going to time everything this spring to check theory against experiment

Richard Jackson wrote:

Yes, I see. Thanks. But with the friction factor taken into account, If everything is correct we should get the same answers, because all the equations are derived from the Bernoulli equation. The answers are different because both our worksheets have problems!

Thanks for finding the unit length thing. I had no idea that was there and causing such trouble.

You calculate the velocity, and then calculate the friction factor. Then the you calculate a new velocity. The problem is, the friction factor you calculated is then not correct. You could iterate everything I guess, but it's not necessary (see the attached worksheet).

That's why the second solve block--let Mathcad find the velocity that fits the flow and the friction factor.

I've tried once more, using the values from your sheet. (For reasons I can't fathom your sheet will not compute in my version 11, so I can't see your answers.) The draining problem is a differential equation--the driving head decreases as the level in the pond drops. I makefilling about 11 minutes, draining about 200 minutes.)

RichardJ
19-Tanzanite
(To:f.kohlhepp)

Well, I finally found time to look at this again. Work is so inconvenient. It gets in the way of everything!

That's why the second solve block--let Mathcad find the velocity that fits the flow and the friction factor.

I'm not sure quite how I missed what you were doing there

I've tried once more, using the values from your sheet. (For reasons I can't fathom your sheet will not compute in my version 11, so I can't see your answers.)

That is because of the temperature units. I have duplicated the new temperature conversion functions in the worksheet, so it now works in MC11.

If I put in exactly the same values for the viscosity, density and roughness then I can get your worksheet to give the same answer as mine if I make a couple of changes. One change is that the Reynolds number is a function of pipe diameter, so that needs to be carried though into the second solve block (otherwise it's always calculated for 20mm). The second change is that I removed the "1+" from your expression in the second solve block. Then the numbers match mine exactly. I believe mine is correct, for the reasons I mentioned earlier: if I take the simple case of laminar flow my derivation (which does not involve a "1+", just a multiplication) results in the Hagen-Poiseuille equation given here:

http://en.wikipedia.org/wiki/Poiseuille%27s_law

The draining problem is a differential equation--the driving head decreases as the level in the pond drops. I makefilling about 11 minutes, draining about 200 minutes.)

Yes, that is a very good point. The pond depth is significant relative to the siphon height, so this can't just be ignored. Draining is about 65 minutes though, not 200, because you are looking at how long it takes for the height to fall to Delta-z, whereas it only needs to fall to the depth of the pond (-2.5ft). I'm a little surprised this doesn't actually make a bigger difference. To get Odesolve to work after the changes mentioned above I had to cut the solution end point back a bit. I reduced the number of steps too, because I don't have enough patience! In my sheet I did use 1000 steps, because odesolve does not call a function based on nested solve blocks, so it calculates much faster.

I have attached two worksheets: a revised version of your sheet, and a revised version of mine.

Richard Jackson wrote:

Well, I finally found time to look at this again. Work is so inconvenient. It gets in the way of everything!

If I put in exactly the same values for the viscosity, density and roughness then I can get your worksheet to give the same answer as mine if I make a couple of changes. One change is that the Reynolds number is a function of pipe diameter, so that needs to be carried though into the second solve block (otherwise it's always calculated for 20mm). The second change is that I removed the "1+" from your expression in the second solve block. Then the numbers match mine exactly. I believe mine is correct, for the reasons I mentioned earlier: if I take the simple case of laminar flow my derivation (which does not involve a "1+", just a multiplication) results in the Hagen-Poiseuille equation given here:

This horse is gasping!

I've attached a one page discussion of why the "1+" is there; simple algebra and the Bernoulli equation. The difference in this case is small, we're spending most of the pressure on friction, very little on accelerating the water. I've reposted the molested corpse one more time to illustrate--23.36 gpm instead of 24.36 gpm on fill requires nearly a half minute more.

This horse is gasping!

That means it's not quite dead yet. Get up horse!

http://communities.ptc.com/servlet/JiveServlet/showImage/38-1286-6326/DeadHorse.gif

Your logic seems sound to me, but then there's something I don't understand. See the attached.

I notice that the derivation of the Hagen-Poiseuille equation on Wikipedia involves a Taylor expansion, and then higher order terms are dropped. If the original equation on the website is not correct, which your derivation certainly indicates, perhaps that's why my derivation can get from that to the Hagen-Poiseuille equation. Seems like quite a coincidence, but it might be the reason.

Poor horse!

I think we're solving the same problem for two different things.

Thanks

I think I finally have it straight. My equation is correct for the filling of the pond, because the water in the house system is not stagnant. The velocity at the inlet and outlet are the same, so the dynamic pressure terms cancel. Your equation is correct when draining the pond, because the water in the pond is stagnant (in my pond, both in the fluid dynamics sense and in the common sense of the word!). Therefore there is no dynamic pressure at the inlet. That is the difference between the derivation of the Hagen-Poiseuille equation and the siphon equation.

I am hoping the horse can have a well deserved rest at this point!

Basically.

If I calculate the dynamic pressure (rho * V_fill^2 /2) at the spigot on the side of your house, I get about 4.3 psi. This is the difference between the pressure in the pipe with the spigot closed (the stagnation pressure) and the "static" pressure in the hose at the spigot with the spigot open and filling the pond.

Both the water in the house and the water in the pond have a "stagnation" pressure--that's the total amount of energy that can be used to push the water through the hose, some going to friction losses and some being converted to kinetic energy in the flowing water (AKA dynamic pressure.)

Total ("stagnation" or "pitot") pressure is always the sum of the static pressure and the dynamic pressure. Friction losses (extracting non-recoverable energy from the flow) subtract from the total pressure (and, therefore, from the static pressure too.) So we're both right.

pressures.jpg

So we're both right

No. There is one physical problem, therefore only one correct (which can allow for stated approximations or assumptions) answer. So we could both be wrong, or one of us could be right (not me, I think), but we cannot both be right. Water does not flow out of the end of a hose at two different velocities simultaneously.

I think I am getting my static and total (stagnation) pressures mixed up, in particular in the Bernoulli equation. If the house water pressure is defined as the stagnation pressure when the spigot is off, not the static pressure when it is on (which seems reasonable), then I get your version of the equation.

But I can then only get the siphon to work if I assume that the stagnation pressure of the pond is zero.

But I can then only get the siphon to work if I assume that the stagnation pressure of the pond is zero.

No, the stagnation pressure of the pond (measured against the end of the hose) is rho x gravity x delta Z, the difference in height between the top of the pond and the discharge end of the hose. That is the total energy available to push water through the hose. Some of that gets dissipated in friction in the hose, the rest gets turned into water velocity at the end of the hose. The static pressure at the hose discharge is equal to the atmospheric pressure at the end of the hose, since that's our reference prressure, static pressure there is 0 psi gage. Total pressure at the top of the pond is rho x g x Z, so that minus the friction factor losses must equal rho x V^2 at the end of the hose. since friction factor losses are f x (rho x V^2), the equation is:

rho x g x Z = rho x V^2 x(1+f).

So the stagnation pressure includes the pressure due to height? Then the house water pressure is not the stagnation pressure. Or the static pressure. It's something else, for which I don't have a term, that is equal to the stagnation pressure minus the pressure due to height.

http://communities.ptc.com/servlet/JiveServlet/showImage/38-1286-7413/confused.gif

Richard Jackson wrote:

So the stagnation pressure includes the pressure due to height? Then the house water pressure is not the stagnation pressure. Or the static pressure. It's something else, for which I don't have a term, that is equal to the stagnation pressure minus the pressure due to height.

http://communities.ptc.com/servlet/JiveServlet/showImage/38-1286-7413/confused.gif

Let us measure ALL pressures relative to the atmospheric pressure at the discharge end of the hose (P_atm).

Then:

The static pressure at the end of the hose is zero. (The pressure water at the end of the hose is the same pressure as the air it's discharging into.

P.static_end = P_atm, but we're measuring relative to that, so P.static_end - P_atm = 0 psig.

The dynamic pressure in the hose is rho/2 x V^2, and it's constant anywhere in the hose. (The hose is constant area, so the water coming out is at the same speed as the water going in and the density is the same (water is virtualy incompressible at these pressures.)

So the stagnation pressure at the end of the hose is: P_o_end = P.static_end + rho/2 x V^2 psi gage (= P.static_end + rho/2 x V^2 +P_atm absolute.) But we said above that P.static end = 0, so P_o_end = rho/2 x V^2 psi gage.

At the spigot (at the house), the dynamic pressure is the same. The stagnation pressure is:

P_o_end +rho x g x Z +rho/2 x V^2 x friction+P_static_in===the stagnation pressure at the discharge plus what was lost to friction plus what we gain from the elevation change plus the static pressure in the hose. Substituting for P_o_end, we get P_o_in = rho x g x Z +(P_o_end = rho/2 x V^2) +rho/2 x V^2 x friction + P_static_in

or P_o_in = rho x g x Z + rho/2 x V^2 x (1+friction) + P_static_in

OK, I guess it's a question of the reference point for the stagnation pressure. For the siphon, I have the pressure due to the height difference as an explicit term, so I should have said the stagnation pressure of the pond relative to the hose inlet is zero". Similarly, if the house water pressure is taken to be the stagnation pressure, that is the stagnation pressure relative to the hose inlet.

Mike Armstrong wrote:

We take the roughness of the pipe/hose material into consideration in our calculations,

Not surprising given that the Reynolds number in your diagram is over 4000000. The transition from laminar to turbulent flow occurs around 2000 to 3000 (see the Moody chart at http://en.wikipedia.org/wiki/File:Moody_diagram.jpg). Surface roughness is only significant in turbulent flow.

Alan

AlanStevens wrote:

Mike Armstrong wrote:

We take the roughness of the pipe/hose material into consideration in our calculations,

Not surprising given that the Reynolds number in your diagram is over 4000000. The transition from laminar to turbulent flow occurs around 2000 to 3000 (see the Moody chart at http://en.wikipedia.org/wiki/File:Moody_diagram.jpg). Surface roughness is only significant in turbulent flow.

Alan

One type of Moody chart - Nikuradse chart on MAS (WebSheet):

http://twt.mpei.ac.ru/MCS/Worksheets/Hydro/La-De-Re.xmcd

MikeArmstrong
5-Regular Member
(To:AlanStevens)

Nicely pointed out.

Mike

http://communities.ptc.com/servlet/JiveServlet/showImage/2-149416-2532/bar_anim.gif

Top Tags