Pages:
1
..
8
9
10
11
12
..
19 |
aga
Forum Drunkard
   
Posts: 7030
Registered: 25-3-2014
Member Is Offline
|
|
Quote: Originally posted by blogfast25  | Did you take into account the surface area of the top and bottom?
I'm curious to see your data and how you will derive your estimate of h from them. I know how I would. |
Erm, i added the area of the middle bit to the last bit but not the other bit, so it could be a way out.
It wasn't clear whether that formula for 'area of a cone' was the entire cone. I assumed it was.
Raw data:-
Secs deg C
0 62
60 60
120 57
180 56
240 55
300 54
360 53.5
420 52
480 50.5
540 49
600 48
I just set the equation up in a spreadsheet, used the calculated values (also the s.h.c. of water @ 4.186) and tried 'h' numbers until the result
roughly agreed with the experimental result.
Quicker for me that way.
|
|
blogfast25
International Hazard
   
Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline
Mood: No Mood
|
|
Quote: Originally posted by aga  |
I just set the equation up in a spreadsheet, used the calculated values (also the s.h.c. of water @ 4.186) and tried 'h' numbers until the result
roughly agreed with the experimental result.
Quicker for me that way. |
The scientific way of doing this would be as follows.
Remember that:
ln[T−TaTi−Ta]=−τt
You would plot in Excell:
ln[T−TaTi−Ta]
versus:
t
That should give something looking like:

The + are the data points, they should roughly be on a straight line. Then make Excell fit a line to the data (linear regression).
The slope of the line is:
−τ
Then calculate h from tau.
|
|
aga
Forum Drunkard
   
Posts: 7030
Registered: 25-3-2014
Member Is Offline
|
|
All day reviving algebra cells did not leave any brain capacity to learn what a linear regression was.
What does a second or third order integrally equation look like ?
|
|
blogfast25
International Hazard
   
Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline
Mood: No Mood
|
|
Checking surface area of the cup, by calculus.
For a revolution body:

... the surface area A (end bits not counted) is given by:
A=2π∫baf(x)√1+[f′(x)]2dx
For the cup, with an x-axis running through it horizontal centre axis:
f(x)=517x+70
f′(x)=517
A=2π∫850√1+52172(517x+70)dx
A=2.085π[534x2+70x]850=45927
Lid:
π4952=7088
Cone + lid:
53015mm2=0.053m2
Quote: Originally posted by aga  | All day reviving algebra cells did not leave any brain capacity to learn what a linear regression was.
What does a second or third order integrally equation look like ? |
If the first derivative is f'(x):
f′(x)=df(x)dx
Then the second derivative is f''(x):
f″(x)=df′(x)dx=f2(x)dx2
y‴=d(y″)dx=d3ydx3
Example:
md2ydt2=mg−ky
or:
my″=mg−ky
Newtonian notation ('overdot notation'):
m¨y=mg−ky
... is the DE of motion of a bob of mass m hanging off a spring with spring constant k.
[Edited on 18-4-2016 by blogfast25]
|
|
aga
Forum Drunkard
   
Posts: 7030
Registered: 25-3-2014
Member Is Offline
|
|
Is 'g' gravitational acceleration ?
Must be, seeing as it popped out of nowhere 
You do realise that i'll soon be able to calculate proof for agaspace ...
Well, maybe if i can finally master those nasty fractions, do 'u' substitutions, drink less, think straight etc.
|
|
blogfast25
International Hazard
   
Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline
Mood: No Mood
|
|
Yes. The bob is subject to it's own weight:
mg
|
|
aga
Forum Drunkard
   
Posts: 7030
Registered: 25-3-2014
Member Is Offline
|
|
Of course, same as everything in agaspace, including Time : all attributes are inseparable.
g′=Σ∫allo′t′otherbits
Still regard these higher equations as 'best fit' or 'best we have found so far' mostly due to things like Limit Theorem : it's all basically an
approximation, although a damned good one.
[Edited on 17-4-2016 by aga]
|
|
blogfast25
International Hazard
   
Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline
Mood: No Mood
|
|
@aga:
I happen to know that for agaspace you need partial derivatives:
For a function f of x AND y:
f(x,y)
Partial derivative to x:
∂f∂x
... and to y:
∂f∂y
You can't calculate agaspace astral mass time differential tensors w/o them! 
But we don't need them for agasyphonTM, so we'll stick to the latter, tomorrow.
[Edited on 17-4-2016 by blogfast25]
|
|
blogfast25
International Hazard
   
Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline
Mood: No Mood
|
|
Here's the plot for:
ln[T−TaTi−Ta]=−τt
[T−TaTi−Ta]
... is also called the "Reduced Temperature", a dimensionless number between 1 and 0.
I used Ta=20 C.

The linear regression equation:
y=−0.0006x−0.021
... gives us an estimate for:
τ=0.0006s−1
τ=hAmcp
With A=0.053 m2, m=0.458 kg, cp=4181 J/kg K then:
h=τmcpA=22Wm−2K−1
Considering there's some noise on the data, this is a crude estimate, of course.
This source puts up a value range of 5 to 37 Wm-2K-1 for free convection in gas (air) or dry vapour, so we're in the right ball
park.
[Edited on 18-4-2016 by blogfast25]
|
|
aga
Forum Drunkard
   
Posts: 7030
Registered: 25-3-2014
Member Is Offline
|
|
22 ?
Doh. I was very uncertain of the Units, so perhaps i got metres/millimetres mixed up somewhere.
f(x,y,z)
now that would definitely come in handy.
Quote: | astral mass time differential tensors |
They sound very cool !
What colour are they ?
|
|
blogfast25
International Hazard
   
Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline
Mood: No Mood
|
|
Quote: Originally posted by aga  | 22 ?
Doh. I was very uncertain of the Units, so perhaps i got metres/millimetres mixed up somewhere.
f(x,y,z)
now that would definitely come in handy.
Quote: | astral mass time differential tensors |
They sound very cool !
What colour are they ? |
Perhaps the most common f(x,y,z) type function is the volume of a right angled box:
V=lhw
They're appleblueseagreen.
In a couple of hours: calculus of the agasyphonTM. Yeehaw!
|
|
aga
Forum Drunkard
   
Posts: 7030
Registered: 25-3-2014
Member Is Offline
|
|
Volume ! Of course !
So if the function f(x,y,z) is simply xyz, then f'(x,y,z) is one of the missing links with t ! (and E, maybe g as well)
|
|
blogfast25
International Hazard
   
Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline
Mood: No Mood
|
|
Syphon: tank emptying time

