gtengine 0.0.2 → 0.0.3

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: 63698ae2d1f94f7ac92937dc0d0219aab3eaa77b
-  data.tar.gz: b26b469da2e306d1d8563b9744f531e9ed513a52
+  metadata.gz: dfce175060b2ea4ef910bfb47e3c33c55668945a
+  data.tar.gz: 084fd7277b4bc64a8621329e09d7ce543dfa614e
 SHA512:
-  metadata.gz: 1d2c9805c67d1de50113b1154b33966d4d054278ebd590e69a0ee02f7d33caf965e5fbae8b0bb063d95366899ba9d05e33cd67815bc062ebe7319737e71bd2de
-  data.tar.gz: 30d3b7e15818c2f0432fd6e029cd296ecc8ca1228a24a38ad5e961241d8114223e7fa178a002720a8584df9ac915231b4e00719222e8ce3caaaad4abc8f8effb
+  metadata.gz: 071a5a798f8547c11cd9eeec5e5b6157c48771cb4bc7d25b17da6828d397486e534b1bb8e7c18b3af608dec4bd4ac7ffbe9a3f350f1b0a9437d64a988823ae61
+  data.tar.gz: a542b14b7845f268c9846f393b7256955573d4860715330553f4b8d00d2ff61639673bc3dfdf01a6410d3da71e098df6a3dad13681168d3a0ec43b9ee3385549

data/Readme.md

@@ -6,11 +6,18 @@ Math models of gas turbine engine and its parts made with ruby language
 
 ## Usage
 
+Add gem to your gemfile
+
+    gem 'gtengine'
+
+and run
+
+    bundle
+
 ### Simple
 
 For simple parts calculation:
 
-    require_relative 'lib/gtengine'
     gas = Gtengine::Gas.new 300, 101325
 
     k = Gtengine::Simple::Compressor.new gas, 4.5, 0.85

data/lib/gtengine/gas.rb

@@ -1,12 +1,9 @@
 class Gtengine::Gas
-
   R = 287.4
   attr_accessor :t, :p, :alfa
 
-  def initialize t = 280.0, p = 101325.0, alfa = 99999999.0
-    @t = t.to_f
-    @p = p.to_f
-    @alfa = alfa
+  def initialize(t = 300.0, p = 101325.0, alfa = 99999999.0)
+    @t, @p, @alfa = t.to_f, p.to_f, alfa.to_f
   end
 
   def density
@@ -14,11 +11,11 @@ class Gtengine::Gas
   end
 
   def k
-    k_from_cp cp
+    k_from(cp)
   end
 
   def average_k
-    k_from_cp average_cp
+    k_from(average_cp)
   end
 
   def cp
@@ -33,23 +30,47 @@ class Gtengine::Gas
 
   private
 
-    def k_from_cp c_p
-      c_p / (c_p - R)
-    end
+  def k_from(c_p)
+    c_p / (c_p - R)
+  end
+
+  def high_temperature_cp
+    (first_high_cp - second_high_cp * third_high_cp) * 4187.0
+  end
+
+  def first_high_cp
+    0.0267 / alfa + 0.26 + (0.032 + 0.0133 / alfa) * (0.001 * 1.176 * t - 0.88235)
+  end
+
+  def second_high_cp
+    0.374 * 0.01 + 0.0094 / (alfa ** 2.0 + 10.0)
+  end
+
+  def third_high_cp
+    5.5556 * 0.000001 * t ** 2.0 - 1.3056 * 0.01 * t + 6.67
+  end
 
-    def high_temperature_cp
-      (0.0267 / alfa + 0.26 + (0.032 + 0.0133 / alfa) * (0.001 * 1.176 * t - 0.88235) - (0.374 * 0.01 + 0.0094 / (alfa ** 2.0 + 10.0)) * (5.5556 * 0.000001 * t ** 2.0 - 1.3056 * 0.01 * t + 6.67)) * 4187.0
-    end
+  def low_temperature_cp
+    (first_low_cp + second_low_cp * third_low_cp) * 4187.0
+  end
 
-    def low_temperature_cp
-      (0.0174 / alfa + 0.2407 + (0.0193 + 0.0093 / alfa) * (0.001 * 2.5 * t - 0.875) + (0.002 - 0.001 * 1.056 / (alfa - 0.2)) * (2.5 * 0.00001 * t ** 2.0 - 0.0275 * t + 6.5625)) * 4187.0
-    end    
+  def first_low_cp
+    0.0174 / alfa + 0.2407 + (0.0193 + 0.0093 / alfa) * (0.001 * 2.5 * t - 0.875)
+  end
 
-    def high_temperature_average_cp
-      (((1.25 + 2.2 * alfa) / (alfa * 100000.0)) * (t + 450.0) + 0.218) * 4187.0
-    end
+  def second_low_cp
+    0.002 - 0.001 * 1.056 / (alfa - 0.2)
+  end
 
-    def low_temperature_average_cp
-      (((2.25 + 1.2 * alfa) / (alfa * 100000.0)) * (t - 70.0) + 0.236) * 4187.0
-    end
+  def third_low_cp
+    2.5 * 0.00001 * t ** 2.0 - 0.0275 * t + 6.5625
+  end
+
+  def high_temperature_average_cp
+    (((1.25 + 2.2 * alfa) / (alfa * 100000.0)) * (t + 450.0) + 0.218) * 4187.0
+  end
+
+  def low_temperature_average_cp
+    (((2.25 + 1.2 * alfa) / (alfa * 100000.0)) * (t - 70.0) + 0.236) * 4187.0
+  end
 end

data/lib/gtengine/simple/burner.rb

@@ -1,44 +1,46 @@
 class Gtengine::Simple::Burner
+  attr_accessor :input, :q_t, :output, :t_g, :options
 
-  attr_accessor :input, :q_t, :output, :t_g
+  DEFAULTS = {
+    t_0: 288.3,
+    q_n: 43000000.0,
+    eta_g: 0.985,
+    l_0: 14.7
+  }
 
-  T0 = 288.3
-  QN = 43000000.0
-  ETA_G = 0.985
-  L0 = 14.7
-
-  def initialize input, t_g
+  def initialize(input, t_g, options = {})
     @input, @t_g = input, t_g
+    @options = DEFAULTS.merge(options)
     cycle
   end
 
   def cp_vh
-    @input.cp
+    input.cp
   end
 
   def cp_vyh
-    @output.cp
+    output.cp
   end
 
   def t_vh
-    @input.t
+    input.t
   end
 
   def p_vh
-    @input.p
+    input.p
   end
 
   def q_ks
-    (cp_vyh * t_g - cp_vh * t_vh - (cp_vyh - cp_vh) * T0) / (QN * ETA_G - (cp_vyh * t_g - cp_mult_t_0))
+    upper_q_ks / lower_q_ks
   end
 
   def cp_mult_t_0
-    Gtengine::Gas.new(T0, 101325.0, 1.0).cp * T0
+    Gtengine::Gas.new(t_0, 101325.0, 1.0).cp * t_0
   end
 
   def alfa
     begin
-      1.0 / (q_ks * L0)
+      1.0 / (q_ks * l_0)
     rescue
       999999999.0
     end
@@ -49,9 +51,35 @@ class Gtengine::Simple::Burner
     5.times { @output.alfa = alfa }
   end
 
+  def t_0
+    options[:t_0]
+  end
+
+  def l_0
+    options[:l_0]
+  end
+
+  def eta_g
+    options[:eta_g]
+  end
+
+  def q_n
+    options[:q_n]
+  end
+
   def info
     puts "== Burner q_ks: #{q_ks}, ALFA: #{alfa}"
     puts "==== Вход T: #{@input.t.to_i} K, P: #{@input.p.to_i} Па"
     puts "==== Выход T: #{@output.t.to_i} K, P: #{@output.p.to_i} Па\n\n"
   end
