Topics created by 李东岳

有几个人联系我了。东西大体上就是网上的东西凑到一起。丝毫看不到用心所在。让我怀疑这种方式领取邀请码是不是浪费我的审核时间。看来为了利益做工作，就是应付。只有真心想写出教程方便大家，才能写出好东西。

如果给我发邮件，没回复。就表示你写的东西我没感到热度。

find src applications -name "*.L" -type f | xargs sed -i -e 's=(YY_FLEX_SUBMINOR_VERSION)=YY_FLEX_MINOR_VERSION < 6 && \1='

ipmitool -I lan -U ADMIN -H 192.168.1.18 sensor thresh FANB lower 100 150 200

各位简直就是大佬中的大佬 厉害厉害 :chitang:

https://askubuntu.com/questions/1143840/downgrade-openmpi-v2-1-1-to-v2-0-2

https://github.com/open-mpi/ompi/issues/6300

有空测试下

可以做如下改动进行输出，至少适用于OpenFOAM-8

打开Smagorinsky.H，添加k_: // Protected data dimensionedScalar Ck_; volScalarField k_;//这是添加的一行 // Protected Member Functions 打开Smagorinsky.C，添加下面代码: : LESeddyViscosity<BasicMomentumTransportModel> ( type, alpha, rho, U, alphaRhoPhi, phi, transport ), Ck_ ( dimensioned<scalar>::lookupOrAddToDict ( "Ck", this->coeffDict_, 0.094 ) ),//这里加个逗号 //下面是添加的 k_ ( IOobject ( IOobject::groupName("k", this->alphaRhoPhi_.group()), this->runTime_.timeName(), this->mesh_, IOobject::MUST_READ, IOobject::AUTO_WRITE ), this->mesh_ ) { if (type == typeName) { this->printCoeffs(type); } }

correctnut()函数通过下面代码替换：

void Smagorinsky<BasicMomentumTransportModel>::correctNut() { k_ = (this->k(fvc::grad(this->U_))); this->nut_ = Ck_*this->delta()*sqrt(k_); this->nut_.correctBoundaryConditions(); fv::options::New(this->mesh_).correct(this->nut_); }

原来这样。

最近在算好多LES算例，直接把服务器当主机用了。日常用没问题。但是发现最大的问题就是word文件的编辑。因为我这面很多都是别人发来的word文件我在上面继续编辑。但是linux下面用各种word编辑软件都会出现格式问题。有些文档就两页，但是导出pdf的时候格式就不对。有些是商务文件，格式对不齐，感觉不是很专业。

我先试试wps，看看会不会好一些。

@李东岳 这个不是，这不是在这边下载BT什么的很危险嘛，搞了个罗马尼亚的抗版权vps 搭了一个aria2用来下载，一个季度8刀，每个月10T流量

说真的，大佬，咱能不能别3分钟限制了。。。。

捕获.JPG

最近一个月浏览量1076980

@cresendo 太复杂了 1 2个小时写不完

这个潘工的联系方式给我呀，500块买不了上当啊

感谢各位大佬！:chitang:

回档了 这是什么意思..

https://www.wikiwand.com/en/Centroid

While in geometry the word barycenter is a synonym for centroid, in astrophysics and astronomy, the barycenter is the center of mass of two or more bodies that orbit each other. In physics, the center of mass is the arithmetic mean of all points weighted by the local density or specific weight. If a physical object has uniform density, its center of mass is the same as the centroid of its shape.

ISB

https://www.isb.bj.edu.cn/admissions/life-in-beijing

不接受中国学生

BSB

https://www.nordangliaeducation.com/en/our-schools/beijing/shunyi/admissions/how-to-apply

暂时不清楚是否接受中国学生

德威

https://beijing.dulwich.org/zh/admissions

不接受中国学生

哈罗

https://www.harrowbeijing.cn/academic-cn/lower-school

貌似接受中国学生

乐成

https://bcis.openapply.cn/

接收中国学生

青苗

https://www.bibs.com.cn/chinese-home/admissions/faq1

招收中国学生

鼎石

https://www.keystoneacademy.cn/

招收中国学生

海嘉

http://www.bibachina.org/index.html?deptId=103

招收中国学生

世青

https://www.ibwya.net/

招收中国学生

其实一楼这个问题我还是有点没想明白。就是一个x方向的一维算例，出现了y方向的速度，应该如何理解..

生活很有规律呀！我就不行了，生物钟紊乱，感觉每天都很疲惫，特“困”生就是我！:haqi:

@李东岳 主要是磁共振成像设计和基于磁共振图像的心脑血管力学分析

