/[ascend]/trunk/models/atoms.a4l
ViewVC logotype

Contents of /trunk/models/atoms.a4l

Parent Directory Parent Directory | Revision Log Revision Log


Revision 288 - (show annotations) (download) (as text)
Fri Feb 10 00:52:07 2006 UTC (16 years, 7 months ago) by johnpye
File MIME type: text/x-ascend
File size: 27045 byte(s)
Changed 'heat_capacity' to 'molar_heat_capacity'
1 REQUIRE "system.a4l";
2 (* => system.a4l, basemodel.a4l *)
3 REQUIRE "measures.a4l";
4 (* => measures.a4l *)
5 PROVIDE "atoms.a4l";
6
7 (*
8 * atoms.a4l
9 * by Joseph J. Zaher, Benjamin A. Allan, Robert S. Huss
10 * Part of the ASCEND Library
11 * $Date: 1998/08/17 13:07:11 $
12 * $Revision: 1.11 $
13 * $Author: ballan $
14 * $Source: /afs/cs.cmu.edu/project/ascend/Repository/models/atoms.a4l,v $
15 *
16 * This file is part of the ASCEND Modeling Library.
17 *
18 * Copyright (C) 1994-1998 Carnegie Mellon University
19 *
20 * The ASCEND Modeling Library is free software; you can redistribute
21 * it and/or modify it under the terms of the GNU General Public
22 * License as published by the Free Software Foundation; either
23 * version 2 of the License, or (at your option) any later version.
24 *
25 * The ASCEND Modeling Library is distributed in hope that it
26 * will be useful, but WITHOUT ANY WARRANTY; without even the implied
27 * warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
28 * See the GNU General Public License for more details.
29 *
30 * You should have received a copy of the GNU General Public License
31 * along with the program; if not, write to the Free Software
32 * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139 USA.
33 *)
34
35 (*============================================================================*
36
37 A T O M S . A 4 L
38 -------------------------
39
40 AUTHOR: Joseph J. Zaher
41
42 DATES: 07/91 - Original code.
43 02/92 - Eliminated "display_units" assignments, added
44 commonly used UNIVERSAL constants and transport
45 quantities.
46 04/92 - Added electro-magnetic quantities.
47 06/92 - Converted by Ben Allan to be compatible with ASCEND3C
48 08/92 - Additional monetary atoms added by Bob Huss
49 02/96 - Brought proper constants into the library. Ben Allan.
50 07/96 - Incorporated Ken Tyner's substitution atoms. BAA.
51 07/96 - Switched to REQUIRE syntax.
52 08/96 - Switched to REQUIRE ivpsystem.a4l.
53 08/97 - Switched to REQUIRE/PROVIDE syntax.
54
55 CONTENTS: ASCEND atom definitions for engineering variable types.
56 Many of the anticipated dimensional variables which occur
57 in engineering design calculations are given to provide a
58 means of standardization. Chosen defaults, nominal, and
59 lower and upper bound values should be re-specified if
60 necessary to enhance the convergence properties of
61 specific models. Units to be displayed are to be controlled
62 using the UNITS tool kit of the environment.
63
64 REQUIRES: see top of file.
65 *===========================================================================*)
66
67 (* G E N E R I C C O N S T A N T S ( C H E M . E . B I A S )
68 -------------------------------------------------------------------- *)
69
70 CONSTANT constant
71 (* any sloppiness about what is a constant will yield a wild real *)
72 REFINES real_constant;
73
74 (* dimensionless *)
75
76 CONSTANT critical_compressibility REFINES real_constant DIMENSIONLESS;
77
78 CONSTANT acentric_factor REFINES real_constant DIMENSIONLESS;
79
80 CONSTANT UNIFAC_size REFINES real_constant DIMENSIONLESS;
81
82 CONSTANT Wilson_constant REFINES real_constant DIMENSIONLESS;
83
84 CONSTANT vapor_pressure_constant REFINES real_constant;
85
86 CONSTANT factor_constant REFINES real_constant DIMENSIONLESS;
87
88 (* molecular weight *)
89 CONSTANT molar_weight_constant
90 REFINES real_constant DIMENSION M/Q;
91
92 (* atomic mass *)
93 CONSTANT atomic_mass_constant
94 REFINES real_constant DIMENSION M;
95
96 (* temperatures *)
97 CONSTANT temperature_constant
98 REFINES real_constant DIMENSION TMP;
99
100 CONSTANT boiling_temperature
101 REFINES temperature_constant;
102
103 CONSTANT critical_temperature
104 REFINES temperature_constant;
105
106 CONSTANT reference_temperature
107 REFINES temperature_constant;
108
109 CONSTANT UNIFAC_a
110 REFINES temperature_constant;
111
112 (* pressures *)
113 CONSTANT pressure_constant
114 REFINES real_constant DIMENSION M/L/T^2;
115
116 CONSTANT critical_pressure
117 REFINES pressure_constant;
118
119 CONSTANT reference_pressure
120 REFINES pressure_constant;
121
122 (* molar volumes *)
123 CONSTANT molar_volume_constant
124 REFINES real_constant DIMENSION L^3/Q;
125
126 CONSTANT critical_volume
127 REFINES molar_volume_constant;
128
129 CONSTANT reference_molar_volume
130 REFINES molar_volume_constant;
131
132 (* mass densities *)
133 CONSTANT reference_mass_density
134 REFINES real_constant DIMENSION M/L^3;
135
136 (* molar energies *)
137 CONSTANT molar_energy_constant
138 REFINES real_constant DIMENSION M*L^2/T^2/Q;
139
140 CONSTANT reference_molar_energy
141 REFINES molar_energy_constant;
142
143 CONSTANT enthalpy_of_formation_constant
144 REFINES molar_energy_constant;
145
146 CONSTANT free_energy_of_formation_constant
147 REFINES molar_energy_constant;
148
149 CONSTANT heat_of_vaporization_constant
150 REFINES molar_energy_constant;
151
152 CONSTANT Wilson_energy_constant
153 REFINES molar_energy_constant;
154
155 (* molar entropies *)
156 CONSTANT molar_entropy_constant
157 REFINES real_constant DIMENSION M*L^2/T^2/Q/TMP;
158
159 CONSTANT reference_molar_entropy
160 REFINES molar_entropy_constant;
161
162 (* other strange correlation coefficients *)
163 CONSTANT heat_capacity_constant
164 REFINES real_constant;
165
166 CONSTANT heat_capacity_a_constant
167 REFINES heat_capacity_constant DIMENSION M*L^2/T^2/Q/TMP;
168
169 CONSTANT heat_capacity_b_constant
170 REFINES heat_capacity_constant DIMENSION M*L^2/T^2/Q/TMP^2;
171
172 CONSTANT heat_capacity_c_constant
173 REFINES heat_capacity_constant DIMENSION M*L^2/T^2/Q/TMP^3;
174
175 CONSTANT heat_capacity_d_constant
176 REFINES heat_capacity_constant DIMENSION M*L^2/T^2/Q/TMP^4;
177
178 (*
179 CONSTANT
180 REFINES real_constant DIMENSION;
181 *)
182
183 (* U N I V E R S A L C O N S T A N T S
184 -------------------------------------- *)
185
186 UNIVERSAL CONSTANT gas_constant
187 (* DIMENSION M*L^2/T^2/Q/TMP *)
188 REFINES real_constant :== 1{GAS_C};
189
190 UNIVERSAL CONSTANT gravity_constant
191 (* DIMENSION L/T^2 *)
192 REFINES real_constant :== 1{EARTH_G};
193
194 UNIVERSAL CONSTANT circle_constant
195 REFINES real_constant :== 1{PI};
196
197 UNIVERSAL CONSTANT speed_of_light
198 REFINES real_constant :== 1{LIGHT_C};
199
200 UNIVERSAL CONSTANT planck_constant
201 (* DIMENSION M*L^2/T *)
202 REFINES real_constant :== 1{PLANCK_C};
203
204 UNIVERSAL CONSTANT avogadro_constant
205 REFINES real_constant :== 1{AVOGADRO_C};
206
207 UNIVERSAL CONSTANT permittivity_constant
208 (* DIMENSION E^2*T^4/M/L^3 *)
209 REFINES real_constant :== 1{EPSILON0};
210
211 UNIVERSAL CONSTANT permeability_constant
212 REFINES real_constant :== 1{MU0};
213
214 UNIVERSAL CONSTANT electron_charge
215 REFINES real_constant :== 1{eCHARGE};
216
217 UNIVERSAL CONSTANT electron_mass
218 REFINES real_constant :== 1{eMASS};
219
220 UNIVERSAL CONSTANT proton_mass
221 REFINES real_constant :== 1{pMASS};
222
223 (* B O O L E A N S *)
224
225 (* use these booleans *)
226 ATOM boolean_start_true REFINES boolean
227 DEFAULT TRUE;
228 END boolean_start_true;
229
230 ATOM boolean_start_false REFINES boolean
231 DEFAULT FALSE;
232 END boolean_start_false;
233
234 (* for backward compatibility *)
235 ATOM start_true REFINES boolean
236 DEFAULT TRUE;
237 END start_true;
238
239 ATOM start_false REFINES boolean
240 DEFAULT FALSE;
241 END start_false;
242
243
244 (* P A R A M E T E R S *)
245
246 UNIVERSAL ATOM bound_width REFINES real
247 (* not really a constant but a parameter to tell us how wide to
248 * put bounds from the current point, in relative terms.
249 * e.g. a.upper_bound := a + a.nominal * bound_width_instance;
250 *)
251 DIMENSIONLESS
252 DEFAULT 1.0e8;
253 END bound_width;
254
255 UNIVERSAL ATOM scaling_constant REFINES real
256 (* not really a constant but a parameter. which
257 * needs problem dependent information to be useful.
258 *)
259 DIMENSIONLESS
260 DEFAULT 1.0;
261 END scaling_constant;
262
263 UNIVERSAL ATOM ode_counter REFINES integer
264 DIMENSIONLESS
265 DEFAULT 1;
266 END ode_counter;
267
268 UNIVERSAL ATOM obs_counter REFINES integer
269 DIMENSIONLESS
270 DEFAULT 1;
271 END obs_counter;
272
273 ATOM real_parameter REFINES real;
274 END real_parameter;
275
276 ATOM length_parameter REFINES real_parameter
277 DIMENSION L
278 DEFAULT 1.0 {m};
279 END length_parameter;
280
281
282 (* D I M E N S I O N L E S S Q U A N T I T I E S
283 ---------------------------------------------- *)
284
285 ATOM positive_variable REFINES solver_var
286 (* one for the gams folks *)
287 DIMENSIONLESS
288 DEFAULT 1.0;
289 lower_bound := 0.0;
290 upper_bound := 1e20;
291 nominal := 1.0;
292 END positive_variable;
293
294 ATOM factor REFINES solver_var
295 DIMENSIONLESS
296 DEFAULT 1.0;
297 lower_bound := -1e50;
298 upper_bound := 1e50;
299 nominal := 1.0;
300 END factor;
301
302 ATOM variable REFINES solver_var
303 DIMENSIONLESS;
304 END variable;
305
306
307 ATOM fraction REFINES solver_var
308 DIMENSIONLESS
309 DEFAULT 0.5;
310 nominal := 1.0;
311 lower_bound := 0.0;
312 upper_bound := 1.0;
313 END fraction;
314
315 ATOM positive_factor REFINES factor;
316 lower_bound := 0.0;
317 END positive_factor;
318
319 ATOM small_factor REFINES factor;
320 lower_bound := -10.0;
321 upper_bound := 10.0;
322 END small_factor;
323
324 ATOM small_positive_factor REFINES factor;
325 lower_bound := 0.0;
326 upper_bound := 10.0;
327 END small_positive_factor;
328
329 ATOM reduced_pressure REFINES factor;
330 END reduced_pressure;
331
332 (* S U B S T I T U T I O N V A R I A B L E S
333 ------------------------------------------ *)
334 ATOM exp_sub REFINES factor ;
335 lower_bound := -1e50;
336 upper_bound := 100;
337 nominal := 1.0;
338 END exp_sub;
339
340 ATOM power_sub REFINES factor;
341 lower_bound := -25;
342 upper_bound := 25;
343 nominal := 1.0;
344 END power_sub;
345
346 (* T E M P E R A T U R E
347 --------------------- *)
348
349 ATOM temperature REFINES solver_var
350 DIMENSION TMP
351 DEFAULT 298.0{K};
352 lower_bound := 1.0e-6{K};
353 upper_bound := 10000{K};
354 nominal := 298.0{K};
355 END temperature;
356
357 ATOM inverse_temperature REFINES solver_var
358 DIMENSION 1/TMP
359 DEFAULT 0.00366099{1/K};
360 lower_bound := 0.0{1/K};
361 upper_bound := 1e50{1/K};
362 nominal := 0.00366099{1/K};
363 END inverse_temperature;
364
365 ATOM delta_temperature REFINES solver_var
366 DIMENSION TMP
367 DEFAULT 0{K};
368 lower_bound := -1000{K};
369 upper_bound := +1000{K};
370 nominal := 0{K};
371 END delta_temperature;
372
373 (* forces *)
374 ATOM force REFINES solver_var
375 DIMENSION M*L/T^2
376 DEFAULT 1.0{N};
377 lower_bound := -1e20{N};
378 upper_bound := 1e20{N};
379 nominal := 1.0{kN};
380 END force;
381
382 (* P R E S S U R E
383 --------------- *)
384
385 ATOM pressure REFINES solver_var
386 DIMENSION M/L/T^2
387 DEFAULT 1.0{atm};
388 lower_bound := 0.001{Pa};
389 upper_bound := 5000{atm};
390 nominal := 1.0{atm};
391 END pressure;
392
393 ATOM vapor_pressure REFINES pressure
394 DIMENSION M/L/T^2
395 DEFAULT 1.0{atm};
396 lower_bound := 0.001{Pa};
397 upper_bound := 5000{atm};
398 nominal := 0.5{atm};
399 END vapor_pressure;
400
401 ATOM k_constant REFINES solver_var
402 (* what IS this ? ? ? *)
403 DIMENSION T^2/L^5
404 DEFAULT 1.0 {s^2/ft^5};
405 lower_bound := 0.001 {s^2/ft^5};
406 upper_bound := 5000 {s^2/ft^5};
407 nominal := 1.0 {s^2/ft^5};
408 END k_constant;
409
410
411 (* M A S S / M O L E Q U A N T I T I E S
412 -------------------------------------- *)
413
414 ATOM molar_mass REFINES solver_var
415 DIMENSION M/Q
416 DEFAULT 100.0{g/g_mole};
417 lower_bound := 0.0{g/g_mole};
418 upper_bound := 1e9{g/g_mole};
419 nominal := 100.0{g/g_mole};
420 END molar_mass;
421
422 ATOM mass REFINES solver_var
423 DIMENSION M
424 DEFAULT 10.0{kg};
425 lower_bound := 0.0{kg};
426 upper_bound := 1e50{kg};
427 nominal := 10.0{kg};
428 END mass;
429
430 ATOM mole_scale REFINES real DIMENSION Q DEFAULT 1 {mole};
431 END mole_scale;
432
433 ATOM mole REFINES solver_var
434 DIMENSION Q
435 DEFAULT 10.0{lb_mole};
436 lower_bound := 0.0{lb_mole};
437 upper_bound := 1e50{lb_mole};
438 nominal := 10.0{lb_mole};
439 END mole;
440
441 ATOM mass_rate REFINES solver_var
442 DIMENSION M/T
443 DEFAULT 50{g/s};
444 lower_bound := 0.0{g/s};
445 upper_bound := 1e50{g/s};
446 nominal := 100.0{g/s};
447 END mass_rate;
448
449 ATOM molar_rate_scale REFINES real DIMENSION Q/T DEFAULT 1 {mole/second};
450 END molar_rate_scale;
451
452 ATOM molar_rate REFINES solver_var
453 DIMENSION Q/T
454 DEFAULT 100.0{lb_mole/hour};
455 lower_bound := 0.0{lb_mole/hour};
456 upper_bound := 1e50{lb_mole/hour};
457 nominal := 100.0{lb_mole/hour};
458 END molar_rate;
459
460 ATOM conc_rate REFINES solver_var
461 DIMENSION Q/L^3/T
462 DEFAULT 100.0{lb_mole/ft^3/hour};
463 lower_bound := 0.0{lb_mole/ft^3/hour};
464 upper_bound := 1e50{lb_mole/ft^3/hour};
465 nominal := 100.0{lb_mole/ft^3/hour};
466 END conc_rate;
467
468
469 ATOM mole_fraction REFINES fraction
470 DIMENSIONLESS
471 DEFAULT 0.5;
472 lower_bound := 0.0;
473 nominal := 0.3;
474 upper_bound := 1.0;
475 END mole_fraction;
476
477 ATOM mass_fraction REFINES fraction
478 DIMENSIONLESS
479 DEFAULT 0.5;
480 lower_bound := 0.0;
481 nominal := 0.3;
482 upper_bound := 1.0;
483 END mass_fraction;
484
485
486 (* V O L U M E Q U A N T I T I E S
487 -------------------------------- *)
488
489 ATOM molar_volume REFINES solver_var
490 DIMENSION L^3/Q
491 DEFAULT 1000.0{cm^3/g_mole};
492 lower_bound := 0.0{cm^3/g_mole};
493 upper_bound := 1e50{cm^3/g_mole};
494 nominal := 1000.0{cm^3/g_mole};
495 END molar_volume;
496
497 ATOM volume_scale REFINES real DIMENSION L^3 DEFAULT 1.0 {m^3};
498 END volume_scale;
499
500 ATOM volume REFINES solver_var
501 DIMENSION L^3
502 DEFAULT 100.0{ft^3};
503 lower_bound := 0.0{ft^3};
504 upper_bound := 1e50{ft^3};
505 nominal := 100.0{ft^3};
506 END volume;
507
508 ATOM volume_rate_scale REFINES real DIMENSION L^3/T DEFAULT 1{m^3/s};
509 END volume_rate_scale;
510
511 ATOM volume_rate REFINES solver_var
512 DIMENSION L^3/T
513 DEFAULT 100.0{gpm};
514 lower_bound := 0.0{gpm};
515 upper_bound := 1e50{gpm};
516 nominal := 100.0{gpm};
517 END volume_rate;
518
519 ATOM volume_rate_square REFINES solver_var
520 DIMENSION L^6/T^2
521 DEFAULT 100.0{ft^6/s^2};
522 lower_bound := 0.0{ft^6/s^2};
523 upper_bound := 1e50{ft^6/s^2};
524 nominal := 100 {ft^6/s^2};
525 END volume_rate_square;
526
527 ATOM volume_expansivity REFINES solver_var
528 DIMENSION 1/TMP
529 DEFAULT 0.001{1/K};
530 lower_bound := 0.0{1/K};
531 upper_bound := 1e50{1/K};
532 nominal := 0.001{1/K};
533 END volume_expansivity;
534
535
536 (* D E N S I T Y Q U A N T I T I E S
537 ---------------------------------- *)
538
539 ATOM molar_density REFINES solver_var
540 DIMENSION Q/L^3
541 DEFAULT 0.1{mole/m^3};
542 lower_bound := 0.0{mole/m^3};
543 upper_bound := 1e50{mole/m^3};
544 nominal := 0.1{mole/m^3};
545 END molar_density;
546
547 ATOM mass_density REFINES solver_var
548 DIMENSION M/L^3
549 DEFAULT 1.0{g/cm^3};
550 lower_bound := 0.0{g/cm^3};
551 upper_bound := 1e50{g/cm^3};
552 nominal := 1.0{g/cm^3};
553 END mass_density;
554
555
556 (* E N E R G Y Q U A N T I T I E S
557 ------------------------------------ *)
558
559 ATOM molar_energy REFINES solver_var
560 DIMENSION M*L^2/T^2/Q
561 DEFAULT 10000.0{BTU/lb_mole};
562 lower_bound := -1e50{BTU/lb_mole};
563 upper_bound := 1e50{BTU/lb_mole};
564 nominal := 10000.0{BTU/lb_mole};
565 END molar_energy;
566
567 ATOM energy_scale REFINES real DIMENSION M*L^2/T^2 DEFAULT 1{joule};
568 END energy_scale;
569
570 ATOM energy REFINES solver_var
571 DIMENSION M*L^2/T^2
572 DEFAULT 100000.0{BTU};
573 lower_bound := -1e50{BTU};
574 upper_bound := 1e50{BTU};
575 nominal := 100000.0{BTU};
576 END energy;
577
578 ATOM energy_rate_scale REFINES real DIMENSION M*L^2/T^3 DEFAULT 1{watt};
579 END energy_rate_scale;
580
581 ATOM energy_rate REFINES solver_var
582 DIMENSION M*L^2/T^3
583 DEFAULT 100000.0{BTU/hour};
584 lower_bound := -1e50{BTU/hour};
585 upper_bound := 1e50{BTU/hour};
586 nominal := 100000.0{BTU/hour};
587 END energy_rate;
588
589 ATOM molar_heat_capacity REFINES solver_var
590 DIMENSION M*L^2/T^2/Q/TMP
591 DEFAULT 1.00e5{J/mole/K};
592 lower_bound := 0.0{J/mole/K};
593 upper_bound := 1e60{J/mole/K};
594 nominal := 1.00e5{J/mole/K};
595 END molar_heat_capacity;
596
597 ATOM molar_energy_rate REFINES solver_var
598 DIMENSION M*L^2/T^3/Q
599 DEFAULT 0 {BTU/lb_mole/hr};
600 lower_bound := -1e50 {BTU/lb_mole/hr};
601 upper_bound := 1e50 {BTU/lb_mole/hr};
602 nominal := 10000.0 {BTU/lb_mole/hr};
603 END molar_energy_rate;
604
605 (* E N T R O P Y Q U A N T I T I E S
606 ---------------------------------- *)
607
608 ATOM molar_entropy REFINES solver_var
609 DIMENSION M*L^2/T^2/Q/TMP
610 DEFAULT 100.0{BTU/lb_mole/R};
611 lower_bound := -1e50{BTU/lb_mole/R};
612 upper_bound := 1e50{BTU/lb_mole/R};
613 nominal := 100.0{BTU/lb_mole/R};
614 END molar_entropy;
615
616 ATOM entropy REFINES solver_var
617 DIMENSION M*L^2/T^2/TMP
618 DEFAULT 1000.0{BTU/R};
619 lower_bound := -1e50{BTU/R};
620 upper_bound := 1e50{BTU/R};
621 nominal := 1000.0{BTU/R};
622 END entropy;
623
624 ATOM entropy_rate REFINES solver_var
625 DIMENSION M*L^2/T^3/TMP
626 DEFAULT 1000.0{BTU/hour/R};
627 lower_bound := -1e50{BTU/hour/R};
628 upper_bound := 1e50{BTU/hour/R};
629 nominal := 1000.0{BTU/hour/R};
630 END entropy_rate;
631
632
633
634 (* E Q U I L I B R I U M Q U A N T I T I E S
635 ------------------------------------------ *)
636
637 ATOM partition_coefficient REFINES factor (* new *)
638 DEFAULT 1.0;
639 lower_bound := 1.0e-10;
640 upper_bound := 30.0;
641 nominal := 1.0;
642 END partition_coefficient;
643
644 ATOM relative_volatility REFINES partition_coefficient; (* new *)
645 END relative_volatility;
646
647
648 (* M O N E T A R Y Q U A N T I T I E S
649 ------------------------------------ *)
650
651 ATOM monetary_unit REFINES solver_var
652 DIMENSION C
653 DEFAULT 100.0{USD};
654 lower_bound := -1e50{USD};
655 upper_bound := 1e50{USD};
656 nominal := 100.0{USD};
657 END monetary_unit;
658
659 ATOM cost_per_volume REFINES solver_var
660 DIMENSION C/L^3
661 DEFAULT 1.0{USD/gallon};
662 lower_bound := 0.0{USD/gallon};
663 upper_bound := 1e50{USD/gallon};
664 nominal := 1.0{USD/gallon};
665 END cost_per_volume;
666
667 ATOM cost_per_mass REFINES solver_var
668 DIMENSION C/M
669 DEFAULT 1.0{USD/lbm};
670 lower_bound := 0.0{USD/lbm};
671 upper_bound := 1e50{USD/lbm};
672 nominal := 1.0{USD/lbm};
673 END cost_per_mass;
674
675 ATOM cost_per_mole REFINES solver_var
676 DIMENSION C/Q
677 DEFAULT 1.0{USD/lb_mole};
678 lower_bound := 0.0{USD/lb_mole};
679 upper_bound := 1e50{USD/lb_mole};
680 nominal := 1.0{USD/lb_mole};
681 END cost_per_mole;
682
683 ATOM cost_per_time REFINES solver_var
684 DIMENSION C/T
685 DEFAULT 1.0{USD/min};
686 lower_bound := 0.0{USD/min};
687 upper_bound := 1e50{USD/min};
688 nominal := 1.0{USD/min};
689 END cost_per_time;
690
691 ATOM cost_per_energy REFINES solver_var
692 DIMENSION C*T^2/M/L^2
693 DEFAULT 1.0{USD/BTU};
694 lower_bound := 0.0{USD/BTU};
695 upper_bound := 1e50{USD/BTU};
696 nominal := 1.0{USD/BTU};
697 END cost_per_energy;
698
699
700 (* S U R V E Y I N G Q U A N T I T I E S
701 --------------------------------------- *)
702
703 ATOM distance REFINES solver_var
704 DIMENSION L
705 DEFAULT 10.0{ft};
706 lower_bound := 0.0{ft};
707 upper_bound := 1e50{ft};
708 nominal := 10.0{ft};
709 END distance;
710
711 ATOM area REFINES solver_var
712 DIMENSION L^2
713 DEFAULT 10.0{ft^2};
714 lower_bound := 0.0{ft^2};
715 upper_bound := 1e50{ft^2};
716 nominal := 10.0{ft^2};
717 END area;
718
719 ATOM inverse_area REFINES solver_var
720 DIMENSION L^-2
721 DEFAULT 1{1/ft^2};
722 lower_bound := 0.0{1/ft^2};
723 upper_bound := 1e50{1/ft^2};
724 nominal := 1.0{1/ft^2};
725 END inverse_area;
726
727 ATOM angle REFINES solver_var
728 DIMENSION P
729 DEFAULT 1 {rad};
730 lower_bound := -1e50 {rad};
731 upper_bound := 1e50 {rad};
732 nominal := 1 {rad};
733 END angle;
734
735 ATOM solid_angle REFINES solver_var
736 DIMENSION S
737 DEFAULT 1 {srad};
738 lower_bound := -1e50 {srad};
739 upper_bound := 1e50 {srad};
740 nominal := 1 {srad};
741 END solid_angle;
742
743 (* M O T I O N Q U A N T I T I E S
744 --------------------------------- *)
745
746 ATOM time REFINES solver_var
747 DIMENSION T
748 DEFAULT 60.0{s};
749 lower_bound := -1e50{s};
750 upper_bound := 1e50{s};
751 nominal := 60.0{s};
752 END time;
753
754 ATOM speed REFINES solver_var
755 DIMENSION L/T
756 DEFAULT 3.0{ft/s};
757 lower_bound := -1e50{m/s};
758 upper_bound := 1e50{m/s};
759 nominal := 1.0{m/s};
760 END speed;
761
762 ATOM acceleration REFINES solver_var
763 DIMENSION L/T^2
764 DEFAULT 9.8{m/s^2};
765 lower_bound := -1e50{m/s^2};
766 upper_bound := 1e50{m/s^2};
767 nominal := 9.8{m/s^2};
768 END acceleration;
769
770 ATOM frequency REFINES solver_var
771 DIMENSION 1/T
772 DEFAULT 60.0{1/s};
773 lower_bound := 0.0{1/s};
774 upper_bound := 1e50{1/s};
775 nominal := 60.0{1/s};
776 END frequency;
777
778
779 (* T R A N S P O R T Q U A N T I T I E S
780 --------------------------------------- *)
781
782 ATOM viscosity REFINES solver_var
783 DIMENSION M/L/T
784 DEFAULT 1.0{cP};
785 lower_bound := 0.0{cP};
786 upper_bound := 1e50{cP};
787 nominal := 1.0{cP};
788 END viscosity;
789
790 ATOM thermal_conductivity REFINES solver_var
791 DIMENSION M*L/T^3/TMP
792 DEFAULT 0.1{BTU/hour/ft/R};
793 lower_bound := 0.0{BTU/hour/ft/R};
794 upper_bound := 1e50{BTU/hour/ft/R};
795 nominal := 0.1{BTU/hour/ft/R};
796 END thermal_conductivity;
797
798 ATOM diffusivity REFINES solver_var
799 DIMENSION L^2/T
800 DEFAULT 1.0{cm^2/s};
801 lower_bound := 0.0{cm^2/s};
802 upper_bound := 1e50{cm^2/s};
803 nominal := 1.0{cm^2/s};
804 END diffusivity;
805
806
807 (* E L E C T R O - M A G N E T I C Q U A N T I T I E S
808 ---------------------------------------------------- *)
809
810 ATOM voltage REFINES solver_var
811 DIMENSION M*L^2/(T^3*E)
812 DEFAULT 10.0{volt};
813 lower_bound := -1e50{volt};
814 upper_bound := 1e50{volt};
815 nominal := 10.0{volt};
816 END voltage;
817
818 ATOM resistance REFINES solver_var
819 DIMENSION M*L^2/(E^2*T^3)
820 DEFAULT 10.0{ohm};
821 lower_bound := -1e50{ohm};
822 upper_bound := 1e50{ohm};
823 nominal := 10.0{ohm};
824 END resistance;
825
826 ATOM current REFINES solver_var
827 DIMENSION E
828 DEFAULT 1.0{amp};
829 lower_bound := -1e50{amp};
830 upper_bound := 1e50{amp};
831 nominal := 1.0{amp};
832 END current;
833
834 ATOM capacitance REFINES solver_var
835 DIMENSION E^2*T^4/M/L^2
836 DEFAULT 1.0{coulomb/volt};
837 lower_bound := -1e50{coulomb/volt};
838 upper_bound := 1e50{coulomb/volt};
839 nominal := 1.0{coulomb/volt};
840 END capacitance;
841
842 ATOM inductance REFINES solver_var
843 DIMENSION M*L^2/E^2/T^2
844 DEFAULT 1.0{volt*s/amp};
845 lower_bound := -1e50{volt*s/amp};
846 upper_bound := 1e50{volt*s/amp};
847 nominal := 1.0{volt*s/amp};
848 END inductance;
849
850 ATOM magnetic_field REFINES solver_var
851 DIMENSION E/L
852 DEFAULT 1.0{amp/m};
853 lower_bound := -1e50{amp/m};
854 upper_bound := 1e50{amp/m};
855 nominal := 1.0{amp/m};
856 END magnetic_field;
857
858 ATOM electric_field REFINES solver_var
859 DIMENSION M*L/E/T^3
860 DEFAULT 1.0{volt/m};
861 lower_bound := -1e50{volt/m};
862 upper_bound := 1e50{volt/m};
863 nominal := 1.0{volt/m};
864 END electric_field;
865
866 (* D I F F E R E N T I A L Q U A N T I T I E S
867 ------------------------------------------- *)
868
869 ATOM delta_distance REFINES solver_var
870 DIMENSION L
871 DEFAULT 10.0{ft};
872 lower_bound := -1e50{ft};
873 upper_bound := 1e50{ft};
874 nominal := 10.0{ft};
875 END delta_distance;
876
877 ATOM temperature_rate REFINES solver_var
878 DIMENSION TMP/T
879 DEFAULT 0{K/s};
880 lower_bound := -100{K/s};
881 upper_bound := 100{K/s};
882 nominal := 298.0{K/s};
883 END temperature_rate;
884
885 ATOM delta_mass REFINES solver_var
886 DIMENSION M
887 DEFAULT 0.0{kg};
888 lower_bound := -1e50{kg};
889 upper_bound := 1e50{kg};
890 nominal := 10.0{kg};
891 END delta_mass;
892
893 ATOM delta_mole REFINES solver_var
894 DIMENSION Q
895 DEFAULT 0{lb_mole};
896 lower_bound := -1e50{lb_mole};
897 upper_bound := 1e50{lb_mole};
898 nominal := 10.0{lb_mole};
899 END delta_mole;
900
901 ATOM delta_mass_rate REFINES solver_var
902 DIMENSION M/T
903 DEFAULT 0{g/s};
904 lower_bound := -1e50{g/s};
905 upper_bound := 1e50{g/s};
906 nominal := 100.0{g/s};
907 END delta_mass_rate;
908
909 ATOM delta_molar_rate REFINES solver_var
910 DIMENSION Q/T
911 DEFAULT 0.0{lb_mole/hour};
912 lower_bound := -1e50{lb_mole/hour};
913 upper_bound := 1e50{lb_mole/hour};
914 nominal := 100.0{lb_mole/hour};
915 END delta_molar_rate;
916
917 ATOM delta_volume_rate REFINES solver_var
918 DIMENSION L^3/T
919 DEFAULT 0.0{gpm};
920 lower_bound := -1e50{gpm};
921 upper_bound := 1e50{gpm};
922 nominal := 100.0{gpm};
923 END delta_volume_rate;
924
925 ATOM delta_energy_rate REFINES solver_var
926 DIMENSION M*L^2/T^3
927 DEFAULT 0.0{BTU/hour};
928 lower_bound := -1e50{BTU/hour};
929 upper_bound := 1e50{BTU/hour};
930 nominal := 100000.0{BTU/hour};
931 END delta_energy_rate;
932
933 ATOM delta_molar_energy_rate REFINES solver_var
934 DIMENSION M*L^2/T^3/Q
935 DEFAULT 0 {BTU/lb_mole/hr};
936 lower_bound := -1e50 {BTU/lb_mole/hr};
937 upper_bound := 1e50 {BTU/lb_mole/hr};
938 nominal := 10000.0 {BTU/lb_mole/hr};
939 END delta_molar_energy_rate;
940
941 ATOM delta_entropy_rate REFINES solver_var
942 DIMENSION M*L^2/T^3/TMP
943 DEFAULT 0.0{BTU/hour/R};
944 lower_bound := -1e50{BTU/hour/R};
945 upper_bound := 1e50{BTU/hour/R};
946 nominal := 1000.0{BTU/hour/R};
947 END delta_entropy_rate;
948
949 (* C O N T R O L L E R Q U A N T I T I E S
950 ---------------------------------------- *)
951
952 ATOM mass_sec REFINES solver_var
953 DIMENSION M*T
954 DEFAULT 0.0{kg*s};
955 lower_bound := -1e50{kg*s};
956 upper_bound := 1e50{kg*s};
957 nominal := 10.0{kg*s};
958 END mass_sec;
959
960 ATOM mole_sec REFINES solver_var
961 DIMENSION Q*T
962 DEFAULT 0.0{lb_mole*s};
963 lower_bound := -1e50{lb_mole*s};
964 upper_bound := 1e50{lb_mole*s};
965 nominal := 10.0{lb_mole*s};
966 END mole_sec;

john.pye@anu.edu.au
ViewVC Help
Powered by ViewVC 1.1.22