Sciencemadness Discussion Board
Not logged in [Login ]
Go To Bottom

Printable Version  
 Pages:  1  ..  8    10    12  ..  19
Author: Subject: Calculus! For beginners, with a ‘no theorems’ approach!
aga
Forum Drunkard
*****




Posts: 7030
Registered: 25-3-2014
Member Is Offline


[*] posted on 17-4-2016 at 08:27


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. :cool: 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.




View user's profile View All Posts By User
blogfast25
International Hazard
*****




Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline

Mood: No Mood

[*] posted on 17-4-2016 at 08:53


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[TTaTiTa]=τt


You would plot in Excell:

ln[TTaTiTa]


versus:

t


That should give something looking like:

line fitted.png - 3kB

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.




View user's profile View All Posts By User
aga
Forum Drunkard
*****




Posts: 7030
Registered: 25-3-2014
Member Is Offline


[*] posted on 17-4-2016 at 10:46


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 ?




View user's profile View All Posts By User
blogfast25
International Hazard
*****




Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline

Mood: No Mood

[*] posted on 17-4-2016 at 11:53


Checking surface area of the cup, by calculus.

For a revolution body:

Revolution body.png - 3kB

... 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π8501+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=mgky

or:
my=mgky


Newtonian notation ('overdot notation'):
m¨y=mgky


... 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]




View user's profile View All Posts By User
aga
Forum Drunkard
*****




Posts: 7030
Registered: 25-3-2014
Member Is Offline


[*] posted on 17-4-2016 at 12:50


Quote: Originally posted by blogfast25  
md2ydt2=mgky

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.




View user's profile View All Posts By User
blogfast25
International Hazard
*****




Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline

Mood: No Mood

[*] posted on 17-4-2016 at 13:39


Quote: Originally posted by aga  
Quote: Originally posted by blogfast25  
md2ydt2=mgky

Is 'g' gravitational acceleration ?



Yes. The bob is subject to it's own weight:

mg




View user's profile View All Posts By User
aga
Forum Drunkard
*****




Posts: 7030
Registered: 25-3-2014
Member Is Offline


[*] posted on 17-4-2016 at 14:07


Of course, same as everything in agaspace, including Time : all attributes are inseparable.
g=Σallototherbits

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]




View user's profile View All Posts By User
blogfast25
International Hazard
*****




Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline

Mood: No Mood

[*] posted on 17-4-2016 at 15:06


@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:
fx

... and to y:
fy


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]




View user's profile View All Posts By User
blogfast25
International Hazard
*****




Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline

Mood: No Mood

[*] posted on 18-4-2016 at 04:58


Here's the plot for:

ln[TTaTiTa]=τt


[TTaTiTa]


... is also called the "Reduced Temperature", a dimensionless number between 1 and 0.

I used Ta=20 C.

Reduced temperature.png - 9kB

The linear regression equation:

y=0.0006x0.021


... gives us an estimate for:

τ=0.0006s1

τ=hAmcp


With A=0.053 m2, m=0.458 kg, cp=4181 J/kg K then:

h=τmcpA=22Wm2K1


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]




View user's profile View All Posts By User
aga
Forum Drunkard
*****




Posts: 7030
Registered: 25-3-2014
Member Is Offline


[*] posted on 18-4-2016 at 09:08


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 ?




View user's profile View All Posts By User
blogfast25
International Hazard
*****




Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline

Mood: No Mood

[*] posted on 18-4-2016 at 09:14


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!




View user's profile View All Posts By User
aga
Forum Drunkard
*****




Posts: 7030
Registered: 25-3-2014
Member Is Offline


[*] posted on 18-4-2016 at 10:45


Quote: Originally posted by blogfast25  


Perhaps the most common f(x,y,z) type function is the volume of a right angled box:
V=lhw

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)




View user's profile View All Posts By User
blogfast25
International Hazard
*****




Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline

Mood: No Mood

[*] posted on 18-4-2016 at 11:39
Syphon: tank emptying time


Syphon problem.png - 8kB

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=c2gy

So the volumetric throughput Qv is:
Qv=π4D2v2=π4D2c2gy


During an infinitesimal amount of time dt, the water level decreases by dy (<0), so:
dV=Ady=Qvdt

Ady=π4D2c2hydt=π4D2c2gydt

Let's call:
α=πD2c2g4A

y1/2dy=αdt

Integrate between relevant boundaries:
h0Hy1/2dy=t0αdt

2[y]Hh0=αt

t=2α(Hh0)

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:
L2(Hh0)+h0

Substituting we get:
12c2v22=gy32μLv2ρD2

12v22+βv2c2gy=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(v22v21)+p0ρ=gh0+p1ρ

A2v2=Av1

v22v21=v22[1(A2A)2]

A2=π4D2

1c2=[1(A2A)2]

12c2v22+p0ρ=gh0+p1ρ

p1=ρg(yh0)+p0

12c2v22=gy

v2=c2gy

Qv=π4D2v2=π4D2c2gy

dV=Ady=Qvdt

Ady=π4D2c2hydt=π4D2c2gydt

α=πD2c2g4A

y1/2dy=αdt

h0Hy1/2dy=t0αdt

2[y]Hh0=αt

