Skip to content
  • 0 赞同
    4 帖子
    4k 浏览
    李东岳

    挺好,老铁

  • 0 赞同
    5 帖子
    8k 浏览

    @zhangxc0223
    Hi
    因为并没有做deterministic的模拟,所以code也基本没看过。。
    stochastic collision 在计算collision probablity的时候是有限制的,具体你可以看一下O‘roucke的论文。具体是怎么定义的,怎么得到一个non-dimensional collision probablity的,还有如何决定这两个parcel发生了碰撞。

    虽然碰撞在这个instantaneously 随机的,但是在统计学上,比如你有1M次碰撞的话,总的碰撞的结果就是相似的。
    这也就是为什么stochastic碰撞会比deterministic相对来说’便宜'的原因,通过统计学的方法达到相似的碰撞结果,但是不需要追踪所有的parcel。

    lagnrangian maxCo 存在的意义,个人认为,在一定程度上可以理解为再一次离散了每个euler time step,所以maxCo一定程度上决定了每个lagrangian time step的大小,和Euler 的Co 定义是相似的。

    碰撞是不会被错过的,因为deterministic理论上来说是追踪所有的parcel,所以下一步这个parcel会在哪,走的路径应该也是计算的,所以不会存在错过的,stochastic collision 会计算两个parcel的碰撞概率,当概率很高的时候,也是不会错过的,而且就算是这次错过了,在茫茫parcel海中,总会碰到的。。统计学上是不会错过的

    以上仅为个人理解

  • Arbitrary-Lagrangian-Eulerian (ALE) 是什么?

    Algorithm
    3
    0 赞同
    3 帖子
    6k 浏览
    小狗狗

    参考文献
    Demirdzic, I.; Peric, M., Space conservation law in finite volume calculations of fluid flow.

  • 0 赞同
    9 帖子
    19k 浏览
    X

    @东岳 流沙老师的教程确实很良心

  • ICEM multi zone 大家用过吗?

    Meshy
    2
    0 赞同
    2 帖子
    5k 浏览
    李东岳

    你可以试试,支持性不是很高。对于简单地几何,玩一下挺有意思的。

  • ICEM非结构网格的Max element问题

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

    请教一个比较初级的问题,查了帮助文档没有找到想得到的答案。
    ICEM采用Octree方法划分非结构网格的时候,帮助文档和网上的一些教程中都写到:“建议全局网格参数设置中的Max element值要是2的幂”
    这里2的幂是指:2的正整数幂?2的正数幂?2的实数幂?
    实际测试过感觉比如填0.05啥的也可以划分,网格也确实要比2的密很多。所以想问下各位高手这里这个数值填写应该有什么经验或者规律呢?尤其是模型尺寸很小的话,有时候确实需要很小的尺寸啊~

  • 关于MRF方法的理解的疑问

    Algorithm
    3
    0 赞同
    3 帖子
    6k 浏览
    Q

    86f6fd36-007a-42a5-b105-de5d4fcefffa-image.png
    cited:http://openfoamwiki.net/index.php/See_the_MRF_development

    图为在旋转坐标系内解绝对速度的控制方程,其中 \vec \omega \times \vec U_i 这项称为科氏力,旋转坐标系内独有的力。他的效果如下面这个视频所示https://haokan.baidu.com/v?vid=4227630901810560602&pd=bjh&fr=bjhauthor&type=video
    我的理解是可以将上面的式子理解为绝对速度产生组成的动量的输运方程,其中的通量是相对通量(\vec U_R \vec U_I),同时由于动量是一个矢量,因此在旋转坐标系中还受到科氏力的作用。等式坐标是压力的贡献,和由于运输的物理量U_I的扩散贡献。

    因此,我在思考是不是诸如湍动能,湍流耗散率这些标量,也可以认为在旋转坐标系内,运输物理量为湍动能,湍流耗散率这些标量,通量为相对速度的通量,但是由于这些量都是标量,因此没有多的项需要添加。另外比如SA模型中的发生项需要计算到voricity的大小,此处的voricity是不是也应该用相对速度?

  • 关于openfoam与结构动力学的问题

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

    @李东岳 在 关于openfoam与结构动力学的问题 中说:

    楼上说的对。小变形比较简单,也可以自己写一个

    好嘞,谢谢李老师!

  • 一颗们组团了

    CFD彩虹条
    1
    0 赞同
    1 帖子
    2k 浏览
    李东岳

    @一颗李子
    @一颗栗子
    @一颗橙子
    @一颗鸭蛋

    :ok: :ok: :ok:

  • AIAA会议找同伴

    CFD彩虹条
    2
    0 赞同
    2 帖子
    4k 浏览
    nuaa_zpN

    我师兄也投了这个会议,正在找小伙伴,我可以给你联系一下。

  • 关于OpenFOAM矩阵求解器

    Algorithm
    9
    0 赞同
    9 帖子
    20k 浏览

    @micro 网格数量越多,网格尺度越小,计算肯定更加不稳定,您可以查一下CFL数

  • 0 赞同
    3 帖子
    4k 浏览
    两月三年

    就是我实际网格模型是左边的哪一个,是右边的模型简化的结果,
    现在想要任意找一个一个过对称轴的平面,不知道如何选取在这个二维fluent中如何选取

  • 关于Fluent PBM模型 UDF的问题

    Fluent
    4
    0 赞同
    4 帖子
    7k 浏览
    waltenW

    @东岳
    然后我在cfd-online上看到有人回答过相似的问题。
    2019-05-09_114138.jpg
    网址:
    https://www.cfd-online.com/Forums/fluent/103074-population-balance-modeling-pbm-ansys-fluent-7.html
    他的意思是不是就是在fluent在计算合并时,都是从最小的气泡尺寸来计算合并的。
    不会跟踪网格里究竟是哪几个直径的气泡在碰撞。
    比如我设定了10个bin。d_1 d_2都是从我的bin9(尺寸最小的bin)开始计算碰撞,得到其他的bin的数密度(bin 8/7/6/5/4/3/2/1/0)。
    我个人认为,这样处理的话,合并会很慢,不符合实际。

  • 博士论文盲审意见

    CFD彩虹条
    1
    1 赞同
    1 帖子
    2k 浏览
    李东岳

    微信图片_20190306174704.png

    微信图片_20190306174712.png

    ....

  • 如何提高混合网格的质量

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

    @simthere 这个非结构网格可以在结构网格上生成吧,不需要强制合并节点,感觉这样会好点,强制合并出来的太差了

  • markdown测试

    CFD彩虹条
    1
    0 赞同
    1 帖子
    3k 浏览
    李东岳

    http://pandao.github.io/editor.md/

    主要特性 支持“标准”Markdown / CommonMark和Github风格的语法,也可变身为代码编辑器; 支持实时预览、图片(跨域)上传、预格式文本/代码/表格插入、代码折叠、搜索替换、只读模式、自定义样式主题和多语言语法高亮等功能; 支持ToC(Table of Contents)、Emoji表情、Task lists、@链接等Markdown扩展语法; 支持TeX科学公式(基于KaTeX)、流程图 Flowchart 和 时序图 Sequence Diagram; 支持识别和解析HTML标签,并且支持自定义过滤标签解析,具有可靠的安全性和几乎无限的扩展性; 支持 AMD / CMD 模块化加载(支持 Require.js & Sea.js),并且支持自定义扩展插件; 兼容主流的浏览器(IE8+)和Zepto.js,且支持iPad等平板设备; 支持自定义主题样式; Editor.md

    x x

    目录 (Table of Contents)

    [TOCM]

    [TOC]

    Heading 1 Heading 2 Heading 3 Heading 4 Heading 5 Heading 6 Heading 1 link Heading link Heading 2 link Heading link Heading 3 link Heading link Heading 4 link Heading link Heading link Heading link Heading 5 link Heading link Heading 6 link Heading link 标题(用底线的形式)Heading (underline) This is an H1 This is an H2 字符效果和横线等

    删除线 <s>删除线(开启识别HTML标签时)</s>
    斜体字 斜体字
    粗体 粗体
    粗斜体 粗斜体

    上标:X<sub>2</sub>,下标:O<sup>2</sup>

    缩写(同HTML的abbr标签)

    即更长的单词或短语的缩写形式,前提是开启识别HTML标签时,已默认开启

    The <abbr title="Hyper Text Markup Language">HTML</abbr> specification is maintained by the <abbr title="World Wide Web Consortium">W3C</abbr>.

    引用 Blockquotes

    引用文本 Blockquotes

    引用的行内混合 Blockquotes

    引用:如果想要插入空白换行即<br />标签,在插入处先键入两个以上的空格然后回车即可,普通链接。

    锚点与链接 Links

    普通链接

    普通链接带标题

    直接链接:https://github.com

    锚点链接

    GFM a-tail link @pandao

    @pandao

    多语言代码高亮 Codes 行内代码 Inline code

    执行命令:npm install marked

    缩进风格

    即缩进四个空格,也做为实现类似<pre>预格式化文本(Preformatted Text)的功能。

    <?php echo "Hello world!"; ?>

    预格式化文本:

    | First Header | Second Header | | ------------- | ------------- | | Content Cell | Content Cell | | Content Cell | Content Cell | JS代码 function test(){ console.log("Hello world!"); } (function(){ var box = function(){ return box.fn.init(); }; box.prototype = box.fn = { init : function(){ console.log('box.init()'); return this; }, add : function(str){ alert("add", str); return this; }, remove : function(str){ alert("remove", str); return this; } }; box.fn.init.prototype = box.fn; window.box =box; })(); var testBox = box(); testBox.add("jQuery").remove("jQuery"); HTML代码 HTML codes <!DOCTYPE html> <html> <head> <mate charest="utf-8" /> <title>Hello world!</title> </head> <body> <h1>Hello world!</h1> </body> </html> 图片 Images

    Image:

    Follow your heart.

    图为:厦门白城沙滩

    图片加链接 (Image + Link):

    图为:李健首张专辑《似水流年》封面

    列表 Lists 无序列表(减号)Unordered Lists (-) 列表一 列表二 列表三 无序列表(星号)Unordered Lists (*) 列表一 列表二 列表三 无序列表(加号和嵌套)Unordered Lists (+) 列表一 列表二 列表二-1 列表二-2 列表二-3 列表三 列表一 列表二 列表三 有序列表 Ordered Lists (-) 第一行 第二行 第三行 GFM task list GFM task list 1 GFM task list 2 GFM task list 3 GFM task list 3-1 GFM task list 3-2 GFM task list 3-3 GFM task list 4 GFM task list 4-1 GFM task list 4-2 绘制表格 Tables 项目 价格 数量 计算机 $1600 5 手机 $12 12 管线 $1 234 First Header Second Header ------------- ------------- Content Cell Content Cell Content Cell Content Cell First Header Second Header ------------- ------------- Content Cell Content Cell Content Cell Content Cell Function name Description ------------- ------------------------------ help() Display the help window. destroy() Destroy your computer! Left-Aligned Center Aligned Right Aligned :------------ :---------------: -----: col 3 is some wordy text $1600 col 2 is centered $12 zebra stripes are neat $1 Item Value --------- -----: Computer $1600 Phone $12 Pipe $1 特殊符号 HTML Entities Codes

    © & ¨ ™ ¡ £
    & < > ¥ € ® ± ¶ § ¦ ¯ « ·

    X² Y³ ¾ ¼ × ÷ »

    18ºC " '

    Emoji表情 :smiley:

    Blockquotes :star:

    GFM task lists & Emoji & fontAwesome icon emoji & editormd logo emoji :editormd-logo-5x: :smiley: @mentions, :smiley: #refs, links, formatting, and <del>tags</del> supported :editormd-logo:; list syntax required (any unordered or ordered list supported) :editormd-logo-3x:; [ ] :smiley: this is a complete item :smiley:; []this is an incomplete item test link :fa-star: @pandao; [ ]this is an incomplete item :fa-star: :fa-gear:; :smiley: this is an incomplete item test link :fa-star: :fa-gear:; :smiley: this is :fa-star: :fa-gear: an incomplete item test link; 反斜杠 Escape

    *literal asterisks*

    科学公式 TeX(KaTeX)

    $$E=mc^2$$

    行内的公式$$E=mc^2$$行内的公式,行内的$$E=mc^2$$公式。

    $$(\sqrt{3x-1}+(1+x)^2)$$

    $$\sin(\alpha)^{\theta}=\sum_{i=0}^{n}(x^i + \cos(f))$$

    多行公式:

    \displaystyle \left( \sum\_{k=1}^n a\_k b\_k \right)^2 \leq \left( \sum\_{k=1}^n a\_k^2 \right) \left( \sum\_{k=1}^n b\_k^2 \right) \displaystyle \frac{1}{ \Bigl(\sqrt{\phi \sqrt{5}}-\phi\Bigr) e^{ \frac25 \pi}} = 1+\frac{e^{-2\pi}} {1+\frac{e^{-4\pi}} { 1+\frac{e^{-6\pi}} {1+\frac{e^{-8\pi}} {1+\cdots} } } } f(x) = \int_{-\infty}^\infty \hat f(\xi)\,e^{2 \pi i \xi x} \,d\xi 绘制流程图 Flowchart st=>start: 用户登陆 op=>operation: 登陆操作 cond=>condition: 登陆成功 Yes or No? e=>end: 进入后台 st->op->cond cond(yes)->e cond(no)->op 绘制序列图 Sequence Diagram Andrew->China: Says Hello Note right of China: China thinks\nabout it China-->Andrew: How are you? Andrew->>China: I am good thanks! End
  • openfoam中cyclic周期性边界的问题

    OpenFOAM
    34
    0 赞同
    34 帖子
    54k 浏览
    W

    @dzw05 请问老师,您这边清楚interface并行边界如何处理吗?我想从openfoam中拿到完整矩阵数据,但是并行的话只能拿到各个线程对应的子矩阵,如何知道并行边界的关系将其组装成一个完整的矩阵

  • CFDer专用春联

    CFD彩虹条
    4
    0 赞同
    4 帖子
    6k 浏览
    I

    psb.jpg

    春风催碧叶
    夜灯盼文章
    横批:诸事如意

    改了一下
    春风催碧叶
    夜灯盼收敛
    横批:哈哈哈哈

    碧叶双关
    收敛双关

  • 0 赞同
    13 帖子
    23k 浏览
    W

    @李东岳 哦哦,没问题了。:142:

  • CFD中文网 定位与规范

    已固定 公告
    101
    0 赞同
    101 帖子
    661k 浏览
    H

    我写 固我在,嘿嘿