+
+  private
+
+  def upper_q_ks
+    (cp_vyh * t_g - cp_vh * t_vh - (cp_vyh - cp_vh) * t_0)
+  end
+
+  def lower_q_ks
+    (q_n * eta_g - (cp_vyh * t_g - cp_mult_t_0))
+  end
 end

data/lib/gtengine/simple/compressor.rb

@@ -1,58 +1,66 @@
 class Gtengine::Simple::Compressor
-  attr_accessor :input, :output, :g, :pi_k, :kpd, :average
+  attr_accessor :input, :output, :g, :pi_k, :average, :options
 
-  def initialize input, pi_k, kpd
-    @input, @pi_k, @kpd = input, pi_k.to_f, kpd.to_f
+  DEFAULTS = {
+    kpd: 0.85
+  }
+
+  def initialize(input, pi_k, options = {})
+    @input, @pi_k = input, pi_k.to_f
+    @options = DEFAULTS.merge(options)
     cycle
   end
 
+  def kpd
+    options[:kpd]
+  end
+
   def t_vyh
-    @input.t * (1.0 + ((pi_k ** ((k - 1.0) / k) - 1.0) / kpd))
+    input.t * (1.0 + ((pi_k ** ((k - 1.0) / k) - 1.0) / kpd))
   end
 
   def p_vyh
-    @input.p * @pi_k
+    input.p * pi_k
   end
 
   def l_k
-    @average.cp * t_vyh * (@pi_k ** ((k - 1.0) / k) - 1.0) / @kpd
+    average.cp * t_vyh * (pi_k ** ((k - 1.0) / k) - 1.0) / kpd
   end
 
   def k
-    @average.average_k
+    average.average_k
   end
 
   def cp
-    @average.average_cp
+    average.average_cp
   end
 
   def alfa
     1.577 * 0.0000001 * t_vyh ** 2.383 + 1.774
   end
 
-  def update_average
-    @average.t = (t_vh + t_vyh) / 2.0
-    @average.p = (p_vh + p_vyh) / 2.0
-  end
-
   def t_vh
-    @input.t
+    input.t
   end
 
   def p_vh
-    @input.p
+    input.p
+  end
+
+  def update_average
+    @average.t = (t_vh + t_vyh) / 2.0
+    @average.p = (p_vh + p_vyh) / 2.0
   end
 
   def cycle
-    @average = Gtengine::Gas.new input.t.to_f, input.p.to_f
+    @average = Gtengine::Gas.new(input.t.to_f, input.p.to_f)
     5.times { update_average }
-    @output = Gtengine::Gas.new t_vyh, p_vyh
+    @output = Gtengine::Gas.new(t_vyh, p_vyh)
   end
 
   def info
-    puts "== Compressor Pi_k: #{@pi_k}, КДП: #{@kpd}, Cp_sr: #{cp}, K: #{k}"
-    puts "==== Вход T: #{@input.t.to_i} K, P: #{@input.p.to_i} Па"
-    puts "==== Выход T: #{@output.t.to_i} K, P: #{@output.p.to_i} Па, ALFA: #{alfa}\n\n"
+    puts "== Compressor Pi_k: #{@pi_k}, КДП: #{kpd}, Cp_sr: #{cp}, K: #{k}"
+    puts "==== Вход T: #{input.t.to_i} K, P: #{input.p.to_i} Па"
+    puts "==== Выход T: #{output.t.to_i} K, P: #{output.p.to_i} Па, ALFA: #{alfa}\n\n"
   end
-
-end
+end

data/lib/gtengine/simple/pik_optimizer.rb

@@ -1,32 +1,31 @@
 class Gtengine::Simple::PikOptimizer
   attr_accessor :start_pik, :end_pik, :t_g, :step, :optimum, :cycles, :air
 
