Hydrostatic problem - Impossible integral

[jaumzaum]

Hydrostatic problem - Impossible integral!

Just studying hidrostatic over the internet and I saw the following problem:

A U-Tube filled with water, initially at rest in a horizontal table, has A1=40cm², A2 = 20cm², A3=30cm², h1 = 80cm, h2 = 20cm and L = 100cm (below piture). It is pressed in A1 by a constant vertical force F of 4N. Find the position x of the surface A1 in function of time (initial position is x=0). Given: g = 10m/s², density of water is 1kg/L

http://img191.imageshack.us/img191/2825/78925268.png

How can I solve this? I tried some stuff but it didn't worked.

A₁.dx₁ = A₂.dx₂ = A₃.dx₃ (where dx is the infinitesimal variation of position for a infinitesimal variation of time).

By energy conservation:

F.dx₁ = Δ(mechanic energy)

Potential Energy -> Now imagine that the dx₁.A₁ volume of water has just raised and become the dx₂.A₂ volume. It raised dx₁.(A₁+A₂)/(2A₂)
Ep = dx₁.(A₁+A₂)/(2A₂).dm.g
But dm = A1.dx₁.ρ
Ep = dx₁².(A₁+A₂)/(2A₂).A₁.ρ.g

Kinetic Enegrgy -> Ek = ∫(dm.v.dv)
dm = A₁.v₁.dt/ρ (v₁ is the instantaneum velocity for the surface A₁)
v = A₁.v₁/A (A is the area of the surface considered)
dv = A₁.dv₁/A

Ek = ∫[(A₁³v1².dt.dv₁)/(ρ.A²)]

But I don't know how to solve this. I don't know even if it is right.
Can anyone help me?

 

[SteamKing]

Remember, water is incompressible.

Apply Pascal's law to the two columns after application of the force.

 

[haruspex]

Apply Pascal's law to the two columns after application of the force.

Not sure that helps. The fluid is accelerating, so it's hydrodynamics, not hydrostatics. jaumzaum, I didn't understand your KE calculation. Can you not just add up the KEs of the three sections (ignoring complications at the corners)? At time t the force has advanced x, the volumes will be (h₁-x)A₁, LA₃ and (h₁+xA₁/A₂))A₂; the speeds v (=dx/dt), vA₁/A₃, and vA₁/A₂.
But the corners bother me. Not sure how to use h₂.

 

[jaumzaum]

I considered the process to be very slow (as the force is 4N). So I think the column of water tend to move orderly. I don't know how to use the h2 too! Actually this problem was from this university: http://www.ita.br/

It was a test applied on the school. You could use computers and you had to plot the graphic position xtime. The other exercises in the same test resulted in integrals that had to be series-aproximated, so I think you have to do some assumptions. But I don't know what!

[]'s
Joao

 

[jaumzaum]

Not sure that helps. The fluid is accelerating, so it's hydrodynamics, not hydrostatics. jaumzaum, I didn't understand your KE calculation. Can you not just add up the KEs of the three sections (ignoring complications at the corners)?

I think you can ignore the corners. I calculated KE for a random "surface" (it is actually a infinitesimal cylinder) of area A. Then I just added the 3 integrals for the 3 surfaces. But I still can't solve it! I mean, inside the integral is there dx and dt. To take the integral signal off I have to derivate. But derivate in relation to what? Time? Velocity? Position? And I also don't know how to derivate F.dx, for example! Another thing I noticed is that dx₁².(A₁+A₂)/(2A₂).A₁.ρ.g is very smaller than F.dx, can we ignore it?

How can I solve this?

 

[haruspex]

Working in terms of energy, you don't need dx's in the PE expression. Just consider that the total advancement of the force at time t is x. Abbreviating h1 to h:
PE = x².(A₁+A₂)/(2A₂).A₁.ρ.g
2*KE = (h-x)A₁ρv² + Lρv²A₁²/A₃ + (h+xA₁/A₂)ρv²A₁²/A₂
PE+KE = Fx
v²{h(1+A₁/A₂)+LA₁/A₃+x((A₁/A₂)²-1)} = 2Fx/(ρA₁) - x².g(A₁+A₂)/A₂

 

[Chestermiller]

