// // all quantities in ISO units m,s,.. // dimension { length <0.7,0.300,0.24> resolution <50,6,44> // grid lines in x direction x 1 0.0000000 2 0.0253562 3 0.0498104 4 0.0733925 5 0.0961298 6 0.1180474 7 0.1391685 8 0.1595147 9 0.1791066 10 0.1979636 11 0.2161046 12 0.2335477 13 0.2503107 14 0.2664112 15 0.2818661 16 0.2966927 17 0.3109079 18 0.3245287 19 0.3375718 20 0.3500543 21 0.3619928 22 0.3734042 23 0.3843051 24 0.3947124 25 0.4046428 26 0.4141196 27 0.4231771 28 0.4318578 29 0.4402088 30 0.4482805 31 0.4561217 32 0.4637720 33 0.4712619 34 0.4786150 35 0.4858482 36 0.4929734 37 0.5000000 38 0.5069093 39 0.5138026 40 0.5207614 41 0.5278965 42 0.5353575 43 0.5433491 44 0.5521566 45 0.5621568 46 0.5738591 47 0.5879867 48 0.6056022 49 0.6283084 50 0.6585519 51 0.7000000 // grid lines in z direction z 1 0.0000000 2 0.0096385 3 0.0185197 4 0.0267125 5 0.0342859 6 0.0413056 7 0.0478335 8 0.0539260 9 0.0596343 10 0.0650045 11 0.0700783 12 0.0748927 13 0.0794812 14 0.0838739 15 0.0880981 16 0.0921789 17 0.0961391 18 0.1000000 19 0.1040837 20 0.1081018 21 0.1120728 22 0.1160146 23 0.1199442 24 0.1238782 25 0.1278336 26 0.1318272 27 0.1358765 28 0.1400000 29 0.1439182 30 0.1479427 31 0.1520938 32 0.1563932 33 0.1608641 34 0.1655319 35 0.1704241 36 0.1755710 37 0.1810056 38 0.1867645 39 0.1928877 40 0.1994192 41 0.2064069 42 0.2139024 43 0.2219612 44 0.2306408 45 0.2400000 } parameter { prstep 100 // #time steps between output deltmax 1.0 // max \Delta t Tfin 0.03 // stop if t>= Tfin eps 1e-3 // stop SOR if residual lower than eps itermax 100 // max. #iterations in SOR omega 1.7 // SOR relaxation parameter alpha 0.0 // blending term for 1st order upwind /central differences tfconv 0.1 // max. CFL number nu 0.0000152 // viscisity of air gx 0.0 // forcing terms gy 0.0 gz 0.0 alphaTC 1.0 // as \alpha for transport of species nuC 1, 0.000007 // diffusion constant of species 1 in air periodboundy } // left wall of box with Dirichler BC for concentration box { south cheminit 1.0 } // lower wall of box box { coords <0,1,1>,<36,6,17> } // upper wall box { coords <0,1,28>,<36,6,44> } // right box{ north slip } // inflow from top box{ top inflow<0.0,0.0,-0.01> cheminit 0.0 } // outflow, natural BC box{ bottom inout 2,0 }