From 68c879ad431061bbb0e53f1e6096c83f09abbac1 Mon Sep 17 00:00:00 2001 From: Michael Krayer Date: Sat, 31 Jul 2021 11:16:44 +0200 Subject: [PATCH] forgot to impose boundary conditions when setting noslip --- parallel.py | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/parallel.py b/parallel.py index a80ff39..5da4329 100644 --- a/parallel.py +++ b/parallel.py @@ -497,11 +497,11 @@ class PPP: import numpy as np self.symmetries[key] = np.zeros((3,3,3),dtype='i') for axis in range(3): - self.impose_mirror_symmetry(key,axis,-1) - self.impose_mirror_symmetry(key,axis,+1) + self.init_mirror_symmetry(key,axis,-1) + self.init_mirror_symmetry(key,axis,+1) return - def impose_mirror_symmetry(self,key,axis,wall): + def init_mirror_symmetry(self,key,axis,wall): '''Applies symmetry: [u,v,w](x,-y,z) -> [u,-v,w](x,y,z) Proposed by Lozano-Duran JFM 2016 and used for free-slip.''' if self.periodicity[axis]: @@ -540,6 +540,7 @@ class PPP: self.symmetries[key][sl] = 1 else: self.symmetries[key][sl] = -1 + self.impose_boundary_conditions(key) return def copy(self,key,key_out,skip_data=False,skip_symmetries=False):