Register to reply 
Diffusion equation and neutron diffusion theory 
Share this thread: 
#1
Mar706, 09:37 AM

Admin
P: 21,865

Basically the steadystate diffusion equation can be written in a form
[itex] \nabla^2\phi\,+\,k^2\phi\,=\,S[/itex] When S = 0, this is just the Helmholtz equation  http://mathworld.wolfram.com/Helmhol...lEquation.html See also  http://www.math.ohiostate.edu/~gerl...t/node107.html To solve it, like any differential equation, one simply applies the boundary conditions. I will add more later. 


#2
Mar1006, 02:55 AM

P: 2

What about the 2 group diffusion equation? If you are familiar, is said equation only solvable by computer program?



#3
Mar1006, 06:18 AM

Admin
P: 21,865

In special cases (e.g. 1D), there can be an analytical solution for two group diffusion theory, basically solving two simultaneous linear differential equations.
However, for most practical (realworld) cases, the two group diffusion theory requires a numerical solution. Most modern nuclear design codes use a modified 2 group approach. There has been some effort at employing transport theory, but it has proved difficult. 


#4
May206, 07:52 AM

Admin
P: 21,865

Diffusion equation and neutron diffusion theory
A brief overview of transport and diffusion theory with respect to neutron propagation.
http://lpsc.in2p3.fr/gpr/PPNPport/node28.html This is not comprehensive  just a summary. The diffusion equation (second Boltzmann equation) is a special case of the transport equation. 


#5
May206, 09:01 AM

P: 395

Just in time for my final today!



#6
May206, 09:36 AM

P: 1,157

Weird  I've been solving this all day in cylindrical polars...
...subject to some nasty conditions 


#7
May206, 09:40 AM

Admin
P: 21,865

Yep  Bessel's functions are the solutions to Bessel's equation, which is the form of the diffusion equation in cylindrical or polar coordinates for radial dependence.



#8
Sep1006, 08:53 PM

P: 6

Hi, I just started graduate level transport theory (sigh, the transport equation) and I'm having a heck of a time with it. Can anyone recommend any fantastic texts or websites? (I very obviously didn't do my undergrad in NE) Thanks!



#9
Sep1106, 12:05 AM

Sci Advisor
PF Gold
P: 1,479




#10
Sep1106, 09:27 PM

P: 395

I think "k" is [tex]\frac {1} {L^2}[/tex] where [tex]L^2[/tex] is diffusion area. I think that makes "S" the source term.
Edit: I am having some problems with latex. Does the above look alright to everyone? 


#11
Sep1206, 08:07 AM

Emeritus
PF Gold
P: 8,147

If you mean the way it displays on white, that is a feature of our installation that showed up the last time we did an upgrade to our main software. It's on the list of problems to be looked at, but other problems have higher priority; people can read this just fine, after all. 


#12
Sep1206, 08:09 AM

Admin
P: 21,865

The diffusion equation is derived from transport theory with several assumptions.
Basically the steadystate neutron diffusion equation can be written as: [tex]D\,\nabla^2\phi\,\,\Sigma_a\,\phi\,+\,S\,=\,0[/tex], where D is the diffusion coefficient, here spatially independent, i.e. constant, and [itex]\Sigma_a[/itex] is the macroscopic absorption coefficient, and [itex]\phi[/itex] is the flux. Buried in here is an assumption that the neutrons are more or less the same energy. If not the case, then one must account for different energy groups and the absorption coefficient becomes a removal coefficient which includes absorption and scattering out of the energy group. Anyway, the above equation becomes, [tex]D\,\nabla^2\phi\,\,\Sigma_a\,\phi\,=\,S[/tex], and dividing by D [tex]\nabla^2\phi\,\,\frac{\Sigma_a}{D}\,\phi\,=\,\frac{S}{D}[/tex], or [tex]\nabla^2\phi\,\,\frac{1}{L^2}\,\phi\,=\,\frac{S}{D}[/tex], where [tex]L^2\,=\,\frac{D}{\Sigma_a}[/tex] Then there is the case where S = a function of the flux [itex]\phi[/itex]. 