Working in terms of energy, you don't need dx's in the PE expression. Just consider that the total advancement of the force at time t is x. Abbreviating h1 to h:
PE = x².(A₁+A₂)/(2A₂).A₁.ρ.g
2*KE = (h-x)A₁ρv² + Lρv²A₁²/A₃ + (h+xA₁/A₂)ρv²A₁²/A₂
PE+KE = Fx
v²{h(1+A₁/A₂)+LA₁/A₃+x((A₁/A₂)²-1)} = 2Fx/(ρA₁) - x².g(A₁+A₂)/A₂

This looks like a very nice development. The OP should also know that

v=dx/dt

This results in a first order ordinary differential equation for calculating x as a function of t. After taking the square root of both sides of the equation, it probably would be best to substitute
x=y² and solve for y.

 

[jaumzaum]

Working in terms of energy, you don't need dx's in the PE expression. Just consider that the total advancement of the force at time t is x. Abbreviating h1 to h:
PE = x².(A₁+A₂)/(2A₂).A₁.ρ.g
2*KE = (h-x)A₁ρv² + Lρv²A₁²/A₃ + (h+xA₁/A₂)ρv²A₁²/A₂
PE+KE = Fx
v²{h(1+A₁/A₂)+LA₁/A₃+x((A₁/A₂)²-1)} = 2Fx/(ρA₁) - x².g(A₁+A₂)/A₂

When I try to solve this (using Wolfram Mathematica of couse), the function does not accept the contour value of 0 for initial value of position.

v² (3.73333 + 3 x) = 2 x - 30 x²

http://www.wolframalpha.com/input/?i=(x')^2+(3.73333+++3+x[t])+==+2+x[t]+-+30+x[t]^2

What this means?

 

[haruspex]

As you can see, plugging in x=0 gives v=0, which is true. You would get the same trying to determine a trajectory under gravity using the energy equation y'² = 2gy. This stops any numerical attempt to produce a plot from there - things never move. You can validly get around it by plugging in a very small nonzero starting value.

 

[Chestermiller]

When I try to solve this (using Wolfram Mathematica of couse), the function does not accept the contour value of 0 for initial value of position.

v² (3.73333 + 3 x) = 2 x - 30 x²

http://www.wolframalpha.com/input/?i=(x')^2+(3.73333+++3+x[t])+==+2+x[t]+-+30+x[t]^2

What this means?

What it means is that you need to follow the advice I gave you in my previous response. Substitute v = dx/dt, divide both sides of the equation by (3.733 + 3x), take the square root of both sides of the equation, substitute x = y², and factor y out from both sides of the equation. Then solve the resulting first order ordinary differential equation for y(t). At short times, the solution should be of the approximate form y ~ kt, or x ~ at²/2, where a is the initial acceleration, and where a can be determined analytically in terms of the constants in your equation.

Wolframalpha won't be able to solve the equation until you work it into the form that I indicated.

 

[Chestermiller]

As you can see, plugging in x=0 gives v=0, which is true. You would get the same trying to determine a trajectory under gravity using the energy equation y'² = 2gy. This stops any numerical attempt to produce a plot from there - things never move. You can validly get around it by plugging in a very small nonzero starting value.

As I said in my previous post, plugging in a small starting value isn't necessary. If

y'² = 2gy

If you take the square root of both sides, you get:

$$y'=\sqrt{2gy}$$

Next, substitute y = z² to get

$$2z\frac{dz}{dt}=z\sqrt{2g}$$
or
$$\frac{dz}{dt}=\frac{\sqrt{2g}}{2}$$

This can now easily be integrated with respect to t. The same type of approach applies to the problem at hand.

 

[haruspex]

As I said in my previous post, plugging in a small starting value isn't necessary.

Sure - an analytic solution gets past the problem. Just wanted to explain that numerical methods can still be made to work when necessary.

 

[Chestermiller]

Sure - an analytic solution gets past the problem. Just wanted to explain that numerical methods can still be made to work when necessary.

Actually, I was alluding to using this approach (substitution of x = y²) to solve the problem numerically, starting out at y = x = 0 at t = 0. I've done this trick many times in obtaining numerical solutions to problems of this type of form.

 

[TSny]

PE+KE = Fx
v²{h(1+A₁/A₂)+LA₁/A₃+x((A₁/A₂)²-1)} = 2Fx/(ρA₁) - x².g(A₁+A₂)/A₂

I believe the factor shown in blue should be 1 + (A₁/A₂)^2. [EDIT: Nevermind, my mistake!]

Anyway, for the numbers given in the problem, I think it's safe to neglect the term in red. Then to a good approximation, the solution is SHM with a small amplitude. [Edit: I'm getting an amplitude of 10/3 cm and a period of about 2.2 s]