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] Remove Flow1D::wdot

Browse Source
This commit is contained in:
Ingmar Schoegl 2023-09-11 19:01:38 -05:00 committed by Ray Speth
parent 0827444065
commit 97efeed2f3
2 changed files with 2 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
include/cantera/oneD/Flow1D.h
View File

@ -396,10 +396,6 @@ protected:
AnyMap getMeta() const override; AnyMap getMeta() const override;
void setMeta(const AnyMap& state) override; void setMeta(const AnyMap& state) override;
double wdot(size_t k, size_t j) const {
return m_wdot(k,j);
}
//! Update the properties (thermo, transport, and diffusion flux). //! Update the properties (thermo, transport, and diffusion flux).
//! This function is called in eval after the points which need //! This function is called in eval after the points which need
//! to be updated are defined. //! to be updated are defined.

4
src/oneD/Flow1D.cpp
View File

@ -569,7 +569,7 @@ void Flow1D::evalEnergy(double* x, double* rsd, int* diag,
double sum = 0.0; double sum = 0.0;
for (size_t k = 0; k < m_nsp; k++) { for (size_t k = 0; k < m_nsp; k++) {
double flxk = 0.5*(m_flux(k,j-1) + m_flux(k,j)); double flxk = 0.5*(m_flux(k,j-1) + m_flux(k,j));
sum += wdot(k,j)*m_hk(k,j); sum += m_wdot(k,j)*m_hk(k,j);
sum += flxk * m_dhk_dz(k,j) / m_wt[k]; sum += flxk * m_dhk_dz(k,j) / m_wt[k];
} }
@ -660,7 +660,7 @@ void Flow1D::evalSpecies(double* x, double* rsd, int* diag,
for (size_t k = 0; k < m_nsp; k++) { for (size_t k = 0; k < m_nsp; k++) {
double convec = rho_u(x,j)*dYdz(x,k,j); double convec = rho_u(x,j)*dYdz(x,k,j);
double diffus = 2.0*(m_flux(k,j) - m_flux(k,j-1)) / (z(j+1) - z(j-1)); double diffus = 2.0*(m_flux(k,j) - m_flux(k,j-1)) / (z(j+1) - z(j-1));
rsd[index(c_offset_Y + k, j)] = (m_wt[k]*(wdot(k,j)) rsd[index(c_offset_Y + k, j)] = (m_wt[k]*m_wdot(k,j)
- convec - diffus)/m_rho[j] - convec - diffus)/m_rho[j]
- rdt*(Y(x,k,j) - Y_prev(k,j)); - rdt*(Y(x,k,j) - Y_prev(k,j));
diag[index(c_offset_Y + k, j)] = 1; diag[index(c_offset_Y + k, j)] = 1;