HydraulicVolume.cpp

Go to the documentation of this file.

// HydraulicVolume.cpp: implementation of the HydraulicVolume class.
//
 
#include "HydraulicVolume.h"
#include <boost/serialization/export.hpp>
#include <boost/archive/binary_oarchive.hpp>
#include <boost/archive/binary_iarchive.hpp>
BOOST_CLASS_EXPORT_IMPLEMENT(HydraulicVolume)
using namespace std;
 
// Construction/Destruction
 
HydraulicVolume::HydraulicVolume()
{
    liquid_temperature = -999.0;
    vapor_temperature = -999.0;
    pressure = -999.0;
    saturation_temperature = -999.0;
    boron_density = -999.0; 
    equilibrium_quality = -999.0; 
    flow_regime_number = -999.0; 
    unit_of_liquid_temperature = "NOT_SPECIFIED";
    unit_of_vapor_temperature = "NOT_SPECIFIED";
    unit_of_pressure = "NOT_SPECIFIED";
    unit_of_saturation_temperature = "NOT_SPECIFIED";
    unit_of_boron_density = "NOT_SPECIFIED"; 
    THElement::setTypeOfCard(Relap5ComponentType::HydraulicVolume);
}
 
HydraulicVolume::~HydraulicVolume()
{
 
}
 
void HydraulicVolume::setLiquidTemperature (double x){liquid_temperature = x;}
void HydraulicVolume::setVaporTemperature(double x){vapor_temperature = x;}
void HydraulicVolume::setPressure(double x){pressure = x;}
void HydraulicVolume::setSaturationTemperature(double x){saturation_temperature = x;}
void HydraulicVolume::setBoronDensity(double x){boron_density = x;} 
void HydraulicVolume::setEquilibriumQuality(double x){equilibrium_quality = x;}
void HydraulicVolume::setFlowRegimeNumber(double x){flow_regime_number = x;}
double HydraulicVolume::getLiquidTemperature(){return liquid_temperature;}
double HydraulicVolume::getVaporTemperature(){return vapor_temperature;}
double HydraulicVolume::getPressure(){return pressure;}
double HydraulicVolume::getSaturationTemperature(){return saturation_temperature;}
double HydraulicVolume::getBoronDensity(){return boron_density;}
double HydraulicVolume::getEquilibriumQuality(){return equilibrium_quality;}
double HydraulicVolume::getFlowRegimeNumber(){return flow_regime_number;}
void HydraulicVolume::setUnits(bool si_unit){
    this->setUnitOfBoronDensity(si_unit);
    this->setUnitOfLiquidTemperature(si_unit);
    this->setUnitOfPressure(si_unit);
    this->setUnitOfSaturationTemperature(si_unit);
    this->setUnitOfVaporTemperature(si_unit);
}
void HydraulicVolume::setUnitOfLiquidTemperature(bool si_unit){
    if(si_unit){
        unit_of_liquid_temperature = "K";}
    else{
        unit_of_liquid_temperature = "F";}
}
void HydraulicVolume::setUnitOfVaporTemperature(bool si_unit){
    if(si_unit){
        unit_of_vapor_temperature = "K";}
    else{
        unit_of_vapor_temperature = "F";}
}
void HydraulicVolume::setUnitOfPressure(bool si_unit){
    if(si_unit){
        unit_of_pressure = "Pa";}
    else{
        unit_of_pressure = "lb/in2";}
};
void HydraulicVolume::setUnitOfSaturationTemperature(bool si_unit){
    if(si_unit){
        unit_of_saturation_temperature = "K";}
    else{
        unit_of_saturation_temperature = "F";}
}
 
