def addFoamSymmTensorField(self, fieldpath, name, time, projectedField=False):
"""
"""
tSurf = TriSurfaceSymmTensor.readFromFoamFile(
fieldpath, self.triSurfaceMesh, time=time, projectedField=projectedField
)
self.fields[name] = tSurf
def addFoamSymmTensorField(self,
fieldpath,
name,
time,
projectedField=False):
'''
'''
tSurf = TriSurfaceSymmTensor.readFromFoamFile(
fieldpath,
self.triSurfaceMesh,
time=time,
projectedField=projectedField)
self.fields[name] = tSurf
def readReStressFromFoamFile(self,
varsFile,
pointsFile,
facesFile,
viewAnchor=(0,0,0),
xViewBasis=(1,0,0),
yViewBasis=(0,1,0),
dx=None,
dy=None,
interpolationMethod='cubic',
kind='min_E'):
'''
'''
tsm = TriSurfaceMesh.readFromFoamFile(pointsFile=pointsFile,
facesFile=facesFile,
viewAnchor=viewAnchor,
xViewBasis=xViewBasis,
yViewBasis=yViewBasis)
tsst = TriSurfaceSymmTensor.readFromFoamFile(varsFile=varsFile,
triSurfaceMesh=tsm,
time=0,
projectedField=False)
points = np.vstack((tsst.x,tsst.y)).T
#if not hasattr(self,'data'):
if not self.data.has_key('dx') or self.data.has_key('dy'):
print 'keys dx and dy does not exist'
if dx==None:
dxlist=[a for a in np.abs(np.diff(points[:,0])) if a>0]
dx=np.min(dxlist)
if dy==None:
dylist=[a for a in np.abs(np.diff(points[:,1])) if a>0]
dy=np.min(dylist)
MaxX=np.max(points[:,0])
MinX=np.min(points[:,0])
MaxY=np.max(points[:,1])
MinY=np.min(points[:,1])
extent=[MinX,MaxX,MinY,MaxY]
#print MinX,MaxX,MinY,MaxY
cellsX=int((MaxX-MinX)/dx)
cellsY=int((MaxY-MinY)/dy)
#print cellsX,cellsY
grid_y, grid_x = np.mgrid[MinY:MaxY:np.complex(0,cellsY),MinX:MaxX:np.complex(0,cellsX)]
triang = tsst.triangulation
uu_bar=self.interpolateField(tsst.txx, grid_x, grid_y, triang, method=interpolationMethod, kind=kind)
uv_bar=self.interpolateField(tsst.txy, grid_x, grid_y, triang, method=interpolationMethod, kind=kind)
uw_bar=self.interpolateField(tsst.tyy, grid_x, grid_y, triang, method=interpolationMethod, kind=kind)
vv_bar=self.interpolateField(tsst.tyy, grid_x, grid_y, triang, method=interpolationMethod, kind=kind)
vw_bar=self.interpolateField(tsst.tyz, grid_x, grid_y, triang, method=interpolationMethod, kind=kind)
ww_bar=self.interpolateField(tsst.tzz, grid_x, grid_y, triang, method=interpolationMethod, kind=kind)
vx_i=np.empty(uu_bar.shape)
vy_i=np.empty(uu_bar.shape)
vz_i=np.empty(uu_bar.shape)
vx_i[:]=np.NAN
vy_i[:]=np.NAN
vz_i[:]=np.NAN
self.vx=np.flipud(vx_i)
self.vy=np.flipud(vy_i)
self.vz=np.flipud(vz_i)
self.extent=extent
self.dx=dx
self.dy=dy
self.minX=MinX
self.maxX=MaxX
self.minY=MinY
self.maxY=MaxY
self.createDataDict()
print 'adding Tensor'
self.data['uu_bar']=np.flipud(uu_bar)
self.data['uv_bar']=np.flipud(uv_bar)
self.data['uw_bar']=np.flipud(uw_bar)
self.data['vv_bar']=np.flipud(vv_bar)
self.data['vw_bar']=np.flipud(vw_bar)
self.data['ww_bar']=np.flipud(ww_bar)
self.data['TKE_bar']=0.5*(self.data['uu_bar']+self.data['vv_bar']+self.data['ww_bar'])
else:
print 'dict exists'
MaxX=self.maxX
MinX=self.minX
MaxY=self.maxY
MinY=self.minY
cellsX=int((MaxX-MinX)/self.dx)
cellsY=int((MaxY-MinY)/self.dy)
#print cellsX,cellsY
grid_y, grid_x = np.mgrid[MinY:MaxY:np.complex(0,cellsY),MinX:MaxX:np.complex(0,cellsX)]
triang = tsst.triangulation
uu_bar=self.interpolateField(tsst.txx, grid_x, grid_y, triang, method=interpolationMethod, kind=kind)
uv_bar=self.interpolateField(tsst.txy, grid_x, grid_y, triang, method=interpolationMethod, kind=kind)
uw_bar=self.interpolateField(tsst.tyy, grid_x, grid_y, triang, method=interpolationMethod, kind=kind)
vv_bar=self.interpolateField(tsst.tyy, grid_x, grid_y, triang, method=interpolationMethod, kind=kind)
vw_bar=self.interpolateField(tsst.tyz, grid_x, grid_y, triang, method=interpolationMethod, kind=kind)
ww_bar=self.interpolateField(tsst.tzz, grid_x, grid_y, triang, method=interpolationMethod, kind=kind)
print 'adding Tensor'
self.data['uu_bar']=np.flipud(uu_bar)
self.data['uv_bar']=np.flipud(uv_bar)
self.data['uw_bar']=np.flipud(uw_bar)
self.data['vv_bar']=np.flipud(vv_bar)
self.data['vw_bar']=np.flipud(vw_bar)
self.data['ww_bar']=np.flipud(ww_bar)
self.data['TKE_bar']=0.5*(self.data['uu_bar']+self.data['vv_bar']+self.data['ww_bar'])