-  def initialize air, start_pik, end_pik, t_g, step=1
+  def initialize(air, start_pik, end_pik, t_g, step = 1)
     @air, @start_pik, @end_pik, @t_g, @step = air, start_pik, end_pik, t_g, step
+    @cycles = []
     optimize
   end
 
   def optimize
-    @cycles = []
     Range.new(start_pik, end_pik).step(step).each do |pik|
       calc_result(pik)
     end
   end
 
   def calc_result(pik)
-    begin
-      c = Gtengine::Simple::Cycle.new(air, pik, t_g) 
-      @cycles << c if c.burner.q_ks
-    rescue
-    end
+    c = Gtengine::Simple::Cycle.new(air, pik, t_g)
+    @cycles << c if c.burner.q_ks
+  rescue
+    false
   end
 
   def optimal
-    @cycles.sort {|a, b| a.burner.q_ks <=> b.burner.q_ks}.first
+    cycles.sort { |a, b| a.burner.q_ks <=> b.burner.q_ks }.first
   end
 
   def info
-    @cycles.each { |c| c.info } 
+    cycles.each { |c| c.info }
     puts "Optimal pi_k: #{optimal.pi_k}, q_ks: #{optimal.q_ks}"
   end
 end

data/lib/gtengine/simple/turbine.rb

@@ -1,18 +1,37 @@
 class Gtengine::Simple::Turbine
+  attr_accessor :burner, :l_k, :average, :output, :options
+
+  DEFAULTS = {
+    eta: 0.9,
+    eta_m: 0.985,
+    g_ohl: 0.001,
+    kpd: 0.9
+  }
+
+  def initialize(burner, l_k, options = {})
+    @burner, @l_k = burner, l_k
+    @options = DEFAULTS.merge(options)
+    cycle
+  end
 
-  attr_accessor :burner, :l_k, :kpd, :average, :output
+  def input
+    burner.output
+  end
 
-  ETA = 0.9
-  ETA_M = 0.985
-  G_OHL = 0.001
+  def kpd
+    options[:kpd]
+  end
 
-  def initialize burner, l_k, kpd = 0.9
-    @burner, @l_k, @kpd = burner, l_k, kpd.to_f
-    cycle
+  def g_ohl
+    options[:g_ohl]
   end
 
-  def input
-    @burner.output
+  def eta_m
+    options[:eta_m]
+  end
+
+  def eta
+    options[:eta]
   end
 
   def t_vh
@@ -24,16 +43,13 @@ class Gtengine::Simple::Turbine
   end
 
   def pi_t
-    pit = l_k / (1.0 + burner.q_ks)
-    pit /= (t_vh * cp * ETA * ETA_M * (1.0 - G_OHL))
-    pit = (1.0 / (1.0 - pit)) ** ((k - 1.0) / k)
-    pit
+    (1.0 / (1.0 - temp_pi_t)) ** ((k - 1.0) / k)
   end
 
   def cycle
-    @average = Gtengine::Gas.new t_vh, p_vh
+    @average = Gtengine::Gas.new(t_vh, p_vh)
     5.times { update_average }
-    @output = Gtengine::Gas.new t_vyh, p_vyh
+    @output = Gtengine::Gas.new(t_vyh, p_vyh)
   end
 
   def update_average
@@ -54,17 +70,22 @@ class Gtengine::Simple::Turbine
   end
 
   def k
-    @average.average_k
+    average.average_k
   end
 
   def cp
-    @average.average_cp
+    average.average_cp
   end
 
   def info
     puts "== Turbine Cp_sr: #{cp}, K: #{k}, Pi_t: #{pi_t}, КДП: #{ETA}"
     puts "==== Вход T: #{t_vh.to_i} K, P: #{p_vh.to_i} Па, ALFA: #{input.alfa}, Cp: #{input.cp}"
