ambu1230的博客

原创 MATLAB非均匀网格梯度计算

在matlab中，gradient函数可以很方便的对均匀网格进行梯度计算，但是对于非均匀网格，但是gradient却无法求解非均匀网格的梯度，这一点我之前犯过错误。我之前以为在gradient函数中指定x，y等坐标，其求解的就是非均匀网格梯度了，然而并不是。于是，今天下午开始写非均匀网格求梯度的函数。首先，函数的要求为：1、边界处采用二阶偏心差分2、内部网格点采用二阶中心差分3、计算三维矩阵的梯度。

原创 Tecplot导出流场Movie

本人最近想利用Tecplot导出流场计算的视频，找了以下两种方法：1、直接一次性打开所有文件，导出视频；2、利用脚本每次打开一个文件，导出其照片，最后合成视频。

原创 MATLAB如何批量更改文件名

在openfoam中，由于timePrecision设置的不合理，常常导致计算得到的文件夹名称的数字末尾多很多零，对于读取数据来说十分的不方便。如何更改为标准的，精确的，writeInterval的整数倍名称，需要一下三步操作。

原创 OF小程序（自用）

读取OpenFOAM中每个文件夹的时刻的函数，每次都用的时候都找不到，这次放在网上。

原创 湍流系统的混沌性质

同样的系统，边界条件一致，初始条件一致，采用不同的核并行，得到的系统的响应参数随着时间推进逐渐开始出现分叉，原因是由于细小的误差被不断放大。But，系统的统计特性依然是不变的。

原创 MO尺度(大气边界层)

但是Liu将上面的热通量Q采用湍流热通量代替，即文中说是温度脉动和速度脉动的协方差，这个操作来源于Kapil Boundary-Layer Meteorol 2013，这个存在两个问题，这里插一句，这里热膨胀系数等于平均温度的倒数，我实验了一下，确实是，但是其原因并不清楚，从课程PPT中也看不出来，太菜了，，，而且，其没有出现热膨胀系数，而是用平均温度来代替，KADER 1990 JFM的说法来看，这两者是等价的。是等于0 的，但是对于大气边界层这类复杂的问题，显然这一项不等于0，那他们为啥还要这样做。

原创 Channel Flow 区域图

原创 Linux终端后台运行Matlab程序

在终端下运行matlab程序时，常见的作法是直接进入matlab，然后run程序，但是一旦电脑断网，终端断联，那么程序就需要重新运行，对于有些非常耗时的程序，这时重新运行的时间成本就很大。还有一点，执行nohup，后台会输出一个nohup.out，这个与log文件内容是一样的。需要注意的是 命令结尾需要加上这一个指令&1与上面的

原创 linux shell 脚本基础知识

解释错误，这是因为脚本格式为dos，需要转化为unix,在vim 里打开，输入：set ff =unix即可。需要对其进行chmod +x xxx是脚本具有执行权限；运行shell脚本时在文件名前加./bash不支持浮点数运算。

原创 Tecplot绘制涡结构（Q准则）

Tecplot绘制涡结构(Q准则判别)并用温度进行染色。

原创 2D槽道流

之前看槽道流时，一直无法在二维槽道流里计算出湍流状态，后来了解到二维槽道流需要额外添加随机扰动，但是这个体积力的植入方式一直不知道。施加定长体积力or固定压差能不能产生2d槽道湍流？

原创 从视频提取音频的方式

缺点：视频超出500MB会得付费升级。

原创 AI生成图片教程（基于DALL-E3）

今天登录GPT时发现openai的官网显示其有生成图片的模型DALL-E3，于是想试一试其效果如何。奈何ChatGPT只能在付费版上使用，但是这个DALL-E3其实免费使用。

原创 瑞丽伯纳德会议视频

链接

原创 cmd命令快速打开MATLAB

唯一的缺点是无法使用plot，这一点比不上linux系统，不过打开速度还是挺快的。

原创 There are not enough slots available in the system to satisfy the 48 slots报错

多核运行时出现这个错误，减少核数运行正常输出命令添加执行命令解决。

原创 Tecplot绘制三维彩色流线

用温度云图对流线进行染色。

原创 Paraview绘制体积渲染云图

原创 槽道流CFD结果与理论解对比

