diff --git a/source/trunk/Model/Equivalence/EQM_POL_PWR_UO2.cxx b/source/trunk/Model/Equivalence/EQM_POL_PWR_UO2.cxx new file mode 100644 index 0000000000000000000000000000000000000000..ebf2dd1f7ca926cf1651ab71648b90850fe4f130 --- /dev/null +++ b/source/trunk/Model/Equivalence/EQM_POL_PWR_UO2.cxx @@ -0,0 +1,34 @@ +#include "EquivalenceModel.hxx" +#include "EQM_POL_PWR_UO2.hxx" +#include "CLASSLogger.hxx" +/*Whatever include you need */ +// ________________________________________________________________________ +// EQM_POL_PWR_UO2 +// +// Brief description +// ________________________________________________________________________ +//Constructor(s) +EQM_POL_PWR_UO2::EQM_POL_PWR_UO2 ( /*parameters*/ ) +{ +//.....Do whatever you want with your parameters +/* +Fill the two isotopic vectors fFissileList and fFertileList +see explanation in the manual */ +// Fertile +ZAI U8(92 ,238 ,0) ; +fFertileList = U8*1; +// Fissile +ZAI U5(92 ,235 ,0) ; +// ... +fFissileList = U5*1; +} +// _______________________________________________________________________ +double EQM_POL_PWR_UO2::GetFissileMolarFraction ( IsotopicVector Fissil , IsotopicVector Fertil , double BurnUp ) +{ +/*Code your Equivalence Model : This function has to return the molar fraction of fissile in +the fuel needed to reach the BurnUp(GWd/tHM) according to +the composition of the Fissil and Fertil vectors +*/ +return 0.0125575 + 0.00053602*BurnUp + 1.80023e-06*BurnUp*BurnUp ; + +} \ No newline at end of file diff --git a/source/trunk/Model/Equivalence/EQM_POL_PWR_UO2.hxx b/source/trunk/Model/Equivalence/EQM_POL_PWR_UO2.hxx new file mode 100644 index 0000000000000000000000000000000000000000..90add1cbd074ebf13f990da83e36fa012b8b9d3e --- /dev/null +++ b/source/trunk/Model/Equivalence/EQM_POL_PWR_UO2.hxx @@ -0,0 +1,26 @@ +#ifndef _EQM_POL_PWR_UO2_HXX +#define _EQM_POL_PWR_UO2_HXX +#include "EquivalenceModel.hxx" +using namespace std; +//−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−// +/*! +Define a EQM_POL_PWR_UO2 +Explain briefly what is it . @author YourName +@version 3.0 +*/ +//________________________________________________________________________ +class EQM_POL_PWR_UO2 : public EquivalenceModel +{ + public: +/*Constructor*/ +EQM_POL_PWR_UO2(/*parameters*/ ); //!< Explain what is the parameters (if any) + +/**This function IS the equivalence model**/ +double GetFissileMolarFraction(IsotopicVector Fissil, IsotopicVector Fertil,double BurnUp) ; // !<Return the molar fraction of fissile element + + private: +/*Your private variables */ + +}; + +#endif \ No newline at end of file diff --git a/source/trunk/src/Makefile b/source/trunk/src/Makefile index 8c2be94fe6104ec8f446c503353f572e0861db4a..68387ed5a8a6525b7b6aeb0c7ace742881c5b0c3 100755 --- a/source/trunk/src/Makefile +++ b/source/trunk/src/Makefile @@ -40,7 +40,7 @@ OBJS = CLASSLogger.o \ CLASSFuelPlan.o\ Scenario.o -OBJMODEL = $(EQM)/EQM_MLP_PWR_MOX.o $(EQM)/EQM_QUAD_PWR_MOX.o $(EQM)/EQM_LIN_PWR_MOX.o \ +OBJMODEL = $(EQM)/EQM_MLP_PWR_MOX.o $(EQM)/EQM_QUAD_PWR_MOX.o $(EQM)/EQM_LIN_PWR_MOX.o $(EQM)/EQM_POL_PWR_UO2.o \ $(XSM)/XSM_MLP.o $(XSM)/XSM_CLOSEST.o \ $(IM)/IM_RK4.o $(IM)/IM_Matrix.o