Algorithm

• 我

  利用交错网格simple算法计算顶盖方腔流发散
  • 我是河滩  

  2
  2
  Posts
  888
  Views

  

  如果能实现顶盖驱动MAC算法，simple也没有问题啊~可以先试试MAC，速度边界条件都是第一类边界条件，其余都为0.只有顶盖有速度。
• 

  CPU在并行时候的指派
  • 李东岳  

  5
  5
  Posts
  2114
  Views

  气

  @程迪 请问linux如何并行计算，在ubuntu下没有问题，但是在linux下并行，显示这样的错误，希望您指教 OpenFOAM-6(23) mpirun -np 2 interFoam parallel -------------------------------------------------------------------------- The value of the MCA parameter "plm_rsh_agent" was set to a path that could not be found: plm_rsh_agent: ssh : rsh Please either unset the parameter, or check that the path is correct
• 

  有限面方法 Finite area method
  • 李东岳  

  2
  2
  Posts
  727
  Views

  W

  以前简单看过，好像是研究表面活性剂的方法。大概是在三维曲面上用类似FVM的方法离散控制方程并求解。
• 

  实验室规模鼓泡床直接模拟
  • 李东岳  

  8
  8
  Posts
  1590
  Views

  星

  Examinination of bubble collisions and coalescence in bubbly flows D. Broder and M. Sommerfeld 还有这一篇，这好像都是小的那个做的，大的具体尺寸我也不知道，不过什么时候能在做实验不好说，实验仪器快被转移了，大的那个。。
• 

  25亿网格直接模拟
  • 李东岳  

  1
  1
  Posts
  1074
  Views

  

  Direct numerical simulation of reactor two-phase flows enabled by high performance computing
• 

  开放合作系列
  • 李东岳  

  3
  3
  Posts
  1584
  Views

  

  Simulations of gas-liquid dispersions by four-way coupled E-QBMM with aggregation and breakage processes Abstract Results Main References Job: Data processing Test cases debugging and running Writing paper in English if you are willing to Requirements: Benefits: You can get around 30000 yuan bonus from your host (if you write the manuscript). You can use this paper to get promoted faster (if you write the manuscript). Another 10000 yuan from my side If I write the paper, I can support you 40000 yuan as well. However, I dont have much time to write. So I prefer you writing the paper and to be the first author. Project duration
• 

  TDMA算法
  • 李东岳  

  3
  3
  Posts
  1364
  Views

  Y

  在isoadevtor的界面重构里面见过，拟合界面等值面和网格体积分数关系的三次函数吧，适合求解小矩阵线性方程
• 

  滑移网格算法相关文献求助
  • 肖艺  

  3
  3
  Posts
  1096
  Views

  

  @东岳 谢谢！
• K

  颗粒表面化学反应的文献
  • kechunhai  

  3
  3
  Posts
  929
  Views

  K

  @ibelief 在 颗粒表面化学反应的文献 中说： 基于简单碰撞理论煤粉燃烧动力学模型的研究 非常感谢
• S

  关于群体平衡方程中sectional method (Gelbard & Seinfeld)的守恒性
  • Shang XP  

  7
  7
  Posts
  1826
  Views

  

  moment method对于nucleation和condensation的模拟简单吗？ 即为简单。但是溶解/蒸发数值上要复杂些。2017年我在HZDR的时候同事用fixed pivot method模拟增长，我用moment method，相比来看moment method要更简单，并且不需要处理相空间的drift：https://www.sciencedirect.com/science/article/pii/S1004954118301174/pdfft?md5=84216a8842ce3fd6d1a99d4f327c1ee1&pid=1-s2.0-S1004954118301174-main.pdf
• Q

  带摩擦的运输方程的矩阵的特征值的虚数和实数部分作图, 为什么呈现一个椭圆?
  • qingdong_wang  

  1
  1
  Posts
  500
  Views

  Q

  带摩擦的运输方程是 u_t = lambda * u_x = eps * u_xx, 通过变化, 我们得到了, ,其中D是导数矩阵, D2 是二阶导数矩阵 然后, 我们对M进行求特征值(EW), 得到了这种图片, 我想知道这是为什么呢?
