Changeset 1951

Timestamp:

12/05/2008 10:23:47 AM (5 weeks ago)

Author:

rfranke

Message:

add nTubes parameter to all pipe models

Location:

Modelica_Fluid/branches/StreamConnector/Modelica_Fluid

Files:

• Examples/AST_BatchPlant.mo (modified) (2 diffs)
• Pipes.mo (modified) (28 diffs)
• Volumes.mo (modified) (1 diff)

Modelica_Fluid/branches/StreamConnector/Modelica_Fluid/Examples/AST_BatchPlant.mo

@@ -1291 +1291 @@
               fillColor={85,170,255},
               fillPattern=FillPattern.Solid),
-            Line(points={{-200,100},{-200,-100},{0,-100},{0,100}}, color={0,0,0}),
+            Line(points={{-200,100},{-200,-100},{0,-100},{0,100}}, color={0,0,0}), 
+
             Text(
               extent={{-198,74},{0,38}},
@@ -2051 +2052 @@
               lineColor={0,0,0},
               textString=DynamicSelect(" ", realString(
-                    level,
-                    1,
+                    level, 
+                    1, 
                     3))),
             Line(

Modelica_Fluid/branches/StreamConnector/Modelica_Fluid/Pipes.mo

@@ -4 +4 @@
 
   model StaticPipe "Basic pipe flow model without storage of mass or energy"
+
     extends Modelica_Fluid.Pipes.BaseClasses.PartialPipe(
-          redeclare model HeatTransfer=BaseClasses.HeatTransfer.PipeHT_ideal);
+          redeclare model HeatTransfer=BaseClasses.HeatTransfer.PipeHT_none(nTubes=nTubes));
 
     replaceable PressureDrop pressureDrop(
@@ -17 +18 @@
             final p_a_start=p_a_start,
             final p_b_start=p_b_start,
-            final m_flow_start=m_flow_start,
+            final m_flow_start=m_flow_start/nTubes,
             final roughness=roughness,
             diameter=4*crossArea/perimeter,
@@ -27 +28 @@
             rotation=0)));
   equation
-    port_a.m_flow = pressureDrop.m_flow[1];
+    port_a.m_flow = pressureDrop.m_flow[1]*nTubes;
     0 = port_a.m_flow + port_b.m_flow;
     port_a.h_outflow = inStream(port_b.h_outflow);
@@ -51 +52 @@
 
     // Extend here to get right ordering in parameter box
-    extends Modelica_Fluid.Pipes.BaseClasses.PartialPipe(
-      final isCircular=true,
-      final perimeter=Modelica.Constants.pi*diameter,
-      final crossArea=Modelica.Constants.pi*diameter*diameter/4);
+    extends Modelica_Fluid.Pipes.BaseClasses.PartialPipe;
 
      //Initialization
@@ -78 +76 @@
       redeclare final package Medium = Medium,
       final n=1,
+      final nTubes=nTubes,
       diameter=4*crossArea/perimeter,
       area=perimeter*length,
@@ -83 +82 @@
       final length=length,
       state={volume.medium.state},
