How do I calculate cascade volume in a steel alloy under ion irradiation?

  • Thread starter Thread starter nukapprentice
  • Start date Start date
  • Tags Tags
    Collision
AI Thread Summary
The discussion focuses on calculating cascade volume in a steel alloy under 1 MeV ion irradiation, specifically using a PKA energy of around 10 keV. The user expresses uncertainty about whether to start calculations from scratch or utilize existing simulation tools like Geant4. They mention deriving a cascade volume analytically but find it unexpectedly small at 50 nm^3. The user also notes that using Averback's equation from TRIM results may not be beneficial for their specific PKA energy of 14 keV. Assistance in finding Matlab code for this calculation is requested, highlighting the need for references on existing methodologies.
nukapprentice
Messages
69
Reaction score
0
So I was wondering if anyone could help me out. I am trying to find a way to calculate the cascade volume in a steel alloy undergoing ~ 1 MeV ion irradiation ( Z ~ 25). For this analysis I will be looking at the electronic stopping regime with PKA energy of ~ 10 keV. Not sure what to do after calculating the PKA energy.
 
Do you have to do it from scratch, or can you check what common simulation tools like Geant4 are doing?
 
I have to do it from scratch since my programming abilities are very limited. I actually came up with a volume analytically using a maximum cascade temperature of ~10000 K, but the volume seems small ( V = 50 nm^3 ). Since I am trying to calculate the volume mid ion range, (about halfway through) I think using Averback's equation from TRIM results will not be helpful. My PKA energy is actually 14 keV which leads to approximately 10 keV for damage energy. You wouldn't happen to know of any Matlab code to handle this would you?
 
I would use some existing software package, or at least see what they did (there should be some references explaining that, I just don't know where).
 
Ok, thanks for the replies.
 
Hello, I'm currently trying to compare theoretical results with an MCNP simulation. I'm using two discrete sets of data, intensity (probability) and linear attenuation coefficient, both functions of energy, to produce an attenuated energy spectrum after x-rays have passed through a thin layer of lead. I've been running through the calculations and I'm getting a higher average attenuated energy (~74 keV) than initial average energy (~33 keV). My guess is I'm doing something wrong somewhere...
Back
Top