flux BC tweak

This commit is contained in:
James E McClure 2018-10-16 11:31:18 -04:00
parent 89aa7e3c87
commit ceb85c2cd0

View File

@ -80,7 +80,7 @@ void ScaLBL_ColorModel::ReadParams(string filename){
outletA=0.f;
outletB=1.f;
if (BoundaryCondition==4) flux = din*rhoA; // mass flux must adjust for density (see formulation for details)
if (BoundaryCondition==4) flux *= rhoA; // mass flux must adjust for density (see formulation for details)
// Read domain parameters
auto L = domain_db->getVector<double>( "L" );