Fixed Rayleigh number

Author: brownbaerchen

pySDC/implementations/problem_classes/RayleighBenard.py

@@ -131,8 +131,10 @@ def __init__(
         self.Dz = S1 @ Dz
         self.Dzz = S2 @ Dzz
 
-        kappa = (Rayleigh * Prandtl) ** (-1 / 2.0)
-        nu = (Rayleigh / Prandtl) ** (-1 / 2.0)
+        # compute rescaled Rayleigh number to extract viscosity and thermal diffusivity
+        Ra = Rayleigh / (abs(BCs['T_top'] - BCs['T_bottom']) * self.axes[1].L ** 3)
+        kappa = (Ra * Prandtl) ** (-1 / 2.0)
+        nu = (Ra / Prandtl) ** (-1 / 2.0)
 
         # construct operators
         L_lhs = {

pySDC/projects/GPU/configs/RBC_configs.py

@@ -385,7 +385,7 @@ def get_controller_params(self, *args, **kwargs):
 class RayleighBenard_large(RayleighBenardRegular):
     # res_per_plume = 256
     # vertical_res = 1024
-    Ra = 2e7
+    Ra = 3.2e8
     relaxation_steps = 5
 
     def get_description(self, *args, **kwargs):

pySDC/projects/Resilience/RBC.py

@@ -14,7 +14,7 @@
 
 import numpy as np
 
-PROBLEM_PARAMS = {'Rayleigh': 2e4, 'nx': 256, 'nz': 128, 'max_cached_factorizations': 30}
+PROBLEM_PARAMS = {'Rayleigh': 3.2e5, 'nx': 256, 'nz': 128, 'max_cached_factorizations': 30}
 
 
 def u_exact(self, t, u_init=None, t_init=None, recompute=False, _t0=None):