是单位化的，写起来忒麻烦了没写一会我改一下。

我试了一下
\begin{equation}
P=
\left(
\begin{matrix}
1&2\\
2&3
\end{matrix}
\right)
\end{equation}
这个就可以:136:

本安装测试的系统为Ubuntu-2004。安装过程需要联网。

在这里下载OpenFOAM-v1912 [版本不能变]，参考 http://dyfluid.com/docs/install.html 进行OpenFOAM-v1912的安装

在这里下载waves2foam [版本不能变]。

创立下图的文件夹路径，注意两点：1）dyfluid-v1912前面的dyfluid是你的linux的用户名，可能是别的。2）把你第二部下载的waves2foam下的所有的文件放在下图中

捕获.JPG

安装必要的库：sudo apt-get install libgsl0-dev

终端切换到上图中的waves2Foam下，./Allwmake编译即可

@李东岳 李老师，可能是您python版本的问题，我用的是python3.8，您应该也有python3但是系统调用的是python2.7，可以试试python3 plot.py运行。
坐等大佬改我的代码:chouchou:

在思考一个审稿人的意见，太难了

漂亮~:xinxin:

测试排版用

引言

湍流运动是目前计算流体力学中困难最多因此也最活跃的领域之一。当湍流存在，则住在其他相关的流动现象，并引致能量耗散、混合以及传热。没有三维的涡，则没有真正的湍流，因为只有在三维的流动中，涡旋才能进行伸展并产生新的涡旋。目前可采用的数值计算方法分为三类：直接模拟（Direct Numerical Simulation，DNS）、大涡模拟（Large Eddy Simulation，LES）和雷诺时均法（Reynolds-averaged Navier–Stokes，RANS）。RANS经过长期的发展，已经非常成熟。但RANS通过将速度进行平均后，并不能捕获湍流中的小涡结构。同时，这些小涡基本是各项同性的。另一方面，从主流中抽取能量的大涡却是各向异性，并且其和计算域的几何、边界、体积力高度关联。在使用RANS的时候，整个流场中必须使用同一个湍流模型对各种尺度下的湍流进行解析，但通常大涡和小涡的表现是不同的。因此研究学者对一种更完善的模型进行了探索。

不同于RANS，LES对大涡进行解析的同时对小涡进行模化。LES认为大涡直接受边界条件的影响因此对其解析，但小涡是各项同性的因此他们表现相同，可以进行模化。由于LES把小涡进行了模化，因此最小的网格单元需要大于Kolmogorov尺度（最小的涡旋尺度）。同时LES的时间步可以比DNS大的多。因此，对于给定的计算资源，相对于DNS，LES可以计算更大雷诺数的算例。另外，不同于RANS中平均的概念，LES使用的是一种空间滤波技术。LES模型的概念如下：

首先要确定一种滤波函数和截止尺度$\Delta$。这样，就可以对所有大于截止尺度的涡进行非稳态计算；

使用滤波函数对依时变量进行空间滤波操作，在这一步，小于截止尺度的涡被过滤掉；

在解析大涡和模化小涡的数学操作中，会产生一个亚格子尺度应力项（Sub-grid-scale Stress，SGS），亚格子尺度应力需要通过SGS模型来模化；

在LES中，截止尺度是用来表明“多大的涡才算大涡”的概念。其可以为任意大小，但是选择比网格还要小的截止尺度是没有意义的。在笛卡尔网格下，最简单的截止尺度这样计算：
\begin{equation}
\Delta=\sqrt[3]{\delta x \delta y \delta z}
\end{equation}
其中$\delta x$等为笛卡尔网格下网格单元的边长。其他不同的截止尺度计算方法还有最大边长法、普朗特混合长法等。

方程与模型 滤波N-S方程

