If you like solving for things like this, I highly encourage you to take up Chemical engineering. You’ll find good jobs too since there is a natural gas boom.
Back to the question.
Here’s what the symbols mean : V = Volume, P = Pressure, M = mass, RHO = density, N = Number of moles, MW = Molecular Weight, R = Universal Gas Constant, T = Absolute temperature
Here’s what the subscripts mean: CO2 = dry ice (assumed pure), AIR = air that fills up the balance volume of the vessel, V = Vessel, amb = Ambient temperature (assumed to be same initially and when the system reaches equilibrium), ATM = atmospheric pressure (assumed initially), FINAL = final stage after equilibrium i.e. the CO2 evaporates and the vessel temperature comes to equilibrium with ambient
So with that
1> Initially (Calculate the number of moles in the system - assume the total volume consists of the volume of solid co2 and the rest is air. Assume the air temperature is ambient and pressure is atmospheric)
V[sub]V[/sub] = V[sub]CO2[/sub] + V[sub]AIR[/sub]
=> V[sub]AIR[/sub] = V[sub]V[/sub] - V[sub]CO2[/sub] …(1)
N[sub]AIR[/sub] = (P[sub]ATM[/sub] x V[sub]AIR[/sub])/ (R x T[sub]AMB[/sub]) …(2)
Using (1), (2) reduces to
N[sub]AIR[/sub] = (P[sub]ATM[/sub]/(R x T[sub]AMB[/sub])) x ( V[sub]V[/sub] - V[sub]CO2[/sub]) …(3)
Now:
Mass[sub]CO2[/sub] = RHO[sub]CO2[/sub] x V[sub]CO2[/sub] …(4)
So moles of CO2, N[sub]CO2[/sub] = Mass[sub]CO2[/sub] / MW[sub]CO2[/sub] …(5)
Using (4), (5) reduces to
N[sub]CO2[/sub] = RHO[sub]CO2[/sub] x V[sub]CO2[/sub] / MW[sub]CO2[/sub] …(6)
Equation (3) gives the moles of air and equation (6) gives the moles of CO2
2> Now look at the final stage when it has equilibriated (assume ideal gas law although CO2 is non ideal)
P[sub]FINAL[/sub] x V[sub]V[/sub] = (N[sub]CO2[/sub] + N[sub]AIR[/sub]) x R x T[sub]AMB[/sub] … (7)
Substituting from Equation N[sub]AIR[/sub] from (3) and N[sub]CO2[/sub] from (6)
P[sub]FINAL[/sub] x V[sub]V[/sub] = (RHO[sub]CO2[/sub] x V[sub]CO2[/sub] / MW[sub]CO2[/sub] + (P[sub]ATM[/sub]/(R x T[sub]AMB[/sub])) x ( V[sub]V[/sub] - V[sub]CO2[/sub])) x R x T …(8)
Dividing both sides of equation (8) by V[sub]V[/sub] and rearranging, we have
P[sub]FINAL[/sub] = (RHO[sub]CO2[/sub]/MW[sub]CO2[/sub])x(V[sub]CO2[/sub]/V[sub]V[/sub]) x R x T[sub]AMB[/sub] + P[sub]ATM[/sub](1 - V[sub]CO2[/sub]/V[sub]V[/sub])
Now - put in the values (Assume T[sub]AMB[/sub] = 300 K or 80 F)
RHO[sub]CO2[/sub] = 1.6 g/cm3 (it can vary between 1.4 and 1.6 per wiki)
MW[sub]CO2[/sub] = 44 44.01 g/mol
V[sub]CO2[/sub] = 10 cm3, V[sub]V[/sub] = 1000 cm3
R = 82 (cm3 atm)/( K mol)
P[sub]ATM[/sub] = 1 atm
T[sub]AMB[/sub] = 300 K
Then
P[sub]FINAL[/sub] = 9.94 atm = 146 PSIA = 1007 kPa
I hope I did not make a mistake in the calcs since it is frustrating to type in code