jay
/
fast-nuclear-sim
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

made tweaks

This commit is contained in:
Jay Robson 2024-01-15 17:15:08 +11:00
parent b9e5b00122
commit 929b14a21d
3 changed files with 12 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
src/reactor/fuel/sample.cpp
View File

@ -66,7 +66,7 @@ void sample::update(double secs)
// sim::reactor::fuelulate fuel use // sim::reactor::fuelulate fuel use
fuel -= neutrons_fuel; fuel -= neutrons_fuel;
energy += neutrons_fuel; energy += neutrons_fuel;
waste.add_fissile(neutrons_fuel * 5); waste.add_fissile(neutrons_fuel * 6);
// do the poison // do the poison
te_135 += neutrons_fuel * (1.0 / 8.0); te_135 += neutrons_fuel * (1.0 / 8.0);

2
src/reactor/fuel/sample.hpp
View File

@ -10,7 +10,7 @@ namespace sim::reactor::fuel
class sample class sample
{ {
constexpr static const double Xe_135_M = 1e3; constexpr static const double Xe_135_M = 1e4;
sim::reactor::fuel::waste waste; sim::reactor::fuel::waste waste;

11
src/reactor/rod.cpp
View File

@ -17,7 +17,16 @@ void rod::add(val_t type, double v)
double rod::extract(val_t type, double s, double k, double o) double rod::extract(val_t type, double s, double k, double o)
{ {
double v = (1 - std::pow(0.5, s * -std::log2(1 - k * get_k(type)))) * 0.5 * (get(type) - o); k *= get_k(type);
double m = 1;
if(k < 1)
{
m = 1 - std::pow(0.5, s * -std::log2(1 - k * get_k(type)));
}
double v = m * 0.5 * (get(type) - o);
vals_in[type] -= v; vals_in[type] -= v;
return v; return v;
} }