OilfieldToolsNet/LBPM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

fix dumb bugs

Browse Source
This commit is contained in:
Zhe Li 2022-04-13 14:53:23 +10:00
parent 0f68c118de
commit c281c0cfda
2 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
models/PoissonSolver.cpp
View File

@ -416,7 +416,7 @@ void ScaLBL_Poisson::Potential_Init(double *psi_init){
auto AffinityList = electric_db->getVector<double>( "InitialValues" ); auto AffinityList = electric_db->getVector<double>( "InitialValues" );
size_t NLABELS = LabelList.size(); size_t NLABELS = LabelList.size();
if (NLABELS != AffinityList.size() || NLABELS != BoundaryConditionSolidList.size()){ if (NLABELS != AffinityList.size()){
ERROR("Error: LB-Poisson Solver: InitialValueLabels and InitialValues must be of the same length! \n"); ERROR("Error: LB-Poisson Solver: InitialValueLabels and InitialValues must be of the same length! \n");
} }

3
tests/lbpm_nernst_planck_simulator.cpp
View File

@ -101,6 +101,9 @@ int main(int argc, char **argv)
timestep++; timestep++;
PoissonSolver.Run(IonModel.ChargeDensity,false,timestep);//solve Poisson equtaion to get steady-state electrical potental PoissonSolver.Run(IonModel.ChargeDensity,false,timestep);//solve Poisson equtaion to get steady-state electrical potental
//StokesModel.Run_Lite(IonModel.ChargeDensity, PoissonSolver.ElectricField);// Solve the N-S equations to get velocity //StokesModel.Run_Lite(IonModel.ChargeDensity, PoissonSolver.ElectricField);// Solve the N-S equations to get velocity
//for debugging
IonModel.getIonConcentration_debug(0);
IonModel.Run(IonModel.FluidVelocityDummy,PoissonSolver.ElectricField); //solve for ion transport and electric potential IonModel.Run(IonModel.FluidVelocityDummy,PoissonSolver.ElectricField); //solve for ion transport and electric potential