If I was still doing consulting engineering I would look for a rather simple practical approach that represents a good proxy for the behavior of the proposed system to achieve reasonable practical relative results within the client’s budget.
Assumptions and approximations (see attached sketch for layout):
1) The area of the round piston Ap is >> Ao (in the range Ao/Ap < 0.1) and Ao is a round orifice area and the cylinder is prismatic container
2) The internal pressure that the piston exerts on the fluid in the cylinder is constant over time and the weight added onto the piston is gently placed on it
3) The flow is 1D and quasi steady state
4) There is no leakage of fluid up along the perimeter sides of the piston. Friction is neglected for a seal on the piston against the cylinder wall
5) The fluid in the cylinder has unit weight w. Ex if water, w = 64.5 lbs/cu. ft and varies with temperature if temp is considered.
6) The control volume region of analysis has streamlines that are parallel to the vertical sides of the container. Therefore near the orifice Ao calculations should refrain from the piston getting to close to the Ao at the bottom where 3D effects start to predominate flow toward the orifice
7) The fluid in the container is incompressible and conservation of volume holds such that the
Ap x dy(t) volume over time t equals the volume that leaves the container through Ao over the same time interval.
8) The discharge through the orifice with assumptions 1 and 3 at any time from the energy and continuity equations is approximated by:
Q = CAo (2gy)^ 0.5 where:
C = Cd / [(1 – n(Ao/Ap)^2]^0.5 where Cd is the discharge coefficient of the orifice with values that
range from about 0.6 for sharp edge to 0.98 for a rounded entrance. Refer to hydraulic text (or
Google) for other orifice conditions and Cd values. It is further assumed that Cd does not vary with
water depth and the orifice is in a thin plate with no additional piping as indicated in the poster’s
sketch. The value of n = 1 just using the energy and continuity equations and n = 2 for Ao/Ap < about
0.5 was derived by Chestermiller in a rather detailed analysis located in his post #12:
The system can reasonably be represented by a tank that drains down through an orifice with the addition of a weighted piston exerting additional internal force on the fluid to drain the cylinder faster than a no piston container with an open top and an same orifice.
The internal imposed pressure of the piston and weight is:
Pp = (W + Wp)/Ap where W = weight put on piston and Wp = piston weight
This can be converted to equivalent fluid head by dividing by fluid unit weight w
Hp = Pp/ w = (W + Wp)/(wAp) = constant uniform pressure head over time
A variable quasi steady draining container can be represented by:
Q = flow = - d(volume)/dt = approx. CAo(2gy)0.5 ; negative sign for time increases as y decreases, using
d(volume) = Ap dy
(Ap dy)/dt = - CAo(2gy)^0.5
dt = - Ap dy / ( CAo(2gy)^ 0.5
Integrating from time t = 0 at (Hp + h) to a later time Tt at (Hp + y) where h is the starting height of the piston in the cylinder and y is the piston position at Tt and y < h:
Tt = from t = 0 to Tt = { (2Ap / {CAo(2g)^0.5) } x [(Hp + h)^ 0.5 – (Hp + y)^0.5]
If y = h, Tt = 0 at start position of piston, if Hp >> h, Tt approaches 0 not desired for Hp value
Select values for Ao, Ap, h, w and Cd, and a y < h and calculate Tt then:
Average Velocity of piston (Vt) over the drop distance (at time Tt):
Vt = (h - y) / Tt with piston velocity = 0 at time = 0
Average Acceleration of piston (At) over the drop distance (at time Tt):
At = Vt / Tt = (h –y) / (Tt)^2 with acceleration = 0 at time = 0.
Note assumption 4): friction of piston seal neglected so the calculated Tt will be greater than Tt. These equations should provide practical relative changes in piston drop times, piston velocity and acceleration for comparisons with different Ao, Ap, h, w, Cd and y < h.
Spread sheet it to play around
