PyFR开发指南（一）：程序框架

最新推荐文章于 2024-05-23 09:57:58 发布
最新推荐文章于 2024-05-23 09:57:58 发布
阅读量1k 收藏 1
点赞数
分类专栏： 编程 文章标签： PyFR
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/iamzhtr/article/details/103114664
版权

本文基本上是官方手册的翻译，原文戳这里。

从哪儿下手

符号链接pyfr指向脚本_main_.py，一切开始于此。具体地说，函数process_run调用_process_common，进而调用get_solver，返回一个Integrator，即Controller和Stepper的复合体。Integrator包含一个方法run，调用它来运行仿真。

Controller

Controller负责沿时间推进仿真。具体地说，Controller包含一个方法advance_to，把仿真推进到一个特定的时间。PyFR 1.9.0提供三种物理时间Controller：

  • pyfr.integrators.std.controllers.StdNoneController
  • pyfr.integrators.std.controllers.StdPIController
  • pyfr.integrators.dual.phys.controllers.DualNoneController

各种物理时间Controller的继承关系如图：
在这里插入图片描述除了物理时间，还有两种伪时间Controller：

  • pyfr.integrators.dual.pseudo.pseudocontrollers.DualNonePseudoController
  • pyfr.integrators.dual.pseudo.pseudocontrollers.DualPIPseudoController

继承关系为：
在这里插入图片描述

Stepper

Stepper负责把仿真推进一个时间步。具体地说，Stepper包含一个方法step，把System推进一个时间步。有八种Stepper：

  • pyfr.integrators.std.steppers.StdEulerStepper

  • pyfr.integrators.std.steppers.StdRK4Stepper

  • pyfr.integrators.std.steppers.StdRK34Stepper

  • pyfr.integrators.std.steppers.StdRK45Stepper

  • pyfr.integrators.std.steppers.StdTVDRK3Stepper

  • pyfr.integrators.dual.phys.steppers.DualBDF2Stepper

  • pyfr.integrators.dual.phys.steppers.DualBDF3Stepper

  • pyfr.integrators.dual.phys.steppers.DualBackwardEulerStepper

继承关系为：
在这里插入图片描述
同样也有伪Stepper，把仿真推进一个伪时间步。通过一个双时间步形式，可以用它们使隐式Stepper时间步收敛。有六种伪Stepper：

  • pyfr.integrators.dual.pseudo.pseudosteppers.DualDenseRKPseudoStepper[source]

  • pyfr.integrators.dual.pseudo.pseudosteppers.DualRK4PseudoStepper[source]

  • pyfr.integrators.dual.pseudo.pseudosteppers.DualTVDRK3PseudoStepper[source]

  • pyfr.integrators.dual.pseudo.pseudosteppers.DualEulerPseudoStepper[source]

  • pyfr.integrators.dual.pseudo.pseudosteppers.DualRK34PseudoStepper[source]

  • pyfr.integrators.dual.pseudo.pseudosteppers.DualRK45PseudoStepper

注意其中的DualDenseRKPseudoStepper包含一类伪Stepper，其系数读取自遵循如下命名规范的.txt文件：

{scheme name}-s{stage count}-p{temporal order}-sp{optimal spatial polynomial order}.txt

继承关系为：
在这里插入图片描述

System

System包含系统的各种信息，包括Elements, Interfaces和后端Backend。System包含一个方法rhs（right-hand-side缩写），用于获取每个solution point上的通量散度（即方程的右端项）。rhs启动预先生成并加载到队列里的各种kernel。System还包含一个方法_gen_kernels，用于生成System需要的所有kernel。kernel是一个"一次性"类的实例，包含一个方法run用于实现需要的功能。每个kernel都由一个kernel提供者生成。PyFR有许多类型的kernel提供者。一个Pointwise Kernel Provider（逐点kernel提供者）生成逐点kernel，比如黎曼求解器和通量函数等。这些逐点kernel是用一种修改版本的Mako，即PyFR-Mako实现的。下文有介绍，更多请参阅Mako官网。有四种System供选择：

  • pyfr.solvers.aceuler.system.ACEulerSystem
  • pyfr.solvers.acnavstokes.system.ACNavierStokesSystem
  • pyfr.solvers.euler.system.EulerSystem
  • pyfr.solvers.navstokes.system.NavierStokesSystem

