Beispiel #1
0
def write_mesh_file_with_fields(field_list, name_=""):
    """
    write .vtu file containing all the cells and field values in the _OUTPUT directory
    
    :param fieldList:   a list of fields to include in the output. *father* of fields defines the mesh
    :type fieldList:    pyCFD_fields.fields.VolumeField
    :param name_:       filename to use
    :type name_:        string
    """
    t_start = time.time()
    print ''
    if name_ == "":
        print 'Writing field file...'
        file_name = __OUTITERDIR__ + 'fields.vtu'
    else:
        print 'Writing field file ' + name_ + ' ...'
        file_name = __OUTITERDIR__ + name_ + '.vtu'
    only_one_field = False
    if isinstance(field_list, fields.Field):
        only_one_field = True
    if only_one_field:
        mesh_object = field_list.father[0]
    else:
        mesh_object = field_list[0].father[0]
    field_list_ = []
    if only_one_field:
        field_list_.append(field_list)
    else:
        field_list_.extend(field_list)
    vtkTools.save_vtu_objects(mesh_object.cells, file_name, field_list_)
    print '   DONE in ' + str(time.time() - t_start) + ' s'
Beispiel #2
0
def write_patch_file_with_fields(mesh_, patch_name, volume_field_list):
    """
    write a patch mesh with field values into separate a .vtu file in the
    :mod:`pyCFD_config.config`.__OUTITERDIR__ directory
    
    :param mesh_:             a mesh object
    :type mesh_:              :class:`pyCFD_mesh.generic_emsh.GenericMesh`
    :param patch_name:        name of geometric patch
    :type patch_name:         string
    :param volume_field_list: list of volume fields
    :type volume_field_list:  :class:`pyCFD_fields.fields.VolumeField`
    """
    t_start = time.time()
    print ''
    print 'Writing patch file for ' + patch_name
    only_one_field = False
    if isinstance(volume_field_list, fields.Field):
        only_one_field = True
    field_list_ = []
    if only_one_field:
        field_list_.append(volume_field_list)
    else:
        field_list_.extend(volume_field_list)
    patch_ = mesh_.get_patch(patch_name)
    for field_ in field_list_:
        field_.update_boundary_values()
    fileName = __OUTITERDIR__ + str(patch_.name) + '.vtu'
    vtkTools.save_vtu_objects(patch_.faces, fileName, field_list_)
    print '   DONE for ' + str(
        patch_.name) + ' in ' + str(time.time() - t_start) + ' s'
Beispiel #3
0
def write_mesh_faces_with_field(surface_field_list, name_=""):
    """
    write .vtu file with all the faces
    
    :param mesh_: a mesh object (GenericMesh or derived)
    :type mesh_: pyCFD_mesh.generic_mesh.GenericMesh
    """
    t_start = time.time()
    print ''
    if name_ == "":
        file_name = __OUTITERDIR__ + 'meshAllFaces.vtu'
        print "Writing mesh file 'meshAllFaces.vtu' with all faces..."
    else:
        file_name = __OUTITERDIR__ + name_ + '.vtu'
        print "Writing mesh file '" + name_ + ".vtu' with all faces..."
    only_one_field = False
    if isinstance(surface_field_list, fields.Field):
        only_one_field = True
    if only_one_field:
        mesh_object = surface_field_list.father[0]
    else:
        mesh_object = surface_field_list[0].father[0]
    field_list = []
    if only_one_field:
        field_list.append(surface_field_list)
    else:
        field_list.extend(surface_field_list)
    vtkTools.save_vtu_objects(mesh_object.faces, file_name, field_list,
                              mesh_object)
    print '   DONE in ' + str(time.time() - t_start) + ' s'
