Changeset 1950

Timestamp:

12/05/2008 07:57:47 AM (5 weeks ago)

Author:

rfranke

Message:

add nTubes parameter to DistributedPipe

Location:

Modelica_Fluid/branches/StreamConnector/Modelica_Fluid

Files:

• Examples/AST_BatchPlant.mo (modified) (2 diffs)
• Pipes.mo (modified) (8 diffs)
• Test/TestOverdeterminedSteadyStateInit.mo (modified) (3 diffs)
• Volumes.mo (modified) (3 diffs)

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}},
@@ -2052 +2051 @@
               lineColor={0,0,0},
               textString=DynamicSelect(" ", realString(
-                    level, 
-                    1, 
+                    level,
+                    1,
                     3))),
             Line(

Modelica_Fluid/branches/StreamConnector/Modelica_Fluid/Pipes.mo

@@ -182 +182 @@
     extends Modelica_Fluid.Volumes.BaseClasses.PartialDistributedVolume(
       final n = nNodes,
-      Qs_flow=heatTransfer.Q_flow);
+      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
@@ -259 +263 @@
       final length=length,
       state=medium.state,
-      m_flow = 0.5*(m_flow[1:n]+m_flow[2:n+1])) "Edit heat transfer parameters"
-              annotation (editButton=true, Placement(transformation(extent={{-20,-5},{20,35}},  rotation=0)));
+      m_flow = 0.5*(m_flow[1:n]+m_flow[2:n+1])/nTubes)
+      "Edit heat transfer parameters" 
+        annotation (editButton=true, Placement(transformation(extent={{-20,-5},{20,35}},  rotation=0)));
 
   equation
@@ -278 +283 @@
 
     // Source/sink terms for mass and energy balances
-    fluidVolume=fill(V/n, n);
+    fluidVolume=fill(V/n*nTubes, n);
     Ws_flow=zeros(n);
     for i in 1:n loop
@@ -317 +322 @@
       fill(medium[1].p, n-1) = medium[2:n].p;
       if modelStructure == ModelStructure.a_v_b then
-        m_flow[1] = pressureDrop.m_flow[1];
+        m_flow[1] = pressureDrop.m_flow[1]*nTubes;
         flowState[1] = state_a;
         flowState[2] = medium[iLumped].state;
         flowState[3] = state_b;
-        m_flow[n+1] = pressureDrop.m_flow[2];
+        m_flow[n+1] = pressureDrop.m_flow[2]*nTubes;
       elseif modelStructure == ModelStructure.av_b then
         port_a.p = medium[1].p;
         flowState[1] = medium[iLumped].state;
         flowState[2] = state_b;
-        m_flow[n+1] = pressureDrop.m_flow[1];
+        m_flow[n+1] = pressureDrop.m_flow[1]*nTubes;
       elseif modelStructure == ModelStructure.a_vb then
-        m_flow[1] = pressureDrop.m_flow[1];
+        m_flow[1] = pressureDrop.m_flow[1]*nTubes;
         flowState[1] = state_a;
         flowState[2] = medium[iLumped].state;
@@ -351 +356 @@
         //m_flow[2:n+1] = pressureDrop.m_flow;
         for i in 2:n+1 loop
-          m_flow[i] = pressureDrop.m_flow[i-1];
+          m_flow[i] = pressureDrop.m_flow[i-1]*nTubes;
         end for;
         port_a.p = medium[1].p;
@@ -359 +364 @@
         //m_flow[1:n] = pressureDrop.m_flow;
         for i in 1:n loop
-          m_flow[i] = pressureDrop.m_flow[i];
+          m_flow[i] = pressureDrop.m_flow[i]*nTubes;
         end for;
         port_b.p = medium[n].p;
@@ -366 +371 @@
         //m_flow[2:n] = pressureDrop.m_flow[1:n-1];
         for i in 2:n loop
-          m_flow[i] = pressureDrop.m_flow[i-1];
+          m_flow[i] = pressureDrop.m_flow[i-1]*nTubes;
         end for;
         port_a.p = medium[1].p;
@@ -463 +468 @@
       //Geometry
       parameter SI.Length length "Length"   annotation(Dialog(tab="General", group="Geometry"));
-      parameter SI.Diameter diameter "Diameter of circular pipe"      annotation(Dialog(group="Geometry", enable=isCircular));
       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" 

Modelica_Fluid/branches/StreamConnector/Modelica_Fluid/Test/TestOverdeterminedSteadyStateInit.mo

@@ -13 +13 @@
     Pipes.DistributedPipe pipe(
       redeclare package Medium = Modelica.Media.Water.StandardWater,
-      nNodes=5,
       h_start=2e6,
       diameter=0.05,
@@ -19 +18 @@
       use_T_start=false,
       modelStructure=Modelica_Fluid.Types.ModelStructure.a_vb,
+      lumpedPressure=true,
       p_a_start=10000000,
       p_b_start=9900000,
-      lumpedPressure=true) 
+      nNodes=5) 
       annotation (Placement(transformation(extent={{-40,-10},{-20,10}})));
     ControlValves.ValveCompressible valve(
@@ -86 +86 @@
                         graphics),
                          Icon(coordinateSystem(preserveAspectRatio=false, extent={{-100,
-              -100},{100,100}}), graphics={Rectangle(extent={{-100,100},{100,-100}}, 
+              -100},{100,100}}), graphics={Rectangle(extent={{-100,100},{100,-100}},
               lineColor={0,0,255}), Text(
             extent={{-60,60},{60,-60}},

Modelica_Fluid/branches/StreamConnector/Modelica_Fluid/Volumes.mo

@@ -332 +332 @@
             lineColor={0,0,0},
             textString=DynamicSelect(" ", realString(
-                level, 
-                1, 
+                level,
+                1,
                 integer(precision)))),
           Line(
@@ -657 +657 @@
             lineColor={0,0,0},
             textString=DynamicSelect(" ", realString(
-                level, 
-                1, 
+                level,
+                1,
                 3))),
           Line(
@@ -1073 +1073 @@
 
     //Source terms, have to be set in inheriting class (to zero if not used)
-    protected
     SI.Volume[n] fluidVolume
         "Discretized volume, determine in inheriting class";