Commit 368b959f authored by Lambert Theisen's avatar Lambert Theisen 🚀

Change names of stab params delta_i to i=(theta,u,p)

parent 0c70313b
......@@ -18,9 +18,9 @@ meshes:
# - stabilization: Must contain cip
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_1: Stabilization of grad(T)*grad(T_test) over edge
# - delta_2: Stabilization of grad(u)*grad(u_test) over edge
# - delta_3: Stabilization of grad(p)*grad(p_test) over edge
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
elements:
theta:
shape: Lagrange
......@@ -40,9 +40,9 @@ elements:
stabilization:
cip:
enable: True
delta_1: 1.0
delta_2: 1.0
delta_3: 0.1
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.1
# Formulation Parameters
# ======================
......
......@@ -18,9 +18,9 @@ meshes:
# - stabilization: Must contain cip
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_1: Stabilization of grad(T)*grad(T_test) over edge
# - delta_2: Stabilization of grad(u)*grad(u_test) over edge
# - delta_3: Stabilization of grad(p)*grad(p_test) over edge
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
elements:
theta:
shape: Lagrange
......@@ -40,9 +40,9 @@ elements:
stabilization:
cip:
enable: True
delta_1: 1.0
delta_2: 1.0
delta_3: 0.1
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.1
# Formulation Parameters
# ======================
......
......@@ -19,9 +19,9 @@ meshes:
# - stabilization: Must contain cip
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_1: Stabilization of grad(T)*grad(T_test) over edge
# - delta_2: Stabilization of grad(u)*grad(u_test) over edge
# - delta_3: Stabilization of grad(p)*grad(p_test) over edge
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
elements:
theta:
shape: Lagrange
......@@ -41,9 +41,9 @@ elements:
stabilization:
cip:
enable: True
delta_1: 1.0
delta_2: 1.0
delta_3: 0.1
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.1
# Formulation Parameters
# ======================
......
......@@ -18,9 +18,9 @@ meshes:
# - stabilization: Must contain cip
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_1: Stabilization of grad(T)*grad(T_test) over edge
# - delta_2: Stabilization of grad(u)*grad(u_test) over edge
# - delta_3: Stabilization of grad(p)*grad(p_test) over edge
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
elements:
theta:
shape: Lagrange
......@@ -40,9 +40,9 @@ elements:
stabilization:
cip:
enable: True
delta_1: 1.0
delta_2: 1.0
delta_3: 0.1
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.1
# Formulation Parameters
# ======================
......
......@@ -48,9 +48,9 @@ class Input:
# - stabilization: Must contain cip
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_1: Stabilization of grad(T)*grad(T_test) over edge
# - delta_2: Stabilization of grad(u)*grad(u_test) over edge
# - delta_3: Stabilization of grad(p)*grad(p_test) over edge
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
elements:
theta:
shape: Lagrange
......@@ -70,9 +70,9 @@ class Input:
stabilization:
cip:
enable: True
delta_1: 1.0
delta_2: 1.0
delta_3: 0.01
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.01
# Formulation Parameters
# ======================
......@@ -321,15 +321,15 @@ class Input:
"type": "boolean",
"required": True
},
"delta_1": {
"delta_theta": {
"type": "float",
"required": True
},
"delta_2": {
"delta_u": {
"type": "float",
"required": True
},
"delta_3": {
"delta_p": {
"type": "float",
"required": True
},
......
......@@ -69,9 +69,9 @@ class Solver:
self.kn = params["kn"]
self.chi_tilde = params["chi_tilde"]
self.use_cip = self.params["stabilization"]["cip"]["enable"]
self.delta_1 = self.params["stabilization"]["cip"]["delta_1"]
self.delta_2 = self.params["stabilization"]["cip"]["delta_2"]
self.delta_3 = self.params["stabilization"]["cip"]["delta_3"]
self.delta_theta = self.params["stabilization"]["cip"]["delta_theta"]
self.delta_u = self.params["stabilization"]["cip"]["delta_u"]
self.delta_p = self.params["stabilization"]["cip"]["delta_p"]
self.write_pdfs = self.params["postprocessing"]["write_pdfs"]
self.massflow = self.params["postprocessing"]["massflow"]
......@@ -483,9 +483,9 @@ class Solver:
bcs = self.bcs
kn = df.Constant(self.kn)
chi_tilde = df.Constant(self.chi_tilde)
delta_1 = df.Constant(self.delta_1)
delta_2 = df.Constant(self.delta_2)
delta_3 = df.Constant(self.delta_3)
delta_theta = df.Constant(self.delta_theta)
delta_u = df.Constant(self.delta_u)
delta_p = df.Constant(self.delta_p)
# Define custom measeasure for boundaries
df.ds = df.Measure("ds", domain=mesh, subdomain_data=boundaries)
......@@ -614,17 +614,17 @@ class Solver:
# 3) CIP Stabilization:
def j_theta():
return (
+ delta_1 * h_avg**3 *
+ delta_theta * h_avg**3 *
df.jump(df.grad(theta), n_vec) * df.jump(df.grad(kappa), n_vec)
) * df.dS
def j_u():
return (
+ delta_2 * h_avg**3 *
+ delta_u * h_avg**3 *
df.dot(df.jump(df.grad(u), n_vec), df.jump(df.grad(v), n_vec))
) * df.dS
def j_p():
return (
+ delta_3 * h_avg *
+ delta_p * h_avg *
df.jump(df.grad(p), n_vec) * df.jump(df.grad(q), n_vec)
) * df.dS
......
......@@ -25,9 +25,9 @@ meshes:
# - stabilization: Must contain cip
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_1: Stabilization of grad(T)*grad(T_test) over edge
# - delta_2: Stabilization of grad(u)*grad(u_test) over edge
# - delta_3: Stabilization of grad(p)*grad(p_test) over edge
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
elements:
theta:
shape: Lagrange
......@@ -47,9 +47,9 @@ elements:
stabilization:
cip:
enable: True
delta_1: 1.0
delta_2: 1.0
delta_3: 0.1
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.1
# Formulation Parameters
# ======================
......
......@@ -25,9 +25,9 @@ meshes:
# - stabilization: Must contain cip
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_1: Stabilization of grad(T)*grad(T_test) over edge
# - delta_2: Stabilization of grad(u)*grad(u_test) over edge
# - delta_3: Stabilization of grad(p)*grad(p_test) over edge
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
elements:
theta:
shape: Lagrange
......@@ -47,9 +47,9 @@ elements:
stabilization:
cip:
enable: False
delta_1: 1.0
delta_2: 1.0
delta_3: 0.1
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.1
# Formulation Parameters
# ======================
......
......@@ -25,9 +25,9 @@ meshes:
# - stabilization: Must contain cip
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_1: Stabilization of grad(T)*grad(T_test) over edge
# - delta_2: Stabilization of grad(u)*grad(u_test) over edge
# - delta_3: Stabilization of grad(p)*grad(p_test) over edge
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
elements:
theta:
shape: Lagrange
......@@ -47,9 +47,9 @@ elements:
stabilization:
cip:
enable: True
delta_1: 1.0
delta_2: 1.0
delta_3: 0.1
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.1
# Formulation Parameters
# ======================
......
......@@ -26,9 +26,9 @@ meshes:
# - stabilization: Must contain cip
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_1: Stabilization of grad(T)*grad(T_test) over edge
# - delta_2: Stabilization of grad(u)*grad(u_test) over edge
# - delta_3: Stabilization of grad(p)*grad(p_test) over edge
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
elements:
theta:
shape: Lagrange
......@@ -48,9 +48,9 @@ elements:
stabilization:
cip:
enable: True
delta_1: 1.0
delta_2: 1.0
delta_3: 0.1
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.1
# Formulation Parameters
# ======================
......
......@@ -65,7 +65,7 @@ class TestHeatConvergence(object):
============= =======================
:math:`Kn` :math:`0.1`
Elements :math:`P_1P_1`
Stabilization CIP, :math:`\delta_1=1`
Stabilization CIP, :math:`\delta_\theta=1`
============= =======================
"""
name = "heat_01_coeffs_p1p1_stab"
......@@ -83,7 +83,7 @@ class TestHeatConvergence(object):
============= =======================
:math:`Kn` :math:`10.0`
Elements :math:`P_2P_2`
Stabilization CIP, :math:`\delta_1=1`
Stabilization CIP, :math:`\delta_\theta=1`
============= =======================
"""
name = "heat_10_coeffs_p2p2_stab"
......@@ -101,7 +101,7 @@ class TestHeatConvergence(object):
============= =======================
:math:`Kn` :math:`0.1`
Elements :math:`P_2P_2`
Stabilization CIP, :math:`\delta_1=1`
Stabilization CIP, :math:`\delta_\theta=1`
============= =======================
"""
name = "heat_01_coeffs_p2p2_stab"
......@@ -119,7 +119,7 @@ class TestHeatConvergence(object):
============= =======================
:math:`Kn` :math:`0.1`
Elements :math:`P_2P_1`
Stabilization CIP, :math:`\delta_1=1`
Stabilization CIP, :math:`\delta_\theta=1`
============= =======================
"""
name = "heat_01_coeffs_p1p2_nostab"
......
......@@ -26,9 +26,9 @@ meshes:
# - stabilization: Must contain cip
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_1: Stabilization of grad(T)*grad(T_test) over edge
# - delta_2: Stabilization of grad(u)*grad(u_test) over edge
# - delta_3: Stabilization of grad(p)*grad(p_test) over edge
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
elements:
theta:
shape: Lagrange
......@@ -48,9 +48,9 @@ elements:
stabilization:
cip:
enable: True
delta_1: 1.0
delta_2: 1.0
delta_3: 0.01
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.01
# Formulation Parameters
# ======================
......
......@@ -25,9 +25,9 @@ meshes:
# - stabilization: Must contain cip
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_1: Stabilization of grad(T)*grad(T_test) over edge
# - delta_2: Stabilization of grad(u)*grad(u_test) over edge
# - delta_3: Stabilization of grad(p)*grad(p_test) over edge
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
elements:
theta:
shape: Lagrange
......@@ -47,9 +47,9 @@ elements:
stabilization:
cip:
enable: True
delta_1: 1.0
delta_2: 1.0
delta_3: 0.01
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.01
# Formulation Parameters
# ======================
......
......@@ -23,9 +23,9 @@ meshes:
# - stabilization: Must contain cip
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_1: Stabilization of grad(T)*grad(T_test) over edge
# - delta_2: Stabilization of grad(u)*grad(u_test) over edge
# - delta_3: Stabilization of grad(p)*grad(p_test) over edge
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
elements:
theta:
shape: Lagrange
......@@ -45,9 +45,9 @@ elements:
stabilization:
cip:
enable: True
delta_1: 1.0
delta_2: 1.0
delta_3: 0.01
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.01
# Formulation Parameters
# ======================
......
......@@ -68,25 +68,25 @@ class TestR13Convergence(object):
Rarefied Gas Flows”. In: Journal of Computational Physics 389
(2019).
========================= ==============================================
==================== ===================================================
Parameter Value
========================= ==============================================
:math:`Kn` :math:`1.0`
:math:`f_{\mathrm{mass}}` :math:`(1-\frac{5R^2}{18{Kn}^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:`\epsilon_w^1` :math:`0`
: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:`\epsilon_w^2` :math:`0`
Elements :math:`P_1P_1P_1P_1P_1`
Stabilization CIP: :math:`\delta_1,\delta_2=1,\delta_3=0.01`
========================= ==============================================
==================== ===================================================
:math:`Kn` :math:`1.0`
:math:`\dot{m}` :math:`(1-\frac{5R^2}{18{Kn}^2})\cos(\phi)`
:math:`r` :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:`\epsilon_w^1` :math:`0`
: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:`\epsilon_w^2` :math:`0`
Elements :math:`P_1P_1P_1P_1P_1`
Stabilization CIP: :math:`\delta_\theta,\delta_u=1,\delta_p=0.01`
==================== ===================================================
"""
name = "r13_1_coeffs_sources_rot_noinflow_p1p1p1p1p1_stab"
self.run_solver("inputs/" + name + ".yml")
......@@ -100,25 +100,25 @@ class TestR13Convergence(object):
Test case is similar to [WES2019]_.
========================= ==============================================
==================== ===================================================
Parameter Value
========================= ==============================================
:math:`Kn` :math:`1.0`
:math:`f_{\mathrm{mass}}` :math:`(1-\frac{5R^2}{18{Kn}^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:`\epsilon_w^1` :math:`0`
: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:`\epsilon_w^2` :math:`0`
Elements :math:`P_2P_2P_2P_2P_2`
Stabilization CIP: :math:`\delta_1,\delta_2=1,\delta_3=0.01`
========================= ==============================================
==================== ===================================================
:math:`Kn` :math:`1.0`
:math:`\dot{m}` :math:`(1-\frac{5R^2}{18{Kn}^2})\cos(\phi)`
:math:`r` :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:`\epsilon_w^1` :math:`0`
: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:`\epsilon_w^2` :math:`0`
Elements :math:`P_2P_2P_2P_2P_2`
Stabilization CIP: :math:`\delta_\theta,\delta_u=1,\delta_p=0.01`
==================== ===================================================
"""
name = "r13_1_coeffs_sources_rot_noinflow_p2p2p2p2p2_stab"
self.run_solver("inputs/" + name + ".yml")
......@@ -136,25 +136,25 @@ class TestR13Convergence(object):
simulations and model error estimation”. In: Journal of
Computational Physics 342 (2017), pp. 66–84.
========================= ==============================================
==================== ===================================================
Parameter Value
========================= ==============================================
:math:`Kn` :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:`\epsilon_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:`\epsilon_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`
========================= ==============================================
==================== ===================================================
:math:`Kn` :math:`1.0`
:math:`\dot{m}` :math:`0`
:math:`r` :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:`\epsilon_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:`\epsilon_w^2` :math:`10^{3}`
Elements :math:`P_1P_1P_1P_1P_1`
Stabilization CIP: :math:`\delta_\theta,\delta_u=1,\delta_p=0.01`
==================== ===================================================
"""
name = "r13_1_coeffs_nosources_norot_inflow_p1p1p1p1p1_stab"
self.run_solver("inputs/" + name + ".yml")
......
......@@ -25,9 +25,9 @@ meshes:
# - stabilization: Must contain cip
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_1: Stabilization of grad(T)*grad(T_test) over edge
# - delta_2: Stabilization of grad(u)*grad(u_test) over edge
# - delta_3: Stabilization of grad(p)*grad(p_test) over edge
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
elements:
theta:
shape: Lagrange
......@@ -47,9 +47,9 @@ elements:
stabilization:
cip:
enable: True
delta_1: 1.0
delta_2: 1.0
delta_3: 0.01
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.01
# Formulation Parameters
# ======================
......
......@@ -26,9 +26,9 @@ meshes:
# - stabilization: Must contain cip
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_1: Stabilization of grad(T)*grad(T_test) over edge
# - delta_2: Stabilization of grad(u)*grad(u_test) over edge
# - delta_3: Stabilization of grad(p)*grad(p_test) over edge
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
elements:
theta:
shape: Lagrange
......@@ -48,9 +48,9 @@ elements:
stabilization:
cip:
enable: True
delta_1: 1.0
delta_2: 1.0
delta_3: 0.01
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.01
# Formulation Parameters
# ======================
......
......@@ -25,9 +25,9 @@ meshes:
# - stabilization: Must contain cip
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_1: Stabilization of grad(T)*grad(T_test) over edge
# - delta_2: Stabilization of grad(u)*grad(u_test) over edge
# - delta_3: Stabilization of grad(p)*grad(p_test) over edge
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
elements:
theta:
shape: Lagrange
......@@ -47,9 +47,9 @@ elements:
stabilization:
cip:
enable: True
delta_1: 1.0
delta_2: 1.0
delta_3: 0.01
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.01
# Formulation Parameters
# ======================
......
......@@ -23,9 +23,9 @@ meshes:
# - stabilization: Must contain cip
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_1: Stabilization of grad(T)*grad(T_test) over edge
# - delta_2: Stabilization of grad(u)*grad(u_test) over edge
# - delta_3: Stabilization of grad(p)*grad(p_test) over edge
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
elements:
theta:
shape: Lagrange
......@@ -45,9 +45,9 @@ elements:
stabilization:
cip:
enable: False
delta_1: 1.0
delta_2: 1.0
delta_3: 0.01
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.01
# Formulation Parameters
# ======================
......
......@@ -24,9 +24,9 @@ meshes:
# - stabilization: Must contain cip
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_1: Stabilization of grad(T)*grad(T_test) over edge
# - delta_2: Stabilization of grad(u)*grad(u_test) over edge
# - delta_3: Stabilization of grad(p)*grad(p_test) over edge
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
elements:
theta:
shape: Lagrange
......@@ -46,9 +46,9 @@ elements:
stabilization:
cip:
enable: True
delta_1: 1.0
delta_2: 1.0
delta_3: 0.01
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.01
# Formulation Parameters
# ======================
......
......@@ -26,9 +26,9 @@ meshes:
# - stabilization: Must contain cip
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_1: Stabilization of grad(T)*grad(T_test) over edge
# - delta_2: Stabilization of grad(u)*grad(u_test) over edge
# - delta_3: Stabilization of grad(p)*grad(p_test) over edge
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
elements:
theta:
shape: Lagrange
......@@ -48,9 +48,9 @@ elements:
stabilization:
cip:
enable: True
delta_1: 1.0
delta_2: 1.0
delta_3: 0.01
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.01
# Formulation Parameters
# ======================
......
......@@ -67,7 +67,7 @@ class TestStressConvergence(object):
:math:`f_{\mathrm{mass}}` :math:`0`
:math:`v_t^1` :math:`10.0`
Elements :math:`P_1P_1P_1`
Stabilization CIP, :math:`\delta_2=1, \delta_3=0.01`
Stabilization CIP, :math:`\delta_u=1, \delta_p=0.01`
========================= =============================================