Beispiel #4
0
def write_patch_files_with_fields(mesh_, volume_field_list):
    """
    write patch meshes with field values into separate .vtu files
    
    :param mesh_: a mesh object (GenericMesh or derived)
    :type mesh_: pyCFD_mesh.generic_mesh.GenericMesh
    :param volume_field_list: list of volume fields
    :type volume_field_list: pyCFD_fields.fields.VolumeField
        
    Output
        separate .vtu files containing patch faces for all the patches in the _OUTPUT directory
    """
    t_start = time.time()
    print ''
    print 'Writing patch files...'
    only_one_field = False
    if isinstance(volume_field_list, fields.Field):
        only_one_field = True
    field_list_ = []
    if only_one_field:
        field_list_.append(volume_field_list)
    else:
        field_list_.extend(volume_field_list)
    for patch_ in mesh_.patches:
        for field_ in field_list_:
            field_.update_boundary_values()
        fileName = __OUTITERDIR__ + str(patch_.name) + '.vtu'
        vtkTools.save_vtu_objects(patch_.faces, fileName, field_list_)
        print '   DONE for ' + str(patch_.name)
    print '   All DONE in ' + str(time.time() - t_start) + ' s'
Beispiel #5
0
def write_mesh_internal_faces_with_field(surface_field_list):
    """
    write .vtu file with all the internal faces
    
    :param mesh_: a mesh object (GenericMesh or derived)
    :type mesh_: pyCFD_mesh.generic_mesh.GenericMesh
    """
    t_start = time.time()
    print ''
    print 'Writing mesh file with internal faces...'
    file_name = __OUTITERDIR__ + 'meshIntFaces.vtu'
    internal_faces = []
    only_one_field = False
    if isinstance(surface_field_list, fields.Field):
        only_one_field = True
    if only_one_field:
        mesh_object = surface_field_list.father[0]
    else:
        mesh_object = surface_field_list[0].father[0]
    for face_ in mesh_object.faces:
        if face_.isBnd == False:
            internal_faces.append(face_)
    internal_face_mesh = sub_mesh.SubMesh(None, internal_faces)
    field_list = []
    if only_one_field:
        field_list.append(surface_field_list)
    else:
        field_list.extend(surface_field_list)
    vtkTools.save_vtu_objects(internal_face_mesh.faces, file_name, field_list,
                              mesh_object)
    print '   DONE in ' + str(time.time() - t_start) + ' s'
Beispiel #6
0
def write_mesh_faces_with_field(surface_field_list, name_=""):
    """
    write .vtu file with all the faces and a list of surface field in the 
    :mod:`pyCFD_config.config`.__OUTITERDIR__ directory
    
    :param mesh_:             a mesh object
    :type mesh_:              :class:`pyCFD_mesh.generic_emsh.GenericMesh`
    """
    t_start = time.time()
    print ''
    if name_ == "":
        file_name = __OUTITERDIR__+'meshAllFaces.vtu'
        print "Writing mesh file 'meshAllFaces.vtu' with all faces..."
    else:
        file_name = __OUTITERDIR__ + name_ + '.vtu'
        print "Writing mesh file '" + name_ + ".vtu' with all faces..."
    only_one_field = False
    if isinstance(surface_field_list, fields.Field):
        only_one_field = True
    if only_one_field:
        mesh_object = surface_field_list.father[0]
    else:
        mesh_object = surface_field_list[0].father[0]
    field_list = []
    if only_one_field:
        field_list.append(surface_field_list)
    else:
        field_list.extend(surface_field_list)
    vtkTools.save_vtu_objects(mesh_object.faces, file_name, field_list, mesh_object)
    print '   DONE in '+str(time.time()-t_start)+' s'
Beispiel #7
0
def write_mesh_file_with_fields(field_list, name_=""):
    """
    write .vtu file containing all the cells and field values in the
    :mod:`pyCFD_config.config`.__OUTITERDIR__ directory
    
    :param fieldList:   a list of fields to include in the output. *father* of fields defines the mesh
    :type fieldList:    :class:`pyCFD_fields.fields.VolumeField`
    :param name_:       filename to use
    :type name_:        string
    """
    t_start = time.time()
    print ''
    if name_ == "":
        print 'Writing field file...'
        file_name = __OUTITERDIR__+'fields.vtu'
    else:
        print 'Writing field file '+name_+' ...'
        file_name = __OUTITERDIR__+name_+'.vtu'
    only_one_field = False
    if isinstance(field_list, fields.Field):
        only_one_field = True
    if only_one_field:
        mesh_object = field_list.father[0]
    else:
        mesh_object = field_list[0].father[0]
    field_list_ = []
    if only_one_field:
        field_list_.append(field_list)
    else:
        field_list_.extend(field_list)
    vtkTools.save_vtu_objects(mesh_object.cells, file_name, field_list_)
    print '   DONE in '+str(time.time()-t_start)+' s'
