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
00029
00030 #include "kunzel.h"
00031 #include "globalt.h"
00032 #include "math.h"
00033 #include "globmatt.h"
00034
00035 dampermeability kunmat::damper;
00036
00037 kunmat::kunmat()
00038 {
00039 long i;
00040
00041
00042 kunzeltype = 0;
00043
00044
00045 kmmtype = 0;
00046
00047
00048 rho_m = 0.0;
00049
00050 rho_w = 1000.0;
00051
00052 M = 0.01801528;
00053
00054 R = 8.314472;
00055
00056
00057 asi = 0.0;
00058
00059 adsi = 0.0;
00060
00061 akappa = 0.0;
00062
00063 adkappa = 0.0;
00064
00065
00066 time_hyst =0.0;
00067
00068
00069 daminfl=off;
00070
00071 for (i=0;i<20;i++){
00072 MatChar[i]=0;
00073 data[i] = NULL;
00074 }
00075 }
00076
00077 kunmat::~kunmat()
00078 {
00079 for (int i=0; i<20; i++)
00080 if (data[i])
00081 delete data[i];
00082 }
00083
00084
00085
00086
00087
00088
00089
00090
00091 void kunmat::read (XFILE *in)
00092 {
00093
00094 xfscanf (in,"%ld %ld",&kunzeltype,&kmmtype);
00095
00096 if (kunzeltype==8) {
00097
00098 rho.read (in);
00099
00100 dpp.read (in);
00101
00102
00103 transpar.read (in);
00104
00105
00106 sorpiso.read (in);
00107
00108
00109 sm.read (in);
00110
00111
00112 c.read (in);
00113
00114
00115 lambda.read (in);
00116
00117 }
00118
00119 else{
00120
00121 rho.read (in);
00122
00123
00124 por.read (in);
00125
00126
00127 mu.read (in);
00128
00129
00130 kappa.read (in);
00131
00132 if (kunzeltype == 11){
00133
00134 xfscanf (in,"%lf",&time_hyst);
00135
00136
00137 xfscanf (in,"%lf",&akappa);
00138
00139
00140 kappadry.read (in);
00141
00142
00143 xfscanf (in,"%lf",&adkappa);
00144 }
00145
00146
00147 sorpiso.read (in);
00148
00149 if (kunzeltype == 11){
00150
00151 xfscanf (in,"%lf",&asi);
00152
00153
00154 desorpiso.read (in);
00155
00156
00157 xfscanf (in,"%lf",&adsi);
00158 }
00159
00160
00161 sm.read (in);
00162
00163
00164 c.read (in);
00165
00166
00167 lambda.read (in);
00168 }
00169
00170 xfscanf (in,"%m",&flagsw_kwdset,&daminfl);
00171 }
00172
00173
00174
00175
00176
00177
00178
00179
00180
00181 void kunmat::print (FILE *out)
00182 {
00183
00184 fprintf (out,"\n %ld %ld\n",kunzeltype,kmmtype);
00185
00186 if (kunzeltype==8)
00187 {
00188
00189 rho.print (out);
00190
00191 dpp.print (out);
00192
00193 transpar.print (out);
00194
00195 sorpiso.print (out);
00196
00197 sm.print (out);
00198
00199 c.print (out);
00200
00201 lambda.print (out);
00202 }
00203 else
00204 {
00205
00206 rho.print (out);
00207
00208 por.print (out);
00209
00210 mu.print (out);
00211
00212 kappa.print (out);
00213
00214 if (kunzeltype == 11){
00215
00216 fprintf (out,"\n %lf",time_hyst);
00217
00218
00219 fprintf (out,"\n %lf",akappa);
00220
00221
00222 kappadry.print (out);
00223
00224 fprintf (out,"\n %lf",adkappa);
00225 }
00226
00227
00228
00229 sorpiso.print (out);
00230
00231 if (kunzeltype == 11){
00232 fprintf (out,"\n %lf",asi);
00233
00234
00235 desorpiso.print (out);
00236
00237 fprintf (out,"\n %lf",adsi);
00238 }
00239
00240
00241 sm.print (out);
00242
00243 c.print (out);
00244
00245 lambda.print (out);
00246 }
00247
00248
00249 }
00250
00251
00252
00253
00254
00255
00256
00257
00258
00259
00260
00261
00262
00263
00264
00265
00266
00267
00268
00269 void kunmat::matcond (matrix &d,long ri,long ci,long ipp)
00270 {
00271 long n;
00272 n = d.n;
00273
00274 switch (n){
00275 case 1:{
00276 matcond1d (d,ri,ci,ipp);
00277 break;
00278 }
00279 case 2:{
00280 matcond2d (d,ri,ci,ipp);
00281 break;
00282 }
00283 case 3:{
00284 matcond3d (d,ri,ci,ipp);
00285 break;
00286 }
00287 default:{
00288 print_err("unknown number of components of conductivity tensor is required",__FILE__,__LINE__,__func__);
00289 }
00290 }
00291 if (daminfl == on){
00292
00293 if((ri == 0) && (ci == 0))
00294 damper.matcond (d,ipp);
00295 }
00296 }
00297
00298
00299
00300
00301
00302
00303
00304
00305
00306
00307
00308
00309 void kunmat::matcond1d (matrix &d,long ri,long ci,long ipp)
00310 {
00311 double k;
00312 k = 0.0;
00313
00314 if((ri == 0) && (ci == 0))
00315 k = kmm(ipp);
00316 if((ri == 0) && (ci == 1))
00317 k = kmt(ipp);
00318 if((ri == 1) && (ci == 0))
00319 k = khm(ipp);
00320 if((ri == 1) && (ci == 1))
00321 k = kht(ipp);
00322
00323 d[0][0] = k;
00324 }
00325
00326
00327
00328
00329
00330
00331
00332
00333
00334
00335
00336 void kunmat::matcond2d (matrix &d,long ri,long ci,long ipp)
00337 {
00338 double k;
00339 k = 0.0;
00340
00341
00342 if((ri == 0) && (ci == 0))
00343 k = kmm(ipp);
00344 if((ri == 0) && (ci == 1))
00345 k = kmt(ipp);
00346 if((ri == 1) && (ci == 0))
00347 k = khm(ipp);
00348 if((ri == 1) && (ci == 1))
00349 k = kht(ipp);
00350
00351 fillm(0.0,d);
00352
00353
00354
00355
00356 d[0][0] = k; d[0][1] = 0.0;
00357 d[1][0] = 0.0; d[1][1] = k;
00358 }
00359
00360
00361
00362
00363
00364
00365
00366
00367
00368
00369
00370
00371 void kunmat::matcond3d (matrix &d,long ri,long ci,long ipp)
00372 {
00373 double k;
00374 k = 0.0;
00375
00376
00377 if((ri == 0) && (ci == 0))
00378 k = kmm(ipp);
00379 if((ri == 0) && (ci == 1))
00380 k = kmt(ipp);
00381 if((ri == 1) && (ci == 0))
00382 k = khm(ipp);
00383 if((ri == 1) && (ci == 1))
00384 k = kht(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
00403 void kunmat::matcap (double &c,long ri,long ci,long ipp)
00404 {
00405 c=0.0;
00406
00407 if((ri == 0) && (ci == 0))
00408 c = cmm(ipp);
00409 if((ri == 0) && (ci == 1))
00410 c = cmt(ipp);
00411 if((ri == 1) && (ci == 0))
00412 c = chm(ipp);
00413 if((ri == 1) && (ci == 1)){
00414 c = cht(ipp);
00415 }
00416 }
00417
00418
00419
00420
00421
00422
00423
00424
00425
00426 void kunmat::values_correction (vector &nv)
00427 {
00428 double rh, tk, pv, pvs;
00429
00430 switch (kunzeltype){
00431 case 1:{
00432
00433 rh=nv[0];
00434
00435 tk=nv[1];
00436
00437 if (rh >= 1.0)
00438 rh = 1.0;
00439 if (rh <= 0.0)
00440 rh = 0.0;
00441
00442 if (tk >= 350.0)
00443 tk = 350.0;
00444 if (tk <= 240.0)
00445 tk = 240.0;
00446
00447 nv[0]=rh;
00448 nv[1]=tk;
00449 break;
00450 }
00451 case 2:
00452 case 3:
00453 case 4:{
00454
00455 pv=nv[0];
00456
00457 tk=nv[1];
00458
00459 pvs = saturated_water_vapor_pressure(tk);
00460
00461 rh = pv/pvs;
00462
00463 if(rh >= 1.0)
00464 pv = pvs;
00465 if(rh <= 0.0)
00466
00467
00468 if (tk >= 350.0)
00469 tk = 350.0;
00470 if (tk <= 240.0)
00471 tk = 240.0;
00472
00473 nv[0]=pv;
00474 nv[1]=tk;
00475 break;
00476 }
00477 case 5:
00478 case 6:
00479 case 7:{
00480
00481 pv=nv[0];
00482
00483 tk=nv[1];
00484
00485
00486 nv[0]=pv;
00487 nv[1]=tk;
00488 break;
00489 }
00490 case 8:{
00491
00492 pv=nv[0];
00493
00494 tk=nv[1];
00495
00496 if(pv >= 3.16995416e+03)
00497
00498 pv = 3.1699e+03;
00499 if(pv <= 3.17014677e-06)
00500 pv = 3.18e-06;
00501
00502 nv[0]=pv;
00503 nv[1]=tk;
00504 break;
00505 }
00506 default:{
00507 print_err("unknown type of Kunzel is required",__FILE__,__LINE__,__func__);
00508 }
00509 }
00510
00511 }
00512
00513
00514
00515
00516
00517
00518
00519
00520
00521 double kunmat::kmm (long ipp)
00522 {
00523 double rh, dp, dmretc, ps, dphi, moist, kk, tk, kapa, pv;
00524 double A,B,r1,r2;
00525
00526 switch (kunzeltype) {
00527 case 1:
00528 case 11:{
00529
00530
00531 rh = Tm->ip[ipp].av[0];
00532
00533
00534 tk = Tm->ip[ipp].av[1];
00535
00536 if (rh>1.0)
00537 rh = 1.0;
00538 if (rh<0.0)
00539 rh = 0.0;
00540
00541
00542
00543
00544 r1 = 0.90;
00545 r2 = 0.976;
00546
00547
00548 moist = Tm->ip[ipp].eqother[0];
00549 A = 0.0;
00550 B = 0.0;
00551
00552 if (rh>=0.0){
00553 if (rh<=r1){
00554 switch (kmmtype){
00555 case 0:{
00556 A=1.0;
00557 B=1.0;
00558 break;
00559 }
00560 case 1:
00561 case 2:
00562 case 3:{
00563 A = 1.0;
00564 B = 0.0;
00565 break;
00566 }
00567 default:{
00568 print_err("unknown type of Kunzel model is required",__FILE__,__LINE__,__func__);
00569 }
00570 }
00571 }
00572 else{
00573 if (rh > r2){
00574 switch (kmmtype){
00575 case 0:{
00576 A=1.0;
00577 B=1.0;
00578 break;
00579 }
00580 case 1:
00581 case 2:
00582 case 3:{
00583 A = 0.0;
00584 B = 1.0;
00585 break;
00586 }
00587 default:{
00588 print_err("unknown type of Kunzel model is required",__FILE__,__LINE__,__func__);
00589 }
00590 }
00591 }
00592 else{
00593 switch (kmmtype){
00594 case 0:{
00595 A=1.0;
00596 B=1.0;
00597 break;
00598 }
00599 case 1:{
00600 B = (1/(r2-r1))*(rh-r1);
00601 A = 1-B;
00602 break;
00603 }
00604 case 2:{
00605 B= 0.5*sin(M_PI/(r2-r1)*(rh-r1)-M_PI/2)+0.5;
00606 A=1-B;
00607 break;
00608 }
00609 case 3:{
00610 if (rh<0.938){
00611 B=32*pow((1/(r2-r1))*(rh-r1),6);
00612 A=1-B;
00613 }
00614 else{
00615 B=1-32*pow((1/(r2-r1))*(r2-rh),6);
00616 A=1-B;
00617 }
00618 break;
00619 }
00620 default:{
00621 print_err("unknown definition of Material KUNZEL is required",__FILE__,__LINE__,__func__);
00622 }
00623 }
00624 }
00625 }
00626 }
00627 else{
00628 print_err("unknown error",__FILE__,__LINE__,__func__);
00629 }
00630
00631
00632
00633
00634
00635 kapa = Tm->ip[ipp].eqother[3];
00636
00637 dp = water_vapour_permeability(tk,ipp);
00638 ps = saturated_water_vapor_pressure( tk);
00639
00640
00641 dmretc = Tm->ip[ipp].eqother[1];
00642
00643 dphi = kapa * dmretc*rho_w;
00644 kk = dphi*B + A*dp*ps;
00645
00646 break;
00647 }
00648
00649
00650
00651
00652
00653
00654
00655
00656
00657
00658
00659
00660
00661
00662
00663
00664
00665
00666
00667
00668
00669
00670
00671
00672
00673
00674
00675
00676
00677
00678
00679
00680
00681
00682
00683
00684
00685
00686
00687
00688
00689
00690
00691
00692
00693
00694
00695
00696
00697
00698
00699
00700
00701
00702
00703
00704
00705
00706
00707
00708
00709
00710
00711
00712
00713
00714
00715
00716
00717
00718
00719
00720
00721
00722
00723
00724
00725
00726
00727
00728 case 8:{
00729
00730
00731 pv = Tm->ip[ipp].av[0];
00732
00733
00734
00735 double transp;
00736 transp = transpar.getval (pv);
00737
00738
00739
00740 kk = transp;
00741 break;
00742 }
00743 default:{
00744 print_err("unknown type of Kunzel is required",__FILE__,__LINE__,__func__);
00745
00746 }
00747 }
00748
00749 return (kk);
00750 }
00751
00752
00753
00754
00755
00756
00757
00758
00759
00760
00761
00762 double kunmat::kmt(long ipp)
00763 {
00764 double kk, tk, rh;
00765 double dp, dpvs;
00766
00767 switch (kunzeltype) {
00768 case 1:
00769 case 11:{
00770
00771
00772 rh = Tm->ip[ipp].av[0];
00773
00774
00775 tk = Tm->ip[ipp].av[1];
00776
00777 if (rh>1.0)
00778 rh = 1.0;
00779 if (rh<0.0)
00780 rh = 0.0;
00781
00782 dpvs = derivative_of_saturated_water_vapor_pressure(tk);
00783 dp = water_vapour_permeability(tk,ipp);
00784
00785
00786 kk = dp * rh * dpvs;
00787 break;
00788 }
00789 case 2:{
00790 double pv,tk, dw;
00791
00792 pv = Tm->ip[ipp].av[0];
00793
00794 tk = Tm->ip[ipp].av[1];
00795
00796
00797 dw = Tm->ip[ipp].eqother[3];
00798
00799 kk = -dw*M*pv/(R*tk*tk);
00800
00801 break;
00802 }
00803
00804 case 3:
00805 case 4:
00806 case 5:
00807 case 6:
00808 case 7:
00809 case 8:{
00810 kk = 0.0;
00811
00812 break;
00813 }
00814 default:{
00815 print_err("unknown type of Kunzel is required",__FILE__,__LINE__,__func__);
00816
00817 }
00818
00819 }
00820
00821
00822
00823
00824 return (kk);
00825 }
00826
00827
00828
00829
00830
00831
00832
00833
00834
00835
00836 double kunmat::kht(long ipp)
00837 {
00838 double kk,tk,rh;
00839 double lv1, lamb, dpvs, dp, moist, rov;
00840
00841 switch (kunzeltype){
00842 case 1:
00843 case 11:{
00844
00845
00846 rh = Tm->ip[ipp].av[0];
00847
00848
00849 tk = Tm->ip[ipp].av[1];
00850
00851 if (rh>1.0)
00852 rh = 1.0;
00853 if (rh<0.0)
00854 rh = 0.0;
00855
00856
00857 moist = Tm->ip[ipp].eqother[0];
00858
00859 lamb = lambda.getval (moist);
00860
00861 dpvs = derivative_of_saturated_water_vapor_pressure(tk);
00862 dp = water_vapour_permeability(tk,ipp);
00863
00864 lv1 = latent_heat_of_evaporation_of_water (tk);
00865
00866
00867 kk = lamb + lv1 * dp * rh * dpvs;
00868 break;
00869 }
00870 case 2:
00871 case 3:
00872 case 4:
00873 case 5:
00874 case 6:{
00875 double moist;
00876
00877 moist = Tm->ip[ipp].eqother[0];
00878
00879
00880
00881 kk = lambda.getval (moist);
00882 break;
00883 }
00884 case 7:{
00885
00886
00887 rov = Tm->ip[ipp].av[0];
00888
00889
00890 tk = Tm->ip[ipp].av[1];
00891
00892
00893 moist = 0.01;
00894
00895
00896 lamb = lambda.getval (moist);
00897 dp = water_vapour_permeability(tk,ipp);
00898
00899 lv1 = latent_heat_of_evaporation_of_water (tk);
00900
00901
00902 kk = lamb + lv1 * dp * R * rov/M;
00903
00904 break;
00905 }
00906 case 8:{
00907
00908
00909
00910
00911
00912
00913
00914
00915
00916
00917 moist = Tm->ip[ipp].eqother[0];
00918
00919 lamb = lambda.getval (moist);
00920
00921
00922 kk = lamb;
00923
00924 break;
00925 }
00926
00927 default:{
00928 print_err("unknown type of Kunzel is required",__FILE__,__LINE__,__func__);
00929
00930 }
00931 }
00932
00933 return (kk);
00934 }
00935
00936
00937
00938
00939
00940
00941
00942
00943
00944
00945 double kunmat::khm (long ipp)
00946 {
00947 double kk,tk;
00948 double lv,dp,ps,moist,D, pv;
00949
00950 switch (kunzeltype){
00951 case 1:
00952 case 11:{
00953
00954
00955 tk = Tm->ip[ipp].av[1];
00956
00957 dp = water_vapour_permeability(tk,ipp);
00958 ps = saturated_water_vapor_pressure (tk);
00959 lv = latent_heat_of_evaporation_of_water (tk);
00960
00961
00962 kk = lv * dp * ps;
00963
00964 break;
00965 }
00966 case 2:
00967 case 3:
00968 case 4:
00969 case 5:{
00970
00971 tk = Tm->ip[ipp].av[1];
00972
00973 dp = water_vapour_permeability(tk,ipp);
00974 lv = latent_heat_of_evaporation_of_water (tk);
00975
00976 kk = lv*dp;
00977
00978 break;
00979 }
00980 case 6:{
00981
00982 tk = Tm->ip[ipp].av[1];
00983
00984 dp = water_vapour_permeability(tk,ipp);
00985 lv = latent_heat_of_evaporation_of_water (tk);
00986
00987 kk = lv*dp;
00988
00989 break;
00990 }
00991 case 7:{
00992
00993
00994
00995
00996
00997 tk = Tm->ip[ipp].av[1];
00998
00999
01000
01001
01002 moist= 0.05;
01003 CorD(18,kd,moist,D,a2,a3);
01004 lv = latent_heat_of_evaporation_of_water (tk);
01005 dp = water_vapour_permeability(tk,ipp);
01006
01007 kk = lv * dp*R*tk/M;
01008
01009 break;
01010 }
01011 case 8:{
01012
01013
01014 tk = Tm->ip[ipp].av[1];
01015
01016
01017 pv = Tm->ip[ipp].av[0];
01018
01019
01020 dp = dpp.getval (pv);
01021
01022
01023 lv = latent_heat_of_evaporation_of_water (tk);
01024
01025 kk = lv*dp;
01026
01027 break;
01028 }
01029 default:{
01030 print_err("unknown type of Kunzel is required",__FILE__,__LINE__,__func__);
01031
01032 }
01033 }
01034
01035 return (kk);
01036 }
01037
01038
01039
01040
01041
01042
01043
01044
01045
01046
01047 double kunmat::cmm(long ipp)
01048 {
01049 double cc,tk;
01050 double moist, pi;
01051
01052 switch (kunzeltype){
01053 case 1:
01054 case 11:{
01055
01056 cc = Tm->ip[ipp].eqother[1]*rho_w;
01057 break;
01058 }
01059 case 2:{
01060 double tk;
01061
01062
01063
01064 tk = Tm->ip[ipp].av[1];
01065
01066 cc = M/(R*tk);
01067
01068 break;
01069 }
01070 case 3:{
01071
01072
01073 moist = Tm->ip[ipp].eqother[0];
01074
01075
01076 pi = por.getval (moist);
01077
01078 cc = (pi-moist);
01079
01080 break;
01081 }
01082 case 4:{
01083
01084
01085
01086 tk = Tm->ip[ipp].av[1];
01087
01088
01089 moist = Tm->ip[ipp].eqother[0];
01090
01091
01092 pi = por.getval (moist);
01093
01094 cc = (pi-moist)*M/(R*tk);;
01095
01096 break;
01097 }
01098 case 5:{
01099
01100
01101
01102 tk = Tm->ip[ipp].av[1];
01103
01104
01105 pi = por.getval (moist);
01106
01107 cc = (pi)*M/(R*tk);;
01108
01109 break;
01110 }
01111 case 6:{
01112
01113
01114
01115 tk = Tm->ip[ipp].av[1];
01116
01117
01118 pi = por.getval (moist);
01119
01120 cc = (pi)*M/(R*tk);
01121
01122 break;
01123 }
01124 case 7:{
01125
01126
01127
01128
01129
01130
01131 pi = por.getval (moist);
01132
01133 cc = pi;
01134
01135 cc = 1.0;
01136
01137 break;
01138 }
01139 case 8:{
01140 double pv, t;
01141
01142 pv = Tm->ip[ipp].av[0];
01143
01144
01145
01146
01147
01148 cc = Tm->ip[ipp].eqother[1]*rho_w;
01149
01150
01151
01152
01153
01154
01155
01156
01157
01158 break;
01159 }
01160 default:{
01161 print_err("unknown type of Kunzel is required",__FILE__,__LINE__,__func__);
01162
01163 }
01164
01165
01166 }
01167
01168 return (cc);
01169 }
01170
01171
01172
01173
01174
01175
01176
01177
01178
01179
01180 double kunmat::cmt(long ipp)
01181 {
01182 double cc,tk,pv;
01183
01184 switch (kunzeltype){
01185 case 1:
01186 case 11:{
01187 cc = 0.0;
01188
01189 break;
01190 }
01191 case 2:{
01192
01193 pv = Tm->ip[ipp].av[0];
01194
01195 tk = Tm->ip[ipp].av[1];
01196
01197 cc = -M*pv/(R*tk*tk);
01198 break;
01199 }
01200 case 3:
01201 case 4:
01202 case 5:
01203 case 6:
01204 case 7:
01205 case 8:{
01206 cc = 0.0;
01207 break;
01208 }
01209
01210 default:{
01211 print_err("unknown type of Kunzel is required",__FILE__,__LINE__,__func__);
01212 }
01213
01214 }
01215
01216 return (cc);
01217 }
01218
01219
01220
01221
01222
01223
01224
01225
01226
01227
01228 double kunmat::chm(long ipp)
01229 {
01230 double cc;
01231
01232 switch (kunzeltype){
01233 case 1:
01234 case 11:{
01235 cc = 0.0;
01236 break;
01237 }
01238 case 2:
01239 case 3:
01240 case 4:
01241 case 5:
01242 case 6:
01243 case 7:
01244 case 8:{
01245 cc = 0.0;
01246
01247 break;
01248 }
01249
01250 default:{
01251 print_err("unknown type of Kunzel is required",__FILE__,__LINE__,__func__);
01252
01253 }
01254
01255 }
01256
01257 return (cc);
01258 }
01259
01260
01261
01262
01263
01264
01265
01266
01267
01268
01269 double kunmat::cht(long ipp)
01270 {
01271 double cc;
01272
01273 switch (kunzeltype) {
01274 case 1:
01275 case 11:{
01276 cc = derivative_of_the_enthalpy_density(ipp);
01277
01278 break;
01279 }
01280 case 2:
01281 case 3:
01282 case 4:
01283 case 5:
01284 case 6:
01285 case 7:
01286 case 8:{
01287 cc = derivative_of_the_enthalpy_density(ipp);
01288
01289 break;
01290 }
01291
01292 default:{
01293 print_err("unknown type of Kunzel is required",__FILE__,__LINE__,__func__);
01294 }
01295
01296 }
01297
01298 return (cc);
01299 }
01300
01301
01302
01303
01304
01305
01306
01307
01308 double kunmat::derivative_of_saturated_water_vapor_pressure(double tk)
01309 {
01310
01311 return ((4042.9 * exp(23.5771-4042.9/(tk - 37.58)))/((tk-37.58)*(tk-37.58)));
01312 }
01313
01314
01315
01316
01317
01318
01319
01320
01321
01322
01323
01324
01325
01326
01327 double kunmat::water_vapour_permeability (double tk,long ipp)
01328 {
01329 double P, Rv, Pa;
01330 double da, dp, mi;
01331
01332
01333 moist = Tm->ip[ipp].eqother[0];
01334
01335
01336
01337 mi = mu.getval (moist);
01338
01339
01340 Rv = R/M;
01341
01342
01343 Pa = 101325.0;
01344 P = 101325.0;
01345
01346 da = (2.306e-5 * Pa)/(Rv * tk * P)*pow((tk/273.15),1.81);
01347
01348 if (mi==0.0)
01349 dp=0.0;
01350 else
01351 dp = da/mi;
01352
01353 return (dp);
01354 }
01355
01356
01357
01358
01359
01360
01361 double kunmat::latent_heat_of_evaporation_of_water(double tk)
01362 {
01363 return (2.5008e6)*pow((273.15/tk),(0.167+tk*3.67e-4));
01364 }
01365
01366
01367
01368
01369
01370
01371
01372
01373
01374
01375
01376
01377
01378 double kunmat::saturated_water_vapor_pressure(double tk)
01379 {
01380 return (exp(23.5771 - 4042.9/(tk - 37.58)));
01381 }
01382
01383
01384
01385
01386
01387
01388
01389
01390
01391
01392
01393 double kunmat::derivative_of_the_enthalpy_density (long ipp)
01394 {
01395 double cap, moist,rho_m;
01396
01397
01398 moist = Tm->ip[ipp].eqother[0];
01399
01400
01401 cap = c.getval (moist);
01402
01403 rho_m = rho.getval (0.0);
01404
01405 return (rho_m * cap);
01406 }
01407
01408
01409
01410
01411
01412
01413
01414
01415
01416
01417
01418
01419
01420
01421
01422
01423
01424
01425
01426 double kunmat::transmission_nodval (double nodval,long ri,long ci,long nid,long bc)
01427 {
01428
01429 long k;
01430 double new_nodval,h,t;
01431 new_nodval = 0.0;
01432
01433
01434 h = nodalval (0,nid,0);
01435
01436 t = nodalval (0,nid,1);
01437
01438 if((ri == 0) && (ci == 0))
01439 new_nodval = get_transmission_nodval_hh (nodval,h,t,bc);
01440 if((ri == 0) && (ci == 1))
01441 new_nodval = 0.0;
01442
01443 if((ri == 1) && (ci == 0))
01444 new_nodval = get_transmission_nodval_th (nodval,bc);
01445 if((ri == 1) && (ci == 1))
01446 new_nodval = get_transmission_nodval_tt (nodval,bc);
01447
01448 return (new_nodval);
01449 }
01450
01451
01452
01453
01454
01455
01456
01457
01458
01459
01460 double kunmat::get_transmission_nodval_hh (double bv,double rh,double t,long bc)
01461 {
01462 double nodval,pgws,moist;
01463
01464 switch (bc){
01465 case 5:{
01466 nodval = bv;
01467 break;
01468 }
01469 case 30:{
01470
01471
01472
01473
01474 nodval = bv;
01475 break;
01476 }
01477 case 31:{
01478
01479 nodval = bv;
01480 break;
01481 }
01482 case 32:{
01483
01484
01485 pgws = saturated_water_vapor_pressure(t);
01486 nodval = pgws*rh;
01487 bv = pgws*bv;
01488
01489 nodval = bv - nodval;
01490
01491 break;
01492 }
01493 case 33:{
01494
01495
01496
01497 rho_m = rho.getval (0.0);
01498
01499 moist = bv/100.0*rho_m/rho_w;
01500
01501
01502 nodval = sorpiso.inverse_isotherm_value (moist);
01503 break;
01504 }
01505 default:{
01506 print_err("no acceptable boundary condition is prescribed",__FILE__,__LINE__,__func__);
01507 exit(0);
01508 }
01509 }
01510 return nodval;
01511 }
01512
01513
01514
01515
01516
01517
01518
01519
01520
01521
01522 double kunmat::get_transmission_nodval_th (double bv,long bc)
01523 {
01524 double nodval,pgws,moist;
01525
01526 switch (bc){
01527 case 4:{
01528 nodval = bv;
01529 break;
01530 }
01531 case 5:{
01532 nodval = bv;
01533 break;
01534 }
01535 case 30:{
01536 nodval = bv;
01537 break;
01538 }
01539 default:{
01540 print_err("no acceptable boundary condition is prescribed",__FILE__,__LINE__,__func__);
01541 exit(0);
01542 }
01543 }
01544 return nodval;
01545 }
01546
01547
01548
01549
01550
01551
01552
01553
01554
01555 double kunmat::get_transmission_nodval_tt (double bv,long bc)
01556 {
01557 double nodval;
01558
01559 switch (bc){
01560 case 5:{
01561 nodval = bv;
01562 break;
01563 }
01564 case 30:{
01565 nodval = bv;
01566 break;
01567 }
01568 default:{
01569 print_err("no acceptable boundary condition is prescribed",__FILE__,__LINE__,__func__);
01570 exit(0);
01571 }
01572 }
01573
01574 return (nodval);
01575 }
01576
01577
01578
01579
01580
01581
01582
01583
01584
01585
01586
01587
01588
01589 double kunmat::transmission_transcoeff (double trc,long ri,long ci,long nid,long bc)
01590 {
01591 long k;
01592 double new_trc,h,t;
01593 new_trc = 0.0;
01594
01595
01596 h = nodalval (0,nid,0);
01597
01598 t = nodalval (0,nid,1);
01599
01600 if((ri == 0) && (ci == 0))
01601 new_trc = get_transmission_transcoeff_hh (trc,bc);
01602 if((ri == 0) && (ci == 1))
01603 new_trc = 0.0;
01604
01605 if((ri == 1) && (ci == 0))
01606 new_trc = get_transmission_transcoeff_th (trc,bc);
01607 if((ri == 1) && (ci == 1))
01608 new_trc = get_transmission_transcoeff_tt (trc,bc);
01609
01610 return (new_trc);
01611 }
01612
01613
01614
01615
01616
01617
01618
01619
01620 double kunmat::get_transmission_transcoeff_hh (double t,long bc)
01621 {
01622 double trc,pgws;
01623
01624 switch (bc){
01625 case 5:{
01626 trc = t;
01627 break;
01628 }
01629 case 30:{
01630
01631
01632
01633 trc = t;
01634 break;
01635 }
01636 case 31:{
01637
01638 pgws = saturated_water_vapor_pressure (t);
01639 trc = pgws;
01640 trc = 1.0;
01641 break;
01642 }
01643 case 32:{
01644
01645 trc = 0.0;
01646 break;
01647 }
01648 case 33:{
01649
01650 trc = 1.0;
01651 break;
01652 }
01653 default:{
01654 print_err("no acceptable boundary condition is prescribed",__FILE__,__LINE__,__func__);
01655 exit(0);
01656 }
01657 }
01658
01659 return trc;
01660 }
01661
01662
01663
01664
01665
01666
01667
01668
01669 double kunmat::get_transmission_transcoeff_th (double t,long bc)
01670 {
01671 double trc;
01672
01673 switch (bc){
01674 case 4:{
01675 trc = t;
01676 break;
01677 }
01678 case 5:{
01679 trc = t;
01680 break;
01681 }
01682 case 30:{
01683 trc = t;
01684 break;
01685 }
01686 default:{
01687 print_err("no acceptable boundary condition is prescribed",__FILE__,__LINE__,__func__);
01688 exit(0);
01689 }
01690 }
01691
01692 return trc;
01693 }
01694
01695
01696
01697
01698
01699
01700
01701
01702 double kunmat::get_transmission_transcoeff_tt (double trcp,long bc)
01703 {
01704 double trc;
01705
01706 switch (bc){
01707 case 5:{
01708 trc = trcp;
01709 }
01710 case 30:{
01711
01712 trc = trcp;
01713 break;
01714 }
01715 default:{
01716 print_err("no acceptable boundary condition is prescribed",__FILE__,__LINE__,__func__);
01717 exit(0);
01718 }
01719 }
01720
01721 return trc;
01722 }
01723
01724
01725
01726
01727
01728
01729
01730
01731
01732
01733
01734
01735
01736
01737 double kunmat::transmission_flux (double nodval,long ri,long ci,long nid,long bc)
01738 {
01739 long k;
01740 double flux,h,t;
01741
01742 flux = 0.0;
01743
01744
01745 h = nodalval (0,nid,0);
01746
01747 t = nodalval (0,nid,1);
01748
01749 if((ri == 0) && (ci == 0))
01750 flux = get_transmission_flux_hh (nodval,h,t,bc);
01751 if((ri == 0) && (ci == 1))
01752 flux = 0.0;
01753
01754 if((ri == 1) && (ci == 0))
01755 flux = 0.0;
01756 if((ri == 1) && (ci == 1))
01757 flux = get_transmission_flux_tt (nodval,t,bc);
01758
01759 return flux;
01760 }
01761
01762
01763
01764
01765
01766
01767
01768
01769
01770
01771 double kunmat::get_transmission_flux_hh (double bv,double rh,double t,long bc)
01772 {
01773 double flux,pgws,moist,nodval;
01774
01775 switch (bc){
01776 case 30:{
01777
01778 pgws = saturated_water_vapor_pressure (t);
01779
01780 flux = pgws*rh;
01781 bv = pgws*bv;
01782 flux = bv - flux;
01783
01784 break;
01785 }
01786 case 31:{
01787
01788 pgws = saturated_water_vapor_pressure (t);
01789 flux = pgws*rh;
01790 flux = bv - flux;
01791 break;
01792 }
01793 case 32:{
01794
01795 pgws = saturated_water_vapor_pressure(t);
01796 flux = pgws*rh;
01797 bv = pgws*bv;
01798 flux = bv - flux;
01799
01800 break;
01801 }
01802 case 33:{
01803
01804
01805
01806 rho_m = rho.getval (0.0);
01807
01808 moist = bv/100*rho_m/rho_w;
01809
01810
01811
01812 nodval = sorpiso.inverse_isotherm_value (moist);
01813
01814 flux = nodval - rh;
01815 break;
01816 }
01817 default:{
01818 print_err("no acceptable boundary condition is prescribed",__FILE__,__LINE__,__func__);
01819 exit(0);
01820 }
01821 }
01822
01823 return flux;
01824 }
01825
01826
01827
01828
01829
01830
01831
01832
01833
01834
01835
01836
01837 double kunmat::get_transmission_flux_tt (double bv,double t,long bc)
01838 {
01839 double flux;
01840
01841 switch (bc){
01842 case 30:{
01843
01844 flux = (bv - t);
01845 break;
01846 }
01847 default:{
01848 print_err("no acceptable boundary condition is prescribed",__FILE__,__LINE__,__func__);
01849 exit(0);
01850 }
01851 }
01852 return(flux);
01853 }
01854
01855
01856
01857
01858
01859
01860
01861
01862
01863
01864
01865 double kunmat::get_othervalue(long compother,double rh,double t, long ipp)
01866 {
01867 double other,w;
01868
01869 switch (compother){
01870 case 0:{
01871 other = rh*100;
01872 break;
01873 }
01874 case 1:{
01875 other = t-273.15;
01876 break;
01877 }
01878 case 2:{
01879
01880 w = sorpiso.isotherm_value (rh);
01881 other = w;
01882 break;
01883 }
01884 case 3:{
01885 other = rh * saturated_water_vapor_pressure (t);
01886 break;
01887 }
01888 case 4:{
01889
01890 rho_m = rho.getval (0.0);
01891 other = w*rho_w/rho_m;
01892 break;
01893 }
01894 case 5:{
01895 t = 298.15;
01896 other = rh/saturated_water_vapor_pressure (t);
01897 break;
01898 }
01899 case 6:{
01900 other = rh*saturated_water_vapor_pressure (t)*M/(R*t);
01901 break;
01902 }
01903 case 7:{
01904 other = rh*R*t/(saturated_water_vapor_pressure (t)*M);
01905 break;
01906 }
01907 case 8:{
01908 other = 0.0;
01909 break;
01910 }
01911 case 9:{
01912 other = 0.0;
01913 break;
01914 }
01915 case 10:{
01916 other = 0.0;
01917 break;
01918 }
01919 default:{
01920 print_err("unknown type of component is required",__FILE__,__LINE__,__func__);
01921 }
01922 }
01923 return (other);
01924
01925 }
01926
01927
01928
01929
01930
01931
01932 void kunmat::print_othervalue_name(FILE *out,long compother)
01933 {
01934 switch (compother){
01935 case 0:{
01936 fprintf (out,"Relative humidity () ");
01937 break;
01938 }
01939 case 1:{
01940 fprintf (out,"Temperature (K) ");
01941 break;
01942 }
01943 case 2:{
01944 fprintf (out,"Moisture content w (m3/m3) ");
01945 break;
01946 }
01947 case 3:{
01948 fprintf (out,"Water vapour pressure (Pa) ");
01949 break;
01950 }
01951 case 4:{
01952 fprintf (out,"Moisture content u (kg/kg) ");
01953 break;
01954 }
01955 default:{
01956 print_err("unknown type of component is required",__FILE__,__LINE__,__func__);
01957 }
01958 }
01959 }
01960
01961
01962
01963
01964
01965
01966
01967
01968
01969
01970
01971
01972 double kunmat::kapa_exp(double a, double b, double rh)
01973 {
01974 return (a * exp(b * rh));
01975 }
01976
01977
01978
01979
01980
01981
01982
01983
01984
01985
01986
01987
01988
01989
01990
01991
01992
01993
01994
01995
01996
01997
01998
01999
02000
02001
02002
02003
02004
02005
02006
02007
02008
02009 void kunmat::CorD(long charid,long &kvyhl, double x, double &y, double &z, double &z2)
02010 {
02011
02012 switch (MatChar[charid]){
02013 case 0:{
02014 kvyhl = 0;
02015 break;
02016 }
02017 case 1:{
02018
02019
02020 y = data[charid]->getval(0.0);
02021 kvyhl = 1;
02022 break;
02023 }
02024 case 2:{
02025
02026
02027 y = data[charid]->getval(x);
02028 kvyhl = 2;
02029
02030
02031 long nr=data[charid]->tabf[0].asize;
02032
02033
02034 z2 = data[charid]->tabf[0].y[nr-1];
02035
02036 break;
02037 }
02038 case 30:{
02039 y = MatFunce[charid][0];
02040 z = MatFunce[charid][1];
02041 kvyhl = 30;
02042 break;
02043 }
02044 case 31:{
02045 y = MatFunce[charid][0];
02046 z = MatFunce[charid][1];
02047 if (charid == 6)
02048 {
02049 z2 = MatFunce[charid][2];
02050 }
02051 kvyhl = 31;
02052 break;
02053 }
02054 case 32:{
02055 y = MatFunce[charid][0];
02056 z = MatFunce[charid][1];
02057 kvyhl = 32;
02058 break;
02059 }
02060 case 33:{
02061 kvyhl = 33;
02062
02063 break;
02064 }
02065 default:{
02066 print_err("unknown definition of material parameter is required",__FILE__,__LINE__,__func__);
02067 }
02068 }
02069 }
02070
02071
02072
02073
02074
02075
02076
02077
02078
02079
02080
02081
02082
02083
02084
02085
02086
02087
02088
02089
02090
02091
02092
02093
02094
02095
02096
02097
02098
02099
02100
02101
02102
02103
02104
02105
02106
02107
02108
02109
02110
02111
02112
02113
02114
02115
02116
02117
02118
02119
02120
02121
02122
02123
02124
02125
02126
02127
02128
02129
02130
02131
02132
02133
02134
02135
02136
02137
02138
02139
02140
02141
02142
02143
02144
02145
02146
02147
02148
02149
02150
02151
02152
02153
02154
02155
02156
02157
02158
02159
02160
02161
02162
02163
02164
02165
02166
02167
02168
02169
02170
02171
02172
02173
02174
02175
02176
02177
02178
02179
02180
02181
02182
02183
02184
02185
02186
02187
02188
02189
02190
02191
02192
02193
02194
02195
02196
02197
02198
02199
02200
02201
02202
02203
02204
02205
02206
02207
02208
02209
02210
02211
02212
02213
02214
02215
02216
02217
02218
02219
02220
02221
02222
02223
02224
02225
02226
02227 void kunmat::aux_values (long ipp,double *inv,double *inp,double *ine,double *out)
02228 {
02229 double rh,tk,w,dwdf, kapaA, smc;
02230
02231 switch (kunzeltype){
02232 case 1:{
02233
02234
02235
02236
02237 rh = inv[0];
02238
02239 tk = inv[1];
02240
02241 if (rh>1.0)
02242 rh = 1.0;
02243 if (rh<0.0)
02244 rh = 0.0;
02245
02246
02247 smc = sm.getval (0.0);
02248
02249 w = sorpiso.isotherm_value (rh);
02250
02251 dwdf = sorpiso.derivative_isotherm_value (rh);
02252
02253 if (w>smc)
02254 w = smc;
02255 if (w<0.0)
02256 w = 0.0;
02257
02258
02259 kapaA = kappa.getval (w);
02260
02261
02262 out[0] = w;
02263
02264 out[1] = dwdf;
02265
02266 out[2] = smc;
02267
02268 out[3] = kapaA;
02269
02270 out[4] = kunzeltype;
02271
02272 break;
02273 }
02274
02275
02276
02277
02278
02279
02280
02281
02282
02283
02284
02285
02286
02287
02288
02289
02290
02291
02292
02293
02294
02295
02296
02297
02298
02299
02300
02301
02302
02303
02304
02305
02306
02307
02308
02309
02310
02311
02312
02313
02314
02315
02316
02317
02318
02319
02320
02321
02322
02323
02324
02325
02326
02327
02328
02329
02330
02331
02332
02333
02334
02335
02336
02337
02338
02339
02340
02341 case 8:{
02342
02343
02344 double pv;
02345
02346 pv = inv[0];
02347
02348 tk = inv[1];
02349
02350
02351 smc = sm.getval (0.0);
02352
02353 w = sorpiso.isotherm_value (pv);
02354
02355 dwdf = sorpiso.derivative_isotherm_value (pv);
02356
02357 if (w>smc)
02358 w = smc;
02359 if (w<0.0)
02360 w = 0.0;
02361
02362
02363
02364
02365
02366 out[0] = w;
02367
02368 out[1] = dwdf;
02369
02370 out[2] = smc;
02371
02372 out[3] = 10000.0;
02373
02374 out[4] = kunzeltype;
02375
02376 break;
02377 }
02378 case 11:{
02379
02380
02381
02382 smc = sm.getval (0.0);
02383
02384 hystereze (ipp,inv,inp,ine,out);
02385
02386 out[2] = smc;
02387 out[4] = kunzeltype;
02388
02389 break;
02390 }
02391 default:
02392 print_err("unknown type of Kunzel model is required",__FILE__,__LINE__,__func__);
02393 }
02394 }
02395
02396
02397
02398
02399
02400
02401
02402
02403
02404 void kunmat::save_values (long ipp,double *out)
02405 {
02406 Tm->ip[ipp].eqother[0]=out[0];
02407 Tm->ip[ipp].eqother[1]=out[1];
02408 Tm->ip[ipp].eqother[2]=out[2];
02409 Tm->ip[ipp].eqother[3]=out[3];
02410 Tm->ip[ipp].eqother[4]=out[4];
02411 }
02412
02413
02414
02415
02416
02417
02418
02419
02420
02421
02422 void kunmat::give_values (long ipp,double *av,double *pv,double *eq)
02423 {
02424 av[0] = Tm->ip[ipp].av[0];
02425 av[1] = Tm->ip[ipp].av[1];
02426
02427 pv[0] = Tm->ip[ipp].pv[0];
02428 pv[1] = Tm->ip[ipp].pv[1];
02429
02430 eq[0] = Tm->ip[ipp].eqother[0];
02431 eq[1] = Tm->ip[ipp].eqother[1];
02432 eq[2] = Tm->ip[ipp].eqother[2];
02433 eq[3] = Tm->ip[ipp].eqother[3];
02434 eq[4] = Tm->ip[ipp].eqother[4];
02435 }
02436
02437
02438
02439
02440
02441
02442
02443
02444 void kunmat::initvalues (long ipp,long ido)
02445 {
02446 Tm->ip[ipp].eqother[4] = kunzeltype;
02447 }
02448
02449
02450
02451
02452
02453
02454
02455
02456
02457
02458
02459
02460
02461
02462
02463
02464
02465
02466
02467
02468
02469
02470
02471
02472
02473
02474
02475
02476 void kunmat::hystereze (long ipp,double *inv,double *inp,double *ine,double *out)
02477 {
02478 double rh,rhp,w,wa,wp,wd,dsi,ddsi,xi,xikappa,kapa,dwdf,kapaa,kapap,kapad,dkdwa,dkdwd;
02479 double logkapap,logkapaa,logkapad,logdkdwa,logdkdwd;
02480 double time,act_time,start_time;
02481
02482
02483 act_time = Tp->time;
02484
02485 start_time = Tp->timecont.start_time;
02486
02487 time = act_time - start_time;
02488
02489
02490 rh = inv[0];
02491
02492 if (time < time_hyst){
02493
02494
02495
02496 w = sorpiso.isotherm_value (rh);
02497
02498 kapa = kappa.getval (w);
02499
02500 dwdf = sorpiso.derivative_isotherm_value (rh);
02501
02502 }
02503 else{
02504
02505
02506
02507 wp = ine[0];
02508
02509 rhp = inp[0];
02510
02511
02512 wa = sorpiso.isotherm_value (rh);
02513
02514 wd = desorpiso.isotherm_value (rh);
02515
02516 dsi = sorpiso.derivative_isotherm_value (rh);
02517
02518 ddsi = desorpiso.derivative_isotherm_value (rh);
02519
02520
02521 dwdf = ((wp-wa)*(wp-wa)*adsi*ddsi + (wp-wd)*(wp-wd)*asi*dsi)/((wd-wa)*(wd-wa));
02522
02523
02524 w = wp + xi*(rh-rhp);
02525
02526
02527
02528 kapap = ine[3];
02529
02530 kapaa = kappa.getval (w);
02531
02532 kapad = kappadry.getval (w);
02533
02534 dkdwa = kappa.getderiv (w);
02535
02536 dkdwd = kappadry.getderiv (w);
02537
02538 logkapap = log (kapap);
02539 logkapaa = log (kapaa);
02540 logkapad = log (kapad);
02541 logdkdwa = log (dkdwa);
02542 logdkdwd = log (dkdwd);
02543
02544
02545 xikappa = ((logkapap-logkapaa)*(logkapap-logkapaa)*adkappa*logdkdwd + (logkapap-logkapad)*(logkapap-logkapad)*akappa*logdkdwa)/((logkapad-logkapaa)*(logkapad-logkapaa));
02546
02547
02548
02549
02550 }
02551
02552
02553 out[0] = w;
02554
02555 out[1] = dwdf;
02556
02557 out[3] = kapa;
02558
02559
02560
02561
02562
02563
02564
02565
02566
02567
02568
02569
02570
02571
02572
02573
02574
02575
02576
02577
02578
02579
02580
02581
02582
02583
02584
02585
02586
02587
02588
02589
02590
02591
02592
02593
02594
02595
02596
02597
02598
02599
02600
02601
02602
02603
02604
02605
02606
02607
02608
02609
02610
02611
02612
02613
02614
02615
02616
02617
02618
02619
02620
02621
02622
02623
02624
02625
02626
02627
02628
02629
02630
02631
02632
02633
02634
02635
02636
02637
02638
02639
02640
02641
02642
02643
02644
02645
02646
02647
02648
02649
02650
02651
02652
02653
02654
02655
02656
02657
02658
02659
02660
02661
02662
02663
02664
02665
02666
02667
02668
02669
02670
02671
02672
02673
02674
02675
02676
02677
02678
02679
02680
02681
02682
02683
02684
02685
02686
02687
02688
02689
02690
02691
02692
02693
02694
02695
02696
02697
02698
02699
02700
02701
02702
02703
02704
02705
02706
02707
02708
02709
02710
02711
02712
02713
02714
02715
02716
02717
02718
02719
02720
02721
02722
02723
02724
02725
02726
02727
02728
02729
02730
02731
02732
02733
02734
02735
02736
02737
02738
02739
02740
02741
02742
02743
02744
02745
02746
02747
02748
02749
02750
02751
02752
02753
02754
02755
02756
02757
02758
02759
02760
02761
02762
02763 }
02764
02765
02766
02767
02768
02769
02770
02771
02772
02773
02774
02775
02776
02777 void kunmat::give_dof_names(namevart *dofname, long ntm)
02778 {
02779 if (ntm < 1)
02780 {
02781 print_err("the model defines %ld unknowns but number of transported media is %ld",
02782 __FILE__, __LINE__, __func__, 1, ntm);
02783 abort();
02784 }
02785 dofname[0] = trf_rel_humidity;
02786 dofname[1] = trf_temperature;
02787 }
02788
02789
02790
02791
02792
02793
02794
02795
02796
02797
02798
02799
02800 double kunmat::give_temperature (long ipp)
02801 {
02802 return Tm->ip[ipp].av[1];
02803 }
02804
02805
02806
02807
02808
02809
02810
02811
02812
02813
02814
02815
02816 double kunmat::give_inittemperature (long ipp)
02817 {
02818 print_err("Not yet implemented", __FILE__, __LINE__, __func__);
02819 abort();
02820 }
02821
02822
02823
02824
02825
02826
02827
02828
02829
02830
02831
02832
02833 double kunmat::give_rel_hum (long ipp)
02834 {
02835 return Tm->ip[ipp].av[0];
02836 }
02837
02838
02839
02840
02841
02842
02843
02844
02845
02846
02847
02848
02849 double kunmat::give_vol_moist (long ipp)
02850 {
02851 return Tm->ip[ipp].eqother[0];
02852 }
02853
02854
02855
02856
02857
02858
02859
02860
02861
02862
02863
02864
02865
02866
02867 void kunmat::give_reqntq(long *antq)
02868 {
02869 if (daminfl == on){
02870
02871 antq[scal_iso_damage-1] = 1;
02872
02873 antq[proc_zone_length-1] = 1;
02874
02875 antq[crack_width-1] = 1;
02876 }
02877 }