polymerViscosity.py

import os
import json
from modena import *


f_polymerViscosity = CFunction(
    Ccode='''
#include "modena.h"
#include "math.h"

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

    const double Aeta = parameters[0];
    const double Eeta = parameters[1];
    const double AA = parameters[2];
    const double B = parameters[3];
    const double Xg = parameters[4];

    const double Rg = 8.31446218;

    if (X<Xg) {
        outputs[0] = Aeta*exp(Eeta/(Rg*T))*pow(Xg/(Xg-X),AA+B*X);
    } else {
        outputs[0] = 1e10;
    }
}
''',
    # These are global bounds for the function
    inputs={
        'T': {'min': 200, 'max': 550 },
        'X': {'min': 0, 'max': 1 },
    },
    outputs={
        'mu': {'min': 0, 'max': +9e99, 'argPos': 0},
    },
    parameters={
        'param1': {'min': -1e9, 'max': 1e9, 'argPos': 0},
        'param2': {'min': -1e9, 'max': 1e9, 'argPos': 1},
        'param3': {'min': -1e9, 'max': 1e9, 'argPos': 2},
        'param4': {'min': -1e9, 'max': 1e9, 'argPos': 3},
        'param5': {'min': -1e9, 'max': 1e9, 'argPos': 4},
    },
)
with open(os.getcwd()+'/inputs/unifiedInput.json') as jsonfile:
    
    inputs=json.load(jsonfile)
    
    X_gel=inputs['kinetics']['gelPoint']


par = [4.1e-8, 38.3e3, 4.0, -2.0, X_gel]


par2 = [10.3e-8, 41.3e3, 1.5, 1.0, X_gel]


par3 = [3.32e-8, 42.9e3, 2.32, 1.4, X_gel]


par4 = [4.1e-8, 38.3e3, 4.0, -2.0, X_gel]


par5 = [3.1e0, 2.24e3, 3.5, -2.0, X_gel]


par6 = [1.6e-7, 44.9e3, 1.29, 1.86, X_gel]


par7 = [2.49707805e-13, 6.88758966e+04, 1.03553549e+01, -1.46554648e+01, 0.5]


par8 = [7.03759789e-17, 8.91382466e+04, 1.32209736e+01, -1.82153953e+01, 0.5]


m_polymerViscosity = ForwardMappingModel(
    _id='polymerViscosity',
    surrogateFunction=f_polymerViscosity,
    substituteModels=[],
    parameters=par,
)