\param map < string , vector <IsotopicVector> > StreamArray, the string is the stream code (fissile fertile ,...) the IsotopicVector the fraction of each IV to take in the (fissile, fertile,..) stock .
IsotopicVectorGetFertileList(){returnfFertileList;}//!<return the fertile list
IsotopicVectorGetFissileList(){returnfFissileList;}//!<return the fissile list
doubleGetBuildFuelFirstGuess(){returnfFirstGuessFissilContent;}//!<Get the initialization value for BuildFuel algorithm
virtualdoubleGetFissileMolarFraction(IsotopicVectorFissil,IsotopicVectorFertil,doubleBurnUp){return0;}//!< Return the molar fraction of fissile element in the fuel according to the burnup, and a given fuel composition (this is the heart of the equivalence model)
map<string,double>GetBuildFuelFirstGuess(){returnfFirstGuessContent;}//!<Get the initialization value for BuildFuel algorithm
virtualmap<string,double>GetMolarFraction(map<string,IsotopicVector>IVStream,doubleBurnUp)=0;//!< Return the molar fractions of each element in the fuel according to the burnup, and a given fuel composition (this is the heart of the equivalence model)
doubleGetRelativMassPrecision()const{returnfRelativMassPrecision;}//!< Mass precision
intGetMaxInterration()const{returnfMaxInterration;}//!< Max iterration in build fueld algorythm
doubleGetActualFissileContent()const{returnfActualFissileContent;}//!< Get the fissile content at the actual algorithm step (usefull for EQM_MLP_Kinf)
voidSetFertileList(IsotopicVectorIV){fFertileList=IV;}//!<set the fertile list
voidSetFissileList(IsotopicVectorIV){fFissileList=IV;}//!<set the fissile list
voidSetBuildFuelFirstGuess(doubleFirstGuess){fFirstGuessFissilContent=FirstGuess;}//!<set the initialization value for BuildFuel algorithm
voidSetBuildFuelFirstGuess(map<string,double>FirstGuess){fFirstGuessContent=FirstGuess;}//!<set the initialization value for BuildFuel algorithm (one FistGuess for each component of the fresh fuel
voidSetRelativMassPrecision(doubleval){fRelativMassPrecision=val;}//!< Mass precision
voidSetMaxInterration(intval){fMaxInterration=val;}//!< Max iterration in build fueld algorythm
voidSetMaxInterration(intval){fMaxInterration=val;}//!< Max iteration in build fuel algorithm
//@}
//@}
/*!
\name Reading NFO related Method
...
...
@@ -140,61 +139,42 @@ class EquivalenceModel : public CLASSObject
*/
voidReadMaximalContent(conststring&line);
//}
voidPrintInfo();//Print the information red in the INFO stream
//{
/// ReadFissil : read the zai name in the TMWA MLP model starts with "k_fissil" keyword
/*!
\param line : line suppossed to contain the fissil list
*/
voidReadFissil(conststring&line);
//}
//{
/// ReadFertil : read the zai name in the TMWA MLP model starts with "k_fertil" keyword
/*!
\param line : line suppossed to contain the fertil list & their default isotopic fraction
*/
voidReadFertil(conststring&line);
//}
voidPrintInfo();//Print the information red in the INFO stream
voidReadList(conststring&line);
//@}
voidReadFirstGuessContent(conststring&line);
boolisIVInDomain(IsotopicVectorIV);
protected:
IsotopicVectorfFertileList;//!< contain the list of zai, needed as fertile, taken in a stock before fabrication
//!< if no stock are provided, will take the isotopic vector in the Park income
IsotopicVectorfFissileList;//!< contain the list of zai, needed as fissile, taken in a stock before fabrication
//!< if no stock are provided the fuel will not be made
doublefFirstGuessFissilContent;//!< fissile content for BuildFuel initialization (in weight proportion)
doublefActualFissileContent;//!< fissile content at the actual dichotomy step (usefull for EQM_MLP_Kinf)
protected:
map<string,IsotopicVector>fStreamList;//!< contains all lists of zai needed to build a fuel (example : 2 -> fissileList+fertileList)
//!< each list is identified by a keyword (example : -> "Fissile" & "Fertile")
map<string,double>fFirstGuessContent;//!< fissile content for BuildFuel initialization (in weight proportion)
map<string,double>fActualMolarContentInFuel;//!< Molar Content in fuel of each list at this step of the calculation
map<string,double>fActualMassContentInFuel;//!< Mass Content in fuel of each list at this step of the calculation
doublefSpecificPower;//!< The specific power in W/gHM (HM: heavy Metal)
doublefMaximalContent;//!< The approx. maximum fissile content reachable by the model
protected:
/*!
\name Algorithm parameters
*/
//@{
doublefRelativMassPrecision;//!< Mass precision
intfMaxInterration;//!< Max iterration in build fueld algorythm
//@}
doubleLambdaCalculation(stringMaterialDenomination,doubleLambdaPreviousStep,doubleMaterialMassNeeded,doubleDeltaMass,vector<IsotopicVector>Stocks);//!< Calculate the proportion of each stocks in StockArray to take in oder to get a mass of MassNeeded (can be Fer(fertile) or Fis(Fissile))
voidSetLambda(vector<double>&lambda,doubleLambda_tot);//!< Set individual lambda according to the LAMBDA_TOT (lambda of all stocks)
intfMaxInterration;//!< Max iterration in build fueld algorythm
/*!
\name Algorithm related functions for default BuildFuel function
*/
//@{
voidSetLambda(vector<double>&lambda,intFirstStockID,intLastStockID,doubleLAMBDA_TOT);//!< Set individual lambda according to the LAMBDA_TOT (lambda of all stocks)
voidSetLambdaToErrorCode(vector<double>&lambda);//!< Set all vector elements to -1
doubleLAMBDA_TOT_FOR(doubleMassNeeded,vector<IsotopicVector>StockArray,stringFisOrFer);//!< Calculate the proportion of each stocks in StockArray to take in oder to get a mass of MassNeeded (can be Fer(fertile) or Fis(Fissile))
boolBuild_Fuel_According_Lambda(vector<double>&lambda,vector<IsotopicVector>FissilArray,vector<IsotopicVector>FertilArray,doubleHMMass,IsotopicVector&Fissile,IsotopicVector&Fertile);//!< Compute the fuel isotopy according to the proportion of each IV taken
//@}
...
...
@@ -203,21 +183,23 @@ class EquivalenceModel : public CLASSObject
#endif
boolfreaded;
map<ZAI,pair<double,double>>fZAILimits;//!< Fresh fuel range : map<ZAI<min edge ,max edge >>
map<ZAI,pair<double,double>>fZAILimits;//!< Fresh fuel range : map<ZAI<min edge ,max edge >>
stringfInformationFile;//!< file containing Reactor Type, Fuel type, HM mass, Power, time vector, and TMVA input variables names (looks the manual for format details)
stringfDBFType;//!< Fuel Type (e.g MOX, UOX, ThU, ThPu ...)
stringfDBRType;//!< Reactor Type (e.g PWR, FBR-Na, ADS..)
doublefSpecificPower;//!< The specific power in W/gHM (HM: heavy Metal)
doublefMaximalContent;//!< The approx. maximum fissile content reachable by the model
stringfInformationFile;//!< file containing Reactor Type, Fuel type, HM mass, Power, time vector, and TMVA input variables names (looks the manual for format details)
stringfDBFType;//!< Fuel Type (e.g MOX, UOX, ThU, ThPu ...)
stringfDBRType;//!< Reactor Type (e.g PWR, FBR-Na, ADS..)
/*!
\name Others
*/
//@{
doublefTotalFissileMassInStocks;//!< Total mass in fissile stocks
doublefTotalFertileMassInStocks;//!< Total mass in fertile stocks
map<string,double>fTotalMassInStocks;//!< Total mass in each vector of stock
map<string,double>fLambdaMax;//!< Total lambda of available stocks
