应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » DPM设置
191/1返回列表
查看: 5509  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

cure727

铜虫 (初入文坛)
已领完
DPM设置
领取红包 (小木虫手机app专属红包)

扫一扫，下载小木虫客户端

各位大佬，冒昧打扰，还请见谅
我想问一下有关DPM设置的问题，首先我想模拟的是人体对着某容器一次性呼出的颗粒物的轨迹以及两个出口的沉积量，模拟思路为稳态流场＋非稳态颗粒追踪，先将流场计算至收敛，再进行DPM设置，具体设置如下:
双向耦合，勾选update DPM source every flow iteration，DPM iteration interval设为20，开启非稳态颗粒追踪，颗粒时间步长为0.001，number of times steps设为45，最大时间步数设的很大(5×10的12次方)，step length factor默认为5，考虑萨夫曼升力，虚拟质量力，压力梯度力
injection设置:面注射，start time为0秒，end time为5秒，粒径为统一粒径5μm，注射速度为2.03米每秒，开启随机轨道模型
边界条件设置:入口为escape，两个出口均为trap，壁面为默认reflect
设置好以上后，直接进行计算
具体问题如下:
(1)颗粒时间步长下面的number of times steps是不是为了确保颗粒在每一时间步里收敛，应该设的大一些?
(2)粒子注射开始时间设为0可以吗？我的模型很小，长度也就20厘米的样子，所以注射时间差设为5秒合理么?
(3)原本流场已收敛，但是加入颗粒物之后，残差曲线突然上升很多，正常吗?是应该再计算算至稳定收敛即计算结束?
烦请各位大佬和同行可以解答一下，小弟在此感激不尽

DPM设置


发自小木虫Android客户端
回复此楼

» 猜你喜欢

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2021-12-01 09:58:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

17854252951

新虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
设置的也不对啊，你这不是非稳态的设置吗

发自小木虫IOS客户端
赞一下 回复此楼
7楼2021-12-01 11:58:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cure727

铜虫 (初入文坛)
引用回帖:
7楼: Originally posted by 17854252951 at 2021-12-01 11:58:56
设置的也不对啊，你这不是非稳态的设置吗

我设置的就是非稳态，稳态的流场+非稳态的颗粒追踪

发自小木虫Android客户端
赞一下 回复此楼
8楼2021-12-01 12:34:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xrs333

木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
（1） 稳态流场＋非稳态颗粒追踪时的颗粒轨迹追踪过程：颗粒轨迹追踪方式为，每隔Number of Continuous Phase Iterations Per DPM Iteration个连续相迭代步数，DPM求解器对每个颗粒进行一轮包含一步或多步的轨迹计算。每一步，DPM求解器计算颗粒从当前状态（位置、速度、尺寸、温度等）起在积分时间（即一个颗粒时间步长）内的运动轨迹以及动量、质量和能量损益，并得到更新的颗粒状态。同时，在每一个颗粒时间步喷射一次颗粒。一轮轨迹计算得到的分散相颗粒的动量、质量和能量损益将在下一个连续相迭代步计入连续相源项。积分时间步长和每一轮的步数由用户给定。这样，随着连续相迭代的进行，颗粒将逐轮、逐步地向前推进。
  Particle Time Step Size（颗粒时间步长）：进行一步颗粒轨迹计算时对运动方程作时间积分的起止时长，即每一步轨迹计算中颗粒运动所持续的时间。
  Number of Time Steps（时间步数）：每一轮轨迹计算的步数。
（2）可以。合理么？不知道。
（3）必然的。正常么？不知道。应该何时结束？不知道。
赞一下 回复此楼
9楼2021-12-01 14:39:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cure727

铜虫 (初入文坛)
引用回帖:
9楼: Originally posted by xrs333 at 2021-12-01 14:39:41
（1） 稳态流场＋非稳态颗粒追踪时的颗粒轨迹追踪过程：颗粒轨迹追踪方式为，每隔Number of Continuous Phase Iterations Per DPM Iteration个连续相迭代步数，DPM求解器对每个颗粒进行一轮包含一步或多步的轨迹计算 ...