The tank's shape is of no consequence as long as its cross-section A is constant.
The outflow speed v2 can be determined by applying Bernoulli's Principle between points 1 and 2:
12v22+gy2+p2ρ=12v21+gy1+p1ρ
The derivation is a bit lengthy, so I've posted it below the second page break.
The result is:
v2=c√2gy
So the volumetric throughput Qv is:
Qv=π4D2v2=π4D2c√2gy
During an infinitesimal amount of time dt, the water level decreases by dy (<0), so:
dV=−Ady=Qvdt
−Ady=π4D2c√2hydt=π4D2c√2g√ydt
Let's call:
α=πD2c√2g4A
−y−1/2dy=αdt
Integrate between relevant boundaries:
−∫h0Hy−1/2dy=∫t0αdt
2[√y]Hh0=αt
t=2α(√H−√h0)
Note that the higher up the tank is placed with respect to 0, the smaller t becomes.
<hr>
But. There's a big but.
The model above doesn't take into account any viscous losses in the syphon's pipe.
To take them into account, with Bernoulli, we get:
12c2v22=gy−Δpsρ
For laminar flow (turbulent flow makes it more complicated) through the pipe:
Δps=32μLv2D2
Where:
L≈2(H−h0)+h0
Substituting we get:
12c2v22=gy−32μLv2ρD2
12v22+βv2−c2gy=0
With:
β=32μLc2ρD2
This quadratic equation has one positive root:
v2=−β+√β2+2c2gy
So:
Qv=π4D2v2=π4D2(√β2+2c2gy−β)
Inserting into the DV above the break:
−∫h0Hdy√β2+2c2gy−β=∫t0αdt
Where:
α=πD2c4A
Although that integral is computable, the result gets very bulky and messy, so I won't go there.
<hr>
Full derivation:
12v22+gy1+p2ρ=12v21+gy2+p1ρ
12(v22−v21)+p0ρ=gh0+p1ρ
A2v2=Av1
⟹v22−v21=v22[1−(A2A)2]
A2=π4D2
1c2=[1−(A2A)2]
12c2v22+p0ρ=gh0+p1ρ
p1=ρg(y−h0)+p0
12c2v22=gy
v2=c√2gy
Qv=π4D2v2=π4D2c√2gy
dV=−Ady=Qvdt
−Ady=π4D2c√2hydt=π4D2c√2g√ydt
α=πD2c√2g4A
−y−1/2dy=αdt
−∫h0Hy−1/2dy=∫t0αdt
2[√y]Hh0=αt
t=2α(√H−√h0)
But.
12c2v22=gy−Δpsρ
Δps=32μLv2D2
L≈2(H−h0)+h0
12c2v22=gy−32μLv2ρD2
12v22+βv2−c2gy=0
β=32μLc2ρD2
v2=−β+√β2+2c2gy
Qv=π4D2v2=π4D2(√β2+2c2gy−β)
−∫h0Hdy√β2+2c2gy−β=∫t0αdt
α=πD2c4A
[Edited on 19-4-2016 by blogfast25]
|
|
aga
Forum Drunkard
   
Posts: 7030
Registered: 25-3-2014
Member Is Offline
|
|
Cor Blimey !
I think a bit of commentary is going to be required to follow that properly !
|
|
blogfast25
International Hazard
   
Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline
Mood: No Mood
|
|
Quote: Originally posted by aga  | Cor Blimey !
I think a bit of commentary is going to be required to follow that properly ! |
Ok. Briefly.
Often the preamble necessary to set up the DV is longer and harder work than actually solving the DV: that was definitely the case here.
The actual DV:
dV=−Ady=Qvdt
... is obtained by looking at an emptying tank during an infinitesimally small t+dt time interval. In that time, a volume:
Qvdt
Is discharged through the syphon. That causes the level in the tank to drop by -dy or a volume dV:
dV=−Ady
So:
−Ady=Qvdt
Then it's a matter of plugging in Qv, rearranging and integrating to get the tank's emptying time t.
Now look at:
Qv(y)=π4D2c√2gy
This implies something anyone with a passing acquaintance of syphons will know: the lower you place the outlet (point 2), the higher the
throughput!
This is also reflected by:
t=2α(√H−√h0)
Take values of H = 9 and h0 = 4 (tank height is 5).
That gives:
t=2α
Now take values of H = 16 and h0 = 11 (tank height is also 5).
That gives:
t=1.37α
32 % shorter!
Now for the case where we try and take pipe viscous losses into account:
Qv=π4D2(√β2+2c2gy−β)
β=32μLc2ρD2
If you look closely at that function, it appears that L on the one hand increases Qv, on the other hand it seems to decrease it.
Could there be an optimum L?
We can find this by Optimum Theory, by evaluating:
dQvdL=0
I did this and found no optimum. That's good news because it means that as in the simpler 'no pipe friction' case, the lower the syphon outlet is, the
higher the output and the smaller the emptying time (there exists no optimum L beyond which that statement is no longer true).
The formula:
t=2α(√H−√h0)
... can be used for a lower limit estimate of the emptying time, bearing in mind that due to losses in the syphon the real emptying time will always
be somewhat higher.
<hr>
Next up (20.00 h UTC): a CHEMICAL example, with 'student participation' (euphemism alert).
[Edited on 19-4-2016 by blogfast25]
|
|
aga
Forum Drunkard
   
