Solving Square Box Reactor Question: Flux General Equation

  • Thread starter Thread starter matt222
  • Start date Start date
  • Tags Tags
    Box Reactor
AI Thread Summary
The discussion focuses on deriving the general equation for flux in a square box, leading to the equation flux = A cos(pi*x/a) cos(pi*y/a) cos(pi*z/a), where 'a' is the side length. It confirms that this equation is valid under the condition of zero flux at the boundaries and maximum flux at the center. The conversation also touches on the neutron diffusion equation, suggesting that the general solution can be expressed as flux(x,y,z) = S/(D[3pi^2/a^2 + 1/L^2]). Additionally, it raises the question of how to express the flux in a steady-state scenario without a source, proposing a solution of flux(x,y,z) = A exp(-x/a) + B exp(x/a). The discussion ultimately aims to clarify the mathematical representation of flux in these contexts.
matt222
Messages
130
Reaction score
0
Iam trying to work out question related to sequare box, my question is what the flux general equation for the sequare box?

I work it out and found it is

flux=Acos(pi*x/a)cos(pi*y/a)cos(pi*z/a)

where a is side lenght, am i right
 
Engineering news on Phys.org
matt222 said:
Iam trying to work out question related to sequare box, my question is what the flux general equation for the sequare box?

I work it out and found it is

flux=Acos(pi*x/a)cos(pi*y/a)cos(pi*z/a)

where a is side length, am i right
Yes - with -a/2 < x < a/2, and similarly for y and z.

This assumes zero flux boundary condition at the edges (flux(-a/2) = flux(a/2) = 0), and max flux at flux(x,y,z) = flux(0,0,0) = A and even symmetry.

Flux is then described by 3 independent functions X, Y, Z which are described by X'' + (pi/a)2X = 0 and X(0) = Y(0) = Z(0) = A1/3
 
do you agree with me when we use neutron diffusion equation, that the general solution will be

flux(x,y,z)=S/(D[3pi^2/a^2+1/L^2])


and what about if we have no source and steady state? How the general equation going to be, should we take it as

flux(x,y,z)=Aexp(-x/a)+Bexp(x/a)
 
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...
Back
Top