In [17]:

```
'''A resistance oven employing nichrome wire is to be operated from 220 V single-phase supply and is to be rated at 16 kW.
If the temperature of the element is to be limited to 1,170°C and average temperature of the charge is 500°C, find the
diameter and length of the element wire.
Radiating efficiency = 0.57, Emmissivity=0.9, Specific resistance of nichrome=(109e–8)ohm-m.'''
P = 16000.0 #W (output power)
V = 220.0 #V (applied voltage)
rho = 109.0e-8 #ohm-m (resistivity)
e = 0.9 # (Emmisivity)
K = 0.57 # (Radiating efficiency)
T1 = 1170.0 + 273.0 #K (Temp of hot body)
T2 = 500.0 + 273.0 #K (Temp of cold body)
#Now , l/d^2 = pi*V^2/4*rho*P = a .Therefore a is
a = 3.14*(V**2)/(4*rho*P) # (a = l/d^2) ------ 1
#Using Stefan's law of radiation
H = 5.72*e*K*((T1/100)**4-(T2/100)**4) #W/m^2
#Total heat dissipated = electrical power input
# (pi*d)*l*H = P . Therefore ,let b = l*d. So,
b = P/(H*3.14)
b2 = b**2 #(b2 = l^2*d^2)------------------------ 2
#Multiplying 1 and 2.
l3 = a*b2 # ( = l^3)
l = l3**(1/3.0) #m (length)
d = b/l*1000 #mm (diameter)
print "Length of wire =",round(l,2),"m."
print "Diameter of wire =",round(d,2),"mm."
```

In [18]:

```
'''A 30-kW, 3-phase, 400-V resistance oven is to employ nickel-chrome strip 0.254 mm thick for the three star-connected heating
elements.If the wire temperature is to be 1,100°C and that of the charge to be 700°C, estimate a suitable width for the strip.
Assume emissivity = 0.9 and radiating efficiency to be 0.5 and resistivity of the strip material is 101.6e-8 ohm- m.What would be
the temperature of the wire if the charge were cold ?'''
import math as m
P = 30.0*1000 *(1/3.0) #W (Power/phase)
V = 400.0/m.sqrt(3) #V (Phase voltage)
rho = 101.6e-8 #ohm-m (resistivity)
e = 0.9 # (emmisivity)
K = 0.5 # (radiating efficiency)
t = 0.254e-3 #m (thickness of strip)
T1 = 1100.0 + 273 #K (Temp. of hot wire)
T2 = 700.0 + 273 #K (Temp. of charge)
R = V*V/P #ohm (Resistance => P = V^2/R)
# R = rho*l/(w*t) l/w = R*t/rho = a.Therefore a is
a = R*t/rho #(=l/w)------------------------------------ 1
#Using stefan's law
H = 5.72*e*K*((T1/100)**4-(T2/100)**4) #W/m^2
#Surface area of strip = 2*w*l .
#Total Heat dissipated = electrical power => wl*2H = P . Let b = wl
b = P/(2*H) #(=wl)----------------------------------------- 2
#Dividing 1 by 2
w2 = b/a #(=w*w)
w = m.sqrt(w2)*1000 #mm (width)
print "Width of strip =",round(w,2),"mm."
```

In [19]:

```
'''A cubic water tank has surface area of 6.0 m^2 and is filled to 90% capacity six times daily. The water is heated from 20°C
to 65°C.The losses per square metre of tank surface per 1°C temperature difference are 6.3 W. Find the loading in kW and the
efficiency of the tank.Assume specific heat of water = 4,200 J/kg/°C and one kWh = 3.6 MJ.'''
import math as m
sa = 6.0 #m^2 (surface area)
T2 = 65.0 #*C (final temp)
T1 = 20.0 #*C (initial temp)
loss = 6.3 #W/*C/m^2 (loss per square metre per 1*C)
s = 4200.0 #J/Kg/*C (Specific heat)
#Now, sa = 6*l^2
l = m.sqrt(sa/6) #m
#Volume = l^3
V = l**3 #m^3
#Volume of water to be heated daily is
V2 = 6*V*0.9 #m^3
#Since 1m^3 = 1000 kg => mass of water to be heated is
mass = V2*1000.0 #kg
#Heat required to raise temp =
H = mass*s*(T2-T1) #MJ
H = H/(3.6*10e+5) #kWh
#Daily loss from surface
L = 6*loss*(T2-T1)*(24.0/1000) #kWh
#Total energy required
Tot = L + H #kWh
#(i) Loading in KW is
load = Tot/24 #kW
#(ii)Efficiency of tank is
eff = (H/Tot)*100.0 # (% efficency)
print "Loading in kW = ",round(load,2),"kW."
print "Efficency of Tank =",round(eff,2),"%"
```