-      m_flow = {0.5*(port_a.m_flow - port_b.m_flow)},
+      m_flow = {0.5*(port_a.m_flow - port_b.m_flow)}/nTubes,
       final useFluidHeatPort=true) "Edit heat transfer parameters" 
       annotation (editButton=true, Placement(transformation(extent={{-11,14},{
@@ -91 +90 @@
       redeclare package Medium = Medium,
       allowFlowReversal=allowFlowReversal,
+      nTubes=nTubes,
       length=length/2,
+      roughness=roughness,
+      diameter=diameter,
+      perimeter=perimeter,
+      crossArea=crossArea,
       height_ab=height_ab/2,
-      roughness=roughness,
       m_flow_start=m_flow_start,
-      diameter=diameter,
       redeclare final model PressureDrop = PressureDrop) 
       annotation (Placement(transformation(extent={{-60,-40},{-40,-20}},
@@ -117 +119 @@
       redeclare package Medium = Medium,
       allowFlowReversal=allowFlowReversal,
+      nTubes=nTubes,
       length=length/2,
+      roughness=roughness,
+      diameter=diameter,
+      perimeter=perimeter,
+      crossArea=crossArea,
       height_ab=height_ab/2,
-      roughness=roughness,
       m_flow_start=m_flow_start,
-      diameter=diameter,
       redeclare final model PressureDrop = PressureDrop)   annotation (Placement(transformation(extent={{40,-40},
               {60,-20}},         rotation=0)));
@@ -184 +189 @@
       Qs_flow = heatTransfer.Q_flow*nTubes);
 
-    // Note: define nTubes as Real to support inverse calculations
-    parameter Real nTubes(min=1)=1 "Number of parallel tubes" 
-     annotation(Dialog(group="Geometry"));
-
     // extending PartialPipe
     extends Modelica_Fluid.Pipes.BaseClasses.PartialPipe(
@@ -231 +232 @@
             final p_a_start=p_a_start,
             final p_b_start=p_b_start,
-            final m_flow_start=m_flow_start,
+            final m_flow_start=m_flow_start/nTubes,
             final roughness=roughness,
             diameter=4*crossArea/perimeter,
@@ -258 +259 @@
       redeclare each final package Medium = Medium,
       final n=n,
+      final nTubes=nTubes,
       diameter=4*crossArea/perimeter,
       area=perimeter*length,
@@ -283 +285 @@
 
     // Source/sink terms for mass and energy balances
-    fluidVolume=fill(V/n*nTubes, n);
+    fluidVolume=fill(V/n, n);
     Ws_flow=zeros(n);
     for i in 1:n loop
@@ -347 +349 @@
         flowState[2:n+1] = medium[1:n].state;
         flowState[n+2] = state_b;
-        //m_flow = pressureDrop.m_flow;
+        //m_flow = pressureDrop.m_flow*nTubes;
         for i in 1:n+1 loop
-          m_flow[i] = pressureDrop.m_flow[i];
+          m_flow[i] = pressureDrop.m_flow[i]*nTubes;
         end for;
       elseif modelStructure == ModelStructure.av_b then
         flowState[1:n] = medium[1:n].state;
         flowState[n+1] = state_b;
-        //m_flow[2:n+1] = pressureDrop.m_flow;
+        //m_flow[2:n+1] = pressureDrop.m_flow*nTubes;
         for i in 2:n+1 loop
           m_flow[i] = pressureDrop.m_flow[i-1]*nTubes;
@@ -362 +364 @@
         flowState[1] = state_a;
         flowState[2:n+1] = medium[1:n].state;
-        //m_flow[1:n] = pressureDrop.m_flow;
+        //m_flow[1:n] = pressureDrop.m_flow*nTubes;
         for i in 1:n loop
           m_flow[i] = pressureDrop.m_flow[i]*nTubes;
@@ -369 +371 @@
       else // avb
         flowState[1:n] = medium[1:n].state;
-        //m_flow[2:n] = pressureDrop.m_flow[1:n-1];
+        //m_flow[2:n] = pressureDrop.m_flow[1:n-1]*nTubes;
         for i in 2:n loop
           m_flow[i] = pressureDrop.m_flow[i-1]*nTubes;
@@ -455 +457 @@
       extends Modelica_Fluid.Interfaces.PartialTwoPort;
 
-       //Initialization
+      // Initialization
       parameter Medium.AbsolutePressure p_a_start=system.p_start
         "Start value of pressure at port a" 
@@ -466 +468 @@
          annotation(Evaluate=true, Dialog(tab = "Initialization"));
 
-      //Geometry
+      // Geometry
+
+      // Note: define nTubes as Real to support inverse calculations
+      parameter Real nTubes(min=1)=1 "Number of parallel tubes" 
+        annotation(Dialog(group="Geometry"));
       parameter SI.Length length "Length"   annotation(Dialog(tab="General", group="Geometry"));
-      parameter SI.Length roughness(min=0)=2.5e-5
-        "Average height of surface asperities (default = smooth steel pipe)" 
-          annotation(Dialog(group="Geometry",enable=WallFriction.use_roughness));
-      parameter SI.Diameter diameter "Diameter of circular pipe"      annotation(Dialog(group="Geometry", enable=isCircular));
       parameter Boolean isCircular=true
         "= true if cross sectional area is circular" 
         annotation (Evaluate, Dialog(tab="General", group="Geometry"));
+      parameter SI.Diameter diameter "Diameter of circular pipe"      annotation(Dialog(group="Geometry", enable=isCircular));
       parameter SI.Length perimeter=Modelica.Constants.pi*diameter
         "Inner perimeter"                                                                                       annotation(Dialog(tab="General", group="Geometry", enable=not isCircular));
       parameter SI.Area crossArea=Modelica.Constants.pi*diameter*diameter/4
         "Inner cross section area"            annotation(Dialog(tab="General", group="Geometry", enable=not isCircular));
-      final parameter SI.Volume V=crossArea*length "volume size";
+      final parameter SI.Volume V=crossArea*length*nTubes "volume size";
+      parameter SI.Length roughness(min=0)=2.5e-5
+        "Average height of surface asperities (default = smooth steel pipe)" 
+          annotation(Dialog(group="Geometry",enable=WallFriction.use_roughness));
 
       // Static head
@@ -785 +791 @@
       replaceable package Medium=Modelica.Media.Interfaces.PartialMedium;
       parameter Integer n=1 "Number of heat transfer segments";
+      parameter Real nTubes "Number of parallel tubes, used for HeatPorts only";
       parameter SI.Length length "Pipe length";
       parameter SI.Length diameter
@@ -796 +803 @@
 
       // Output defined by heat transfer model
-      output SI.HeatFlowRate[n] Q_flow "Heat flow rates";
+      output SI.HeatFlowRate[n] Q_flow "Heat flow rates per tube";
 
       // Heat ports
@@ -821 +828 @@
     equation
       T = Medium.temperature(state);
+      wallHeatPort.Q_flow = Q_flow*nTubes;
       if useFluidHeatPort then
-        prescribedHeatFlow.Q_flow = Q_flow;
+        prescribedHeatFlow.Q_flow = Q_flow*nTubes;
       else
         prescribedHeatFlow.port.T = T;
       end if;
+
+      connect(prescribedHeatFlow.port, fluidHeatPort) annotation (Line(
+          points={{-1.83697e-015,-20},{0,-20},{0,-60}},
+          color={191,0,0},
+          smooth=Smooth.None));
+
       annotation (Icon(coordinateSystem(preserveAspectRatio=true,  extent={{-100,
                   -100},{100,100}}), graphics={Ellipse(
@@ -840 +854 @@
 Base class for heat transfer models that can be used in distributed pipe models.
 </html>"),
-        Diagram(coordinateSystem(preserveAspectRatio=true, extent={{-100,-100},
-                  {100,100}}),
+        Diagram(coordinateSystem(preserveAspectRatio=true, extent={{-100,-100},{
+                100,100}}),
                         graphics));
-      connect(prescribedHeatFlow.port, fluidHeatPort) annotation (Line(
-          points={{-1.83697e-015,-20},{0,-20},{0,-60}},
-          color={191,0,0},
-          smooth=Smooth.None));
     end PartialPipeHeatTransfer;
 
@@ -867 +877 @@
       Re = CharacteristicNumbers.ReynoldsNumber(m_flow, diameter, crossArea, eta);
       Nu = CharacteristicNumbers.NusseltNumber(alpha, diameter, lambda);
-      wallHeatPort.Q_flow=Q_flow;
-      wallHeatPort.Q_flow=alpha.*area.*(wallHeatPort.T - T);
+      Q_flow=alpha.*area.*(wallHeatPort.T - T);
         annotation (Documentation(info="<html>
 Base class for heat transfer models that are expressed in terms of the Nusselt number and which can be used in distributed pipe models.
 </html>"));
     end PartialPipeHT_Nu;
+
+    model PipeHT_none
+        "PipeHT_none: No heat transfer assuming perfect isolation"
+      extends PartialPipeHeatTransfer;
+    equation
+      Q_flow = zeros(n);
+      annotation(Documentation(info="<html>
+Ideal heat transfer without thermal resistance.
+</html>"));
+    end PipeHT_none;
 
     model PipeHT_ideal
@@ -879 +898 @@
     equation
       T = wallHeatPort.T;
-      Q_flow = wallHeatPort.Q_flow;
       annotation(Documentation(info="<html>
 Ideal heat transfer without thermal resistance.
@@ -894 +912 @@
 </html>"));
     equation
-      wallHeatPort.Q_flow = alpha0*area/n*(wallHeatPort.T - T);
-      wallHeatPort.Q_flow = Q_flow;
+      Q_flow = alpha0*area/n*(wallHeatPort.T - T);
     end PipeHT_constAlpha;
     annotation (Documentation(info="<html>
@@ -1046 +1063 @@
         redeclare final package Medium = Medium,
         final n=1,
+        final nTubes=1,
         diameter=4*crossArea/perimeter,
         area=perimeter*length,
@@ -1188 +1206 @@
        redeclare final package Medium = Medium,
        final n=nNodes,
+       final nTubes=1,
        diameter=4*crossArea/perimeter,
        area=perimeter*length,
@@ -1387 +1406 @@
       redeclare each final package Medium = Medium,
       final n=n,
+      final nTubes=1,
       diameter=4*crossArea/perimeter,
       area=perimeter*length,
@@ -1649 +1669 @@
               fillColor={0,0,0},
               fillPattern=FillPattern.Solid)}),
-    Diagram(coordinateSystem(preserveAspectRatio=true,  extent={{-100,-100},{100,
-                100}},
+    Diagram(coordinateSystem(preserveAspectRatio=true,  extent={{-100,-100},{
+                100,100}},
             grid={1,1}),
             graphics),

Modelica_Fluid/branches/StreamConnector/Modelica_Fluid/Volumes.mo

@@ -657 +657 @@
             lineColor={0,0,0},
             textString=DynamicSelect(" ", realString(
-                level,
-                1,
+                level, 
+                1, 
                 3))),
           Line(