Eksempel 15.6 volumarbeid_isoterm.py

# Program volumarbeid_isoterm.py

# Definerer størrelser
def p(V): # trykk [Pa]
  n = 1.00 # stoffmengde [mol]
  R = 8.31 # gasskonstant [J/(molK)]
  T = 722 # temperatur [K]
  return n*R*T/V

V1 = 0.015 # startverdi volum [m^3]
V2 = 0.030 # sluttverdi volum [m^3]
V = V1 # aktuell verdi volum [m^3]
N = 10000 # antall intervaller
dV = (V2-V1)/N # stegverdi [m^3]
W = 0 # startverdi arbeid [J]

# Summerer opp ΔW-rektanglene
while V <= V2:
  W = W + p(V)*dV # summering av arbeid
  V = V + dV # neste rektangelposisjon

print(f'W = {W:.1f} J') # utskrift med 1 desimal

W = 4158.9 J