• L

  到底静压值是怎么得来的？
  • l.j刘侃  

  1
  1
  Posts
  719
  Views

  L

  压差是分析很多工程设备的关键参数，在之前，主要都是针对多孔介质的Ergun模式来处理这一块。 也跟试验对比很好。但是，最近新开发的一个模型，虽然与试验的结果对比比较好，但感觉总是有问题。 于是不知道CFD中静压值到底是怎么得出来的？ 对于一块多孔介质好像比较好处理，但是如果这个时候多孔介质表面沉积了灰，这个时候这块的整个压降会升高，然后灰越堆越多整个表面的压降都会升高。然而，这个局部压降和表面灰的分布十分相关，造成整个表面的这个灰分的压降相差一千多Pa。而最后一般的检测压降点，都是压差变松器的两个接口，一般都是设置在容器外壁的。但是如果我模拟检测点和试验一样的时候，这个时候的压差比我的表面计算压差平均值要小得多。（压差是1200多Pa，但是我的多孔介质有800多Pa，加上平均粉尘层的压降有640多Pa，合起来应该有1400多Pa，而不是1200多Pa.） 这就产生了一个问题，这个通过N-S方程求解的压力项到底是怎么处理的? 当他通过一层多孔介质的时候，多孔介质内部有压力源项，而在外部没有，那么当他气体从多孔介质外向内流通的时候，他的这层压力和速度是怎么计算的？如果这个时候我再添加一个源项在多孔介质层，而这个源项在多孔介质外是没有的，那么这层额外的源项是不是也会通过增加流体的静压来表达呢？（通过伯努利方程似乎很好理解，但是N-S是怎么实现的，他的这个\Delta P 是怎么获得的？）
• 

  GPU计算: 多相流应用
  • 李东岳  

  3
  3
  Posts
  1552
  Views

  L

  DEM加速比很高 CFD好像一般
• 

  如何理解CFD问题只有三维
  • 凤雏岂能消逝  

  3
  3
  Posts
  1327
  Views

  L

  二维的vorticity只有一个方向 实际上从涡的形成来看 二维似乎并不适合一些涡量的模拟 但是对于一些简单的流场 实际上2D具备非常大的优势 他的处理简单 很多算法很好自编程实现 并且网格可以完全mapped 最大的优势 还是在于后处理的一目了然 这就是为什么大多数理论研究还是集中2D 但是2D可能不是太适用Fluent Fluent主要接轨工程实际问题
• L

  关于多孔介质近壁面层的流体速度拟合问题
  多孔介质；近壁面流动；rans • • l.j刘侃  

  1
  1
  Posts
  635
  Views

  L

  各位专家好： 在多孔介质模拟的过程中，我发现了一个问题，虽然在发表文章中省略了。但是确实是一个直观的问题，我的模型，RANS方法中，不管采用哪一种，在多孔介质的近壁面处的拟合的速度与实测值总是具备一定的差距。这种差距我分析认为是在较高阻力降的多孔介质环境下，气流在流经多孔介质表面时，一部分流体会因为高阻降而产生流向变化，形成比较复杂的近壁面流场。而在这个近壁面流场中，因为Fluent默认是采用的多孔介质域或者是多孔跳跃面，应该是对这层近壁面没有处理的，依然按照充分发展层处理，造成这层近壁面的流速（甚至是流向）与实测值都有较大区别。 虽然能够较好的拟合实时压降（Endo based equations），但是如果在近壁面流场不符合实际的话，那么很显然对于最后的表面覆尘的预算难以符合实际情况。不知道大家对于这种多孔介质的近壁面流是否有过研究，或者能够有比较好的解决思路？欢迎大家提出宝贵的想法和建议，感谢，再次感谢！ 下图为提出新模型对于滤管外壁的模拟和实测结果对比：
