Updated 2016-05-15 17:48:56 by gold

One Dimension Heat Flow Model and eTCL Slot Calculator Demo Example, numerical analysis edit

This page is under development. Comments are welcome, but please load any comments in the comments section at the bottom of the page. Please include your wiki MONIKER in your comment with the same courtesy that I will give you. Its very hard to reply intelligibly without some background of the correspondent. Thanks,gold

gold Here is some eTCL code for One Dimension Heat Flow Model.

Testcases Section

In planning any software, it is advisable to gather a number of testcases to check the results of the program. The math for the testcases can be checked by pasting statements in the TCL console. Aside from the TCL calculator display, when one presses the report button on the calculator, one will have console show access to the capacity functions (subroutines).

Testcase 1

table 1printed in tcl wiki format
quantity value comment, if any
testcase number:1
150.0 :initial temperature t1 degrees
640.0 :final temperature t2 degrees
0.3048 :length meters
10.0 :boundary length meters
0.3 :temperature constant
0.3 :temperature coefficient
481.98346456692906 :heat flux
164.9352 :temperature degrees


Testcase 2

table 2printed in tcl wiki format
quantity value comment, if any
testcase number:2
100.0 :initial temperature t1 degrees
600.0 :final temperature t2 degrees
1.0 :length meters
10.0 :boundary length meters
1.0 :temperature constant
1.0 :temperature coefficient
499.0 :heat flux
150.0 :temperature degrees

Testcase 3

table 3printed in tcl wiki format
quantity value comment, if any
testcase number:3
120.0 :initial temperature t1 degrees
530.0 :final temperature t2 degrees
5.0 :length meters
10.0 :boundary length meters
1.0 :temperature constant
1.0 :temperature coefficient
81.0 :heat flux
325.0 :temperature degrees


Screenshots Section

figure 1.


References:

  • Joran

Appendix Code edit

appendix TCL programs and scripts

        # pretty print from autoindent and ased editor
        # One Dimensional Heat Flow calculator
        # written on Windows XP on eTCL
        # working under TCL version 8.5.6 and eTCL 1.0.1
        # gold on TCL WIKI, 10may2016
        package require Tk
        namespace path {::tcl::mathop ::tcl::mathfunc}
        frame .frame -relief flat -bg aquamarine4
        pack .frame -side top -fill y -anchor center
        set names {{} {initial temperature t1 degrees:} }
        lappend names {final temperature t2 degrees :}
        lappend names {length meters: }
        lappend names {boundary length meters:}
        lappend names {temperature constant:}
        lappend names {temperature coeffient: }
        lappend names {answers:heat flux: }
        lappend names {temperature :}
        foreach i {1 2 3 4 5 6 7 8} {
            label .frame.label$i -text [lindex $names $i] -anchor e
            entry .frame.entry$i -width 35 -textvariable side$i
            grid .frame.label$i .frame.entry$i -sticky ew -pady 2 -padx 1 }
        proc about {} {
            set msg "Calculator for One Dimensional Heat Flow 
            from TCL WIKI,
            written on eTCL "
            tk_messageBox -title "About" -message $msg } 
       proc pi {} {expr acos(-1)}       
       proc calculate {     } {
            global answer2
            global side1 side2 side3 side4 side5
            global side6 side7 side8
            global testcase_number
            incr testcase_number 
            set side1 [* $side1 1. ]
            set side2 [* $side2 1. ]
            set side3 [* $side3 1. ]
            set side4 [* $side4 1. ]
            set side5 [* $side5 1. ]
            set side6 [* $side6 1. ]
            set side7 [* $side7 1. ]
            set side8 [* $side8 1. ]
            set temperature1 $side1  
            set temperature2 $side2 
            set length1 $side3
            set lengthx $side4
            set constant1 $side5
            set coefficient2 $side6
            set heat_flow [- [* $coefficient2 [/ [- $temperature2 $temperature1 ] $length1] ] $constant1 ]
            #set heat_flow [* $coefficient2 [/ [- $temperature2 $temperature1 ] $length1 ] ]
            set side7 $heat_flow 
            #set result1 [+ $temperature1 [* [/ [- $temperature2 $temperature1 ] $length1] $lengthx ]]
            set result1 [+ $temperature1 [* [/ [- $temperature2 $temperature1 ] $lengthx] $length1 ]]
            set side8 $result1       }               
        proc fillup {aa bb cc dd ee ff gg hh} {
            .frame.entry1 insert 0 "$aa"
            .frame.entry2 insert 0 "$bb"
            .frame.entry3 insert 0 "$cc"
            .frame.entry4 insert 0 "$dd"
            .frame.entry5 insert 0 "$ee"
            .frame.entry6 insert 0 "$ff" 
            .frame.entry7 insert 0 "$gg"
            .frame.entry8 insert 0 "$hh" 
             }
        proc clearx {} {
            foreach i {1 2 3 4 5 6 7 8 } {
                .frame.entry$i delete 0 end } }
        proc reportx {} {
            global side1 side2 side3 side4 side5
            global side6 side7 side8 
            global testcase_number
            console show;
            puts "%|table $testcase_number|printed in| tcl wiki format|% "
            puts "&| quantity| value| comment, if any|& "
            puts "&| testcase number:|$testcase_number | |&"
            puts "&| $side1 :|initial temperature t1 degrees|   |&"
            puts "&| $side2 :|final temperature t2 degrees| |& "
            puts "&| $side3 :|length meters| |& "
            puts "&| $side4 :|boundary  length meters |      |&"
            puts "&| $side5 :|temperature constant|     |&"
            puts "&| $side6 :|temperature coefficient |  |&"
            puts "&| $side7 :|heat flux|  |&"
            puts "&| $side8 :|temperature degrees|  |&" 
            }
        frame .buttons -bg aquamarine4
        ::ttk::button .calculator -text "Solve" -command { calculate   }
        ::ttk::button .test2 -text "Testcase1" -command {clearx;fillup 150. 640.0 .3048 10.0  0.3 0.3 472. 164.09}
        ::ttk::button .test3 -text "Testcase2" -command {clearx;fillup 100. 600.0 1.0 10.0   1.0  1.0 450.5 165.09 }
        ::ttk::button .test4 -text "Testcase3" -command {clearx;fillup 120. 530.0 5.0 10.0   1.0 1.0  420.5 170.09 }
        ::ttk::button .clearallx -text clear -command {clearx }
        ::ttk::button .about -text about -command about
        ::ttk::button .cons -text report -command { reportx }
        ::ttk::button .exit -text exit -command {exit}
        pack .calculator  -in .buttons -side top -padx 10 -pady 5
        pack  .clearallx .cons .about .exit .test4 .test3 .test2   -side bottom -in .buttons
        grid .frame .buttons -sticky ns -pady {0 10}
        . configure -background aquamarine4 -highlightcolor brown -relief raised -border 30
        wm title . "One Dimensional Heat Flow Calculator"

Pushbutton Operation

For the push buttons, the recommended procedure is push testcase and fill frame, change first three entries etc, push solve, and then push report. Report allows copy and paste from console.

For testcases in a computer session, the eTCL calculator increments a new testcase number internally, eg. TC(1), TC(2) , TC(3) , TC(N). The testcase number is internal to the calculator and will not be printed until the report button is pushed for the current result numbers. The current result numbers will be cleared on the next solve button. The command { calculate; reportx } or { calculate ; reportx; clearx } can be added or changed to report automatically. Another wrinkle would be to print out the current text, delimiters, and numbers in a TCL wiki style table as
  puts " %| testcase $testcase_number | value| units |comment |%"
  puts " &| volume| $volume| cubic meters |based on length $side1 and width $side2   |&"  

gold This page is copyrighted under the TCL/TK license terms, this license.

Comments Section edit

Please place any comments here, Thanks.