decay_heat_with_nuclear 0.0.2

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: cefcf0711745f3546c8b0965529aa4ea97a0c9f1
+  data.tar.gz: e505da9206e1da0d7746dc3ad9f330ecf8bf4144
+SHA512:
+  metadata.gz: 719356ded9a0b1b548164f6fea5deb2b7e6cdbf7b7a5df41750c14b245694dde809b9d3f99860b52a6a6710d001e408ce17732933992ab6da727d68f9400fe22
+  data.tar.gz: 645af1f55684aa918b610cdecca5d0bed1b4f25f5273b22d356ecbd483a612fe26534fddca67b81583738a0542a13f2f8d25adc388d9f94e2a21ca1ce589667e

data/.gitignore

@@ -0,0 +1,35 @@
+*.gem
+*.rbc
+/.config
+/coverage/
+/InstalledFiles
+/pkg/
+/spec/reports/
+/test/tmp/
+/test/version_tmp/
+/tmp/
+
+## Specific to RubyMotion:
+.dat*
+.repl_history
+build/
+
+## Documentation cache and generated files:
+/.yardoc/
+/_yardoc/
+/doc/
+/rdoc/
+
+## Environment normalisation:
+/.bundle/
+/vendor/bundle
+/lib/bundler/man/
+
+# for a library or gem, you might want to ignore these files since the code is
+# intended to run in multiple environments; otherwise, check them in:
+# Gemfile.lock
+# .ruby-version
+# .ruby-gemset
+
+# unless supporting rvm < 1.11.0 or doing something fancy, ignore this:
+.rvmrc

data/Gemfile

@@ -0,0 +1,4 @@
+source 'https://rubygems.org'
+
+# Specify your gem's dependencies in decay_heat_with_nuclear.gemspec
+gemspec

data/LICENSE.txt

@@ -0,0 +1,166 @@
+                   GNU LESSER GENERAL PUBLIC LICENSE
+                       Version 3, 29 June 2007
+
+ Copyright (C) 2007 Free Software Foundation, Inc. <http://fsf.org/>
+ Everyone is permitted to copy and distribute verbatim copies
+ of this license document, but changing it is not allowed.
+
+
+  This version of the GNU Lesser General Public License incorporates
+the terms and conditions of version 3 of the GNU General Public
+License, supplemented by the additional permissions listed below.
+
+  0. Additional Definitions.
+
+  As used herein, "this License" refers to version 3 of the GNU Lesser
+General Public License, and the "GNU GPL" refers to version 3 of the GNU
+General Public License.
+
+  "The Library" refers to a covered work governed by this License,
+other than an Application or a Combined Work as defined below.
+
+  An "Application" is any work that makes use of an interface provided
+by the Library, but which is not otherwise based on the Library.
+Defining a subclass of a class defined by the Library is deemed a mode
+of using an interface provided by the Library.
+
+  A "Combined Work" is a work produced by combining or linking an
+Application with the Library.  The particular version of the Library
+with which the Combined Work was made is also called the "Linked
+Version".
+
+  The "Minimal Corresponding Source" for a Combined Work means the
+Corresponding Source for the Combined Work, excluding any source code
+for portions of the Combined Work that, considered in isolation, are
+based on the Application, and not on the Linked Version.
+
+  The "Corresponding Application Code" for a Combined Work means the
+object code and/or source code for the Application, including any data
+and utility programs needed for reproducing the Combined Work from the
+Application, but excluding the System Libraries of the Combined Work.
+
+  1. Exception to Section 3 of the GNU GPL.
+
+  You may convey a covered work under sections 3 and 4 of this License
+without being bound by section 3 of the GNU GPL.
+
+  2. Conveying Modified Versions.
+
+  If you modify a copy of the Library, and, in your modifications, a
+facility refers to a function or data to be supplied by an Application
+that uses the facility (other than as an argument passed when the
+facility is invoked), then you may convey a copy of the modified
+version:
+
+   a) under this License, provided that you make a good faith effort to
+   ensure that, in the event an Application does not supply the
+   function or data, the facility still operates, and performs
+   whatever part of its purpose remains meaningful, or
+
+   b) under the GNU GPL, with none of the additional permissions of
+   this License applicable to that copy.
+
+  3. Object Code Incorporating Material from Library Header Files.
+
+  The object code form of an Application may incorporate material from
+a header file that is part of the Library.  You may convey such object
+code under terms of your choice, provided that, if the incorporated
+material is not limited to numerical parameters, data structure
+layouts and accessors, or small macros, inline functions and templates
+(ten or fewer lines in length), you do both of the following:
+
+   a) Give prominent notice with each copy of the object code that the
+   Library is used in it and that the Library and its use are
+   covered by this License.
+
+   b) Accompany the object code with a copy of the GNU GPL and this license
+   document.
+
+  4. Combined Works.
+
+  You may convey a Combined Work under terms of your choice that,
+taken together, effectively do not restrict modification of the
+portions of the Library contained in the Combined Work and reverse
+engineering for debugging such modifications, if you also do each of
+the following:
+
+   a) Give prominent notice with each copy of the Combined Work that
+   the Library is used in it and that the Library and its use are
+   covered by this License.
+
+   b) Accompany the Combined Work with a copy of the GNU GPL and this license
+   document.
+
+   c) For a Combined Work that displays copyright notices during
+   execution, include the copyright notice for the Library among
+   these notices, as well as a reference directing the user to the
+   copies of the GNU GPL and this license document.
+
+   d) Do one of the following:
+
+       0) Convey the Minimal Corresponding Source under the terms of this
+       License, and the Corresponding Application Code in a form
+       suitable for, and under terms that permit, the user to
+       recombine or relink the Application with a modified version of
+       the Linked Version to produce a modified Combined Work, in the
+       manner specified by section 6 of the GNU GPL for conveying
+       Corresponding Source.
+
+       1) Use a suitable shared library mechanism for linking with the
+       Library.  A suitable mechanism is one that (a) uses at run time
+       a copy of the Library already present on the user's computer
+       system, and (b) will operate properly with a modified version
+       of the Library that is interface-compatible with the Linked
+       Version.
+
+   e) Provide Installation Information, but only if you would otherwise
+   be required to provide such information under section 6 of the
+   GNU GPL, and only to the extent that such information is
+   necessary to install and execute a modified version of the
+   Combined Work produced by recombining or relinking the
+   Application with a modified version of the Linked Version. (If
+   you use option 4d0, the Installation Information must accompany
+   the Minimal Corresponding Source and Corresponding Application
+   Code. If you use option 4d1, you must provide the Installation
+   Information in the manner specified by section 6 of the GNU GPL
+   for conveying Corresponding Source.)
+
+  5. Combined Libraries.
+
+  You may place library facilities that are a work based on the
+Library side by side in a single library together with other library
+facilities that are not Applications and are not covered by this
+License, and convey such a combined library under terms of your
+choice, if you do both of the following:
+
+   a) Accompany the combined library with a copy of the same work based
+   on the Library, uncombined with any other library facilities,
+   conveyed under the terms of this License.
+
+   b) Give prominent notice with the combined library that part of it
+   is a work based on the Library, and explaining where to find the
+   accompanying uncombined form of the same work.
+
+  6. Revised Versions of the GNU Lesser General Public License.
+
+  The Free Software Foundation may publish revised and/or new versions
+of the GNU Lesser General Public License from time to time. Such new
+versions will be similar in spirit to the present version, but may
+differ in detail to address new problems or concerns.
+
+  Each version is given a distinguishing version number. If the
+Library as you received it specifies that a certain numbered version
+of the GNU Lesser General Public License "or any later version"
+applies to it, you have the option of following the terms and
+conditions either of that published version or of any later version
+published by the Free Software Foundation. If the Library as you
+received it does not specify a version number of the GNU Lesser
+General Public License, you may choose any version of the GNU Lesser
+General Public License ever published by the Free Software Foundation.
+
+  If the Library as you received it specifies that a proxy can decide
+whether future versions of the GNU Lesser General Public License shall
+apply, that proxy's public statement of acceptance of any version is
+permanent authorization for you to choose that version for the
+Library.
+

