00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028 #include "kunzel.h"
00029 #include "globalt.h"
00030
00031
00032 kunmat::kunmat()
00033 {
00034
00035 rho_m = 0.0;
00036
00037 rho_w = 1000.0;
00038 }
00039
00040 kunmat::~kunmat()
00041 {
00042
00043 }
00044
00045
00046
00047
00048
00049
00050 void kunmat::read(XFILE *in)
00051 {
00052 long i,nr,nc;
00053
00054 for (i = 2; i<20; i++){
00055 xfscanf(in, "%d",&MatChar[i]);
00056 }
00057
00058 for (i = 2;i <= 19; i++){
00059
00060 switch (MatChar[i]){
00061 case 0:{
00062
00063 break;
00064 }
00065 case 1:{
00066
00067
00068 data[i]=new gfunct ();
00069
00070
00071 data[i]->tfunc=stat;
00072
00073 xfscanf(in, "%lf",&data[i]->f);
00074
00075 break;
00076 }
00077 case 2:{
00078
00079
00080
00081
00082 xfscanf(in, "%ld %ld",&nr,&nc);
00083
00084
00085 data[i]=new gfunct (tab,nr);
00086
00087 data[i]->tabf[0].itype = piecewiselin;
00088
00089 data[i]->tabf[0].readval (in);
00090
00091
00092 xfscanf(in, "%d",&nr);
00093
00094 break;
00095 }
00096 case 3:{
00097
00098 break;
00099 }
00100 case 30:{
00101 if( MatChar[i-1] ==1)
00102 {
00103 xfscanf(in, "%lf %lf",&MatFunce[i][0],&MatFunce[i][1]);
00104 }
00105 else
00106 {
00107 xfscanf(in, "%lf %lf",&MatFunce[i][0],&MatFunce[i][1]);
00108 }
00109 break;
00110 }
00111 case 31:{
00112 if( MatChar[i-1] ==1)
00113 {
00114 if (i == 6)
00115 {
00116 xfscanf(in, "%lf %lf %lf",&MatFunce[i][0],&MatFunce[i][1], &MatFunce[i][2]);
00117 }
00118 else
00119 {
00120 xfscanf(in, "%lf %lf",&MatFunce[i][0],&MatFunce[i][1]);
00121 }
00122 }
00123 else
00124 {
00125 if (i == 6)
00126 {
00127 xfscanf(in, "%lf %lf %lf",&MatFunce[i][0],&MatFunce[i][1], &MatFunce[i][2]);
00128 }
00129 else
00130 {
00131 xfscanf(in, "%lf %lf",&MatFunce[i][0],&MatFunce[i][1]);
00132 }
00133 }
00134
00135 break;
00136 }
00137 case 32:{
00138 if( MatChar[i-1] ==1)
00139 {
00140 xfscanf(in, "%lf %lf",&MatFunce[i][0],&MatFunce[i][1]);
00141 }
00142 else
00143 {
00144 xfscanf(in, "%lf %lf",&MatFunce[i][0],&MatFunce[i][1]);
00145 }
00146 break;
00147 }
00148 default:{
00149 print_err("unknown definition of Material KUNZEL is required",__FILE__,__LINE__,__func__);
00150 }
00151 }
00152 }
00153
00154 CorD(2,kd,0,rho_m,a2,a3);
00155 }
00156
00157
00158
00159
00160
00161
00162
00163 void kunmat::print(FILE *out)
00164 {
00165 long i,j;
00166
00167 fprintf (out,"\n");
00168 for (i = 2; i<19; i++)
00169 {
00170 fprintf(out, " %ld ",MatChar[i]);
00171 }
00172 fprintf(out, " %ld ",MatChar[19]);
00173 for (i = 2;i <= 19; i++)
00174 {
00175
00176 switch (MatChar[i]){
00177 case 0:{
00178
00179 break;
00180 }
00181 case 1:{
00182
00183 fprintf(out, "\n %e ",data[i]->f);
00184 break;
00185 }
00186 case 2:{
00187
00188
00189 fprintf(out, "\n %ld %d ",data[i]->tabf[0].asize,2);
00190
00191 for (j=0;j<data[i]->tabf[0].asize;j++){
00192 fprintf(out, "\n %e %e ",data[i]->tabf[0].x[j],data[i]->tabf[0].y[j]);
00193 }
00194
00195 fprintf(out, "\n %ld ",data[i]->tabf[0].asize);
00196 break;
00197 }
00198 case 3:{
00199
00200 break;
00201 }
00202 case 30:{
00203 if( MatChar[i-1] ==1)
00204 {
00205 fprintf(out, "\n %e %e ",MatFunce[i][0],MatFunce[i][1]);
00206 }
00207 else
00208 {
00209 fprintf(out, "\n %e %e ",MatFunce[i][0],MatFunce[i][1]);
00210 }
00211 break;
00212 }
00213 case 31:{
00214 if( MatChar[i-1] ==1)
00215 {
00216 if (i == 6)
00217 {
00218 fprintf(out, "\n %e %e %e ",MatFunce[i][0],MatFunce[i][1], MatFunce[i][2]);
00219 }
00220 else
00221 {
00222 fprintf(out, "\n %e %e ",MatFunce[i][0],MatFunce[i][1]);
00223 }
00224 }
00225 else
00226 {
00227 if (i == 6)
00228 {
00229 fprintf(out, "\n %e %e %e",MatFunce[i][0],MatFunce[i][1], MatFunce[i][2]);
00230 }
00231 else
00232 {
00233 fprintf(out, "\n %e %e ",MatFunce[i][0],MatFunce[i][1]);
00234 }
00235 }
00236
00237 break;
00238 }
00239 case 32:{
00240 if( MatChar[i-1] ==1)
00241 {
00242 fprintf(out, "\n %e %e ",MatFunce[i][0],MatFunce[i][1]);
00243 }
00244 else
00245 {
00246 fprintf(out, "\n %e %e ",MatFunce[i][0],MatFunce[i][1]);
00247 }
00248 break;
00249 }
00250 default:{
00251 print_err("unknown definition of Material KUNZEL is required",__FILE__,__LINE__,__func__);
00252 fprintf (out,"\n");
00253 }
00254 }
00255 }
00256 fprintf (out,"\n");
00257 }
00258
00259
00260
00261
00262
00263
00264
00265
00266
00267
00268
00269
00270
00271 void kunmat::matcond (matrix &d,long ri,long ci,long ipp)
00272 {
00273 long n;
00274 n = d.n;
00275
00276 switch (n){
00277 case 1:{
00278 matcond1d (d,ri,ci,ipp);
00279 break;
00280 }
00281 case 2:{
00282 matcond2d (d,ri,ci,ipp);
00283 break;
00284 }
00285 case 3:{
00286 matcond3d (d,ri,ci,ipp);
00287 break;
00288 }
00289 default:{
00290 print_err("unknown number of components of conductivity tensor is required",__FILE__,__LINE__,__func__);
00291 }
00292 }
00293 }
00294
00295
00296
00297
00298
00299
00300
00301
00302
00303
00304 void kunmat::matcond1d (matrix &d,long ri,long ci,long ipp)
00305 {
00306 double k;
00307 double phi,t;
00308 k = 0.0;
00309
00310 phi = Tm->ip[ipp].av[0];
00311 t = Tm->ip[ipp].av[1];
00312
00313 if((ri == 0) && (ci == 0))
00314 k = tokJ1(phi,t,ipp);
00315 if((ri == 0) && (ci == 1))
00316 k = tokJ2(phi,t,ipp);
00317 if((ri == 1) && (ci == 0))
00318 k = tokJ3(phi,t,ipp);
00319 if((ri == 1) && (ci == 1))
00320 k = tokJ4(phi,t,ipp);
00321
00322 d[0][0] = k;
00323 }
00324
00325
00326
00327
00328
00329
00330
00331
00332
00333 void kunmat::matcond2d (matrix &d,long ri,long ci,long ipp)
00334 {
00335 double k;
00336 double phi,t;
00337 k = 0.0;
00338
00339 phi = Tm->ip[ipp].av[0];
00340 t = Tm->ip[ipp].av[1];
00341
00342 if((ri == 0) && (ci == 0))
00343 k = tokJ1(phi,t,ipp);
00344 if((ri == 0) && (ci == 1))
00345 k = tokJ2(phi,t,ipp);
00346 if((ri == 1) && (ci == 0))
00347 k = tokJ3(phi,t,ipp);
00348 if((ri == 1) && (ci == 1))
00349 k = tokJ4(phi,t,ipp);
00350
00351 fillm(0.0,d);
00352
00353
00354
00355 d[0][0] = k; d[0][1] = 0.0;
00356 d[1][0] = 0.0; d[1][1] = k;
00357 }
00358
00359
00360
00361
00362
00363
00364
00365
00366
00367
00368 void kunmat::matcond3d (matrix &d,long ri,long ci,long ipp)
00369 {
00370 double k;
00371 double phi,t;
00372 k = 0.0;
00373
00374 phi = Tm->ip[ipp].av[0];
00375 t = Tm->ip[ipp].av[1];
00376
00377 if((ri == 0) && (ci == 0))
00378 k = tokJ1(phi,t,ipp);
00379 if((ri == 0) && (ci == 1))
00380 k = tokJ2(phi,t,ipp);
00381 if((ri == 1) && (ci == 0))
00382 k = tokJ3(phi,t,ipp);
00383 if((ri == 1) && (ci == 1))
00384 k = tokJ4(phi,t,ipp);
00385
00386 fillm(0.0,d);
00387
00388 d[0][0]=k; d[0][1]=0.0; d[0][2]=0.0;
00389 d[1][0]=0.0; d[1][1]=k; d[1][2]=0.0;
00390 d[2][0]=0.0; d[2][1]=0.0; d[2][2]=k;
00391 }
00392
00393
00394
00395
00396
00397
00398
00399
00400
00401
00402 void kunmat::matcap (double &c,long ri,long ci,long ipp)
00403 {
00404 double phi,t;
00405 c=0.0;
00406
00407 if((ri == 0) && (ci == 0))
00408 c = Tm->ip[ipp].eqother[1];
00409 if((ri == 0) && (ci == 1))
00410 c = 0.0;
00411 if((ri == 1) && (ci == 0))
00412 c = 0.0;
00413 if((ri == 1) && (ci == 1)){
00414 phi = Tm->ip[ipp].av[0];
00415 t = Tm->ip[ipp].av[1];
00416 c = DerivaceHustotyEntalpiePodleTeploty(phi,t,ipp);
00417 }
00418 }
00419
00420
00421
00422
00423
00424
00425
00426
00427
00428 void kunmat::values_correction (vector &nv)
00429 {
00430 double rh,t;
00431
00432
00433 rh=nv[0];
00434
00435 t=nv[1];
00436
00437 if (rh >= 1.0)
00438 rh = 1.0;
00439 if (rh <= 0.0)
00440 rh = 0.0;
00441
00442
00443 if (t >= 350.0)
00444 t = 350.0;
00445 if (t <= 240.0)
00446 t = 240.0;
00447
00448 nv[0]=rh;
00449 nv[1]=t;
00450 }
00451
00452
00453
00454
00455
00456
00457
00458
00459
00460 double kunmat::DerivaceTlakuNasycenychParNaTeplote(double rh, double tk)
00461 {
00462
00463 return ((4042.9 * exp(23.5771-4042.9/(tk - 37.58)))/((tk-37.58)*(tk-37.58)));
00464 }
00465
00466
00467
00468
00469
00470
00471
00472
00473
00474
00475
00476
00477 double kunmat::PermeabilitaVodniPary(double rh, double tk,long ipp)
00478 {
00479 double P, Rv, Pa;
00480 double da, dp, mi, smc;
00481
00482 if (rh >1) rh = 1.0;
00483 if (rh < 0) rh = 0.0;
00484
00485
00486
00487 moist = Tm->ip[ipp].eqother[0];
00488 smc = Tm->ip[ipp].eqother[2];
00489
00490 CorD(4,kd,moist,mi,a2,a3);
00491
00492 Rv = 461.5;
00493 Pa = 101325;
00494 P = 101325;
00495
00496 da = (2.306e-5 * Pa)/(Rv * tk * P)*pow((tk/273.15),1.81);
00497
00498 if (mi==0.0)
00499 dp=0.0;
00500 else
00501 dp = da/mi;
00502
00503 return (dp);
00504 }
00505
00506
00507
00508
00509
00510
00511
00512
00513
00514
00515 double kunmat::tokJ2(double rh, double tk, long ipp)
00516 {
00517 double dp, dpvs;
00518
00519 if (rh>1.0) rh = 1.0;
00520 if (rh<0.0) rh = 0.0;
00521
00522 dpvs = DerivaceTlakuNasycenychParNaTeplote(rh,tk);
00523 dp = PermeabilitaVodniPary(rh, tk,ipp);
00524
00525 J2 = dp * rh * dpvs;
00526
00527
00528
00529 return (J2);
00530 }
00531
00532
00533
00534
00535
00536
00537
00538
00539
00540
00541
00542 double kunmat::tokJ4(double rh, double tk, long ipp)
00543 {
00544 double lv1, lambdaakt, dpvs, dp, moist,smc;
00545
00546 if (rh >1) rh = 1;
00547 if (rh < 0) rh = 0;
00548
00549
00550
00551 moist = Tm->ip[ipp].eqother[0];
00552 smc = Tm->ip[ipp].eqother[2];
00553
00554 CorD(10,kd,moist,lambdaakt,a2,a3);
00555
00556 dpvs = DerivaceTlakuNasycenychParNaTeplote(rh, tk);
00557 dp = PermeabilitaVodniPary(rh, tk,ipp);
00558
00559 lv1 = LatentHeatofEvaporationOfWater (tk);
00560
00561 J4 = lambdaakt + lv1 * dp * rh * dpvs;
00562
00563 return (J4);
00564 }
00565
00566
00567
00568
00569
00570
00571
00572
00573
00574
00575
00576
00577 double kunmat::tokJ3(double rh, double tk, long ipp)
00578 {
00579 double lv1;
00580 double dp, ps;
00581
00582 if (rh >1) rh = 1;
00583 if (rh < 0) rh = 0;
00584
00585 dp = PermeabilitaVodniPary(rh, tk,ipp);
00586 ps = TlakNasycenychVodnichParNaTeplote (rh, tk);
00587 lv1 = LatentHeatofEvaporationOfWater (tk);
00588
00589 J3 = lv1 * dp * ps;
00590
00591 return (J3);
00592 }
00593
00594
00595
00596
00597
00598
00599
00600 double kunmat::LatentHeatofEvaporationOfWater(double tk)
00601 {
00602 return (2.5008e6)*pow((273.15/tk),(0.167+tk*3.67e-4));
00603 }
00604
00605
00606
00607
00608
00609
00610
00611
00612
00613
00614
00615
00616
00617 double kunmat::TlakNasycenychVodnichParNaTeplote(double rh, double tk)
00618 {
00619 return (exp(23.5771 - 4042.9/(tk - 37.58)));
00620 }
00621
00622
00623
00624
00625
00626
00627
00628
00629 double kunmat::DerivaceHustotyEntalpiePodleTeploty(double rh, double tk, long ipp)
00630 {
00631 double rom, c, moist;
00632 if (rh >1) rh = 1;
00633 if (rh < 0) rh = 0;
00634
00635
00636 moist = Tm->ip[ipp].eqother[0];
00637
00638 CorD(2,kd,0,rom,a2,a3);
00639 CorD(9,kd,moist,c,a2,a3);
00640 rho_m = rom;
00641
00642 return (rom * c);
00643 }
00644
00645
00646
00647
00648
00649
00650
00651 double kunmat::HygroscopicMoisture(double rh, double tk)
00652 {
00653
00654 double mhmc, mhrh,a3;
00655 CorD(6,kd,rh,mhmc,mhrh,a3);
00656
00657 return (mhmc/(1-sqrt(1-mhrh)));
00658
00659 }
00660
00661
00662
00663
00664
00665
00666
00667
00668 double kunmat::DerivativeOfTheSorptionIsotherm(double rh, double tk)
00669 {
00670
00671 double mhv;
00672 if (rh >1) rh = 1;
00673 if (rh < 0) rh = 0;
00674
00675 mhv = HygroscopicMoisture (rh,tk);
00676
00677 return (1000 * mhv/(2*sqrt(1-rh)));
00678
00679 }
00680
00681
00682
00683
00684
00685
00686
00687 double kunmat::DerivativeOfTheRetentionCurve(double rh, double tk, long ipp)
00688 {
00689
00690 double smc, mhmc, mhrh;
00691
00692 smc = Tm->ip[ipp].eqother[2];
00693 CorD(6,kd,rh,mhmc,mhrh,a3);
00694
00695
00696 return (1000 * (smc - mhmc)/(1 - mhrh));
00697 }
00698
00699
00700
00701
00702
00703
00704
00705 double kunmat::tokJ1(double rh, double tk, long ipp)
00706 {
00707 double kapa1, mi, dp, dmretc, ps, drh, moist, smc;
00708
00709
00710
00711 if (rh < 0) rh = 0;
00712 if (rh >1) rh = 1;
00713
00714
00715 moist = Tm->ip[ipp].eqother[0];
00716
00717 smc = Tm->ip[ipp].eqother[2];
00718 CorD(4,kd,moist,mi,a2,a3);
00719 kapa1 = Tm->ip[ipp].eqother[3];
00720
00721
00722
00723
00724
00725
00726
00727
00728
00729
00730
00731
00732
00733
00734
00735
00736
00737
00738
00739
00740
00741
00742
00743
00744 dp = PermeabilitaVodniPary(rh, tk,ipp);
00745 ps = TlakNasycenychVodnichParNaTeplote (rh, tk);
00746
00747
00748
00749 dmretc = Tm->ip[ipp].eqother[1];
00750
00751 drh = kapa1 * dmretc;
00752
00753 J1 = drh + dp * ps;
00754
00755
00756
00757 return (J1);
00758 }
00759
00760
00761
00762
00763
00764
00765
00766
00767
00768 double kunmat::sorptionizothermDerivation(double rh, double tk)
00769 {
00770 double a, b, n;
00771
00772 CorD(6,kd,rh,a,b,n);
00773
00774
00775 if (rh >1) rh = 1;
00776 if (rh < 0) rh = 0;
00777
00778
00779 return (1000 * (a/(b*rh*n))*pow((1-(log(rh))/b),(-1-(1/n))));
00780
00781 }
00782
00783
00784
00785
00786
00787
00788
00789
00790
00791
00792
00793
00794
00795
00796
00797
00798
00799 double kunmat::transmission_transcoeff(double trc,long ri,long ci,long nn,long bc,long ipp)
00800 {
00801 long k;
00802 double new_trc,h,t;
00803 new_trc = 0.0;
00804
00805 k=Gtt->give_dof(nn,0);
00806 if (k>0) {h = Lsrst->lhs[k-1]+Lsrst->lhsi[k-1];}
00807 if (k==0) {h = 0.0;}
00808 if (k<0) {h = Tb->lc[0].pv[0-k-1].getval();}
00809 k=Gtt->give_dof(nn,1);
00810 if (k>0) {t = Lsrst->lhs[k-1]+Lsrst->lhsi[k-1];}
00811 if (k==0) {t = 0.0;}
00812 if (k<0) {t = Tb->lc[0].pv[0-k-1].getval();}
00813
00814 if((ri == 0) && (ci == 0))
00815 new_trc = get_transmission_transcoeff_hh(h,t,bc,ipp);
00816 if((ri == 0) && (ci == 1))
00817 new_trc = 0.0;
00818
00819 if((ri == 1) && (ci == 0))
00820 new_trc = 0.0;
00821 if((ri == 1) && (ci == 1))
00822 new_trc = get_transmission_transcoeff_tt(h,t,bc,ipp);
00823
00824 new_trc = new_trc*trc;
00825
00826 return (new_trc);
00827 }
00828
00829
00830
00831
00832
00833
00834
00835
00836
00837
00838
00839
00840
00841
00842 double kunmat::transmission_nodval(double nodval,double trc2,long ri,long ci,long nn,long bc,long ipp)
00843 {
00844 long k;
00845 double new_nodval,h,t;
00846 new_nodval = 0.0;
00847
00848 k=Gtt->give_dof(nn,0);
00849 if (k>0) {h = Lsrst->lhs[k-1]+Lsrst->lhsi[k-1];}
00850 if (k==0) {h = 0.0;}
00851 if (k<0) {h = Tb->lc[0].pv[0-k-1].getval();}
00852 k=Gtt->give_dof(nn,1);
00853 if (k>0) {t = Lsrst->lhs[k-1]+Lsrst->lhsi[k-1];}
00854 if (k==0) {t = 0.0;}
00855 if (k<0) {t = Tb->lc[0].pv[0-k-1].getval();}
00856
00857 if((ri == 0) && (ci == 0))
00858 new_nodval = get_transmission_nodval_hh(nodval,h,t,bc,ipp);
00859 if((ri == 0) && (ci == 1))
00860 new_nodval = 0.0;
00861
00862 if((ri == 1) && (ci == 0))
00863 new_nodval = 0.0;
00864 if((ri == 1) && (ci == 1))
00865 new_nodval = get_transmission_nodval_tt(nodval,h,t,bc,ipp);
00866
00867 return (new_nodval);
00868 }
00869
00870
00871
00872
00873
00874
00875
00876
00877
00878
00879
00880
00881
00882
00883
00884 double kunmat::transmission_flux(double nodval,double trc2,long ri,long ci,long nn,long bc,long ipp)
00885 {
00886 long k;
00887 double flux,h,t;
00888 flux = 0.0;
00889
00890 k=Gtt->give_dof(nn,0);
00891 if (k>0) {h = Lsrst->lhs[k-1]+Lsrst->lhsi[k-1];}
00892 if (k==0) {h = 0.0;}
00893 if (k<0) {h = Tb->lc[0].pv[0-k-1].getval();}
00894 k=Gtt->give_dof(nn,1);
00895 if (k>0) {t = Lsrst->lhs[k-1]+Lsrst->lhsi[k-1];}
00896 if (k==0) {t = 0.0;}
00897 if (k<0) {t = Tb->lc[0].pv[0-k-1].getval();}
00898
00899 if((ri == 0) && (ci == 0))
00900 flux = get_transmission_flux_hh(nodval,h,t,bc,ipp);
00901 if((ri == 0) && (ci == 1))
00902 flux = 0.0;
00903
00904 if((ri == 1) && (ci == 0))
00905 flux = 0.0;
00906 if((ri == 1) && (ci == 1))
00907 flux = get_transmission_flux_tt(nodval,h,t,bc,ipp);
00908
00909 return (flux);
00910 }
00911
00912
00913
00914
00915
00916
00917
00918
00919
00920 double kunmat::get_transmission_transcoeff_hh(double rh,double t,long bc,long ipp)
00921 {
00922 double trc,pgws;
00923
00924 switch (bc){
00925 case 30:{
00926 pgws = TlakNasycenychVodnichParNaTeplote(rh,t);
00927 trc = pgws;
00928 break;
00929 }
00930 case 31:{
00931 pgws = TlakNasycenychVodnichParNaTeplote(rh,t);
00932 trc = pgws;
00933 break;
00934 }
00935 case 32:{
00936 trc = 0.0;
00937 break;
00938 }
00939 default:{
00940 print_err("no real boundary condition is prescribed",__FILE__,__LINE__,__func__);
00941 exit(0);
00942 }
00943 }
00944 return(trc);
00945 }
00946
00947
00948
00949
00950
00951
00952
00953
00954
00955 double kunmat::get_transmission_nodval_hh(double bv,double rh,double t,long bc,long ipp)
00956 {
00957 double nodval,pgws, pgwspred;
00958
00959 switch (bc){
00960 case 30:{
00961
00962 pgws = TlakNasycenychVodnichParNaTeplote(rh,t);
00963 pgwspred = TlakNasycenychVodnichParNaTeplote(rh,t);
00964
00965 nodval = bv*pgwspred;
00966
00967 break;
00968 }
00969 case 31:{
00970 nodval = bv;
00971 break;
00972 }
00973 case 32:{
00974 pgws = TlakNasycenychVodnichParNaTeplote(rh,t);
00975
00976 nodval = pgws*rh;
00977 bv = pgws*bv;
00978 nodval = bv - nodval;
00979
00980 break;
00981 }
00982 default:{
00983 print_err("no real boundary condition is prescribed",__FILE__,__LINE__,__func__);
00984 exit(0);
00985 }
00986 }
00987 return(nodval);
00988 }
00989
00990
00991
00992
00993
00994
00995
00996
00997
00998
00999 double kunmat::get_transmission_flux_hh(double bv,double rh,double t,long bc,long ipp)
01000 {
01001 double flux,pgws;
01002
01003 switch (bc){
01004 case 30:{
01005 pgws = TlakNasycenychVodnichParNaTeplote(rh, t);
01006
01007 flux = pgws*rh;
01008 bv = pgws*bv;
01009 flux = bv - flux;
01010
01011 break;
01012 }
01013 case 31:{
01014 pgws = TlakNasycenychVodnichParNaTeplote(rh, t);
01015 flux = pgws*rh;
01016
01017 flux = bv - flux;
01018 break;
01019 }
01020 case 32:{
01021 pgws = TlakNasycenychVodnichParNaTeplote(rh, t);
01022
01023 flux = pgws*rh;
01024 bv = pgws*bv;
01025 flux = bv - flux;
01026
01027 break;
01028 }
01029 default:{
01030 print_err("no real boundary condition is prescribed",__FILE__,__LINE__,__func__);
01031 exit(0);
01032 }
01033 }
01034 return(flux);
01035 }
01036
01037
01038
01039
01040
01041
01042
01043
01044
01045 double kunmat::get_transmission_transcoeff_tt(double h,double t,long bc,long ipp)
01046 {
01047 double trc;
01048
01049 switch (bc){
01050 case 30:{
01051 trc = 1.0;
01052 break;
01053 }
01054 default:{
01055 print_err("no real boundary condition is prescribed",__FILE__,__LINE__,__func__);
01056 exit(0);
01057 }
01058 }
01059
01060 return(trc);
01061 }
01062
01063
01064
01065
01066
01067
01068
01069
01070
01071 double kunmat::get_transmission_nodval_tt(double bv,double h,double t,long bc,long ipp)
01072 {
01073 double nodval;
01074
01075 switch (bc){
01076 case 30:{
01077 nodval = bv;
01078 break;
01079 }
01080 default:{
01081 print_err("no real boundary condition is prescribed",__FILE__,__LINE__,__func__);
01082 exit(0);
01083 }
01084 }
01085 return(nodval);
01086 }
01087
01088
01089
01090
01091
01092
01093
01094
01095
01096
01097 double kunmat::get_transmission_flux_tt(double bv,double h,double t,long bc,long ipp)
01098 {
01099 double flux;
01100
01101 switch (bc){
01102 case 30:{
01103 flux = (bv - t);
01104 break;
01105 }
01106 default:{
01107 print_err("no real boundary condition is prescribed",__FILE__,__LINE__,__func__);
01108 exit(0);
01109 }
01110 }
01111 return(flux);
01112 }
01113
01114
01115
01116
01117
01118
01119
01120
01121
01122
01123
01124 double kunmat::get_othervalue(long compother,double rh,double t, long ipp)
01125 {
01126 double other;
01127
01128 switch (compother){
01129 case 0:{
01130 other = rh;
01131 break;
01132 }
01133 case 1:{
01134 other = t;
01135 break;
01136 }
01137 case 2:{
01138 double w, a3;
01139 sorption_izoterms_giva_data(0,rh,t,w,a3, ipp);
01140 other = w;
01141 break;
01142 }
01143 case 3:{
01144 other = rh * TlakNasycenychVodnichParNaTeplote (rh, t);
01145 break;
01146 }
01147 case 4:{
01148 double w, a3;
01149 sorption_izoterms_giva_data(0,rh,t,w,a3, ipp);
01150 other = w*rho_w/rho_m;
01151 break;
01152 }
01153 case 5:{
01154 other = 0.0;
01155 break;
01156 }
01157 case 6:{
01158 other = 0.0;
01159 break;
01160 }
01161 case 7:{
01162 other = 0.0;
01163 break;
01164 }
01165 case 8:{
01166 other = 0.0;
01167 break;
01168 }
01169 case 9:{
01170 other = 0.0;
01171 break;
01172 }
01173 case 10:{
01174 other = 0.0;
01175 break;
01176 }
01177 default:{
01178 print_err("unknown type of component is required",__FILE__,__LINE__,__func__);
01179 }
01180 }
01181 return (other);
01182
01183 }
01184
01185
01186
01187
01188
01189
01190 void kunmat::print_othervalue_name(FILE *out,long compother)
01191 {
01192 switch (compother){
01193 case 0:{
01194 fprintf (out,"Relative humidity () ");
01195 break;
01196 }
01197 case 1:{
01198 fprintf (out,"Temperature (K) ");
01199 break;
01200 }
01201 case 2:{
01202 fprintf (out,"Moisture content w (m3/m3) ");
01203 break;
01204 }
01205 case 3:{
01206 fprintf (out,"Water vapour pressure (Pa) ");
01207 break;
01208 }
01209 case 4:{
01210 fprintf (out,"Moisture content u (kg/kg) ");
01211 break;
01212 }
01213 default:{
01214 print_err("unknown type of component is required",__FILE__,__LINE__,__func__);
01215 }
01216 }
01217 }
01218
01219
01220 void kunmat::give_data(double rh,double Mhmc, double Smc, double Mhrh, double & moistakt)
01221 {
01222 double hvezdicka;
01223 hvezdicka = Mhmc/(1-sqrt(1-Mhrh));
01224
01225 if (rh < Mhrh)
01226 {
01227 moistakt = (1-sqrt(1-rh)) * hvezdicka ;
01228 }
01229 else {
01230 if (rh > 1)
01231 {
01232 moistakt = Smc;
01233 }
01234 else
01235 {
01236 moistakt = Mhmc + (rh - Mhrh )/(1-Mhrh )* (Smc - Mhmc );
01237 }
01238 }
01239
01240 }
01241
01242
01243
01244
01245
01246
01247
01248
01249
01250
01251 double kunmat::kapa_exp(double a, double b, double rh)
01252 {
01253 return (a * exp(b * rh));
01254 }
01255
01256
01257
01258
01259
01260
01261
01262
01263
01264
01265
01266
01267
01268
01269
01270
01271
01272
01273
01274
01275
01276
01277
01278
01279
01280
01281
01282
01283
01284
01285
01286
01287
01288 void kunmat::CorD(long charid,long &kvyhl, double x, double &y, double &z, double &z2)
01289 {
01290
01291 switch (MatChar[charid]){
01292 case 0:{
01293 kvyhl = 0;
01294 break;
01295 }
01296 case 1:{
01297
01298
01299 y = data[charid]->getval(0.0);
01300 kvyhl = 1;
01301 break;
01302 }
01303 case 2:{
01304
01305
01306 y = data[charid]->getval(x);
01307 kvyhl = 2;
01308
01309
01310 long nr=data[charid]->tabf[0].asize;
01311
01312
01313 z2 = data[charid]->tabf[0].y[nr-1];
01314
01315 break;
01316 }
01317 case 30:{
01318 y = MatFunce[charid][0];
01319 z = MatFunce[charid][1];
01320 kvyhl = 30;
01321 break;
01322 }
01323 case 31:{
01324 y = MatFunce[charid][0];
01325 z = MatFunce[charid][1];
01326 if (charid == 6)
01327 {
01328 z2 = MatFunce[charid][2];
01329 }
01330 kvyhl = 31;
01331 break;
01332 }
01333 case 32:{
01334 y = MatFunce[charid][0];
01335 z = MatFunce[charid][1];
01336 kvyhl = 32;
01337 break;
01338 }
01339 case 33:{
01340 kvyhl = 33;
01341
01342 break;
01343 }
01344 default:{
01345 print_err("unknown definition of material parameter is required",__FILE__,__LINE__,__func__);
01346 }
01347 }
01348 }
01349
01350 void kunmat::sorption_izoterms_giva_data(long kod,double rh, double tk, double & moist, double & dmoistdrh, long ipp)
01351 {
01352
01353
01354
01355 long s;
01356 double hmrh, mhmc, smc, u, a, n;
01357
01358 if (rh < 0) rh = 0;
01359
01360 CorD(6,s,0,a1,a2,a3);
01361
01362 switch (s){
01363 case 2:{
01364 if(rh>0.976)
01365 {
01366
01367 CorD(7,kd,rh,smc,a2,a3);
01368 CorD(6,kd,rh,mhmc,a2,mhmc);
01369 give_data(rh,mhmc, smc, 0.976,moist);
01370 dmoistdrh =1000 * (smc - mhmc)/(1 - 0.976);
01371 }
01372 else
01373 {
01374 CorD(6,s,rh,moist,dmoistdrh,a3);
01375 dmoistdrh =1000* derivative_sorption_izoterm_data(rh);
01376 }
01377 break;
01378 }
01379 case 30:{
01380 CorD(6,s,rh,a1,hmrh,a3);
01381 CorD(7,s,rh,smc,a2,a3);
01382 CorD(6,s,rh,mhmc,hmrh,a3);
01383 give_data(rh,mhmc, smc, hmrh,moist);
01384 if (rh < hmrh)
01385 {
01386 dmoistdrh = DerivativeOfTheSorptionIsotherm (rh, tk);
01387 }
01388 else {
01389 dmoistdrh = DerivativeOfTheRetentionCurve(rh, tk, ipp);
01390 }
01391
01392 break;
01393 }
01394 case 31:{
01395 switch (kod){
01396 case 0:{
01397 CorD(6,kod,0,u,a,n);
01398 moist = si_kk_hansen(rh, 0,u,a,n);
01399 dmoistdrh = sorptionizothermDerivation(rh, tk);
01400 break;
01401 }
01402 case 1:{
01403 dmoistdrh = sorptionizothermDerivation(rh, tk);
01404 break;
01405 }
01406 default:{
01407 print_err("unknown definition of Sorption isotherm HANSEN is required",__FILE__,__LINE__,__func__);
01408 }
01409 }
01410 break ;
01411 }
01412 default:{
01413 print_err("unknown definition of Sorption isotherm is required",__FILE__,__LINE__,__func__);
01414 }
01415 }
01416
01417
01418 }
01419
01420 double kunmat::si_kk_hansen(double rh, double tk,double u, double a, double n)
01421 {
01422 return (u*pow((1-(log(rh))/a),(-1/n)));
01423 }
01424
01425
01426
01427
01428
01429
01430
01431
01432 double kunmat::derivative_sorption_izoterm_data(double x1)
01433 {
01434 long charid;
01435 double derbi;
01436
01437
01438 charid=6;
01439 derbi = data[charid]->getderiv (x1);
01440
01441 return derbi;
01442 }
01443
01444 void kunmat::aux_values (long ipp,double *in,double *inp, double *ineq,double *out)
01445 {
01446 double x1,x2,x1p, x2p, ypv1,ypv2,ypv3,ypv4,ypv5, ypv0;
01447
01448 double w,dwdf, kapaA;
01449 double t;
01450 t = Tp->time;
01451 x1 = in[0];
01452 x2 = in[1];
01453
01454 x1p = inp[0];
01455 x2p = inp[1];
01456
01457
01458 ypv0 =ineq[0];
01459 ypv1 =ineq[1];
01460 ypv2 =ineq[2];
01461 ypv3 =ineq[3];
01462 ypv4 =ineq[4];
01463 ypv5 =ineq[5];
01464
01465 sorption_izoterms_giva_data(0,x1,x2,w,dwdf, ipp);
01466
01467 kapa_values(MatChar[5],ipp,w,ypv0,ypv2, kapaA);
01468 out[3] = kapaA;
01469
01470
01471 out[5] = -10000.0;
01472
01473 out[0] = w;
01474 out[1] = dwdf;
01475
01476 out[4] = -10000.0;
01477
01478
01479
01480 out[2] = data[7]->getval(0.0);
01481
01482
01483
01484
01485 }
01486
01487
01488
01489
01490
01491
01492
01493
01494
01495 void kunmat::save_values (long ipp,double *out)
01496 {
01497 Tm->ip[ipp].eqother[0]=out[0];
01498 Tm->ip[ipp].eqother[1]=out[1];
01499 Tm->ip[ipp].eqother[2]=out[2];
01500 Tm->ip[ipp].eqother[3]=out[3];
01501 }
01502
01503
01504
01505
01506
01507
01508
01509
01510
01511
01512 void kunmat::give_values (long ipp,double *av,double *eq)
01513 {
01514 av[0] = Tm->ip[ipp].av[0];
01515 av[1] = Tm->ip[ipp].av[1];
01516
01517 eq[0] = Tm->ip[ipp].eqother[0];
01518 eq[1] = Tm->ip[ipp].eqother[1];
01519 eq[2] = Tm->ip[ipp].eqother[2];
01520 eq[3] = Tm->ip[ipp].eqother[3];
01521 }
01522
01523 void kunmat::kapa_values (long kod, long ipp,double x1, double xpv, double ineq1, double &kapa)
01524 {
01525 double kapa1, kapa2, kapa3;
01526 double smc, a, n;
01527 long kd1;
01528
01529
01530 CorD(7,kd1,x1,smc,a,n);
01531
01532 if( x1 > smc)
01533 {
01534 x1 = smc;
01535 }
01536
01537 switch (kod){
01538 case 0: kapa = 0.0;
01539 break;
01540 case 1: CorD(5,kod,x1,kapa,kapa1,kapa2);
01541 break;
01542 case 2: CorD(5,kod,x1,kapa,kapa1,kapa2);
01543 break;
01544 case 30:{
01545 CorD(5,kod,x1,kapa1,kapa2,kapa3);
01546 kapa = kapa_exp(kapa1,kapa2,x1);
01547 break;
01548 }
01549 default:{
01550 print_err("unknown kod is required",__FILE__,__LINE__,__func__);
01551 }
01552 }
01553 }
01554
01555 void kunmat::initvalues (long ipp,long ido)
01556 {
01557 double x1,x2;
01558
01559 x1 = Tm->ip[ipp].av[0];
01560 x2 = Tm->ip[ipp].av[1];
01561
01562 }
01563
01564