• 

  欧拉双流体模型中应力张量的表示问题
  • Alex  

  3
  3
  Posts
  1082
  Views

  

  @东岳 老师，没理解您的意思。后面的$\alpha_q\lambda_q\nabla·\vec{v_q}I$中单位二阶张量是表示体形变率仅出现在法相应力中，红框内等于$-\frac{2}{3}\lambda_q$时才能与法相应力对应啊
• 

  关于“CFD中的能量方程”中最后的方程（16）
  • hurricane007  

  10
  10
  Posts
  4919
  Views

  

  fvc::laplacian(M,phi)和(fvc::grad(M)&fvc::grad(phi)) + M*fvc::laplacian(phi)是等价的吗？ 数学上是等价的。数值上后者第一项可能会出现震荡问题，尤其在网格比较糙的情况下。
• 

  二维轴对称动量方程
  • Alex  

  3
  3
  Posts
  937
  Views

  

  谢谢老师，我找出来看一下
• J

  目前计算可压缩流体中，NS方程和能量方程到底如何耦合求解？
  • Juanyong  

  4
  4
  Posts
  1260
  Views

  

  https://www.ansys.com/-/media/ansys/corporate/resourcelibrary/article/aa-v2-i2-coupling-momentum-and-continuity.pdf 刚看了下，Fluent那面也只是速度和压力耦合求解。能量抛在了外面。 但是另一个文档速度压力能量放在了一起 ： http://www.bakker.org/dartmouth06/engs150/05-solv.pdf 速度压力耦合求解在高超那面用的多，如果分析N-S方程的数学特征（双曲抛物），你需要把u, v, w, p当成矢向量，然后求特征值特征向量。接下来求解这个矩阵就是耦合算法。我目前只是了解到思路，没亲自植入过。可
• G

  dyfluid.com上“icoFoam解析”中“方程(3)中的...” 应该是方程（6）吧
  • ghoulyy  

  2
  2
  Posts
  925
  Views

  

  多谢，已更新
• G

  dyfluid.com上“CFD中的壁面函数“”公式（3）（5）（6）有出入
  • ghoulyy  

  2
  2
  Posts
  837
  Views

  

  5右边多写了个$\rho$吧，谢谢！回北京我更新一下
• 云

  超声速湍流平板计算问题
  • 云飞  

  3
  3
  Posts
  1050
  Views

  云

  网格加密到Y+<1，计算结果变化不大，对流通量用二阶MUSCL, AUSM格式
• 

  关于bubble flow中的通量计算问题
  • yhdthu  

  5
  5
  Posts
  1415
  Views

  

  @东岳 如此说来，除了加密网格，没有好的解决方案啦
• 迅

  连续性方程是双曲型方程，还是椭圆型方程
  • 迅哥的三哼经  

  2
  2
  Posts
  917
  Views

  

  从特征来看，是双曲形方程。是步进形式的问题。
• 

  一个小bug调试了一个月
  • 李东岳  

  3
  3
  Posts
  1154
  Views

  

  对，然后不计算边界的话就是0
• 偷

  求助贴！关于无稳态流场的CFD计算
  • 偷懒的小熊猫  

  5
  5
  Posts
  1698
  Views

  偷

  @him 谢谢解答，另外一般商业软件的定常求解器是用的时间推进来求解的吗？如果是时间推进的定常求解器，这种不稳定的振荡可否理解为流场随时间的变化？
• 

  LES和RANS中，离散的对象是相同的吗？
  • .J.  

  31
  31
  Posts
  20943
  Views

  H

  如果使用standard KE，计算二维流场，在某个高雷诺数下把Cd, Cl, St都算对，那么说明碰对了，意思就是好运。这并不能说明standard KE可以算对漩涡脱离。
• 程

  Rhie-Chow插值有没有对应的连续的微分方程/积分方程？
  • 程迪  

  21
  21
  Posts
  9736
  Views

  H

  但是RC插值没有离散，因为它是显式计算源项 -- “速度插值+高阶压力导数项”。
