Criticality calculation of an homogeneous finite reactor

Click For Summary

Discussion Overview

The discussion focuses on the criticality calculation of a homogeneous finite cylindrical reactor core using four-group diffusion equations. Participants are exploring the application of iterative methods to solve linear equations derived from discretized multigroup diffusion equations, specifically in a 2D geometry.

Discussion Character

  • Technical explanation
  • Mathematical reasoning
  • Homework-related

Main Points Raised

  • One participant is working on criticality calculations using four-group diffusion equations and has discretized these equations using the finite difference method but is struggling with the iterative method in MATLAB.
  • Another participant assumes the calculations are being done in a 2D geometry and inquires about the existence of examples for a 2D 2-group system.
  • A participant confirms the 2D geometry involves r and z coordinates but notes that most available examples are 1D or limited 2D systems.
  • One participant expresses difficulty in uploading their set of discretized equations to the forum.
  • Another participant suggests using LaTeX or creating an image or PDF to upload the equations.
  • A participant mentions they have successfully attached a document containing their equations.
  • One participant seeks guidance on using an iterative method (Gauss-Seidel or SOR) in MATLAB to solve the linear equations related to their attached document.
  • Another participant recalls that programming the Gauss-Seidel method was straightforward in FORTRAN and suggests it should be doable in MATLAB.
  • A participant requests an explicit algorithm or calculation strategy for implementing the Gauss-Seidel method for criticality calculations.

Areas of Agreement / Disagreement

Participants do not appear to reach a consensus on the best approach for implementing the iterative method in MATLAB, and there are varying levels of experience and resources shared regarding the programming of the Gauss-Seidel method.

Contextual Notes

Participants express uncertainty regarding the availability of 2D examples and the specific details of the iterative methods, indicating a reliance on personal experience and shared resources.

Babatunde22
Messages
7
Reaction score
0
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
 
Engineering news on Phys.org
I'm assuming one is doing this in a 2D geometry.

Has one studied an example of a 2D 2-group system?
 
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
 
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
 
PLEASE FIND ATTACHED
 

Attachments

thank you sir, I have been able to attached the document
 
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
 
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.
 
  • #10
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
 

Similar threads

Replies
3
Views
3K
Replies
4
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 31 ·
2
Replies
31
Views
4K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 3 ·
Replies
3
Views
5K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 0 ·
Replies
0
Views
2K