Otto cycle mathematical model - Heat added by combustion?

AI Thread Summary
The discussion focuses on improving a spreadsheet model of the Otto cycle by updating the heat added during combustion calculations. The user is trying to replace the pressure input at ignition (P4) with a calculated value based on air-fuel ratio (AFR) and compression ratio. Two methods for calculating the heat added (Qin) are presented, with the first yielding a high ignition temperature and torque figure, while the second produces a lower ignition temperature but a more favorable torque outcome. The user is seeking guidance on which method is more accurate, as the results vary significantly with changes in AFR. The goal is to enhance the model's accuracy by incorporating additional functions and thermodynamic principles.
Kozy
Messages
140
Reaction score
7
I am attempting to improve a spreadsheet I created back in uni some 5-6 years ago, which models all the velocities and accelerations of the pistons and rods, as well calculating pressure, temperature and resulting torque. It is only a simple model so far, originally it started simply by calculating the temperature and pressure prior to ignition, then had an input for pressure at ignition (P4) from which the PV through power stroke was then calculated.

I have dug this out of my hard drive and would like to update it, replacing the P4 input with a calculated value taking into account AFR and compression ratio. So far I am having some difficulty with the 'Qin' heat added portion of this.

I have two methods for calculating this and I am unsure which, if any is correct.

Version 1:

Qin = V*1/(R*FAR)

Where:
V= V1+V2
R = 8.314
AFR = 1/AFR

This gives me 45951 of undisclosed units, I suspect joules.

Ignition temperature is then:

T4 = T3+(Qin/Cv)

Where:
T3 = 770°K
Qin = 45951
Cv = 3R/2 =12.471

Which gives me 4455°K

This method results in a torque figure with a BMEP of 291psi, however a leaner AFR does result in a hotter ignition temperature which would appear to be correct.

The second method was a bit more simple:

T4 = T3 + (HHVfuel * Mfuel)

Where:
T3 = 770°K
HHVfuel = 47.3Mj/kg
Mfuel = Mair * FAR = 4.6e-5

This resulted in a T4 of 2935°K which produced a better torque figure (BMEP of 170psi), but which become hotter with a richer AFR, which does not appear to be correct.

Can someone point me in the right direction with this? I've been chasing my tail for a while and in the absence of any thermodynamics textbooks I've not been able to find a decent answer to my problem.
 
Last edited:
Engineering news on Phys.org
Currently in the process of adding Weiber functions into the model, hopefully that should give a more realistic result.
 
Posted June 2024 - 15 years after starting this class. I have learned a whole lot. To get to the short course on making your stock car, late model, hobby stock E-mod handle, look at the index below. Read all posts on Roll Center, Jacking effect and Why does car drive straight to the wall when I gas it? Also read You really have two race cars. This will cover 90% of problems you have. Simply put, the car pushes going in and is loose coming out. You do not have enuff downforce on the right...
I'm trying to decide what size and type of galvanized steel I need for 2 cantilever extensions. The cantilever is 5 ft. The space between the two cantilever arms is a 17 ft Gap the center 7 ft of the 17 ft Gap we'll need to Bear approximately 17,000 lb spread evenly from the front of the cantilever to the back of the cantilever over 5 ft. I will put support beams across these cantilever arms to support the load evenly
Thread 'Physics of Stretch: What pressure does a band apply on a cylinder?'
Scenario 1 (figure 1) A continuous loop of elastic material is stretched around two metal bars. The top bar is attached to a load cell that reads force. The lower bar can be moved downwards to stretch the elastic material. The lower bar is moved downwards until the two bars are 1190mm apart, stretching the elastic material. The bars are 5mm thick, so the total internal loop length is 1200mm (1190mm + 5mm + 5mm). At this level of stretch, the load cell reads 45N tensile force. Key numbers...

Similar threads

Back
Top