Mathcad - Example 2 - 2 Finite Heat Release
Mathcad - Example 2 - 2 Finite Heat Release
Mathcad - Example 2 - 2 Finite Heat Release
Given:
Asked :
(a) the net work, the thermal efficiency of the cycle, and the mean effective pressure, and
(b) the pressure, temperature and work profiles versus crank angle. Assume that the ideal gas
specific heat ratio is 1.4, the gas molecular weight is 29, the combustion duration is constant at
40o , and that the Weibe parameters are a=5 and n=3.
Solution:
a)
Numerical Solution:
γ− 1
ηt := 1 − 1 = 0.602 R := 2
L
r S
c
m
Up := 2 ⋅ S⋅ N = 62.832
s
π 2
Vd := ⋅ B ⋅ S = 0.785 L
4
Vd Vd cos( 2 ⋅ θ)
⋅ R + 1 − cos( θ) − R − 1 −
2
V( θ) := +
rc − 1 2 2
−3
8× 10
−3
6× 10
V( θ) −3
4× 10
L
−3
2× 10
0
− 100 0 100
θ
deg
0 otherwise 0 otherwise
Heat Release Example Quinn Straub
3 2
2× 10
3 1.5
1.5× 10
Q( θ) 3 xb( θ) 1
1× 10
J
500 0.5
0 0
− 100 0 100 − 100 0 100
θ θ
deg deg
Boundary Conditions:
5
θ0 := −180 deg P0 := P1 = 1 × 10 Pa θf := 180deg n := 128
θ := for i ∈ 0 .. n
(θf − θ0)
M ← θ0 + i ⋅
i n
M
P d γ − 1 d P d γ − 1 d
f ( θ , P) := −γ⋅ ⋅
V( θ) + ⋅ Q( θ) f2 ( θ , P) := −γ⋅ ⋅ V( θ) + ⋅ Q2 ( θ)
V( θ) dθ V( θ) dθ V( θ) dθ V( θ) dθ
(
k1 ← f θ0 + i ⋅ h , M
i) (
k1 ← f2 θ0 + i ⋅ h , N
i)
k2 ← f ( θ0 + i ⋅ h + 0.5⋅ h , M + 0.5⋅ k1⋅ h ) k2 ← f2 ( θ0 + i ⋅ h + 0.5⋅ h , N + 0.5⋅ k1⋅ h )
i i
k3 ← f ( θ0 + i ⋅ h + 0.5⋅ h , M + 0.5⋅ k2⋅ h ) k3 ← f2 ( θ0 + i ⋅ h + 0.5⋅ h , N + 0.5⋅ k2⋅ h )
i i
k4 ← f ( θ0 + i ⋅ h + h , M + k3⋅ h ) k4 ← f2 ( θ0 + i ⋅ h + h , N + k3⋅ h )
i i
1 1
M ←M + ( k1 + 2 ⋅ k2 + 2 ⋅ k3 + k4) ⋅ h N ←N + ( k1 + 2 ⋅ k2 + 2 ⋅ k3 + k4) ⋅ h
i+ 1 i 6 i+ 1 i 6
M N
Heat Release Example Quinn Straub
3
8× 10
3
6× 10
P 3
4× 10 Engine 1
kPa
3
2× 10
0
− 200 − 100 0 100 200
θ
deg
3
8× 10
3
6× 10
Engine 1
P Engine 2
kPa
3
4× 10
P2
kPa
3
2× 10
0
− 200 − 100 0 100 200
θ
deg
The solution of just the Pressure vs. Crank Angle is shown. The question ask for other parameters
which can be solved for given the necessary equations.
Heat Release Example Quinn Straub
d
Heat Release Example Quinn Straub