Skip to content

Latest commit

 

History

History
99 lines (61 loc) · 2.43 KB

readme.md

File metadata and controls

99 lines (61 loc) · 2.43 KB
Raw

2-D Neutron Diffusion Equation Solver in Cartesian Coordinate Geometry

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

Build Status License Badges

Table of Contents (Optional)

Installation

None

Features

Any Core Size

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.

Addendum

  • Added reactivity ρ calculation in the code.
  • Added criticality declaration in the code.

Compile

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

Plotting Of Results

$ gnuplot < plot.gp

Label in the plot "thermal.txt" corrected from the earlier code where incorrect y-label is generated.

Support

Reach out to me at the following place!

License

License