Skip to content
  • C++11数组初始化

    Algorithm
    4
    0 赞同
    4 帖子
    7k 浏览
    R

    @wwzhao

    简单的感受了一下两种赋值方式,效率确实很大提高。

    #include <stdio.h> #include <string> #include <chrono> #include <iostream> using namespace std ; class Timer { private: // Type aliases to make accessing nested type easier using clock_t = std::chrono::high_resolution_clock; using second_t = std::chrono::duration<double, std::ratio<1> >; std::chrono::time_point<clock_t> m_beg; public: Timer() : m_beg(clock_t::now()) { } void reset() { m_beg = clock_t::now(); } double elapsed() const { return std::chrono::duration_cast<second_t>(clock_t::now() - m_beg).count(); } }; class SomethingBefore11 { private: int m_array[5]; public: SomethingBefore11() // zero the member array { // If we want the array to have values, we'll have to use assignment here m_array[0] = 1; m_array[1] = 2; m_array[2] = 3; m_array[3] = 4; m_array[4] = 5; } }; class SomethingAfter11 { private: int m_array[5]; public: SomethingAfter11(): m_array { 1, 2, 3, 4, 5 } //zero the member array { } }; int main(){ Timer tBefore11; SomethingBefore11 m_array_before11; std::cout << "Time elapsed: " << tBefore11.elapsed() << ‘n’; Timer tAfter11; SomethingAfter11 m_array_after11; std::cout << "Time elapsed: " << tAfter11.elapsed() << ‘n’; return 0; } [xx OFtutorial0_helloWorld]$ whatAboutThisGuy Time elapsed: 7.506e-06 Time elapsed: 1.47e-07 [xx OFtutorial0_helloWorld]$ whatAboutThisGuy Time elapsed: 8.664e-06 Time elapsed: 1.9e-07 [xx OFtutorial0_helloWorld]$ whatAboutThisGuy Time elapsed: 7.646e-06 Time elapsed: 1.89e-07
  • buoyantBoussinesqSimpleFoam层流求解问题?

    OpenFOAM
    1
    0 赞同
    1 帖子
    2k 浏览
    D

    大家好,请教一个问题,我想用buoyantBoussinesqSimpleFoam求解器求解层流问题。将湍流参数字典设置为层流,我发现仍然需要在transportProperties字典中提供湍流Prt数?这里的Prt数参与计算吗?
    同时,0字典中,也仍需要提供alphat文件。

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

    0_1516841984913_捕获.PNG
    动图1:河床沙石冲刷运动

    河流的流动会改变河床的形状。底部的沙子以及大石头通常会分为两个区:

    顶层的bed-load zone(有研究河海的么,这个怎么翻译?),其中的颗粒会随之移动; 底层:沙子以及石头会被困住,但是会缓慢的爬行;

    在碎石较多的河床中,较大的沙石倾向于堆积在顶层,这种现象被称之为armoring。个人感觉之所以称之为armoring是因为较大的沙石趋向于在表面,类似一种武装。

    实验表明,在这个区域,大沙石有一个净得向上的运动速度,小颗粒具有一个小的向下的移动速度。 目前对于这种现象的原因尚不清楚。不过现存了一些潜在的理论。例如支持程度比较高的巴西果效应(Brazil Nut Effect)。

    0_1516841984913_捕获.PNG
    动图2:有人晃荡大米里面的密度差不多的小石头

    替代文字
    动图3:进行的计算模拟

    巴西果效应是指如果把两种颗粒的混合物置于容器中,然后施加外加的振荡,体积比较大的颗粒会上升到表层,而较小的颗粒会沉降到底部。关于这种古老的效应的动力学机制至今仍众说纷纭。

    DEM模拟:文中作者采用LIGGGHTS对这个现象进行了模拟,下图可以看出,随着时间的推移,大颗粒漂浮出来。

    0_1516841984913_捕获.PNG
    图4:Nature Communication中的LIGGGHTS模拟

    看起来很有意思 :cheeky:

    本文来自Nature Communication

  • 0 赞同
    18 帖子
    23k 浏览

    @tens 有可能下午来的 或者晚上来的有问题。
    咱讨论下一话题吧:mihu:
    :tishizi:

  • 求助:FLUENT模拟二维机翼不收敛

    Fluent
    7
    0 赞同
    7 帖子
    10k 浏览
    N

    @ildcfder 最后选用了耦合算法收敛了,谢谢提醒,一开始是密度基,后来改成压力基啦

  • 0 赞同
    6 帖子
    8k 浏览
    C

    @bestucan 嗯嗯,感谢大佬。检查过,check和Quallity其实都没有问题,再次感谢~

  • PBM之QMOM中最后一步。。。

    Algorithm
    11
    0 赞同
    11 帖子
    17k 浏览
    AstremonA

    @李东岳 http://dx.doi.org/10.1016/j.compchemeng.2015.09.008
    这篇文章对三种重构方法做了对比。

  • 0 赞同
    3 帖子
    5k 浏览
    Y

    已经解决。
    (在一个视频中提到这一点)网格统计问题有三种情况:
    第一种,是显示问题的网格数目很少,就几个,占总体的0.0002%。对于这种情况直接更改限定值。
    第二种,是显示问题的网格数目很多,占总体的超过10%,甚至更多,对于这种情况需要更改MeshData的参数设置,重新对网格进行绘制。
    必须保证,不能出现负网格。

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

    @yingqing 不知道为啥,有一次成功了,在弄就不能成功了。节点的位置进行了设置,但网格的节点与我设置的节点不重合。就像这样
    177924ad-d438-468f-8cf5-a612fe049465-image.png

  • DPM入射质量问题

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

    @东岳 谢谢李博士,我感觉好像是不算的,因为在后处理的过程中,MassFlow的进口和出口的所有流体的质量流量值和液相的相等,但是固相的就是固相的质量流量,好像不包括,好像DPM模型中的液相和固相是采用分离流计算的,那我们在设置出口质量流量的时候就只考虑液相的质量流量,不过可能这个是这个模型只用于计算固相体积分数较少的情况

  • Latex表格自动换行

    Algorithm
    6
    0 赞同
    6 帖子
    9k 浏览
    李东岳

    谢谢各位,解决了:mihu:

  • 超算安装openfoam问题

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

    @东岳 感谢东岳老师指点:chitang:

  • Icem CFD 几何拓扑报错

    Meshy
    2
    0 赞同
    2 帖子
    4k 浏览
    bestucanB

    应该和左下角的 tolerance 相关。
    这tolerance都那么小了,还能12点覆盖

    估计是把一个几何体修复成一个点了,所以崩了

  • 有限体积法的强弱形式

    Algorithm
    7
    0 赞同
    7 帖子
    11k 浏览
    H

    @李东岳 大神,这个我真没注意。好像不多。我是土木专业的,不是搞固体力学的。您可以再找找。

  • 并行计算中求最大值

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

    思路:如果是内部场的值,调取a最大值的cell label,然后拿出这个label的b值就可以了。对于边界场同理。

  • 0 赞同
    3 帖子
    5k 浏览
    A

    谢谢您的解答!目前我做的是流固多场耦合,固体是阀门(上百个零件),流体是阀门及附属管道内区域(开度约10~20°),介质是可压缩空气,目前已知的只有质量流量入口,换算成管道入口平均流速约400~500m/s,若把流体当不可压计算,那么在阀板附近将达到一万多米每秒。如果当做可压缩计算,结果总是发散:我把流体密度设置为idealgas,其他参数恒定,入口用质量流量,入口恒温,压力出口为大气压。目前的发散状态是阀板附近流体温度不断降低到1K以下,超过fluent最低温度限制了。初步怀疑是阀板附近网格问题,当然,这样的设置是否合理,还未经过严格的论证,请问您还有什么建议么?

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

    any official links?

  • fluent模拟气泡/液滴聚并破碎求助

    Fluent
    18
    0 赞同
    18 帖子
    30k 浏览
    P

    厉害了!我好奇问一下,这篇文章里的控制方程一共有多少个?对于三维的情况,下面的方程总数是8个吗?(连续方程2个,动量方程6个)也就是说作者用求解的是双流体模型吗?是不是还缺一个相分数的输运方程?
    07c4a5d3-1477-404c-b074-5afecb87284b-image.png

  • Tecplot中涡环作图问题

    CFD彩虹条
    3
    0 赞同
    3 帖子
    3k 浏览
    凤雏岂能消逝

    @zhouxu 感谢您的回复,我尝试了一下,的却可以做出来。这部分本就有代码在计算。

  • 搞CFD的都能玩点啥?

    Algorithm
    4
    0 赞同
    4 帖子
    7k 浏览
    李东岳

    @youmengtian
    你用的什么显示器 :cheeky: