1d diffusion equation solution for slab with non symmetric source

Click For Summary
SUMMARY

The discussion focuses on solving the 1D diffusion equation for a slab with a non-symmetric source distribution defined as Se^(-k(x+a/2)). The user has derived both homogeneous and particular solutions but struggles with boundary conditions, specifically setting the neutron flux to zero at -a/2 and a/2. The conversation highlights the importance of applying the condition that "there always exists a physically realizable solution" and suggests exploring Green's functions for convolution with the source term as a potential method for finding a solution.

PREREQUISITES
  • Understanding of 1D diffusion equations
  • Familiarity with boundary conditions in differential equations
  • Knowledge of Green's functions and their applications
  • Basic concepts of neutron flux in nuclear physics
NEXT STEPS
  • Research the application of Green's functions in solving differential equations
  • Study the implications of boundary conditions in diffusion problems
  • Explore the concept of physically realizable solutions in nuclear reactor physics
  • Review analytical techniques for solving non-homogeneous differential equations
USEFUL FOR

Students and professionals in nuclear engineering, particularly those dealing with diffusion equations and neutron transport theory, will benefit from this discussion.

Mojo
Messages
1
Reaction score
0
Thread moved from the technical forums, so no HH Template is shown.
Disclaimer: This is a homework problem

I need to analytically solve the diffusion equation for a 1d 1 group slab with width a, and source distribution Se^(-k(x+a/2))

I've gone through the math, and come up with my homogeneous and particular solution and attempted to apply the boundary conditions of the flux set to 0 at -a/2 and a/2 but my constants are functions of each other. I was thinking of solving the source for the average, and setting the neutron current at that point to 0 but I didn't think that would work. Is there any other boundary condition I can use? I read in another text (I went through about 4 texts and 7 power points to no avail) "There always exists a physically realizable solution (no critical buckling!)") (http://www.mit.edu/~lululi/school/22.211_Nuclear_Reactor_Physics_I/notes/__all__.pdf) But I am not sure how to use this condition, or if it applies to this problem. I would greatly appreciate any insight into this.
 
Physics news on Phys.org
If you haven't tried it already, I would try finding a Green's function and taking it's convolution with the source term.
 

Similar threads

  • · Replies 11 ·
Replies
11
Views
3K
Replies
0
Views
1K
  • · Replies 2 ·
Replies
2
Views
4K
  • · Replies 3 ·
Replies
3
Views
3K
Replies
1
Views
2K
  • · Replies 8 ·
Replies
8
Views
4K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 9 ·
Replies
9
Views
3K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K