void
AddCoupledSolidKinSpeciesAuxKernelsAction::act()
{
std::vector<std::string> reactions = getParam<std::vector<std::string> >("kin_reactions");
std::vector<Real> logk = getParam<std::vector<Real> >("log10_keq");
std::vector<Real> r_area = getParam<std::vector<Real> >("specific_reactive_surface_area");
std::vector<Real> ref_kconst = getParam<std::vector<Real> >("kinetic_rate_constant");
std::vector<Real> e_act = getParam<std::vector<Real> >("activation_energy");
Real gas_const = getParam<Real>("gas_constant");
std::vector<Real> ref_temp = getParam<std::vector<Real> >("reference_temperature");
std::vector<Real> sys_temp = getParam<std::vector<Real> >("system_temperature");
// NEED TO ADD AN ERROR MESSAGE IF THE SIZES OF ABOVE ARRAYS ARE NOT THE SAME //
for (unsigned int j=0; j < reactions.size(); j++)
{
std::vector<std::string> tokens;
std::vector<std::string> solid_kin_species(reactions.size());
// Parsing each reaction
MooseUtils::tokenize(reactions[j], tokens, 1, "+=");
std::vector<Real> stos(tokens.size()-1);
std::vector<VariableName> rxn_vars(tokens.size()-1);
for (unsigned int k=0; k < tokens.size(); k++)
{
Moose::out << tokens[k] << "\t";
std::vector<std::string> stos_vars;
MooseUtils::tokenize(tokens[k], stos_vars, 1, "()");
if (stos_vars.size() == 2)
{
Real coef;
std::istringstream iss(stos_vars[0]);
iss >> coef;
stos[k] = coef;
rxn_vars[k] = stos_vars[1];
}
else
{
solid_kin_species[j] = stos_vars[0];
}
}
void
AddCoupledSolidKinSpeciesKernelsAction::act()
{
std::vector<NonlinearVariableName> vars = getParam<std::vector<NonlinearVariableName> >("primary_species");
std::vector<std::string> reactions = getParam<std::vector<std::string> >("kin_reactions");
Moose::out<< "Solid kinetic reaction list:" << "\n";
for (unsigned int i=0; i < reactions.size(); i++)
{
Moose::out<< reactions[i] << "\n";
}
for (unsigned int i=0; i < vars.size(); i++)
{
Moose::out << "primary species - " << vars[i] << "\n";
std::vector<bool> primary_participation(reactions.size(), false);
std::vector<std::string> solid_kin_species(reactions.size());
std::vector<Real> weight;
for (unsigned int j=0; j < reactions.size(); j++)
{
std::vector<std::string> tokens;
// Parsing each reaction
MooseUtils::tokenize(reactions[j], tokens, 1, "+=");
std::vector<Real> stos(tokens.size()-1);
std::vector<std::string> rxn_vars(tokens.size()-1);
for (unsigned int k=0; k < tokens.size(); k++)
{
Moose::out << tokens[k] << "\t";
std::vector<std::string> stos_vars;
MooseUtils::tokenize(tokens[k], stos_vars, 1, "()");
if (stos_vars.size() == 2)
{
Real coef;
std::istringstream iss(stos_vars[0]);
iss >> coef;
stos[k] = coef;
rxn_vars[k] = stos_vars[1];
Moose::out << "stochiometric: " << stos[k] << "\t";
Moose::out << "reactant: " << rxn_vars[k] << "\n";
// Check the participation of primary species
if (rxn_vars[k] == vars[i]) primary_participation[j] = true;
}
else
{
solid_kin_species[j] = stos_vars[0];
}
}
// Done parsing, recorded stochiometric and variables into separate arrays
Moose::out << "whether primary present (0 is not): " << primary_participation[j] << "\n";
if (primary_participation[j])
{
// Assigning the stochiometrics based on parsing
for (unsigned int m=0; m < rxn_vars.size(); m++)
{
if (rxn_vars[m] == vars[i])
{
weight.push_back(stos[m]);
Moose::out << "weight for " << rxn_vars[m] <<" : " << weight[weight.size()-1] << "\n";
}
}
Moose::out << "solid kinetic species: " << solid_kin_species[j] << "\n";
std::vector<VariableName> coupled_var(1);
coupled_var[0] = solid_kin_species[j];
// Building kernels for solid kinetic species
InputParameters params_kin = _factory.getValidParams("CoupledBEKinetic");
params_kin.set<NonlinearVariableName>("variable") = vars[i];
params_kin.set<std::vector<Real> >("weight") = weight;
params_kin.set<std::vector<VariableName> >("v") = coupled_var;
_problem->addKernel("CoupledBEKinetic", vars[i]+"_"+solid_kin_species[j]+"_kin", params_kin);
Moose::out << vars[i]+"_"+solid_kin_species[j]+"_kin" << "\n";
params_kin.print();
}
}