Hi Jake,
I took another look at your problem and I think it is fundamentally ill-posed. Your reaction equations need to be linearly independent. You have 5 reactions, but only 3 are independent. You can check my reasoning below.
I am assuming that I can extract a stoichiometric matrix,
$\nu$, from you combustion energy equations or
(*combustion energy*)
DELCP1[T_] := CPB[T] + CPD[T] - CPA[T] - 0.5*CPC[T]
DELCP2[T_] := 2*CPE[T] + 3*CPD[T] - CPA[T] - 3.5*CPC[T]
DELCP3[T_] := 2*CPF[T] + 3*CPD[T] - CPA[T] - 2.5*CPC[T]
DELCP4[T_] := 2*CPE[T] + 2*CPD[T] - CPB[T] - 3*CPC[T]
DELCP5[T_] := 2*CPF[T] + 2*CPD[T] - CPB[T] - 2*CPC[T]
This should lead to the following
$\nu$ matrix
nu = {
{-1, 1 , -1/2, 1, 0, 0},
{-1, 0, -7/2, 3, 2, 0},
{-1, 0, -5/2, 3, 0, 2},
{0, -1, -3, 2, 2, 0},
{0, -1, -2, 2, 0, 2}
};
The function MatrixRank[] returns the number of linearly independent rows or columns. When applied to
$\nu$, we obtain
In[7]:= MatrixRank[nu]
Out[7]= 3
Therefore, only three of you reactions are independent and you have over specified your system with 5 reactions and over constrained systems are not solvable.