data/README.md

@@ -0,0 +1,44 @@
+# decay_heat_with_nuclear
+
+It's a easy way to get nuclear fuel decay heat when it remove from reactor core.
+Using 3 method `ASB9-2`, `ANS-1973` and `ANS-1979` to calculate it.
+
+## Install this gem
+
+```
+gem install decay_heat_with_nuclear
+```
+
+
+## Yor need a HASH
+
+The hash will like this:
+
+```
+hash = {
+          ts: array,
+          t0: array
+       }
+```
+
+where
+ts is `Time after remove (sec)`,
+t0 is `Cumulative reactor operating time (sec)`
+
+
+## Calculate from this gem
+
+```
+DecayHeatWithNuclear.run(hash)
+```
+
+and you will get the hash data like:
+
+```
+output_hash = {
+                ans1979: { ts: array, P/P0: array },
+                ans1973: { ts: array, P/P0: array },
+                asb9_2:  { ts: array, P/P0: array },
+              }
+
+```

data/Rakefile

@@ -0,0 +1,2 @@
+require "bundler/gem_tasks"
+

data/decay_heat_with_nuclear.gemspec

@@ -0,0 +1,23 @@
+# coding: utf-8
+lib = File.expand_path('../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'decay_heat_with_nuclear/version'
+
+Gem::Specification.new do |spec|
+  spec.name          = 'decay_heat_with_nuclear'
+  spec.version       = DecayHeatWithNuclear::VERSION
+  spec.authors       = ['C.lin']
+  spec.email         = ['demeter.yeh@gmail.com']
+  spec.summary       = %q{It's a easy way to get nuclear fuel decay heat when it remove from reactor core.}
+  spec.description   = %q{It's a easy way to get nuclear fuel decay heat when it remove from reactor core.}
+  spec.homepage      = ''
+  spec.license       = 'GPL3.0'
+
+  spec.files         = `git ls-files -z`.split("\x0")
+  spec.executables   = spec.files.grep(%r{^bin/}) { |f| File.basename(f) }
+  spec.test_files    = spec.files.grep(%r{^(test|spec|features)/})
+  spec.require_paths = ['lib']
+
+  spec.add_development_dependency 'bundler', '~> 1.6'
+  spec.add_development_dependency 'rake', '~> 10.0'
+end

data/lib/decay_heat_with_nuclear/main_run.rb

@@ -0,0 +1,338 @@
+module MainRun
+  class RunInit
+    attr_accessor :ts, :t0, :power
+
+    def initialize(hash_data, option = 0, file_name = 'default.txt')
+      @ts = hash_data[:ts]
+      @t0 = hash_data[:t0]
+      @power = []
+      @option = option
+      @file_name = file_name
+      @output_hash = {
+                        ts:             [],
+                        p_p0:           [],
+                        p_p0_without_k: []
+                     }
+    end
+
+    def dataout(ts, p_p0, p_p0_without_k = nil)
+      @output_hash[:ts] << ts
+      @output_hash[:p_p0] << p_p0
+      @output_hash[:p_p0_without_k] << p_p0_without_k if p_p0_without_k
+    end
+
+    # day2sec(day)
+    def day2sec(day)
+      if day < 86400 * 10
+        day * 24 * 3600.0
+      end
+    end
+
+    # sec2day(sec)
+    def sec2day(sec)
+      sec / 24 / 3600.0
+    end
+  end
+
+  class RunAns1973 < RunInit
+    def run(ts = @ts, t0 = @t0, power = @power)
+      read_data = ThermalData::DataForANS_5_1_1973.new
+      ts.each_index do |i|
+        ts_to_f = ts[i].to_f
+        t0_to_f = t0[i].to_f
+        # power_to_f = power[i].to_f
+
+        t0_add_ts   = ts_to_f + t0_to_f
+
+        p_p0_source = calc_thermal_fission_functions(t0_to_f, ts_to_f, read_data)
+        p_p0        = calc_sum_thermal_fission(p_p0_source)
+
+        p_p0_U239   = calc_thermal_fission_functions_with_U239(t0_to_f, ts_to_f)
+        p_p0_Np239  = calc_thermal_fission_functions_with_Np239(t0_to_f, ts_to_f)
+
+        p_p0_tatal    = p_p0 + p_p0_U239 + p_p0_Np239
+
+        dataout(ts_to_f, p_p0_tatal)
+      end
+      @output_hash
+    end
+
+    private
+
+    # Calculate thermal fission functions from ASB9-2.
+    # need t0, ts and data with thermal fission(from class DataForASB_9_2)
+    # ts: Time after remove (sec)
+    # t0: Cumulative reactor operating time (sec)
+    #
+    # calc_thermal_fission_functions(t0, ts, read_data)
+    #
+    # return { :ts          => P/P0(t_inf, ts)
+    #          :ts_add_t0   => P/P0(t_inf, ts + t0) }
+    #
+    def calc_thermal_fission_functions(t0, ts, read_data)
+      total_times = t0 + ts
+      ff = { ts: 0.0, ts_add_t0: 0.0 }
+
+      (0..read_data.theAn.size - 1).each do |i|
+        p_p0_tinf2ts        = (1.0 / 200.0) * read_data.theAn[i] * Math.exp(-read_data.thean[i] * ts)
+
+        p_p0_tinf2ts_add_t0 = (1.0 / 200.0) * read_data.theAn[i] * Math.exp(-read_data.thean[i] * total_times)
+
+        ff[:ts] = ff[:ts] + p_p0_tinf2ts
+        ff[:ts_add_t0] = ff[:ts_add_t0] + p_p0_tinf2ts_add_t0
+      end
+      ff
+    end
+
+    def calc_sum_thermal_fission(p_p0_source)
+      p_p0_source[:ts] - p_p0_source[:ts_add_t0]
+    end
+
+    def calc_thermal_fission_functions_with_U239(t0, ts)
+      0.00228 * 0.7 * (1 - Math.exp(-0.000491 * t0)) * Math.exp(-0.000491 * ts)
+    end
+
+    def calc_thermal_fission_functions_with_Np239(t0, ts)
+      a1 = 4.91E-4
+      a2 = 3.41E-6
+      0.00217 * 0.7 * ((a1 / (a1 + a2)) * (1 - Math.exp(-a2 * t0)) * Math.exp(-a2 * ts) - (a2 / (a1 + a2)) * (1 - Math.exp(-a1 * t0)) * Math.exp(-a1 * ts))
+    end
+  end
+
+  class RunAns1979 < RunInit
+    def run(ts = @ts, t0 = @t0, power = @power)
+      read_data = ThermalData::DataForANS_5_1_1979.new
+      ts.each_index do |i|
+        ts_to_f = ts[i].to_f
+        t0_to_f = t0[i].to_f
+        # power_to_f = power[i].to_f
+
+        t0_add_ts = ts_to_f + t0_to_f
+        f_ts2t0 = calc_thermal_fission_functions(t0_to_f, ts_to_f, read_data)
+        thePd_apostrophe = calc_sum_thermal_fission(f_ts2t0, read_data)
+        thePd = calc_total_fission_product(thePd_apostrophe, t0_to_f, ts_to_f, read_data)
+        fU239 = calc_thermal_fission_functions_with_U239(t0_to_f, ts_to_f, read_data)
+        fNp239 = calc_thermal_fission_functions_with_Np239(t0_to_f, ts_to_f, read_data)
+        thePd_all = thePd + ((fU239 + fNp239)/read_data.theQ[:others])
+        un = calc_Un(ts_to_f)
+
+        if ts_to_f > 1.0E+03 && ts_to_f < 1.0E+04
+          thePd_all_un = thePd_all * 1.03
+        elsif ts_to_f >= 1.0E+04
+          thePd_all_un = thePd_all * (1.0 + un)
+        else
+          thePd_all_un = thePd_all * (1 + 3.0E-06 * ts_to_f)
+        end
+
+        dataout(ts_to_f, thePd_all_un)
+      end
+      @output_hash
+    end
+
+    private
+
+    # Calculate thermal fission functions from ANS-5.1-1979 Table 7~9 include U235, Pu239 and U238.
+    # need t0, ts and data with thermal fission(from class DataForANS_5_1_1979)
+    # ts: Time after remove (sec)
+    # t0: Cumulative reactor operating time (sec)
+    #
+    # calc_thermal_fission_functions(t0, ts, read_data)
+    #
+    # return { :U235  => f_U235(ts, t0)
+    #          :Pu239 => f_Pu235(ts, t0)
+    #          :U238  => f_U238(ts, t0) }
+    #
+    def calc_thermal_fission_functions(t0, ts, read_data)
+      total_times = t0 + ts
+      ff = ThermalData::HashWithThermalFission.new
+      (0..read_data.theU235_alpha.size-1).each do |i|
+        f_U235_ts2tinf        = read_data.theU235_alpha[i] / read_data.theU235_lamda[i] *
+                                Math.exp(-read_data.theU235_lamda[i] * ts) *
+                                (1.0 - Math.exp(-read_data.theU235_lamda[i] * read_data.tinf))
+
+        f_U235_ts_add_t02tinf = read_data.theU235_alpha[i] / read_data.theU235_lamda[i] *
+                                Math.exp(-read_data.theU235_lamda[i] * total_times) *
+                                (1.0 - Math.exp(-read_data.theU235_lamda[i] * read_data.tinf))
+
+        ff.thermal_fission[:U235] = ff.thermal_fission[:U235] + f_U235_ts2tinf - f_U235_ts_add_t02tinf
+
+        f_Pu239_ts2tinf         = read_data.thePu239_alpha[i] / read_data.thePu239_lamda[i] *
+                                  Math.exp(-read_data.thePu239_lamda[i] * ts) *
+                                  (1.0 - Math.exp(-read_data.thePu239_lamda[i] * read_data.tinf))
+
+        f_Pu239_ts_add_t02tinf  = read_data.thePu239_alpha[i] / read_data.thePu239_lamda[i] *
+                                  Math.exp(-read_data.thePu239_lamda[i] * total_times) *
+                                  (1.0 - Math.exp(-read_data.thePu239_lamda[i] * read_data.tinf))
+
+        ff.thermal_fission[:Pu239] = ff.thermal_fission[:Pu239] + f_Pu239_ts2tinf - f_Pu239_ts_add_t02tinf
+
+        f_U238_ts2tinf          = read_data.theU238_alpha[i] / read_data.theU238_lamda[i] *
+                                  Math.exp(-read_data.theU238_lamda[i] * ts) *
+                                  (1.0 - Math.exp(-read_data.theU238_lamda[i] * read_data.tinf))
+
+        f_U235_ts_add_t02tinf   = read_data.theU238_alpha[i] / read_data.theU238_lamda[i] *
+                                  Math.exp(-read_data.theU238_lamda[i] * total_times) *
+                                  (1.0 - Math.exp(-read_data.theU238_lamda[i] * read_data.tinf))
+
+        ff.thermal_fission[:U238] = ff.thermal_fission[:U238] + f_U238_ts2tinf - f_U235_ts_add_t02tinf
+      end
+      ff.thermal_fission
+    end
+
+    # Calculate P'di (the uncorrected decay heat power) from ANS-5.1-1979 Eq.6
+    # need f_ts2t0, Pi(thePi), Qi(theQ) from class DataForANS_5_1_1979.
+    #
+    # calc_sum_thermal_fission(f_ts2t0, read_data)
+    #
+    # return P'd
+    #
+    def calc_sum_thermal_fission(f_ts2t0, read_data)
+      prd = ThermalData::HashWithThermalFission.new
+      pd = 0
+      prd.thermal_fission.each do |key, value|
+        value = value + read_data.thePi[key] * f_ts2t0[key] / read_data.theQ[key]
+        pd += value
+      end
+      pd
+    end
+
+    # Calculate total fission product decay heat power
+    # at t(ts) sec after shutdown from an operating history of T(t0) sec duration. (from ANS-5.1-1979 Eq.1 and 11)
+    # need t0, ts and data with thermal fission(from class DataForANS_5_1_1979)
+    # ts: Time after remove (sec)
+    # t0: Cumulative reactor operating time (sec)
+    #
+    # calc_total_fission_product(thePd_apostrophe, t0, ts, read_data)
+    #
+    # return Pd
+    #
+    def calc_total_fission_product(thePd_apostrophe, t0, ts, read_data)
+      g = 1.0 + ((3.24E-06 + 5.23E-10 * ts) * (t0 ** (4.0E-01)) * read_data.phi)
+      if g >= 1.1
+        g = 1.1
+      end
+      thePd_apostrophe * g
+    end
+
+    # Calculate U239 fission product decay heat power from ANS-5.1-1979 Eq.14
+    # ts: Time after remove (sec)
+    # t0: Cumulative reactor operating time (sec)
+    #
+    # calc_thermal_fission_functions_with_U239(t0, ts, read_data)
+    #
+    # return fU239
+    #
+    def calc_thermal_fission_functions_with_U239(t0, ts, read_data)
+      read_data.theEU239 * read_data.theR * (1 - Math.exp(-read_data.lamda1 * t0)) * Math.exp(-read_data.lamda1 * ts)
+    end
+
+    # Calculate Np239 fission product decay heat power from ANS-5.1-1979 Eq.15
+    # ts: Time after remove (sec)
+    # t0: Cumulative reactor operating time (sec)
+    #
+    # calc_thermal_fission_functions_with_Np239(t0, ts, read_data)
+    #
+    # return fNp239
+    #
+    def calc_thermal_fission_functions_with_Np239(t0, ts, read_data)
+      read_data.theENp239 * read_data.theR * ((read_data.lamda1 / (read_data.lamda1 - read_data.lamda2)) *
+      (1 - Math.exp(-read_data.lamda2 * t0)) * Math.exp(-read_data.lamda2 * ts) -
+      (read_data.lamda2 / (read_data.lamda1 - read_data.lamda2)) * (1 - Math.exp(-read_data.lamda1 * t0)) * Math.exp(-read_data.lamda1 * ts))
+    end
+
+    # Calculate Un is a factor which accounts for the SIL-636 additional terms effect.
+    # ts: Time after remove (sec)
+    #
+    # calc_Un(ts)
+    #
+    # return un
+    #
+    def calc_Un(ts)
+      un = (6.0E-02 - 3.0E-02) / (1.0E+06 - 1.0E+04) * (ts - 1.0E+04) + 3.0E-02
+      if un >= 6.0E-02
+        un = 6.0E-02
+      end
+      un
+    end
+  end
+
+  class RunASB9_2 < RunInit
+    def run(ts = @ts, t0 = @t0, power = @power)
+      read_data = ThermalData::DataForASB_9_2.new
+      ts.each_index do |i|
+        ts_to_f = ts[i].to_f
+        t0_to_f = t0[i].to_f
+        # power_to_f = power[i].to_f
+
+        t0_add_ts   = ts_to_f + t0_to_f
+
+        p_p0_source = calc_thermal_fission_functions(t0_to_f, ts_to_f, read_data)
+        p_p0        = calc_sum_thermal_fission(p_p0_source, ts_to_f)
+
+        p_p0_U239   = calc_thermal_fission_functions_with_U239(t0_to_f, ts_to_f)
+        p_p0_Np239  = calc_thermal_fission_functions_with_Np239(t0_to_f, ts_to_f)
+
+        p_p0_tatal = { with_k: 0, without_k: 0 }
+        p_p0_tatal[:with_k]    = p_p0[:with_k] + p_p0_U239 + p_p0_Np239
+        p_p0_tatal[:without_k] = p_p0[:without_k] + p_p0_U239 + p_p0_Np239
+
+        dataout(ts_to_f, p_p0_tatal[:with_k], p_p0_tatal[:without_k])
+      end
+      @output_hash
+    end
+
+    private
+
+    # Calculate thermal fission functions from ASB9-2.
+    # need t0, ts and data with thermal fission(from class DataForASB_9_2)
+    # ts: Time after remove (sec)
+    # t0: Cumulative reactor operating time (sec)
+    #
+    # calc_thermal_fission_functions(t0, ts, read_data)
+    #
+    # return { :P/P0(t_inf, ts)    => value
+    #          :P/P0(t_inf, ts+t0) => f_Pu235(ts, t0) }
+    #
+    def calc_thermal_fission_functions(t0, ts, read_data)
+      total_times = t0 + ts
+      ff = { ts: 0.0, ts_add_t0: 0.0 }
+
+      (0..read_data.theAn.size - 1).each do |i|
+        p_p0_tinf2ts        = (1.0 / 200.0) * read_data.theAn[i] * Math.exp(-read_data.thean[i] * ts)
+
+        p_p0_tinf2ts_add_t0 = (1.0 / 200.0) * read_data.theAn[i] * Math.exp(-read_data.thean[i] * total_times)
+
+        ff[:ts] = ff[:ts] + p_p0_tinf2ts
+        ff[:ts_add_t0] = ff[:ts_add_t0] + p_p0_tinf2ts_add_t0
+      end
+      ff
+    end
+
+    def calc_sum_thermal_fission(p_p0_source, ts)
+      p_p0 = { with_k: 0.0, without_k: 0.0 }
+      if ts >= 0 && ts < 1.0E+03
+        k1 = 0.2
+        k2 = 0.2
+      elsif ts >= 1.0E+03 && ts <= 1.0E+07
+        k1 = 0.1
+        k2 = 0.1
+      else
+        k1 = 0.1
+        k2 = 0.0
+      end
+      p_p0[:with_k] = p_p0[:with_k] + ((1 + k1) * p_p0_source[:ts]) - p_p0_source[:ts_add_t0]
+      p_p0[:without_k] = p_p0[:without_k] + ((1 + k2) * p_p0_source[:ts]) - p_p0_source[:ts_add_t0]
+      p_p0
+    end
+
+    def calc_thermal_fission_functions_with_U239(t0, ts)
+      0.00228 * 0.7 * (1 - Math.exp(-0.000491 * t0)) * Math.exp(-0.000491 * ts)
+    end
+
+    def calc_thermal_fission_functions_with_Np239(t0, ts)
+      a1 = 4.91E-4
+      a2 = 3.41E-6
+      0.00217 * 0.7 * ((a1 / (a1 + a2)) * (1 - Math.exp(-a2 * t0)) * Math.exp(-a2 * ts) - (a2 / (a1 + a2)) * (1 - Math.exp(-a1 * t0)) * Math.exp(-a1 * ts))
+    end
+  end
+end

data/lib/decay_heat_with_nuclear/thermal_data.rb

@@ -0,0 +1,328 @@
+module ThermalData
+  class HashWithThermalFission
+    attr_accessor :thermal_fission
+
+    def initialize(init_data = 0)
+      @thermal_fission = {}
+      ary = [:U235, :Pu239, :U238]
+      ary.each { |a| @thermal_fission[a] = init_data }
+    end
+  end
+
+  # ANSI/ANS-5.1-1979, American National Standard for Decay
+  # Heat Power in Light Water Reactors.
+  class DataForANS_5_1_1979
+    # phi is the number of fissions per initial fissile atoms.
+    # tinf is t = inf, setting it 1.0E+13.
+    # R(theR) is the nuber of U-239 atoms produced per fission.
+    # lamda1 is decay constant for U239 and Np239.
+    # lamda2 is decay constant for Np239.
+    # theEU239 is average energy from decay of on U239 atom (0.474 MeV)
+    # theENp239 is average energy from decay of on Np239 atom (0.419 MeV)
+    attr_reader :phi, :tinf, :theR, :lamda1, :lamda2, :theEU239, :theENp239
+
+    def initialize
+      @phi = 3.0
+      @tinf = 1.0E+13
+      @theR = 6.20E-01
+      @lamda1 = 4.91E-04
+      @lamda2 = 3.41E-06
+      @theEU239 = 0.474
+      @theENp239 = 0.419
+    end
+
+    # thePi is Pi in ANS-5.1-1979.
+    def thePi
+      hash = HashWithThermalFission.new
+
+      hash.thermal_fission[:U235]  = 0.98
+      hash.thermal_fission[:Pu239] = 0.01
+      hash.thermal_fission[:U238]  = 0.01
+
+      hash.thermal_fission
+    end
+
+    # theQ is Q in ANS-5.1-1979.
+    # Total recoverable energy assicuated with on fission of nuclide.
+    # Default is 200 MeV/sec in U235, Pu239, U238 and others
+    def theQ
+      hash = HashWithThermalFission.new(200.0)
+      hash.thermal_fission[:others] = 200.0
+      hash.thermal_fission
+    end
+
+    # theU235_alpha is alpha in ANS-5.1-1979 Table 7.
+    def theU235_alpha
+      array = Array.new(23)
+
+      array[0] = 6.5057E-01
+      array[1] = 5.1264E-01
+      array[2] = 2.4384E-01
+      array[3] = 1.3850E-01
+      array[4] = 5.544E-02
+      array[5] = 2.2225E-02
+      array[6] = 3.3088E-03
+      array[7] = 9.3015E-04
+      array[8] = 8.0943E-04
+      array[9] = 1.9567E-04
+      array[10] = 3.2535E-05
+      array[11] = 7.5595E-06
+      array[12] = 2.5232E-06
+      array[13] = 4.9948E-07
+      array[14] = 1.8531E-07
+      array[15] = 2.6608E-08
+      array[16] = 2.2398E-09
+      array[17] = 8.1641E-12
+      array[18] = 8.7797E-11
+      array[19] = 2.5131E-14
+      array[20] = 3.2176E-16
+      array[21] = 4.5038E-17
+      array[22] = 7.4791E-17
+
+      array
+    end
+
+    # theU235_lamda is lamda in ANS-5.1-1979 Table 7.
+    def theU235_lamda
+      array = Array.new(23)
+
+      array[0] = 2.2138E+01
+      array[1] = 5.1587E-01
+      array[2] = 1.9594E-01
+      array[3] = 1.0314E-01
+      array[4] = 3.3656E-02
+      array[5] = 1.1681E-02
+      array[6] = 3.5870E-03
+      array[7] = 1.3930E-03
+      array[8] = 6.2630E-04
+      array[9] = 1.8906E-04
+      array[10] = 5.4988E-05
+      array[11] = 2.0958E-05
+      array[12] = 1.0010E-05
+      array[13] = 2.5438E-06
+      array[14] = 6.6361E-07
+      array[15] = 1.2290E-07
+      array[16] = 2.7213E-08
+      array[17] = 4.3714E-09
+      array[18] = 7.5780E-10
+      array[19] = 2.4786E-10
+      array[20] = 2.2384E-13
+      array[21] = 2.4600E-14
+      array[22] = 1.5699E-14
+
+      array
+    end
+
+    # thePu239_alpha is alpha in ANS-5.1-1979 Table 8.
+    def thePu239_alpha
+      array = Array.new(23)
+
+      array[0] = 2.083E-01
+      array[1] = 3.853E-01
+      array[2] = 2.213E-01
+      array[3] = 9.460E-02
+      array[4] = 3.531E-02
+      array[5] = 2.292E-02
+      array[6] = 3.946E-03
+      array[7] = 1.317E-03
+      array[8] = 7.052E-04
+      array[9] = 1.432E-04
+      array[10] = 1.765E-05
+      array[11] = 7.347E-06
+      array[12] = 1.747E-06
+      array[13] = 5.481E-07
+      array[14] = 1.671E-07
+      array[15] = 2.112E-08
+      array[16] = 2.996E-09
+      array[17] = 5.107E-11
+      array[18] = 5.730E-11
+      array[19] = 4.138E-14
+      array[20] = 1.088E-15
+      array[21] = 2.454E-17
+      array[22] = 7.557E-17
+
+      array
+    end
+
+    # thePu239_lamda is lamda in ANS-5.1-1979 Table 8.
+    def thePu239_lamda
+      array = Array.new(23)
+
+      array[0] = 1.002E+01
+      array[1] = 6.433E-01
+      array[2] = 2.186E-01
+      array[3] = 1.004E-01
+      array[4] = 3.728E-02
+      array[5] = 1.435E-02
+      array[6] = 4.549E-03
+      array[7] = 1.328E-03
+      array[8] = 5.356E-04
+      array[9] = 1.730E-04
+      array[10] = 4.881E-05
+      array[11] = 2.006E-05
+      array[12] = 8.319E-06
+      array[13] = 2.358E-06
+      array[14] = 6.450E-07
+      array[15] = 1.278E-07
+      array[16] = 2.466E-08
+      array[17] = 9.378E-09
+      array[18] = 7.450E-10
+      array[19] = 2.426E-10
+      array[20] = 2.210E-13
+      array[21] = 2.640E-14
+      array[22] = 1.380E-14
+
+      array
+    end
+
+    # theU238_alpha is alpha in ANS-5.1-1979 Table 9.
+    def theU238_alpha
+      array = Array.new(23)
+
+      array[0] = 1.2311E+0
+      array[1] = 1.1486E+0
+      array[2] = 7.0701E-01
+      array[3] = 2.5209E-01
+      array[4] = 7.187E-02
+      array[5] = 2.8291E-02
+      array[6] = 6.8382E-03
+      array[7] = 1.2322E-03
+      array[8] = 6.8409E-04
+      array[9] = 1.6975E-04
+      array[10] = 2.4182E-05
+      array[11] = 6.6356E-06
+      array[12] = 1.0075E-06
+      array[13] = 4.9894E-07
+      array[14] = 1.6352E-07
+      array[15] = 2.3355E-08
+      array[16] = 2.8094E-09
+      array[17] = 3.6236E-11
+      array[18] = 6.4577E-11
+      array[19] = 4.4963E-14
+      array[20] = 3.6654E-16
+      array[21] = 5.6293E-17
+      array[22] = 7.1602E-17
+
+      array
+    end
+
+    # theU238_lamda is lamda in ANS-5.1-1979 Table 9.
+    def theU238_lamda
+      array = Array.new(23)
+
+      array[0] = 3.2881E+0
+      array[1] = 9.3805E-01
+      array[2] = 3.7073E-01
+      array[3] = 1.1118E-01
+      array[4] = 3.6143E-02
+      array[5] = 1.3272E-02
+      array[6] = 5.0133E-03
+      array[7] = 1.3655E-03
+      array[8] = 5.5158E-04
+      array[9] = 1.7873E-04
+      array[10] = 4.9032E-05
+      array[11] = 1.7058E-05
+      array[12] = 7.0465E-06
+      array[13] = 2.3190E-06
+      array[14] = 6.4480E-07
+      array[15] = 1.2649E-07
+      array[16] = 2.5548E-08
+      array[17] = 8.4782E-09
+      array[18] = 7.5130E-10
+      array[19] = 2.4188E-10
+      array[20] = 2.2739E-13
+      array[21] = 9.0536E-14
+      array[22] = 5.6098E-15
+
+      array
+    end
+  end
+
+  # ANS/ANSI-5.1-1973, Standard for Decay Heat.
+  class DataForANS_5_1_1973
+    def initialize
+    end
+
+    # theAn is An in ANS/ANSI-5.1-1973.
+    def theAn
+      array = Array.new(11)
+
+      array[0] = 0.598
+      array[1] = 1.65
+      array[2] = 3.1
+      array[3] = 3.87
+      array[4] = 2.33
+      array[5] = 1.29
+      array[6] = 0.462
+      array[7] = 0.328
+      array[8] = 0.17
+      array[9] = 0.0865
+      array[10] = 0.114
+
+      array
+    end
+
+    # thean is an in ANS/ANSI-5.1-1973.
+    def thean
+      array = Array.new(10)
+
+      array[0] = 1.772E+00
+      array[1] = 5.774E-01
+      array[2] = 6.743E-02
+      array[3] = 6.214E-03
+      array[4] = 4.739E-04
+      array[5] = 4.810E-05
+      array[6] = 5.344E-06
+      array[7] = 5.716E-07
+      array[8] = 1.036E-07
+      array[9] = 2.959E-08
+      array[10] = 7.585E-10
+
+      array
+    end
+  end
+
+  # ASB9-2, Standard for Decay Heat.
+  class DataForASB_9_2
+    def initialize
+    end
+
+    # theAn is An in ASB9-2.
+    def theAn
+      array = Array.new(11)
+
+      array[0] = 0.598
+      array[1] = 1.65
+      array[2] = 3.1
+      array[3] = 3.87
+      array[4] = 2.33
+      array[5] = 1.29
+      array[6] = 0.462
+      array[7] = 0.328
+      array[8] = 0.17
+      array[9] = 0.0865
+      array[10] = 0.114
+
+      array
+    end
+
+    # thean is an in ASB9-2.
+    def thean
+      array = Array.new(10)
+
+      array[0] = 1.772E+00
+      array[1] = 5.774E-01
+      array[2] = 6.743E-02
+      array[3] = 6.214E-03
+      array[4] = 4.739E-04
+      array[5] = 4.810E-05
+      array[6] = 5.344E-06
+      array[7] = 5.716E-07
+      array[8] = 1.036E-07
+      array[9] = 2.959E-08
+      array[10] = 7.585E-10
+
+      array
+    end
+  end
+end

data/lib/decay_heat_with_nuclear/tools.rb

@@ -0,0 +1,21 @@
+module Tools
+  class TextToArray
+    attr_accessor :output
+
+    def initialize(filename)
+      @output = {}
+      ts = []
+      t0 = []
+      power = []
+      f = File.open("#{filename}", 'r')
+      f.each_line do |line|
+        a, b, c = line.split("\t")
+        ts << a.to_f
+        t0 << b.to_f
+        power << c.to_f
+      end
+      f.close
+      @output = { ts: ts, t0: t0, power: power }
+    end
+  end
+end

data/lib/decay_heat_with_nuclear/version.rb

@@ -0,0 +1,3 @@
+module DecayHeatWithNuclear
+  VERSION = '0.0.2'
+end

data/lib/decay_heat_with_nuclear.rb

@@ -0,0 +1,18 @@
+require 'decay_heat_with_nuclear/version'
+require 'decay_heat_with_nuclear/thermal_data'
+require 'decay_heat_with_nuclear/tools'
+require 'decay_heat_with_nuclear/main_run'
+
+module DecayHeatWithNuclear
+  def self.run(hash)
+    test1 = MainRun::RunAns1979.new(hash)
+    test2 = MainRun::RunASB9_2.new(hash)
+    test3 = MainRun::RunAns1973.new(hash)
+
+    {
+      ans1979: test1.run,
+      ans1973: test3.run,
+      asb9_2:  test2.run
+    }
+  end
+end

data/lib/run_this.rb

@@ -0,0 +1,69 @@
+require './text_to_array.rb'
+require './run_ans_1979.rb'
+require './run_ans_1973.rb'
+require './run_asb_9_2.rb'
+
+# 執行此程式需要一個hash
+# 裡面必須含有ts (sec)、t0 (sec)、power (MW)
+#
+# hash =  {ts: array,
+#          t0: array,
+#          power: array}
+#
+# 抑或可以利用TextToArray.new(相對路徑)讀取文字檔
+# 文字檔格式如下，其中(TAB)表示輸入Tab鍵、(ENTER)表示輸入Enter鍵
+#
+# ts(TAB)t0(TAB)power(ENTER)
+#
+# ks_ts_t0_power.txt為範例檔
+
+output = TextToArray.new('./ks_ts_t0_power.txt')
+hash1 = output.output
+
+# ts = []
+# times = 1
+# t_index = 0
+# while t_index <= 1 * 365 * 86400
+#   index = 1
+#   while index < 10
+#     t_index = index * times
+#     ts << t_index
+#     index = index + 1
+#   end
+#   times = times * 10
+# end
+# ts << 365 * 86400
+ts = Array.new(20) { |i| (i + 1) * 365 * 24 * 3600 }
+hash2 = { ts: ts,
+          t0: Array.new(ts.size) { 63 * 30 * 24 * 3600 },
+          power: Array.new(ts.size) { 2943 } }
+
+# 執行RunAns1979.new(hash, case)來配置前置檔案
+# hash即為上面說明的hash
+# case為輸出的類型
+#
+# case = 0
+#  ts(year)     t(sec)           P/P0
+#     1        31536000      0.000281683414
+#
+# case = 1
+#       ts             P/P0
+# ts =    1 years P/P0= 0.000281683414
+#
+# case = 2
+#         ts                     t0                 un                P/P0 without un/K            P/P0 with un/K           power
+# ts = 31536000.0 sec, t0 = 141912000.0 sec, un = 0.06000000, P/P0(without un) = 0.00028168, P/P0(with un) = 0.00029858, power = 0.90471 MW
+# P/P0 with un = (P/P0 without un) * (1 + un)
+#
+# P.S RunASB9_2的K值比較特殊，without K 的意思是指當 ts > 10^7 時將 K = 0
+#
+# 或是可以編輯'run_ans_1979.rb'第37行，增加自訂的輸出方式
+
+test2 = RunASB9_2.new(hash2, 3, 'ASB9-2_1to20years.txt')
+test2.run
+
+test3 = RunAns1973.new(hash2, 3, 'Ans1973_1to20years.txt')
+test3.run
+
+test1 = RunAns1979.new(hash2, 3, 'Ans1979_1to20years.txt')
+test1.run

metadata

@@ -0,0 +1,86 @@
+--- !ruby/object:Gem::Specification
+name: decay_heat_with_nuclear
+version: !ruby/object:Gem::Version
+  version: 0.0.2
+platform: ruby
+authors:
+- C.lin
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2015-07-09 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.6'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.6'
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '10.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '10.0'
+description: It's a easy way to get nuclear fuel decay heat when it remove from reactor
+  core.
+email:
+- demeter.yeh@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- Gemfile
+- LICENSE.txt
+- README.md
+- Rakefile
+- decay_heat_with_nuclear.gemspec
+- lib/decay_heat_with_nuclear.rb
+- lib/decay_heat_with_nuclear/main_run.rb
+- lib/decay_heat_with_nuclear/thermal_data.rb
+- lib/decay_heat_with_nuclear/tools.rb
+- lib/decay_heat_with_nuclear/version.rb
+- lib/run_this.rb
+homepage: ''
+licenses:
+- GPL3.0
+metadata: {}
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 2.2.2
+signing_key: 
+specification_version: 4
+summary: It's a easy way to get nuclear fuel decay heat when it remove from reactor
+  core.
+test_files: []