pulo 0.1.1 → 0.1.2

data/lib/pulo/machine/steam/steam_process.rb

@@ -0,0 +1,27 @@
+
+module Pulo
+  module Steam
+    class SteamProcess
+
+      attr_reader :supply, :condensate
+      attr_reader :supply_pressure, :supply_quality, :process_power, :condensate_recovery
+      attr_reader :supply_massflow, :condensate_massflow, :condensate_power
+      def initialize(supply_pressure: nil, supply_quality:nil, process_power: nil, condensate_recovery: nil)
+        raise "Need all parameters" unless supply_pressure && supply_quality && process_power &&  condensate_recovery
+
+        @supply_pressure=supply_pressure
+        @supply_quality=supply_quality
+        @process_power=process_power
+        @condensate_recovery=condensate_recovery
+
+        @supply=WaterSteam.new(pressure: @supply_pressure, quality: @supply_quality)
+        @condensate=WaterSteam.new(pressure: supply_pressure, quality: Dimensionless.n(0))
+        @evap_energy=@supply.specific_enthalpy-@condensate.specific_enthalpy
+
+        @supply_massflow=@process_power/@evap_energy
+        @condensate_massflow=@supply_massflow*@condensate_recovery
+        @condensate_power=@condensate.specific_enthalpy*@condensate_massflow
+      end
+    end
+  end
+end

data/lib/pulo/machine/steam/steam_turbine.rb

@@ -0,0 +1,42 @@
+module Pulo
+  module Steam
+    class SteamTurbine
+      attr_reader :inlet_steam,:outlet_steam
+      attr_reader :isentropic_efficiency,:mechanical_efficiency
+      attr_reader :electrical_power
+
+      def initialize(inlet_pressure: nil, inlet_temperature:nil,outlet_pressure: nil,isentropic_efficiency: nil, mechanical_efficiency: nil,
+                    electrical_power: nil, mass_flow:nil)
+
+        raise "Need all parameters" unless isentropic_efficiency && mechanical_efficiency &&
+            inlet_pressure && outlet_pressure && inlet_temperature &&
+            (electrical_power || mass_flow)
+
+        @electrical_power=electrical_power
+        @isentropic_efficiency=isentropic_efficiency
+        @mechanical_efficiency=mechanical_efficiency
+
+        @inlet_steam=WaterSteam.new(pressure: inlet_pressure, temperature: inlet_temperature)
+        outlet_ideal=WaterSteam.new(pressure: outlet_pressure, specific_entropy: @inlet_steam.specific_entropy)
+        outlet_enthalpy=outlet_enthalpy_from_isentropic(@inlet_steam.specific_enthalpy,outlet_ideal.specific_enthalpy,@isentropic_efficiency)
+        @outlet_steam=WaterSteam.new(pressure: outlet_pressure,specific_enthalpy: outlet_enthalpy)
+        if @electrical_power
+          energy_output=@electrical_power/@mechanical_efficiency
+          outlet_steam.mass_flow=energy_output/(@inlet_steam.specific_enthalpy-outlet_enthalpy)
+        else
+          outlet_steam.mass_flow=mass_flow
+          energy_output=mass_flow*(@inlet_steam.specific_enthalpy-outlet_enthalpy)
+          @electrical_power=energy_output*@mechanical_efficiency
+        end
+        inlet_steam.mass_flow=outlet_steam.mass_flow
+      end
+
+      def isentropic_efficiency_from_enthalpy inlet_enthalpy,outlet_enthalpy,ideal_outlet_enthalpy
+        (inlet_enthalpy-outlet_enthalpy)/(inlet_enthalpy-ideal_outlet_enthalpy)
+      end
+      def outlet_enthalpy_from_isentropic inlet_enthalpy,ideal_outlet_enthalpy,isentropic_efficiency
+        inlet_enthalpy-isentropic_efficiency*(inlet_enthalpy-ideal_outlet_enthalpy)
+      end
+    end
+  end
+end

data/lib/pulo/machine/steam/water_steam.rb

@@ -0,0 +1,229 @@
+# encoding: utf-8
+
+require_relative '../../../../lib/pulo/machine/steam/if97'
+
+module Pulo
+  module Steam
+    class WaterSteam
+      attr_accessor :pressure,:temperature
+      attr_accessor :specific_volume,:specific_volume_liquid
+      attr_accessor :specific_internal_energy,:specific_entropy,:specific_enthalpy
+      attr_accessor :saturation_temperature, :saturation_pressure
+      attr_accessor :state,:if97_region
+      attr_accessor :quality
+      attr_reader :mass_flow, :energy_flow
+
+      def to_s(us: false)
+        r=case if97_region
+             when "1"
+               "Water "
+             when "2" || "3"
+               "Superheated "
+            when "4"
+              case quality
+                when Dimensionless.n(1)
+                  "Dry steam "
+                when Dimensionless.n(0)
+                  "Water "
+                else
+                  "Quality:#{quality} "
+              end
+
+             else
+               if97_region.to_s
+           end.ljust(12)
+        if us
+          r+="mf:#{mass_flow.pounds_per_hour.to_s(0).rjust(15)} ef:#{energy_flow.millions_btu_per_hour.to_s.rjust(17)} p:#{pressure.psig.to_s.rjust(9)}  t:#{temperature.fahrenheit.to_s(0).rjust(8)}  h:#{specific_enthalpy.btu_per_pound.to_s(0).rjust(15)} s:#{specific_entropy.btu_per_pound_rankine}"
+        else
+          r+="mf:#{mass_flow.kilograms_per_second.to_s(0).rjust(15)} ef:#{energy_flow.megawatts.to_s.rjust(17)} p:#{pressure.megapascals} t:#{temperature.kelvin} #{temperature.celsius} h:#{specific_enthalpy.kilojoules_per_kilogram} s:#{specific_entropy.kilojoules_per_kilogram_kelvin}"
+        end
+
+      end
+
+      def mass_flow=(val)
+        @mass_flow=val
+        if @specific_enthalpy
+          @energy_flow=@mass_flow*@specific_enthalpy
+        end
+      end
+      def energy_flow=(val)
+        @energy_flow=val
+        if @specific_enthalpy
+          @mass_flow=@energy_flow/@specific_enthalpy
+        end
+      end
+
+      def initialize(pressure: nil,temperature: nil,specific_entropy: nil, specific_enthalpy: nil, quality: nil)
+
+        return if !pressure && !temperature && !specific_entropy && !specific_enthalpy && !quality
+
+        @temperature=temperature.kelvin if temperature
+        @pressure=pressure.megapascals if pressure
+        @specific_entropy=specific_entropy.kilojoules_per_kilogram_kelvin if specific_entropy
+        @specific_enthalpy=specific_enthalpy.kilojoules_per_kilogram if specific_enthalpy
+        @quality=quality.n if quality
+        @mass_flow=MassFlow.new(0)
+        @energy_flow=Power.new(0)
+
+        if @temperature && (@temperature>Temperature.kelvin(1073.15) || @temperature<Temperature.kelvin(273.15))
+          raise "Temperature outside valid range 273.15K to 1073.15K"
+          #TODO: Implement R5
+        end
+        if @pressure && (@pressure<Pressure.megapascals(0) || @pressure>Pressure.megapascals(100))
+          raise "Pressure outside valid range 0MPa to 100MPa"
+        end
+
+        if (quality && (pressure || temperature))
+          @if97_region="4"
+          if temperature
+            @pressure=IF97.pressure_from_temperature_r4 @temperature
+          else
+            @temperature=IF97.temperature_from_pressure_r4 @pressure
+          end
+          fluid=WaterSteam.new
+          IF97.other_props_r1(@temperature,@pressure,fluid)
+
+          gas=WaterSteam.new
+          IF97.other_props_r2(@temperature,@pressure,gas)
+
+          @specific_enthalpy=fluid.specific_enthalpy+(gas.specific_enthalpy-fluid.specific_enthalpy)*@quality
+          @specific_entropy=fluid.specific_entropy+(gas.specific_entropy-fluid.specific_entropy)*@quality
+          @specific_volume=fluid.specific_volume+(gas.specific_volume-fluid.specific_volume)*@quality
+          @specific_internal_energy=fluid.specific_internal_energy+(gas.specific_internal_energy-fluid.specific_internal_energy)*@quality
+          return
+        end
+
+        if (temperature && !pressure && !specific_entropy  && !specific_enthalpy) || (!temperature && pressure && !specific_entropy && !specific_enthalpy)
+          #if97 region 4 saturated steam
+          @if97_region="4"
+          if temperature
+            @pressure=IF97.pressure_from_temperature_r4 @temperature
+          else
+            @temperature=IF97.temperature_from_pressure_r4 @pressure
+          end
+          @saturation_temperature=@temperature
+          @saturation_pressure=@pressure
+          return
+        end
+
+        if @temperature && @pressure
+          #determine region
+          if @temperature<Temperature.kelvin(623.15)
+            @saturation_pressure=IF97.pressure_from_temperature_r4 @temperature
+            @saturation_temperature=IF97.temperature_from_pressure_r4 @pressure
+            if @pressure<@saturation_pressure
+              @if97_region="2"
+              IF97.other_props_r2(@temperature,@pressure,self)
+            else
+              @if97_region="1"
+              IF97.other_props_r1(@temperature,@pressure,self)
+            end
+          else
+            p_b23=IF97.pressure_from_b23(@temperature)
+            if @pressure>p_b23
+              @if97_region="3"
+              #TODO: Implement R3 forward
+            else
+              @if97_region="2"
+              IF97.other_props_r2(temperature,@pressure,self)
+            end
+          end
+          return
+        end
+
+        if @pressure && @specific_enthalpy
+          region=IF97.get_region_ph(@pressure,@specific_enthalpy)
+          case region
+            when 1
+              #raise "Reverse for region 1 not implemented"
+              IF97.temperature_from_pressure_enthalpy_r1(@pressure,@specific_enthalpy,self)
+              IF97.other_props_r1(@temperature,@pressure,self)
+              @if97_region="1"
+            when 2
+              IF97.temperature_from_pressure_enthalpy_r2(@pressure,@specific_enthalpy,self)
+              IF97.other_props_r2(@temperature,@pressure,self)
+              @if97_region="2"
+            when 3
+              raise "Reverse for region 3 not implemented"
+              @if97_region="3"
+            when 4
+              @temperature=IF97.temperature_from_pressure_r4 @pressure
+              @if97_region="4"
+
+              fluid=WaterSteam.new
+              IF97.other_props_r1(@temperature,@pressure,fluid)
+
+              gas=WaterSteam.new
+              IF97.other_props_r2(@temperature,@pressure,gas)
+
+              @quality=(self.specific_enthalpy-fluid.specific_enthalpy)/(gas.specific_enthalpy-fluid.specific_enthalpy)
+              @specific_entropy=fluid.specific_entropy+(gas.specific_entropy-fluid.specific_entropy)*@quality
+              @specific_volume=fluid.specific_volume+(gas.specific_volume-fluid.specific_volume)*@quality
+              @specific_internal_energy=fluid.specific_internal_energy+(gas.specific_internal_energy-fluid.specific_internal_energy)*@quality
+          end
+          return
+        end
+
+        if @pressure && @specific_entropy
+          #get the r1,r2 boundary saturation entropy
+          region=IF97.get_region_ps(@pressure,@specific_entropy)
+
+          case region
+            when 1
+              IF97.temperature_from_pressure_entropy_r1(@pressure,@specific_entropy,self)
+              IF97.other_props_r1(@temperature,@pressure,self)
+              @if97_region="1"
+            when 2
+              IF97.temperature_from_pressure_entropy_r2(@pressure,@specific_entropy,self)
+              IF97.other_props_r2(@temperature,@pressure,self)
+              @if97_region="2"
+            when 3
+              raise "Reverse for region 3 not implemented"
+              @if97_region="3"
+            when 4
+              @temperature=IF97.temperature_from_pressure_r4 @pressure
+              @if97_region="4"
+
+              fluid=WaterSteam.new
+              IF97.other_props_r1(@temperature,@pressure,fluid)
+
+              gas=WaterSteam.new
+              IF97.other_props_r2(@temperature,@pressure,gas)
+
+              @quality=(self.specific_entropy-fluid.specific_entropy)/(gas.specific_entropy-fluid.specific_entropy)
+              @specific_enthalpy=fluid.specific_enthalpy+(gas.specific_enthalpy-fluid.specific_enthalpy)*@quality
+              @specific_volume=fluid.specific_volume+(gas.specific_volume-fluid.specific_volume)*@quality
+              @specific_internal_energy=fluid.specific_internal_energy+(gas.specific_internal_energy-fluid.specific_internal_energy)*@quality
+          end
+
+          case region
+            when 1
+              IF97.temperature_from_pressure_entropy_r1(@pressure,@specific_entropy,self)
+              @if97_region="1"
+            when 2
+              IF97.temperature_from_pressure_entropy_r2(@pressure,@specific_entropy,self)
+              @if97_region="2"
+            when 3
+              raise "Reverse for region 3 not implemented"
+              @if97_region="3"
+            when 4
+              @temperature=IF97.temperature_from_pressure_r4 @pressure
+              @if97_region="4"
+
+              fluid=WaterSteam.new
+              IF97.other_props_r1(@temperature,@pressure,fluid)
+
+              gas=WaterSteam.new
+              IF97.other_props_r2(@temperature,@pressure,gas)
+
+              @quality=(self.specific_entropy-fluid.specific_entropy)/(gas.specific_entropy-fluid.specific_entropy)
+              @specific_enthalpy=fluid.specific_enthalpy+(gas.specific_enthalpy-fluid.specific_enthalpy)*@quality
+              @specific_volume=fluid.specific_volume+(gas.specific_volume-fluid.specific_volume)*@quality
+              @specific_internal_energy=fluid.specific_internal_energy+(gas.specific_internal_energy-fluid.specific_internal_energy)*@quality
+            return
+          end
+        end
+      end
+    end
+  end
+end

data/lib/pulo/material/water.rb

@@ -0,0 +1,34 @@
+# encoding: utf-8
+
+module Pulo
+
+  class Water
+  #def self.bulk_modulus
+    #BulkModulus.new(2.15*10**9)
+  #end
+  #def self.expansion_coefficient
+  #  Dimensionless.n(0.000088)
+  #end
+  #def self.enthalpy_of_vaporization
+  #  40680
+  #end
+  def self.standard_density
+    Density.kilograms_per_cubic_meter(1000)
+  end
+  #def self.density(temperature,pressure)
+  #  tmp=temperature.is_a?(Temperature) ? temperature.value : temperature
+  #  pres=pressure.is_a?(Pressure) ? pressure.value : pressure
+  #  Density.new(999.8 / (1 + Water.expansion_coefficient * tmp) / (1 - (pres - 1105) / Water.bulk_modulus.value))
+  #end
+  #def self.boiling_point(pressure)
+  #  pres=pressure.is_a?(Pressure) ? pressure.value : pressure
+  #  k=GAS_CONSTANT/Water.enthalpy_of_vaporization
+  #  Temperature.new(1/((1/Temperature.new(100).kelvin)-k*Math.log(pres/Atmosphere.standard_pressure.value)),:kelvin)
+  #end
+  #def self.melting_point(pressure)
+  #  pres=pressure.is_a?(Pressure) ? pressure.value : pressure
+  #  Temperature.new(Math::E**(Math.log(Temperature.new(0).kelvin)+((-1.56*10**-6)/6030)*(pres-Atmosphere.standard_pressure.value)),:kelvin)
+  #end
+  end
+
+end

data/lib/pulo/quantity/dimension.rb

@@ -0,0 +1,63 @@
+# encoding: utf-8
+
+module Pulo
+  class Dimension
+  include Comparable
+  def initialize(spec)
+    @spec=spec #hash keyed on dimension name as a symbol eg :L
+    clean
+  end
+  def spec
+    @spec
+  end
+
+  #def is_base?
+  #  @spec.count==1 && @spec.first[1]==1
+  #end
+
+  def +(other_spec)
+    Dimension.new @spec.merge(other_spec.spec) {|_key,val1,val2| val1+val2}
+  end
+  def -(other_spec)
+    new_hash = Hash[other_spec.spec.map{|k,itm| [k,-itm] } ]
+    Dimension.new @spec.merge(new_hash) {|_key,val1,val2| val1+val2}
+  end
+  def *(val)
+    Dimension.new Hash[@spec.map{|itm| [itm[0],itm[1]*val]}]
+  end
+  def /(val)
+    Dimension.new Hash[@spec.map{|itm| [itm[0],itm[1]/val]}]
+  end
+
+  def ==(other_spec)
+    val=(@spec.merge(other_spec.spec) {|_key,val1,val2| val1-val2}).inject(0) {|sum,n| sum+n[1].abs}
+    val==0
+  end
+  def eql?(other_spec)
+    self==(other_spec)
+  end
+  def hash
+    @spec.to_s.hash
+  end
+  def to_s(supress_pretty=false)
+    sort_spec=@spec.sort_by { |item| item[1]*-1}
+
+    if supress_pretty
+      sort_spec.reduce('') do |ret,item|
+        ret+item[0].to_s + item[1].to_s
+      end
+    else
+      sort_spec.reduce('') do |ret,item|
+        ret+=item[0].to_s + (item[1]==1 ? '' : Pulo.super_digit(item[1])) + '.'
+      end[0..-2]
+    end
+
+  end
+
+  private
+
+  def clean
+    @spec=Hash[@spec.delete_if {|_key,value| value==0}.sort]
+  end
+end
+end

data/lib/pulo/quantity/numeric_overloads.rb

@@ -0,0 +1,217 @@
+# encoding: utf-8
+
+#Cos, Sin and Tan of angles
+class Numeric
+  def clamp(min, max)
+    if self < min
+      min
+    else
+      self > max ? max : self
+    end
+  end
+end
+
+module Math
+  class << self
+    alias :old_sqrt :sqrt
+    def sqrt(val)
+      if val.is_a?(Pulo::Quantity)
+        val.rt(2)
+      else
+        self.old_sqrt(val)
+      end
+    end
+
+    alias :old_cos :cos
+    def cos(angle)
+      if angle.is_a?(Pulo::Angle)
+        Pulo::Dimensionless.n(self.old_cos(angle.radians.value))
+      else
+        self.old_cos(angle)
+      end
+    end
+
+    alias :old_acos :acos
+    def acos(dimless)
+      if dimless.is_a?(Pulo::Dimensionless)
+        Pulo::Angle.radians(self.old_acos(dimless.n.value.clamp(-1,1)))
+      else
+        self.old_acos(dimless)
+      end
+    end
+
+    alias :old_sin :sin
+    def sin(angle)
+      if angle.is_a?(Pulo::Angle)
+        Pulo::Dimensionless.n(self.old_sin(angle.radians.value))
+      else
+        self.old_sin(angle)
+      end
+    end
+
+    alias :old_asin :asin
+    def asin(dimless)
+      if dimless.is_a?(Pulo::Dimensionless)
+        Pulo::Angle.radians(self.old_asin(dimless.n.value.clamp(-1,1)))
+      else
+        self.old_asin(dimless)
+      end
+    end
+
+    alias :old_tan :tan
+    def tan(angle)
+      if angle.is_a?(Pulo::Angle)
+        Pulo::Dimensionless.n(self.old_tan(angle.radians.value))
+      else
+        self.old_tan(angle)
+      end
+    end
+
+    alias :old_atan :atan
+    def atan(dimless)
+      if dimless.is_a?(Pulo::Dimensionless)
+        Pulo::Angle.radians(self.old_tan(dimless.n.value))
+      else
+        self.old_atan(dimless)
+      end
+    end
+
+  end
+end
+
+#Add and subtract work for Dimensionless
+class BigDecimal
+  alias :old_minus :-
+  def -(other)
+    if other.is_a?(Pulo::Dimensionless)
+      Pulo::Dimensionless.new(self-other.to_base_unit.value)
+    else
+      self.old_minus(other)
+    end
+  end
+end
+
+class Integer
+  alias :old_minus :-
+  def -(other)
+    if other.is_a?(Pulo::Dimensionless)
+      Pulo::Dimensionless.new(self-other.to_base_unit.value)
+    else
+      self.old_minus(other)
+    end
+  end
+end
+
+class Float
+  alias :old_minus :-
+  def -(other)
+    if other.is_a?(Pulo::Dimensionless)
+      Pulo::Dimensionless.new(self-other.to_base_unit.value)
+    else
+      self.old_minus(other)
+    end
+  end
+end
+
+class BigDecimal
+  alias :old_plus :+
+  def +(other)
+    if other.is_a?(Pulo::Dimensionless)
+      Pulo::Dimensionless.new(self+other.to_base_unit.value)
+    else
+      self.old_plus(other)
+    end
+  end
+end
+
+class Integer
+  alias :old_plus :+
+  def +(other)
+    if other.is_a?(Pulo::Dimensionless)
+      Pulo::Dimensionless.new(self+other.to_base_unit.value)
+    else
+      self.old_plus(other)
+    end
+  end
+end
+
+class Float
+  alias :old_plus :+
+  def +(other)
+    if other.is_a?(Pulo::Dimensionless)
+      Pulo::Dimensionless.new(self+other.to_base_unit.value)
+    else
+      self.old_plus(other)
+    end
+  end
+end
+
+#Division (Scalar/Quantity) gives (1/Quantity) * Scalar
+class BigDecimal
+  alias :old_div :/
+  def /(other)
+    if other.is_a?(Pulo::Quantity)
+      other.inverse*self
+    else
+      self.old_div(other)
+    end
+  end
+end
+
+class Integer
+  alias :old_div :/
+  def /(other)
+    if other.is_a?(Pulo::Quantity)
+      other.inverse*self
+    else
+      self.old_div(other)
+    end
+  end
+end
+
+class Float
+  alias :old_div :/
+  def /(other)
+    if other.is_a?(Pulo::Quantity)
+      other.inverse*self
+    else
+      self.old_div(other)
+    end
+  end
+end
+
+#Multiplication is a reversal for Quantity * Scalar
+class BigDecimal
+  alias :old_times :*
+  def *(other)
+    if other.is_a?(Pulo::Quantity)
+      other*self
+    else
+      self.old_times(other)
+    end
+  end
+end
+
+class Integer
+  alias :old_times :*
+  def *(other)
+    if other.is_a?(Pulo::Quantity)
+      other*self
+    else
+      self.old_times(other)
+    end
+  end
+end
+
+
+class Float
+  alias :old_times :*
+  def *(other)
+    if other.is_a?(Pulo::Quantity)
+      other*self
+    else
+      self.old_times(other)
+    end
+  end
+end
+