Rheology_Arrhenius.py

import os
import modena
import SurfaceTension, polymerViscosity
import Rheology
from modena import ForwardMappingModel, BackwardMappingModel, SurrogateModel, CFunction, ModenaFireTask
import modena.Strategy as Strategy
from fireworks.user_objects.firetasks.script_task import FireTaskBase, ScriptTask
from fireworks import Firework, Workflow, FWAction
from fireworks.utilities.fw_utilities import explicit_serialize
from jinja2 import Template

f = CFunction(
    Ccode= '''
#include "modena.h"
#include "math.h"
#include <stdio.h>

void Arrhenius
(
    const modena_model_t* model,
    const double* inputs,
    double *outputs
)
{
    {% block variables %}{% endblock %}

    const double A_mu = 0.0387; // could be loaded from some library for consistence
    const double E_mu = 10000;
    const double R_rh = 8.314;

    double f_t, mu_ap;

    f_t = A_mu * exp(E_mu / R_rh / T );

    mu_ap = mu_car * f_t;
    // printf("apparent viscosity %f", mu_ap);
    outputs[0] = mu_ap;
}
''',
   inputs={
       'T': {'min': 0, 'max': 9e99 },
       'shear': {'min': 0, 'max': 9e99 },
       'X': {'min': 0, 'max': 1 },
       'm0' : {'min': 0, 'max' : 9e99},
       'm1' : {'min': 0, 'max' : 9e99},
       'mu' : {'min': 0, 'max' : 9e99},
       'ST' : {'min': 0, 'max' : 9e99},
       'mu_car' : {'min': 0, 'max' : 9e99},
   },
   outputs={
       'mu_ap': { 'min': 0, 'max': 9e99, 'argPos': 0 },
   },
   parameters={
       'Rgas' : {'min': 8.31, 'max': 8.32, 'argPos': 0},
   }
)


m = ForwardMappingModel(
    _id='Rheology_Arrhenius',
    surrogateFunction=f,
    substituteModels=[ Rheology.m ],
    parameters=[ 8.314 ],
)