Neutron Diffusion Equation Solver in Cartesian Coordinate Geometry
- DiffSolve-2D.cpp
- DESCRIPTION:steady state testbench
- REQUIREMENTS:DiffSolve-2D.hpp,ccomp,xsec
- BUGS: None
- AUTHOR: Alok_Jha
- E-mail: jhaalok1984_at_gmail.com
- ORGANIZATION:NPCIL
- VERSION:1.1
- REVISION:0.0
- CREATED:Thursday, 17 January 2019. Released on 18-06-2019
None
Any 2-D Core size input will work. Modify in the header file. However, for WINDOWS machine of 32 bits, the core size smaller than 12 rows by 16 columns, the results are not correct due to 'double' and 'long double' limitations. Please check by reducing the
#define F_max 1e5 //maximum iteration
#define K_max 1e5 //maximum iteration
in the header file. It may yield reasonably acceptable results.
- Added reactivity ρ calculation in the code.
- Added criticality declaration in the code.
Use makefile
$ make -k makefile
or, at the terminal
Call
$ g++ DiffSolve-2D.cpp -o DiffSolve-2D && ./DiffSolve-2D
or,
$ c++ DiffSolve-2D.cpp -o DiffSolve-2D && ./DiffSolve-2D
$ gnuplot < plot.gp
Label in the plot "thermal.txt" corrected from the earlier code where incorrect y-label is generated.
Reach out to me at the following place!
- Twitter at
@amigo_2222