继承关系为：
在这里插入图片描述

Elements

Elements（注意s）包含一组单元信息。有四种Elements：

  • pyfr.solvers.aceuler.elements.ACEulerElements
  • pyfr.solvers.acnavstokes.elements.ACNavierStokesElements
  • pyfr.solvers.euler.elements.EulerElements
  • pyfr.solvers.navstokes.elements.NavierStokesElements

继承关系为：
在这里插入图片描述

Interfaces

Interfaces（注意s）包含一组界面信息。有八种（非边界）界面：

  • pyfr.solvers.aceuler.inters.ACEulerIntInters

  • pyfr.solvers.aceuler.inters.ACEulerMPIInters

  • pyfr.solvers.acnavstokes.inters.ACNavierStokesIntInters

  • pyfr.solvers.acnavstokes.inters.ACNavierStokesMPIInters

  • pyfr.solvers.euler.inters.EulerIntInters

  • pyfr.solvers.euler.inters.EulerMPIInters

  • pyfr.solvers.navstokes.inters.NavierStokesIntInters

  • pyfr.solvers.navstokes.inters.NavierStokesMPIInters

继承关系为：
在这里插入图片描述

Backend

Backend包含了后端信息。有三种（原文是四种，疑似笔误）后端可选：

  • pyfr.backends.cuda.base.CUDABackend
  • pyfr.backends.opencl.base.OpenCLBackend
  • pyfr.backends.openmp.base.OpenMPBackend

继承关系为：
在这里插入图片描述

逐点kernel提供者

逐点kernel提供者（简写为PKP）负责逐点生成kernel。具体地说，PKP包含一个方法register，负责把一个新方法添加到一个PKP实例里。这个新方法被调用时会返回一个kernel。所谓的kernel，是一个“一次性”类的实例，其包含的方法run用于实现需要的功能。kernel的功能由PyFR-Mako实现，因此PKP还包含一个方法_render_kernel，用于把PyFR-Mako渲染为底层代码。_render_kernel首先为Mako设定上下文（例如Backend细节等），然后使用Mako渲染PyFR-Mako内容。当Mako遇到一个pyfr:kernel时，会创建一个kernel生成器的实例，用于渲染pyfr:kernel的主体内容。有三种（原文写作四种，疑似笔误）PKP可选：

  • pyfr.backends.cuda.provider.CUDAPointwiseKernelProvider
  • pyfr.backends.opencl.provider.OpenCLPointwiseKernelProvider
  • pyfr.backends.openmp.provider.OpenMPPointwiseKernelProvider

继承关系为：
在这里插入图片描述

Kernel生成器

Kernel生成器把pyfr:kernel中的PyFR-Mako渲染为底层代码。具体地说，kernel生成器包含一个方法render。有三种（原文写作四种，疑似笔误）kernel生成器：

  • pyfr.backends.cuda.generator.CUDAKernelGenerator
  • yfr.backends.opencl.generator.OpenCLKernelGenerator
  • pyfr.backends.openmp.generator.OpenMPKernelGenerator

继承关系为：
在这里插入图片描述

后记

该网页最后有关于PyFR-Mako的介绍，放在下一篇翻译。

