Pygmsh生成有限元网格数据应用到fipy中

最近我在研究基于Python有限体积求解微分方程。

首先通过Pygmsh生成一个正方体有限元网格数据文档,代码如下:

import pygmsh
with pygmsh.occ.Geometry() as geom:
    geom.characteristic_length_max = 0.1
    ellipsoid = geom.add_box([0.0, 0.0, 0.0], [1, 1, 1])
    mesh = geom.generate_mesh()
    mesh.write("D:\\desktop\\tiji\\zfx.vtk")

生成一个有限元网格数据文件zfx.vtk。将数据转成有限体积网格数据,利用有限体积求解微分方程,具体程序:

from fipy import CellVariable, Gmsh3D, TransientTerm, DiffusionTerm, Viewer, MayaviClient
import pygmsh

with pygmsh.occ.Geometry() as geom:
    geom.characteristic_length_max = 0.1
    ellipsoid = geom.add_box([0.0, 0.0, 0.0], [1, 1, 1])
    mesh = geom.generate_mesh()
    mesh.write("D:\\desktop\\tiji\\zfx1.vtk")
mesh = Gmsh3D('''        
              D:\\desktop\\tiji\\zfx1.vtk'''
              )
print(mesh.cellCenters)

但是

  • 4
    点赞
  • 10
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值