from __future__ import division
# Soltion
print "Example: 1.1 - Page: 4\n\n"
#*****Data*****#
g_Earth = 9.83 # [m/square s]
F_Earth = 800 # [N]
g_Moon = 3.2 # [m/square s]
#************#
# From the expression of force, the force on the man on the Eath's surface is given by:
# F = m*g_Earth
m = F_Earth/g_Earth # [kg]
# On the moon, the weight of the mass is equal to the force acting on the mass on the moon and is given by
F_Moon = m*g_Moon # [N]
print "Weight of the man on the moon is %0.2f N\n"%(F_Moon)#
# Solution
print "Example: 1.2 - Page: 5\n\n"
#*****Data*****#
m1 = 1.5 # [mass of the body, kg]
m2 = 6*10**(24) # [mass of the Earth, kg]
G = 6.672*10**(-11) # [N.square m/square.kg]
r = 6000*10**(3) # [m]
#************#
# According to Newton's universal law of gravity:
F = G*m1*m2/r**2 # [N]
print "Gravitational force on the body is %.2f N\n"%(F)#
# Solution
print "Example: 1.3 - Page: 5\n\n"
#*****Data*****#
r_Moon = 0.3 # [km]
r_Earth = 1 # [km]
m2 = 1 # [mass of body, kg]
mMoon_By_mEarth = 0.013 # [kg/kg]
#***************#
# According to the Newton's universal law of gravitation:
Fe_By_Fm = (1/mMoon_By_mEarth)*(r_Moon/r_Earth)**2#
print "Mass of 1 kg will weigh %.2f kg on moon\n"%(Fe_By_Fm)#
# Solution
print "Example: 1.4 - Page: 6\n\n"
#*****Data*****#
h = 40 # [cm]
density = 14.02 # [g/cubic cm]
g = 9.792 # [m/square s]
#*************#
P = h*density*g/1000 # [N/square cm]
P = P*10 # [kPa]
print "The absolute pressure is %.3f kPa\n"%(P)#
# Solution
print "Example: 1.5 - Page: 7\n\n"
#*****Data*****#
Patm = 112 # [kPa]
density = 1200 # [kg/cubic m]
g = 9.81 # [m/sqaure s]
h = 0.62 # [m]
#**************#
P = Patm + (density*g*h/1000) # [kPa]
print "The absolute pressure within the container is %.3f kPa\n"%(P)#
# Solution
print "Example: 1.6 - Page: 9\n\n"
#*****Data*****#
F = 150 # [N]
Displacement = 10 # [m]
#**************#
W = F*Displacement # [J]
print "Work done by the system is %d J\n"%(W)#
# Solution
print "Example: 1.7 - Page: 9\n\n"
#*****Data*****#
P = 560*10**3 # [Pa]
Vinit = 3 # [cubic m]
Vfinal = 5 # [cubic m]
Wext = 210*10**3 # [J]
#*************#
W = P*(Vfinal - Vinit) # [J]
# Again the system receives 210 kJ of work from the external agent.
W = W - Wext # [J]
print "Actual Work done by the system is %.1e J\n"%(W)#
# Solution
print "Example: 1.8 - Page: 11\n\n"
#*****Data*****#
g = 9.81 # [m/square s]
Z = 100 # [m]
#***************#
# Basis: 1 kg of water
m = 1 # [kg]
Ep = m*g*Z # [J]
print "Change in potential Energy is %d J\n"%(Ep)
from __future__ import division
# Solution
print "Example: 1.9 - Page: 11\n\n"
#*****Data*****#
m = 15# # [kg]
g = 9.81 # [m/square s]
V1 = 0 # [m/square s]
Z1 = 12 # [m]
Z2 = 0 # [m]
#***************#
# At initial condition, V1 = 0, so kinetic energy is zero.
# At final condition, Z2 = 0, so potential energy is zero.
# Ep1 + Ek1 = Ep2 + Ek2
#deff('[y] = f(V2)','y = ((1/2)*m*V1**2) + (m*g*Z1) - (((1/2)*m*V2**2) + (m*g*Z2))')#
def f(V2):
y = ((1/2)*m*V1**2) + (m*g*Z1) - (((1/2)*m*V2**2) + (m*g*Z2))
return y
from scipy.optimize import fsolve
V2 = fsolve(f,7)#
print "The velocity of the metal block is %.2f m/s\n"%(V2)#
Ek2 = (1/2)*m*V2**2 # [J]
print "The final Kinetic Energy is %.1f J\n"%(Ek2)#
from __future__ import division
# Solution
print "Example: 1.10 - Page: 12\n\n"
#*****Data*****#
m = 1200 # [kg]
v1 = 10 # [km/h]
v2 = 100 # [km/h]
time = 1 # [min]
#***************#
v1 = 10*1000/3600 # [m/s]
v2 = 100*1000/3600 # [m/s]
W = (1/2)*m*(v2**2 - v1**2) # [J]
time = time*60 # [s]
P = W/time # [W]
print "Power required is %.2f kW\n"%(P/1000)#
from __future__ import division
from math import pi
print "Example: 1.11 - Page: 13\n\n"
#*****Data*****#
dia = 0.3 # [m]
m = 100 # [kg]
P_atm = 1.013*10**5 # [N/square m]
g = 9.792 # [m/square s]
#**************#
Area = (pi/4)*dia**2 # [square m]
#Solution (a)(i)
# Force exerted by the atmosphere:
F_atm = P_atm*Area # [N]
# Force exerted by piston & metal block:
F_mass = m*g # [N]
# Total force acting upon the gas:
F = F_atm + F_mass # [N]
print "Total Force eacting upon the gas is %.1f N\n"%(F)#
# Solution (a)(ii)
Pressure = F/Area # [N/square m]
print "Pressure exerted is %.3f kPa\n\n"%(Pressure/1000)#
# Solution (b)
# The gas expands on application of heat, the volume of the gas goes on increasing and the piston moves upward.
Z = 0.5 # [m]
# Work done due to expansion of gas:
W = F*Z # [J]
print "Work due to expansion by the gas is %.3f kJ\n\n"%(W/1000)#
# Solution (c)
# Change in potential energy of piston and weight after expansion process:
Ep = m*g*Z # [J]
print "Change in Potential Energy is %.1f J\n"%(Ep)#
from __future__ import division
print "Example: 1.12 - Page: 24\n\n"
# Solution
# The relation is:
# (C/5) = ((F - 32)/9)
# For C = F
C = - (32*5/4)# # [degree Celsius]
print "The temperature which has the same value on both the centigrade and Fahrenheit scales is %d degree Celsius or %d degree Fahrenheit\n"%(C,C)#
from __future__ import division
print "Example: 1.13 - Page: 24\n\n"
# Solution
#*****Data*****#
delta_T_C = 30 # [OC]
#*************#
# The relation between the Kelvin temperature scale and the Celsius temperature scale:
# T(K) = T(OC) + 273.15
# Here, the temperature rise is to be expressed in terms of K, but the difference in temperature will be the same in the Kelvin and Celsius scales of temperature:
delta_T_K = delta_T_C # [K]
print "The rise in temperature in the Kelvin scale is %d K\n"%(delta_T_K)#
# The emperical relationship between the Rankine and Kelvin scales is given by:
delta_T_R = 1.8*delta_T_K # [R]
print "The rise in temperature in the Rankine scale is %d R\n"%(delta_T_R)#
delta_T_F = delta_T_R # [OF]
print "The rise in temperature in the Fahrenheit scale is %d OF\n"%(delta_T_F)#