Commit 930a9e70 authored by Lambert Theisen's avatar Lambert Theisen 🚀

Revert back to positive inflow (+sig_nn) for reason:

1) Fully symmetric system ala:
  |A   B^T| |x| = |f|
  |B       0| |y|    |g|
2) Still have positive-definite system when stabilized in almost anti-sym formulation because:
  => d + f1 + f2 + h > 0 (quadratic term)
parent e2798235
......@@ -583,7 +583,7 @@ class Solver:
(tt(psi_) + (1/2) * nn(psi_))
)
+ (1/chi_tilde) * nt(sigma_) * nt(psi_)
) * df.ds - sum([ # Changed inflow condition => minus
) * df.ds + sum([ # Changed inflow condition => minus
bcs[bc]["epsilon_w"] * chi_tilde * nn(sigma_) * nn(psi_) *
df.ds(bc)
for bc in bcs.keys()
......@@ -639,7 +639,7 @@ class Solver:
lhs[1] = +1*b(kappa, s)+0 +0 +0 +0
lhs[2] = +1*c(s, psi) +0 +d(sigma, psi)-e(u, psi)+f(p, psi)
lhs[3] = +1*0 +0 +e(v, sigma) +0 +g(p, v)
lhs[4] = +1*0 +0 -f(q, sigma) -g(q, u) +h(p, q)
lhs[4] = +1*0 +0 +f(q, sigma) -g(q, u) +h(p, q)
# 2) Right-hand sides:
rhs[0] = - sum([
n(r) * bcs[bc]["theta_w"] * df.ds(bc)
......
......@@ -102,7 +102,7 @@ bcs:
# - relative_errors: Use relative errors. If exact sol. is zero, use absolute.
convergence_study:
enable: True
exact_solution: esols/1_coeffs_nosources_norot_inflow.cpp
exact_solution: esols/1_coeffs_nosources_norot_inflow_positive.cpp
plot: False
write_systemmatrix: False
rescale_pressure: False
......
......@@ -102,7 +102,7 @@ bcs:
# - relative_errors: Use relative errors. If exact sol. is zero, use absolute.
convergence_study:
enable: True
exact_solution: esols/1_coeffs_nosources_norot_inflow.cpp
exact_solution: esols/1_coeffs_nosources_norot_inflow_positive.cpp
plot: False
write_systemmatrix: False
rescale_pressure: False
......
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.0749854935308,0.0787613855322,0.508527345579,0.343945391205,0.370708290954,0.257777630188,1.50582801652,2.76383201873,0.697087612134,0.607698803914,0.687752348432,0.862015454556,1.21709005197,0.592311682958,1.98517604408,1.27582426608,2.24662155114,1.11030513218
0.6340332990709842,0.0763355660647,0.0701057613416,0.916193173135,0.785916331974,0.618094593262,0.323854879025,3.6513380924,5.49636989361,5.49013492407,2.28584893436,2.15378978384,1.54260735924,5.62586360461,3.49915548567,6.6525810967,3.10873999144,8.46182462893,4.98845940255
0.32904683851469807,0.0261031827137,0.0468979908421,0.0907310150632,0.0846420173058,0.0952112618586,0.113078818103,0.234634421851,0.274321184496,0.94151114893,0.670842834426,0.343243144001,0.720540819098,0.463283077117,0.493052871279,0.572795386851,0.587533226053,0.51950745745,0.501755329362
0.16754966839339377,0.0102785848201,0.0240550863362,0.0267326911862,0.0303909118002,0.0257865756701,0.0267569797033,0.0717261830671,0.122041120931,0.230837129,0.280498266541,0.10060831498,0.314397821549,0.112691328922,0.135349547413,0.13702471256,0.172639836751,0.128410213898,0.13352358799
0.08734460120995041,0.00334923619914,0.0113326215569,0.00724853402965,0.0100002493692,0.00717263422709,0.00822905805696,0.0226609861373,0.0304988807026,0.0532194851553,0.0643177706557,0.0220932195036,0.0812448628194,0.0370049830733,0.0307155223862,0.0462839205587,0.0471886097176,0.0439334788056,0.0324014370782
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
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.0143768913243,0.0167763774855,0.154427796117,0.212905569859,0.170021307007,0.189184423016,0.222665842777,0.38785328282,0.946625073133,0.567225546841,0.310747128191,0.414217805868,0.376671752432,0.462493535845,0.833076883105,0.836036693911,0.546864985397,0.607810788706
0.6340332990709842,0.0207579137652,0.0367173799369,0.169037418712,0.248093709767,0.163463969828,0.130693891918,0.665027312211,0.8072723145,2.01920543274,0.91507796314,0.701235155631,0.781210911785,1.47078769847,1.04351764263,1.49146487718,1.08353617431,2.00889778669,1.3755469062
0.32904683851469807,0.0130451011078,0.0214463940876,0.11719988638,0.135578567567,0.0945258805968,0.0635837546267,0.478092455792,0.458272915175,1.30107452488,0.646522741078,0.454477192923,0.649877565913,0.875180306235,0.718588084973,1.07272217111,0.790755661626,1.03959737869,0.753488608423
0.16754966839339377,0.00244245700739,0.00484554351159,0.0224211862256,0.0282299306717,0.0191812627504,0.0141904437469,0.0888811739689,0.0840260609066,0.220940884681,0.113340275219,0.0826896049548,0.151911670261,0.155970621267,0.128324961683,0.195690972003,0.153115738912,0.190054619914,0.139421512179
0.08734460120995041,0.000595068041661,0.00121440933827,0.0055147176246,0.00700562569891,0.00471691302173,0.00347684971707,0.0217752311678,0.0216173663209,0.0531285069094,0.0331765062864,0.0204957947248,0.0441530177482,0.0379894104042,0.0316570167234,0.0479717016566,0.0387767133382,0.0457767021787,0.0339001090382
0.9886573325052778,0.0129233871354,0.0122078711717,0.142606796966,0.181065193865,0.171232198277,0.183332375248,0.153825868414,0.251222821161,0.467199059836,0.351578955684,0.230834001404,0.320524042489,0.208790031071,0.243407603825,0.319992672675,0.381318318486,0.290442860302,0.293072931386
0.6340332990709842,0.00895852938885,0.0123205599413,0.0888620443435,0.119606329269,0.0939259762802,0.109788786856,0.113074392385,0.197491883103,0.371038603195,0.243923106839,0.147943964597,0.182169123232,0.246500049281,0.170876040277,0.294962464201,0.277948273655,0.340638217095,0.233412652185
0.32904683851469807,0.00410802891336,0.00934057522246,0.0379399997286,0.0455212061874,0.0335501009603,0.0359811705404,0.0540253959852,0.0772723743642,0.246188980709,0.127028969543,0.081590098508,0.131052472565,0.166863682621,0.113594364991,0.230218538928,0.160644473446,0.224852901943,0.128644181635
0.16754966839339377,0.00144444138091,0.00370496500738,0.0135958948481,0.0154303118198,0.0102535216357,0.00986740957961,0.021418211269,0.0247389660103,0.089924516004,0.056207161841,0.0312738073508,0.064344793798,0.0680554247232,0.054310404741,0.0930171192734,0.0500502953316,0.0923992386119,0.0519505190851
0.08734460120995041,0.000397238867704,0.000821517015023,0.00379286200302,0.00419956890815,0.002671688347,0.00241044530682,0.00616089809359,0.00768745653881,0.0249099236635,0.0231595940585,0.00908472799792,0.0230123785227,0.0197002729691,0.0137748045052,0.0271022278468,0.0121428020021,0.0268543719238,0.0143682011848
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment