Criticality calculation of an homogeneous finite reactor

Babatunde22

Please,I am working on the criticality calculation of an homogeneous finite cylindrical reactor core using four-group diffusion equations. I have been able to discretize the multigroup diffusion equations using the finite difference method(FED). But

I am stocked on the iterative method to solve the linear equations using matlab. Thanks

Related Nuclear Engineering News on Phys.org

Astronuc

Staff Emeritus
I'm assuming one is doing this in a 2D geometry.

Has one studied an example of a 2D 2-group system?

Babatunde22

yes sir, it is a 2D geometry which entails r and z, but most examples i have come by are mostly 1D and very few 2D(x and y) 2-group system

Babatunde22

please, I wanted to upload my sets of equations on the forum, but I couldn't. How can I display my set of discretized equations

Staff Emeritus

Attachments

• 392 KB Views: 110
• 392 KB Views: 126

Babatunde22

thank you sir, I have been able to attached the document

Babatunde22

please I am stocked on how to use an iterative method(gauss-seidel or SOR method) to solve the linear equations using matlab with the above attached pdf document. Thanks

NukeProf

I think you will have to program the Gauss-Seidel method yourself. I recall that it was straightforward in FORTRAN, but should be doable in MATLAB.

Babatunde22

Yes I want to programme it using the Gauss- Seidel method, but i need an algorithm or the calculation strategy well explicit that I can use to calculate the criticality calculation. Thanks

Physics Forums Values

We Value Quality
• Topics based on mainstream science
• Proper English grammar and spelling
We Value Civility
• Positive and compassionate attitudes
• Patience while debating
We Value Productivity
• Disciplined to remain on-topic
• Recognition of own weaknesses
• Solo and co-op problem solving