t=2α(Hh0)





But.
12c2v22=gyΔpsρ

Δps=32μLv2D2

L2(Hh0)+h0

12c2v22=gy32μLv2ρD2

12v22+βv2c2gy=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]




View user's profile View All Posts By User
aga
Forum Drunkard
*****




Posts: 7030
Registered: 25-3-2014
Member Is Offline


[*] posted on 19-4-2016 at 01:34


Cor Blimey !

I think a bit of commentary is going to be required to follow that properly !




View user's profile View All Posts By User
blogfast25
International Hazard
*****




Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline

Mood: No Mood

[*] posted on 19-4-2016 at 05:44


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)=π4D2c2gy


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α(Hh0)


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α(Hh0)


... 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]




View user's profile View All Posts By User
aga
Forum Drunkard
*****




Posts: 7030
Registered: 25-3-2014
Member Is Offline


[*] posted on 19-4-2016 at 08:09


Jeez. All that algebra is seriously heavy stuff.

Are all the students going to draw lots to decide who does the experiment(s) ? ;)




View user's profile View All Posts By User
blogfast25
International Hazard
*****




Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline

Mood: No Mood

[*] posted on 19-4-2016 at 09:35


Quote: Originally posted by aga  

Are all the students going to draw lots to decide who does the experiment(s) ? ;)


Of course. But I've got some Darwins on a dead cert! :D

Re. the algebra, the next one is much lighter on it. :)

[Edited on 19-4-2016 by blogfast25]




View user's profile View All Posts By User
blogfast25
International Hazard
*****




Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline

Mood: No Mood

[*] posted on 19-4-2016 at 12:18
Hydrogen generator: (student problem)


Hydrogen generator.png - 6kB

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=c0V2NV




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!
:D

[Edited on 19-4-2016 by blogfast25]




View user's profile View All Posts By User
blogfast25
International Hazard
*****




Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline

Mood: No Mood

[*] posted on 21-4-2016 at 13:27
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

=1aduu

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+b1ax+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=(5x)6dx=17(5x)7+C

Then:
y=[17(5x)7+C]=17×7(5x)6×(5x)+0=(5x)6×(1)

=(5x)6

If:
y=e32xdx=12e32x+C

Then:
y=[12e32x+C]=12e32x×(32x)+0

=12e32x×(2)

=e32x


[Edited on 22-4-2016 by blogfast25]




View user's profile View All Posts By User
aga
Forum Drunkard
*****




Posts: 7030
Registered: 25-3-2014
Member Is Offline


[*] posted on 22-4-2016 at 06:59


This sodding 'u' thing is doing my head in.

Let's see if i got it yet :

dNdt=kc

c=c0V2NV

dNdt=kc0V2NV

dNc0V2N=kVdt

let
u=c0V2N

Right. Here is where it gets messy.
du=u=(c0V2N)

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
(c0V2N)=(02)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

12duu=kVdt

12ln(u)=kVt

substitute for u
12ln(c0V2N)=kVt

ln(c0V2N)=2kVt

c0V2N=e2ktV

2N=c0Ve2ktV

N=12(c0Ve2ktV)


i'm assuming the C generated in each integration is the same C, so cancel out.




View user's profile View All Posts By User
aga
Forum Drunkard
*****




Posts: 7030
Registered: 25-3-2014
Member Is Offline


[*] posted on 22-4-2016 at 07:13


I guess the rate of gas generation is the first derivative of that, which i think is
12e2kV


[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]




View user's profile View All Posts By User
blogfast25
International Hazard
*****




Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline

Mood: No Mood

[*] posted on 22-4-2016 at 08:05


Quote: Originally posted by aga  

i'm assuming the C generated in each integration is the same C, so cancel out.



Heck. Everything was 100 % correct up to that point.

N=12(c0Ve2ktV)


... is incorrect because:

12ln(c0V2N)=kVt


... should have been:


12ln(c0V2N)=kVt+C

ln(c0V2N)=2kVt2C

c0V2N=e2kVt2C=e2Ce2kVt

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:
c0V2N=Ke2kVt

Now with the initial condition:
t=0,N=0

c0V=K

So:
N=c0V2(1e2kVt)

And:
dNdt=kc0e2kVt

We can see that these are correct because:
tN=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
12e2kV


[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]




View user's profile View All Posts By User
aga
Forum Drunkard
*****




Posts: 7030
Registered: 25-3-2014
Member Is Offline


[*] posted on 22-4-2016 at 08:49


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]




View user's profile View All Posts By User
blogfast25
International Hazard
*****




Posts: 10562
Registered: 3-2-2008
Location: Neverland
Member Is Offline

Mood: No Mood

[*] posted on 22-4-2016 at 10:07


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. :D

After that: partial derivatives!

[Edited on 22-4-2016 by blogfast25]




View user's profile View All Posts By User
aga
Forum Drunkard
*****




Posts: 7030
Registered: 25-3-2014
Member Is Offline


[*] posted on 22-4-2016 at 10:11


Quote: Originally posted by blogfast25  
Hold beer till about 8 to 8.30 UTC. :D

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]




View user's profile View All Posts By User
 Pages:  1  ..  8    10    12  ..  19

  Go To Top