-    puts "==== Выход T: #{@output.t.to_i} K, P: #{@output.p.to_i} Па, , Cp: #{@output.cp}\n\n"
+    puts "==== Выход T: #{output.t.to_i} K, P: #{output.p.to_i} Па, , Cp: #{output.cp}\n\n"
   end
 
-end
+  private
+
+  def temp_pi_t
+    (l_k / (1.0 + burner.q_ks)) / (t_vh * cp * eta * eta_m * (1.0 - g_ohl))
+  end
+end

data/lib/gtengine/simple_cycle.rb

@@ -8,25 +8,23 @@ module Gtengine
 
       attr_accessor :air, :pi_k, :t_g, :compressor, :turbine, :burner
 
-      def initialize air=Gas.new(300, 101325), pi_k, t_g
-        @air = air
-        @t_g = t_g
-        @pi_k = pi_k
+      def initialize(air = Gas.new, pi_k, t_g)
+        @air, @t_g, @pi_k = air, t_g, pi_k
         init_cycle
       end
 
       def init_cycle
-        @compressor = Compressor.new air, pi_k, 0.85
-        @burner = Burner.new @compressor.output, t_g
-        @turbine = Turbine.new burner, compressor.l_k, 0.9
+        @compressor = Compressor.new(air, pi_k)
+        @burner = Burner.new(@compressor.output, t_g)
+        @turbine = Turbine.new(burner, compressor.l_k)
       end
 
       def q_ks
-        @burner.q_ks
+        burner.q_ks
       end
 
       def info
-        puts "== Cycle info: pi_k: #{@compressor.pi_k}, q_ks: #{q_ks}"
+        puts "== Cycle info: pi_k: #{compressor.pi_k}, q_ks: #{q_ks}"
       end
     end
   end

data/lib/gtengine/version.rb

@@ -1,3 +1,3 @@
 module Gtengine
-  VERSION = "0.0.2"
+  VERSION = "0.0.3"
 end

data/spec/models/{simple/gas_spec.rb → gas_spec.rb}

@@ -4,24 +4,24 @@ describe Gtengine::Gas do
   let(:gas) { Gtengine::Gas.new }
 
   it "Считает плотность" do
-    expect(gas.density).to eq 1.259133611691023
+    expect(gas.density).to eq 1.1751913709116215
   end
 
   it "Считает k" do
-    expect(gas.k).to eq 1.4029968590096924
+    expect(gas.k).to eq 1.4019467210274585
   end
 
   it "Считает k за процесс" do
-    expect(gas.average_k).to eq 1.404058445015607
+    expect(gas.average_k).to eq 1.4034884084723653
   end
 
   context "LOW TEMPERATURE" do
     it "Считает теплоемкость" do
-      expect(gas.cp).to eq 1000.556923124028
+      expect(gas.cp).to eq 1002.4201381549933
     end
 
     it "Считает теплоемкость за процесс" do
-      expect(gas.average_cp).to eq 998.6832401978357
+      expect(gas.average_cp).to eq 999.6881202166772
     end
   end
 

data/spec/models/simple/burner_spec.rb

@@ -2,14 +2,14 @@ require 'spec_helper'
 
 describe Gtengine::Simple::Burner do
   let(:gas) { Gtengine::Gas.new }
-  let(:compressor) { Gtengine::Simple::Compressor.new gas, 10, 0.87 }
-  let(:burner) { Gtengine::Simple::Burner.new compressor.output, 1400 }
+  let(:compressor) { Gtengine::Simple::Compressor.new(gas, 10) }
+  let(:burner) { Gtengine::Simple::Burner.new(compressor.output, 1400) }
 
   it "считает расход топлива" do
-    expect(burner.q_ks).to eq 0.026948712923834613
+    expect(burner.q_ks).to eq 0.0255346596182821
   end
 
   it "считает альфа" do
-    expect(burner.alfa).to eq 2.524321331286569
+    expect(burner.alfa).to eq 2.664112696283923
   end
 end

data/spec/models/simple/compressor_spec.rb

