GreenPolyProducts/cantera
5
0
Fork 0
mirror of https://github.com/Cantera/cantera.git synced 2025-02-25 18:55:29 -06:00
Code Issues Packages Projects Releases Wiki Activity

[oneD] Drop terms that are zero by definition

Browse Source
This commit is contained in:
Ingmar Schoegl 2023-08-13 07:20:26 -05:00 committed by Ray Speth
parent 2650acc48d
commit 349316b19f
1 changed files with 5 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
src/oneD/StFlow.cpp
View File

@ -1054,24 +1054,23 @@ void StFlow::evalContinuity(size_t j, double* x, double* rsd, int* diag, double
-(rho_u(x,j+1) - rho_u(x,j))/m_dz[j]
-(density(j+1)*V(x,j+1) + density(j)*V(x,j));
} else if (m_isFree) {
// terms involving V are zero as V=0 by definition
if (grid(j) > m_zfixed) {
rsd[index(c_offset_U,j)] =
- (rho_u(x,j) - rho_u(x,j-1))/m_dz[j-1]
- (density(j-1)*V(x,j-1) + density(j)*V(x,j));
- (rho_u(x,j) - rho_u(x,j-1))/m_dz[j-1];
} else if (grid(j) == m_zfixed) {
if (m_do_energy[j]) {
rsd[index(c_offset_U,j)] = (T(x,j) - m_tfixed);
} else {
rsd[index(c_offset_U,j)] = (rho_u(x,j)
- m_rho[0]*0.3);
- m_rho[0]*0.3); // why 0.3?
}
} else if (grid(j) < m_zfixed) {
rsd[index(c_offset_U,j)] =
- (rho_u(x,j+1) - rho_u(x,j))/m_dz[j]
- (density(j+1)*V(x,j+1) + density(j)*V(x,j));
- (rho_u(x,j+1) - rho_u(x,j))/m_dz[j];
}
} else {
// unstrained
// unstrained with fixed mass flow rate
rsd[index(c_offset_U, j)] = rho_u(x, j) - rho_u(x, j - 1);
}
}