Beispiel #8
0
def write_mesh_internal_faces(mesh_):
    """
    write .vtu file with all the internal faces in the
    :mod:`pyCFD_config.config`.__OUTITERDIR__ directory
    
    :param mesh_:             a mesh object
    :type mesh_:              :class:`pyCFD_mesh.generic_emsh.GenericMesh`
    """
    t_start = time.time()
    print ''
    print 'Writing mesh file with internal faces...'
    fileName = __OUTITERDIR__ + 'meshIntFaces.vtu'
    internal_faces = []
    for face_ in mesh_.faces:
        if face_.isBnd == False:
            internal_faces.append(face_)
    internal_face_mesh = sub_mesh.SubMesh(None, internal_faces)
    surf_normal = fields.SurfaceVectorField(internal_face_mesh, 'Sf')
    for i in range(len(surf_normal.A)):
        surf_normal.A[i] = surf_normal.father[0].faces[i].Sf
    fieldList = []
    fieldList.append(surf_normal)
    vtkTools.save_vtu_objects(internal_face_mesh.faces, fileName, fieldList,
                              mesh_)
    print '   DONE in ' + str(time.time() - t_start) + ' s'
Beispiel #9
0
def write_patch_file_with_fields(mesh_, patch_name, volume_field_list):
    """
    write a patch mesh with field values into separate a .vtu file in the
    :mod:`pyCFD_config.config`.__OUTITERDIR__ directory
    
    :param mesh_:             a mesh object
    :type mesh_:              :class:`pyCFD_mesh.generic_emsh.GenericMesh`
    :param patch_name:        name of geometric patch
    :type patch_name:         string
    :param volume_field_list: list of volume fields
    :type volume_field_list:  :class:`pyCFD_fields.fields.VolumeField`
    """
    t_start = time.time()
    print ''
    print 'Writing patch file for '+patch_name
    only_one_field = False
    if isinstance(volume_field_list, fields.Field):
        only_one_field = True
    field_list_ = []
    if only_one_field:
        field_list_.append(volume_field_list)
    else:
        field_list_.extend(volume_field_list)
    patch_ = mesh_.get_patch(patch_name)
    for field_ in field_list_:
        field_.update_boundary_values()
    fileName = __OUTITERDIR__+str(patch_.name)+'.vtu'
    vtkTools.save_vtu_objects(patch_.faces, fileName, field_list_)
    print '   DONE for '+str(patch_.name)+' in '+str(time.time()-t_start)+' s'
Beispiel #10
0
def write_mesh_internal_faces_with_field(surface_field_list):
    """
    write .vtu file with all the internal faces in the
    :mod:`pyCFD_config.config`.__OUTITERDIR__ directory
    
    :param mesh_:             a mesh object
    :type mesh_:              :class:`pyCFD_mesh.generic_emsh.GenericMesh`
    """
    t_start = time.time()
    print ''
    print 'Writing mesh file with internal faces...'
    file_name = __OUTITERDIR__+'meshIntFaces.vtu'
    internal_faces = []
    only_one_field = False
    if isinstance(surface_field_list, fields.Field):
        only_one_field = True
    if only_one_field:
        mesh_object = surface_field_list.father[0]
    else:
        mesh_object = surface_field_list[0].father[0]
    for face_ in mesh_object.faces:
        if face_.isBnd == False:
            internal_faces.append(face_)
    internal_face_mesh = sub_mesh.SubMesh(None,internal_faces)
    field_list = []
    if only_one_field:
        field_list.append(surface_field_list)
    else:
        field_list.extend(surface_field_list)
    vtkTools.save_vtu_objects(internal_face_mesh.faces, file_name, field_list, mesh_object)
    print '   DONE in '+str(time.time()-t_start)+' s'
