Forse ho capito dov' era l' inghippo, per cui i conti non mi tornavano mai ..e cambiavo formule all' impazzata
l' errore grosso che facevo, prendendo seppur il funzionamento più semplice, che riepilogo, 1) immissione del gas spingendo il pistone alla pressione costante del serbatoio grande,
2) sottrazione di calore al punto morto inferiore a volume costante riportando i valori di temperatura e pressione a quelli esterni 3) compressione isoterma fino al volume v3=v1*P1/p2,
4) re-immissione del gas per una trasformazione isobara nel serbatoio..
L' errore stava proprio che re-immettendo il gas nel serbatoio l' isobara per riportarlo alle condizioni di temperatura del serbatoio sarebbe venuta di energia us=n*cp*(t2-t1)
ma non ci contavo il fatto di aver svolto un lavoro sul gas per reimmetterlo li dentro ..uguale al volume, una volta raggiunta la stessa pressione che nel serbatoio, completata l' isoterma,
chiamato v3 ..che pertanto la fiamma ad ogni giro doveva fornire una energia data dalla differenza us-v3*(p2-p1)
..che re-impostando questo conto al PC in Basiic online ..questo..
https://www.tutorialspoint.com/execute_basic_online.phpbasterà col copia incolla riportarci il file che metto sotto, dove è possibili modificare i valori di temperatura t1 e t2 (con t2>t1) di pressione p1 e di cilindrata v1 a piacimento,
per ottenere dal programma le potenze e i rendimenti premendo execute in alto a sinistra.
Il rendimento viene maggiore e allora mi propongo di provare altre combinazioni migliorandolo ancora..
P.s molti comandi del basic sul link a me non funzionano però quel che ho fatto mi pare vada ..provatelo.. facendo copia incolla così com' è questo scritto sotto
..dopo aver cancellato sulla prima linea dove trovate scritto Print "Hello World"
p1=5*(10^4)*9.81
'pressione esterna la prima cifra nella formula sono le atmosfere
t1=100+273
'temperatura esterna la prima cifra nella formula sono i gradi centigradi
t2=300+273
'temperatura nel serbatoio la prima cifra nella formula sono i gradi centigradi
ng=500/60
'numero di giri dichiarati la prima cifra sono i giri al minuto
p2=p1*t2/t1
'ricavo la pressione nel riscaldatore
print "Pressione nel serbatoio riscaldato in atmosfere =", p2/10^4/9.81
print
'stampo la pressione a gas riscaldato
v1=2000*10^-6
'volume cilindro la prima cifra sono i cm^3
'____________________________________________________________________________
'sopra questa linea ci sono i valori che potete scegliere di cambiare
'sotto il programma va avanti da se per un gas monoatomico
'____________________________________________________________________________
r=8.31
'costante dei gas perfetti
cv=3/2*r
' ricavo calore specifico a volume costante per gas monoatomici es. idrogeno
cp=cv+r
'ricavo calore specifico sempre per gas monoatomici
n=p1*v1/r/t1
'ricavo numero di moli per il volume dato in base a temperatura e pressione iniziali, ammesse uguali a quelle esterne
pd=p2-p1
'ricavo la differenza tra pressione ottenuta riscaldando il serbatoio e quella iniziale come l' esterna
e=pd*v1
'ricavo il lavoro della forza agente sul pistone fino al punto morto inferiore
ee=n*cv*(t2-t1)
'ricavo il calore dissipato raggiunto il punto morto inferiore, ammesso in un istante a volume costante
v3=v1*p1/p2
' ricavo il volume della iniziale compressione isoterma fino al raggiungimento della pressione p2 uguale a quella del sebatoio
v2=v1-v3
'ricavo il volume rimanente
w=pd*v2/2
'ricavo il lavoro necessario per la compressione isoterma che dunque coincide con il calore disperso, anch' esso però rimandato ai radiatori
w1=pd*v3
'ricavo il lavoro per ri-pompare il gas nel serbatoio
eu=e-v3*pd-w
'ricavo il lavoro utile del motore sottraendo al lavoro ottenuto in fase di espansione quello speso durante le due fasi durante la compressione
us=n*cp*(t2-t1)
'ricavo il lavoro necessario per riportare il gas alle condizioni iniziali con una trasformazione isobara
us=us-pd*v3
'ci calcolo l' energia necessaria alla fiamma sottraendoci, il lavoro apportato comprimendo il gas, per entrare nel serbatoio
print "Rendimento %=",eu/us*100
'calcolo e stampo il rendimento del motore
print "Rendimento ideale %=",(t2-t1)/t2*100
'calcolo e stampo il rendimento per il secondo principio della termodinamica
print
print "Potenza resa kW=", eu*ng/1000
'calcolo e stampo la potenza resa in kW
print "Calore fiamma KW=",us*ng/1000
'stampo il calore dato dalla fiamma in kW
print "Calore dissipato kW=", (ee+w)*ng/1000
'stampo il calore dissipato dai radiatori in kW
print
print "rapporto dei rendimenti %=", eu/us/((t2-t1)/t2)*100
'stampo il rapporto tra rendimento del modello teorico del motore e il rendimento ideale
end