next up previous
Naprej: Prikaz rezultatov Gor: Laboratorijske vaje pri predmetu Nazaj: Iteracija

Program

Program prebere podatke iz tekstne datoteke, ki jo uporabnik napise z editorjem v predpisani obliki. Pravilnost podatkov se sproti kontrolira in izpisujejo se napake ter opozorila pri interpretiranju vhodne datoteke. Vhodna datoteka naj ima koncnico .TXI. Po preracunu pa naj program izpise na zaslon rezultate obdelave, kot so: vhodni podatki, pretoki, hitrosti, Reynoldsovo stevilo, stevilo iteracij ... Kasneje pa bo potrebno rezultate zapisati se v binarni obliki z direktno datoteko za posredovanje podatkov postprocesorju. Ne trudite se izpisovati rezultate se v tekstovno datoteko, saj to omogoca ze sam operacijski sistem z ukazom kot je

echo podatki.txi | flow > rezultati.txt

Vhodna datoteka naj ima sledec format:

constants 1000 1.01e-6  # gostota, viskoznost
junctions 6             # i, x, y, z, Qi
1   0.0 100.0 1.0 -0.040
2  50.0 100.0 2.0 -0.020
3   0.0  50.0 0.5  0.025
4  50.0  50.0 1.5  0.000
5   0.0   0.0 1.7 -0.005
6   0.0  50.0 0.0  0.040
pipes 7 # i, za\v{c}, kon, l[m], d[m], k[m]
1 1 2 50.0 0.1  0.0005
2 2 4 50.0 0.1  0.0005
3 1 3 50.0 0.1  0.0005
4 4 3 50.0 0.1  0.0005
5 5 3 50.0 0.1  0.0005
6 4 6 50.0 0.1  0.0005
7 5 6 50.0 0.1  0.0005
loops 2 # \v{s}t. cevi {smer cevi ...}
4  1  2  4 -3
4 -4  6 -7  5
solve

Zacetne delovne tocke se izracuna kot povprecje vseh pretokov kar je osnova za zacetne priblizke pretokov.

Celoten zacetni sistem enacb z upostevanjem enacb (10 in 11) ima obliko

Konstanta R = 85927 ima to vrednost le pri prvem izracunu sistema enacb, kasneje pa se matrika pri iteraciji spreminja (le zancni del enacb). Sistem enacb je linearen. Neznane vrednosti so , ki jih izracunamo z eno od metod za resevanje sistema linearnih enacb.

Rezultat sedmih iteracij so pretoki delovne tocke z relativno natancnostjo 1%





Leon Kos
Mon Oct 9 07:22:25 GMT+0100 1995