Beispiel #11
0
def write_patch_files_with_fields(mesh_, volume_field_list):
    """
    write patch meshes with field values into separate .vtu files
    
    :param mesh_: a mesh object (GenericMesh or derived)
    :type mesh_: pyCFD_mesh.generic_mesh.GenericMesh
    :param volume_field_list: list of volume fields
    :type volume_field_list: pyCFD_fields.fields.VolumeField
        
    Output
        separate .vtu files containing patch faces for all the patches in the _OUTPUT directory
    """
    t_start = time.time()
    print ''
    print 'Writing patch files...'
    only_one_field = False
    if isinstance(volume_field_list, fields.Field):
        only_one_field = True
    field_list_ = []
    if only_one_field:
        field_list_.append(volume_field_list)
    else:
        field_list_.extend(volume_field_list)
    for patch_ in mesh_.patches:
        for field_ in field_list_:
            field_.update_boundary_values()
        fileName = __OUTITERDIR__+str(patch_.name)+'.vtu'
        vtkTools.save_vtu_objects(patch_.faces, fileName, field_list_)
        print '   DONE for '+str(patch_.name)
    print '   All DONE in '+str(time.time()-t_start)+' s'
Beispiel #12
0
def write_mesh_file(Mesh):
    """
    write .vtu file containing all the cells in the _OUTPUT directory
    
    :param mesh_:   mesh object
    :type mesh_:    pyCFD_mesh.generic_emsh.GenericMesh
    """
    t_start = time.time()
    print ''
    print 'Writing mesh file...'
    fileName = __OUTITERDIR__ + 'mesh.vtu'
    vtkTools.save_vtu_objects(Mesh.cells, fileName)
    print '   DONE in ' + str(time.time() - t_start) + ' s'
Beispiel #13
0
def write_mesh_file(Mesh):
    """
    write .vtu file containing all the cells in the _OUTPUT directory
    
    :param mesh_:   mesh object
    :type mesh_:    pyCFD_mesh.generic_emsh.GenericMesh
    """
    t_start = time.time()
    print ''
    print 'Writing mesh file...'
    fileName = __OUTITERDIR__+'mesh.vtu'
    vtkTools.save_vtu_objects(Mesh.cells,fileName)
    print '   DONE in '+str(time.time()-t_start)+' s'
Beispiel #14
0
def write_patch_files(mesh_):
    """
    write separate .vtu files containing patch faces for all the patches in the
    _OUTPUT directory
    
    :param mesh_:   mesh object
    :type mesh_:    pyCFD_mesh.generic_emsh.GenericMesh
    """
    t_start = time.time()
    print ''
    print 'Writing patch files...'
    for patch in mesh_.patches:
        fileName = __OUTITERDIR__ + str(patch.name) + '.vtu'
        vtkTools.save_vtu_objects(patch.faces, fileName)
        print '   DONE for ' + str(patch.name)
    print '   All DONE in ' + str(time.time() - t_start) + ' s'
Beispiel #15
0
def write_patch_files(mesh_):
    """
    write separate .vtu files containing patch faces for all the patches in the
    :mod:`pyCFD_config.config`.__OUTITERDIR__
    
    :param mesh_:   mesh object
    :type mesh_:    :class:`pyCFD_mesh.generic_emsh.GenericMesh`
    """
    t_start = time.time()
    print ''
    print 'Writing patch files...'
    for patch in mesh_.patches:
        fileName = __OUTITERDIR__+str(patch.name)+'.vtu'
        vtkTools.save_vtu_objects(patch.faces,fileName)
        print '   DONE for '+str(patch.name)
    print '   All DONE in '+str(time.time()-t_start)+' s'
Beispiel #16
0
def write_mesh_faces(mesh_):
    """
    write .vtu file with all the faces (internal + boundary)
    
    :param mesh_: a mesh object (GenericMesh or derived)
    :type mesh_: pyCFD_mesh.generic_mesh.GenericMesh
    """
    t_start = time.time()
    print ''
    print 'Writing mesh file with all faces...'
    fileName = __OUTITERDIR__ + 'meshAllFaces.vtu'
    surf_normal = fields.SurfaceVectorField(mesh_, 'Sf')
    for i in range(len(surf_normal.A)):
        surf_normal.A[i] = surf_normal.father[0].faces[i].Sf
    fieldList = []
    fieldList.append(surf_normal)
    vtkTools.save_vtu_objects(mesh_.faces, fileName, fieldList)
    print '   DONE in ' + str(time.time() - t_start) + ' s'
