### Modify inflow BC

```- Only affect R13 inflow test case (adjusted referrors, qualitative convergence remains)
- System is now fully antisymmetric```
parent 76a63aab
 ... ... @@ -583,7 +583,7 @@ class Solver: (tt(psi_) + (1/2) * nn(psi_)) ) + (1/chi_tilde) * nt(sigma_) * nt(psi_) ) * df.ds + sum([ # TODO: Fix inflow with 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() ... ... @@ -634,12 +634,12 @@ class Solver: lhs = [None] * 5 rhs = [None] * 5 # 1) Left-hand sides # TODO: Fix inflow with minus in lhs4 before f # Changed inflow condition => minus before f(q, sigma) lhs = +1*a(s, r) -b(theta, r)-c(r, sigma) +0 +0 lhs = +1*b(kappa, s)+0 +0 +0 +0 lhs = +1*c(s, psi) +0 +d(sigma, psi)-e(u, psi)+f(p, psi) lhs = +1*0 +0 +e(v, sigma) +0 +g(p, v) lhs = +1*0 +0 +f(q, sigma) -g(q, u) +h(p, q) lhs = +1*0 +0 -f(q, sigma) -g(q, u) +h(p, q) # 2) Right-hand sides: rhs = - sum([ n(r) * bcs[bc]["theta_w"] * df.ds(bc) ... ...
 ... ... @@ -13,114 +13,136 @@ namespace py = pybind11; double kn = 1.0; // // Constants eps=10^-5 // double C_0 = 1.088418574133493; // double C_1 = -0.4708784837854026; // double C_2 = -0.1060780078877844; // double C_3 = 0.6397363738946663; // double C_4 = -0.04268373030990083; // double C_5 = 0.1027291801998076; // double C_6 = -0.1867020667903523; // double C_7 = 2.864923360570221e-7; // double C_8 = 1.963244537638472; // double C_9 = 0.007704644954804194; // double CK_1 = -0.01841555365364274; // double CK_2 = 0.02647436311982035; // double CK_3 = -0.0722116290730638; // double CK_4 = 0.01263372155827893; // double CK_5 = 0.03846910453592139; // double CK_6 = 0.003998659067908465; // double CK_7 = 0.05469703377538136; // double CK_8 = 0.04448263813229264; // double CK_9 = -0.001607567940383826; // double CK_10 = 0.1781491089766695; // double C_0 = 1.337436195119152; // double C_1 = -0.6064772869609508; // double C_2 = -0.1726980022374072; // double C_3 = 0.7901754852718373; // double C_4 = -0.08396168002205009; // double C_5 = 0.1252124495285327; // double C_6 = -0.1867020765534401; // double C_7 = 4.568817251907919e-7; // double C_8 = 1.963244517881108; // double C_9 = -0.03843441110894395; // double CK_1 = -0.03180963566002944; // double CK_2 = 0.03411564204278963; // double CK_3 = -0.07221148250544122; // double CK_4 = 0.01263372282796629; // double CK_5 = 0.03846898363886858; // double CK_6 = 0.003998656953005047; // double CK_7 = 0.1003792356748235; // double CK_8 = 0.05030863773067603; // double CK_9 = -0.001651068357981648; // double CK_10 = 0.3495787712905514; // // Constants eps=10^-4 // double C_0 = 1.088458465830652; // double C_1 = -0.4708795265589713; // double C_2 = -0.1060844454754048; // double C_3 = 0.6397268622237638; // double C_4 = -0.0426888328939176; // double C_5 = 0.1027285697835575; // double C_6 = -0.1867022145254646; // double C_7 = 2.864825757859002e-6; // double C_8 = 1.963244238669897; // double C_9 = 0.007703775736883896; // double CK_1 = -0.01841695056397639; // double CK_2 = 0.02647443224729624; // double CK_3 = -0.07220941121077909; // double CK_4 = 0.01263374077119675; // double CK_5 = 0.03846727512087732; // double CK_6 = 0.003998627065175332; // double CK_7 = 0.0547022288973836; // double CK_8 = 0.04447888927384885; // double CK_9 = -0.001604798505932295; // double CK_10 = 0.178170030606713; // double C_0 = 1.337429881309671; // double C_1 = -0.6064620302518795; // double C_2 = -0.17269389293132; // double C_3 = 0.7901524616369147; // double C_4 = -0.0839597704511531; // double C_5 = 0.1252094721825813; // double C_6 = -0.1867023121536032; // double C_7 = 4.568671844035206e-6; // double C_8 = 1.963244041101796; // double C_9 = -0.03843414811236764; // double CK_1 = -0.03180886812453335; // double CK_2 = 0.03411478821870104; // double CK_3 = -0.0722079455756748; // double CK_4 = 0.01263375346771417; // double CK_5 = 0.03846606618426868; // double CK_6 = 0.003998605916734596; // double CK_7 = 0.1003768640953838; // double CK_8 = 0.05030581411268845; // double CK_9 = -0.001649480514012068; // double CK_10 = 0.3495706866039064; // // Constants eps=10^-3 double C_0 = 1.088856712167725; double C_1 = -0.4708896102041703; double C_2 = -0.1061486462553457; double C_3 = 0.6396313746480878; double C_4 = -0.04273974910778836; double C_5 = 0.102722409367338; double C_6 = -0.1867036913227729; double C_7 = 0.00002863849458707334; double C_8 = 1.963241250104882; double C_9 = 0.007695099574563038; double CK_1 = -0.01843088435918791; double CK_2 = 0.02647510412125701; double CK_3 = -0.07218724090202373; double CK_4 = 0.01263393282835227; double CK_5 = 0.03844898782835787; double CK_6 = 0.003998307157812554; double CK_7 = 0.05475405925519971; double CK_8 = 0.04444138977514516; double CK_9 = -0.001577106876225382; double CK_10 = 0.178378791898016; double C_0 = 1.337366717236668; double C_1 = -0.6063094894924456; double C_2 = -0.1726528012092394; double C_3 = 0.7899222753828902; double C_4 = -0.08394067339153539; double C_5 = 0.1251797046222099; double C_6 = -0.1867046673301068; double C_7 = 0.00004567217257213145; double C_8 = 1.963239274978468; double C_9 = -0.03843149872072692; double CK_1 = -0.0318011928378749; double CK_2 = 0.03410625144135048; double CK_3 = -0.07217258866517415; double CK_4 = 0.012634059757885; double CK_5 = 0.03843690185588499; double CK_6 = 0.003998095732771025; double CK_7 = 0.1003531476907441; double CK_8 = 0.05027758648967144; double CK_9 = -0.00163360750417993; double CK_10 = 0.3494898346560165; // // Constants eps=10^-2 // double C_0 = 1.09277278324226; // double C_1 = -0.4709565408029582; // double C_2 = -0.1067733521869535; // double C_3 = 0.6386400415463009; // double C_4 = -0.04323806992636355; // double C_5 = 0.1026552684836559; // double C_6 = -0.1867184037607001; // double C_7 = 0.0002854059606868768; // double C_8 = 1.963211476840119; // double C_9 = 0.007609935703161957; // double CK_1 = -0.0185667325855968; // double CK_2 = 0.02647991106816222; // double CK_3 = -0.07196637153183492; // double CK_4 = 0.0126358461775725; // double CK_5 = 0.03826680259919404; // double CK_6 = 0.003995120114343059; // double CK_7 = 0.05526041387534725; // double CK_8 = 0.04406535024733006; // double CK_9 = -0.001300480668267124; // double CK_10 = 0.180421409499547; // double C_0 = 1.336732533609719; // double C_1 = -0.6047867380273104; // double C_2 = -0.1722420311465265; // double C_3 = 0.7876254443155443; // double C_4 = -0.08374957633036394; // double C_5 = 0.1248826224847917; // double C_6 = -0.1867281362946452; // double C_7 = 0.0004552621117990246; // double C_8 = 1.963191781306897; // double C_9 = -0.03840306975572652; // double CK_1 = -0.03172444950554837; // double CK_2 = 0.03402103132100606; // double CK_3 = -0.07182026259642857; // double CK_4 = 0.01263711189140766; // double CK_5 = 0.03814628389686272; // double CK_6 = 0.003993011829577151; // double CK_7 = 0.1001159319865724; // double CK_8 = 0.04999616566544083; // double CK_9 = -0.001475419473718654; // double CK_10 = 0.3486808426586021; // // Constants eps=10^-1 // double C_0 = 1.125921604594353; // double C_1 = -0.4686941217654443; // double C_2 = -0.1114820681020402; // double C_3 = 0.6256505450040041; // double C_4 = -0.04724969749859539; // double C_5 = 0.1015107112011764; // double C_6 = -0.1868598401676765; // double C_7 = 0.002753811809744579; // double C_8 = 1.962925254827631; // double C_9 = 0.006913904797024802; // double CK_1 = -0.01961423092232162; // double CK_2 = 0.0263626002451694; // double CK_3 = -0.06984306808715408; // double CK_4 = 0.01265423995026528; // double CK_5 = 0.03651538489775871; // double CK_6 = 0.003964481821802226; // double CK_7 = 0.05925618682225351; // double CK_8 = 0.04024172550289172; // double CK_9 = 0.001420920976733996; // double CK_10 = 0.1968170158183649; // double C_0 = 1.330185932037747; // double C_1 = -0.5898443496048875; // double C_2 = -0.1681607691585836; // double C_3 = 0.7651780840964286; // double C_4 = -0.08183362422804108; // double C_5 = 0.1219741008236003; // double C_6 = -0.1869542885544411; // double C_7 = 0.004402163569437981; // double C_8 = 1.962734121526709; // double C_9 = -0.0379330049303497; // double CK_1 = -0.03096036349825703; // double CK_2 = 0.03318469878613935; // double CK_3 = -0.06842516876936609; // double CK_4 = 0.01266652294157217; // double CK_5 = 0.03534582339667038; // double CK_6 = 0.003944022186510511; // double CK_7 = 0.09774692911694641; // double CK_8 = 0.04726696889230132; // double CK_9 = 0.0000532462109035967; // double CK_10 = 0.3405754068324816; // // Constants eps=10^-0 // double C_0 = 1.26608703305865; // double C_1 = -0.4723447945692413; // double C_2 = -0.1340987755805701; // double C_3 = 0.5922083059862845; // double C_4 = -0.06517741879346212; // double C_5 = 0.09936322021509666; // double C_6 = -0.1883345440537925; // double C_7 = 0.02849094538051043; // double C_8 = 1.959940926229848; // double C_9 = -0.02236820885608115; // double CK_1 = -0.02452200335350193; // double CK_2 = 0.02660474893156432; // double CK_3 = -0.04770418667594679; // double CK_4 = 0.01284602485664462; // double CK_5 = 0.01825401538380585; // double CK_6 = 0.003645027941151212; // double CK_7 = 0.07750866083895219; // double CK_8 = 0.02706654817009965; // double CK_9 = 0.01115352331657436; // double CK_10 = 0.2703240873448836; double lambda_1 = sqrt(5.0/9.0); double lambda_2 = sqrt(5.0/6.0); ... ...
 ... ... @@ -13,117 +13,114 @@ namespace py = pybind11; double kn = 1.0; // // Constants eps=10^-5 // double C_0 = 1.337436195119152; // double C_1 = -0.6064772869609508; // double C_2 = -0.1726980022374072; // double C_3 = 0.7901754852718373; // double C_4 = -0.08396168002205009; // double C_5 = 0.1252124495285327; // double C_6 = -0.1867020765534401; // double C_7 = 4.568817251907919e-7; // double C_8 = 1.963244517881108; // double C_9 = -0.03843441110894395; // double CK_1 = -0.03180963566002944; // double CK_2 = 0.03411564204278963; // double CK_3 = -0.07221148250544122; // double CK_4 = 0.01263372282796629; // double CK_5 = 0.03846898363886858; // double CK_6 = 0.003998656953005047; // double CK_7 = 0.1003792356748235; // double CK_8 = 0.05030863773067603; // double CK_9 = -0.001651068357981648; // double CK_10 = 0.3495787712905514; // double C_0 = 1.088418574133493; // double C_1 = -0.4708784837854026; // double C_2 = -0.1060780078877844; // double C_3 = 0.6397363738946663; // double C_4 = -0.04268373030990083; // double C_5 = 0.1027291801998076; // double C_6 = -0.1867020667903523; // double C_7 = 2.864923360570221e-7; // double C_8 = 1.963244537638472; // double C_9 = 0.007704644954804194; // double CK_1 = -0.01841555365364274; // double CK_2 = 0.02647436311982035; // double CK_3 = -0.0722116290730638; // double CK_4 = 0.01263372155827893; // double CK_5 = 0.03846910453592139; // double CK_6 = 0.003998659067908465; // double CK_7 = 0.05469703377538136; // double CK_8 = 0.04448263813229264; // double CK_9 = -0.001607567940383826; // double CK_10 = 0.1781491089766695; // // Constants eps=10^-4 // ... // double C_0 = 1.088458465830652; // double C_1 = -0.4708795265589713; // double C_2 = -0.1060844454754048; // double C_3 = 0.6397268622237638; // double C_4 = -0.0426888328939176; // double C_5 = 0.1027285697835575; // double C_6 = -0.1867022145254646; // double C_7 = 2.864825757859002e-6; // double C_8 = 1.963244238669897; // double C_9 = 0.007703775736883896; // double CK_1 = -0.01841695056397639; // double CK_2 = 0.02647443224729624; // double CK_3 = -0.07220941121077909; // double CK_4 = 0.01263374077119675; // double CK_5 = 0.03846727512087732; // double CK_6 = 0.003998627065175332; // double CK_7 = 0.0547022288973836; // double CK_8 = 0.04447888927384885; // double CK_9 = -0.001604798505932295; // double CK_10 = 0.178170030606713; // // Constants eps=10^-3 double C_0 = 1.337366717236668; double C_1 = -0.6063094894924456; double C_2 = -0.1726528012092394; double C_3 = 0.7899222753828902; double C_4 = -0.08394067339153539; double C_5 = 0.1251797046222099; double C_6 = -0.1867046673301068; double C_7 = 0.00004567217257213145; double C_8 = 1.963239274978468; double C_9 = -0.03843149872072692; double CK_1 = -0.0318011928378749; double CK_2 = 0.03410625144135048; double CK_3 = -0.07217258866517415; double CK_4 = 0.012634059757885; double CK_5 = 0.03843690185588499; double CK_6 = 0.003998095732771025; double CK_7 = 0.1003531476907441; double CK_8 = 0.05027758648967144; double CK_9 = -0.00163360750417993; double CK_10 = 0.3494898346560165; double C_0 = 1.088856712167725; double C_1 = -0.4708896102041703; double C_2 = -0.1061486462553457; double C_3 = 0.6396313746480878; double C_4 = -0.04273974910778836; double C_5 = 0.102722409367338; double C_6 = -0.1867036913227729; double C_7 = 0.00002863849458707334; double C_8 = 1.963241250104882; double C_9 = 0.007695099574563038; double CK_1 = -0.01843088435918791; double CK_2 = 0.02647510412125701; double CK_3 = -0.07218724090202373; double CK_4 = 0.01263393282835227; double CK_5 = 0.03844898782835787; double CK_6 = 0.003998307157812554; double CK_7 = 0.05475405925519971; double CK_8 = 0.04444138977514516; double CK_9 = -0.001577106876225382; double CK_10 = 0.178378791898016; // // Constants eps=10^-2 // double C_0 = 1.336732533609719; // double C_1 = -0.6047867380273104; // double C_2 = -0.1722420311465265; // double C_3 = 0.7876254443155443; // double C_4 = -0.08374957633036394; // double C_5 = 0.1248826224847917; // double C_6 = -0.1867281362946452; // double C_7 = 0.0004552621117990246; // double C_8 = 1.963191781306897; // double C_9 = -0.03840306975572652; // double CK_1 = -0.03172444950554837; // double CK_2 = 0.03402103132100606; // double CK_3 = -0.07182026259642857; // double CK_4 = 0.01263711189140766; // double CK_5 = 0.03814628389686272; // double CK_6 = 0.003993011829577151; // double CK_7 = 0.1001159319865724; // double CK_8 = 0.04999616566544083; // double CK_9 = -0.001475419473718654; // double CK_10 = 0.3486808426586021; // double C_0 = 1.09277278324226; // double C_1 = -0.4709565408029582; // double C_2 = -0.1067733521869535; // double C_3 = 0.6386400415463009; // double C_4 = -0.04323806992636355; // double C_5 = 0.1026552684836559; // double C_6 = -0.1867184037607001; // double C_7 = 0.0002854059606868768; // double C_8 = 1.963211476840119; // double C_9 = 0.007609935703161957; // double CK_1 = -0.0185667325855968; // double CK_2 = 0.02647991106816222; // double CK_3 = -0.07196637153183492; // double CK_4 = 0.0126358461775725; // double CK_5 = 0.03826680259919404; // double CK_6 = 0.003995120114343059; // double CK_7 = 0.05526041387534725; // double CK_8 = 0.04406535024733006; // double CK_9 = -0.001300480668267124; // double CK_10 = 0.180421409499547; // // Constants eps=10^-1 // double C_0 = 1.330185932037747; // double C_1 = -0.5898443496048875; // double C_2 = -0.1681607691585836; // double C_3 = 0.7651780840964286; // double C_4 = -0.08183362422804108; // double C_5 = 0.1219741008236003; // double C_6 = -0.1869542885544411; // double C_7 = 0.004402163569437981; // double C_8 = 1.962734121526709; // double C_9 = -0.0379330049303497; // double CK_1 = -0.03096036349825703; // double CK_2 = 0.03318469878613935; // double CK_3 = -0.06842516876936609; // double CK_4 = 0.01266652294157217; // double CK_5 = 0.03534582339667038; // double CK_6 = 0.003944022186510511; // double CK_7 = 0.09774692911694641; // double CK_8 = 0.04726696889230132; // double CK_9 = 0.0000532462109035967; // double CK_10 = 0.3405754068324816; // // Constants eps=10^-0 // double C_0 = 1.26608703305865; // double C_1 = -0.4723447945692413; // double C_2 = -0.1340987755805701; // double C_3 = 0.5922083059862845; // double C_4 = -0.06517741879346212; // double C_5 = 0.09936322021509666; // double C_6 = -0.1883345440537925; // double C_7 = 0.02849094538051043; // double C_8 = 1.959940926229848; // double C_9 = -0.02236820885608115; // double CK_1 = -0.02452200335350193; // double CK_2 = 0.02660474893156432; // double CK_3 = -0.04770418667594679; // double CK_4 = 0.01284602485664462; // double CK_5 = 0.01825401538380585; // double CK_6 = 0.003645027941151212; // double CK_7 = 0.07750866083895219; // double CK_8 = 0.02706654817009965; // double CK_9 = 0.01115352331657436; // double CK_10 = 0.2703240873448836; // double C_0 = 1.125921604594353; // double C_1 = -0.4686941217654443; // double C_2 = -0.1114820681020402; // double C_3 = 0.6256505450040041; // double C_4 = -0.04724969749859539; // double C_5 = 0.1015107112011764; // double C_6 = -0.1868598401676765; // double C_7 = 0.002753811809744579; // double C_8 = 1.962925254827631; // double C_9 = 0.006913904797024802; // double CK_1 = -0.01961423092232162; // double CK_2 = 0.0263626002451694; // double CK_3 = -0.06984306808715408; // double CK_4 = 0.01265423995026528; // double CK_5 = 0.03651538489775871; // double CK_6 = 0.003964481821802226; // double CK_7 = 0.05925618682225351; // double CK_8 = 0.04024172550289172; // double CK_9 = 0.001420920976733996; // double CK_10 = 0.1968170158183649; double lambda_1 = sqrt(5.0/9.0); double lambda_2 = sqrt(5.0/6.0); ... ...
 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.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 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
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!