gtengine 0.0.1

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: d7a6124fbc51b2021b795bdaeeddf672e29bbc1f
+  data.tar.gz: 0e3e1fc6b8a5e2ce233fb83aeb16b319252abe99
+SHA512:
+  metadata.gz: ae805ae344231992d2974dccd6724bfaa38905048d129b2eb49d811ee0e1eaa36ec0a073f1c4f941a9104c953da6c11126f0faf3c3c76895f7c1ceaece386a9b
+  data.tar.gz: 030452f07c07cbf8878c756739b9f0e610c492c575443ab75b56af22baf4169afb4a5e0b5b794fde417b3ad280b4700cd5134b224195b59311f23aa5091b5fa6

data/.gitignore

@@ -0,0 +1,21 @@
+*.gem
+*.rbc
+.bundle
+.config
+.yardoc
+Gemfile.lock
+InstalledFiles
+_yardoc
+coverage
+doc/
+lib/bundler/man
+pkg
+rdoc
+spec/reports
+test/tmp
+test/version_tmp
+tmp
+
+# Custom ignores
+/.rvmrc
+/gogo.rb

data/.rspec

@@ -0,0 +1 @@
+--colour

data/.travis.yml

@@ -0,0 +1,7 @@
+language: ruby
+rvm:
+  - 2.0.0
+  - 2.1.0
+notifications:
+  email: false
+script: bundle exec rspec

data/Gemfile

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

data/Guardfile

@@ -0,0 +1,24 @@
+# A sample Guardfile
+# More info at https://github.com/guard/guard#readme
+
+guard :rspec do
+  watch(%r{^spec/.+_spec\.rb$})
+  watch(%r{^lib/(.+)\.rb$})     { |m| "spec/lib/#{m[1]}_spec.rb" }
+  watch('spec/spec_helper.rb')  { "spec" }
+
+  # Rails example
+  watch(%r{^app/(.+)\.rb$})                           { |m| "spec/#{m[1]}_spec.rb" }
+  watch(%r{^app/(.*)(\.erb|\.haml|\.slim)$})          { |m| "spec/#{m[1]}#{m[2]}_spec.rb" }
+  watch(%r{^app/controllers/(.+)_(controller)\.rb$})  { |m| ["spec/routing/#{m[1]}_routing_spec.rb", "spec/#{m[2]}s/#{m[1]}_#{m[2]}_spec.rb", "spec/acceptance/#{m[1]}_spec.rb"] }
+  watch(%r{^spec/support/(.+)\.rb$})                  { "spec" }
+  watch('config/routes.rb')                           { "spec/routing" }
+  watch('app/controllers/application_controller.rb')  { "spec/controllers" }
+
+  # Capybara features specs
+  watch(%r{^app/views/(.+)/.*\.(erb|haml|slim)$})     { |m| "spec/features/#{m[1]}_spec.rb" }
+
+  # Turnip features and steps
+  watch(%r{^spec/acceptance/(.+)\.feature$})
+  watch(%r{^spec/acceptance/steps/(.+)_steps\.rb$})   { |m| Dir[File.join("**/#{m[1]}.feature")][0] || 'spec/acceptance' }
+end
+

data/LICENSE.md

@@ -0,0 +1,22 @@
+Copyright (c) 2014 Oleg Bovykin
+
+MIT License
+
+Permission is hereby granted, free of charge, to any person obtaining
+a copy of this software and associated documentation files (the
+"Software"), to deal in the Software without restriction, including
+without limitation the rights to use, copy, modify, merge, publish,
+distribute, sublicense, and/or sell copies of the Software, and to
+permit persons to whom the Software is furnished to do so, subject to
+the following conditions:
+
+The above copyright notice and this permission notice shall be
+included in all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
+EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
+MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
+NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
+LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
+WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.

data/Rakefile

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

data/Readme.md

