/************************************************************/ // DESCRIPTION // // // Simple Scenario with one PWR UOX and one PWR MOX // // // _______ ____ _______ ___________ // | | | | | | | | // |Reactor| =>|Pool|=>|Storage|=>|Fabrication| // | UOx | |UOX | | UOX | | Plant | // |_______| |____| |_______| |___________| // || // \/ // _______ ____ _______ // | | | | | | // |Storage|<= |Pool|<==|Reactor| // | MOX | |MOX | | MOX | // |_______| |____| |_______| // // //@author FrA /***********************************************************/ #include "CLASSHeaders.hxx" #include #include #include #include #include "XS/XSM_MLP.hxx" //Load the include for Neural network cross section predictor #include "Irradiation/IM_RK4.hxx" //Load the include for Runge Kutta 4 resolution #include "Equivalence/EQM_PWR_MLP_MOX.hxx" //Load the include for Neural Network Equivalence Model (PWRMOX) using namespace std; int main(int argc, char** argv) { //seconds in one year cSecond year = 3600*24.*365.25; /******LOG MANAGEMENT**********************************/ //Definition of the Log file : CLASS messages output int Std_output_level = 0; // Only error are shown in terminal int File_output_level = 2; // Error + Warning + Info are shown in the file CLASS_OUTPUT.log CLASSLogger *Logger = new CLASSLogger("CLASS_OUTPUT.log",Std_output_level,File_output_level); /******SCENARIO****************************************/ // The scenario start at year 2015 Scenario *gCLASS=new Scenario(2015*year,Logger); gCLASS->SetStockManagement(true); //If false all the IsotopicVector in stocks are mixed together. gCLASS->SetTimeStep(year/4.); //the scenario calculation is updated every 3 months cSecond EndOfScenarioTime=2060*year; //Scenario ends in year 2060 gCLASS->SetOutputFileName("Scenario6.root");//Set the name of the output file /******DATA BASES**************************************/ //Geting CLASS to path string CLASS_PATH = getenv("CLASS_PATH"); if (CLASS_PATH=="") { cout<<" Please setenv CLASS_PATH to your CLASS installation folder in your .bashs or .tcshrc"<GetLog(), PATH_TO_DATA + "DECAY/ALL/Decay.idx"); //you may have to open this file and do the proper changes according your path gCLASS->SetDecayDataBase(DecayDB); //This decay data base will be used for all the decay calculations in this Scenario /*===Reactor data base===*/ // Reprocessed fuel PWR MOX XSM_MLP* XSMOX = new XSM_MLP(gCLASS->GetLog(), PATH_TO_DATA + "PWR/MOX/XSModel/30Wg_FullMOX"); //Defining the XS Predictor IM_RK4 *IMRK4 = new IM_RK4(gCLASS->GetLog()); //Bateman's equation solver method (RungeKutta4) //NO PATH FOR THE RK SOLVER??? EQM_PWR_MLP_MOX* EQMMLPPWRMOX = new EQM_PWR_MLP_MOX(gCLASS->GetLog(),PATH_TO_DATA + "PWR/MOX/EQModel/MLP/EQM_MLP_PWR_MOX_3batch.xml");//Defining the EquivalenceModel PhysicsModels* PHYMOD = new PhysicsModels(XSMOX, EQMMLPPWRMOX, IMRK4);//The PhysicsModels containing the 3 object previously defined //Fixed fuel : PWR UOX EvolutionData *CYCLADE =new EvolutionData(gCLASS->GetLog(), PATH_TO_DATA + "PWR/UOX/FixedFuel/CYCLADES.dat"); EvolutionData *GARANCE = new EvolutionData(gCLASS->GetLog(), PATH_TO_DATA + "PWR/UOX/FixedFuel/GARANCE.dat"); EvolutionData *STD900 = new EvolutionData(gCLASS->GetLog(), PATH_TO_DATA + "PWR/UOX/FixedFuel/STD900.dat"); /******FACILITIES*************************************/ /*=== Stock===*/ //Storage for infinite U and Pu Storage *Stock = new Storage(gCLASS->GetLog()); //Definition of the stock Stock->SetName("Stock"); //Its name Stock->AddToStock(ZAI(92,238,0) * 3e31); Stock->AddToStock(ZAI(92,235,0) * 0.04375* 3e31); /* Stock->AddToStock(ZAI(94,238,0) * 10e28); Stock->AddToStock(ZAI(94,239,0) * 10e28); Stock->AddToStock(ZAI(94,240,0) * 10e28); Stock->AddToStock(ZAI(94,241,0) * 10e28); Stock->AddToStock(ZAI(94,242,0) * 10e28);*/ gCLASS->Add(Stock); //Adding the stock to the Scenario //Storage for MOX Storage *StockMOX = new Storage(gCLASS->GetLog());//Definition of the stock StockMOX->SetName("StockMOX"); //Its name gCLASS->Add(StockMOX); //Adding the stock to the Scenario //Storage for UOX Storage *StockUOX = new Storage(gCLASS->GetLog());//Definition of the stock StockUOX->SetName("StockUOX"); //Its name gCLASS->Add(StockUOX); //Adding the stock to the Scenario /*===Pool===*/ //Pool for UOX Pool *Cooling_UOX = new Pool(gCLASS->GetLog(),StockUOX, 5*year); //After 5 years of cooling, the pool sends its content to "StockUOX" Cooling_UOX->SetName("Pool_UOX"); gCLASS->Add(Cooling_UOX); //Pool for MOX Pool *Cooling_MOX = new Pool(gCLASS->GetLog(),StockMOX, 5*year); //After 5 years of cooling, the pool sends its content to "StockMOX" Cooling_MOX->SetName("Pool_MOX"); gCLASS->Add(Cooling_MOX); /*===A FabricationPlant===*/ FabricationPlant *FP_MOX = new FabricationPlant(gCLASS->GetLog(), 3*year); //Declare a FabricationPlant. After the build of the fuel, it decays during 3years before to be loaded in Reactor FP_MOX->SetFiFo(false); //The latest isotopicVector to enter in "Stock" will be used to build the fuel (Opposite of First In First Out) FP_MOX->SetName("Fab_MOX"); FP_MOX->AddFissileStorage(StockUOX); //Tell the FP to look in Stock for fissionable material FP_MOX->AddFertileStorage(Stock);//Tell the FP to look in Stock for fertile material //If fertile stock is not defined (like here), CLASS get fertile from nature (OUTCOMING vector) //FP_MOX->SetReUsableStorage(wastestock);//By default the fabricationplant get the list of nuclei defined in the EquivalenceModel (here EQM_MLP_MOX) from stock and send the others nuclei in WASTE. If user want these nuclei to go in another stock he can use the SetReUsableStorage function gCLASS->AddFabricationPlant(FP_MOX); /*=== Reactor===*/ double Power_CP0 = 2.66e10; //Thermal power (in W) //Combustibles type : CYCLADE double BurnUp_Cyclade = 47; // GWd/tHM double HMMass_Cyclade = 72.3; //heavy metal mass (in tons) //REACTOR PWR UOX cSecond StartingTime = 2016*year; cSecond LifeTime = (EndOfScenarioTime - StartingTime); Reactor* PWR_UOX = new Reactor(gCLASS->GetLog(),// Log CYCLADE, // The DataBase used Cooling_UOX, // Connected Backend facility : The reactor discharge its fuel into the Pool "Cooling_UOX" StartingTime,// Starting time LifeTime, // time of reactor life time Power_CP0, // Power HMMass_Cyclade,// HM mass BurnUp_Cyclade,// BurnUp 0.8); // Load Factor PWR_UOX->SetName("PWR_UOX");// name of the reactor (as it will show up in the CLASSGui) gCLASS->AddReactor(PWR_UOX);//Add this reactor to the scenario //POWER reactor MOX double Power_CPY = 2.785e9; //Thermal power (in W) //Combustibles type : STANDARD 900 double BurnUpMOX = 35; //GWd/tHM //double BurnUp_STD900 = 33; // GWd/tHM double HMMass_STD900 = 72.3; //heavy metal mass (in tons) //the PWR MOX StartingTime = 2022*year; LifeTime = EndOfScenarioTime - StartingTime; Reactor* PWR_MOX = new Reactor(gCLASS->GetLog(),// Log PHYMOD, // The models used to build the fuel & to calculate its evolution FP_MOX, // The FabricationPlant Cooling_MOX, // Connected Backend StartingTime, // Starting time LifeTime, // time of reactor l Power_CPY, // Power HMMass_STD900, // HM mass BurnUpMOX, // BurnUp 0.8); // Load Factor PWR_MOX->SetName("PWR_MOX");// name of the reactor (as it will appear in the CLASSGui) gCLASS->AddReactor(PWR_MOX);//Add this reactor to the scenario gCLASS->Evolution((double)EndOfScenarioTime);//Perform the calculation from year 2015(defined in Scenario declaration) to year 2060 delete gCLASS; } //========================================================================================== // Compilation //========================================================================================== /* \rm CLASS* ; g++ -o CLASS_Exec Scenario6.cxx -I $CLASS_include -L $CLASS_lib -lCLASSpkg `root-config --cflags` `root-config --libs` -fopenmp -lgomp -Wunused-result */