Beispiel #17
0
def write_mesh_faces(mesh_):
    """
    write .vtu file with all the faces (internal + boundary)
    
    :param mesh_: a mesh object (GenericMesh or derived)
    :type mesh_: pyCFD_mesh.generic_mesh.GenericMesh
    """
    t_start = time.time()
    print ''
    print 'Writing mesh file with all faces...'
    fileName = __OUTITERDIR__+'meshAllFaces.vtu'
    surf_normal = fields.SurfaceVectorField(mesh_,'Sf')
    for i in range(len(surf_normal.A)):
        surf_normal.A[i] = surf_normal.father[0].faces[i].Sf
    fieldList = []
    fieldList.append(surf_normal)
    vtkTools.save_vtu_objects(mesh_.faces,fileName,fieldList)
    print '   DONE in '+str(time.time()-t_start)+' s'
Beispiel #18
0
def write_mesh_faces(mesh_):
    """
    write .vtu file with all the faces in the 
    :mod:`pyCFD_config.config`.__OUTITERDIR__ directory
    
    :param mesh_:             a mesh object
    :type mesh_:              :class:`pyCFD_mesh.generic_emsh.GenericMesh`
    """
    t_start = time.time()
    print ''
    print 'Writing mesh file with all faces...'
    fileName = __OUTITERDIR__+'meshAllFaces.vtu'
    surf_normal = fields.SurfaceVectorField(mesh_,'Sf')
    for i in range(len(surf_normal.A)):
        surf_normal.A[i] = surf_normal.father[0].faces[i].Sf
    fieldList = []
    fieldList.append(surf_normal)
    vtkTools.save_vtu_objects(mesh_.faces,fileName,fieldList)
    print '   DONE in '+str(time.time()-t_start)+' s'
Beispiel #19
0
def write_mesh_faces(mesh_):
    """
    write .vtu file with all the faces in the 
    :mod:`pyCFD_config.config`.__OUTITERDIR__ directory
    
    :param mesh_:             a mesh object
    :type mesh_:              :class:`pyCFD_mesh.generic_emsh.GenericMesh`
    """
    t_start = time.time()
    print ''
    print 'Writing mesh file with all faces...'
    fileName = __OUTITERDIR__ + 'meshAllFaces.vtu'
    surf_normal = fields.SurfaceVectorField(mesh_, 'Sf')
    for i in range(len(surf_normal.A)):
        surf_normal.A[i] = surf_normal.father[0].faces[i].Sf
    fieldList = []
    fieldList.append(surf_normal)
    vtkTools.save_vtu_objects(mesh_.faces, fileName, fieldList)
    print '   DONE in ' + str(time.time() - t_start) + ' s'
Beispiel #20
0
def write_mesh_internal_faces(mesh_):
    """
    write .vtu file with all the internal faces
    
    :param mesh_: a mesh object (GenericMesh or derived)
    :type mesh_: pyCFD_mesh.generic_mesh.GenericMesh
    """
    t_start = time.time()
    print ''
    print 'Writing mesh file with internal faces...'
    fileName = __OUTITERDIR__+'meshIntFaces.vtu'
    internal_faces = []
    for face_ in mesh_.faces:
        if face_.isBnd == False:
            internal_faces.append(face_)
    internal_face_mesh = sub_mesh.SubMesh(None,internal_faces)
    surf_normal = fields.SurfaceVectorField(internal_face_mesh,'Sf')
    for i in range(len(surf_normal.A)):
        surf_normal.A[i] = surf_normal.father[0].faces[i].Sf
    fieldList = []
    fieldList.append(surf_normal)
    vtkTools.save_vtu_objects(internal_face_mesh.faces,fileName,fieldList,mesh_)
    print '   DONE in '+str(time.time()-t_start)+' s'