dv_chi_dmf.cc

Go to the documentation of this file.

 
#include "dv_chi_dmf.h"
#include "domain.h"
#include <cstdlib>
#include <cmath>
 
 
dv_chi_dmf::dv_chi_dmf(domain  *line,
                 const         string s,
                 const bool    Lt,
                 const bool    Lo) {
 
    domn          = line;
    var_name      = s;
    L_transported = Lt;
    L_output      = Lo;
    d             = vector<double>(domn->ngrd, 0.0);
 
    Dmf           = domn->io->streamProps["Dmf"] ? domn->io->streamProps["Dmf"].as<double>() : 0.0;
    if(Dmf == 0.0) {
        cout << endl << "ERROR: if you are defining dv_chi_dmf for chi, you need to set Dmf";
        exit(0);
    }
 
}
 
 
void dv_chi_dmf::setVar(const int ipt){
 
    if(ipt != -1) {
        cout << endl << "ERROR in setVar: ipt must be -1" << endl;
        exit(0);
    }
 
    d.resize(domn->ngrd);
 
    domn->mixf->setVar();     // this may be redundant
 
    vector<double> gradZ(domn->ngrd);
 
    //------------- Compute chi
 
    gradZ.at(0) = (domn->mixf->d.at(1)-domn->mixf->d.at(0))/(domn->pos->d.at(1)-domn->pos->d.at(0));
    double d1, d2;
    for(int i=1; i<domn->ngrd-1; i++) {
        d1 = 0.5*(domn->pos->d.at(i)  -domn->pos->d.at(i-1));
        d2 = 0.5*(domn->pos->d.at(i+1)-domn->pos->d.at(i));
        gradZ.at(i) = (d2*d2*(domn->mixf->d.at(i)  -domn->mixf->d.at(i-1)) +
                    d1*d1*(domn->mixf->d.at(i+1)-domn->mixf->d.at(i)))/(2*d1*d2*(d1+d2));
    }
    gradZ.at(domn->ngrd-1) = (domn->mixf->d.at(domn->ngrd-1)-domn->mixf->d.at(domn->ngrd-2)) /
                          (domn->pos->d.at(domn->ngrd-1)-domn->pos->d.at(domn->ngrd-2));
 
    for(int i=0; i<domn->ngrd; i++)
        d.at(i) = 2.0*Dmf*gradZ.at(i)*gradZ.at(i);
 
}