00001 #include "globalt.h"
00002 #include "saltmat1.h"
00003 #include "globmatt.h"
00004
00005 dampermeability saltmat1::damper;
00006
00007 saltmat1::saltmat1 ()
00008 {
00009
00010 daminfl = off;
00011 }
00012
00013 saltmat1::~saltmat1 ()
00014 {}
00015
00016
00017
00018
00019
00020
00021
00022
00023 void saltmat1::read (XFILE *in)
00024 {
00025
00026 kappa.read (in);
00027
00028 sm.read (in);
00029
00030 dcoef.read (in);
00031
00032 bindiso.read (in);
00033 }
00034
00035
00036
00037
00038
00039
00040
00041
00042 void saltmat1::print (FILE *out)
00043 {
00044
00045 kappa.print (out);
00046
00047 sm.print (out);
00048
00049 dcoef.print (out);
00050
00051 bindiso.print (out);
00052 }
00053
00054
00055
00056
00057
00058
00059
00060
00061
00062
00063
00064
00065
00066 void saltmat1::matcond (matrix &d,long ri,long ci,long ipp)
00067 {
00068 long n;
00069 n = d.n;
00070
00071 switch (n){
00072 case 1:{
00073 matcond1d (d,ri,ci,ipp);
00074 break;
00075 }
00076 case 2:{
00077 matcond2d (d,ri,ci,ipp);
00078 break;
00079 }
00080 case 3:{
00081 matcond3d (d,ri,ci,ipp);
00082 break;
00083 }
00084 default:{
00085 print_err("unknown number of space dimension is required",__FILE__,__LINE__,__func__);
00086 }
00087 }
00088
00089 if (daminfl == on){
00090
00091 damper.matcond (d,ipp);
00092 }
00093 Tm->ip[ipp].eqother[0] = d[0][0];
00094 }
00095
00096
00097
00098
00099
00100
00101
00102
00103
00104
00105 void saltmat1::matcond1d (matrix &d,long ri,long ci,long ipp)
00106 {
00107 double k;
00108 double w,cf;
00109 k = 0.0;
00110
00111 w = Tm->ip[ipp].av[0];
00112 cf = Tm->ip[ipp].av[1];
00113
00114 if((ri == 0) && (ci == 0))
00115 k = k11 (w);
00116 if((ri == 0) && (ci == 1))
00117 k = k12 ();
00118 if((ri == 1) && (ci == 0))
00119 k = k21 (w,cf);
00120 if((ri == 1) && (ci == 1))
00121 k = k22 (w,cf);
00122
00123 d[0][0] = k;
00124 }
00125
00126
00127
00128
00129
00130
00131
00132
00133
00134 void saltmat1::matcond2d (matrix &d,long ri,long ci,long ipp)
00135 {
00136 double k;
00137 double w,cf;
00138 k = 0.0;
00139
00140 w = Tm->ip[ipp].av[0];
00141 cf = Tm->ip[ipp].av[1];
00142
00143 if((ri == 0) && (ci == 0))
00144 k = k11 (w);
00145 if((ri == 0) && (ci == 1))
00146 k = k12 ();
00147 if((ri == 1) && (ci == 0))
00148 k = k21 (w,cf);
00149 if((ri == 1) && (ci == 1))
00150 k = k22 (w,cf);
00151
00152 fillm(0.0,d);
00153
00154 d[0][0] = k; d[0][1] = 0.0;
00155 d[1][0] = 0.0; d[1][1] = k;
00156 }
00157
00158
00159
00160
00161
00162
00163
00164
00165
00166
00167 void saltmat1::matcond3d (matrix &d,long ri,long ci,long ipp)
00168 {
00169 double k;
00170 double w,cf;
00171 k = 0.0;
00172
00173 w = Tm->ip[ipp].av[0];
00174 cf = Tm->ip[ipp].av[1];
00175
00176 if((ri == 0) && (ci == 0))
00177 k = k11 (w);
00178 if((ri == 0) && (ci == 1))
00179 k = k12 ();
00180 if((ri == 1) && (ci == 0))
00181 k = k21 (w,cf);
00182 if((ri == 1) && (ci == 1))
00183 k = k22 (w,cf);
00184
00185 fillm(0.0,d);
00186
00187 d[0][0]=k; d[0][1]=0.0; d[0][2]=0.0;
00188 d[1][0]=0.0; d[1][1]=k; d[1][2]=0.0;
00189 d[2][0]=0.0; d[2][1]=0.0; d[2][2]=k;
00190 }
00191
00192
00193
00194
00195
00196
00197
00198
00199
00200
00201 void saltmat1::matcap (double &c,long ri,long ci,long ipp)
00202 {
00203 double w,cf;
00204
00205 if((ri == 0) && (ci == 0))
00206 c = c11 ();
00207 if((ri == 0) && (ci == 1))
00208 c = c12 ();
00209 if((ri == 1) && (ci == 0)){
00210 cf = Tm->ip[ipp].av[1];
00211 c = c21 (cf);
00212 }
00213 if((ri == 1) && (ci == 1)){
00214 w = Tm->ip[ipp].av[0];
00215 cf = Tm->ip[ipp].av[1];
00216 c = c22 (w,cf);
00217 }
00218 }
00219
00220
00221
00222
00223
00224
00225 double saltmat1::c11 ()
00226 {
00227 double c11;
00228 c11 = 1.0;
00229 return c11;
00230 }
00231
00232 double saltmat1::c12 ()
00233 {
00234 double c12;
00235 c12 = 0.0;
00236 return c12;
00237 }
00238
00239 double saltmat1::c21 (double cf)
00240 {
00241 double c21;
00242 c21 = cf;
00243 return c21;
00244 }
00245
00246 double saltmat1::c22 (double w,double cf)
00247 {
00248 double c22,dCbDcf;
00249
00250 if (cf <= 0.0)
00251 cf = 1.0e-10;
00252 if (w < 0.0)
00253 w = 0.0;
00254
00255 dCbDcf = bindiso.derivative_isotherm_value (cf);
00256
00257 c22 = w + dCbDcf;
00258
00259 return c22;
00260 }
00261
00262
00263
00264 double saltmat1::k11 (double w)
00265 {
00266 double k11,wsat;
00267
00268 if (w < 0.0)
00269 w = 0.0;
00270
00271 wsat = sm.getval (w);
00272 if (w > wsat)
00273 w = wsat;
00274
00275 k11 = kappa.getval (w);
00276
00277 return k11;
00278 }
00279
00280 double saltmat1::k12 ()
00281 {
00282 double k12;
00283
00284 k12 = 0.0;
00285
00286 return k12;
00287 }
00288
00289 double saltmat1::k21 (double w,double cf)
00290 {
00291 double k21,wsat;
00292
00293 if (w < 0.0)
00294 w = 0.0;
00295 wsat = sm.getval (w);
00296 if (w > wsat)
00297 w = wsat;
00298
00299 k21 = w*dcoef.getval (cf);
00300
00301 return k21;
00302 }
00303
00304 double saltmat1::k22 (double w,double cf)
00305 {
00306 double k22,wsat;
00307
00308 if (w < 0.0)
00309 w = 0.0;
00310 wsat = sm.getval (w);
00311 if (w > wsat)
00312 w = wsat;
00313
00314
00315 k22 = cf*kappa.getval (w);
00316
00317 return k22;
00318 }
00319
00320
00321
00322
00323
00324
00325
00326
00327
00328
00329
00330
00331
00332
00333
00334
00335
00336
00337
00338 double saltmat1::transmission_transcoeff (double trc,long ri,long ci,long nid,long bc)
00339 {
00340 double new_trc,w,cf;
00341 new_trc = 0.0;
00342
00343
00344 w = nodalval (0,nid,0);
00345
00346 cf = nodalval (0,nid,1);
00347
00348 if((ri == 0) && (ci == 0))
00349 new_trc = trc;
00350 if((ri == 0) && (ci == 1))
00351 new_trc = 0.0;
00352
00353 if((ri == 1) && (ci == 0))
00354 new_trc = 0.0;
00355 if((ri == 1) && (ci == 1))
00356 new_trc = trc;
00357
00358 return new_trc;
00359 }
00360
00361
00362
00363
00364
00365
00366
00367
00368
00369
00370
00371
00372
00373
00374
00375
00376
00377
00378
00379
00380
00381
00382
00383
00384
00385
00386
00387
00388
00389
00390
00391
00392
00393
00394
00395
00396
00397
00398
00399
00400
00401
00402 double saltmat1::transmission_nodval (double nodval,long ri,long ci,long nid,long bc)
00403 {
00404 double new_nodval,w,cf;
00405 new_nodval = 0.0;
00406
00407
00408 w = nodalval (0,nid,0);
00409
00410 cf = nodalval (0,nid,1);
00411
00412 if((ri == 0) && (ci == 0))
00413
00414 new_nodval = nodval;
00415 if((ri == 0) && (ci == 1))
00416 new_nodval = 0.0;
00417
00418
00419 if((ri == 1) && (ci == 0))
00420 new_nodval = 0.0;
00421 if((ri == 1) && (ci == 1))
00422 new_nodval = nodval;
00423
00424 return new_nodval;
00425 }
00426
00427
00428
00429
00430
00431
00432
00433
00434
00435
00436
00437
00438
00439
00440
00441
00442
00443
00444
00445
00446
00447
00448
00449
00450
00451
00452
00453
00454
00455
00456
00457
00458
00459
00460
00461
00462
00463
00464
00465 double saltmat1::transmission_flux (double nodval,long ri,long ci,long nid,long bc)
00466 {
00467 double flux,w,cf;
00468 flux = 0.0;
00469
00470
00471 w = nodalval (0,nid,0);
00472
00473 cf = nodalval (0,nid,1);
00474
00475
00476 if((ri == 0) && (ci == 0))
00477 flux = get_transmission_flux_ww (nodval,w,bc);
00478 if((ri == 0) && (ci == 1))
00479 flux = 0.0;
00480
00481 if((ri == 1) && (ci == 0))
00482 flux = 0.0;
00483 if((ri == 1) && (ci == 1))
00484 flux = 0.0;
00485
00486 return (flux);
00487 }
00488
00489
00490
00491
00492
00493
00494
00495
00496
00497 double saltmat1::get_transmission_flux_ww (double bv,double w,long bc)
00498 {
00499 double flux;
00500
00501 switch (bc){
00502 case 30:{
00503 flux = (bv - w);
00504 break;
00505 }
00506 default:{
00507 print_err("no acceptable boundary condition is prescribed",__FILE__,__LINE__,__func__);
00508 exit(0);
00509 }
00510 }
00511
00512 return(flux);
00513 }
00514
00515
00516
00517
00518
00519
00520
00521
00522
00523
00524
00525
00526
00527
00528 double saltmat1::get_othervalue(long compother,long ipp, double x1,double x2,double x3)
00529 {
00530 double other;
00531
00532 switch (compother){
00533 case 0:{
00534 other = x1;
00535 break;
00536 }
00537 default:{
00538 fprintf (stderr,"\n\n unknown type of component is required in function (%s, line %d).\n",__FILE__,__LINE__);
00539 }
00540 }
00541 return (other);
00542 }
00543
00544
00545
00546
00547
00548
00549 void saltmat1::print_othervalue_name(FILE *out,long compother)
00550 {
00551 switch (compother){
00552 case 0:{
00553 fprintf (out,"First unknown (Units) ");
00554 break;
00555 }
00556 default:{
00557 fprintf (stderr,"\n\n unknown type of component is required in function (%s, line %d).\n",__FILE__,__LINE__);
00558 }
00559 }
00560 }
00561
00562
00563
00564
00565
00566
00567
00568
00569
00570
00571
00572
00573 void saltmat1::give_dof_names (namevart *dofname, long ntm)
00574 {
00575 if (ntm < 1){
00576 print_err("the model defines %ld unknowns but number of transported media is %ld",
00577 __FILE__, __LINE__, __func__, 1, ntm);
00578 abort();
00579 }
00580 dofname[0] = trf_moisture;
00581 dofname[1] = trf_salt_conc;
00582 }