Index: mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/CarnotCycle/CarnotCycle_abg.fig ================================================================== --- mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/CarnotCycle/CarnotCycle_abg.fig +++ mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/CarnotCycle/CarnotCycle_abg.fig @@ -1,26 +1,36 @@ -#FIG 3.1 +#FIG 3.2 Portrait Center Metric +A4 +100.00 +Single +-2 1200 2 2 1 0 2 -1 7 0 0 -1 0.000 0 0 -1 0 0 3 3827 3152 5402 3152 5177 3377 2 1 0 2 -1 7 0 0 -1 0.000 0 0 -1 0 0 3 6977 3152 8552 3152 8327 3377 2 1 0 2 -1 7 0 0 -1 0.000 0 0 -1 0 0 2 8552 2927 8552 3377 2 1 0 2 -1 7 0 0 -1 0.000 0 0 -1 0 0 3 - 1352 3150 2927 3150 2702 3375 -2 4 0 1 31 7 1 0 -1 0.000 0 0 7 0 0 5 - 9900 5850 9900 450 0 450 0 5850 9900 5850 + 902 3150 2477 3150 2252 3375 +2 1 0 2 -1 7 0 0 -1 0.000 0 0 -1 0 0 2 + 2475 2925 2475 3375 2 1 0 2 -1 7 0 0 -1 0.000 0 0 -1 0 0 2 - 2925 2925 2925 3375 -4 1 -1 0 0 0 20 0.0000000 4 270 1440 6212 3242 Cycle:cycle\001 -4 0 -1 0 0 0 20 0.0000000 4 210 1110 8687 3242 SS:Work\001 -4 1 1 1 0 3 20 0.0000000 4 210 735 4456 3466 dS/dT\001 -4 1 1 1 0 3 20 0.0000000 4 210 180 4501 2926 T\001 -4 1 1 1 0 3 20 0.0000000 4 210 180 7921 2926 P\001 -4 1 1 1 0 3 20 0.0000000 4 210 765 7786 3466 dV/dT\001 -4 2 -1 0 0 0 20 0.0000000 4 210 990 1170 3285 SS:Heat\001 -4 1 1 1 0 3 30 0.0000000 4 390 2310 5085 4635 Carnot cycle\001 -4 1 -1 0 0 0 20 0.0000000 4 210 570 3375 3240 RT:r\001 + 3825 2925 3825 3375 +2 1 0 2 -1 7 0 0 -1 0.000 0 0 -1 0 0 3 + 3152 1125 3152 2700 2927 2475 +2 4 0 2 31 7 1 0 -1 0.000 0 0 7 0 0 5 + 9900 5175 -450 5175 -450 450 9900 450 9900 5175 +4 1 -1 0 0 0 20 0.0000 4 285 1440 6212 3242 Cycle:cycle\001 +4 0 -1 0 0 0 20 0.0000 4 210 1050 8687 3242 Sf:Work\001 +4 1 1 1 0 3 20 0.0000 4 210 735 4456 3466 dS/dT\001 +4 1 1 1 0 3 20 0.0000 4 195 180 4501 2926 T\001 +4 1 1 1 0 3 20 0.0000 4 195 180 7921 2926 P\001 +4 1 1 1 0 3 20 0.0000 4 210 765 7786 3466 dV/dT\001 +4 1 1 1 0 3 30 0.0000 4 405 2325 5085 4635 Carnot cycle\001 +4 2 -1 0 0 0 20 0.0000 4 210 975 810 3240 Se:Heat\001 +4 1 -1 0 0 0 20 0.0000 4 195 1005 3150 3240 FMRT:r\001 +4 1 4 100 0 18 14 0.0000 4 195 600 3150 2925 [mod]\001 +4 1 -1 0 0 0 20 0.0000 4 210 1170 3150 990 Sf:switch\001 Index: mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/CarnotCycle/CarnotCycle_input.txt ================================================================== --- mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/CarnotCycle/CarnotCycle_input.txt +++ mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/CarnotCycle/CarnotCycle_input.txt @@ -1,42 +1,42 @@ -# Numerical parameter file (CarnotCycle_input.txt) +# Input file (CarnotCycle_input.txt) # Generated by MTT at Mon Dec 8 20:05:30 GMT 1997 # %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% # %% Version control history # %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% # %% $Id$ # %% $Log$ +# %% Revision 1.1 2000/12/28 18:14:40 peterg +# %% To RCS +# %% # %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% alpha = 0.553; # Set the inputs if ((t>=0.0)&&(t<1.0)) #Cooling - u(1) = 300.0; # Source temperature - r = 1e-5; # Small r - isothermal - u(2) = -alpha; # Volume rate-of-change + carnotcycle__heat__u= 300.0; # Source temperature + carnotcycle__switch__u = 1e10; # Large conductance - isothermal + carnotcycle__work__u = -alpha; # Volume rate-of-change endif; if ((t>=1.0)&&(t<2.0)) #Compression - u(1) = 300.0; # Source temperature - r = 1e5; # Large resistance -- isentropic - u(2) = -(0.8-alpha); # Volume rate-of-change - TopTemp = x(1)/(m_g*c_v); + carnotcycle__heat__u = 300.0; # Source temperature + carnotcycle__switch__u = 0; # Small conductance -- isentropic + carnotcycle__work__u = -(0.8-alpha); # Volume rate-of-change endif; if ((t>=2.0)&&(t<3.0)) #Heating - u(1) = TopTemp; # Source temperature - r = 1e-5; # Small r - isothermal - u(2) = (0.8-alpha); # Volume rate-of-change + TopTemp = x(1)/(m_g*c_v); + carnotcycle__heat__u = TopTemp; # Source temperature + carnotcycle__switch__u = 1e10; # Large conductance - isothermal + carnotcycle__work__u = (0.8-alpha); # Volume rate-of-change endif; if (t>=3.0) #Cooling -% Temp = x(1)/(m_g*c_v); # - u(1) = 300.0; # Source temperature - r = 1e5; # Large resistance -- no entropy flow - u(2) = alpha; # Volume rate-of-change + carnotcycle__heat__u = 300.0; # Source temperature + carnotcycle__switch__u = 0; # Small conductance -- isentropic + carnotcycle__work__u = alpha; # Volume rate-of-change endif; - - Index: mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/CarnotCycle/CarnotCycle_lbl.txt ================================================================== --- mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/CarnotCycle/CarnotCycle_lbl.txt +++ mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/CarnotCycle/CarnotCycle_lbl.txt @@ -1,42 +1,51 @@ -%SUMMARY CarnotCycle: a simple closed thermodynamic cycle -%DESCRIPTION The Carnot cycle is a simple closed thermodynamic cycle -%DESCRIPTION with four parts: -%DESCRIPTION o Isentropic compression -%DESCRIPTION o Heat injection at constant temperature -%DESCRIPTION o Isentropic expansion -%DESCRIPTION o Heat extraction at constant temperature - -%VAR P_0 -%VAR T_0 -%VAR V_0 -%VAR TopTemp -%VAR alpha - -%% Label file for system CarnotCycle (CarnotCycle_lbl.txt) - -% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -% %% Version control history -% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -% %% $Id$ -% %% $Log$ -% %% Revision 1.1 1998/07/21 15:18:18 peterg -% %% Initial revision -% %% -% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% - -%% Each line should be of one of the following forms: -% a comment (ie starting with %) -% Component-name CR_name arg1,arg2,..argn -% blank - - - -% Component type Cycle +#SUMMARY CarnotCycle: a simple closed thermodynamic cycle +#DESCRIPTION The Carnot cycle is a simple closed thermodynamic cycle +#DESCRIPTION with four parts: +#DESCRIPTION o Isentropic compression +#DESCRIPTION o Heat injection at constant temperature +#DESCRIPTION o Isentropic expansion +#DESCRIPTION o Heat extraction at constant temperature + +#PAR P_0 +#PAR T_0 +#PAR V_0 +#PAR alpha + +#NOTPAR ideal_gas +#NOTPAR alpha +#NOTPAR TopTemp + +## Label file for system CarnotCycle (CarnotCycle_lbl.txt) + +# ############################################################### +# ## Version control history +# ############################################################### +# ## $Id$ +# ## $Log$ +# ## Revision 1.2 1998/08/10 16:40:07 peterg +# ## Added VARs and parametrs +# ## +# ## Revision 1.1 1998/07/21 15:18:18 peterg +# ## Initial revision +# ## +# ############################################################### + +## Each line should be of one of the following forms: +# a comment (ie starting with #) +# Component-name CR_name arg1,arg2,..argn +# blank + + + +# Component type Cycle cycle none ideal_gas;c_v;gamma_g;m_g -% Component type RT - r lin flow,r +# Component type RT + r lin flow,1 + +# Component type Se + Heat SS external -% Component type SS - Heat SS external,internal - Work SS internal,external +# Component type Sf + Work SS external + switch SS external Index: mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/CarnotCycle/CarnotCycle_numpar.txt ================================================================== --- mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/CarnotCycle/CarnotCycle_numpar.txt +++ mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/CarnotCycle/CarnotCycle_numpar.txt @@ -4,10 +4,13 @@ # %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% # %% Version control history # %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% # %% $Id$ # %% $Log$ +# %% Revision 1.2 2000/12/28 18:14:40 peterg +# %% To RCS +# %% # %% Revision 1.1 1998/03/04 11:49:01 peterg # %% Initial revision # %% # %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% @@ -15,13 +18,9 @@ P_0 = 1e5; V_0 = 1; T_0 = 300; # Parameters -c_v = 718.0; # Parameter c_v for CU -gamma_g = 1.4; # Parameter gamma for CU +c_v = 718.0; # Parameter c_v for CU +gamma_g = 1.4; # Parameter gamma for CU m_g = P_0*V_0/(T_0*(gamma_g-1)*c_v); # Parameter m for CU - - -r = 1.0; # Parameter r for CarnotCycle - Index: mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/CarnotCycle/CarnotCycle_sspar.r ================================================================== --- mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/CarnotCycle/CarnotCycle_sspar.r +++ mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/CarnotCycle/CarnotCycle_sspar.r @@ -4,41 +4,45 @@ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % Version control history % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % $Id$ % % $Log$ +% % Revision 1.2 2000/12/28 18:14:40 peterg +% % To RCS +% % % % Revision 1.1 1998/03/04 11:49:14 peterg % % Initial revision % % % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Set a pressure of 1 bar -P_0 := 10^5; +%P_0 := 10^5; % Unit initial volume -V_0 := 1; +%V_0 := 1; % Internal energy U_0 := P_0*V_0/(gamma_g-1); % Set initial temperature of 300k -T_0 := 300; +%T_0 := 300; % Deduce the mass of gas -m := U_0/(T_0*c_v); +m_g := U_0/(T_0*c_v); % Entropy S_0 := U_0/T_0; % Steady-state states MTTX1 := U_0; % CarnotCycle_cycle_gas (c) MTTX2 := V_0; % CarnotCycle_cycle_gas (c) -MTTX3 := S_0; % CarnotCycle_cycle_entropy (3) -MTTX4 := V_0; % CarnotCycle_cycle_volume (3) +MTTX3 := S_0; % CarnotCycle_cycle_volume (3) +MTTX4 := V_0; % CarnotCycle_cycle_entropy (3) % Steady-state inputs -MTTU1 := 0; % CarnotCycle (Heat) +MTTU1 := T_0; % CarnotCycle (Heat) MTTU2 := 0; % CarnotCycle (Work) +MTTU3 := 0; % Switch (conduction) ;;END; Index: mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/DieselCycle/DieselCycle_abg.fig ================================================================== --- mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/DieselCycle/DieselCycle_abg.fig +++ mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/DieselCycle/DieselCycle_abg.fig @@ -1,33 +1,41 @@ -#FIG 3.1 +#FIG 3.2 Portrait Center Metric +A4 +100.00 +Single +-2 1200 2 2 1 0 2 -1 7 0 0 -1 0.000 0 0 -1 0 0 3 2026 3151 3601 3151 3376 3376 2 1 0 2 -1 7 0 0 -1 0.000 0 0 -1 0 0 3 5176 3151 6751 3151 6526 3376 2 1 0 2 -1 7 0 0 -1 0.000 0 0 -1 0 0 2 6751 2926 6751 3376 2 1 0 2 -1 7 0 0 -1 0.000 0 0 -1 0 0 2 2026 2926 2026 3376 -2 4 0 1 31 7 1 0 -1 0.000 0 0 7 0 0 5 - 10350 5850 10350 450 675 450 675 5850 10350 5850 2 1 0 2 -1 7 0 0 -1 0.000 0 0 -1 0 0 2 7425 2925 7425 3375 2 1 0 2 -1 7 0 0 -1 0.000 0 0 -1 0 0 3 7426 3150 9001 3150 8776 3375 2 1 0 2 -1 7 0 0 -1 0.000 0 0 -1 0 0 3 - 7065 2924 7065 1349 7290 1574 -4 1 4 0 0 0 20 0.0000 4 240 735 3601 2926 [Heat]\001 -4 1 -1 0 0 0 20 0.0000 4 255 1290 4411 3241 Cycle:cycle\001 -4 1 4 0 0 0 20 0.0000 4 240 825 5176 2926 [Work]\001 -4 2 -1 0 0 0 20 0.0000 4 195 960 1936 3241 SS:Heat\001 -4 1 1 1 0 3 20 0.0000 4 210 690 2655 3465 dS/dT\001 -4 1 1 1 0 3 20 0.0000 4 195 165 2700 2925 T\001 -4 1 1 1 0 3 20 0.0000 4 195 165 6120 2925 P\001 -4 1 1 1 0 3 20 0.0000 4 210 720 5985 3465 dV/dT\001 -4 0 -1 0 0 0 20 0.0000 4 195 1050 9091 3241 SS:Work\001 -4 1 1 1 0 3 30 0.0000 4 405 2085 5535 4500 Diesel cycle\001 -4 1 -1 0 0 0 20 0.0000 4 195 135 7066 3241 1\001 -4 1 -1 0 0 0 20 0.0000 4 195 375 7021 1261 R:r\001 + 6975 2924 6975 1349 7200 1574 +2 1 0 2 -1 7 0 0 -1 0.000 0 0 -1 0 0 3 + 5176 1125 6751 1125 6526 1350 +2 4 0 2 31 7 1 0 -1 0.000 0 0 7 0 0 5 + 10350 4950 675 4950 675 450 10350 450 10350 4950 +4 1 4 0 0 0 20 0.0000 4 240 795 3601 2926 [Heat]\001 +4 1 -1 0 0 0 20 0.0000 4 285 1440 4411 3241 Cycle:cycle\001 +4 1 4 0 0 0 20 0.0000 4 255 900 5176 2926 [Work]\001 +4 2 -1 0 0 0 20 0.0000 4 210 945 1936 3241 Sf:Heat\001 +4 1 1 1 0 3 20 0.0000 4 210 735 2655 3465 dS/dT\001 +4 1 1 1 0 3 20 0.0000 4 195 180 2700 2925 T\001 +4 1 1 1 0 3 20 0.0000 4 195 180 6120 2925 P\001 +4 1 1 1 0 3 20 0.0000 4 210 765 5985 3465 dV/dT\001 +4 0 -1 0 0 0 20 0.0000 4 210 1080 9091 3241 Se:Work\001 +4 1 1 1 0 3 30 0.0000 4 405 2160 5535 4500 Diesel cycle\001 +4 1 -1 0 0 0 20 0.0000 4 195 150 7066 3241 1\001 +4 0 -1 0 0 0 20 0.0000 4 195 825 6930 1215 FMR:r\001 +4 1 4 0 0 0 20 0.0000 4 255 750 6705 945 [mod]\001 +4 2 -1 0 0 0 20 0.0000 4 210 1170 5085 1215 Sf:Switch\001 Index: mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/DieselCycle/DieselCycle_input.txt ================================================================== --- mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/DieselCycle/DieselCycle_input.txt +++ mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/DieselCycle/DieselCycle_input.txt @@ -11,32 +11,31 @@ # %% # %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% # Set the inputs - if ((t>=0.0)&&(t<1.0)) #Compression - dieselcycle__heat = 0.0; # Entropy flow - r = 1e10; # Big R to make a flow source - dieselcycle__work = 0.8*r; #- Volume rate-of-change - TopPressure = (gamma_g-1)*(x(1)/x(2)) + dieselcycle__Heat__u = 0.0; # Entropy flow + dieselcycle__switch__u = 1/big; # Large R to make a flow source + dieselcycle__Work__u = 0.8/big; #- Volume rate-of-change endif; if ((t>=1.0)&&(t<2.0)) #Heating - dieselcycle__heat = 1000; # Entropy flow - r = 1e4; # small r constant pressure - dieselcycle__work = TopPressure; # Pressure source + TopPressure = (gamma_g-1)*(x(1)/x(2)) + dieselcycle__Heat__u = 1000; # Entropy flow + dieselcycle__switch__u = big; # small r constant pressure + dieselcycle__Work__u = TopPressure; # Pressure source Volume = x(4); endif; if ((t>=2.0)&&(t<3.0)) #Expansion - dieselcycle__heat = 0.0; # Entropy flow - r = 1e10; # Big R to make a flow source - dieselcycle__work = -(1-Volume)*r; # Volume rate-of-change + dieselcycle__Heat__u = 0.0; # Entropy flow + dieselcycle__switch__u = 1/big; # Large R to make a flow source + dieselcycle__Work__u = -(1-Volume)*big; # Volume rate-of-change endif; if (t>=3.0) #Cooling Pressure = (gamma_g-1)*x(1)/x(2); - dieselcycle__heat = -1000*(Pressure>1e5); # Entropy flow - r = 1e10; # Big R to make a flow source - dieselcycle__work = 0.0*r; # Volume rate-of-change + dieselcycle__Heat__u = -1000*(Pressure>1e5); # Entropy flow + dieselcycle__switch__u = 1/big; # Large R to make a flow source + dieselcycle__Work__u = 0.0; # Zero volume rate-of-change endif; Index: mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/DieselCycle/DieselCycle_lbl.txt ================================================================== --- mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/DieselCycle/DieselCycle_lbl.txt +++ mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/DieselCycle/DieselCycle_lbl.txt @@ -9,10 +9,11 @@ #PAR P_0 #PAR T_0 #PAR V_0 #PAR S_0 #PAR U_0 +#PAR big ### NB TopPressure should be computed within input.txt !!! #PAR TopPressure #PAR Volume #PAR Pressure @@ -43,10 +44,15 @@ # Component type Cycle cycle none ideal_gas;c_v;gamma_g;m_g # Component type R - r lin flow,r + r lin flow,1 + +# Component type Sf + Heat SS external + +# Component type Se + Work SS external -# Component type SS - Heat SS internal,external - Work SS external,internal +# Component type Sf + Switch SS external Index: mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/DieselCycle/DieselCycle_numpar.txt ================================================================== --- mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/DieselCycle/DieselCycle_numpar.txt +++ mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/DieselCycle/DieselCycle_numpar.txt @@ -11,10 +11,12 @@ # %% # %% Revision 1.1 1998/03/04 11:45:49 peterg # %% Initial revision # %% # %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + +big = 1e5; # Large number # Initial states -- needed to choose an approppriate mass P_0 = 1e5; V_0 = 1; T_0 = 300; @@ -23,13 +25,12 @@ c_v = 718.0; # Parameter c_v for CU gamma_g = 1.4; # Parameter gamma_g for CU m_g = P_0*V_0/(T_0*(gamma_g-1)*c_v);# Parameter m for CU ## Removed by MTT on Wed Aug 6 17:32:59 BST 2003: pressure = 1.0; # Added by MTT on Wed Aug 06 17:31:21 BST 2003 -r = 1.0; # Added by MTT on Wed Aug 06 17:31:21 BST 2003 s_0 = 1.0; # Added by MTT on Wed Aug 06 17:31:21 BST 2003 ## Removed by MTT on Wed Aug 6 17:32:59 BST 2003: ## Removed by MTT on Wed Aug 6 17:32:59 BST 2003: toppressure = 1.0; # Added by MTT on Wed Aug 06 17:31:21 BST 2003 u_0 = 1.0; # Added by MTT on Wed Aug 06 17:31:21 BST 2003 ## Removed by MTT on Wed Aug 6 17:32:59 BST 2003: volume = 1.0; # Added by MTT on Wed Aug 06 17:31:21 BST 2003 pressure = 1.0; # Added by MTT on Wed Aug 06 17:36:47 BST 2003 toppressure = 1.0; # Added by MTT on Wed Aug 06 17:36:47 BST 2003 volume = 1.0; # Added by MTT on Wed Aug 06 17:36:47 BST 2003 Index: mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/DieselCycle/DieselCycle_sspar.r ================================================================== --- mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/DieselCycle/DieselCycle_sspar.r +++ mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/DieselCycle/DieselCycle_sspar.r @@ -4,27 +4,30 @@ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % Version control history % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % $Id$ % % $Log$ +% % Revision 1.1 2000/12/28 18:15:21 peterg +% % To RCS +% % % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Set a pressure of 1 bar -P_0 := 10^5; +%P_0 := 10^5; % Unit initial volume -V_0 := 1; +%V_0 := 1; % Internal energy U_0 := P_0*V_0/(gamma_g-1); % Set initial temperature of 300k -T_0 := 300; +%T_0 := 300; % Deduce the mass of gas -m := U_0/(T_0*c_v); +m_g := U_0/(T_0*c_v); % Entropy S_0 := U_0/T_0; % Steady-state states Index: mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/OttoCycle/OttoCycle_abg.fig ================================================================== --- mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/OttoCycle/OttoCycle_abg.fig +++ mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/OttoCycle/OttoCycle_abg.fig @@ -1,9 +1,13 @@ -#FIG 3.1 +#FIG 3.2 Portrait Center Metric +A4 +100.00 +Single +-2 1200 2 2 1 0 2 -1 7 0 0 -1 0.000 0 0 -1 0 0 3 2026 3151 3601 3151 3376 3376 2 1 0 2 -1 7 0 0 -1 0.000 0 0 -1 0 0 3 5176 3151 6751 3151 6526 3376 @@ -11,13 +15,13 @@ 6751 2926 6751 3376 2 1 0 2 -1 7 0 0 -1 0.000 0 0 -1 0 0 2 2026 2926 2026 3376 2 4 0 1 31 7 1 0 -1 0.000 0 0 7 0 0 5 8325 5850 8325 450 675 450 675 5850 8325 5850 -4 1 -1 0 0 0 20 0.0000000 4 270 1440 4411 3241 Cycle:cycle\001 -4 0 -1 0 0 0 20 0.0000000 4 210 1110 6886 3241 SS:Work\001 -4 2 -1 0 0 0 20 0.0000000 4 210 990 1936 3241 SS:Heat\001 -4 1 1 1 0 3 20 0.0000000 4 210 735 2655 3465 dS/dT\001 -4 1 1 1 0 3 20 0.0000000 4 210 180 2700 2925 T\001 -4 1 1 1 0 3 20 0.0000000 4 210 180 6120 2925 P\001 -4 1 1 1 0 3 20 0.0000000 4 210 765 5985 3465 dV/dT\001 -4 1 1 1 0 3 30 0.0000000 4 390 1800 4365 4770 Otto cycle\001 +4 1 -1 0 0 0 20 0.0000 4 285 1440 4411 3241 Cycle:cycle\001 +4 0 -1 0 0 0 20 0.0000 4 210 1050 6886 3241 Sf:Work\001 +4 2 -1 0 0 0 20 0.0000 4 210 945 1936 3241 Sf:Heat\001 +4 1 1 1 0 3 20 0.0000 4 210 735 2655 3465 dS/dT\001 +4 1 1 1 0 3 20 0.0000 4 195 180 2700 2925 T\001 +4 1 1 1 0 3 20 0.0000 4 195 180 6120 2925 P\001 +4 1 1 1 0 3 20 0.0000 4 210 765 5985 3465 dV/dT\001 +4 1 1 1 0 3 30 0.0000 4 405 1815 4365 4770 Otto cycle\001 Index: mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/OttoCycle/OttoCycle_input.txt ================================================================== --- mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/OttoCycle/OttoCycle_input.txt +++ mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/OttoCycle/OttoCycle_input.txt @@ -19,27 +19,27 @@ # %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% # Set the inputs if ((t>=0.0)&&(t<1.0)) #Compression - OttoCycle__Heat = 0.0; # Entropy flow - OttoCycle__Work = -0.8; # Volume rate-of-change + OttoCycle__Heat__u = 0.0; # Entropy flow + OttoCycle__Work__u = -0.8; # Volume rate-of-change endif; if ((t>=1.0)&&(t<2.0)) #Heating - OttoCycle__Heat = 1000; # Entropy flow - OttoCycle__Work = 0.0; # Volume rate-of-change + OttoCycle__Heat__u = 1000; # Entropy flow + OttoCycle__Work__u = 0.0; # Volume rate-of-change endif; if ((t>=2.0)&&(t<3.0)) #Expansion - OttoCycle__Heat = 0.0; # Entropy flow - OttoCycle__Work = 0.8; # Volume rate-of-change + OttoCycle__Heat__u = 0.0; # Entropy flow + OttoCycle__Work__u = 0.8; # Volume rate-of-change endif; if (t>=3.0) #Cooling - OttoCycle__Heat = -1000; # Entropy flow - OttoCycle__Work = 0.0; # Volume rate-of-change + OttoCycle__Heat__u = -1000; # Entropy flow + OttoCycle__Work__u = 0.0; # Volume rate-of-change endif; Index: mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/OttoCycle/OttoCycle_lbl.txt ================================================================== --- mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/OttoCycle/OttoCycle_lbl.txt +++ mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/OttoCycle/OttoCycle_lbl.txt @@ -42,8 +42,8 @@ # Component type Cycle cycle none ideal_gas;c_v;gamma_g;m_g -# Component type SS - Heat SS internal,external - Work SS internal,external +# Component type Sf + Heat SS external + Work SS external Index: mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/OttoCycle/OttoCycle_sspar.r ================================================================== --- mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/OttoCycle/OttoCycle_sspar.r +++ mttroot/mtt/lib/examples/Thermal/ThermodynamicCycles/OttoCycle/OttoCycle_sspar.r @@ -4,27 +4,30 @@ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % Version control history % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % $Id$ % % $Log$ +% % Revision 1.1 2000/12/28 18:15:52 peterg +% % To RCS +% % % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Set a pressure of 1 bar -P_0 := 10^5; +%%P_0 := 10^5; % Unit initial volume -V_0 := 1; +%%V_0 := 1; % Internal energy U_0 := P_0*V_0/(gamma_g-1); % Set initial temperature of 300k -T_0 := 300; +%%T_0 := 300; % Deduce the mass of gas -m := U_0/(T_0*c_v); +m_g := U_0/(T_0*c_v); % Entropy S_0 := U_0/T_0; % Steady-state states