Solve Diffusion Equation for Neutron Flux in Multiplying Sphere

  • Thread starter Thread starter EngNewbit
  • Start date Start date
  • Tags Tags
    Diffusion Theory
AI Thread Summary
The discussion focuses on solving the diffusion equation for neutron flux in a multiplying sphere with a constant source strength. The equation must be expressed in spherical coordinates, incorporating boundary conditions where the flux is finite at the origin and zero at the sphere's surface. The solution for the neutron flux when k inf equals 1 is derived as (Q/6D)(R^2 − r^2). Additionally, it is emphasized that the current should also be zero at the origin. Understanding these concepts is crucial for progressing in the set of problems presented.
EngNewbit
Messages
1
Reaction score
0
This is in the beginning of a long set of problems, and I am lost. I don't get anything like this answer. Any guidance? I have a feeling its simple but haven't done much of these.

Write down the diffusion equation for the neutron flux in a multiplying sphere of radius R containing a constant distributed source of strength Q neutrons/cm3/s. Assuming that the flux vanishes at the sphere surface and that it remains finite at the origin:
Solve the diffusion equation when k inf = 1 and show that the neutron flux is given by:
(Q/6D)(R^2 − r^2)
 
Engineering news on Phys.org
One has to write the diffusion equation (in spherical coordinates) with a constant, distributed source, and apply the boundary conditions, that the flux \phi is finite at r = 0, and 0 at r = R.

The current should also be 0 at r = 0.
 
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