void HydraulicVolume::setUnitOfBoronDensity(bool si_unit){
        if(si_unit){
        unit_of_boron_density = "kg/m3";}
    else{
        unit_of_boron_density = "lb/ft3";}
}
string HydraulicVolume::getUnitOfLiquidTemperature(){return unit_of_liquid_temperature;}
string HydraulicVolume::getUnitOfVaporTemperature(){return unit_of_vapor_temperature;}
string HydraulicVolume::getUnitOfPressure(){return unit_of_pressure;}
string HydraulicVolume::getUnitOfSaturationTemperature(){return unit_of_saturation_temperature;}
string HydraulicVolume::getUnitOfBoronDensity(){return unit_of_boron_density;}
void HydraulicVolume::readInfoFromFile(ifstream &fin){
    THElement::readInfoFromFile(fin);
}
void HydraulicVolume::saveHydraulicVolumeToFile(ofstream &fout){
    double ltemp, vtemp, xpress, tsat, bor, qual;
    //
    ltemp = liquid_temperature;
    if (ltemp > -1e-30 && ltemp < 1e-30) {
        ltemp = 0.0;
    }
    vtemp = vapor_temperature;
    if (vtemp > -1e-30 && vtemp < 1e-30) {
        vtemp = 0.0;
    }
    xpress = pressure;
    if (xpress > -1e-30 && xpress < 1e-30) {
        xpress = 0.0;
    }
    tsat = saturation_temperature;
    if (tsat > -1e-30 && tsat < 1e-30) {
        tsat = 0.0;
    }
    bor = boron_density;
    if (bor > -1e-30 && bor < 1e-30) {
        bor = 0.0;
    }
    qual = equilibrium_quality;
    if (qual > -1e-30 && qual < 1e-30) {
        qual = 0.0;
    }
    //
    fout << ltemp << ' '
    << vtemp  << ' '
    << xpress << ' '
    << tsat   << ' '
    << bor    << ' '
    << qual   << ' '
    
    //fout << liquid_temperature << ' '
    //<< vapor_temperature << ' '
    //<< pressure << ' '
    //<< saturation_temperature << ' '
    //<< boron_density << ' '
    //<< equilibrium_quality << ' '
    //
    << flow_regime_number << ' '
    << unit_of_liquid_temperature << ' '
    << unit_of_vapor_temperature << ' '
    << unit_of_pressure << ' '
    << unit_of_saturation_temperature << ' '
    << unit_of_boron_density << '\n';
    this->saveTHElementToFile(fout);
}
void HydraulicVolume::retrieveHydraulicVolumeFromFile(ifstream &fin){
    fin >> liquid_temperature
    >> vapor_temperature
    >> pressure
    >> saturation_temperature
    >> boron_density
    >> equilibrium_quality
    >> flow_regime_number
    >> unit_of_liquid_temperature
    >> unit_of_vapor_temperature
    >> unit_of_pressure
    >> unit_of_saturation_temperature
    >> unit_of_boron_density;
    this->retrieveTHElementFromFile(fin);
}
void HydraulicVolume::copyHydraulicVolume(HydraulicVolume &_hv){
    this->copyTHElement(_hv);
    liquid_temperature = _hv.getLiquidTemperature();//(tempf)
    vapor_temperature = _hv.getVaporTemperature();//(tempg)
    pressure = _hv.getPressure();//(p)
    saturation_temperature = _hv.getSaturationTemperature();//(satt) based on the steam partial pressure
    boron_density = _hv.getBoronDensity();//(boron) mass of boron per celll volume 
    equilibrium_quality = _hv.getEquilibriumQuality(); //(quale)
    flow_regime_number = _hv.getFlowRegimeNumber(); //(floreg)
    unit_of_liquid_temperature = _hv.getUnitOfLiquidTemperature();//(tempf)
    unit_of_vapor_temperature = _hv.getUnitOfVaporTemperature();//(tempg)
    unit_of_pressure = _hv.getUnitOfPressure();//(p)
    unit_of_saturation_temperature = _hv.getUnitOfSaturationTemperature();//(satt) based on the steam partial pressure
    unit_of_boron_density = _hv.getUnitOfBoronDensity();//(boron) mass of boron per celll volume 
}