- 博客(19)
- 资源 (2)
- 收藏
- 关注
RSS订阅
原创 OpenFOAM理论笔记 目录
第一章 有限体积离散1.1 计算域的空间离散1.2系统的控制方程1.2.1流动与传热的通用控制方程1.2.2 通用的标量输运方程1.2.3 控制方程的分类1.3 控制方程的离散1.3.1 高斯定理1.3.2 标量输运方程的离散1.3.3 离散后方程的求解第二章 压力速度耦合算法2.1 压力泊松方程2.2 SIMPLE系列算法2.2.1标准SIMPLE算法2.2.2 SIMPLE算法的亚松弛2.2.3 SIMPLEC算法2.2.4 OpenFOAM中..
2020-08-04 21:31:48
2606
3
原创 OpenFOAM分段边界条件的设置
问题描述在实际使用OpenFOAM中，我们经常会遇到根据平面上点的坐标给出边界条件的情况，例如一个三维方腔，在距离墙壁底端3m的位置上有一个直径0.5m的管道入口，流体从其中以恒定的流量流入计算域。在这样的背景下，我们就需要对这个壁面给定分段的边界条件。在本文中我们就具体讨论一下如何设置这样的边界条件。具体步骤我们用到的两个工具是creatPatch和topoSet。 creatPatch需要在算例的system文件下调用creatPatchDict文件，其原本用来调整周期性边界条件的位置，在这里我
2020-10-23 17:50:12
4184
1
原创 3.2.1 对流项离散格式的基本介绍 | 3.2.2 中心离散格式与一阶迎风格式（OpenFOAM理论笔记系列）
3.2 对流项的离散格式3.2.1 对流项离散格式的基本介绍在开始本节的讨论开始，笔者首先要说明，本节所介绍的内容相对于整个对流项离散格式的开发的历史可以说是极其简略的。可以这么说，整个计算流体动力学发展的历史有一半以上是由对流项离散格式写成的。本节的内容以Versteeg和Malalasekera的著作《An Introduction to Computational Fluid Dynamic: The Finite Volume Method Second edition》中的对流项离散格式介绍内
2020-08-28 21:52:44
5263
原创 3.1 扩散项的离散格式（OpenFOAM理论笔记系列）
第三章 离散格式3.1 扩散项的离散格式3.1.1 面法向梯度的计算有限体积法中，物理量都存储在控制体体心上，而由于我们在进行有限体积离散时使用了高斯定理，因此在最终组建代数方程时，我们往往都用的是控制面上的值。在这里，就牵扯倒一个如何根据体心值去计算面上值的过程，这就是离散格式要解决的问题。对于一个标量ϕ\phiϕ,按照1.3.2节的方法，我们有：∫VP∇⋅(Γ∇ϕ)dV=∑f[Sf⃗⋅(Γ∇ϕ)f]=∑f[Sf⃗⋅Γf(∇ϕ)f](1.29)\int_{V_P}\nabla\cdot(\Gam
2020-08-13 16:31:24
4691
9
原创 2.4 PIMPLE算法 | 2.5 附加显式力的压力速度耦合（OpenFOAM理论笔记系列）
2.4 PIMPLE算法2.4.1速度的非线性耦合在2.1节中，我们曾对流动系统的压力速度耦合问题进行过讨论，实际上，在流动系统中除了压力和速度的耦合关系，我们还存在着另一种重要的耦合关系，速度的非线性耦合关系。当我们对速度进行求解时（例如求解动量预测方程），速度是未知量，因此对流项∇⋅(U⃗U⃗)\nabla\cdot(\vec U\vec U)∇⋅(UU)是一个未知量乘以未知量的结果，这样会使得原来呈线性的方程变为非线性方程（即未知数的最高次数从1变为2）。直接求解非线性系统的矩阵方程需要调用非线性
2020-08-04 21:28:04
6147
2
原创 2.2 SIMPLE系列算法 | 2.3 PISO算法（OpenFOAM理论笔记系列）
2.2 SIMPLE系列算法2.2.1标准SIMPLE算法SIMPLE算法(Semi-Implicit Method for PressureLinked Equations)1最初被设计用来求解稳态问题，即控制方程中不包含瞬态项的计算。按照1.3.3节的约定，我们假设计算开始的时候有初始的压力和速度值Po,Uo⃗P^o,\vec{U^o}Po,Uo,待求的真值为Pn,Un⃗P^n,\vec{U^n}Pn,Un。实际上，速度的初始值可以完全随意给定，因为我们完全首先可以将初始压力带入式2.1解出一个速度
2020-07-31 21:46:16
15826
4
原创 2.1 压力泊松方程（OpenFOAM理论笔记系列）
第二章 压力速度耦合算法2.1 压力泊松方程在第一章中，我们推导出了不可压缩牛顿流体的控制方程：连续性方程∇⋅U⃗=0(1.8)\nabla \cdot \vec U =0\tag{1.8}∇⋅U=0(1.8)动量方程∂U⃗∂t+∇⋅(U⃗U⃗)=∇⋅ν∇U⃗−∇pρ+g⃗(1.15){{\partial \vec U}\over{\partial t}}+\nabla\cdot(\vec U\vec U)=\nabla\cdot\nu\nabla\vec U-\nabla {p\o
2020-07-25 21:28:10
11013
1
原创 1.3 控制方程的离散（OpenFOAM理论笔记系列）
1.3 控制方程的离散1.3.1 高斯定理在有限体积法中，离散的核心数学公式为高斯定理：∫V(∇⋅a⃗)dV=∮∂VdS⃗⋅a⃗(1.20.a)\int_V(\nabla\cdot\vec a )dV=\oint_{\partial V}d\vec S\cdot\vec a\tag{1.20.a}∫V(∇⋅a)dV=∮∂VdS⋅a(1.20.a)∫V(∇ϕ)dV=∮∂VdS⃗ϕ(1.20.b)\int_V(\nabla\phi )dV=\oint_{\partial V}d\vec S\
2020-07-19 20:15:29
4471
7
原创 1.2.2 通用的标量输运方程|1.2.3 控制方程的分类（OpenFOAM理论笔记系列）
1.2.2 通用的标量输运方程1.2.1节涉及的动量方程和能量方程都可以整理为如下的标准的标量输运方程：∂(ρϕ)∂t+∇⋅(ρϕu⃗)=∇⋅(Γ∇ϕ)+Sϕ(1.16)\frac{\partial(\rho\phi)}{\partial t}+\nabla\cdot(\rho\phi\vec u)=\nabla\cdot(\Gamma\nabla\phi)+S_\phi\tag{1.16}∂t∂(ρϕ)+∇⋅(ρϕu)=∇⋅(Γ∇ϕ)+Sϕ(1.16)式(1.16)从左至右依次是表示控制体
2020-07-19 20:10:50
3815
3
原创 1.2.1流动与传热的通用控制方程（OpenFOAM理论笔记系列）
1.2系统的控制方程1.2.1流动与传热的通用控制方程通用的流体流动与传热方程如下：连续性方程：∂ρ∂t+∇⋅(ρU⃗)=0(1.4){{\partial \rho}\over{\partial t}}+\nabla\cdot(\rho \vec U)=0\tag {1.4}∂t∂ρ+∇⋅(ρU)=0(1.4)动量方程：∂(ρU⃗)∂t+∇⋅(ρU⃗U⃗)=∇⋅τ−∇p+ρg⃗(1.5){{\partial (\rho\vec U)}\over{\partial t}}+\nabl
2020-07-07 21:42:30
6230
1
原创 1.1 计算域的空间离散（OpenFOAM理论笔记系列）
1.1 计算域的空间离散导语： 如果有人问你有限体积法的精度是几阶，你也许会不假思索地回答：“二阶!” ，但是为什么是二阶呢？使用更高阶的离散格式能不能使得有限体积法的精度提升呢？认识有限体积法，我们先从认识有限体积入手。相比有限体积这个说法，更多的人可能更熟悉“网格”这个词。在一张白纸上画几道横线，再画几道竖线，你就能得到一个最简单的网格系统。在你的第一印象里，也许有横线和竖线构成的一个小四方块就是一个网格单元，也就是我们说的一个有限体积（即控制体，之后均采用控制体一词）。然而事实并非如此。图1.
2020-07-07 21:37:35
3600
8
原创 gnuplot常用命令(随时更新)
最后更新日期：202006141.已有数据的点线图plot '.\line_rhoe1.xy' u 1:2 w lp lw 2 lt 2 pt 1 lc rgb "black" ,\'.\line_rhoe2.xy' u 1:2 w l lw 2 lc rgb "black" #u 1:2 使用第一列数据为x轴，第二列数据为y轴 还可以进行运算 例如 u 1:($2)/($3)#w lp 绘制数据点和数据线 w l 仅绘制数据线 w p仅绘制数据点#lt/pt 2 点和线的绘制风格 lt
2020-06-14 20:34:32
1927
转载 OpenFOAM 工具列表
OpenFOAM工具可以在终端输入命令 cd $FOAM_UTILITIES查看，标准小工具列表如下： Pre-processing 前处理 applyBoundaryLayer Apply a simplified boundary-layer model to the velocity and turbulence fields based on the 1/7th power-law applyWallF...
2020-05-29 21:36:59
4118
原创 interFoam进行两相流模拟的一些典型参数设置对比¶view绘制相界面的方法
1 背景我的个人课题是基于interFoam的一些二次开发,然而我在使用自己的求解器进行模拟的时候往往会因为模拟参数配置的不合适导致这样或那样的一些问题。这些问题本质上是自己对interFoam进行模拟时的一些关键参数不太了解导致的。前几天看到一本资料上关于interFOAM的讲解，其中的拓展练习中包括一些不同模拟参数的对比，借此机会我就使用溃坝算例对一些典型参数对模拟的影响进行了实际的试算分析,也将之前积累的一些心得一同分享出来。2 算例二维溃坝算例，来自OpenFOAM7的自带的教程。将模拟参数延
2020-05-29 17:21:06
5388
1
原创 icoFoam代码解析
icoFoam代码解析icoFoam是OpenFOAM最基础的求解器之一，想要以OpenFOAM为基础进行二次开发的同学必须对icoFOAM的代码有所了解。关于此部分内容，由于篇幅限制，更多内容可以参考本博文对应的资源"OpenFOAM编程日志，单相不可压缩流动"或去本人的CSDN资源空间搜索下载。/*----------------------------------------------...
2020-04-10 20:21:59
3411
3
原创 简单自编程实现手动划分OpenFOAM并行分区（以manual方式进行划分）
最近在进行计算的时候，为了提高并行效率，我决定手动进行OpenFOAM并行网格的划分。一般来说，提高并行效率可以从两个方面着手：使每个CPU负担的网格数量尽可能相等，避免出现限制计算的短板。使每个CPU负担的网格之间的公用面尽可能少，减少CPU之间的信息交互。一般来说，如果网格比较简单，可以轻易使用simple方法划分，如果网格比较复杂，可以用scotch方法划分（关于各划分方法可以查看...
2020-04-06 20:50:59
3688
7
原创 一些OpenFOAM自带的有用小工具
前言这篇博文总结自OpenFOAM的用户手册。OpenFOAM的用户手写得很好，初学OpenFOAM的时候没觉得什么，等差不多入门之后，回头再看一遍，能发现其中很多亮点，很多能解决实际使用痛点的操作。有空多读手册，多回头看看。（本文长期更新）场映射 mapFields插值得到中间步的场 temporalInterpolate重新进行并行区域划分redistributeParOpenF...
2020-04-05 12:55:08
5593
1
原创 OpenFOAM后台运行求解器的多种方式
OpenFoam后台运行求解器的多种方式严格说，这不是OpenFOAM的使用技巧，而是Linux系统的运行程序的技巧。最简单的后台运行一个求解器可以通过在终端输入求解器名 &的方式进行，例如：icoFoam &注意，后台运行的命令不是不在控制台输出信息，而是指在求解器运行的时候还允许用户在同一个终端输入其它命令。这对于一般的求解器运行没有什么意义，因为OpenFOAM的计...
2020-04-04 20:25:13
2057
1
原创 OpenFOAM中重力的植入方式
OpenFOAM中重力的植入方式考虑重力的NS方程可以写为：ρ∂u⃗∂t+ρ∇(u⃗u⃗)=∇(μ∇u⃗)−∇P+ρg⃗(1)\rho \frac{\partial \vec u}{\partial t}+\rho \nabla(\vec u \vec u)=\nabla(\mu\nabla \vec u)-\nabla P +\rho \vec g \tag {1}ρ∂t∂u+ρ∇(uu)...
2020-03-25 18:59:16
3179
OpenFOAM标量场转化为labelList代码.tar.xz
2020-04-06
openFOAM编程日志_单相不可压缩流动2020edition2.5.pdf
2020-03-23
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人