@@ -2,21 +2,21 @@ require 'spec_helper'
 
 describe Gtengine::Simple::Compressor do
   let(:gas) { Gtengine::Gas.new }
-  let(:compressor) { Gtengine::Simple::Compressor.new gas, 10, 0.87 }
+  let(:compressor) { Gtengine::Simple::Compressor.new(gas, 10) }
 
   it "считает давление выхода" do
     expect(compressor.p_vyh).to eq 1013250.0
   end
 
   it "считает температуру выхода" do
-    expect(compressor.t_vyh).to eq 579.410983909093
+    expect(compressor.t_vyh).to eq 627.5847203290919
   end
 
   it "считает работу компрессора" do
-    expect(compressor.l_k).to eq 631083.3734847266
+    expect(compressor.l_k).to eq 701724.7666737017
   end
 
   it "считает альфа" do
-    expect(compressor.alfa).to eq 2.37938137624045
+    expect(compressor.alfa).to eq 2.5062930805003583
   end
 end

data/spec/models/simple/turbine_spec.rb

@@ -2,24 +2,24 @@ require 'spec_helper'
 
 describe Gtengine::Simple::Turbine do
   let(:gas) { Gtengine::Gas.new }
-  let(:compressor) { Gtengine::Simple::Compressor.new gas, 10, 0.87 }
-  let(:burner) { Gtengine::Simple::Burner.new compressor.output, 900 }
-  let(:turbine) { Gtengine::Simple::Turbine.new burner, compressor.l_k }
+  let(:compressor) { Gtengine::Simple::Compressor.new(gas, 10) }
+  let(:burner) { Gtengine::Simple::Burner.new(compressor.output, 900) }
+  let(:turbine) { Gtengine::Simple::Turbine.new(burner, compressor.l_k) }
 
   it "считает степень понижения давления" do
-    expect(turbine.pi_t).to eq 1.4858595860294537
+    expect(turbine.pi_t).to eq 1.6842480477100261
   end
 
   it "считает температуру выхода" do
-    expect(turbine.t_vyh).to eq 815.5179703599931
+    expect(turbine.t_vyh).to eq 790.5487885221904
   end
 
   it "считает давление выхода" do
-    expect(turbine.p_vyh).to eq 681928.5008670494
+    expect(turbine.p_vyh).to eq 601603.7847736602
   end
 
   it "считает работу турбины l_t" do
-    expect(turbine.l_t).to eq 87289.27555567237
+    expect(turbine.l_t).to eq 112962.28663065929
   end
-  
+
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: gtengine
 version: !ruby/object:Gem::Version
-  version: 0.0.2
+  version: 0.0.3
 platform: ruby
 authors:
 - Oleg Bovykin
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2014-01-12 00:00:00.000000000 Z
+date: 2014-09-04 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler
@@ -91,9 +91,9 @@ files:
 - lib/gtengine/simple/turbine.rb
 - lib/gtengine/simple_cycle.rb
 - lib/gtengine/version.rb
+- spec/models/gas_spec.rb
 - spec/models/simple/burner_spec.rb
 - spec/models/simple/compressor_spec.rb
-- spec/models/simple/gas_spec.rb
 - spec/models/simple/turbine_spec.rb
 - spec/spec_helper.rb
 homepage: https://github.com/arrowcircle/engine
@@ -116,13 +116,13 @@ required_rubygems_version: !ruby/object:Gem::Requirement
       version: '0'
 requirements: []
 rubyforge_project: 
-rubygems_version: 2.2.0.rc.1
+rubygems_version: 2.2.2
 signing_key: 
 specification_version: 4
 summary: Gas Turbine engine and its components math models written in ruby
 test_files:
+- spec/models/gas_spec.rb
 - spec/models/simple/burner_spec.rb
 - spec/models/simple/compressor_spec.rb
-- spec/models/simple/gas_spec.rb
 - spec/models/simple/turbine_spec.rb
 - spec/spec_helper.rb