在OpenFOAM中获取网格详细信息

转载自: https://jibranhaider.weebly.com/blog/mesh-information-in-openfoam。

Updated: 29/04/2019
When developing your own code in OpenFOAM, it is almost certain that access to mesh information will be required in order to evaluate various parameters. In this post, I have provided a short description on how to access some of this information which could come in handy while programming. These have been characterized into mesh connectivities, mesh coordinates and some other useful mesh data.

Connectivities

// CELLS
const labelListList& cellPoints = mesh.cellPoints();        // Cell to point connectivity
const labelListList& cellEdges =  mesh.cellEdges();         // Cell to edge connectivity
const cellList& cells = mesh.cells();                       // Cell to face connectivity
const labelListList& cellCells =  mesh.cellCells();         // Cell to cell connectivity
// FACES
const faceList& faces = mesh.faces();                       // Face to node connectivity
const labelListList& faceEdges = mesh.faceEdges();          // Face to edge connectivity
const labelList& faceOwner = mesh.faceOwner();              // Face to owner cell
const labelList& faceNeighbour = mesh.faceNeighbour();      // Face to neighbour cell
// POINTS
const labelListList& pointPoints = mesh.pointPoints();      // Point to point connectivity
const labelListList& pointEdges = mesh.pointEdges();        // Point to edge connectivity
const labelListList& pointFaces = mesh.pointFaces();        // Point to face connectivity
const labelListList& pointCells = mesh.pointCells();        // Point to cell connectivity
// EDGES
const edgeList& edges = mesh.edges();                       // Edge to point connectivity
const labelListList& edgeFaces = mesh.edgeFaces();          // Edge to face connectivity
const labelListList& edgeCells = mesh.edgeCells();          // Edge to cell connectivity

Coordinates

const volVectorField& C = mesh.C();                         // Cell center coordinates
const surfaceVectorField& Cf = mesh.Cf();                   // Face center coordinates
const pointField& points = mesh.points();                   // Point coordinates// EXAMPLE: COMPUTE EDGE CENTRE COORDINATES
#include "vectorList.H"
vectorList Ce(edges.size(), vector::zero);                  // Edge centre coordinates
forAll(edges, edge)
{const label& own = edges[edge][0];                      // Label of owner pointconst label& nei = edges[edge][1];                      // Label of neighbour pointCe[edge] = 0.5*(points[own] + points[nei]);             // Compute edge centre coordinates
}

Boundaries

const polyBoundaryMesh& boundaryMesh = mesh.boundaryMesh();           // Boundary mesh// EXAMPLE: LOOP OVER ALL FACES
forAll(faces, face)
{if (mesh.isInternalFace(face))                                     // Check for an internal face{// Do your calculations}else                                                               // Boundary face{const label& patch = boundaryMesh.whichPatch(face);            // Boundary patchIDconst label& facei = boundaryMesh[patch].whichFace(face);      // Local boundary faceID// Do your calculations}
}// EXAMPLE: LOOP OVER FACES OF BOUNDARY PATCHES
forAll(mesh.boundary(), patch)
{const word& patchName = mesh.boundary()[patch].name();             // Boundary patch nameforAll(mesh.boundary()[patch], facei){const label& bCell = boundaryMesh[patch].faceCells()[facei];   // Boundary cellIDconst label& face = boundaryMesh[patch].start() + facei;       // FaceID// Do your calculations// U.boundaryField()[patch][facei] = vector::zero;}
}// EXAMPLE: LOOP OVER POINTS OF BOUNDARY PATCHES
forAll(mesh.boundary(), patch)
{forAll(boundaryMesh[patch].meshPoints(), pointi){const label& point = boundaryMesh[patch].meshPoints()[pointi];  // PointID// Do your calculations}
}

Other useful parameters

const scalarField& V = mesh.V();                            // Cell volumes
const surfaceVectorField& Sf = mesh.Sf();                   // Face area normal vectors
const surfaceScalarField& magSf = mesh.magSf();             // Face areas
const surfaceVectorField& N = Sf/magSf;                     // Face normal vectors
const label& nCells = mesh.nCells();                        // Total number of cells in the mesh
const label& nPoints = mesh.nPoints();                      // Total number of points in the mesh
const label& nInternalFaces = mesh.nInternalFaces();        // Number of internal faces in the mesh
const label& nInternalPoints = mesh.nInternalPoints();      // Number of internal points in

在OpenFOAM中获取网格详细信息相关推荐

PHP 获取服务器详细信息
PHP 获取系统信息,PHP 获取服务器详细信息获取系统类型及版本号: php_uname() (例:Windows NT ...
淘宝/天猫API接口，item_sku - 淘宝商品SKU详细信息查询，淘宝/天猫获取sku详细信息 API 返回值说明
欢迎使用淘宝API接口(item_sku - 淘宝商品SKU详细信息查询 ) 你好! 这是你使用我们的淘宝API接口获取商品详细接口说明. 测试地址:获取Key和secret 测试结果: Reques ...
淘宝/天猫如何获取sku API接口，item_sku - 获取sku详细信息
一.接口参数说明: 通过淘宝/天猫app中分享至PC端的短链接商品ID,通过调用item_sku - 获得sku接口,可以得到接口要的商品sku详情数据. 点击获取测试key和secret 请求地址: ...
Android获取歌曲详细信息
Android获取歌曲详细信息一.目标二.实现方案三.HTML方案四.JSON方案五.接下来六.Finally 在<Android实现歌词笔记构思>时,构思了如何在笔记中添加歌 ...
Python3 | 通过百度地图API获取商家详细信息（包括店名，地址，经纬度，电话）
一.打开百度地图开放平台,选择Web服务API . 官方文档对API接口的调用描述的比较详细,在调用之前先要申请百度地图的账号,然后申请密钥,获取密钥的步骤官方也有说明. 二.通过百度地图 ...
Lucene学习-深入Lucene分词器,TokenStream获取分词详细信息
Lucene学习-深入Lucene分词器,TokenStream获取分词详细信息在此回复牛妞的关于程序中分词器的问题,其实可以直接很简单的在词库中配置就好了,Lucene中分词的所有信息我们都可以从 ...
python查看系统进程_在Python中获取操作系统的进程信息
本文主要介绍在 Python 中使用 psutil 获取系统的进程信息. 1 概述 psutil 是 Python 的一个进程和系统工具集模块,通过使用 psutil,我们可以在 Python 中获取 ...
淘宝API接口：获取sku详细信息
今天分享的是淘宝平台API,item_sku - 获取sku详细信息点击注册测试地址(获取Key和secret) { "item": { "num_iid": ...
微信小程序控制台报错对应的服务器证书无效控制台输入 showRequestInfo() 可以获取更详细信息原因是ssl证书过期重新申请即可
微信小程序控制台报错对应的服务器证书无效.控制台输入 showRequestInfo() 可以获取更详细信息报错原因:域名下的SSL证书过期我这边报错的原因是ssl证书过期重新申请后就好了, ...
NX二次开发-UFUN获取工程图详细信息UF_DRAW_ask_drawing_info
NX二次开发-UFUN获取工程图详细信息UF_DRAW_ask_drawing_info NX9+VS2012#include <uf.h> #include <uf_draw.h& ...

在OpenFOAM中获取网格详细信息

Connectivities

Coordinates

Boundaries

Other useful parameters

在OpenFOAM中获取网格详细信息相关推荐

最新文章

热门文章