ReactorMesh.hxx

Go to the documentation of this file.

#ifndef _ReactorMesh_
#define _ReactorMesh_

#include "MureHeaders.hxx"
using namespace std;



class ReactorMesh
{
 public:
    ReactorMesh(
    int NZones=1,           
    int NLevels=1,          
    string AssemblyName="A", 
    bool Hex=false          
    ); 
    virtual ReactorMesh* Clone(){return new ReactorMesh(*this);}
    virtual ~ReactorMesh(); 



    void SetExternalShape(
    Shape_ptr TheShape      
    );
    void SetSpecialExternalShape(
    Shape_ptr TheShape          
    ) {fSAssemblyZone = TheShape;}      
    void Translate(
    double dx,          
    double dy,          
    double dz           
    );
    void SetDimensions(
    double FuelRadius,      
    double CladThickness,   
    double Pitch,           
    int Zone=-1             
    );                   
    void SetSpecialDimensions(
    double FuelPinRadius,   
    double CladThickness,   
    double Pitch,           
    double R1,              
    int Zone=-1             
    );                   
    void SetGTDimensions(
    double InnerCladRadius, 
    double CladThickness,   
    int Zone=-1             
    );                   
    void SetCRDimensions(
    double InnerCRRadius,       
    double CRCladThickness,     
    double CoolantGapThickness, 
    double CladThickness,       
    int Zone=-1             
    );                   
    void SetCenteredPosition(
    double dx=0.,                       
    double dy=0.                        
    ) {fXDeviation+=dx; fYDeviation+=dy;fXCenterDeviationOfMesh=dx;fYCenterDeviationOfMesh=dy;} 
    void SetGasSpaceData(
    double GasThickness=0,  
    double GasPressure=0,   
    double HeFraMol=0,      
    double XeFraMol=0,      
    int Zone=-1             
    );



    void AddMaterials(
    Material *Fuel,         
    Material *Cladding,     
    Material *Coolant,      
    int Zone=-1,            
    int Level=-1            
    ); 
    void AddSpecialMaterials(
    Material *Fuel,         
    Material *Cladding,     
    Material *Coolant,      
    Material *Spec1,        
    int Zone=-1             
    ); 
    void AddGTMaterials(
    Material *GTInside,     
    Material *GTCladding,   
    int Zone=-1             
    ); 
    void AddCRMaterials(
    Material *CRMat,        
    Material *CRCladding,   
    Material *CoolantGap,   
    Material *ExtCladding,  
    int Zone=-1             
    ); 



    void AddCircleZoneRadius(
    double Radius       
    );

    void AddPinPos(
    int X,      
    int Y,      
    int Zone    
    );


    void AddGuideTube(
    int X,      
    int Y,      
    int Z=-1    
    );

    void AddControlRod(
    int X,      
    int Y,      
    int Z=-1    
    );

    void AddSpecialPin(
    int X,          
    int Y,          
    int Zone=-1     
    );

    void ForcedPinCloseToBorder(int PinZoneForItsDimension=-1){fForcedPinCloseToBorder=true;fZoneOfThePinsCloseToBorder=PinZoneForItsDimension;}
    void CreateGeometry();
//  void 



    void SetSource(
    int NPS=5000,               
    int ActivesCycles=100,      
    int InactivesCycles=100     
    );
    void UpdatePointer();
    void SetPlenum(double PlenumThickness){fPlenumThickness=PlenumThickness;fIsPlenum=true;}
    bool IsACOBRAObject(){return fIsCOBRA;}
    vector<vector<Cell*> > & GetFuelCell(){return fCFuels;} 



 protected:
  
    void Copy(const ReactorMesh &m);
    ReactorMesh(const ReactorMesh &m);
    void CreateShapes();
    void CreateCells();
    void FillSquaredCells();
    void FillTriangularCells();
    void CreatePlenum();
    bool ParTest(int n);

    int ZoneDiscretization(
    int i,      
    int j,      
    int k       
    );
    void Init();
    void PrintData();
    void Closing();
    void InitializeLatticeDiagram();
    void CreateGuideTubesCells();
    void CreateControlRodsCells();
    int CheckFillZoneOfPin(
    int i,          
    int j,          
    int k,          
    int FillZone    
    );