• 韬

  自适应网格优化 求助！！!
  • 韬光养晦  

  1
  1
  Posts
  842
  Views

  韬

  基于2011年发表在Journal of Computational Physics的文章Adaptive mesh refinement based on high order finite difference WENO scheme for multi-scale simulations，并根据您文章的网格自适应算法编写了一套计算程序，在计算双马赫反射算例过程中，遇到了一些问题，受困许久，希望大神指点迷津~ 在此诚挚感谢 遇到的问题具体问题如下： 在计算双马赫反射问题时，在Kelvin-Helmholtz rolls两侧流场出现较大的震荡，局部密度云图及网格自适应加密情况如图1所示，震荡区域存在粗细网格过渡,故怀疑在粗细网格过渡中存在算法错误或者代码错误（我已多次检查代码，未发现代码错误，但不能完全排除仍有代码错误存在）。 （a）局部密度云图（红框处为存在流场震荡处） （b）局部网格生成 （c）总体密度云图 （d）总体网格生成（黑色区域表示网格加密） 图1 局部密度云图及网格自适应加密情况 该程序采用C++语言编写，流场计算使用有限差分法，五阶WENO格式，Steger-Warming通量分裂，2D欧拉方程，时间推进采用3阶Runge-Kutta方法；网格自适应加密方法和加密准则与您的论文有所不同，采用基于原始根网格进行加密，每个网格单元生成3×3个子网格单元，加密准则使用密度梯度(∂ρ⁄(∂x+∂ρ⁄∂y))。一共使用三层网格，在计算域[0,3]×[0,1]中，原始根网格单元数量为300×100。Layer 0表示原始根网格，Layer 2表示最密的网格。Layer 2比Layer 1略大，覆盖Layer 1，两层网格的间隙为3个根网格单元的尺度，如图2所示。以此希望精准捕捉激波，并适当减少网格生成数量。子网格层的时间步长为父网格层的1/3。 父网格层每计算完一步，对父网格节守恒量（ρ，ρu，ρv，E）和守恒量对时间的导数[ρ’，（ρu）’，(ρv)’，E’]进行五阶加权(类似WENO加权)空间插值，求解出子网格层边界的虚网格节点的相应守恒量及其导数。并将前一时间步和本时间步插值出的通量及通量导数使用Hermite插值，以获得三个子时间步的相应变量，作为加密网格层的边界条件，并使用式（1）作为Runge-Kutta时间推进的中间步的边界条件。 在子网格层三个时间步计算完成后，父网格层被子网格层覆盖的节点数值被子网格对应同位置节点更新。每一或两个根网格时间步重新生成一次加密网格。 （a）粗网格和细网格 b）不同层次网格 图2 网格加密示意图 其他相关的情况介绍： 在一次测试中使网格只加密，不减少，结果显示流场基本没有前述震荡，结果如图3所示。此算例应该可以证明该程序计算格式是正确的。然而，由于生成网格实在太多，不具实际应用意义。 （a）密度云图 （b）网格生成情况（黑色为加密区域） 图3 网格只加密不减少的计算结果 不知大家在调试程序过程中，是否遇到过类似的情况。 希望大家能够指出 （1）我对算法的理解是否存在错误，特别是粗细网格过渡部分的算法； （2）解决该问题有何建议； （3）能够告知粗细网格过渡区域，造成的流场震荡的原因可能有哪些; 再次感谢大家！
• 一

  OpenFOAM里面有热通量heat flux的边界条件吗？
  • 一二  

  4
  4
  Posts
  1518
  Views

  一

  @东岳 这个边界条件已经整合到exteranlHeatFluxTemperature里面了，算的也是温度梯度，也就是$q=-\kappa \frac{\partial T}{\partial n}$，这种边界按照FDS技术手册里面应该是DNS的计算方法。 LES计算时给定的方法是$q=h \Delta T$,$h$是对流换热系数。openfoam里面的边界条件要么是固定值要么是梯度，这种差值不知道在of里面能不能实现
• C

  关于气液分离的采出分率问题
  • charles1992  

  2
  2
  Posts
  838
  Views

  C

  补充一下，不知道这些图中的这些数据点是如何提取的，谢谢大家