鼠标君
关注
专栏目录
服务器应用程序无法使用,使用通道2运行服务器时无法导入ASGI_应用程序模块
weixin_36429576的博客
08-02 1442
我已经遵循了频道教程，但是在运行这些错误抛出时包的版本是频道==2.1.2Django==2.0.4我错过了什么？在settings.py中INSTALLED_APPS = ["channels"....]ROOT_URLCONF = 'myapp.urls'ASGI_APPLICATION = "myapp.routing.application"添加了文件mayapp/routing.pyfr...
PyFR上手指南
iamzhtr的博客
11-17 2511
通量重构(FR)格式是一种高精度紧致型流体数值格式。PyFR 是一个开源 CFD 包，采用非常精确的FR格式来处理世界上最具挑战性的一些流体流动问题。尤其是涉及不平稳湍流的问题。相关的博客可见Nvidia博客。PyFR的官网上提供了下载和手册。本文是基于以上内容的一篇入门指南，操作平台为Ubuntu 19.10。 安装 安装基本的依赖，命令为 sudo pip3 install xxx 其中...
PyFR:吡咯
03-20
吡咯 概述 PyFR是一个基于Python的开源框架，用于使用Huynh的Flux Reconstruction方法解决流传输架构上的对流扩散类型问题。该框架旨在解决包含各种元素类型的混合非结构化网格上的一系列控制系统。它还设计为通过使用从Mako模板引擎派生的内置领域特定语言来针对各种硬件平台。 例子 examples目录中提供了使用PyFR的examples 。当前可用的示例有： 2D Couette流 二维Euler涡 2D不可压缩气缸流量 贡献 要为PyFR做出贡献，请遵循以下步骤： 将此存储库分叉到您的GitHub帐户 在分支的存储库中创建一个新分支 在您的新分支中进行更改 通过创建“拉取请求”将更改提交到原始PyFR存储库的develop分支 最终，对develop分支的修改将合并到master分支中以发布新版本。 帮助与社区 有关PyFR的用法和结构的信息可以在和。 在开发
推荐开源项目：PyFR——流体动力学的Python框架
最新发布
gitblog_00093的博客
05-23 464
推荐开源项目：PyFR——流体动力学的Python框架 去发现同类优质开源项目:https://gitcode.com/ 项目介绍 PyFR 是一款基于 Python 的开源框架，用于解决在流架构上的扩散-对流类型问题，其采用了 Huynh 提出的 Flux Reconstruction 方法。设计上，PyFR 能处理混合无结构网格中的各种元素类型的广义守恒方程系统，并能针对多种硬件平台进行优化，...
开发功能需求（FR）组件的技术指南
N201871643的博客
05-11 502
在开发功能需求组件时，遵循上述注意事项和技术指南是至关重要的。这不仅能够帮助团队高效地交付符合预期的产品，还能够确保产品的长期成功和可持续发展。通过对需求的深入理解、精心设计、周密实施和严格测试，可以构建出既健壮又灵活的软件系统。
Julia调用PyFR
iamzhtr的博客
11-19 509
PyFR由Python3写成，并且不提供接口，所以我们需要自己动手，改造一番，才能在Julia里调用它的内部函数。上回我们介绍了PyFR的运行示例：二维欧拉涡。在该示例中，有一个转换.msh文件的步骤。如何在Julia调用PyFR实现这个步骤呢？一切开始前，假定你已经安装了Julia的PyCall包，并且读过PyCall的Github项目的readme文档。 首先，我们要安装PyFR。这次不是在P...
StageFright框架流程解读
menuconfig的专栏
06-23 8104
1、    StageFright介绍    Androidfroyo版本多媒体引擎做了变动，新添加了stagefright框架，并且默认情况android选择stagefright，并没有完全抛弃opencore，主要是做了一个OMX层，仅仅是对opencore的omx-component部分做了引用。stagefright是在MediaPlayerService这一层加入的，和
angular-bootstrap-include-j7pyfr：用StackBlitz:high_voltage:创建
02-10
标题 "angular-bootstrap-include-j7pyfr：用StackBlitz:high_voltage:创建" 暗示了这是一个关于使用Angular框架结合Bootstrap库进行前端开发的项目，并且该项目可能是在StackBlitz在线代码编辑器上创建的。...
PyPI 官网下载 | pyfr-1.10.0.tar.gz
01-29
- **Python库**：PyFR作为一个Python库，可以无缝融入Python开发环境，利用Python的丰富生态进行扩展和集成。 总体来说，PyFR 1.10.0是一个强大的CFD求解器，它集成了现代计算技术，如分布式计算和云原生，使得复杂...
双曲柱面matlab代码-CFD-Simulation:用MATLAB和基于Python的PyFR编写的流体求解器代码的数值模拟。方法包括有限
05-25
一旦要研究阻力效应，这就是一个重要的假设。 当前涡流强度设置为R = 1.5，马赫数为0.4，特征半径gamma为1.4。 有关更多信息，请咨询。 保守变量的存储方式如下： 字段初始化完成后，我们就可以开始采用集中方案及时...
fr8016 sdk源码结构和开发方式
weixin_43810563的博客
03-30 948
下面的代码中，user_task_func是一个任务，需要被其他任务通过事件驱动，当其他任务调用os_msg_post向user_task_func发送os_event_t时，user_task_func得到执行。sdk开发中，任务采用事件驱动，不同任务间通过消息os_event_t *进行通信，当任务A产生的事件可以驱动任务B时，任务A通过os_msg_post向任务B发送消息，任务B接收消息，并从消息中获取事件id和数据，根据事件id做相应的处理。下面按照协议栈处理流程给出接口的说明。
python global用法_用Python设计PCR引物: Primer3py
weixin_39717026的博客
12-07 1545
PCR引物设计应该算是生物实验基本技能吧，工具也非常多。Primer-BLAST、Primer Premier 都是比较经典的软件，除此之外还有很多新的在线设计软件，也很方便。不过今天不是要讲怎么用这些现成软件来设计引物，而是要讲讲怎么用Python来做引物设计。Python拥有众多的第三方包，其中有个叫: Primer3-pyhttps://pypi.org/project/prim...
FR技术简介
将勤补拙
12-14 4773
FR：帧中继 DLCL：标识符 帧中继( Frame Relay,FR)网络用虚电路来连接网络两端的帧中继设备。每条虚电路用数据链路连接标 识符( Data Link Connection Identifier,DLCL)定义了一条帧中继连接通道。帧中继协议是一种统计复用协 议,它能够在单一物理传输线路上提供多条虚电路。帧中继用户的接入速率在64kbis~2 Mbit/s,甚至可达...
Julia安装与运行
iamzhtr的博客
08-20 4610
（1）Julia的安装很简单。在官网https://julialang.org/下载exe直接安装即可。 （2）运行方法： 　　（2.1）点击Julia图标，打开所谓的REPL界面（其实就是Julia自己的命令行），像Matlab一样使用。 　　 　　（2.2）用任意文本编辑器创建一个源代码文件test.jl（注意后缀是.jl）。在REPL中用include命令运行，注意路径层级用\代替\...
Atom配置Julia/latex
iamzhtr的博客
08-23 1630
在Atom运行Julia/latex
在VS Code中配置码云进行代码管理（解决你的报错）
iamzhtr的博客
06-23 1110
在网上下载安装Git最新版本，安装时有个选项可选择Visual Studio Code作为默认编辑器，其余全部默认。 安装后，从开始菜单打开Git Bash，依次输入以下命令： ssh-keygen -t rsa -C "gitee.com" cat ~/.ssh/id_rsa.pub 注意第一个命令会提示你输入名称，千万不要填，直接连按三下回车跳过。第二个命令会返回一串SSH公钥。如果你填了名...
getfem++ 热弹性-电耦合算例 Python 实现
iamzhtr的博客
10-22 851
原文戳这里。 问题设定 考虑一个厚度为 ε\varepsilonε 的二维平板，受到外力、电势和加热作用。令 Ω⊂IR2\Omega \subset \rm I\hspace{-0.15em}R^2Ω⊂IR2为平板计算域，记 θ:Ω→IR\theta : \Omega \rightarrow \rm I\hspace{-0.15em}Rθ:Ω→IR 为温度场（摄氏度），V:Ω→IRV : \Omega \rightarrow \rm I\hspace{-0.15em}RV:Ω→IR 为电势场，u:Ω→IRu
方腔驱动流问题python
05-09
3. PyFR：是一个用于求解Navier-Stokes方程的高阶有限元方法库，可高效地求解方腔驱动流问题。 4. Fluidity：是一个用于求解多相流和流变学问题的计算流体力学软件，可使用有限元和有限体积方法求解Navier-Stokes...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值