这个图来自JFM论文，其中ST是标准槽道流，可以看出其与直线的误差也挺大，但是这里的k和B似乎是两组参数的组合。不知道是笔误还是就是这样。CFD有两组数据，一组是pisoFoam跑的，摩擦雷诺数精确为180.可以看出蓝线更接近实验标定值，红线有点偏离，不知道这个误差能不能接受。对于槽道流，理论上近壁面时均速度为线性分布，远离壁面处为对数分布。在对数区，有两个未知参数，两个人标定了不同的两组数据。图中的线采取的参数为k=0.40，B=5.5.可以看出还是有一定区别的。两组数据网格不太一样。

原创 buoyantBoussinesqPimpleFoam代码解析

buoyantBoussinesqPimpleFoam求解器只是在pimpleFoam求解器的基础上添加了布西内斯克近似的浮力项，感觉只要把pimpleFoam求解器搞明白，在pimpleFoam求解器的基础很容易修改处理，由于该求解器只在低版本的OF中存在，高版本的OF中已经融合到buoyantPimpleFoam中，但是由于该求解器比较复杂，所以我们还是得到buoyantBoussinesqPimpleFoam比较简单。

原创 OpenFOAM运动壁面边界条件

链接

原创 如何通过图片提取文献图里的数据点

即可显示具体坐标，这里我们查看原始的数据点为（0.1，1.14977）精度还可以。点击图片中你想知道的点，选中后会显示红色圆点，然后点击左侧的view data。标定是按照X1,X2,Y1,Y2的顺序来的，完成之后点击complete!点击右侧的Add Point可以新增数据点。可以勾选对数坐标，完成之后点击OK。

原创 OpenFOAM颗粒求解器中的integrationSchemes

隐式欧拉：

原创 CFDEM学习笔记

虚拟机：链接：https://pan.baidu.com/s/1MPMTJQfl76mW0H5bbT_rAg提取码：rqli开机密码：530944988知乎博客：作者说明了如何关闭颗粒碰撞计算来达到提升计算速度。Github：这个是github上的CFDEM，在原来的基础上增加了许多求解器。官网：官网教程，但是比较鸡肋。

原创 解决OpenFOAM颗粒计算输出文件Paraview无法打开问题（二）

第二个方案的源是在CFD中文网上看到的一篇帖子，其具体链接忘了。这个帖子给了一个github的链接，就是将OpenFOAM输出的颗粒位置信息转变为真实的位置信息的脚本。。

原创 解决OpenFOAM颗粒计算输出文件Paraview无法打开问题（一）

长话短说，关于这个问题，有两个解决方案，一是将文件类型转换为VTK格式，二是采用github上的脚本，将OpenFOAM输出的位置文件转换成真实的颗粒的真实位置文件。

原创 CFDEM里coupingProperties文件解析

CFDEM里，在进行两相耦合求解时，需要在CFD/constant文件夹下设置coupingProperties文件，本文给出这个文件里一些参数的意义。

原创 模拟退火算法，遗传算法，禁忌搜索算法的特点

（1）借助物理学中退火的思想，从某一高温出发，随着温度参数不断下降，在解空间中寻找目标函数的全局最优解，温度影响着当新解不优于当前解时，接受新解的概率，温度越高，接受新解的概率越高。（2）基于概率的算法（3）需要设置如何产生新解。

原创 OpenFOAM的fvOptions

对此，我在进行槽道流仿真时，进行了验证。对比了中间截面的平均速度和出口截面的平均速度，发现它们几乎一致。这个fvOptions影响的是方程中的源项，类似重力项，通过调整体积力的大小来设定速度的大小。采用OpenFoam中的fvOptions。进行如下设置，可以大致设定流动的平均速度。

原创 CFDEM-OpenFOAM-Yade安装教程

在网上搜索与OpenFOAM相关的颗粒两相流计算资料时，发现了一个。在此之前，我学习过OpenFOAM自带的颗粒计算求解器，但是自带的求解器有很多缺点，最大的缺点就是颗粒运动方程的求解精度只有一阶。还学习过CFDEM，是OF+liggghts做联合仿真，这个liggghts学习起来很简单，很多教程可以看，但是CFDEM就没这么多教程了，唯二的一个交流群里的几节入门课程买到了1k，壁垒拉的满满的。这个我简单的跑过一个案例，但是由于参数太多，最终没有继续学下去。

原创 集群安装OpenFOAM

在个人电脑上安装比较简单，在服务器上安装多出几个步骤，需要首先安装boost库以及openmpi库下面步骤主要参考，这里我自己重复一遍。

原创 MATLAB程序初始化OpenFOAM颗粒位置