In [20]:

```
'''A 4-phase electric arc furnace has the following data :
Current drawn = 5000 A ; Arc voltage = 50 V
Resistance of transformer referred to secondary = 0.002 ohm
Resistance of transformer referred to secondary = 0.004 ohm
(i) Calculate the power factor and kW drawn from the supply.
(ii) If the overall efficiency of the furnace is 65%, find the time required to melt 2 tonnes of steel if
latent heat of steel = 8.89 kcal/kg, specific heat of steel = 0.12, melting point of steel = 1370°C and
initial temperature of steel = 20°C.'''
import math as m
I = 5000.0 #A (current drawn)
V = 50.0 #V (Arc Voltage)
Rs = 0.002 #ohm (transformer resistance on secondary)
Xs = 0.004 #ohm (transformer reactance on secondary)
T2 = 1370.0 #*C (final temp)
T1 = 20.0 #*C (initial temp)
# Voltage drop due to resistance =
Vr = I*Rs #V
# Voltage drop due to reactance =
Vx = I*Xs #V
#Total Voltage is (Using vector sum)
V_tot = m.sqrt((V+Vr)**2 + Vx**2) #V
#(i) Supply power factor is
pf = (V+Vr)/V_tot
#Total Power drawn => P = 3*VI*(power factor)
P = 3*V_tot*I*pf/1000 #kW
print "Power factor =",round(pf,4)
print "Power drawn from supply =",round(P,2),"kW."
#Energy required to melt 2 tonnes of steel
m = 2000.0 #kg
s = 0.12 # (specific heat of steel)
L = 8.89 #kcal/kg (latent heat of steel)
enrgy = m*s*(T2-T1) + m*L #kcal
enrgy = enrgy/860.0 #kWh
#Utilised power
P = 0.65*P #kW
#Time required for melting steel
Time = enrgy/P #hr
Time = Time*60 #min
Sec = (round(Time,2) - round(Time,-1) )*60 #sec
print "Time required for melting steel =",round(Time,-1),"minutes. and ",round(Sec),"seconds."
```

In [16]:

```
'''If a 3-phase arc furnace is to melt 10 tonne steel in 2 hours, estimate the average input to the furnace if overall
efficiency is 50%. If the current input is 9,000 A with the above kW input and the resistance and reactance of furnace leads
(including transformer) are 0.003 ohm and 0.005 ohm respectively, estimate the arc voltage and total kVA taken from the supply
Specific heat of steel = 444 J /kg/°C ,Latent heat of fusion of steel = 37.25 kJ/kg , Melting point of steel = 1,370 °C.'''
from sympy import Symbol,solve,Eq,sqrt
import math as m
mass = 10000.0 #kg (mass in kg)
t = 2.0 #hr (time taken to melt)
eff = 50.0 #% (overall efficiency)
I = 9000.0 #A (current input)
Rs = 0.003 #ohm (secondary resistance)
Xs = 0.005 #ohm (secondary reactance)
s = 444.0 #J/kg/*C (specific heat of steel)
L = 37250 #J/kg (latent heat of fusion)
T2 = 1370.0 #*C (final temp)
T1 = 20.0 #*C (initial temp)
#Energy required to melt 10 tonnes of steel
enrgy = mass*s*(T2-T1) + mass*L #J
enrgy = enrgy/(1000*3600) #kWh
#Avg output power = energy/time
P = enrgy/t #kW
#Avg input power =
Pin = P/eff*100
#Voltage drop due to resistance
Vr = I*Rs #V
#Voltage drop due to reactance
Vx = I*Xs #V
#Now,Let Va is arc drop voltage
Va = Symbol('Va') #V
Vt = sqrt((Va+Vr)**2 + Vx**2)
pf = (Va+Vr)/Vt
#Total power input = 3*(Vt*It*pf)
eq = Eq(Pin*1000,3*Vt*I*pf)
Va = solve(eq) #V
Va1 = Va[0]
Vt = sqrt((Va1+Vr)**2 + Vx**2)
#Total KVA taken from supply line =
power = 3*Vt*I/1000
print "#Total KVA taken from supply line =",round(power,2),"KVA ."
```

