[Transport] Update test of the power law adding missing molecular diff

parent 46f5d3bf
......@@ -151,8 +151,9 @@ void test_new_parameters (const SoluteParameters& sparam,
auto alpha = parameters.find("alpha")->second;
auto gamma = parameters.find("gamma")->second;
auto mol_diff = parameters.find("mol_diff")->second;
auto _disp = gamma*std::pow(pe,alpha);
auto _disp = mol_diff*gamma*std::pow(pe,alpha);
for (int i = 0; i < dim; ++i)
for (int j = 0; j < dim; ++j)
......@@ -162,7 +163,6 @@ void test_new_parameters (const SoluteParameters& sparam,
assert(FloatCmp::eq(disp[i][j],0.));
auto char_len = parameters.find("char_length")->second;
auto mol_diff = parameters.find("mol_diff")->second;
assert(FloatCmp::eq(pe,char_len*fluxn/(mol_diff*wc)));
}
else if (name == "material_2") // const diff + const hm disp
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment