firedrake求解NS方程

已于 2023-04-12 19:37:47 修改
阅读量1.1k 收藏 1
点赞数
分类专栏： 传统数值方法 NS方程求解 文章标签： 概率论 几何学 线性代数
于 2022-05-06 19:50:15 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/forrestguang/article/details/124604503
版权
传统数值方法 同时被 2 个专栏收录
35 篇文章 20 订阅 ¥199.90 ¥299.90
订阅专栏
NS方程求解
8 篇文章 28 订阅
订阅专栏
本文介绍了如何用firedrake解决流体力学中的Bump和Recirculating Flow问题，详细阐述了弱形式的推导过程，并探讨了出流边界的复杂边界条件处理。
摘要由CSDN通过智能技术生成

博客内容参考本人知乎添加链接描述

Bump问题

在这里插入图片描述

使用firedrake求解正问题，
在这里插入图片描述
在这里插入图片描述

下面列出三个bump.geo分别对应三个边界
bump1,geo

lr1 = 5e-1;
lr2 =
了解本专栏 订阅专栏 解锁全文
firedrake的安装（docker环境）
forrestguang的博客
04-19 1732
注意以下安装过程都是默认在docker环境下的服务器进行的，用户需要被服务器管理员添加锦docker组 1:百度搜索firedrake官网，进去，然后Ctrl + f，搜索docker，找到下面这个部分 2：点击进入上面提到的docker image，复制右边的命令：docker pull firedrakeproject/firedrake，把这段命令复制到服务器终端，然后服务器会开始下载跟firedrake相关的库，大概有两三个G. 3：下载完了以后，键入命令： docker run -it – p
firedrake创建mesh和求解具体问题
forrestguang的博客
05-06 435
界面问题 −∇⋅(σ∇u)+u=5,Ω-\nabla \cdot(\sigma \nabla u) + u = 5,\Omega−∇⋅(σ∇u)+u=5,Ω u∣∂Ω1=0u|_{\partial \Omega_1} = 0u∣∂Ω1​​=0 [σ∇u⋅n]=3,in∂Ω2[\sigma\nabla u \cdot n] = 3,in \partial \Omega_2[σ∇u⋅n]=3,in∂Ω2​ ∂Ω1\partial \Omega_1∂Ω1​为矩形边界 ∂Ω2\partial \Omega_2∂Ω2
求解2维NS方程的Matlab源程序
07-01
2D NS方程求解Matlab源程序<br>有用户界面<br>非结构化网格<br>Matlab与C混合编程
CFD求解NS方程代码
12-26
求解NS方程会碰到一个问题就是CFL问题，该代码研究时间参数设置对结果影响
有限元开源代码firedrake_Python编写的有限元
05-18
Firedrake is an automated system for the portable solution of partial differential equations using the finite element method (FEM). Firedrake enables users to employ a wide range of discretisations to an infinite variety of PDEs and employ either conventional CPUs or GPUs to obtain the solution.
fireshape求解NS方程约束问题
forrestguang的博客
05-19 96
代码主要包括三个部分 第一：定义PDE约束 第二：定义目标函数 第三：编写优化过程 管道流问题 初始形状参考链接的geo文件。 参考链接添加链接描述 不知道为什么原链接代码无法收敛，会报错，因此本人做了一些修改 本人把u∇uu\nabla uu∇u去掉了 import firedrake as fd import fireshape as fs import ROL import fireshape.zoo as fsz from fireshape import PdeConstraint from
seigen:基于有限元法的地震波弹性方程求解器
05-14
是基于有限元框架的弹性波方程求解器。 它是地震成像项目的一部分。 快速开始 Seigen需要安装Firedrake，并且必须在Firedrake虚拟环境中运行。 要首先安装Firedrake，请按照以下。 一旦安装了Firedrake并激活了虚拟...
基于FreeFEM++的有限元编程--3
love code, love life.
11-14 5079
目录 前言 1、泊松方程的描述 2、结果保存和读取 3、更多网格生成方法 4、边界条件介绍 5、总结， 前言 FreeFem++是通过有限元方法求解PDE的免费软件，该软件由法国第六大学开发，可以运行在windows、linux和Mac系统上。可以进行网格生成、自动建立有限元模型、整合了各类线性求解器、能适用于与二维和三维问题，可以生成图片、文字和文件作为计算结果的输出。使用方法是编写脚本文件、保存为后缀为.edp为后缀的文本文件，运行软件调用该脚本文件即可。 用网上找到的材料，通过六类方
有限元开源代码FreeFEM++
05-18
有限元开源代码FreeFEM++,个人感觉是基于OOFEM开发的开源代码。
matlab 有限元工具箱
chenyuanyuan的专栏
07-21 2559
开源工具箱
用Python有限元框架Feon求解弹簧支座二维杆问题
Iplay4FuN的博客
11-23 5013
1、Feon简介 Feon是湖北工业大学土木建筑与环境学院教师裴尧尧基于Python开发的一个开源免费的有限元计算框架。这是一个致力于有限元编程教学和有限元理论研究的框架，上手快，可扩展性强，代码开源。 2、安装 使用pip安装 pip install feon，或者源码安装 python setup.py install 3、使用feon求解如下图所示二维杆系统 Feon使用非常简单，只包含...
通用有限元软件简介及比较
zhouschina的专栏
07-04 8245
一般而言,有限元软件属于科技计算软件,评价标准是为3个:(1)可移植性,能够在各个平台之间传输;(2)可使用性——使用简便,可靠;(3)可维护性——可读,易检验,可扩充。当今世界有限元软件舞台上,最为著名的三个通用软件是:ANSYS,ADINA, ABAQUS,下面进行一一介绍以及之间的比较。 三大通用软件 ANSYS: ANSYS是最早进入国内市场的商业化软件之一,也是现下在国内知名度最高的
计算不可压缩流体- NS方程求解算法
热门推荐
lizhengjiang的专栏
03-20 1万+
在上一篇介绍差分格式可以看到求解不可压缩NS方程，主要是求解速度的动量方程和压力的波松方程。本文介绍一些历史中的算法，下一文介绍目前常用的算法。 一  Marker and Celler （MAC) 算法 在自然差分格式中给出PPE如下：    laplace ( p) = laplace ( body_force) - rho * div ( u * grad_u)    （1） 。
针对neumann边界条件的差分法代码
forrestguang的博客
10-11 4931
需要添加人工的Dirichlet边界 考虑的问题仍然是矩形区域的泊松方程。 如果整个区域都采用Neumann边界条件，即∂ΩN=∂Ω\partial\Omega_{N}=\partial\Omega∂ΩN​=∂Ω，则这个问题没有定解，因为对于某个精确解uuu而言，针对任意的常数C，u+cu+cu+c都是解 为此，我们常常需要添加一个人工提供的Dirichlet边界条件，比如说在这里，我们选取{−1}×[−1,1]\{-1\}\times[-1,1]{−1}×[−1,1]为Dirichlet边界，也就是说矩
五点差分法求解泊松方程-共轭梯度法和G-S迭代法
forrestguang的博客
12-30 3062
差分法求解矩形区域椭圆方程 Δu=f,x∈Ω=[−1,1]∗[−1,1]\Delta u = f, x \in \Omega = [-1,1]*[-1,1]Δu=f,x∈Ω=[−1,1]∗[−1,1] u=g,x∈∂Ωu = g, x \in \partial \Omegau=g,x∈∂Ω 在x轴和y轴上分别取分划长度为dx,dy，令M=int(2/dx)+1,N=int(2/dx)+1M = int(2/dx) + 1,N = int(2/dx) + 1M=int(2/dx)+1,N=int(2/dx)+
一维完全气体的Riemann问题求解
forrestguang的博客
05-04 2446
一维完全气体欧拉方程的引入 在这里我们将利用LF格式，MacCormack格式和Roe格式对一维完全气体欧拉方程进行数值求解。 {(ρρuE)t+(ρuρu2+pu(E+p))x=0,p=(γ−1)(E−12ρu2),γ=1.4.\left\{\begin{array}{c} \left(\begin{array}{c} \rho \\ \rho u \\ E \end{array}\right)_{t}+\left(\begin{array}{c} \rho u \\ \rho u^{2}+p \\ u(
区域分解与并行计算
forrestguang的博客
05-27 1983
区域分解算法 以下都以泊松方程来举例 问题引入 给定二维空间的泊松方程： −Δu=f,x∈Ω-\Delta u = f,x \in \Omega−Δu=f,x∈Ω u=g,x∈∂Ωu = g,x \in \partial \Omegau=g,x∈∂Ω 我们把区域分成若干部分分别求解，然后相邻区域做一个通信，这样就可以得到一个并行算法。下面我们先把区域分成两部分做一个介绍： 考虑区域Ω1,Ω2\Omega_1,\Omega_2Ω1​,Ω2​ 其中红色的是Ω1\Omega_1Ω1​，蓝色的是Ω2\Omega
有限元方法求解二维矩形区域椭圆方程
forrestguang的博客
09-06 1915
利用Galerkin方法 def UU(X, order): temp = 10*(X[:,0]+X[:,1])**2 + (X[:,0]-X[:,1])**2 + 0.5 if order[0]==0 and order[1]==0: return torch.log(temp) if order[0]==1 and order[1]==0: return temp**(-1) * (20*(X[:,0]+X[:,1]) + 2*(X[:,0]-X
有限元法流体力学python
最新发布
09-28
有限元法是一种常用于数值求解流体力学问题的方法。它通过将连续领域离散化为有限数量的子域，并在每个子域上使用适当的基函数来逼近解。在有限元法中，流体力学问题的域被分割成小的几何单元，如三角形或四边形。然后，在每个单元上，通过构造适当的基函数来表示未知解的近似值。这些基函数通常是多项式函数。通过将这些基函数的线性组合代入原始的流体力学方程中，并确定基函数系数的值，可以得到一个代数系统，这个代数系统可以进一步用于求解流体力学问题的解。 在Python中，有许多开源的库可以用于实现有限元方法求解流体力学问题。其中一些库包括FEniCS、DOLFIN、Firedrake等。这些库提供了丰富的功能和工具，可以帮助用户构建和求解流体力学问题的有限元离散化。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Galerkin码农选手

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值