Add Delta/R contribution

- Needs slightly modification of referrors, but trend stays the same

fenicsR13/solver.py

@@ -519,7 +519,7 @@ class Solver:
                 )
                 - 24/75 * kn * df.div(s_) * df.div(r_)
                 # For Delta-term, works for R13 but fails for heat:
-                # + 4/5 * kn * df.div(s_) * df.div(r_)
+                + 4/5 * cpl * kn * df.div(s_) * df.div(r_)
                 + 4/15 * (1/kn) * df.inner(s_, r_)
             ) * df.dx + (
                 + 1/(2*xi_tilde) * n(s_) * n(r_)

tests/r13/referrors/r13_1_coeffs_nosources_norot_inflow_p1p1p1p1p1_stab/errors.csv

 h,theta_L_2,theta_l_inf,sx_L_2,sx_l_inf,sy_L_2,sy_l_inf,p_L_2,p_l_inf,ux_L_2,ux_l_inf,uy_L_2,uy_l_inf,sigmaxx_L_2,sigmaxx_l_inf,sigmaxy_L_2,sigmaxy_l_inf,sigmayy_L_2,sigmayy_l_inf
-0.9886573325052778,0.0586870540823,0.0676503244787,0.329696368224,0.199785942014,0.309522209551,0.191084853197,0.800387820221,1.48974508465,1.28875169247,0.769191646313,0.596420728898,0.918128560147,0.657181039184,0.508388607105,0.942099171231,0.696258201166,1.0170505251,0.566855115942
-0.6340332990709842,0.0437897459673,0.0558994267328,0.292129271704,0.176316844509,0.191328945359,0.116255847771,0.568763832117,1.08762985584,0.721169499849,0.731448278431,0.483342193529,0.821899359254,0.846780349993,0.375249672323,1.06481525539,0.642386852352,1.33471169629,0.571848421672
-0.32904683851469807,0.0192119603565,0.0324129668358,0.128726432938,0.093244665703,0.076416618118,0.0553658475547,0.183222887021,0.390491659321,0.211253809652,0.434764609749,0.24136899917,0.515274169846,0.436346032961,0.221314372068,0.630026442467,0.371620754375,0.649296254338,0.319473520575
-0.16754966839339377,0.00694051547747,0.0157948502449,0.0304363993486,0.0287934223037,0.0205467775845,0.0160415867278,0.0419775349096,0.106868753733,0.0766925874737,0.174209454888,0.0739343482418,0.212297022033,0.118284867953,0.071854044975,0.165833774216,0.115059880779,0.170685083055,0.0956612639106
-0.08734460120995041,0.00216991291744,0.0075396687322,0.00521800488363,0.00609255611535,0.0046040070456,0.00481737873745,0.006312864438,0.0113328146225,0.0258526595398,0.0283359807048,0.0122057677818,0.0453378900007,0.0202453041047,0.0175683708925,0.0280622231826,0.0189196691116,0.0279502614323,0.0194507271278
+0.9886573325052778,0.0681867538775,0.0744539102566,0.376530880493,0.24884567619,0.363972780949,0.261175535864,0.801010940085,1.4955964945,1.28894577031,0.769154212937,0.596419281717,0.91806487552,0.6584749777,0.505522473004,0.943493976447,0.695665581149,1.01895680449,0.572317265128
+0.6340332990709842,0.0563199044772,0.0664024797673,0.316590229902,0.215754837442,0.235908092371,0.168467984865,0.56913807613,1.09243767221,0.721347887552,0.73148393092,0.483286600964,0.821814217016,0.847254469345,0.375308208233,1.06545502655,0.642082935475,1.33553232386,0.578535273719
+0.32904683851469807,0.0263501516916,0.0424544126873,0.136617355955,0.100549312594,0.0962316643263,0.0945520415061,0.183371466493,0.392314115421,0.211273788406,0.43487341045,0.241364752444,0.515290568846,0.436581878873,0.223060236894,0.630309389112,0.371877961948,0.649563290376,0.322130350885
+0.16754966839339377,0.0102906267322,0.0220800835808,0.034229331923,0.0350341879696,0.0257316468831,0.0233166421871,0.0420244617947,0.107524381474,0.0766733862495,0.174234155942,0.0739334595795,0.21237321596,0.118361130309,0.0723062928699,0.165931947258,0.115090003178,0.170777578084,0.0964224216797
+0.08734460120995041,0.00332107495712,0.0105706355644,0.00629371863111,0.00773578935345,0.0060944254744,0.00632790922008,0.00633021225488,0.0114770503456,0.025844870766,0.0283573560354,0.0122049218696,0.0453076116213,0.0202692438884,0.0174774026492,0.0280936390693,0.0189922287296,0.0279809707499,0.0194754324376

tests/r13/referrors/r13_1_coeffs_sources_rot_noinflow_p1p1p1p1p1_stab/errors.csv

 h,theta_L_2,theta_l_inf,sx_L_2,sx_l_inf,sy_L_2,sy_l_inf,p_L_2,p_l_inf,ux_L_2,ux_l_inf,uy_L_2,uy_l_inf,sigmaxx_L_2,sigmaxx_l_inf,sigmaxy_L_2,sigmaxy_l_inf,sigmayy_L_2,sigmayy_l_inf
-0.9886573325052778,0.16035552156,0.124845108097,0.431581477045,0.280727348439,0.375452001809,0.26823762042,3.43434620181,3.80916463169,0.634754458214,0.696037494878,0.794142372165,0.835325265391,0.397123703791,0.397278399293,0.330409813609,0.366490080602,0.38002082965,0.398994331601
-0.6340332990709842,0.0981425236808,0.0757081274074,0.312881560923,0.179081792245,0.30126728478,0.206039076705,1.92897949058,2.20850722268,0.373575105683,0.652062186827,0.517223581063,0.793872461571,0.236027817072,0.234983180651,0.227370825237,0.240312230842,0.223063649476,0.233430971298
-0.32904683851469807,0.0294465940072,0.0566183297865,0.0940610773307,0.0759668908424,0.0823585227834,0.0577914752096,0.298267114799,0.383599562406,0.190401003388,0.550888314039,0.311831688207,0.74822378413,0.0933013644238,0.090132029723,0.0909167726683,0.0958021911184,0.0897557756587,0.0845623131466
-0.16754966839339377,0.0113750937644,0.0329618204475,0.0228522613663,0.0258030628389,0.0210836152103,0.0237901500964,0.0355675752413,0.0775969119984,0.0840986393902,0.343378651656,0.145026026069,0.477215173129,0.030469079139,0.042346579513,0.0291424062688,0.0369222545308,0.0290337337309,0.0370293192259
-0.08734460120995041,0.00355815829269,0.0149426179682,0.0047393746931,0.00565353285106,0.00483782644175,0.00644081420721,0.00282114157963,0.0102867614987,0.0283054809269,0.144814385734,0.0510019972174,0.230149103304,0.00747478428943,0.0120302297447,0.00703714030599,0.0104552330516,0.00717958701229,0.0112495987517
+0.9886573325052778,0.182351535382,0.14076396052,0.448005026494,0.333048810364,0.373635108049,0.244188750391,3.44561693309,3.82139145711,0.633901565655,0.695535690804,0.793350509306,0.835155259264,0.397412108941,0.397960622392,0.330874986524,0.366342944567,0.380546525621,0.398702952703
+0.6340332990709842,0.133801149839,0.103286607921,0.327144099552,0.21380140044,0.306717019209,0.190397483259,1.93465922204,2.21516348477,0.373113555102,0.651741272663,0.516840870497,0.793701456383,0.236149547516,0.235122742916,0.227566098121,0.240340419061,0.223342718558,0.233263878085
+0.32904683851469807,0.0425610827068,0.0691291340179,0.102795838958,0.0853369124106,0.0968061169972,0.0748375167591,0.298594597354,0.384467541498,0.190404674516,0.550939889058,0.311852999463,0.748365920545,0.093275336168,0.0902217705246,0.0909420799688,0.0958402132679,0.0897348336527,0.0845141677969
+0.16754966839339377,0.0165773552561,0.0420803547265,0.0269617502395,0.0292648191225,0.0247391924075,0.0332452844217,0.035636739436,0.0779090175274,0.084110231319,0.343410338765,0.14502819864,0.477193331598,0.0304649026046,0.042348490009,0.0291474819284,0.0369085150459,0.0290293002466,0.0370417377333
+0.08734460120995041,0.00543472071402,0.0202277761477,0.0057269591288,0.00737481073138,0.00593992934544,0.00883109775742,0.00281806970715,0.0103402010456,0.0283047289264,0.144806489883,0.0510047804751,0.230206532554,0.00747485872226,0.0120282294928,0.0070357642113,0.0104607724696,0.00717957303192,0.0112525092712

tests/r13/referrors/r13_1_coeffs_sources_rot_noinflow_p2p2p2p2p2_stab/errors.csv

 h,theta_L_2,theta_l_inf,sx_L_2,sx_l_inf,sy_L_2,sy_l_inf,p_L_2,p_l_inf,ux_L_2,ux_l_inf,uy_L_2,uy_l_inf,sigmaxx_L_2,sigmaxx_l_inf,sigmaxy_L_2,sigmaxy_l_inf,sigmayy_L_2,sigmayy_l_inf
-0.9886573325052778,0.0496740788862,0.0396821486168,0.203797597346,0.202943920328,0.173253875471,0.179980827012,0.66927242935,0.329803235806,0.314874282264,0.508438646953,0.468689465066,0.588220435812,0.309544263748,0.965443887003,0.267008301967,0.940441696786,0.29862018587,0.938652897221
-0.6340332990709842,0.0148084989825,0.0199122153764,0.0961336112955,0.0982214142727,0.0987467799292,0.105063764991,0.156682147585,0.113687716538,0.11428919551,0.317372947244,0.166845825588,0.346794758727,0.138063227948,0.457540885406,0.134849106181,0.455047465514,0.132906095532,0.440531903568
-0.32904683851469807,0.00336013596909,0.00884253282265,0.0315595158733,0.0384697233064,0.0324785982145,0.0350718023168,0.0195660651934,0.0332944992568,0.0279925920888,0.115727218391,0.0396120591337,0.147670415333,0.0465790672836,0.134122648167,0.045683498084,0.129641822966,0.0452670047363,0.127475547996
-0.16754966839339377,0.00104772582563,0.00356184215138,0.00934247774584,0.0116966428567,0.00942873490841,0.0112787510004,0.00550470207638,0.0130369190563,0.00650116266157,0.0541222266242,0.00906469680812,0.0479083943841,0.0136167591163,0.030468413604,0.0134499886918,0.0281816990923,0.0130249729767,0.0284316989442
+0.9886573325052778,0.0461927691018,0.0344434354208,0.197613772515,0.19275229515,0.167509030969,0.173966240854,0.669547382307,0.330031241428,0.31503782588,0.50874376628,0.468991166778,0.588484810943,0.30964169231,0.965845136824,0.267110740203,0.940886851784,0.298707457316,0.939018437795
+0.6340332990709842,0.0154715541439,0.0219757212353,0.0957342734793,0.0994403436492,0.0982394905764,0.104432748819,0.156715928541,0.113692060823,0.114309344292,0.317480776344,0.166890707415,0.346932390922,0.138065528782,0.457600966042,0.134853451324,0.455112745034,0.132907675389,0.440590849737
+0.32904683851469807,0.00367649110892,0.0106253994445,0.0316157554632,0.0394624981612,0.032710534065,0.0344316992028,0.0195656669564,0.0333103800285,0.0279930836005,0.115719663763,0.0396104611193,0.147651953195,0.0465778517115,0.134133652206,0.0456840036123,0.129647480917,0.0452660401111,0.127484012579
+0.16754966839339377,0.0012171822591,0.00455428138134,0.0093593528868,0.0117920563772,0.00944368709076,0.0114312916128,0.00550464444808,0.0130355892877,0.00650094042838,0.0541210605118,0.0090643315239,0.0479081568585,0.0136166796242,0.0304685124146,0.0134500201023,0.0281813412898,0.0130248923172,0.0284313656633