Skip to content
  • 求解颗粒沉降的方法选择问题

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

    @hurricane007 好的,非常感谢

  • 0 赞同
    1 帖子
    2k 浏览
    thegameT

    00d4f1b7-3cbe-4541-8fc1-e899bf8ec13a-image.png
    代码路径:
    channelIndex.H
    channelIndexTemplates.C

    想请教各位,他在channelIndex.H这里传的是const bool asymmetric = false,但是在模板实现channelIndexTemplates.C中又是const bool asymmetric。那么在初始化的时候,到底是用的哪个呢?

  • 颗粒液固两相流的界面怎么确定呢?

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

    那个有alpha场,可以用来定义界面。类似interFoam

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

    压力没有一个绝对值的话,必须需要给一个参考值。但是你是自己写的代码,不好debug

  • 0 赞同
    7 帖子
    8k 浏览
    chengan.wangC

    @星星星星晴 感觉nSector不能设为1,最后输出的sum(total mass),sum(average mass flow rate)是nSector>1的4倍左右。

    if (nSector_ > 1) { refDir = this->coeffDict().lookup("refDir"); refDir -= normal_[0]*(normal_[0] & refDir); refDir /= mag(refDir); } else { // set 4 quadrants for single sector cases nS = 4; vector tangent = vector::zero; scalar magTangent = 0.0; Random rnd(1234); while (magTangent < SMALL) { vector v = rnd.vector01(); tangent = v - (v & normal_[0])*normal_[0]; magTangent = mag(tangent); } refDir = tangent/magTangent; }

    可能跟源代码中nS = 4;有关系

  • 美国底层打工人的2021年终总结

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

    我是19年国庆后来的美国,基本上一直在干餐馆,中间短暂送过三周外卖。最初第一年干了11家餐馆,3家自助餐,5家日餐,3家中餐。2021年1月29号开始在这家中餐馆做企台的,就是做服务员,俗称端盘子。包吃住,底薪750$,剩下全靠小费,扣除跑菜的12%,和付给早上替我做side work的(摆桌子拖地打扫厕所),剩下的小费都是我的。

    上班时间是上午11点半左右到店里简单准备一下就开始接客打台,营业时间是11点半到晚上10点,打烊后还要收拾桌椅,录入小费,洗杯子,因为疫情,市政府允许室外就餐,摆了很多桌椅在门口大马路上,我要把椅子搬进来,把桌子归拢好,拿铁链子锁起来,最后负责锁门。离开店里大概要晚上11点左右了,所以工作时间大概是12个小时,每周6天,周二休息一天,旺季的时候只休半天。类似于996,不过我是11-11-6。

    工作内容就是打台,就是给客人点餐,上饮料,添茶续水,吃完买单,收桌子,然后翻桌接下一桌,很多时候是同时进行的。还包括接电话,把4个网络外卖平台的来单录入电脑,送到后厨。后来老板请了人负责跑菜,就是传菜,我就不用端盘子了,扣15%小费给跑菜的,后来有了寿司需要我自己跑菜,扣除的部分经常争吵降到了12%。旺季的时候请过一个老墨amigo负责带位,点饮料,上饮料,最后收桌子。那段时间我的工作就简化到点菜,买单,收桌子上的小费单。

    工作强度大概只能用每天的步数衡量了,我每天要走2万到3万步。旺季和周末五六日大概2万5到3万步,平均7天每天2万步出头。

    收入和支出情况,大概淡季6000多刀,旺季有几个月9000多。因为包吃住,每月的开销就是休息的那天去趟超市,买点吃的喝的,大概花三四十$。手机话费每月40,加油一次三十多,能跑三个多星期,上下班大概9个迈。车险合每月不到60,给儿子每月打生活费1000,偶尔每月去按摩一次加小费60块。

    这样一年下来,收入合计90341,支出19095,结余71246。跟在国内的时候收入几乎一样,不同的是单位从人民币变成了美元,从月光变成了能结余一大笔,这也是我在出国前了解了一些情况后,最终决定来美国的原因。

    21年3月底收到了移民局寄来的社安号和工卡(EAD,Employment Authorization Document)。面谈因为疫情还是等待中,所以也不算打黑工了。过两天找个会计师报税。

    为什么要来美国?我编辑了一大段,还是决定删掉了,因为涉及个人隐私,另外带有敏感的话题,不适合在公开论坛说。

    以后的计划。说实话没想好,走一步看一步吧。干几年存点钱,干不动了就去送外卖或者开卡车也说不准。实话说,不太想回去,没绿卡前回去了也就回不来了,当初的铁饭碗也没了,回去也找不到像样的工作了,因为学的文科,也没什么技能。

    另外,我的收入情况有点特殊,不是每个餐馆的企台都能挣这么多,正常情况是每月4000-5000。为什么我能挣的比别人多,因为店里就我一个全职企台,只要我能忙的过来,老板会一直让我干,我接不过来的桌子他才出手,也就再去请别人,因为再来一个人,都挣的少了,都不会继续干下去,都留不住。因为我能干,肯定,他也比较照顾我,忙的时候,有大桌他会让给我,他去接一个两个人的小桌,因为美国餐馆惯例,超过6个人就可以强制收取18%的服务费,运气好或者客人大方的,有的会在18%的基础上再签18%或者20%,多的时候一桌收到过一百四五十的小费。因为有人替我跑菜,我同时可以照顾7桌客人,平均下来每桌10刀,一个人头5刀。最多一天打了快600块。

    算下来,每天2万步,一年大概走了730万步,每天12个小时,一周6天,一年3700小时,平均一小时24.4刀,比加州最多时薪15刀还是高出不少,只有洛杉矶新毕业护士时薪的6成,所以一度我想存钱考护士,因为听说湾区护士平均年薪14万刀,有一阵很心动。就算读本科加速班,也要脱产读三年差不多,时间和学费上又让我打了退堂鼓,加上年纪也大了,马上40了。

    跟隔壁意大利餐馆的墨西哥小妹也交流过,他们生意好,他们单位时间可能挣的比我多,但是不是每天都排班,所以有的人一周要打三份工。因为老外老板会避免安排太多工时,这样雇佣的人少雇佣的时间短就不用给工人什么福利。

    其他底层人不同工种的收入情况也说一下,我做过送外卖,在洛杉矶跑doordash,类似于国内的某团,20年圣诞到21年元旦那周我还是生手(刚学会开车,刚拿到驾照刚买的车),第二周一周就挣了快1400,熟手使劲干一周能挣2000,加过华人外卖群了解的。在大城市开Uber比送外卖收入只高不低。

    了解过按摩,按摩店男的四五千,女的六七千。我去的那家床脚店,男工按脚一小时35,女工按全身一小时50,所以女的挣得多。夫妻工12000左右。

    看过一则旧金山搬家招聘月薪8000,淡季6000+,旺季10000+,不过一般人没那个力气和身子骨是干不了的。流调最辛苦那位大哥估计可以。

    根据华人职介所的招工广告和我了解的实际情况,中餐馆洗碗打杂3500,帮炒3500-4000,炒锅(厨师)4000-5000,前台接电话打带位包3500-3800,日餐寿司学徒3000-3500,助手4000,二手4500,寿司师傅5000。企台4000-5000。加州老墨洗碗15/小时。

    家政:住家阿姨,打扫卫生煮饭,15刀/小时,100-120/天,3000-3500/月

    物流公司叉车工人18-20刀/小时,经常加班5000+,送货司机200-300/天。

    装修小工150-170/天,中工200。

    结语,前几天跟朋友调侃道,早知道美国人傻钱多速来,还读什么研究生,上什么大学,跟福建人一样,中学毕业就来这边干就完了,现在不知道攒多少钱了。我也见过几个来美国赚了几套房甚至开过跑车的厨师,因为吃喝嫖赌抽,到60来岁还在干餐馆的。

    https://weibo.com/u/3764674343

  • 虚拟机双版本OpenFOAM怎么切换使用

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

    @李东岳谢谢东岳老师,已经搞定了在求解器窗口输入 of2106 就切换回V2106版本了:high:

  • 有关泡沫金属与石蜡传热传质问题

    Fluent
    1
    0 赞同
    1 帖子
    2k 浏览
    Z

    替代文字 我模拟的是一个矩形泡沫铜浸润石蜡的传热传质过程,现在模型验证一直没弄好。首先是编写了一个简单的UDF,采用热平衡模型,导入fluent后提醒我定义比热的那一块有错误,求各位前辈帮我看一下;其次,我添加了Boussinesq项后,优化了一点点液相云图中熔融线一直平行的问题,但是与论文中的还是有差距,且我的熔融线上下厚度保持一致(实际上应该下面厚很多),UDF中还需要添加什么项来解决这个问题呢?以下三张图分别是我的UDF、我跑出来的液相云图以及论文的液相云图,求各位帮我解惑,不甚感激。2.png 3.pn替代文字g 4.png

  • 关于QMOM矩初值的给定

    Algorithm
    26
    0 赞同
    26 帖子
    32k 浏览
    waltenW

    @东岳 我现在接触的还很浅。。以后碰到这方面问题再来请教~~

  • Icem CFD 几何拓扑报错

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

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

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

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

    我也发现了这个情况,并且写在了ns方程笔记里面。不过没有详细研究。我觉得是因为多相流密度差异导致的压差计算不同导致的。有空是应该研究一下。

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

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

  • InterFoam设置壁面Wall上不同的润湿性

    OpenFOAM
    7
    0 赞同
    7 帖子
    8k 浏览
    D

    @wsxfyy 谢谢,我试一下

  • 保持登录信息时发帖会显示Forbidden

    CFD彩虹条
    7
    0 赞同
    7 帖子
    9k 浏览
    李东岳

    我最近出现了这个情况。还有人有这个情况么,很奇怪

  • DM塔式服务器 cpu 参数

    CFD彩虹条
    1
    0 赞同
    1 帖子
    2k 浏览
    Y

    @李东岳 李博士,

    请问东岳流体官网的塔式工作站(DM0812, DM3.1)cpu 具体参数,像型号,主频,睿频等可否简要分享下。

    谢谢!

  • 0 赞同
    15 帖子
    17k 浏览
    I

    大家好,我来更新这个问题的相关进展。最近两天我查了些资料,基本上解决了这个问题,感谢大家的帮助,尤其是bestucan大佬。
    参考帖子关于相邻cell我写了个遍历进行暴力求解。
    (1)代码思路:对所有网格做遍历,找到气液界面网格,针对每一个气液界面网格都对所有网格再做一次遍历,计算两网格之间的距离,对满足距离要求的网格在markRegionNearInterface这个场上+1。最后,可以利用markRegionNearInterface>=1来筛选出所有满足条件的网格。
    (2)这段代码计算起来实在是费时间,因为我的需求(见贴名)不需要每一个时间步都运行这段代码,所以我在代码最外层套了一个时间步判断,可以实现多个时间步运行一次该代码。
    (3)我从其他地方抄了点代码写了两个接口,分别从transportProperties读取 加密时间步间隔 和 加密范围,其中加密范围是依据气泡半径写的无量纲范围,所以需要配合前面代码(未在本段代码中给出)求出的气泡半径radius.value()定义标记区域。

    代码如下

    //read timestep interval required for the following code, as the following code is very time consuming so I needn't it run for every timestep. const label calculateTimestepInterval = transportProperties.get<label>("calculateTimestepInterval"); //read information from transportProperties so that we can modify the size of marked region easily. scalar markRegionNearInterfaceSize = transportProperties.get<scalar>("markRegionNearInterfaceSize"); if (mesh.time().timeIndex() % calculateTimestepInterval == 0)//judge whether the timestep satisfy the timestep interval criterion. { markRegionNearInterface.field()=0;//reset the field to be zero first. Of course such field is claimed in createFields.H first. forAll(mesh.C(), i)//Cycle all meshes { if (alpha1[i]>0.0001&&alpha1[i]<0.9999)//judge whether the cell contains gas-liquid interface { vector centerOfSurfacePoint = mesh.C()[i];//store center information of such a interface cell forAll(mesh.C(), cellI)//Cycle all meshes again { scalar offset = mag(centerOfSurfacePoint - mesh.C()[cellI]);//distance between any cells and the interface cell if (offset <= markRegionNearInterfaceSize*radius.value())//judge the criterion, note the radius is calculated before with code which is not shown here { markRegionNearInterface.field()[cellI] += 1;//any cells which meet the criterion will be marked with value>0 } } } } }

    运行结果如下
    10652d23-2e1a-470d-a7ef-a8637883685a-image.png

    BTW
    (1) 因为我是个beginner,代码是东抄抄西抄抄弄出来的,上述代码可能有些冗余,如果大佬们有时间,想麻烦大佬们帮我看看我这些代码能从哪些角度优化优化。代码中的变量类型scalar vector label之类都是我一次次试错试出来的,也不知道这些变量类型用的合理不合理。
    (2)我在代码中写了两个链接到transportProperties的接口读取 加密时间步间隔 和 加密范围,测试中发现程序只会在刚开始运行时读取这些常量。请问各位大佬,有没有什么办法能否让程序运行中也能读取这些常量呢?因为假设运行过程中我想改动这两个常量,目前只能把计算停下来改了后重新继续算,不能边算边改。

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

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

  • 网格+粒子的两相流仿真

    Algorithm
    12
    0 赞同
    12 帖子
    16k 浏览
    JerryJ

    @东岳 谢谢李老师o(^o^)o

  • 0 赞同
    7 帖子
    9k 浏览

    @FluidGao 哦,我下来看看,谢谢

  • 导出网格数据

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

    @禾斗石开 fluent