• S

  均匀各项同性湍流的初场设置问题
  • ShiFC  

  4
  4
  Posts
  1719
  Views

  B

  @shifc 不知道你现在解决了没有，我最近做格栅湍流也碰到了同样的问题，我的解决方法是谱空间的复速度对三个坐标轴必须都是共轭对称的。
• P

  Hill vortex是怎么产生的？
  • Phoenix_DXY  

  2
  2
  Posts
  1053
  Views

  S

  你好👋
• 

  数值传热学P94，关于内节点法的离散过程。
  • 羽之下  

  4
  4
  Posts
  1593
  Views

  

  额，自问自答？
• C

  计算对流扩散方程没有扩散效果？
  • calrel  

  8
  8
  Posts
  3179
  Views

  

  不好意思，那个算例里面求解器我这面编译不了（openfoam-4.x and 5.x）.. 估计后来问题已经解决了吧
• 

  压力方程和解
  • 李东岳  

  2
  2
  Posts
  1483
  Views

  程

  这个是理论大坑啊。 如果压力基的求解器，求解的是速度和压力场，那么得用连续性方程来求压力。 如果是密度基的，连续性可以求密度，然后用状态方程求压力。
• 

  有意思的CFD算法类研究
  • 李东岳  

  9
  9
  Posts
  9076
  Views

  

  可压缩PBE方程形式的比较 咦、之前写的这个研究方向已经发出去了 SuperBee好像就会锐化界面。 Yes, SUPERBEE加了下风格式，是叫下风么？Down-wind 哈哈，东岳大神，感谢你做了这个平台 ... 感觉学习CFD就是要不断地重复才能领悟其中的真谛，当然有通俗易懂的资料会更为便捷。 @mohui 多谢支持，CFD任重道远比较耐玩！只要有兴趣，是能研究一辈子的，不像王者荣耀，玩几年到头了 哈哈 资料以后会持续更新，精力有限，稍微慢点
• 

  LES一个方程
  • 李东岳  

  5
  5
  Posts
  1651
  Views

  

  此问题同类于湍流动能耗散率的定义。湍流动能耗散率的定义为： \begin{equation} \varepsilon=2\nu_t\overline{S_{ij}S_{ij}} \end{equation} 经过各向同性假定后有： \begin{equation} \varepsilon=\nu_t\overline{\frac{\p u_i'}{\p x_k}\frac{\p u_i'}{\p x_k}} \end{equation}
• 

  大气泡小气泡不同的移动方向
  • 李东岳  

  1
  1
  Posts
  825
  Views

  

  Measurement of Lift Force of Air Bubble in Linear Shear Flow of Water Using L-Channel 在这篇文章中，作者研究了单一的不同粒径的气泡在水中的移动方向，主要研究升力系数问题。感兴趣可以看看。 下图就是实验，高速相机监控的不同直径气泡不同的切向移动。
• 

  高等数学的CFD相关问题
  • 李东岳  

  1
  1
  Posts
  1212
  Views

  

  页数按照pdf的页数，不是书上的页数 上册： https://www.jianguoyun.com/p/DVcobNoQ9s3ZBhjejV8 73页【直线运动的速度】：导数的直观描述，CFD方程这面全是导数，大部分为偏导数，后面会介绍 75页、76页【公式1-4】 138页【公式3-3】：泰勒公式，CFD很多方程都是基于泰勒公式推导出来的 下册 https://www.jianguoyun.com/p/DQnanmcQ9s3ZBhjijV8 23页：流量计算公式 71-75页：偏导数，CFD这面都是偏导，重点看一下，理解一下 77页提到了个拉普拉斯方程，随便扫一眼 78页全微分，跟CFD的物质导数关联 79页3-3 81页有个例子 112页介绍了梯度，但是感觉描述的不是很简单，可以看看 141-144页 重积分，为之后理解有限体积法做准备 224-225页 曲面积分的定义，为之后理解有限体积法做准备。其中225页有个图，比较形象 230-231页，看高数中对通量的定义 243-245页，高数中明确的对速度散度的定义
