Skip to content
GitLab
Menu
Projects
Groups
Snippets
Loading...
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Sign in
Toggle navigation
Menu
Open sidebar
Lambert Theisen
fenicsR13
Commits
bd0beb79
Commit
bd0beb79
authored
Aug 11, 2019
by
Lambert Theisen
Browse files
Add full r13 testcase with inflow and change description of others
parent
8a55f146
Changes
3
Hide whitespace changes
Inline
Sidebyside
tests/r13/inputs/1_coeffs_nosources_norot_inflow_p1p1p1p1p1_stab.yml
View file @
bd0beb79
...
...
@@ 11,8 +11,8 @@ meshes:

../mesh/ring2.h5

../mesh/ring3.h5

../mesh/ring4.h5

../mesh/ring5.h5

../mesh/ring6.h5
# fails on Macbook 8GB
#
 ../mesh/ring5.h5
#
 ../mesh/ring6.h5 # fails on Macbook 8GB
#  ../mesh/ring7.h5 # works on iMac 48GB
#  ../mesh/ring8.h5
...
...
@@ 86,6 +86,6 @@ mass_source: 0
convergence_study
:
enable
:
True
exact_solution
:
esols/1_coeffs_nosources_norot_inflow.cpp
plot
:
Tru
e
plot
:
Fals
e
write_systemmatrix
:
False
rescale_pressure
:
False
tests/r13/referrors/1_coeffs_nosources_norot_inflow_p1p1p1p1p1_stab/errors.csv
0 → 100644
View file @
bd0beb79
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.11932437570021297,0.13595652832606953,0.1357094242477064,0.0822180095307,0.10512602608540127,0.0648883609082,0.40328949518145724,0.7582104995759621,1.849066997296789,1.10362282931,0.4938638585189256,0.760237819678,0.1301514425663452,0.101086048134,0.14351293987429084,0.106079819283,0.15130276131828505,0.0846972440835
0.6340332990709842,0.08837093734127678,0.11230519628506075,0.1202583618245216,0.072578267753,0.07144135832949701,0.0433922431076,0.2914934835196426,0.553613110345462,1.3793925027742815,1.39911865559,0.4621462427166392,0.785815559565,0.16765378623708282,0.0742363947388,0.18338301820895725,0.110647248116,0.2092287142596212,0.0896092176968
0.32904683851469807,0.038878913285623366,0.06513265378196054,0.052991357165989954,0.0383836531617,0.028442759582438318,0.0206055361137,0.09329458776349556,0.19874056809990387,0.38408651429722457,0.790456990514,0.22499621738432946,0.48030615872,0.09589310806745673,0.0486442593952,0.10748624465986259,0.0634059777899,0.11944866340176732,0.0587848984789
0.16754966839339377,0.013990056923050439,0.03174718715099378,0.012528960652537962,0.0118524978057,0.007708859068859351,0.00601859213766,0.02163501715131492,0.05437817015116003,0.14457204454593733,0.328392916107,0.07017533091090264,0.201500212761,0.026236999466893236,0.0159430939322,0.02863576480445802,0.0198678252547,0.031559409826889574,0.0176864517838
0.08734460120995041,0.004367241782151902,0.015156170732704366,0.0021479302498918066,0.00250793213301,0.0017273481915284917,0.00180739630822,0.003295234168161358,0.005766112224362463,0.04945190648447543,0.0542011784537,0.011585256605152479,0.0430322264031,0.004516352703061235,0.0039184868827,0.0048457828444671116,0.00326680175412,0.005183908326272927,0.0036084964106
tests/r13/test_r13_convergence.py
View file @
bd0beb79
...
...
@@ 50,7 +50,7 @@ class TestR13Convergence(object):
def
test_1_coeffs_sources_rot_noinflow_p1p1p1p1p1_stab
(
self
):
r
"""
Executes
westerkamp2019 coupled
system test and check with
Executes
full linear R13
system test and check with
reference errors.
========================= ==============================================
...
...
@@ 59,8 +59,16 @@ class TestR13Convergence(object):
:math:`\tau` :math:`1.0`
:math:`f_{\mathrm{mass}}` :math:`(1\frac{5R^2}{18\tau^2})\cos(\phi)`
:math:`f_{\mathrm{heat}}` :math:`0`
:math:`\theta_w^1` :math:`1.0`
:math:`v_t^1` :math:`10.0`
:math:`v_n^1` :math:`0`
:math:`p_w^1` :math:`0`
:math:`\gamma_w^1` :math:`10^{3}`
:math:`\theta_w^2` :math:`0.5`
:math:`v_t^2` :math:`0.0`
:math:`v_n^2` :math:`0`
:math:`p_w^2` :math:`0`
:math:`\gamma_w^2` :math:`10^{3}`
Elements :math:`P_1P_1P_1P_1P_1`
Stabilization CIP: :math:`\delta_1,\delta_2=1,\delta_3=0.01`
========================= ==============================================
...
...
@@ 73,7 +81,7 @@ class TestR13Convergence(object):
def
test_1_coeffs_sources_rot_noinflow_p2p2p2p2p2_stab
(
self
):
r
"""
Executes
westerkamp2019 coupled
system test and check with
Executes
full linear R13
system test and check with
reference errors.
========================= ==============================================
...
...
@@ 82,8 +90,16 @@ class TestR13Convergence(object):
:math:`\tau` :math:`1.0`
:math:`f_{\mathrm{mass}}` :math:`(1\frac{5R^2}{18\tau^2})\cos(\phi)`
:math:`f_{\mathrm{heat}}` :math:`0`
:math:`\theta_w^1` :math:`1.0`
:math:`v_t^1` :math:`10.0`
:math:`v_n^1` :math:`0`
:math:`p_w^1` :math:`0`
:math:`\gamma_w^1` :math:`10^{3}`
:math:`\theta_w^2` :math:`0.5`
:math:`v_t^2` :math:`0.0`
:math:`v_n^2` :math:`0`
:math:`p_w^2` :math:`0`
:math:`\gamma_w^2` :math:`10^{3}`
Elements :math:`P_2P_2P_2P_2P_2`
Stabilization CIP: :math:`\delta_1,\delta_2=1,\delta_3=0.01`
========================= ==============================================
...
...
@@ 93,3 +109,34 @@ class TestR13Convergence(object):
errors
=
name
+
"/"
+
"errors.csv"
referrors
=
"referrors/"
+
name
+
"/errors.csv"
self
.
compare_errors
(
errors
,
referrors
)
def
test_1_coeffs_nosources_norot_inflow_p1p1p1p1p1_stab
(
self
):
r
"""
Executes full linear R13 system test and check with reference errors.
Test case is similar to Torrilhon2017.
========================= ==============================================
Parameter Value
========================= ==============================================
:math:`\tau` :math:`1.0`
:math:`f_{\mathrm{mass}}` :math:`0`
:math:`f_{\mathrm{heat}}` :math:`0`
:math:`\theta_w^1` :math:`1.0`
:math:`v_t^1` :math:`0`
:math:`v_n^1` :math:`0`
:math:`p_w^1` :math:`0`
:math:`\gamma_w^1` :math:`10^{3}`
:math:`\theta_w^2` :math:`2.0`
:math:`v_t^2` :math:`1.00 \sin(\phi)`
:math:`v_n^2` :math:`+1.00 \cos(\phi)`
:math:`p_w^2` :math:`0.27 \cos(\phi)`
:math:`\gamma_w^2` :math:`10^{3}`
Elements :math:`P_1P_1P_1P_1P_1`
Stabilization CIP: :math:`\delta_1,\delta_2=1,\delta_3=0.01`
========================= ==============================================
"""
name
=
"1_coeffs_nosources_norot_inflow_p1p1p1p1p1_stab"
self
.
run_solver
(
"inputs/"
+
name
+
".yml"
)
errors
=
name
+
"/"
+
"errors.csv"
referrors
=
"referrors/"
+
name
+
"/errors.csv"
self
.
compare_errors
(
errors
,
referrors
)
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment