QJ Meng 01-优快云博客

原创简单的MATLAB输出Excel&TXT方法

MATLAB中计算出一组结果，比如F, V, S等，如何写入电子表格中？F_h=0; %水平力Vx1=0; %x速度Sx1=0; %位移 F_h(cycle/SaveStep)=Fx; % cycle是总循环次数，SaveStep是存储间隔 Vx1(cycle/SaveStep)=Vx; Sx1(cycle/SaveStep)=(Sx-Sx01); data1 = [F_h', Vx1',

2020-12-26 16:41:09 1185

原创 MATLAB的动态图，伪色彩图、矢量图、流线图，散点图折线图的绘制例子--流体力学；LBM

还是老规矩先宣传一下QQ群群：格子玻尔兹曼救星：293267908。免费群！一切为了早日毕业。最近群友问画图的挺多，伪色彩图、矢量图、流线图，我在这里贡献一下自己的画图代码算是给大家提供参考。1 流线图。用任何语言生成的xy坐标和uv速度场，怎么画流线图和矢量图呢，MATLAB提供streamslice函数：%From https://ww2.mathworks.cn/help/matlab/ref/streamslice.html hh=streamslice(ux',u

2020-05-22 12:57:55 10949 14

原创 TECPLOT画流固耦合的流场图--学习笔记

TECPLOT画流固耦合图片–学习笔记（参考[link]https://www.cnblogs.com/zhubinglong/p/8735426.html）还是老规矩先进下QQ群：格子玻尔兹曼救星 293267908。一 .dat文件Title可以不写。通常流体作图的抬头F=POINT，I和J分别是点的个数与元的个数，LBM的流场中，IJ分别写你的xy轴的流体网格数量。zone t=“1”，zone t=“2”…的设定非常有用，可以用来设定不同的区用不同色彩显示，这里分别用来设置流场和固体场，然后

2021-05-08 21:01:34 3538 1

原创涡线图用tecplot咋画

步骤：Zone style-contour-contour type-lines直接这样操作是一个颜色的。还可用颜色区分不同等值线：

2021-04-11 15:43:19 3269

原创非平衡反弹格式更正-[王富海. 基于MRT-LBM的流场与声场仿真计算[D]. 华中科技大学, 2017.]

先说下我们qq群：格子玻尔兹曼救星：293267908。之前分享过这个文章的代码，而论文中有一些问题更正如下。感兴趣的可以推一下，推导的方法论文中写的很详细了。

2021-03-20 18:13:12 1073 1

原创关于反弹操作

obst = (x-obst_x).^2 + (y-obst_y).^2 <= obst_r.^2;%圆球obst(:,[1,ly]) = 1;%边界处理bbRegion = find(obst);fOut(i,bbRegion) = fIn(opp(i),bbRegion);%反弹式边界这样的写法中，bbRegion指的是标记了格子空间重新排列为一维数组后，障碍物的格子的下标。除了圆柱边界，圆柱内部的格子，也执行了回弹边界的操作。...

2020-07-09 10:47:17 826 1

原创 MATLAB图像二值化

还是老规矩先宣传一下QQ群群：格子玻尔兹曼救星：293267908。一切为了早日毕业。听说可以在PPT作图然后用MATLAB二值化把复杂边界用01矩阵表达出来。上网看了下：I = imread('18_4_25.jpg');imshow(I);thresh = graythresh(I); %自动确定二值化阈值I2 = im2bw(I,thresh); %对图像二值化imshow(I2);图片是我18年拍摄的北洋园校区：二值化之后让我有一种想裱起来的冲动

2020-06-29 15:09:08 19999 4

原创敲黑板！半步(half-way)反弹格式和标准反弹格式的随堂笔记

先说下我们qq群：格子玻尔兹曼救星：293267908。每个人都是自己的救星。铁打的群流水的学生，总有人问的问题，还是得记录一下。郭照立：形式上完全相同，是否编程完全相同？实际上half_way是碰撞之后立马反弹，也就是写入反弹的时候坐标不变。在边界上，直接令 f(k,x,y) = fpost(inv(k),x,y); 其中fpost(inv(k),x,y)是上一步走向固体的分布函数，f(k,x,y) 则是当前由固体走向流体的分布函数，也就是一步就弹回来了，你的xy坐标没有体现固体的位置

2020-06-25 16:26:58 2370 2

原创 Tecplot导出来的图片空白太大了怎么办？

有没有发现，tecplot调整好宽之后，高度方向总是留着很多空白？然后你是不是打开AI开始手动裁剪了？下面告诉大家调整的方法。哦对了老规矩，先宣传一下QQ群群：格子玻尔兹曼救星：293267908。我们的宗旨是：一切为了早日毕业。上面是计算结果的图。留白难看。打开frame：这个frame可以保存下来以后再用。调整轴：大功告成v~...

2020-06-19 21:25:45 10973

原创 MATLAB，几种让你储存计算数据的方法

大家好，又是我。老规矩先宣传一下QQ群群：格子玻尔兹曼救星：293267908。我们的宗旨是：一切为了早日毕业。1、最简单的还是生成一串0，然后挨个把每个时间步的值存进去。%初始设置：Vx1=zeros(1,timesteps);Vy1=zeros(1,timesteps);%程序运行到cycle步Vx1(cycle)=Vx;Vy1(cycle)=Vy;%Vx是当前速度2、不妨写出到TXT：%设置位置： xy coordinate datadata_x1 = fope

2020-06-19 20:16:40 3413

原创 SPALART–ALLMARAS TURBULENCE MODEL--Sa大涡的方腔流代码--感谢Robert Lee

还是老规矩先宣传一下QQ群群：格子玻尔兹曼救星：293267908。不收费的哦，就是为了早点毕业建的群。%文章来源：AN IMPLEMENTATION OF THE SPALART–ALLMARAS TURBULENCE MODEL IN A MULTI-DOMAIN LATTICE BOLTZMANN METHOD FOR SOLVING TURBULENT AIRFOIL FLO...

2020-03-02 21:07:54 1138

原创 Splart-Allmaras湍流模型及MATLAB编程~

还是老规矩先宣传一下QQ群群：格子玻尔兹曼救星：293267908。不收费的哦，就是为了早点毕业建的群。代码完整版在群里lbm_matlab-master。文章请搜索：(2015) Nicolas Pellerin, Sebastien Leclaire, and Marcelo Reggio. AN IMPLEMENTATION OF THE SPALART–ALLMARAS TURBUL...

2020-03-02 16:28:53 3371 3

原创物理参数到LBM格子参数的转换方法--注意雷诺数的L取值与粘度单位换算！

还是老规矩先宣传一下QQ群群：格子玻尔兹曼救星：293267908。不收费的哦，就是为了早点毕业建的群。%参数计算: P=C*d, d是dimensionless无量纲值,P是physical参数，C是转化参数 %离散空间和时间尺度 δx和δt ,lb是离散数目. CC=P/lb是从物理参数到格子参数的转换参数lb。deltas=1;deltat=1;lattc=deltas/delt...

2020-02-25 21:37:11 7084 16

原创 Fortran语言--自由界面程序。

还是老规矩先宣传一下QQ群：格子玻尔兹曼救星：293267908。 PROGRAM dambreak USE IFPORT implicit none! //LBM model Real*8, parameter:: w(0:8) = (/4.0d0/9.0d0,1.0d0/9.0d0,1.0d0/9.0d0,1.0d0/9.0d0,1.0d0/9.0d0,1...

2020-01-05 09:12:09 3007

原创加重力的方腔流！

%还是老规矩先宣传一下QQ群：格子玻尔兹曼救星：293267908。%加了重力之后上侧的密度会有逐渐的减轻，因为LBM中，压力和密度是耦合的！% ^y % 7 3 6 % \|/ % 4-1-2 ---> x % /|\ % 8 5 9 %%% INPUTS: rho The flui...

2019-12-22 22:37:00 1012

原创外力项来自Guo Phy.2002 2D immersed boundary lattice Boltzmann method code. Code from Timm Kruger.

// 还是老规矩先宣传一下QQ群：格子玻尔兹曼救星：293267908。// It uses the D2Q9 lattice with Guo's forcing term. // Rigid bottom and top walls are parallel to the x-axis (channel). // The flow is periodic along the ...

2019-12-22 15:54:27 1099 1

原创无滑移边界的水冲小球例子，用浸没式边界（速度差法）求曳力--修改自王富海硕士论文

%这个例子采用 MRT-LBM 模拟旋转圆柱绕流%左边速度边界-泊肃叶流，右边压力边界，上下无滑移壁面（全部用非平衡外推格式）%力的计算来自：基于 MRT-LBM 的流场与声场仿真计算 --王富海2017%对照例题来自：P33，王露. 一种浸没边界格子Boltzmann方法的算法改进及在CFD中的应用[D].% 还是老规矩先宣传一下QQ群：格子玻尔兹曼救星：293267908。...

2019-11-05 08:57:34 1649 2

原创采用 MRT-LBM 模拟旋转圆柱绕流2---MATLAB代码--王富海2017--基于 MRT-LBM 的流场与声场仿真计算

%这段代码之前发过，结束后生成图也都贴出来了，但是很多地方没有写出详细的说明，加上王富海的3.2图做的一塌糊涂，力的计算引用自王星，但是王星的学位论文画图也是字母全标错了，当时看到这里也是欲哭无泪。拜托你们毕业的能不能认真一点。所以在这再发一次，对一些内容进行补充。%还是老规矩先宣传一下QQ群：格子玻尔兹曼救星：293267908。不收费的哦，就是为了早点毕业建的群。%这个例子采...

2019-11-01 10:27:08 3691 2

原创采用 MRT-LBM 模拟旋转圆柱绕流---MATLAB代码--王富海2017--基于 MRT-LBM 的流场与声场仿真计算

%这个例子采用 MRT-LBM 模拟旋转圆柱绕流%基于 MRT-LBM 的流场与声场仿真计算 --王富海2017%左边速度边界-泊肃叶流，右边压力边界，上下无滑移壁面（全部用非平衡外推格式）%还是老规矩先宣传一下QQ群：格子玻尔兹曼救星：293267908。不收费的哦，就是为了早点毕业建的群。clcclearclose all%% 设置仿真参数uM...

2019-10-28 12:02:36 6277 1

原创这个例子采用 MRT-LBM 模拟矩形腔绕流--基于 MRT-LBM 的流场与声场仿真计算 --王富海2017

%这个例子采用 MRT-LBM 模拟矩形腔绕流%基于 MRT-LBM 的流场与声场仿真计算 --王富海2017%上边界速度边界，其它边界-上下非平衡反弹格式-无滑移壁面%还是老规矩先宣传一下QQ群：格子玻尔兹曼救星：293267908。不收费的哦，就是为了早点毕业建的群。clcclearclose all%设置仿真参数uMax=0.1;%中间最大速度xLen=151...

2019-10-25 16:40:18 2112 2

原创力驱动泊萧叶流，基于 MRT-LBM 的流场与声场仿真计算 --王富海2017

%这个例子采用 MRT-LBM 模拟体积力驱动的二维泊肃叶流，%左右周期边界格式，上下非平衡反弹格式.%基于 MRT-LBM 的流场与声场仿真计算 --王富海2017%还是老规矩先宣传一下QQ群：格子玻尔兹曼救星：293267908。不收费的哦，就是为了早点毕业建的群。clcclearclose all%设置仿真参数xLen=60; %水平方向格子数yLen=...

2019-10-25 16:07:45 3015 1

原创方腔流--格子玻尔兹曼方法，基于MATLAB，带画图程序。

还是老规矩先宣传一下QQ群：格子玻尔兹曼救星：293267908。不收费的哦，就是为了早点毕业建的群。% To prove the credibility of fluid part of solver. %方腔流体--速度计算。%来自书--莫哈默德中文 A.A.Mohamad P141%Qiujie Meng 2019.10.24.% ^y % ...

2019-10-24 17:03:09 6432 8

原创郭照立P23页公式推导

感谢云龙同志的贡献。

2019-10-18 16:39:34 1029

原创 D2Q9模型--基于MATLAB的LBM代码

哇咔咔又扒代码啦，感谢我的难兄老田同志贡献了自己的账号。这个代码上面没署名，不知道是哪个朋友写的，侵权立刻删哈。还是老规矩先宣传一下QQ群群：格子玻尔兹曼救星：293267908。不收费的哦，就是为了早点毕业建的群。% ------------------% -> 二维通道，左进右出，上下为壁面% -----------...

2019-09-28 20:42:15 9027 23

原创改变 STREAMING 部分--基于MATLAB的LBM代码

有啥LBM相关问题Feel free to tell us~~毕竟群众力量大嘛~格子玻尔兹曼救星QQ群：293267908。不收费的哦，就是为了早点毕业建的群。之前上传的代码有震荡问题，而且运行时间久了就爆炸。。。昨天重写传播部分，不再使用circshift函数之后就解决了。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%...

2019-09-20 08:19:24 849 1

原创关于LBM-IBM方法中的无滑移边界--Niu X D A momentum exchange-based immersed boundary-lattice Boltzmann method.

感谢群里的朋友们提供帮助！还是老样子，有啥问题Feel free to tell us~毕竟群众力量大嘛~格子玻尔兹曼救星QQ群：293267908。这两天看了大佬Niu X D 的文章（Niu X D , Shu C , Chew Y T , et al. A momentum exchange-based immersed boundary-lattice Boltzmann met...

2019-09-19 09:16:03 1931

原创 An example： 2D immersed boundary lattice Boltzmann method code--By Timm Krüger.2011.

//谢谢Timm Krüger大佬的代码。还是老样子，有啥问题Feel free to tell us~毕竟群众力量大嘛~QQ群：293267908。 // 话不多说，上干货！ // ISBN 978-3-319-44649-3 (Electronic) // 978-3-319-44647-9 (Print) // http://www.springer.c...

2019-09-18 10:07:54 1264 1

原创基于LBM-IBM-DEM的圆形颗粒在粘性流体中沉降的耦合模型及数值模拟[D]. 湘潭大学, 2013. 的解读.

感谢李浩同学的文章。在这里记录几个问题。一、这个公式中，第二个fβ是上一步中的分布函数，而fα与fβ的区别是，以BounceBack边界为例：% 7 3 6 ^y % \|/ | x % 4-1-2 ---> % /|\ % 8 5 9cx= [ 0, 1, 0, -1, 0, 1, -1, -1...

2019-09-06 20:40:53 2093 1

原创基于MATLAB的LBM代码： Rough jet model

%By liu-2017.0403. 谢谢这位没写名字的大佬。%又上网扒代码了。感人的效果。流体流入水腔。%还是老样子，有啥问题Feel free to tell us~毕竟群众力量大嘛~QQ群：293267908。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clear% GENERAL FLOW CONSTANTSlx =...

2019-07-10 15:31:59 2425 5

原创 MATLAB写的LBM管道小球例子

%%%这个程序是有问题的。也就是借鉴一下，用的话需要自己改！！！%%还是老样子, QQ群：293267908。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% cylinder.m: Channel flow past a cylinderical % obstacle, using ...

2019-06-26 10:04:51 2270 5

原创 The Immeased boundary Method示例程序第二弹~~~

%%%一个兄弟找的，还没来得及调试。先分享出来~有需要研究这个的朋友大家一起讨论啊~%%%最近发现了个好东西叫MATDEM，还没来得及详细研究。http://www.matdem.com/list/?94_1.html，开源软件，可做后%%%处理。%%%还是老样子，有啥问题Feel free to tell me, 我会的一定解答！ QQ群：293267908。global...

2019-06-26 08:56:49 564 1

原创基于The Immeased boundary Method~Charles S.Peskin的IBsolver，很好的示例程序。MATLAB写的。

2019-06-25 21:50:52 702

原创矩形障碍算例（附Fortran计算代码及MATLAB后处理代码）

! 我又来啦，做了个小算例，在100*100个格子的空腔中加了几个矩形障碍!上代码！西边界为水流进入区parameter (n=100,m=100)real f(0:8,0:n,0:m)real feq(0:8,0:n,0:m),rho(0:n,0:m)real w(0:8), cx(0:8),cy(0:8)real u(0:n,0:m), v(0:n,0:m)integer i...

2019-04-05 18:08:57 1570

原创如何把笛卡尔坐标系的矢量数据导入paraview？

由我的第一篇博客我们可以知道，在Fortran中可以生成矢量数据，包括每个点的坐标以及该点上xy方向的速度分量。然后的一周我都在探寻如何将这个数据画出矢量图最后画出流线图。今天在网上搜索到一种方法，可以在MATLAB中写出.raw文件，然后再导入paraview，我刚刚试了一下不知道是数据问题还是使用问题，反正显示还不太对。先记录下，下午再回来研究问题的解决方法。来源：https://ww2...

2019-04-02 11:19:17 1656

原创格子玻尔兹曼法学习记录（附MATLAB画图源程序）

感谢群里的朋友们提供帮助！还是老样子，有啥问题Feel free to tell us~毕竟群众力量大嘛~格子玻尔兹曼救星QQ群：293267908。流体计算领域中，LBM还是个比较新的思想，最近宝宝正在尝试性的进行研究。首先说一下教材吧。我在某宝买了中文版，杨大勇翻译的，中国工信出版。还可以看英文原版哈，英文教材链接：https://pan.baidu.c...

2019-03-30 16:59:37 9750 21

weixin_37783345的博客