在OpenFOAM的颗粒两相流求解器中，我们可以采用manualInjection的方式进行自定义颗粒的初始位置，这个命令十分方便，在CFDEM中也有类似的命令，不过CFDEM中的命令更加强大，我们不仅可以定义颗粒的初始位置，而且还可以定义颗粒的初始速度,唯一的缺点就是CFDEM学起来资料比较少。于是自定义粒子的位置的思路就非常简单，我们利用MATLAB创建上面的文件内容即可。我们把上面的坐标粘贴到kinematicCloudPositions文件里即可。这里不对程序做解说，直接放在这里。

原创 DPMFoam给粒子添加周期性边界条件

DPMFoam是OpenFOAM中的一个颗粒两相流求解器，在流体的运算中，我们经常采用周期性边界条件，在OpenFOAM中即为cyclic，那么对于颗粒，能否在边界也设置成周期性边界呢，在CFDEM中是可以的，因为liggghts的设置项中即有颗粒的周期性边界条件设置，p代表周期性边界。其中rebund为反弹边界，stick为吸附边界，escape为逃逸边界，这几个都有比较明确的意义，但是没有出现周期性边界。首先我们需要设置流体区域的上下表面为cyclic边界条件，然后设置粒子的相应边界为none。

原创 小白电脑装机(自用)

我这边曾将内存条和主板一块寄给维修人员，得到的回复是有一个内存条无法点亮，这里我们认为主板和CPU是没问题的。但是我拿到之后，仅有那个好的内存条还是无法点亮，一直卡在DRAM灯这里，令人崩溃，目前我已经排除掉第一种可能，因为内存条已经反复试过，并且也用橡皮擦擦过。这里的第二点目前尚不明确，因为商家都点亮了，说明还是适配的，但是我目前还无法查证哪里可以看到主板和内存条是否适配的资料。目前是等着新买的内存条到货，然后测试，没办法，目前只能这样了，后续更新进度。按照主板简介给出的信息应该是适配内存条的。

原创 Paraview调出左侧工具栏

勾选这两项即可调出。

原创 VM里ubuntu虚拟无法启动

和我的错误一模一样，都是只有一个光标在闪。于是按照这个操作进行了一遍，发现是home文件满了，于是删除了一些东西，然后reboot，成功启动。关闭的3d加速图像显示，以及那个cmd命令，但是没什么用。开始认为是VM的设置问题，按照。

原创 paraview显示颗粒

后处理方法

原创 Edge浏览器和Google浏览器占用内存情况

说edge浏览器确实占用内存较多。但是给出的解决方案是更新windows系统，这个方法显然不能接受，采用google浏览器即可减小内存占用。最近发现适用edge浏览器看视频时，或者打开多个网页，电脑就会变得非常卡顿，打开任务管理器发现内存占用较多，google浏览器内存占比比edge浏览器少百分之十左右，因此建议采用google浏览器看视频。开始之前，内存占用率稳定在55%-57%开始之前，内存占用率稳定在55%-57%1、Edge浏览器占用内存情况。2、Edge浏览器占用内存情况。下面进行一组对比实验。

原创 2、瑞丽-伯纳德对流的拉格朗日拟序结构（FTLE场结果对比）

在中，我简单比较了瑞丽伯纳德对流的FTLE场，但是因为粒子追踪采用的是欧拉方法，所以精度不是很高， 因此与文献中的结果还是有些差别。下面放一张文献中的FTLE场，参数与上篇文章是一致的，Ra = 1e8;Pr = 1;本文将详细的介绍如何得到这样的图片，好吧，其实也不是很详细。总之，先展示以下我个人的结果。可以看出，上图基本上与文献是一致的，因为这里展示的是瞬时场，所以当然不可能是完全一致，至于提取的LCS结构，我在下面也会提及。本着详细介绍的原则，所以这里我很良心的从流场计算开始介绍。

原创 麻了。。。

我目前跑算例发现长时间步保存数据得到的结果不太对，跟文献差的有点多，其实也就不到百分之十，但是不同网格之间的差距虽然小一点，但是也还是有，细网格更接近文献，我感觉有种肯是时间推进步长的问题，再小些似乎会更好，但是。可以看出其实是不太一样的，一个取平均是24.72，一个是25.34，差距其实有的。同样的初始条件，一个每隔0.05s保存一步数据，一个每隔1.13s保存一步数据。有点希望是这个原因，又不是很希望是这个原因。上面横轴代表时间，纵轴代表Nu。湍流真是一个混沌的系统。