SystemDisplay.cpp

Go to the documentation of this file.

// SystemDisplay.cpp: implementation of the SystemDisplay class.
//
 
#include "SystemDisplay.h"
#include "../System/SystemModule.h"
#include "../Branch/BranchRecord.h"
#include <boost/serialization/export.hpp>
#include "../Global/SimulationContext.h"
#include <boost/archive/binary_oarchive.hpp>
#include <boost/archive/binary_iarchive.hpp>
BOOST_CLASS_EXPORT_IMPLEMENT(SystemDisplay)
using namespace std;
 
 
// Construction/Destruction
 
SystemDisplay::SystemDisplay()
{
    functional_state = sim::State::Succeed;
    corresponding_system_parameter_name_in_number = -1;
    ptr_to_corresponding_TH_element = 0;
}
 
SystemDisplay::~SystemDisplay()
{
 
}
void SystemDisplay::setTypeOfDisplay(DisplayType n){type_of_display = n;}
DisplayType SystemDisplay::getTypeOfDisplay(){return type_of_display;}
void SystemDisplay::setFunctionalState(int n) {functional_state = n;}
int SystemDisplay::getFunctionalState(){return functional_state;}
void SystemDisplay::readInfoFromFile(ifstream &fin){
    DisplayElement::readInfoFromFile(fin);
}
void SystemDisplay::readInfoFromFileWithoutDisplayName(ifstream &fin){
    DisplayElement::readInfoFromFileWithoutDisplayName(fin);
}
double SystemDisplay::getCorrespondingHeatStructureValue(int _sequence_of_mesh_point){
    HeatStructure *hs_ptr;
    double r_value;
    string err_str;
    r_value = -99999.0;
//
    if(ptr_to_corresponding_TH_element == 0){
        ptr_to_corresponding_TH_element = g_sim.Plant.getComponentFromSystemMap(SystemDisplay::getCorrespondingSystemName());
    }
    if(ptr_to_corresponding_TH_element != 0){
        switch(ptr_to_corresponding_TH_element->getTypeOfCard()){
        case Relap5ComponentType::HeatStructure:
            hs_ptr = (HeatStructure *) ptr_to_corresponding_TH_element;
            r_value = hs_ptr->getMeshPointTemperature(_sequence_of_mesh_point);
            break;
        default:
            cout << "SystemDisplay::getCorrespondingSystemParameterValue Heat Structure has incorrect type of card: "
                << static_cast<int>(ptr_to_corresponding_TH_element->getTypeOfCard()) << endl;
            ofstream fout("err messages.txt", ios::app);
            fout << r5pars_.timesno << "  SystemDisplay::getCorrespondingSystemParameterValue Heat Structure has incorrect type of card: "
                << static_cast<int>(ptr_to_corresponding_TH_element->getTypeOfCard()) << '\n';
            fout.close();
            r5pars_.sim_err = true;
            g_sim.Branch_center.setTempESMessage("SystemDisplay::getCorrespondingHeatStructureValue err.");
            g_sim.Branch_center.setTempTypeOfConsequence(Consequence::AdsSimErr);
            break;
        }
    }
    else{//could not find the corresponding component
        cout << "SystemDisplay::getCorrespondingSystemParameterValue Heat Structure can not find the corresponding component: "
            << SystemDisplay::getCorrespondingSystemName() << endl;
        ofstream fout("err messages.txt", ios::app);
        fout << r5pars_.timesno << "  SystemDisplay::getCorrespondingSystemParameterValue Heat Structure can not find the corresponding component: "
            << SystemDisplay::getCorrespondingSystemName() << '\n';
        fout.close();
        r5pars_.sim_err = true;
        g_sim.Branch_center.setTempESMessage("SystemDisplay::getCorrespondingHeatStructureValue err.");
        g_sim.Branch_center.setTempTypeOfConsequence(Consequence::AdsSimErr);
    }
    return r_value;
}
 