笛卡尔坐标下的连续性方程为：
\begin{equation}
\frac{\partial \rho}{\partial t}+\nabla \cdot(\rho \mathbf{U})=0
\end{equation}
对$\nabla \cdot(\rho \mathbf{U})$进行滤波后有滤波连续性方程：
\begin{equation}
\frac{\partial \rho}{\partial t}+\nabla \cdot(\rho \mathbf{\bar{U}})=0
\end{equation}
其中$\mathbf{\bar{U}}$为滤波后的速度。且有：
\begin{equation}
\mathbf{U}=\mathbf{\bar{U}}+\mathbf{U'}
\end{equation}
其中$\mathbf{U'}$为残余速度，且$\mathbf{U'} \neq 0$。下图表示的是一个对$x$分量速度进行高斯滤波之后的滤波速度分量和残余速度分量的示意图。
U.jpg
另外，有不可压缩动量方程为：
\begin{equation}
\frac{\partial \mathbf{U}}{\partial t}+\nabla \cdot (\mathbf{U} \mathbf{U})=-\nabla \frac{p}{\rho}+\nabla \cdot(\nu \nabla \mathbf{U})
\label{mom}
\end{equation}
同理有：
\begin{equation}
\frac{\partial \mathbf{\bar{U}}}{\partial t}+\nabla \cdot (\overline{\mathbf{U} \mathbf{U}})=-\nabla \frac{\bar{p}}{\rho}+\nabla \cdot(\nu \nabla \mathbf{\bar{U}})
\label{momF}
\end{equation}
在方程\eqref{momF}中，除了待求的$\bar{\mathbf{U}}$和$\bar{p}$外增加了一个未知量$\overline{\mathbf{U} \mathbf{U}}$。为了将问题简化，把方程\eqref{momF}的第二项进行变化：
\begin{equation}
\nabla \cdot (\overline{\mathbf{U} \mathbf{U}})=\nabla \cdot (\mathbf{\bar{U}} \mathbf{\bar{U}})+\left(\nabla \cdot (\overline{\mathbf{U} \mathbf{U}})-\nabla \cdot (\mathbf{\bar{U}} \mathbf{\bar{U}})
\right)
\label{W}
\end{equation}
将方程\eqref{W}带入到方程\eqref{momF}中有：
\begin{equation}
\frac{\partial \mathbf{\bar{U}}}{\partial t}+\nabla \cdot (\mathbf{\bar{U}} \mathbf{\bar{U}})=-\nabla \frac{\bar{p}}{\rho}+\nabla \cdot(\nu \nabla \mathbf{\bar{U}})
-\left(\nabla \cdot (\overline{\mathbf{U} \mathbf{U}})-\nabla \cdot (\mathbf{\bar{U}} \mathbf{\bar{U}})
\right)
\label{momFF}
\end{equation}
对比最初的的N-S方程\eqref{mom}，方程\eqref{momFF}中的最后一项$-\left(\nabla \cdot (\overline{\mathbf{U} \mathbf{U}})-\nabla \cdot (\mathbf{\bar{U}} \mathbf{\bar{U}})
\right)$为滤波操作产生的特殊项。对其展开有：
\begin{equation}
\nabla \cdot (\mathbf{\bar{U}} \mathbf{\bar{U}})=\nabla \cdot \left[\begin{matrix}
\bar{u}_1\
\bar{u}_2\
\bar{u}_3
\end{matrix}\right][\bar{u}_1, \bar{u}_2, \bar{u}_3]=\nabla \cdot \left[
\begin{matrix}
\bar{u}_1 \bar{u}_1 & \bar{u}_1 \bar{u}_2 & \bar{u}_1 \bar{u}_3\
\bar{u}_2 \bar{u}_1 & \bar{u}_2 \bar{u}_2 & \bar{u}_2 \bar{u}_3\
\bar{u}_3 \bar{u}_1 & \bar{u}_3 \bar{u}_2 & \bar{u}_3 \bar{u}_3
\end{matrix}
\right]
\label{T1}
\end{equation}

\begin{equation}
\nabla \cdot (\overline{\mathbf{U} \mathbf{U}})=\nabla \cdot \overline{\left[\begin{matrix}
u_1\
u_2\
u_3
\end{matrix}\right][u_1, u_2, u_3]}=\nabla \cdot \left[
\begin{matrix}
\overline{u_1 u_1} & \overline{u_1 u_2} & \overline{u_1 u_3}\
\overline{u_2 u_1} & \overline{u_2 u_2} & \overline{u_2 u_3}\
\overline{u_3 u_1} & \overline{u_3 u_2} & \overline{u_3 u_3}
\end{matrix}
\right]
\label{T2}
\end{equation}

@wangfei9088 哈哈 多谢支持多谢支持 :chigua:

学校有大礼堂可以放电影，春节档电影均已刷完~

在breaking wave领域，这种海浪结构，会将空气卷吸并形成很多小气泡。最近提交一个稿件，审稿人推荐我们讨论一下这方面的研究。我详细看了一下，确实有点意思。

看下图，一些实验研究已经证明，这种海浪拍打过程卷吸进的气泡具有尺寸分布。见过很多文献，都是用的DNS直接模拟研究的，结果也都能证实。

但目前我想的是，能否用PBE模型，或者其他类似模型做相关的研究。

捕获.JPG

气泡直径和分布的关系
捕获.JPG

一些参考文献

Deike, L., Melville, W. K. & Popinet, S. 2016. Air entrainment and bubble statistics in breaking waves. J. Fluid Mech. 801, 91–129.
Wang, Z., Yang, J., Stern, F., 2016. High-fidelity simulations of bubble, droplet and spray formation in breaking waves. J. Fluid Mech. 792, 307–327.
G. Soligo, A. Roccon, and A. Soldati, Breakage, coalescence and size distribution of surfactant-laden droplets in turbulent flow, J. Fluid Mech. 881, 244 (2019).
Ahmed, Z., Izbassarov, D., Costa, P., Muradoglu, M., Tammisola, O. 2020. Turbulent bubbly channel flows: Effects of soluble surfactant and viscoelasticity, Computers & Fluids 212, 104717.

@李东岳 用的是真系统

也是实测。OpenFOAM自带算例。256万网格。虚拟机16核跑了1772s，真系统跑了548s。

虚拟机大约要慢3倍

@西湖冷月 这就是天秀

没币发帖

感觉这个问题是计算流体力学本身不同方向的割裂造成的，流体力学可以不十分懂算法和底层的方程原理也能把流动算出来，做算法的可以不用管具体流动应用场景，直管把方程解出来。但其实这两者是有桥梁的，粘性流体力学基础之类的相对底层一点的东西才是涵盖算法和工程的，算是个纲领。

@李东岳 wifi用5G就不会干扰了啊

TEMOM只展开$(v+v_1)^{1/2}$，任意阶矩展开任意阶都是泰勒的截断误差，是个原式的高阶无穷小，原式也就是$1/2$次方。
但DEMM把$(v+v_1)^k-v^k-v_1^k$也带上展开了，这个截断误差是变的，k越大，误差越大。越高阶矩，误差的量级越大，虽然是都是无穷小。

所以我觉得越高阶这个方法的优势越不明显。当然都是我猜的:135:
在提出DEMM的那篇文章里的数据也有这样的苗头，低阶矩符合很好，高阶的出点问题又修复了。
我把相关资料发给老师。

:143:

好像有个wsl和wsl2，我也没太研究，我直接关了..

捕获.PNG

看方程20，看起来植入还挺容易的，类似这种模型在OpenFOAM中植入做个验证研究都可以发

李老师您好，
我按照您的提示：

sudo update-alternatives --install /usr/bin/gcc gcc /usr/bin/gcc-5 100

目前默认gcc-5

gcc (Ubuntu 5.3.1-14ubuntu2) 5.3.1 20160413 Copyright (C) 2015 Free Software Foundation, Inc. This is free software; see the source for copying conditions. There is NO warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.

我想切换回g9，用命令

sudo update-alternatives --config gcc

得到：

There is only one alternative in link group gcc (providing /usr/bin/gcc): /usr/bin/gcc-5 Nothing to configure.

没有切换提示，咋办？

我用命令

ls /usr/bin/gcc* /usr/bin/gcc /usr/bin/gcc-ar /usr/bin/gcc-nm /usr/bin/gcc-ranlib /usr/bin/gcc-5 /usr/bin/gcc-ar-5 /usr/bin/gcc-nm-5 /usr/bin/gcc-ranlib-5 /usr/bin/gcc-9 /usr/bin/gcc-ar-9 /usr/bin/gcc-nm-9 /usr/bin/gcc-ranlib-9

应该是g9也在的，不知道为啥不能切换，跟百度不一样

点赞关注

ugks, dugks和LBM都是解得动力学方程，但是算法还是有极大区别的

虽然我知道我弄不成（硬盘没有没分区的剩余空间了），但我还是想试试，然后根目录炸了:jingya:
我说我的根目录怎么一点空间都不剩了，还报警。生成了一个 swapfile 占满了。操作有风险，尝试需谨慎。

这是占根目录的空间当作swap用，绕过了重新分区，挺好，就是空间没的有点……突然 _(┐「ε:)_

为什么要配台丽台的卡？

595269271@qq.comqq.com

https://packages.vmware.com/tools/esx/latest/windows/

@东岳
不不不，您才是大佬，我是仰慕小菜鸡:xiexie:

@sibo 并行限制2核，网格限制是52万

http://dyfluid.com/rhocentralfoam.html 参考这个 把Lax-Friedrichs方法植入进去了 耗散的很 不过是个OKS好材料

捕获.JPG

@walten 在 如何CFD自学成才 中说：

@闻久STU
乐观点， 六年:high: 这样才学得深入

不能提:135: 长期抗战是肯定的了