您好，谢谢您的回复
(1)相关理论看过您发的这个更加了解了一些，因为我前些天在网上看到一片有关时间步长的介绍，里面写到微米级别的颗粒运动时间步长一般为1－10毫秒级别，所以我设的1毫秒，但number of times steps还是不太懂，名义上是说每个颗粒运动时间步长里的迭代次数，是不是应该大一些保证收敛？
(2)懂了，谢谢您的回答
(3)既然合理的话，那我想问一下，我之前稳态流场残差设置的10的负6，必须将残差设置的更小才能继续计算(10的负7)，残差显示加入颗粒之后突然上升，残差曲线很奇怪，而且算的也很慢，几乎半天才动一点点，所以一直在犹豫是不是该继续算下去，或者若是设置的不对就可以及时止损

发自小木虫Android客户端
赞一下 回复此楼
10楼2021-12-01 15:32:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xrs333

木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
10楼: Originally posted by cure727 at 2021-12-01 15:32:11
您好，谢谢您的回复
(1)相关理论看过您发的这个更加了解了一些，因为我前些天在网上看到一片有关时间步长的介绍，里面写到微米级别的颗粒运动时间步长一般为1－10毫秒级别，所以我设的1毫秒，但number of times st ...

求解设置可能对DPM模型的计算支出有影响。你的问题需要考虑萨夫曼升力，虚拟质量力，压力梯度力这几个力吗？他们的量级有多大？
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

cure727
11楼2021-12-01 15:55:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cure727

铜虫 (初入文坛)
送红花一朵
引用回帖:
11楼: Originally posted by xrs333 at 2021-12-01 15:55:54
求解设置可能对DPM模型的计算支出有影响。你的问题需要考虑萨夫曼升力，虚拟质量力，压力梯度力这几个力吗？他们的量级有多大？...

你好，再次感谢您的回复，其实这三个力也是看师兄师姐总选就跟着选上了，具体量级并未进行详细计算，刚刚看了help文档和一些有介绍到粒子受力方面的文献，发现对理论知识还缺乏一定的理解，其中三个公式中都投提到流体速度和颗粒速度，按理说我设的入口速度和颗粒的注射速度是一样的啊，所以按照公式计算岂不都等于零?还请大佬不吝赐教，指导一下小弟

发自小木虫Android客户端
赞一下 回复此楼
12楼2021-12-01 19:55:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xrs333

木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
12楼: Originally posted by cure727 at 2021-12-01 19:55:26
你好，再次感谢您的回复，其实这三个力也是看师兄师姐总选就跟着选上了，具体量级并未进行详细计算，刚刚看了help文档和一些有介绍到粒子受力方面的文献，发现对理论知识还缺乏一定的理解，其中三个公式中都投提到 ...

通常拽力是主要的。你可以比较一下选不选那三个力时的颗粒轨迹，可以在稳态颗粒追踪方式下比较，但比较时要把湍流随机模型关掉。
非稳态颗粒追踪时要注意颗粒包（parcel）数，它太大时计算支出会非常大。
赞一下 回复此楼
13楼2021-12-01 20:32:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cure727

铜虫 (初入文坛)
引用回帖:
13楼: Originally posted by xrs333 at 2021-12-01 20:32:14
通常拽力是主要的。你可以比较一下选不选那三个力时的颗粒轨迹，可以在稳态颗粒追踪方式下比较，但比较时要把湍流随机模型关掉。
非稳态颗粒追踪时要注意颗粒包（parcel）数，它太大时计算支出会非常大。...

您好，谢谢您的回复
我最近大概确定了一下，由于我的颗粒密度远大于流场(空气)密度，所以虚拟质量力暂不考虑，由于入口处的速度较大，并且速度衰减很明显，所以考虑了压力梯度力和萨夫曼升力。
至于您说的注射的parcel数，因为我选择的是面注射，面网格数量为2000多，每次注射2000多个parcel，而我若想模拟2秒的总时间，颗粒注射时间步长暂时打算取0.0秒，追踪步数的确不能取太大，否则parcel数量太过庞大

发自小木虫Android客户端
赞一下 回复此楼
14楼2021-12-06 10:54:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xrs333