double SystemDisplay::getCorrespondingSystemParameterValue(){
    ControlVariable *cv_ptr;
//  HeatStructure *hs_ptr;
    HydraulicJunction *hj_ptr;
    HydraulicVolume *hv_ptr;
    InteractiveVariable *iv_ptr;
    LogicalTrip *lt_ptr;
    VariableTrip *vt_ptr;
    double r_value;
    string err_str, str;
    r_value = -99999.0;
//  map<string, THElement *> system_map; 
    if(ptr_to_corresponding_TH_element == 0){
        ptr_to_corresponding_TH_element = g_sim.Plant.getComponentFromSystemMap(SystemDisplay::getCorrespondingSystemName());
    }
    if(ptr_to_corresponding_TH_element != 0){
        switch(ptr_to_corresponding_TH_element->getTypeOfCard()){
        case Relap5ComponentType::HydraulicVolume:
            str = SystemDisplay::getCorrespondingSystemParameterName();
            if (corresponding_system_parameter_name_in_number == -1){
                if(str.compare("Pressure") == 0){
                    corresponding_system_parameter_name_in_number = 0;
                }
                else if(str.compare("Liquid_Temperature") == 0){
                    corresponding_system_parameter_name_in_number = 1;
                }
                else if(str.compare("Vapor_Temperature") == 0){
                    corresponding_system_parameter_name_in_number = 2;
                }
                else if(str.compare("Saturation_Temperature") == 0){
                    corresponding_system_parameter_name_in_number = 3;
                }
                else if(str.compare("Boron_Density") == 0){
                    corresponding_system_parameter_name_in_number = 4;
                }
                else if(str.compare("Equilibrium_Quality") == 0){
                    corresponding_system_parameter_name_in_number = 5;
                }
                else if(str.compare("Flow_Regime_Number") == 0){
                    corresponding_system_parameter_name_in_number = 6;
                }
                else {
                    cout << "SystemDisplay::getCorrespondingSystemParameterValue (SystemDisplay::getCorrespondingSystemName()): "
                        << SystemDisplay::getCorrespondingSystemName() << " can not be recognized" << endl;
                    ofstream fout("err messages.txt", ios::app);
                    fout << r5pars_.timesno << "  SystemDisplay::getCorrespondingSystemParameterValue (SystemDisplay::getCorrespondingSystemName()): "
                        << SystemDisplay::getCorrespondingSystemName() << " can not be recognized\n";
                    fout.close();
                    r5pars_.sim_err = true;
                    g_sim.Branch_center.setTempESMessage("SystemDisplay::getCorrespondingSystemParameterValue err.");
                    g_sim.Branch_center.setTempTypeOfConsequence(Consequence::AdsSimErr);
                }
            }
            if(!r5pars_.sim_err){
                hv_ptr = (HydraulicVolume *) ptr_to_corresponding_TH_element;
                switch (corresponding_system_parameter_name_in_number){
                case 0:
                        r_value = hv_ptr->getPressure();
                    break;
                case 1:
                        r_value = hv_ptr->getLiquidTemperature();
                    break;
                case 2:
                    r_value = hv_ptr->getVaporTemperature();
                    break;
                case 3:
                    r_value = hv_ptr->getSaturationTemperature();
                    break;
                case 4:
                    r_value = hv_ptr->getBoronDensity();
                    break;
                case 5:
                    r_value = hv_ptr->getEquilibriumQuality();
                    break;
                case 6:
                    r_value = hv_ptr->getFlowRegimeNumber();
                    break;
                default:
                    break;
                }
            }
            break;
        case Relap5ComponentType::HydraulicJunction:
            if (corresponding_system_parameter_name_in_number == -1){
                str = SystemDisplay::getCorrespondingSystemParameterName(); 
                if(str.compare("Mass_Flow_Rate") == 0){
                    corresponding_system_parameter_name_in_number = 0;
                }
                else if(str.compare("Liquid_Velocity") == 0){
                    corresponding_system_parameter_name_in_number = 1;
                }
                else if(str.compare("Vapor_Velocity") == 0){
                    corresponding_system_parameter_name_in_number = 2;
                }
                else {
                    cout << "SystemDisplay::getCorrespondingSystemParameterValue (SystemDisplay::getCorrespondingSystemName()): "
                        << (SystemDisplay::getCorrespondingSystemName()) << " can not be recognized" << endl;
                    ofstream fout("err messages.txt", ios::app);
                    fout << r5pars_.timesno << "  SystemDisplay::getCorrespondingSystemParameterValue (SystemDisplay::getCorrespondingSystemName()): "
                        << (SystemDisplay::getCorrespondingSystemName()) << " can not be recognized\n";
                    fout.close();
                    r5pars_.sim_err = true;
                    g_sim.Branch_center.setTempESMessage("SystemDisplay::getCorrespondingSystemParameterValue err.");
                    g_sim.Branch_center.setTempTypeOfConsequence(Consequence::AdsSimErr);
                }
            }
            if(!r5pars_.sim_err){
                hj_ptr = (HydraulicJunction *) ptr_to_corresponding_TH_element;
                switch (corresponding_system_parameter_name_in_number){
                case 0:
                        r_value = hj_ptr->getMassFlowRate();
                    break;
                case 1:
                        r_value = hj_ptr->getLiquidVelocity();
                    break;
                case 2:
                    r_value = hj_ptr->getVaporVelocity();
                    break;
                default:
                    break;
                }   
            }
            break;
        case Relap5ComponentType::VariableTrip: 
            vt_ptr = (VariableTrip*) ptr_to_corresponding_TH_element;
            r_value = vt_ptr->getTripTime();
            break;
        case Relap5ComponentType::LogicalTrip:
            lt_ptr = (LogicalTrip*) ptr_to_corresponding_TH_element;
            r_value = lt_ptr->getTripTime();
            break;
        case Relap5ComponentType::InteractiveVariable:
            iv_ptr = (InteractiveVariable *) ptr_to_corresponding_TH_element;
            r_value = iv_ptr->getDisplayValue();
            break;
        case Relap5ComponentType::ControlVariable:
            cv_ptr = (ControlVariable *) ptr_to_corresponding_TH_element;
            r_value = cv_ptr->getControlValue();
            break;
//      case Relap5ComponentType::HeatStructure:
//          hs_ptr = (HeatStructure *) ptr_to_corresponding_TH_element;
//          r_value = hs_ptr->getMeshPointTemperature(sequence_of_mesh_point);
//          break;
        default:
            cout << "SystemDisplay::getCorrespondingSystemParameterValue can not find the corresponding type of component: "
                << SystemDisplay::getCorrespondingSystemName() << endl;
            ofstream fout("err messages.txt", ios::app);
            fout << r5pars_.timesno << "  SystemDisplay::getCorrespondingSystemParameterValue can not find the corresponding type of component: "
                << SystemDisplay::getCorrespondingSystemName() << '\n';
            fout.close();
            r5pars_.sim_err = true;
            g_sim.Branch_center.setTempESMessage("SystemDisplay::getCorrespondingSystemParameterValue err.");
            g_sim.Branch_center.setTempTypeOfConsequence(Consequence::AdsSimErr);
            break;
        }
    }
    else{//could not find the corresponding component
        cout << "SystemDisplay::getCorrespondingSystemParameterValue can not find the corresponding component: "
            << SystemDisplay::getCorrespondingSystemName() << endl;
        ofstream fout("err messages.txt", ios::app);
        fout << r5pars_.timesno << "  SystemDisplay::getCorrespondingSystemParameterValue can not find the corresponding component: "
            << SystemDisplay::getCorrespondingSystemName() << '\n';
        fout.close();
        r5pars_.sim_err = true;
        g_sim.Branch_center.setTempESMessage("SystemDisplay::getCorrespondingSystemParameterValue err.");
        g_sim.Branch_center.setTempTypeOfConsequence(Consequence::AdsSimErr);
    }
    return r_value;
}
 