Posts: 7030
Registered: 25-3-2014
Member Is Offline
|
|
Jeez. All that algebra is seriously heavy stuff.
Are all the students going to draw lots to decide who does the experiment(s) ?
|
|
blogfast25
International Hazard
   
Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline
Mood: No Mood
|
|
Of course. But I've got some Darwins on a dead cert! 
Re. the algebra, the next one is much lighter on it. 
[Edited on 19-4-2016 by blogfast25]
|
|
blogfast25
International Hazard
   
Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline
Mood: No Mood
|
|
Hydrogen generator: (student problem)

In the hydrogen gas generator above, a (constant) volume V of dilute hydrochloric acid of initial concentration c0 reacts with zinc
granules acc.:
Zn(s) + 2 HCl(aq) === > ZnCl2(aq) + H2(g)
We assume the amount of Zn (number of moles) present to be much, much higher than that of HCl and the temperature to be constant. In those
circumstances kinetics allows us to write:
dNdt=kc
N is the number of moles of hydrogen formed, t is time, k the reaction rate constant and c the hydrogen chloride concentration.
The concentration of HCl is obviously NOT constant, as each mol of hydrogen generated requires 2 mol of HCl to react away. Stoichiometrically we can
write:
c=c0V−2NV
We start off at t=0 with N=0.
Develop expressions for:
1. N(t): the number of moles hydrogen generated as a function of time.
2. The rate of hydrogen generation dN(t)/dt as a function of time.
Heigh-ho, heigh-ho!
Off to work you go! 
[Edited on 19-4-2016 by blogfast25]
|
|
blogfast25
International Hazard
   
Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline
Mood: No Mood
|
|
Mini-recap on substitutions and differentiation
Suppose we have an integral of the type:
∫dxax+b
This can't be integrated directly because the stuff behind the d isn't the same as ax+b (differential doesn't match integrand).
To remedy this we try a substitution:
u=ax+b
Now derive:
dudt=adxdt+0=adxdt
Multiply both sides by dt to obtain the differential:
du=adx
So:
dx=1adu
Now substitute:
∫1aduu
=1a∫duu
Integrate this, then substitute back with u=ax+b.
Why this u-kerfuffle? Well, had we said that:
∫dxax+b=ln(ax+b) , that would have been incorrect because, if we take the first derivative, with the chain rule:
[ln(ax+x)]′=1ax+b(ax+b)′=aax+b≠1ax+n
But:
[1aln(ax+x)]′=1a×1ax+b(ax+b)′=1a×aax+b
[1aln(ax+x)]′=1ax+n
... which is correct.
<hr>
Here's a few more examples: remember that derivation and integration are anti-operators.
If:
y=∫cos3xdx=13sin3x+C
Then:
y′=[13sin3x+C]′=13×(sin3x)′+0
=13×cos3x×(3x)′=13×cos3x×3=cos3x
If:
y=∫(5−x)6dx=−17(5−x)7+C
Then:
y′=[−17(5−x)7+C]′=−17×7(5−x)6×(5−x)′+0=−(5−x)6×(−1)
=(5−x)6
If:
y=∫e3−2xdx=−12e3−2x+C
Then:
y′=[−12e3−2x+C]′=−12e3−2x×(3−2x)′+0
=−12e3−2x×(−2)
=e3−2x
[Edited on 22-4-2016 by blogfast25]
|
|
aga
Forum Drunkard
   
Posts: 7030
Registered: 25-3-2014
Member Is Offline
|
|
This sodding 'u' thing is doing my head in.
Let's see if i got it yet :
dNdt=kc
c=c0V−2NV
dNdt=kc0V−2NV
dNc0V−2N=kVdt
let
u=c0V−2N
Right. Here is where it gets messy.
du=u′=(c0V−2N)′
Left-hand part first
u′=1du=du
In the right-hand part only N is a variable, so we get the derivative of that, i.e. dN
(c0V−2N)′=(0−2)dN=−2dN
overall, we get
du=−2dN
therefore
dN=−12du
so now we can substitute for dN in the next step and get all integrobblers the same species
dNu=kVdt
and finally integrate
∫−12du1u=∫kVdt
−12∫duu=kV∫dt
−12ln(u)=kVt
substitute for u
−12ln(c0V−2N)=kVt
ln(c0V−2N)=−2kVt
c0V−2N=e−2ktV
2N=c0V−e−2ktV
N=12(c0V−e−2ktV)
i'm assuming the C generated in each integration is the same C, so cancel out.
|
|
aga
Forum Drunkard
   
Posts: 7030
Registered: 25-3-2014
Member Is Offline
|
|
I guess the rate of gas generation is the first derivative of that, which i think is
−12e−2kV
[Edited on 22-4-2016 by aga]
Hang on, surely the Rate will depend on the concentration as well ?
Doh.
[Edited on 22-4-2016 by aga]
|
|
blogfast25
International Hazard
   
Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline
Mood: No Mood
|
|
Heck. Everything was 100 % correct up to that point.
N=12(c0V−e−2ktV)
... is incorrect because:
−12ln(c0V−2N)=kVt
... should have been:
−12ln(c0V−2N)=kVt+C
⟹ln(c0V−2N)=−2kVt−2C
c0V−2N=e−2kVt−2C=e−2Ce−2kVt
Now I'm gonna use a little trick that isn't indispensable but just helps a bit. e<sup>-2C</sup> is a constant too, so I'll call it K:
c0V−2N=Ke−2kVt
Now with the initial condition:
t=0,N=0
c0V=K
So:
N=c0V2(1−e−2kVt)
And:
dNdt=kc0e−2kVt
We can see that these are correct because:
t→∞⟹N=c0V2
c0V mol of HCl does indeed generate (c0V)/2 mol of H2.
Also:
(dNdt)t=0=kc0
which is also correct.
That your solution was incorrect could be seen because you're subtracting a dimensionless number from moles. "apples minus oranges"
error!
Quote: Originally posted by aga  | I guess the rate of gas generation is the first derivative of that, which i think is
−12e−2kV
[Edited on 22-4-2016 by aga]
Hang on, surely the Rate will depend on the concentration as well ?
|
Indeed, Sir. Indeed. QED.
[Edited on 22-4-2016 by blogfast25]
|
|
aga
Forum Drunkard
   
Posts: 7030
Registered: 25-3-2014
Member Is Offline
|
|
Oh bollocks.
Big fat fecking C bollocks !
Ah well, did i get into the top 30% for effort, or do i remain in the bottom 20% ?
[Edited on 22-4-2016 by aga]
|
|
blogfast25
International Hazard
   
Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline
Mood: No Mood
|
|
Quote: Originally posted by aga  | Oh bollocks.
Big fat fecking C bollocks !
Ah well, did i get into the top 30% for effort, or do i remain in the bottom 20% ?
[Edited on 22-4-2016 by aga] |
Top 40 %!
One more problem coming up tonite, similar to this one (mathematically speaking). Hold beer till about 8 to 8.30 UTC. 
After that: partial derivatives!
[Edited on 22-4-2016 by blogfast25]
|
|
aga
Forum Drunkard
   
Posts: 7030
Registered: 25-3-2014
Member Is Offline
|
|
Erm, er ...
Edit:
Adding that impossible task into the equation makes this possibly the most difficult SM page of all time !
[Edited on 22-4-2016 by aga]
|
|
Pages:
1
..
8
9
10
11
12
..
19 |