Neutron Flux in a sub-critical reactor

AI Thread Summary
The discussion focuses on deriving the neutron flux equation for a sub-critical reactor with a specific configuration of nested spheres. The user is struggling with the derivation of the equation Φ(r) = S/4πD [1/r - 1/R2], particularly the transition from the diffusion equation to the desired form. The reactor consists of a lead target, actinide fuel, and a reflector, with neutron flux vanishing at the outer radius of the reflector. The user notes that their attempts result in exponential terms, indicating a challenge in applying the point source assumption correctly. Assistance is sought to clarify the derivation process while considering symmetry and current continuity at the interfaces.
Oxlade
Messages
5
Reaction score
0
Hi everyone,

I am supposed to derive the neutron flux equation provided for region A of my reactor. Just wondering if anyone can help me out since I stuck on the derivation for [1/r - 1/R2]; S/4πD aspect is very similar to a solving the constant for a point source spherical reactor

Here is my reactor diagram consisting of nested spheres (sorry, this is my first time posting):

Sphere 1 --> Lead target with radius of R0 (Target-region)
Sphere 2 --> Actinide Fuel surrounding lead target; has outer radius of R1 (A-region)
Sphere 3 --> Reflector surrounding the fuel with outer radius of R2 (E-region)

Neutron Flux vanished at R2 and the flux at A-region is given as:

Φ(r) = S/4πD [1/r - 1/R2]My attempts at deriving this equation always leaves me with exponential terms if I'm assuming I'm solving the diffusion equation for a finite spherical reactor with a point source. Can anyone help me?
 
Engineering news on Phys.org
Assuming symmetry, the current at the center of the Pb is zero. The current at the boundaries/interfaces on one side must equal that on the other side.
 
Hello everyone, I am currently working on a burnup calculation for a fuel assembly with repeated geometric structures using MCNP6. I have defined two materials (Material 1 and Material 2) which are actually the same material but located in different positions. However, after running the calculation with the BURN card, I am encountering an issue where all burnup information(power fraction(Initial input is 1,but output file is 0), burnup, mass, etc.) for Material 2 is zero, while Material 1...
Hi everyone, I'm a complete beginner with MCNP and trying to learn how to perform burnup calculations. Right now, I'm feeling a bit lost and not sure where to start. I found the OECD-NEA Burnup Credit Calculational Criticality Benchmark (Phase I-B) and was wondering if anyone has worked through this specific benchmark using MCNP6? If so, would you be willing to share your MCNP input file for it? Seeing an actual working example would be incredibly helpful for my learning. I'd be really...
Back
Top