void SystemDisplay::setCorrespondingSystemParameterUnit(){
    ControlVariable *cv_ptr;
    HeatStructure *hs_ptr;
    HydraulicJunction *hj_ptr;
    HydraulicVolume *hv_ptr;
    string r_str = "NOT_SPECIFIED";
    string str;
//  map<string, THElement *> system_map; 
    if(ptr_to_corresponding_TH_element == 0){
        ptr_to_corresponding_TH_element = g_sim.Plant.getComponentFromSystemMap(SystemDisplay::getCorrespondingSystemName());
    }
    if(ptr_to_corresponding_TH_element != 0){
        switch(ptr_to_corresponding_TH_element->getTypeOfCard()){
        case Relap5ComponentType::HydraulicVolume:
            if (corresponding_system_parameter_name_in_number == -1){
                str = SystemDisplay::getCorrespondingSystemParameterName();
                if(str.compare("Pressure") == 0){
                    corresponding_system_parameter_name_in_number = 0;
                }
                else if(str.compare("Liquid_Temperature") == 0){
                    corresponding_system_parameter_name_in_number = 1;
                }
                else if(str.compare("Vapor_Temperature") == 0){
                    corresponding_system_parameter_name_in_number = 2;
                }
                else if(str.compare("Saturation_Temperature") == 0){
                    corresponding_system_parameter_name_in_number = 3;
                }
                else if(str.compare("Boron_Density") == 0){
                    corresponding_system_parameter_name_in_number = 4;
                }
                else if(str.compare("Equilibrium_Quality") == 0){
                    corresponding_system_parameter_name_in_number = 5;
                }
                else if(str.compare("Flow_regime_Number") == 0){
                    corresponding_system_parameter_name_in_number = 6;
                }
                else {
                    cout << "SystemDisplay::getCorrespondingSystemParameterValue (SystemDisplay::getCorrespondingSystemName()): "
                        << (SystemDisplay::getCorrespondingSystemName()) << " can not be recognized" << endl;
                    ofstream fout("err messages.txt", ios::app);
                    fout << r5pars_.timesno << "  SystemDisplay::getCorrespondingSystemParameterValue (SystemDisplay::getCorrespondingSystemName()): "
                        << (SystemDisplay::getCorrespondingSystemName()) << " can not be recognized\n";
                    fout.close();
                    r5pars_.sim_err = true;
                    g_sim.Branch_center.setTempESMessage("SystemDisplay::setCorrespondingSystemParameterUnit err.");
                    g_sim.Branch_center.setTempTypeOfConsequence(Consequence::AdsSimErr);
                }
            }
            if(!r5pars_.sim_err){
                hv_ptr = (HydraulicVolume *) ptr_to_corresponding_TH_element;
                switch (corresponding_system_parameter_name_in_number){
                case 0:
                    str = hv_ptr->getUnitOfPressure();
                    if(type_of_display == DisplayType::ParameterTrend) str += "_per_min";                       
                    DisplayElement::setCorrespondingSystemParameterUnit(str);
                    break;
                case 1:
                    str = hv_ptr->getUnitOfLiquidTemperature();
                    if(type_of_display == DisplayType::ParameterTrend) str += "_per_min";                       
                    DisplayElement::setCorrespondingSystemParameterUnit(str);
                    break;
                case 2:
                    str = hv_ptr->getUnitOfVaporTemperature();
                    if(type_of_display == DisplayType::ParameterTrend) str += "_per_min";                       
                    DisplayElement::setCorrespondingSystemParameterUnit(str);
                    break;
                case 3:
                    str = hv_ptr->getUnitOfSaturationTemperature();
                    if(type_of_display == DisplayType::ParameterTrend) str += "_per_min";                       
                    DisplayElement::setCorrespondingSystemParameterUnit(str);
                    break;
                case 4:
                    str = hv_ptr->getUnitOfBoronDensity();
                    if(type_of_display == DisplayType::ParameterTrend) str += "_per_min";                       
                    DisplayElement::setCorrespondingSystemParameterUnit(str);
                    break;
                    //case 5 and case 6 do not have units
                default:
                    break;
                }
            }
            break;
        case Relap5ComponentType::HydraulicJunction:
            if (corresponding_system_parameter_name_in_number == -1){
                str = SystemDisplay::getCorrespondingSystemParameterName();
                if(str.compare("Mass_Flow_Rate") == 0){
                    corresponding_system_parameter_name_in_number = 0;
                }
                else if(str.compare("Liquid_Velocity") == 0){
                    corresponding_system_parameter_name_in_number = 1;
                }
                else if(str.compare("Vapor_Velocity") == 0){
                    corresponding_system_parameter_name_in_number = 2;
                }
                else {
                    cout << "SystemDisplay::getCorrespondingSystemParameterValue (SystemDisplay::getCorrespondingSystemName()): "
                        << (SystemDisplay::getCorrespondingSystemName()) << " can not be recognized" << endl;
                    ofstream fout("err messages.txt", ios::app);
                    fout << r5pars_.timesno << "  SystemDisplay::getCorrespondingSystemParameterValue (SystemDisplay::getCorrespondingSystemName()): "
                        << (SystemDisplay::getCorrespondingSystemName()) << " can not be recognized\n";
                    fout.close();
                    r5pars_.sim_err = true;
                    g_sim.Branch_center.setTempESMessage("SystemDisplay::setCorrespondingSystemParameterUnit err.");
                    g_sim.Branch_center.setTempTypeOfConsequence(Consequence::AdsSimErr);
                }
            }
            if(!r5pars_.sim_err){
                hj_ptr = (HydraulicJunction *) ptr_to_corresponding_TH_element;
                switch (corresponding_system_parameter_name_in_number){
                case 0:
                    str = hj_ptr->getUnitOfMassFlowRate();
                    if(type_of_display == DisplayType::ParameterTrend) str += "_per_min";                       
                    DisplayElement::setCorrespondingSystemParameterUnit(str);
                    break;
                case 1:
                    str = hj_ptr->getUnitOfLiquidVelocity();
                    if(type_of_display == DisplayType::ParameterTrend) str += "_per_min";                       
                    DisplayElement::setCorrespondingSystemParameterUnit(str);
                    break;
                case 2:
                    str = hj_ptr->getUnitOfVaporVelocity();
                    if(type_of_display == DisplayType::ParameterTrend) str += "_per_min";                       
                    DisplayElement::setCorrespondingSystemParameterUnit(str);
                    break;
                default:
                    break;
                }   
            }
            break;
        case Relap5ComponentType::VariableTrip: 
            DisplayElement::setCorrespondingSystemParameterUnit("sec");
            break;
        case Relap5ComponentType::LogicalTrip:
            DisplayElement::setCorrespondingSystemParameterUnit("sec");
            break;
        case Relap5ComponentType::ControlVariable:
            cv_ptr = (ControlVariable *) ptr_to_corresponding_TH_element;
            str = cv_ptr->getUnitOfControlValue();
            if(type_of_display == DisplayType::ParameterTrend) str += "_per_min";                       
            DisplayElement::setCorrespondingSystemParameterUnit(str);
            break;
        case Relap5ComponentType::HeatStructure:
            hs_ptr = (HeatStructure *) ptr_to_corresponding_TH_element;
            str = hs_ptr->getUnitOfMeshPointTemperature();
            if(type_of_display == DisplayType::ParameterTrend) str += "_per_min";                       
            DisplayElement::setCorrespondingSystemParameterUnit(str);
            break;
        default:
//      case Relap5ComponentType::InteractiveVariable does not have unit
            break;
        }
    }
    else{//could not find the corresponding component
        cout << "DisplayElement::getCorrespondingSystemParameterValue can not find the corresponding component: "
            << DisplayElement::getCorrespondingSystemName() << endl;
        ofstream fout("err messages.txt", ios::app);
        fout << r5pars_.timesno << "  DisplayElement::getCorrespondingSystemParameterValue can not find the corresponding component: "
            << DisplayElement::getCorrespondingSystemName() << '\n';
        fout.close();
        r5pars_.sim_err = true;
        g_sim.Branch_center.setTempESMessage("SystemDisplay::setCorrespondingSystemParameterUnit err.");
        g_sim.Branch_center.setTempTypeOfConsequence(Consequence::AdsSimErr);
    }
}
void SystemDisplay::setCorrespondingSystemParameterNameInNumber(int n){corresponding_system_parameter_name_in_number = n;}
int SystemDisplay::getCorrespondingSystemParameterNameInNumber(){return corresponding_system_parameter_name_in_number;}
void SystemDisplay::setPtrToCorrespondingTHElement(THElement *_th){ptr_to_corresponding_TH_element = _th;}
THElement* SystemDisplay::getPtrToCorrespondingTHElement(){return ptr_to_corresponding_TH_element;}
 
void SystemDisplay::copySystemDisplay(SystemDisplay &_sdp){
    this->copyDisplayElement(_sdp);
    type_of_display = _sdp.getTypeOfDisplay();
    corresponding_system_parameter_name_in_number = _sdp.getCorrespondingSystemParameterNameInNumber();
    functional_state = _sdp.getFunctionalState();
}
void SystemDisplay::saveSystemDisplayToFile(ofstream &fout){
    this->saveDisplayElementToFile(fout);
    fout << static_cast<int>(type_of_display) << ' ' << corresponding_system_parameter_name_in_number << ' ' << functional_state << '\n';
}
void SystemDisplay::retrieveSystemDisplayFromFile(ifstream &fin){
    this->retrieveDisplayElementFromFile(fin);
    { int td; fin >> td >> corresponding_system_parameter_name_in_number >> functional_state; type_of_display = static_cast<DisplayType>(td); }
}