• 

  升力模型
  • 李东岳  

  1
  1
  Posts
  1017
  Views

  

  为什么我这个推导的反了？应该是大于0 。
• 

  ANSYS Fluent中QMOM的植入缺陷
  • 李东岳  

  4
  4
  Posts
  1931
  Views

  S

  @东岳 谢谢
• 

  一些CFD的方程
  • 李东岳  

  1
  1
  Posts
  1078
  Views

  

  张量标识预习测验 下面的方程取自各个不同的CFD教材。方程并没有采用代码的形式而是直接截图，这样更加直观的可以感受到文献里面方程这种多样性的书写方式。各位老师和同学需要能够拆分为分量形式。
• H

  动量方程为什么要忽略重力项？
  • huadian001  

  2
  2
  Posts
  1345
  Views

  

  大都忽略重力源项，哪位大神给解释下为什么忽略重力源项 重力对流体的作用就是引起水动力学压力，如果是密闭的流体，其不会引起流动变化。因此可以忽略重力项，这样压力边界条件也可以好给出。 一般考虑重力的主要是重力为一种重要的体积力，如密度不均引起的浮力，其中密度不均进一步由温度导致。 还有一种情况要考虑重力是模拟流体的自由落体。 其他情况可能跟具体的应用有关，如旋转系、海洋热风那种普遍都添加重力。 如果不忽略，怎样恰当的处理？ 如果考虑重力，通常求解的压力不是先前的压力p，而是压力p减去水力学压力$\rho g h$。从算法层面理解，重力引致下落，压力对这种效果抵消，采用这种压力计算的时候，可以省去这种来回下落-抵消的求解过程。求解出$p-\rho gh$之后，真实的压力即为$p-\rho gh+\rho gh$
• 小

  求助帖：关于高阶Boussinesq方程的推导
  • 小o0男  

  5
  5
  Posts
  1784
  Views

  小

  @东岳 好的，麻烦您啦！我在一点一点仔细研究下。如果能研究明白，万一以后也有人问，我再做个帖子出来♪(^∇^*)
• 

  i7-5820k超频（九州风神散热器）+ 使用OpenFOAM做CPU超频测试
  • 李东岳  

  7
  7
  Posts
  5646
  Views

  

  今天上到了4.4G，也还行
• B

  湍流的尺度的计算
  • benqing  

  2
  2
  Posts
  1674
  Views

  

  是否是真实的流场中湍流长度尺度呢？ Kolmogorov是一种假说：所有雷诺数足够高的湍流中，最小的尺度由一个普适性方程来决定，其参数为粘度和耗散率。这个假说建立在若干的假定基础之上，如各向同性。随后，有了Kolmogorov尺度。所以针对你的问题，至少各大教材都认为Kolmogorov尺度是最小尺度。 这个泰勒湍流尺度能否代表RANS求解中的尺度信息呢？ 是的 在大涡模拟中的湍流长度如何让计算？ LES里面的湍流长度通过$\Delta$控制。大于$\Delta$被解析，小于$\Delta$被模化。不确定是否可以通过$\Delta$定义LES的湍流长度。
• L

  关于努塞尔数问题请教
  • lsprxd  

  10
  10
  Posts
  3685
  Views

  

  @lsprxd $h=\frac{q}{T_w - T_b}$ 其中：$q$为热流密度，$T_w$为壁面温度，$T_b$为参考温度 这样可以吗？
• L

  在方腔自然对流中为什么大家都不分析压力的变化
  • lsprxd  

  4
  4
  Posts
  1382
  Views

  L

  @东岳 是的，我看的文献大多数是关注速度，温度和涡量场变化，几乎没有涉及压力场变化，于是我编程计算的时候也没有考虑，结果老师问我为什么，我都不知道如何回答了
• 

  层流模型指的是什么？
  • .J.  

  5
  5
  Posts
  3261
  Views

  

  @winsway_zero 有空整理一下文献，感觉发在CFD界比较合适 类似还有这个帖子 http://www.cfd-china.com/topic/1056/3d-or-2d-les