    int fAxialLevels;               
    int fRadialZones;               
    bool fIsHexagon;                
    bool fIsGasSpace;               
    bool fIsPinPos;                 
    bool fIsGuideTubes;             
    bool fIsControlRods;            
    bool fIsSpecialPins;            
    vector<double> fFuelRadius;     
    vector<double> fCladThickness;  
    vector<double> fGTRadius;       
    vector<double> fGTCladThickness;
    vector<double> fRadius1;        
    vector<double> fCRRadius;       
    vector<double> fCRCladThickness;
    vector<double> fCRCoolantGapThickness; 
    vector<double> fCRExtCladThickness;     
    double fPitch;                  
    double fChannelHeight;          
    double fAssemblyWidth;          
    double fAssemblyLength;         
    double fAssemblySide;           
    double fAssemblyHeight;         
    int fXNumberOfPins;             
    int fYNumberOfPins;             
    double fXDeviation;             
    double fYDeviation;             
    double fZDeviation;             
    double fXCenterDeviationOfMesh; 
    double fYCenterDeviationOfMesh; 
    vector<double> fGasThickness;   
    double fGasPressure;            
    double fHeFraMol;               
    double fXeFraMol;               
    int fXmin;                      
    int fXmax;                      
    int fYmin;                      
    int fYmax;                      
    int fZmin;                      
    int fZmax;                      
    vector<int> fXpos;              
    vector<int> fYpos;              
    vector<int> fZonepos;           
    vector<int> fGTXpos;            
    vector<int> fGTYpos;            
    vector<int> fGTZpos;            
    vector<int> fCRXpos;            
    vector<int> fCRYpos;            
    vector<int> fCRZpos;            
    vector<bool> fZonesOfGT;        
    vector<bool> fZonesOfCR;        
    vector<int> fSPXpos;            
    vector<int> fSPYpos;            
    vector<int> fZSP;               
    vector<bool> fZonesOfSP;        
    double fXTrans;                     
    double fYTrans;                     
    double fZTrans;                     



    Shape_ptr fSAssemblyZone;   
    Shape_ptr fSLatticeGen;     
    Cell* fCAssemblyZone;       
    Cell* fCLatticeGen;         
    vector<Material*> fMWhole;  
    Shape_ptr* fSWhole;         
    vector<Cell*> fCWhole;      
    vector<double> fCZR;        
    Shape_ptr fSPlenumUp;       
    Shape_ptr fSPlenumDown;     
    Material* fMPlenumUp;       
    Material* fMPlenumDown;     
    Cell* fCPlenumUp;           
    Cell* fCPlenumDown;         
    double fPlenumThickness;    



    vector<vector<Material*> > fMFuels;             
    vector<vector<Material*> > fMCladdings;         
    vector<vector<Material*> > fMCoolants;          
    vector<vector<Material*> > fMSpec1;             
    vector<vector<Material*> > fMGTInside;          
    vector<vector<Material*> > fMGTCladdings;       
    vector<vector<Material*> > fMGTOutside;         
    vector<vector<Material*> > fMCRMat;             
    vector<vector<Material*> > fMCRCladding;        
    vector<vector<Material*> > fMCRCoolantGap;      
    vector<vector<Material*> > fMCRExtCladding;     
    vector<vector<Material*> > fMCROutside;         
    vector<Shape_ptr*> fSFuels;                     
    vector<Shape_ptr*> fSCladding;                  
    vector<Shape_ptr*> fSCoolants;                  
    vector<Shape_ptr*> fSSpec1;                     
    vector<Shape_ptr*> fSGTInside;                  
    vector<Shape_ptr*> fSGTCladdings;               
    vector<Shape_ptr*> fSGTOutside;                 
    vector<Shape_ptr*> fSCRMat;           
    vector<Shape_ptr*> fSCRCladding;      
    vector<Shape_ptr*> fSCRCoolantGap;    
    vector<Shape_ptr*> fSCRExtCladding;       
    vector<Shape_ptr*> fSCROutside;       
    vector<vector<Cell*> > fCFuels;       
    vector<vector<Cell*> > fCCladdings;     
    vector<vector<Cell*> > fCCoolants;      
    vector<vector<Cell*> > fCSpec1;         
    vector<vector<Cell*> > fCGTInside;      
    vector<vector<Cell*> > fCGTCladdings;   
    vector<vector<Cell*> > fCGTOutside;     
    vector<vector<Cell*> > fCCRMat;           
    vector<vector<Cell*> > fCCRCladding;      
    vector<vector<Cell*> > fCCRCoolantGap;    
    vector<vector<Cell*> > fCCRExtCladding;       
    vector<vector<Cell*> > fCCROutside;       
    bool fInitVectors;                      


 
    string fAssemblyName;                           
    vector<int> fNumberOfPinInZone;                 
    vector<int> fNumberOfGTInZone;                  
    vector<int> fNumberOfCRInZone;                  
    vector<ofstream*> fFiles;                       
    vector<vector<vector<int> > > fLatticeZ;        
    string fExternalShapeType;                      
    bool fIsPlenum;                                 
    bool fIsCOBRA;                                  
    vector<bool> fAxialCrushSecurityHeterogeneity;  
    bool fForcedPinCloseToBorder;                   
    int fZoneOfThePinsCloseToBorder;                



};

#endif