research/cantera
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

[Samples/Matlab] Fix bugs in conhp and conuv Matlab samples

Browse source
This commit is contained in:
Bryan W. Weber 2014-07-18 14:35:22 +00:00
parent 6cd78c0233
commit a067fcbc37
2 changed files with 5 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
samples/matlab/conhp.m
View file

@ -1,4 +1,4 @@
function dydt = conhp(t,y,gas,mw) %#ok<INUSL>
function dydt = conhp(t, y, gas, mw) %#ok<INUSL>
% CONHP ODE system for a constant-pressure, adiabatic reactor.
%
% Function CONHP evaluates the system of ordinary differential
@ -10,6 +10,7 @@ function dydt = conhp(t,y,gas,mw) %#ok<INUSL>
% Set the state of the gas, based on the current solution vector.
setMassFractions(gas, y(2:end), 'nonorm');
set(gas, 'T', y(1), 'P', pressure(gas));
nsp = nSpecies(gas);
% energy equation
wdot = netProdRates(gas);
@ -18,11 +19,10 @@ tdot = - temperature(gas) * gasconstant * enthalpies_RT(gas)' ...
% set up column vector for dydt
dydt = [ tdot
zeros(53,1) ];
zeros(nsp, 1) ];
% species equations
rrho = 1.0/density(gas);
nsp = nSpecies(gas);
for i = 1:nsp
dydt(i+1) = rrho*mw(i)*wdot(i);
end

4
samples/matlab/conuv.m
View file

@ -1,4 +1,4 @@
function dydt = conuv(t,y,gas,mw) %#ok<INUSL>
function dydt = conuv(t, y, gas, mw) %#ok<INUSL>
% CONUV ODE system for a constant-volume, adiabatic reactor.
%
% Function CONUV evaluates the system of ordinary differential
@ -19,7 +19,7 @@ tdot = - temperature(gas) * gasconstant * (enthalpies_RT(gas) - ones(nsp,1))' ..
% set up column vector for dydt
dydt = [ tdot
zeros(53,1) ];
zeros(nsp, 1) ];
% species equations
rrho = 1.0/density(gas);