@@ -0,0 +1,14 @@
+[![Build Status](https://travis-ci.org/arrowcircle/engine.png?branch=master)](https://travis-ci.org/arrowcircle/engine)
+[![Code Climate](https://codeclimate.com/github/arrowcircle/engine.png)](https://codeclimate.com/github/arrowcircle/engine)
+
+# GTEngine
+
+Math models of gas turbine engine and its parts made with ruby language
+
+## Requirements
+
+	* Ruby 2.0.0 or newer
+
+### Copyright
+
+Copyright © __2014__ __Oleg Bovykin__. See [LICENSE]() for details.

data/gtengine.gemspec

@@ -0,0 +1,25 @@
+# coding: utf-8
+lib = File.expand_path('../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'gtengine/version'
+
+Gem::Specification.new do |spec|
+  spec.name          = "gtengine"
+  spec.version       = Gtengine::VERSION
+  spec.authors       = ["Oleg Bovykin"]
+  spec.email         = ["oleg.bovykin@gmail.com"]
+  spec.description   = %q{Gas Turbine engine and its components math models written in ruby}
+  spec.summary       = %q{Gas Turbine engine and its components math models written in ruby}
+  spec.homepage      = "https://github.com/arrowcircle/engine"
+  spec.license       = "MIT"
+
+  spec.files         = `git ls-files`.split($/)
+  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.3"
+  spec.add_development_dependency "rake"
+  spec.add_development_dependency "rspec"
+  spec.add_development_dependency "guard-rspec"
+end

data/lib/gtengine.rb

@@ -0,0 +1,7 @@
+require 'gtengine/version'
+require 'gtengine/gas'
+require 'gtengine/simple_cycle'
+
+module Gtengine
+  # Your code goes here...
+end

data/lib/gtengine/gas.rb

@@ -0,0 +1,55 @@
+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
+  end
+
+  def density
+    p / (t * R)
+  end
+
+  def k
+    k_from_cp cp
+  end
+
+  def average_k
+    k_from_cp average_cp
+  end
+
+  def cp
+    return low_temperature_cp if t < 750.0
+    high_temperature_cp
+  end
+
+  def average_cp
+    return low_temperature_average_cp if t < 700.0
+    high_temperature_average_cp
+  end
+
+  private
+
+    def k_from_cp c_p
+      c_p / (c_p - R)
+    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
+      (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 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

@@ -0,0 +1,50 @@
+class Gtengine::Simple::Burner
+
+  attr_accessor :input, :q_t, :output, :t_g
+
+  T0 = 288.3
+  QN = 43000000.0
+  ETA_G = 0.985
+  L0 = 14.7
+
+  def initialize input, t_g
+    @input = input
+    @t_g = t_g
+    @output = Gtengine::Gas.new(t_g, @input.p, @input.alfa)
+    cycle
+  end
+
+  def cp_vh
+    @input.cp
+  end
+
+  def cp_vyh
+    @output.cp
+  end
+
+  def t_vh
+    @input.t
+  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))
+  end
+
+  def cp_mult_t_0
+    Gtengine::Gas.new(T0, 101325.0, 1.0).cp * T0
+  end
+
+  def alfa
+    begin
+      1.0 / (q_ks * L0)
+    rescue
+      999999999.0
+    end
+  end
+
+  def cycle
+    5.times do
+      @output.alfa = alfa
+    end
+  end
+end

data/lib/gtengine/simple/compressor.rb

@@ -0,0 +1,65 @@
+class Gtengine::Simple::Compressor
+  attr_accessor :input, :output, :g, :pi_k, :kpd, :average
+
+  def initialize input, pi_k, kpd
+    @input = input
+    @pi_k = pi_k.to_f
+    @kpd = kpd.to_f
+    @average = Gtengine::Gas.new input.t.to_f, input.p.to_f
+    cycle
+  end
+
+  def t_vyh
+    @input.t * (1.0 + ((pi_k ** ((k - 1.0) / k) - 1.0) / kpd))
+  end
+
+  def p_vyh
+    @input.p * @pi_k
+  end
+
+  def l_k
+    @average.cp * t_vyh * (@pi_k ** ((k - 1.0) / k) - 1.0) / @kpd
+  end
+
+  def k
+    @average.average_k
+  end
+
+  def 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
+  end
+
+  def p_vh
+    @input.p
+  end
+
+  def cycle
+    5.times do
+      update_average
+    end
+    @output = Gtengine::Gas.new t_vyh, p_vyh
+  end
+
+  def info
+    puts "== Compressor"
+    puts "==== Pi_k: #{@pi_k}, КДП: #{@kpd}"
+    puts "==== Вход T: #{@input.t.to_i} K, P: #{@input.p.to_i} Па"
+    puts "====== Cp_sr: #{cp}, K: #{k}"
+    puts "==== Выход T: #{@output.t.to_i} K, P: #{@output.p.to_i} Па, ALFA: #{alfa}"
+    puts
+  end
+
+end

data/lib/gtengine/simple/pik_optimizer.rb

@@ -0,0 +1,35 @@
+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
+    @air = air
+    @start_pik = start_pik
+    @end_pik = end_pik
+    @t_g = t_g
+    @step = step
+    optimize
+  end
+
+  def optimize
+    @cycles = []
+    Range.new(start_pik, end_pik).step(step).each do |pik|
+      begin
+        c = Cycle.new(air, pik, t_g) 
+        @cycles << c if c.burner.q_ks
+      rescue
+        
+      end
+    end
+  end
+
+  def optimal
+    @cycles.sort {|a, b| a.burner.q_ks <=> b.burner.q_ks}.first
+  end
+
+  def info
+    @cycles.each do |c|
+      puts "pi_k: #{c.pi_k}, q_ks: #{c.q_ks}"
+    end
+    puts "Optimal pi_k: #{optimal.pi_k}, q_ks: #{optimal.q_ks}"
+  end
+end

data/lib/gtengine/simple/turbine.rb

@@ -0,0 +1,77 @@
+class Gtengine::Simple::Turbine
+
+  attr_accessor :burner, :l_k, :kpd, :average, :output
+
+  ETA = 0.9
+  ETA_M = 0.985
+  G_OHL = 0.001
+
+  def initialize burner, l_k, kpd = 0.9
+    @burner = burner
+    @l_k = l_k
+    @kpd = kpd.to_f
+    @average = Gtengine::Gas.new t_vh, p_vh
+    cycle
+  end
+
+  def input
+    @burner.output
+  end
+
+  def t_vh
+    burner.output.t
+  end
+
+  def p_vh
+    burner.output.p
+  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
+  end
+
+  def cycle
+    5.times do
+      update_average
+    end
+    @output = Gtengine::Gas.new t_vyh, p_vyh
+  end
+
+  def update_average
+    @average.t = (t_vh + t_vyh) / 2.0
+    @average.p = (p_vh + p_vyh) / 2.0
+  end
+
+  def t_vyh
+    t_vh * (1.0 - ((1.0 - (pi_t ** ((1.0 - k) / k))) * kpd))
+  end
+
+  def p_vyh
+    p_vh / pi_t
+  end
+
+  def l_t
+    cp * t_vh * ((1.0 - (pi_t ** ((1.0 - k) / k))) * kpd)
+  end
+
+  def k
+    @average.average_k
+  end
+
+  def cp
+    @average.average_cp
+  end
+
+  def info
+    puts "== Turbine"
+    puts "==== Pi_t: #{pi_t}, КДП: #{ETA}"
+    puts "==== Вход T: #{t_vh.to_i} K, P: #{p_vh.to_i} Па, ALFA: #{input.alfa}, Cp: #{input.cp}"
+    puts "====== Cp_sr: #{cp}, K: #{k}"
+    puts "==== Выход T: #{@output.t.to_i} K, P: #{@output.p.to_i} Па, , Cp: #{@output.cp}"
+    puts
+  end
+
+end

data/lib/gtengine/simple_cycle.rb

@@ -0,0 +1,29 @@
+module Gtengine
+  module Simple
+    class Cycle
+      require_relative 'simple/turbine'
+      require_relative 'simple/compressor'
+      require_relative 'simple/burner'
+      require_relative 'simple/pik_optimizer'
+
+      attr_accessor :air, :pi_k, :t_g, :compressor, :turbine, :burner
+
+      def initialize air=Gas.new(300, 101325), pi_k, t_g
+        @air = air
+        @pi_k = pi_k
+        @compressor = Compressor.new air, pi_k, 0.85
+        @burner = Burner.new @compressor.output, t_g
+        @turbine = Turbine.new @burner, @compressor.l_k, 0.9
+      end
+
+      def q_ks
+        @burner.q_ks
+      end
+
+      def info
+        puts "== Cycle info"
+        puts "pi_k: #{@compressor.pi_k}, q_ks: #{q_ks}"
+      end
+    end
+  end
+end

data/lib/gtengine/version.rb

@@ -0,0 +1,3 @@
+module Gtengine
+  VERSION = "0.0.1"
+end

data/spec/models/simple/burner_spec.rb

@@ -0,0 +1,15 @@
+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 }
+
+  it "считает расход топлива" do
+    expect(burner.q_ks).to eq 0.026948712923834613
+  end
+
+  it "считает альфа" do
+    expect(burner.alfa).to eq 2.524321331286569
+  end
+end

data/spec/models/simple/compressor_spec.rb

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

data/spec/models/simple/gas_spec.rb

@@ -0,0 +1,39 @@
+require 'spec_helper'
+
+describe Gtengine::Gas do
+  let(:gas) { Gtengine::Gas.new }
+
+  it "Считает плотность" do
+    expect(gas.density).to eq 1.259133611691023
+  end
+
+  it "Считает k" do
+    expect(gas.k).to eq 1.4029968590096924
+  end
+
+  it "Считает k за процесс" do
+    expect(gas.average_k).to eq 1.404058445015607
+  end
+
+  context "LOW TEMPERATURE" do
+    it "Считает теплоемкость" do
+      expect(gas.cp).to eq 1000.556923124028
+    end
+
+    it "Считает теплоемкость за процесс" do
+      expect(gas.average_cp).to eq 998.6832401978357
+    end
+  end
+
+  context "HIGH TEMPERATURE" do
+    let(:gas) { Gtengine::Gas.new 900, 1013251 }
+
+    it "Считает теплоемкость" do
+      expect(gas.cp).to eq 1121.2960248302861
+    end
+
+    it "Считает теплоемкость за процесс" do
+      expect(gas.average_cp).to eq 1037.119900706556
+    end
+  end
+end

data/spec/models/simple/turbine_spec.rb

@@ -0,0 +1,25 @@
+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 }
+
+  it "считает степень понижения давления" do
+    expect(turbine.pi_t).to eq 1.4858595860294537
+  end
+
+  it "считает температуру выхода" do
+    expect(turbine.t_vyh).to eq 815.5179703599931
+  end
+
+  it "считает давление выхода" do
+    expect(turbine.p_vyh).to eq 681928.5008670494
+  end
+
+  it "считает работу турбины l_t" do
+    expect(turbine.l_t).to eq 87289.27555567237
+  end
+  
+end

data/spec/spec_helper.rb

@@ -0,0 +1,6 @@
+require 'bundler/setup'
+require 'rspec'
+
+require 'gtengine'
+Dir[Pathname.new(File.expand_path('../', __FILE__)).join('support/**/*.rb')].each {|f| require f}
+

metadata

@@ -0,0 +1,128 @@
+--- !ruby/object:Gem::Specification
+name: gtengine
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+platform: ruby
+authors:
+- Oleg Bovykin
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2014-01-12 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.3'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.3'
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: rspec
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: guard-rspec
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+description: Gas Turbine engine and its components math models written in ruby
+email:
+- oleg.bovykin@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- ".rspec"
+- ".travis.yml"
+- Gemfile
+- Gemfile.lock
+- Guardfile
+- LICENSE.md
+- Rakefile
+- Readme.md
+- gtengine.gemspec
+- lib/gtengine.rb
+- lib/gtengine/gas.rb
+- lib/gtengine/simple/burner.rb
+- lib/gtengine/simple/compressor.rb
+- lib/gtengine/simple/pik_optimizer.rb
+- lib/gtengine/simple/turbine.rb
+- lib/gtengine/simple_cycle.rb
+- lib/gtengine/version.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
+licenses:
+- MIT
+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.0.rc.1
+signing_key: 
+specification_version: 4
+summary: Gas Turbine engine and its components math models written in ruby
+test_files:
+- 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