By Ingo Dahn, (Koblenz, Germany [email protected]) and Ambjörn Naeve (Stockholm, Sweden [email protected])
This notebook is part of an endeavour, to let the learner discover scientific concepts. Our tools shall be observations, models describing these operations and the critical inspection of these models in the light of new observations.
We take the historic emergence of these concepts in history as guidance, but may deviate from that route, follow dead ends or allow us be carried away by curiosity. We emphasize, that this is not a well-planned travel and by no means we claim that this is the only way to go. Quite the opposite - we are looking for travel mates, opening up alternative views, bridging gaps and pointing us, and the learner, to discoveries in the gardens and open landscapes aside of the road. We hope, that the interested reader/actor may take this as a travel guide, but more importantly we encourage walking ones' own roads.
We try to be highly interactive, supporting experiments, discovery and critical inspection, enjoying many flowers from the garden of open educational resources, into which we plant this notebook as well.
The first target of our travel is the concept of force as it emerged in science in ancient Greece in the form of buoyancy.
This notebook is made of a sequence of cells:
Text cells, explaining what we do. The text cells analyze a specific default balloon flight. They do not change and, depending on your working environment, you may not be able to change the content of text cells.
Hovering the mouse over the SageMath icon will reveal some additional explanations on the code - try it
Don't worry, if you don't understand the content of the code cells! We start simple and so you don't need to know a lot of Math, Physics or Programming - just go and experiment! Whenever you run into trouble - have a look at our .
License: This document is made available under a Creative Commons Attribution Share Alike License 4.0.
Observations:
Hypothesis (Aristotle, 384 BC - 322 BC): The natural place of a stone is the earth. Therefore no force is needed to let a stone fall. Moving a stone upwards is unnatural, force is only needed for unnatural motions.
Criticism: An arrow is flying through the air though apparently no force is acting on the arrow.
Improved Hypothesis: The air displaced through the arrow is providing the force that prevents the arrow from falling
Archimedes (287 BC - 212 BC):
Observation: Buoyancy lifts objects immersed in a fluid upward https://en.wikipedia.org/wiki/Buoyancy
Hypothesis (Archimedean Principle): A body immersed in a fluid experiences a buoyant force equal to the weight of the fluid it displaces.
Note: Unlike Aristotle, Archimedes doesn't just talk about the nature of Buoyancy, he gives a way to calculate it - Math enters Physics!
Example: A stone with a weight of 20g displaces a water volume of wight 5g and loses a weight of 5g.
Stone and water are connected by the volume: The volume of the stone equals the volume of the displaced water.
Volume and weight are connected by the concept of density. This concept is based on the observation, that the weight of a compact object is proportional to it's volume - if you double (triple) it's volume, the weight is doubled (tripled) as well. The proportionality factor of this relation is nothing else but the density ρ:weight=ρ⋅volume resp. Density ρ is the quotient of weight and volume and is measured in cm3g or m3kg. Can you translate between these units?
Note: You may learn later, that Archimedes should perhaps have used mass instead of weight but he had no reason to and therefore we shall also not distinguish between mass and weight in this notebook. This distinction must be discussed elsewhere.
If we assume that water has a density of 1cm3g, the water displaced and the stone must have a volume of 5cm3. The density of the stone will be 520cm3g=4cm3g.
If we denote the weight of the object (stone) and of the displaced fluid (water) by mo,mf respectively and the density of the object respectively fluid by ρo,ρf, then ρomo=ρfmf. According to the definition of density, this equation just states that the volume of the stone equals the volume of the water displaced.
Now Archimedes' principle says that the difference in the measurement of weight of the object outside and inside the fluid is equal to the weight of the fluid displaced-
Question (King Hiero II): King Hiero II from Syracuse provided a certain amount of gold to a goldsmith for the delivery of a votive pure crown. The king asks you (Archimedes) to determine whether some silver had been substituted by the dishonest goldsmith and, if so, how much gold has been displaced with silver.
Data: Density of water: ρf=0.997cm3g, Density of pure gold: ρg=19.3cm3g. mass of the crown: mo=400g, volume of the water displaced by the crown: vf=22.3cm3.
m_o=400 v_f=22.3 rho_g=19.3 rho_o=m_o/v_f print("The density of the crown is",rho_o,"g/ccm which is less than the density of gold of",rho_g,"g/ccm.")
But how much of the gold has the goldsmith taken?
Suppose the crown with a weight of 400g contains s gram of silver. The density of silver is ρs=10.49cm3g.
Then the crown must have a volume of 22.3=ρg400−s+ρss. Solving this equation for s tells us, how much the goldsmith has taken.
rho_s=10.49 var('s') eq=22.3 == (400-s)/rho_g +s/rho_s sol=solve(eq,s) show("Solution:",sol) r=sol[0].rhs().n() print("The goldsmith has replaced",r,"g of gold by silver.")
Fooling Archimedes
Had the goldsmith a better knowledge of physics, he might have fooled archimedes by adding a 3rd substance which is more dense than gold. A compilation of the density of various materials is available at Wikipedia. Unfirtunately, such materials tend to have unwanted properties like a high price or being radioactive. Anyhow, as a thought experiment, let's find out how much of the densest, Osmium, the goldsmith would have to add to compensate for the loss of that amount of gold.
The density of Osmium is ρos=22.57 cm3g. A crown consisting of 400−r g gold, s g silver and o g Osmium would have a weight of 400−r+g+o g and a volume of ρg400−r+ρss+ρoso, while a pure gold crown of 400 g would have a volume of ρg400. Thus we need to solve the equation
ρg400−r+ρss+ρoso = ρg400
which is equivalent with
ρss+ρoso = ρgr
We also know that
(400−r)+s+o = 400,
hence
o+s = r.
We must solve a linear equation system.
var('o') rho_os=22.57 eqs=[s/rho_s + o/rho_os == r/rho_g, o+s == r] sol2=solve(eqs,(s,o)) show('Solution:',sol2) ss=sol2[0][0].rhs().n() oo=sol2[0][1].rhs().n() print("To fool Archimedes, replacing",r,"g of gold by silver and osmium, the goldsmith would need",ss,"gram of silver and",oo,"g of osmium.")
According to Archimedes' principle, the weight of the object in the fluid is reduced by the weight of the displaced fluid. This reduction is called buoyancy.
The buoyancy can be calculated as the weight of the fluid displaced by the volume of the object, i.e. b=vo⋅ρf=ρomo⋅ρf=moρfρo. It is worth noting that the buoyancy does not depend on the form of the object
Experience buoyancy in the following virtual lab from Dimitris Mirogiannis.
Note: This lab measures the buoyant force in Newton (N). We have not discussed the reasons for using the Newton unit. For now you may assume that 1 N = 102 g and 1 kg = 9.81 N.
Question: Can you calculate the buoyancy of objects made from various materials?
When the object has the same density as the fluid (ρo=ρf), the buoyancy compensates the mass of the object - the object floats in the fluid. If the object's density is less than that of the fluid, buoyancy is greater than the mass of the object, buoyancy provides an upthrust and the object drifts upwards until it reaches the surface of the fluid and swims.
In the latter case, the calculation of the weight of the object immersed into the fluid yields a negative result, which means that the buoyancy acts as an upward directed force on the object.
Objects that are hollow (like ships) have a lower density than the materials they consist of. In the following Question you can explore the effect. If d) yields a negative value, the ball will float to the surface and swim.
We have seen that an object will float in a fluid if it's density equals the density of the fluid.
Question: Can you calculate, how large a ball must be inorder to float or swim?
If a hollow ball has a radius of r cm and a hull thickness of h cm, it's volume is vouter=34π⋅r3 cm3 The volume of the empty space inside is vinner=34π⋅(r−h)3 cm3. The volume of the hull alone is vouter−vinner and the mass of the hollow ball is (vouter−vinner)ρm, where ρm is the density of the material of the hull. Hence, the density of the hollow ball is ρo=vouter(vouter−vinner)ρm (we ignore the weight of air inside).
The question is, how big r must be such that the ball begins to float in the water, i.e. such that ρo=ρw.
We calculate this for a ball of glass (ρm=2.5) with a hull thickness of 2 mm, but feel free to change these values.
h=0.2 rho_w=0.997 rho_m=2.5 var('r') v_outer(r) = (4/3)*pi*r^3 v_inner(r) = (4/3)*pi*(r-h)^3 rho_o(r) = (v_outer(r)-v_inner(r))*rho_m/v_outer(r) rr=solve(rho_o(r) == rho_w,r) rrr=rr[2].rhs().n() print("The ball will swim in water if it's radius is greater than",rrr,"cm.") p1=plot(rho_o,(1,5),legend_label="Density of object") p2=plot(rho_w,(1,5),color='red',legend_label="Density of water") show(p1+p2,axes_labels=["Radius","Density"])
The following diagram shows a section of a hollow glass ball, immersed in water, and the buoyancy. You can vary the radius and the thickness of the hull of the ball with the sliders, either with your mouse or with your arrow keys. The cell reuses values from previous cells. Since the thickness h of the hull now varies, the volume vinner of the inner ball and the density of the ball ρo now depend also on h. w(r,h)=ρo(r,h)vouter(r)−ρwvouter is the loss of weight by immersion into water.
Questions: Try to find values for the radius and the thickness of the hull such that the ball floats in the water. Set, at the start of the cell, new values for the densities rho_m
and rho_w
of glass and water to experiment with other materials.
var('r,h') v_inner(r,h) = (4/3)*pi*(r-h)^3 rho_o(r,h) = (v_outer(r)-v_inner(r,h))*rho_m/v_outer(r) w(r,h)=v_outer(r)/rho_o(r,h)-v_outer(r)/rho_w @interact def _(r=slider(1,5,step_size=0.05),h=slider(0.2,1,step_size=0.05)): Outer=circle((0,0),r,fill=True,zorder=1) Inner=circle((0,0),r-h,fill=True,color='white',zorder= 2) B=arrow((0,0),(0,w(r,h)),color='black',zorder=3) show(Outer+Inner+B,xmin=-5,xmax=5,ymin=-5,ymax=5)
The following image shows the density of the hollow ball, depending on the radius as well as on the thickness of the hull, where the thickness h
of the hull varies between 0.2 and 1. You can drag the diagram around with the mouse. Where the blue surface is below the red plane, the ball swims. You should recognize the curve above in the foreground when the label Radius=3.00
is shown at the lower axis in the background.
Note that the thickness of the hull, stretching from the foreground to the background in the view described above, has a bigger effect than the radius!
Plot_o=plot3d(rho_o,(r,1,5),(h,0.2,1),opacity=0.7) Plot_w=plot3d(rho_w,(r,1,5),(h,0.2,1),color='red',opacity=0.3) show(Plot_o+Plot_w,axes_labels=["Radius","Hull","Density"])
Question: Experiment!
(h,0.2,1)
in lines 4 an 5 with (h,0.2,5)
to inspect the whole range from a small hull to a compact ball.Think: How can the consideration in this notebook be applied to calculate the buoyancy of balloons floating in the air, when the density of the air decreases when the balloon gets higher? Experiment with our notebooks on weather balloons!
Bonus: Did you know, the the United States built swimming ships from concrete during the second world war? (Note: The density of concrete is 2.4cm3g, that of iron is 7.87cm3g):
You may use the following cell for your own experiments.