mwtkrayer
/
suspendtools
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

forgot to impose boundary conditions when setting noslip

This commit is contained in:
Michael Krayer 2021-07-31 11:16:44 +02:00
parent 1e2488a9bf
commit 68c879ad43
1 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
parallel.py
View File

@ -497,11 +497,11 @@ class PPP:
import numpy as np import numpy as np
self.symmetries[key] = np.zeros((3,3,3),dtype='i') self.symmetries[key] = np.zeros((3,3,3),dtype='i')
for axis in range(3): for axis in range(3):
self.impose_mirror_symmetry(key,axis,-1) self.init_mirror_symmetry(key,axis,-1)
self.impose_mirror_symmetry(key,axis,+1) self.init_mirror_symmetry(key,axis,+1)
return 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) '''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.''' Proposed by Lozano-Duran JFM 2016 and used for free-slip.'''
if self.periodicity[axis]: if self.periodicity[axis]:
@ -540,6 +540,7 @@ class PPP:
self.symmetries[key][sl] = 1 self.symmetries[key][sl] = 1
else: else:
self.symmetries[key][sl] = -1 self.symmetries[key][sl] = -1
self.impose_boundary_conditions(key)
return return
def copy(self,key,key_out,skip_data=False,skip_symmetries=False): def copy(self,key,key_out,skip_data=False,skip_symmetries=False):