#13
Nov1007, 08:47 AM

P: 15

[tex]D_{2}d^{2}\phi_{1}[/tex]/dx[tex]^{2}[/tex]+[tex]\Sigma_{R1}\phi_{1}[/tex]=1/k(v[tex]_{1}[/tex][tex]\Sigma_{f1}\phi_{1}[/tex]+v[tex]_{2}[/tex][tex]\Sigma_{f2}\phi_{2}[/tex]) [tex]D_{2}d^{2}\phi_{2}[/tex]/dx[tex]^{2}[/tex]+[tex]\Sigma_{a2}\phi_{2}[/tex]=[tex]\Sigma_{s12}\phi_{1}[/tex] This is solving a two group diffusion in a slab reactor with plane source at the center of the subcritical slab. So analytical or numerical? Any guess on the method to solve for flux? Sorry, but v1 and v2 are not squared. They are just v_1 and v_2. Some reason it puts as a square. 


#14
Nov1007, 08:56 AM

Admin
P: 21,865

Certainly the problem can be solved numerically (FD or FE), and I believe analytically, but I'd have to dig back in my archives for that.
one could write \nu_1 and \nu_2 in thex LaTeX expressions before \Sigma. I think this is how the equations are supposed to look: [tex]{D_1}\frac{{d^2}\phi_1}{dx^2}\,+\,\Sigma_{R1}\phi_1\,=\,\frac{1}{k}({\n u_1}{\Sigma_{f1}\phi_1}\,+\,{\nu_2}{\Sigma_{f2}\phi_2})[/tex] [tex]{D_2}\frac{{d^2}\phi_2}{dx^2}\,+\,\Sigma_{a2}\phi_2\,=\,{\Sigma_{s12}}\ phi_1[/tex] Just looking these, one could collect coefficents and rewrite the equations as: [itex]\phi_1[/itex]'' + A [itex]\phi_1[/itex] = B [itex]\phi_2[/itex] [itex]\phi_2[/itex]'' + C [itex]\phi_2[/itex] = D [itex]\phi_1[/itex] 


#15
Nov1207, 06:31 PM

Admin
P: 21,865




#16
Oct2210, 04:41 PM

P: 1

this forum has helped with this challenge problem but I am still having a hard time getting started, the problem is:
"There is a waiting room on the opposite side of a very large wall adjacent to a proton therapy treatment room at the local hospital. Compute the neutron flux [tex]\phi[/tex](x) into the waiting room using a diffusion theory approximation. Assume the neutrons are emitted from the wall surface via a uniform planar surface source emitting Snot neutrons/cm^2/s. ( At wall surface, x=0) The diffusion equation is D*(d^2[tex]\phi[/tex]/dx^2)[tex]\sigma_{a}[/tex]=0 for x not equal to 0. Assuming a 1dimensional flux approximation, and other dimensions of the room (y,z) relative to the wall (at x=0) are infinite, use the following conditions: (i)Flux must remain finite as x> infinity (ii)The Xray all current(Neutron coming out of the wall) has a limit as x>0, where: limx>0 J(x)*i^>limx>0(D(d[tex]\phi[/tex]/dx))=(Snot/2) (iii) D, [tex]\sigma_{a}[/tex], Snot are all constants" any help is appreciated! 


Register to reply 
Related Discussions  
What is neutron diffusion?  High Energy, Nuclear, Particle Physics  15  
Neutron diffusion problem.  Nuclear Engineering  16  
Diffusion equation in 1D  Introductory Physics Homework  4  
Diffusion equation  Calculus & Beyond Homework  0  
Neutron diffusion and nucleosynthesis  Astronomy & Astrophysics  2 