- 博客(31)
- 收藏
- 关注
RSS订阅
原创 OpenFOAM 使用技巧总结
OpenFOAM使用技巧总结（持续更新2018.12.10）矢量数组可以新建为矢量列表形式，初始化方法为：List<vector> a = (n,vector::zero)生成不同数据类型的链表并初始化：#include"Tuple2.H" List<Tuple
2018-12-10 11:07:56
8740
原创 solidwork导出Ascii格式的STL文件
1. 打开*.STEP文件如图点击如下图标选择STEP格式文件并打开2. 选择ASCii格式如上图所示，默认的输出格式为二进制，若需要Ascii格式，需要点击步骤1中设置，然后搜索Ascii，选择Ascii，也可修改单位3. 选择stl格式另存为即可...
2021-02-19 22:39:27
3755
3
原创 OpenFOAM中进程之间的信息交流
应用场景：在OpenFOAM中利用Gussian插值等插值方法将欧拉网格中流场信息插值到某一点过程中需要寻找该点周围网格信息以及流场信息，如果该点落在进程交接面处那么就需要进程进行信息传递，让计算进程可以看到相邻进程的网格信息以及流场信息。解决方案：如上图所示计算区域被分为4个进程，主要解决思路为由于其他进程的流场信息以及网格拓扑信息在当前进程是看不到的，所以需要一个二维数组来保存其他进程发送的信息。以processor0二维数组为例，如下表所示processor_Idcells
2020-12-10 16:10:03
1208
1
原创 vscode+github进行版本控制
参考自该链接1. 本地初始化gitgit init ：创建一个.git目录git status -sb可以查看git下文件的状态git add 文件路径 :添加文件到暂存区/git add .git commit -m “版本1”:添加该版本说明2.上传到github生成RAS keygit remote add orgin —（---为SSH内容）git push -u origin master出错 ! [rejected] master -> mast
2020-12-10 16:08:42
254
原创 OverLeaf编辑latex投稿Physics of FLuids参考文献导入
问题描述：OverLeaf编辑latex投稿Physics of FLuids时不能上传Reference.bib文件，所以不能在线编译出带参考文献的PDF解决方案：如下图代码块展示的样子，首先在overLeaf中带着Reference.bib完整的编译并成功生成PDF全文。\section*{References}\bibliography{Reference}% Produces the bibliography via BibTeX.如图所示点击`logs and outputs f
2020-10-20 15:10:05
1240
原创 latex以及perl安装
latex以及perl安装前言安装前言latex安装版本为（CTeX_2.9.2.164_Full.exe）,ActivePerl主要是在生成审阅模式下使用其安装版本最好为:ActivePerl_5.16.2.3010812913安装首先安装latex看Micktex里latexdiff安装包安装没有具体做法：在开始菜单—CTeX—Package Manager (Admin)，然后在name中搜索latexdiff,如有只有减号没有加号，说明已经安装上了。安装Aciveperl
2020-10-20 11:29:03
950
原创 latex审阅模式出现Missing \begin{document}错误
项目场景：windows 10中使用latexdiff生成审阅版本tex文件问题描述：latex生成审阅版本的.tex文件编译过程中出现Missing \begin{document}，出现该错误并不一定说明你的.tex文件缺少\begin{document}，而是说明你的\begin{document}之前出现了一些编辑器（我用的是WinEdt）不认识的字符@OverrideThis is pdfTeX, Version 3.1415926-2.3-1.40.12 (MiKTeX 2.9) (p
2020-10-20 11:18:48
19378
原创 paraview生成颗粒在床内运动的迹线
如图所示，采用TemporalParticlesToPathlines Feilter可以获取颗粒的path lines 并且对path lines使用tube Feilter可以改变迹线的宽度。
2020-10-15 14:57:56
1010
原创 使用Box–Muller transform随机生成服从Gaussian以及log-normal分布的点
随机生成服从Gaussian分布点的方法有很多种，具体可以参考链接；本文采用的是Box-Muller transform方法。#include<stdio.h>#include<stdlib.h>#include<math.h>#define M_PI 3.14159265358979323846double generateGaussianNoise(double mu, double sigma){ double twopi = 2.0 * M_PI;
2020-10-14 19:22:18
750
原创 计算概率密度分布函数后处理程序
最近需要批量计算部分数据的概率密度分布函数，肝了两天总算搞定，先整理如下：具体思路分为3步：从原始数据中爬取数据计算概率密度分布画图画图过程中由于横纵坐标都是对数坐标系下完成，所以首先需要将数据转换到对数域，然后划分区间，最后再求指数转换到正常域；绘制曲线时选择对数坐标系即可。具体实现细节请点击链接...
2020-09-15 17:29:35
313
原创 为什么C++交换函数传入数据无法完成交换
#include <iostream>using namespace std;void swap(int a, int b) { int x = a; int y = b; a = y; b = x; cout<<"swap "<<&a<<" "<<&b<<endl;}int main(){ int a = 2; int b = 8; swap
2020-07-23 19:53:06
329
原创 APPLICATION TERMINATED WITH THE EXIT STRING: Floating point exception (signal 8)
其实这个问题发现好久了，做一下笔记：自己在超算上运行OpenFOAM中solver时候会碰见这种错误：APPLICATION TERMINATED WITH THE EXIT STRING: Floating point exception (signal 8)发现解决办法就是自己重新写一个Makefile编译一下；具体做法为：将Make/options中项添加到Makefile的CINC...
2020-04-02 17:17:55
3298
原创 LaTex参考文献由author-year变为numbers
在撰写科研文章时不同期刊对文章引用格式不同。最近需要将authoryear模式的参考文献转换为numbers（即,gas-solid fluidization (Grace, 1996; Zhao et al., 2015; Hou et al., 2016)转换为gas-solid fluidization [14–16]）我需要修改的地方有：将\documentclass[authorye...
2020-04-01 14:35:48
3500
原创 Sutherland's law动力学粘度关联式
Sutherland’s law是关于动力学粘度 μ\muμ 的经验关联式：μ(T)=T3/2(T0+ST+S)\mu(T)=T^{3 / 2}\left(\frac{T_{0}+S}{T+S}\right)μ(T)=T3/2(T+ST0+S)其中T0=273.15KT_0=273.15KT0=273.15K,S=110.4KS=110.4KS=110.4K为经验参数。1更详细的内...
2020-03-17 15:55:30
5844
原创 reactingParcleFileFoam解析
离散粒子相求解：parcels.evolve()parcels的原型为creatClouds.H中的basicReactingCloud parcels( "reactingCloud1", rho, U, g, slgThermo);basicReactingCloud ,该原型定义位于：/opt/OpenFOAM/OpenFOAM-2...
2020-03-17 12:32:29
1298
原创 OpenFOAM中的particle,parcel和cloud
/src/lagrangian/中包括/solidParticle/库以及/intermediate/其中solidParticle更新颗粒位置与与速度但是流固之间为one way couple即只包括气固之间的相互影响？ 。而intermediate则包括更多的不同复杂度的parcel和cloud库。当然/spray/中也包含parcel和cloud库，与/intermediate/中不同的时，...
2020-03-13 16:54:58
2218
原创 python读取txt数据处理后.csv格式输出
本代码读取了一个两列的txt文本文件，对其中数据进行处理画图后，并用csv格式输出：#!/usr/bin/pythonimport numpy as npimport matplotlib.pyplot as plt import csv#plt.rc('font',family='Times New Roman')x=[]y=[]xcorr=[]ycorr=[]with...
2020-03-07 14:02:01
605
原创 一个openfoam后处理APP（offline）
本文介绍一个openfoam的后处理APP，主要作用是读取并输出轴向颗粒浓度，输出到文件，并用python直接绘图，主要是对数组的一些处理，该链接附带详细代码及算例，具体如下：后处理代码/*---------------------------------------------------------------------------*\ ========= ...
2020-01-13 10:35:09
1060
原创 利用python中plt绘图技巧（持续更新）
plt.rc('text',usetex=True) #使用latexplt.rc('font',family='serif')#修改字体plt.xlabel(r"$\theta^ \circ$",size = 14)#latex中含有希腊字母plt.savefig('post_sum.png', format='png', dpi = 300)#修改图表分辨率...
2020-01-08 19:38:40
1661
原创 将stl文件中vertex顶点由顺时针（clockwise）排列改为逆时针（counterclockwise）
最近需要读取solidWork软件生成的复杂几何体（stl格式），用于识别一系列点位于复杂几何体内部还是外部，发现识别函数需要stl文件符合右手规则，即点的排布要逆时针顺序，而且法线方向朝外，但是solidWork生成的文件不符合这个特征，所以写了一个c++的小程序，用来实现该过程，具体内容如下：1. stl文件格式： facet normal -0.996195 0.0871562 ...
2020-01-08 19:35:53
852
原创 OpenFOAM中使用subsetMesh挖走一个box
OpenFOAM中使用subsetMesh挖走一个box利用`subsetMesh`程序挖走一个box，具体做法是：利用subsetMesh程序挖走一个box，具体做法是：利用topoSet选中要挖出的box c0并删除-(system/topoSetDict)actions( { name c0; type cellSet; ...
2020-01-03 20:34:10
1956
1
原创 坐标点转换为paraview形式
**应用场景：**将含有xyz坐标的文件转换为可以在paraview中可视化的文件**缺点：**需要制定点的个数**格式：**如下为wallParticleThermo.out文件的分别为点的xyz坐标0.00577497 0.0339325 0.0437007代码如下：int nParticles = 16514;float x = 0.0;float y = 0.0;floa...
2019-09-07 17:10:28
997
2
原创 openfoam中获得不同方向上的网格长度
获取不同方向上的网格长度cell size(x,y,z)const faceList & ff = mesh.faces();const pointField & pp = mesh.points();forAll ( mesh.C(), celli){ const cell & cc = mesh.cells()[celli]; labelLis...
2019-09-06 09:14:20
1567
原创 openfoam寻找与某一边界相邻的一层网格
基本思路定位到边界的faceID在owner中根据faceID选择网格cellvolScalarField post ( IOobject ( "post", runTime.timeName(), mesh, IOobject::NO_READ, ...
2019-09-06 09:12:12
1134
原创 虚拟机上在线安装OpenFOAMv6网络连接问题
虚拟机上网络配置在虚拟机虚拟机上安装CentOs7系统教程很多，在此不多做赘述，详情可参考.不过需要注意的是，如果您想最小安装（无图像化界面，通过命令来操控）则选择下图的软件选择中最小安装即可如果要图像化界面，则需要点击软件选择，进去选择GNOME界面：接下来安装教程完成所有的安装以及网络设置。可以通过在终端键入ping baidu.com来检查网络是否配置成功（如下图即为成功），...
2019-06-13 10:07:15
1066
原创 OpenFOAM按进程输出和输入文件
在某些前处理过程中需要按进程输入或者输出某些文件，例如输出网格中心高度大于0.2m的网格中心点： const volVectorField& C = mesh_.C(); string nProc = std::to_string(Pstream::myProcNo());//不同编译器形式不同 fileName gamma = "gamma" + "-" + n...
2019-03-19 16:52:04
1966
原创 openfoam中Eular网格上的流场值插值到某一个点
openfoam中Eular网格上的流场值插值到某一个点在模拟气固两相鼓泡床过程中需要将场量（流体速度，压力梯度）插值到相应颗粒位置，openfoam中自带了这种插值函数源文件可点击总共有六种插值方法cell, cellPoint, cellPointFace, cellPointWallModified, cellPatchConstrained 和pointMVC详细的插值方法介绍。以cel...
2018-11-15 18:50:52
1814
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人