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

Correctly size grad_lnAC so that we don't write past the end of it

Browse source 
for 2D and 3D problems.

Fix the creation of the Stefan-Maxwell matrix coefficients.
This commit is contained in:
Victor Brunini 2012-11-26 20:20:52 +00:00
parent c3233c1355
commit 985b1f07a8
1 changed files with 7 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

14
src/transport/LiquidTransport.cpp
View file

@ -416,7 +416,7 @@ bool LiquidTransport::initLiquid(LiquidTransportParams& tr)
m_chargeSpecies[i] = m_thermo->charge(i);
}
m_volume_spec.resize(m_nsp, 0.0);
m_Grad_lnAC.resize(m_nsp, 0.0);
m_Grad_lnAC.resize(m_nDim * m_nsp, 0.0);
m_spwork.resize(m_nsp, 0.0);
// resize the internal gradient variables
@ -1692,9 +1692,9 @@ void LiquidTransport::stefan_maxwell_solve()
//if ( !( m_bdiff(i,j) > 0.0 ) )
//throw CanteraError("LiquidTransport::stefan_maxwell_solve",
// "m_bdiff has zero entry in non-diagonal.");
tmp = m_molefracs_tran[j] * m_bdiff(i,j);
tmp = m_molefracs_tran[i] * m_molefracs_tran[j] * m_bdiff(i,j);
m_A(i,i) -= tmp;
m_A(i,j) = + tmp;
m_A(i,j) = tmp;
}
}
}
@ -1756,9 +1756,9 @@ void LiquidTransport::stefan_maxwell_solve()
//if ( !( m_bdiff(i,j) > 0.0 ) )
//throw CanteraError("LiquidTransport::stefan_maxwell_solve",
// "m_bdiff has zero entry in non-diagonal.");
tmp = m_molefracs_tran[j] * m_bdiff(i,j);
tmp = m_molefracs_tran[i] * m_molefracs_tran[j] * m_bdiff(i,j);
m_A(i,i) -= tmp;
m_A(i,j) = + tmp;
m_A(i,j) = tmp;
}
}
}
@ -1801,9 +1801,9 @@ void LiquidTransport::stefan_maxwell_solve()
//if ( !( m_bdiff(i,j) > 0.0 ) )
//throw CanteraError("LiquidTransport::stefan_maxwell_solve",
// "m_bdiff has zero entry in non-diagonal.");
tmp = m_molefracs_tran[j] * m_bdiff(i,j);
tmp = m_molefracs_tran[i] * m_molefracs_tran[j] * m_bdiff(i,j);
m_A(i,i) -= tmp;
m_A(i,j) = + tmp;
m_A(i,j) = tmp;
}
}
}