Skip to content
  • 0 赞同
    1 帖子
    2k 浏览
    T

    teoria.png
    如图所示,这是非常经典的Le Méhauté (1976)波浪理论适用图,但我有点迷糊,主要有两方面疑问:

    图中红色部分看起来是Stokes三阶波/四阶波与五阶波的交接处,那么我应该选择用三阶波/四阶波还是五阶波呢?同样地,图中浅蓝色和五阶波也有交集,那么在既是浅蓝色又是五阶波的区域应该用cnoidal还是五阶波呢? 我还注意到图中五阶波下面还有一个流函数理论,但我看olaFlow的用户手册,用流函数的话需要输入一大堆参数(如下图所示),这些参数应该怎么选取呢?本菜鸡在此请教一下各位:135:
    2024-07-08 21-55-00 的屏幕截图.png
    2024-07-08 22-01-59 的屏幕截图.png
  • 输入公式方法

    已固定 Algorithm
    9
    0 赞同
    9 帖子
    19k 浏览
    李东岳

    @coolhhh 厉害厉害!!

  • 0 赞同
    18 帖子
    17k 浏览

    三维的类似情况,我最后一个打错字了。

  • 0 赞同
    5 帖子
    6k 浏览
    V

    @李东岳 感谢回复。我是想先写个简单的特定算例学习一下算法,格式边界条件啥的都直接写死,然后再慢慢做的通用一点,能搞出来肯定开源。

  • 关于dpm的udf.

    Fluent
    6
    0 赞同
    6 帖子
    9k 浏览
    L

    这个属于dpm的particle的fate了 但是fate里面没有trapped 哈哈
    所以实际上ABORT了也不会被dpm output header统计进去呀
    这个时候不妨用用tp->gvtp g开头的意思是global
    最开始的错误也很明显的 你不管怎么样都是返回的ACTIVE 自然也就不行了
    实际上还能PATH_END也可以做
    但是ABORT并没有移除particle的链表这条记录 所以有的时候可能你需要 可能你不需要
    最好还是再MARK一下remove吧

  • 论坛搜索功能下线了?

    CFD彩虹条
    3
    0 赞同
    3 帖子
    3k 浏览
    D

    @李东岳 我发了帖子之后才发现你之前发过通知了,不是催更,抱歉抱歉~:chitang:

  • icem非结构网格如何划分block

    Meshy
    16
    0 赞同
    16 帖子
    21k 浏览
    Y

    结构化网格非结构化网格按照网格是六面体还是体面体区分好像是不完全正确的。我看一个老师讲解说,结构化和非结构化区分是文件保存格式的区分,结构化是使用索引,非结构化是把每个网格都区分出来了。我看各位大佬上面回答可能是直接讨论的四面体网格,因为我最近也在学openfoam和icem,icem可以通过自动化分成四面体非结构网格,也可以分块化成6面体网格,但是我通常画完6面体网格后,icem有可以转化成非结构网格按钮,将画好的网格转化成非结构体网格(转化后依然是6面体)。
    所以如果你要话非结构化网格也可以分块(画了6面体转化一下),如果你说的是画四面体非结构化网格分块的话,这方面知识我就不知道了

  • 说说关于dpm系统的半透膜的问题

    Fluent
    3
    0 赞同
    3 帖子
    5k 浏览
    L

    因为时间比较紧,叙述不是很明了,现在按照李老师的建议,加上标点。

    一般dpm的碰撞判断基本是BC宏,这个不多说。

    碰撞算法fluent未透露,但是在ANSYS博客中间有一部分提及了,是采用了颗粒轨迹线与面相交的办法。这个办法是属于碰撞检测算法 line-surface。 这个方法非常的复杂。同时在CFD-ONLINE中间,也有人讨论过这个问题。 从我观点,这点Fluent的保密算法不可能是采用最原始的intersection算法;如果是的话,那么巨量的颗粒沉降计算很明显没有fluent这样的高效。

    好了,这是说了题外话,但是也跟接下来的相关。

    也就是当你不利用BC宏的时候,你需要写这样一个算法(line-surface)来判断 dpm 的截留。而这显然是在半透处理上是较为简单的。因为这样的话,你只要一个在某一个面上,甚至是一个解析几何的构成虚拟面上,就可以构造你的半透面。只要颗粒符合你的条件。 当然这种复杂在于你的颗粒将在第几步碰撞? 他能够在途中的每个网格的步数是多少?(假设颗粒在第一次进入一个cell内被捕获,这一般是现实的,因为你基本上只要处理碰撞表面网格的数据)。 那么又有一个很关键的问题: 你需要重新编制一个颗粒更新步的计算。 假设你是最简单的欧拉蛙跳,也显然你需要了解一个parcel的合外力,所以你必须了解整个dpm.h和其他几个.h的头文件。 那么,你将可以直接解决这个半透的问题。

    而如果你不写trap算法。那么很明显,你需要使用fluent所谓的DPM体系的BC宏。这个宏在颗粒碰撞你挂载的壁面的时候执行,记住这个宏执行将中断颗粒时间步。 而他本身是不出网格的。 因为出网格需要赋值新的P0值。(P0属于particle结构体而p属于tracked particle )。也就是说,当颗粒碰到这个面的时候,碰撞数据被计算,碰撞点被获取,颗粒必须更新他的P_v0值,而particle依然是在当前cell,且颗粒的fate只有active, end, abort, escape 四种。 你的trap可以使用abort,但是他的显示数值决于于p->gvtp, 如果要他在head中显示正常, 你还得做一番处理。

    最要的问题是,颗粒将在当前cell被分割,但是这个这个新的P0值依然是在当前cell内。 那么就是说,这些颗粒继续会被这个BC执行。 最终的结果是所有颗粒都被在封闭在这个cell内来回碰壁, 截留, 最后造成完全的incomplete。
    那么如何半透? 现在已经说明了在fluent实现这个半透可能的方法。

    设想一个现实的条件。颗粒一部分会被反弹(他的vo重新被设置) 一部分被捕获(path_abort) 一部分透过(记住你的时间步已经被分割, 而透过是不分割时间步的)。但就是这样最简单的日常物理现象,因为时间步被分割,碰撞算法被黑盒子,依然是十分困难的。(并且fluent的更新过程并不是属于经典的更新过程,时间步的更新会被更小的切分,根据你的设定,并且每一个切分步,你的UPDATA_SCALAR宏都会被执行。也就是说一个理论上的碰撞算法,很可能并不具备真正的可实现性。)

  • ICEM版本问题

    Meshy
    1
    0 赞同
    1 帖子
    3k 浏览
    Y

    请问大家,有个ICEM19.4导出的.blk文件,在新版本(2022)中导入显示版本不匹配,但是.tin文件就可以,有没有什么好办法可以查看.blk文件?

  • 0 赞同
    1 帖子
    2k 浏览

    3d408cf3-836a-4aef-bb0b-a370d462f7e5-image.png 如图,左边的上部是大气边界,右边是封闭壁面。如果不在左边画和右边一样的细网格,得到的网格质量很差。如果画了细网格,则计算速度很慢。有更好的处理方式吗

  • 关于LES模拟中采用大时间步长的问题

    Algorithm
    5
    0 赞同
    5 帖子
    6k 浏览
    T

    @李东岳 在 关于LES模拟中采用大时间步长的问题 中说:

    你可以放大库朗数,或者你可以用URANS更快。这些都没问题。真正分析的时候库朗数放小就行。并且库朗数一般是局部库朗数的最大值。如果1000个网格里面有1个网格特别小。会出现999个网格的库朗数都在1以下,但这1个网格在10左右。都不是问题。

    好的,谢谢李老师!我明白了

  • 0 赞同
    1 帖子
    2k 浏览
    李东岳

    背景:我要做自己的CFD软件,有生之年,跟商软来一次硬刚!

    我这面目前要做的东西是基于OpenFOAM的CFD软件。OpenFOAM是一个国外的开源的CFD求解器,类似linux内核或安卓。我这面就是要基于国外的OpenFOAM内核,做一个中文的界面出来。就类似小米的MIUI手机系统。软件做出来之后就是国产CFD软件。但是内核是基于OpenFOAM。具体为什么要基于OpenFOAM而不是自己写。主要有2个原因:1)自己写出来的CFD代码没人信,2)所需要的时间要三五倍以上。这也就类似目前的C++,C++是国外提出来的,目前大可以放心使用,没必要自己从头开发一个编程语言。

    在有了高度开发化的CFD求解器核心之后,在顶层GUI设计的过程中,要充分的做一个新东西出来:1)不能基于老套路。要按照CFD用户最方便的角度去重新设计。2)要把CFD相关的经验融合进去。做智能化,最易用的操作。目前的CFD软件还不够智能,比如商软,她会给你提供一系列的湍流模型可以选择,但是不会给你提供一个最优解。那这个CFD软件在做的时候,要克服相关的缺点。

    软件做出来之后,基于我对产品的前沿性的设计。我有充分的信心和十足的把握去推动。软件要按照GPL协议,必然是要提供免费试用版。但同时,可以提供定制版,这部分可以获取利润。目前国内CFD软件的打法,都是直接推收费版软件。我觉得在没有人用的前提下,这些软件是推不出去的。可以赚些小钱,但是不会做成一个国家层面的具有影响力的工业软件。

    我对这个产品具有充分的信心把它推出去。到时候只要国人用CFD,那除了商软,开源openfoam,我希望大家能想起这个。这么多年来国家一直在推国产工业软件,我国确实也有,但是没有一个让大家能说得出手。我国有些软件做的不错了,比如金山软件可以抗衡office,中望CAD可以抗衡SolidWorks。但CFD领域却没有人能说我国的xxx能与国外抗衡。我要做的,就是在这个空白的地方做出一个产品。这个软件对标的市场就是西方国家在国内售卖的CFD商业软件。

    对合作方要说的:

    要充分理解我做CFD软件的目的。我要做出一个最适合国内工业界的普适性CFD软件。希望真正的能让大家用起来。这建立在软件充分的“好用”的基础上。用户占有率在第一位。 也希望能认识到,免费产品未必不可盈利。如微信、抖音、快手、QQ、淘宝等。另一方面。工业软件投入的周期非常长,资本可以做一定程度的催化剂,但底层算法开发以及软件的更新迭代需要长期的坚守。

    捐赠

    我个人认为这个事情很难做成。主要就是缺少一个前端,普通前端的工资不低,一年30万打底。好一点的一年得更多。咱也没有这个钱一直去投入。目前尚且没有国家基金以及民间资本对这个项目进行支持。于是我决定尝试最后一条路,虽然希望也很渺茫,那就是捐赠。

    我个人的想法是在最短的时间内,募集最少40万元,来支付一个前端的工资。当然在法律层面可能需要斟酌(比如一年期聘用关系?)同时,为了保证资金的透明。我这面可以做以下承诺:

    2024年9月1日,捐赠达到40万元,我自己将出10万元,凑够50万元,用于第一年的业务支出(40万元用于发工资,10万元用于成立公司组办公室等杂项)。同时继续募集资金用于下一年的支出。

    2024年9月1日,不足40万元。项目终止。资金退回给捐赠者。

    如果在此期间中彩票拿到天使投,资金退回给捐赠者。

    如果在此期间有合作方可以投人,那表明不需要额外资金招聘,资金退回给捐赠者。

    捐赠期间,款项进出都将在本帖公示。

    捐赠方式

    如果是一些个人的捐赠,请参考下方的微信支付二维码。

    如果是一些机构的捐赠。可以通过合作的方式来进行。例如,我这面有服务器售价5万元,合作方若愿意捐赠1万元,可走6万元的合同。其中差价1万元则进入本项目捐赠资金。其他的与我这面的其他项目合作都可以,如CFD项目、CFD课程之类。

    也希望大家多多宣传。CFD这面我可以弄,最需要的就是一个能做事情的前端。谢谢大家。:xiexie:

    资金列表 日期 机构/个人/匿名 金额 备注
  • snappyHexMesh 绘制网格报错

    OpenFOAM
    4
    0 赞同
    4 帖子
    4k 浏览
    Y

    好的
    我试试

  • 0 赞同
    5 帖子
    6k 浏览
    L

    @李东岳 这个案例我用的商软,发现空隙率0.01和0.5差别不大,影响最大的是阻力系数,我试了这两个空隙率,水都能很顺畅的流出来

  • 0 赞同
    4 帖子
    6k 浏览
    凤雏岂能消逝

    @ibelief 非常感谢

  • N-S方程前传:普罗米修斯

    Algorithm
    5
    0 赞同
    5 帖子
    8k 浏览
    B

    居然还发公众号了

  • Fluent报错,求帮忙看一下

    Fluent
    4
    0 赞同
    4 帖子
    6k 浏览
    西湖冷月西

    一方面是动网格变形或更替方式不合适,另外我觉得是中间细小的几何部位网格有点粗糙,可以试着加密局部网格,毕竟2D计算,网格资源消耗也不是很大。

  • mesh画的翼型的网格求助

    Meshy
    1
    0 赞同
    1 帖子
    4k 浏览
    L

    烦请各位大佬指点一下,workbench里面mesh画的带角度的翼型的网格,为什么边界层厚度不一样啊?怎么处理一下呢?上面的感觉太窄了

    e9d97488-fa0f-453d-9f2b-1884d5136695-image.png

  • 0 赞同
    2 帖子
    3k 浏览
    C

    您好,我也在研究topset和refinemesh的问题,您知道在blockmesh生成一个梯形(具有斜面)的网格时,使用refinemesh会出现报错是什么原因吗?

  • snappyHexMesh画网格如何处理尖角和狭缝

    Meshy
    15
    0 赞同
    15 帖子
    15k 浏览

    @李东岳 谢谢东哥,就是我把resolveFeatrueAngle调大之后画网格没有捕捉这条特征线,画几何的时候因为在两个不同平面上,这条特征线是存在的