Commit d38021c1 authored by Timo Koch's avatar Timo Koch
Browse files

[tracer] Fix analytic storage derivative

Update 2ptracer reference solution after correcting derivative
parent 4bbb4f5a
......@@ -192,9 +192,10 @@ public:
const VolumeVariables& curVolVars,
const SubControlVolume& scv) const
{
const auto saturation = curVolVars.saturation(phaseIdx);
const auto porosity = curVolVars.porosity();
const auto rho = useMoles ? curVolVars.molarDensity() : curVolVars.density();
const auto d_storage = scv.volume()*porosity*rho/this->timeLoop().timeStepSize();
const auto d_storage = scv.volume()*porosity*rho*saturation/this->timeLoop().timeStepSize();
for (int compIdx = 0; compIdx < numComponents; ++compIdx)
partialDerivatives[compIdx][compIdx] += d_storage;
......
......@@ -105,33 +105,33 @@
3.4718e-11 8.36407e-11 1.43353e-10 1.82696e-10 1.82164e-10 1.47469e-10 9.95053e-11 5.70289e-11 2.81472e-11 1.0375e-11 3.65714e-12 1.19036e-11
4.35136e-11 9.73039e-11 1.58954e-10 2.08124e-10 2.29209e-10 2.18307e-10 1.82971e-10 1.36522e-10 9.14155e-11 5.52475e-11 3.02579e-11 1.50593e-11
0 0 0 0 0 0 0 0 0 0 0 1.01428e-12
5.00547e-12 2.17914e-11 5.88986e-11 1.04228e-10 1.32025e-10 1.28213e-10 1.00404e-10 6.56164e-11 3.66233e-11 1.77438e-11 7.54994e-12 8.20426e-12
5.00547e-12 2.17914e-11 5.88986e-11 1.04228e-10 1.32025e-10 1.28214e-10 1.0041e-10 6.56239e-11 3.66294e-11 1.77477e-11 7.55187e-12 8.20432e-12
3.06799e-11 8.1233e-11 1.42983e-10 1.83849e-10 1.83881e-10 1.49003e-10 1.00544e-10 5.75409e-11 2.81902e-11 1.16441e-11 4.25205e-12 9.92947e-12
3.85413e-11 9.06378e-11 1.51632e-10 1.99906e-10 2.19421e-10 2.07042e-10 1.71344e-10 1.26021e-10 8.31155e-11 4.94633e-11 2.6675e-11 1.30732e-11
0 0 0 0 0 0 0 0 0 0 0 8.2919e-13
7.69297e-12 3.28126e-11 7.94951e-11 1.27196e-10 1.47817e-10 1.33529e-10 9.88514e-11 6.18329e-11 3.33287e-11 1.57056e-11 6.54202e-12 7.62727e-12
3.01163e-11 8.09088e-11 1.42944e-10 1.83901e-10 1.83931e-10 1.4913e-10 1.00825e-10 5.79636e-11 2.87146e-11 1.23664e-11 4.73966e-12 8.10985e-12
7.69297e-12 3.28126e-11 7.94951e-11 1.27196e-10 1.47806e-10 1.33608e-10 9.90402e-11 6.20414e-11 3.3489e-11 1.5803e-11 6.59056e-12 7.62484e-12
3.01124e-11 8.09075e-11 1.42944e-10 1.83901e-10 1.83931e-10 1.4913e-10 1.00825e-10 5.79636e-11 2.87146e-11 1.23664e-11 4.73966e-12 8.10985e-12
3.21574e-11 7.99268e-11 1.38821e-10 1.87135e-10 2.07745e-10 1.96753e-10 1.62571e-10 1.18939e-10 7.7829e-11 4.58685e-11 2.44633e-11 1.18446e-11
0 0 0 0 0 0 0 0 0 0 0 2.11613e-12
1.66468e-11 5.66721e-11 1.14547e-10 1.59059e-10 1.639e-10 1.34541e-10 9.18954e-11 5.36905e-11 2.73357e-11 1.22897e-11 4.9305e-12 6.98314e-12
2.984e-11 8.08563e-11 1.43052e-10 1.83933e-10 1.83882e-10 1.49204e-10 1.01164e-10 5.85235e-11 2.93506e-11 1.29416e-11 5.10442e-12 7.0024e-12
1.66468e-11 5.66721e-11 1.14547e-10 1.59056e-10 1.64006e-10 1.35263e-10 9.30114e-11 5.47338e-11 2.80657e-11 1.27057e-11 5.13078e-12 6.93322e-12
2.97509e-11 8.08164e-11 1.43037e-10 1.83929e-10 1.83881e-10 1.49204e-10 1.01164e-10 5.85235e-11 2.93506e-11 1.29416e-11 5.10442e-12 7.0024e-12
2.75698e-11 7.15275e-11 1.28457e-10 1.77084e-10 1.99242e-10 1.89935e-10 1.57144e-10 1.14674e-10 7.46289e-11 4.36481e-11 2.30647e-11 1.10512e-11
0 0 0 0 0 0 0 0 0 0 0 7.38702e-12
3.94511e-11 1.00922e-10 1.64919e-10 1.94327e-10 1.75736e-10 1.29142e-10 7.9844e-11 4.26268e-11 2.00228e-11 8.38726e-12 3.16607e-12 6.74943e-12
3.09933e-11 8.23612e-11 1.44074e-10 1.84209e-10 1.83726e-10 1.49104e-10 1.01375e-10 5.89976e-11 2.98945e-11 1.33876e-11 5.37793e-12 6.43047e-12
3.94511e-11 1.00922e-10 1.64919e-10 1.94317e-10 1.76682e-10 1.3148e-10 8.25816e-11 4.48275e-11 2.14079e-11 9.11375e-12 3.49422e-12 6.51806e-12
3.04868e-11 8.20282e-11 1.43904e-10 1.84136e-10 1.837e-10 1.49095e-10 1.01373e-10 5.8997e-11 2.98943e-11 1.33876e-11 5.37793e-12 6.43047e-12
2.48511e-11 6.62647e-11 1.21834e-10 1.70756e-10 1.94107e-10 1.85991e-10 1.54041e-10 1.1217e-10 7.26678e-11 4.22304e-11 2.2143e-11 1.05167e-11
0 0 0 0 0 0 0 0 0 0 0 2.47819e-11
8.87271e-11 1.70043e-10 2.24253e-10 2.24151e-10 1.7859e-10 1.17817e-10 6.61452e-11 3.22979e-11 1.39532e-11 5.40507e-12 1.89764e-12 7.11272e-12
3.36789e-11 8.64466e-11 1.47348e-10 1.85666e-10 1.83858e-10 1.48864e-10 1.01341e-10 5.92413e-11 3.02482e-11 1.369e-11 5.56671e-12 6.13933e-12
2.33074e-11 6.31698e-11 1.17893e-10 1.67029e-10 1.91163e-10 1.83774e-10 1.52271e-10 1.10675e-10 7.14387e-11 4.1309e-11 2.15297e-11 1.01564e-11
8.87271e-11 1.70043e-10 2.24253e-10 2.24233e-10 1.80852e-10 1.21846e-10 7.01695e-11 3.52095e-11 1.56389e-11 6.22739e-12 2.24606e-12 6.53712e-12
3.22392e-11 8.51398e-11 1.46456e-10 1.85155e-10 1.83607e-10 1.48757e-10 1.013e-10 5.92277e-11 3.0244e-11 1.36889e-11 5.56642e-12 6.13926e-12
2.33073e-11 6.31698e-11 1.17893e-10 1.67029e-10 1.91163e-10 1.83774e-10 1.52271e-10 1.10675e-10 7.14387e-11 4.1309e-11 2.15297e-11 1.01564e-11
0 0 0 0 0 0 0 0 0 0 0 6.7536e-11
1.70756e-10 2.54439e-10 2.78503e-10 2.41919e-10 1.7329e-10 1.04908e-10 5.4716e-11 2.49826e-11 1.01198e-11 3.67821e-12 1.21148e-12 8.01687e-12
3.77554e-11 9.32148e-11 1.53821e-10 1.89623e-10 1.85296e-10 1.48989e-10 1.01195e-10 5.92348e-11 3.0373e-11 1.38354e-11 5.66911e-12 5.98646e-12
2.24576e-11 6.14316e-11 1.15662e-10 1.64931e-10 1.89533e-10 1.82559e-10 1.51281e-10 1.09807e-10 7.07001e-11 4.07427e-11 2.1148e-11 9.93072e-12
1.70756e-10 2.54439e-10 2.78503e-10 2.42217e-10 1.7661e-10 1.09834e-10 5.91625e-11 2.79857e-11 1.17659e-11 4.44284e-12 1.52073e-12 6.95282e-12
3.48707e-11 8.9868e-11 1.50953e-10 1.87581e-10 1.84056e-10 1.48336e-10 1.00893e-10 5.91103e-11 3.0327e-11 1.382e-11 5.66438e-12 5.98512e-12
2.24573e-11 6.14315e-11 1.15662e-10 1.64931e-10 1.89533e-10 1.82559e-10 1.51281e-10 1.09807e-10 7.07001e-11 4.07427e-11 2.1148e-11 9.93072e-12
0 0 0 0 0 0 0 0 0 0 0 1.3129e-10
2.60275e-10 3.25822e-10 3.14309e-10 2.49124e-10 1.66908e-10 9.58216e-11 4.79147e-11 2.11077e-11 8.27045e-12 2.90767e-12 9.24814e-13 9.45164e-12
4.3415e-11 1.02615e-10 1.63728e-10 1.96639e-10 1.88537e-10 1.49642e-10 1.00817e-10 5.87615e-11 3.00922e-11 1.37211e-11 5.63641e-12 5.88735e-12
2.20571e-11 6.06205e-11 1.14617e-10 1.63951e-10 1.88779e-10 1.82001e-10 1.50824e-10 1.09399e-10 7.03477e-11 4.04699e-11 2.09633e-11 9.82138e-12
2.60275e-10 3.25822e-10 3.14311e-10 2.49582e-10 1.70814e-10 1.01043e-10 5.23489e-11 2.39851e-11 9.8013e-12 3.6009e-12 1.19849e-12 7.76038e-12
3.86932e-11 9.6138e-11 1.57183e-10 1.91239e-10 1.84777e-10 1.47385e-10 9.9632e-11 5.82103e-11 2.98626e-11 1.36347e-11 5.60673e-12 5.87796e-12
2.20543e-11 6.06197e-11 1.14617e-10 1.63951e-10 1.88779e-10 1.82001e-10 1.50824e-10 1.09399e-10 7.03477e-11 4.04699e-11 2.09633e-11 9.82138e-12
</DataArray>
<DataArray type="Float32" Name="X^tracer_0" NumberOfComponents="1" format="ascii">
0 0 0 0 0 0 0 0 0 0 0 2.03283e-22
......@@ -235,33 +235,33 @@
5.78633e-10 1.39401e-09 2.38921e-09 3.04493e-09 3.03607e-09 2.45781e-09 1.65842e-09 9.50481e-10 4.69119e-10 1.72917e-10 6.09524e-11 1.98394e-10
7.25227e-10 1.62173e-09 2.64923e-09 3.46874e-09 3.82016e-09 3.63845e-09 3.04951e-09 2.27536e-09 1.52359e-09 9.20792e-10 5.04298e-10 2.50988e-10
0 0 0 0 0 0 0 0 0 0 0 1.69047e-11
8.34246e-11 3.6319e-10 9.81643e-10 1.73713e-09 2.20042e-09 2.13688e-09 1.6734e-09 1.09361e-09 6.10388e-10 2.9573e-10 1.25832e-10 1.36738e-10
8.34246e-11 3.6319e-10 9.81643e-10 1.73713e-09 2.20042e-09 2.1369e-09 1.6735e-09 1.09373e-09 6.10491e-10 2.95795e-10 1.25865e-10 1.36739e-10
5.11332e-10 1.35388e-09 2.38304e-09 3.06415e-09 3.06468e-09 2.48338e-09 1.67573e-09 9.59016e-10 4.69837e-10 1.94068e-10 7.08676e-11 1.65491e-10
6.42356e-10 1.51063e-09 2.5272e-09 3.33177e-09 3.65702e-09 3.4507e-09 2.85574e-09 2.10035e-09 1.38526e-09 8.24389e-10 4.44584e-10 2.17887e-10
0 0 0 0 0 0 0 0 0 0 0 1.38198e-11
1.28216e-10 5.46876e-10 1.32492e-09 2.11993e-09 2.46362e-09 2.22549e-09 1.64752e-09 1.03055e-09 5.55478e-10 2.6176e-10 1.09034e-10 1.27121e-10
5.01938e-10 1.34848e-09 2.38241e-09 3.06502e-09 3.06552e-09 2.4855e-09 1.68042e-09 9.6606e-10 4.78576e-10 2.06107e-10 7.89943e-11 1.35164e-10
1.28216e-10 5.46876e-10 1.32492e-09 2.11993e-09 2.46343e-09 2.2268e-09 1.65067e-09 1.03402e-09 5.5815e-10 2.63383e-10 1.09843e-10 1.27081e-10
5.01874e-10 1.34846e-09 2.3824e-09 3.06501e-09 3.06552e-09 2.4855e-09 1.68042e-09 9.6606e-10 4.78576e-10 2.06107e-10 7.89943e-11 1.35164e-10
5.35956e-10 1.33211e-09 2.31368e-09 3.11891e-09 3.46241e-09 3.27922e-09 2.70952e-09 1.98232e-09 1.29715e-09 7.64475e-10 4.07722e-10 1.9741e-10
0 0 0 0 0 0 0 0 0 0 0 3.52688e-11
2.77446e-10 9.44536e-10 1.90912e-09 2.65098e-09 2.73167e-09 2.24235e-09 1.53159e-09 8.94842e-10 4.55595e-10 2.04829e-10 8.2175e-11 1.16386e-10
4.97334e-10 1.3476e-09 2.38419e-09 3.06555e-09 3.0647e-09 2.48674e-09 1.68607e-09 9.75391e-10 4.89177e-10 2.15693e-10 8.50736e-11 1.16707e-10
2.77446e-10 9.44536e-10 1.90912e-09 2.65094e-09 2.73344e-09 2.25439e-09 1.55019e-09 9.1223e-10 4.67762e-10 2.11762e-10 8.5513e-11 1.15554e-10
4.95848e-10 1.34694e-09 2.38396e-09 3.06548e-09 3.06468e-09 2.48674e-09 1.68607e-09 9.75391e-10 4.89177e-10 2.15693e-10 8.50736e-11 1.16707e-10
4.59497e-10 1.19213e-09 2.14095e-09 2.9514e-09 3.3207e-09 3.16558e-09 2.61907e-09 1.91123e-09 1.24382e-09 7.27468e-10 3.84412e-10 1.84186e-10
0 0 0 0 0 0 0 0 0 0 0 1.23117e-10
6.57518e-10 1.68203e-09 2.74864e-09 3.23879e-09 2.92893e-09 2.15236e-09 1.33073e-09 7.10447e-10 3.33713e-10 1.39788e-10 5.27678e-11 1.12491e-10
5.16554e-10 1.37269e-09 2.40123e-09 3.07015e-09 3.06211e-09 2.48506e-09 1.68959e-09 9.83293e-10 4.98241e-10 2.23127e-10 8.96322e-11 1.07175e-10
6.57518e-10 1.68203e-09 2.74864e-09 3.23861e-09 2.9447e-09 2.19134e-09 1.37636e-09 7.47125e-10 3.56798e-10 1.51896e-10 5.82371e-11 1.08634e-10
5.08113e-10 1.36714e-09 2.39841e-09 3.06894e-09 3.06166e-09 2.48492e-09 1.68955e-09 9.83283e-10 4.98239e-10 2.23127e-10 8.96321e-11 1.07175e-10
4.14185e-10 1.10441e-09 2.03057e-09 2.84593e-09 3.23512e-09 3.09984e-09 2.56735e-09 1.8695e-09 1.21113e-09 7.0384e-10 3.69049e-10 1.75279e-10
0 0 0 0 0 0 0 0 0 0 0 4.13032e-10
1.47878e-09 2.83406e-09 3.73755e-09 3.73585e-09 2.97649e-09 1.96361e-09 1.10242e-09 5.38298e-10 2.32553e-10 9.00845e-11 3.16273e-11 1.18545e-10
5.61314e-10 1.44078e-09 2.45579e-09 3.09443e-09 3.0643e-09 2.48107e-09 1.68902e-09 9.87356e-10 5.04136e-10 2.28167e-10 9.27784e-11 1.02322e-10
1.47878e-09 2.83406e-09 3.73755e-09 3.73722e-09 3.01421e-09 2.03076e-09 1.16949e-09 5.86826e-10 2.60649e-10 1.0379e-10 3.74343e-11 1.08952e-10
5.37319e-10 1.419e-09 2.44093e-09 3.08592e-09 3.06012e-09 2.47928e-09 1.68834e-09 9.87128e-10 5.04067e-10 2.28148e-10 9.27736e-11 1.02321e-10
3.88456e-10 1.05283e-09 1.96488e-09 2.78381e-09 3.18605e-09 3.06291e-09 2.53785e-09 1.84458e-09 1.19064e-09 6.88484e-10 3.58829e-10 1.69273e-10
0 0 0 0 0 0 0 0 0 0 0 1.1256e-09
2.84594e-09 4.24065e-09 4.64171e-09 4.03199e-09 2.88816e-09 1.74847e-09 9.11933e-10 4.16377e-10 1.68663e-10 6.13035e-11 2.01913e-11 1.33614e-10
6.29256e-10 1.55358e-09 2.56369e-09 3.16038e-09 3.08827e-09 2.48316e-09 1.68659e-09 9.87246e-10 5.06217e-10 2.30589e-10 9.44851e-11 9.97743e-11
3.74294e-10 1.02386e-09 1.92769e-09 2.74885e-09 3.15889e-09 3.04265e-09 2.52135e-09 1.83011e-09 1.17833e-09 6.79046e-10 3.52467e-10 1.65512e-10
2.84594e-09 4.24065e-09 4.64171e-09 4.03696e-09 2.94351e-09 1.83057e-09 9.86042e-10 4.66429e-10 1.96098e-10 7.40473e-11 2.53455e-11 1.1588e-10
5.81178e-10 1.4978e-09 2.51588e-09 3.12634e-09 3.06761e-09 2.47227e-09 1.68155e-09 9.85172e-10 5.05451e-10 2.30333e-10 9.44064e-11 9.9752e-11
3.74288e-10 1.02386e-09 1.92769e-09 2.74885e-09 3.15889e-09 3.04265e-09 2.52135e-09 1.83011e-09 1.17833e-09 6.79046e-10 3.52467e-10 1.65512e-10
0 0 0 0 0 0 0 0 0 0 0 2.18817e-09
4.33792e-09 5.43037e-09 5.23849e-09 4.15207e-09 2.7818e-09 1.59703e-09 7.98578e-10 3.51795e-10 1.37841e-10 4.84612e-11 1.54136e-11 1.57527e-10
7.23584e-10 1.71025e-09 2.7288e-09 3.27732e-09 3.14229e-09 2.49404e-09 1.68029e-09 9.79358e-10 5.01536e-10 2.28685e-10 9.39402e-11 9.81225e-11
3.67618e-10 1.01034e-09 1.91028e-09 2.73252e-09 3.14631e-09 3.03334e-09 2.51373e-09 1.82331e-09 1.17246e-09 6.74499e-10 3.49388e-10 1.6369e-10
4.33792e-09 5.43037e-09 5.23852e-09 4.1597e-09 2.8469e-09 1.68405e-09 8.72482e-10 3.99752e-10 1.63355e-10 6.0015e-11 1.99749e-11 1.2934e-10
6.44887e-10 1.6023e-09 2.61971e-09 3.18732e-09 3.07961e-09 2.45642e-09 1.66053e-09 9.70172e-10 4.97711e-10 2.27245e-10 9.34454e-11 9.7966e-11
3.67572e-10 1.01033e-09 1.91028e-09 2.73252e-09 3.14631e-09 3.03334e-09 2.51373e-09 1.82331e-09 1.17246e-09 6.74499e-10 3.49388e-10 1.6369e-10
</DataArray>
<DataArray type="Float32" Name="rho" NumberOfComponents="1" format="ascii">
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
......
Supports Markdown
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