In [27]:

```
'''Determine the efficiency of a high-frequency induction furnace which takes 10 minutes to melt 2 kg of a aluminium initially
at a temperature of 20°C. The power drawn by the furnace is 5 kW, specific heat of aluminium = 0.212, melting point of
aluminium = 660° C and latent heat of fusion of aluminium. = 77 kcal/kg.'''
m = 2.0 #kg (mass of alluminium)
L = 77.0 #kcal/kg (Latent heat of fusion)
T2 = 660.0 #*C (final temp)
T1 = 20.0 #*C (initial temp)
s = 0.212 # (specific heat of alluminium)
Pin = 5.0 #kW (input power)
#Heat required to melt alluminium
H1 = m*L #kcal
#Heat required to raise the temperature
H2 = m*s*(T2 - T1) #kcal
#Total heat
heat_tot = H1 + H2 #kcal
#Heat required per hour
enrgy = heat_tot/(10.0/60) #kcal
#Power delivered to alluminium
Power = enrgy/860 #kW
eff = Power/Pin*100 #(% efficiency)
print "Efficiency of induction furnace = ",round(eff,2),"%."
```

In [21]:

```
'''A low-frequency induction furnace has a secondary voltage of 20V and takes 600 kW at 0.6 p.f. when the hearth is full. If the
secondary voltage is kept constant, determine the power absorbed and the p.f. when the hearth is half-full. Assume that the
resistance of the secondary circuit is doubled but the reactance remains the same.'''
import math as m
V = 20.0 #V (secondary voltage)
P = 600*1000 #W (Input Power)
pf = 0.6 # (power factor)
#Inital secondary current using P = VI*pf
I = P/(V*pf) #A (secondary current)
#Now, pf = cosQ , .'. sinQ = sqrt(1-pf^2)
Vr = V*pf #V (Voltage across resistance)
Vx = V*m.sqrt(1-pf**2) #V (Voltage across reactance)
#As, Vr = I*R and Vx = I*X
R = Vr/I #ohm
X = Vx/I #ohm
#When hearth is half-full
R2 = 2*R
X2 = X
pf = R2/m.sqrt(R2**2 + X2**2) #(new power factor)
Vr = V*pf #V (Voltage across resistance)
#As, Vr = I*R
I = Vr/R2 #A
power = V*I*pf/1000 #kW
print "Power absorbed =",round(power),"kW."
print "power factor =",round(pf,2)
```

In [22]:

```
'''Estimate the energy required to melt 0.5 tonne of brass in a single-phase induction furnace. If the melt is to be carried out
in 0.5 hour, what must be the average power input to the furnace?
Specific heat of brass = 0.094
Latent heat of fusion of brass = 39 kilocal/kg
Melting point of brass = 920°C
Furnace efficiency = 60.2%
The temperature of the cold charge may be taken as 20°C.'''
m = 0.5*1000 #kg (mass of brass)
s = 0.094 # (specific heat)
T2 = 920.0 #*C (final temp)
T1 = 20.0 #*C (initial temp)
L = 39.0 #kcal/kg(Latent heat of fusion)
eff = 60.2 # (% efficiency)
#Total amount of heat req to melt 0.5 kg brass
H = m*L + m*s*(T2-T1) #kcal
H = H/860 #kWh
H_tot = H/(eff)*100 #kWh (input energy required)
print "Total furnace input =",round(H_tot,2),"kWh ."
```

In [23]:

```
'''A low-frequency induction furnace whose secondary voltage is maintained constant at 10 V, takes 400 kW at 0.6 p.f. when the
hearth is full.Assuming the resistance of the secondary circuit to vary inversely as the height of the charge and reactance to
remain constant,find the height upto which the hearth should be filled to obtain maximum heat.'''
import math as m
V2 = 10.0 #V (secondary voltage)
P = 400.0*1000 #kW (power)
pf = 0.6 # (power factor)
#Secondary current is (Using P = VI*cosQ)
I = P/(V2*pf) #A
#Impedance of secondary circuit is
Z = V2/I #ohm
#Now, R = Z*cosQ X = Z*sinQ
R = Z*pf #ohm (resistance)
X = Z*m.sqrt(1-pf**2) #ohm (reactance)
#Let height of charge be 'x' times of the full hearth h = x*H
#Resistance varies inersely as height .Therefore,
# R' = R/x
#Now ,for max heat resistance should be equal to reactance.Therefore,
x = R/X
print "height for maximum heat = ",round(x,2),"*H."
```

In [24]:

```
'''A slab of insulating material 150 cm^2 in area and 1 cm thick is to be heated by dielectric heating. The power required is
400 W at 30 MHz.Material has relative permittivity of 5 and p.f. of 0.05. Determine the necessary voltage. Absolute
permittivity = 8.854e - 12 F/m.'''
import math as m
P = 400.0 #W (power)
f = 30.0e+6 #Hz (frequency)
A = 150.0e-4 #m^2 (area)
d = 1.0e-2 #m (thickness)
er = 5.0 # (relative permitivity)
e0 = 8.89e-12#F/m (absolute permitivity)
pf = 0.05 # (power factor)
#Capacitance
C = (A/d)*(er*e0) #F
#Now, P = (2*pi*f)*(C*V^2)*pf .Therefore V is
V = m.sqrt(P/(2*3.14*f*C*pf)) #V
print "Voltage =",round(V,2),"V."
```

In [25]:

```
'''An insulating material 2 cm thick and 200 cm^2 in area is to be heated by dielectric heating. The material has relative
permitivity of 5 and power factor of 0.05.Power required is 400 W and frequency of 40 MHz is to be used. Determine the necessary
voltage and the current that will flow through the material.If the voltage were to be limited to 700 V, what will be the
frequency to get the same loss? '''
import math as m
d = 2.0e-2 #m (Thickness)
A = 200e-4 #m^2 (Area)
er = 5 # (relative permitivity)
pf = 0.05 # (power factor)
f = 40.0e+6 #Hz (frequency)
P = 400.0 #W (power)
e0 = 8.89e-12 # (absolute permitivity)
C = (A/d)*(e0*er) #F (Capacitance)
#Now, P = 2*pi*f*C*V^2*pf
V = m.sqrt(P/(2*3.14*f*C*pf)) #V
#Also, P = VI*cosQ .Therefore,current through material
I = P/(V*pf) #A
#Heat produced is propotional to V^2*f. (V2/V1)^2 = (f1/f2)
f2 = f*(V/700)**2 #Hz
f2 = f2/(10e+5) #MHz
print "Voltage = ",round(V,2),"V."
print "Current through material =",round(I,2),"A."
print "Frequency = ",round(f2,2),"MHz."
```

In [26]:

```
'''A plywood board of 0.5*0.25*0.02 metre is to be heated from 25 to 125°C in 10 minutes by dielectric heating employing a
frequency of 30 MHz. Determine the power required in this heating process. Assume specific heat of wood 1500/J/kg/°C;
weight of wood 600 kg/m3 and efficiency of process 50%.'''
Vol = 0.5*0.25*0.02 #m^3 (Volume of plywood to be heated)
f = 30.0e+6 #Hz (frequency)
t = 10.0 #minutes (time)
T2 = 125.0 #*C (final temperature)
T1 = 25.0 #*C (initial temperature)
s = 1500.0 #J/kg/*C (specific heat of wood)
den = 600.0 #kg/m^3 (weight of wood)
eff = 50.0 #% (efficiency of process)
wt = den*Vol #kg (weight of plywood)
#Heat required to raise the temp is
H = wt*s*(T2-T1) #J
H = H/3600 #Wh
#As P = H/t .Therfore power required for heating
P = H/(10.0/60) #W
#As efficiency is 50%
Pin = P/eff*100 #W
print "Power input =",round(Pin,2),"W ."
```

In [ ]:

```
```