木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
14楼: Originally posted by cure727 at 2021-12-06 10:54:41
您好，谢谢您的回复
我最近大概确定了一下，由于我的颗粒密度远大于流场(空气)密度，所以虚拟质量力暂不考虑，由于入口处的速度较大，并且速度衰减很明显，所以考虑了压力梯度力和萨夫曼升力。
至于您说的注射的pa ...

每次2000个parcel是很多的，非稳态追踪parcel数很快就很庞大。建议surface injection不直接使用face zone，而是另外建一个surface，用 bounded surface with sample points方式。这样parcel数才好控制。
赞一下 回复此楼
15楼2021-12-06 18:29:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dpmer

金虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
10楼: Originally posted by cure727 at 2021-12-01 15:32:11
您好，谢谢您的回复
(1)相关理论看过您发的这个更加了解了一些，因为我前些天在网上看到一片有关时间步长的介绍，里面写到微米级别的颗粒运动时间步长一般为1－10毫秒级别，所以我设的1毫秒，但number of times st ...

Number of Time Steps（时间步数）：每一轮轨迹计算的步数。
这个参数还是不太懂，名义上是说每个颗粒运动时间步长里的迭代次数，是不是应该大一些保证收敛？
请问楼主有理解了吗
赞一下 回复此楼
16楼2022-02-23 16:49:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dpmer

金虫 (初入文坛)
引用回帖:
12楼: Originally posted by cure727 at 2021-12-01 19:55:26
你好，再次感谢您的回复，其实这三个力也是看师兄师姐总选就跟着选上了，具体量级并未进行详细计算，刚刚看了help文档和一些有介绍到粒子受力方面的文献，发现对理论知识还缺乏一定的理解，其中三个公式中都投提到 ...

我觉得在流动过程中，两相的速度不一定相等，存在相对速度，这时这些力就不为零了
赞一下 回复此楼
17楼2022-02-23 16:50:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang--pen

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
加入颗粒之后，残差曲线上升很正常。
赞一下 回复此楼
18楼2022-03-07 13:59:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dpmer

金虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
18楼: Originally posted by zhang--pen at 2022-03-07 13:59:19
加入颗粒之后，残差曲线上升很正常。

Number of Time Steps（时间步数）：每一轮轨迹计算的步数。
这个参数还是不太懂，名义上是说每个颗粒运动时间步长里的迭代次数，是不是应该大一些保证收敛？
请问你有什么见解吗
赞一下 回复此楼
19楼2022-03-08 08:39:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
6号枫木塔2楼
2021-12-01 10:00   回复  
cure727(金币+1): 谢谢参与
已获得1个金币 发自小木虫IOS客户端
zhuhuawe3楼
2021-12-01 10:01   回复  
cure727(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
tzynew4楼
2021-12-01 10:02   回复  
cure727(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
卓耿1235楼
2021-12-01 10:03   回复  
cure727(金币+1): 谢谢参与
已获得1个金币 发自小木虫IOS客户端
冰山雪莲灯6楼
2021-12-01 10:06   回复  
cure727(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
相关版块跳转 我要订阅楼主 cure727 的主题更新
191/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 建议基金发布提前给出明确的时间点 +5 kulium 2026-08-21 6/300 2026-08-21 18:15 by alongwaytogo
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] filecode，4个jtjc了 +11 ziyangfang 2026-08-19 14/700 2026-08-21 16:51 by applepetty
[基金申请] 让我中一个面上吧！ +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +15 医学老男孩 2026-08-20 18/900 2026-08-21 13:46 by 医学老男孩
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 只有每年这种时候来逛逛小木虫 +23 yaoyewhu2008 2026-08-20 25/1250 2026-08-21 12:19 by Poppy1104
[基金申请] 时间戳又变了8-15 +15 archvillain 2026-08-15 27/1350 2026-08-21 11:32 by tim76
[基金申请] 科研孤儿太难了 +15 我4大白菜 2026-08-20 16/800 2026-08-21 11:07 by yming1319
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今天系统多次维护，明天很可能放榜！ +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜吗？ +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
信息提示
关闭
请填处理意见
关闭