24小时热门版块排行榜    

查看: 2274  |  回复: 11

winterwac

铁虫 (小有名气)

[求助] Helmholtz方程仿真光传输问题(边界条件设置,3D简化2D) 已有2人参与

初学者,试图用Helmholtz仿真光在生物组织内的传输问题,最后观察变量为组织表面的u(fluence rate),通过查阅文献,查到了可以使用Helmholtz作为控制方程,Robin boundary condition作为边界条件,但根据文献所述的设置利用2D axisymmetric模型得出的结果与Monte Carlo得出的结果相差甚远,然后又试着用3D模型进行仿真,在使用点光源(边界内部一个自由光程处设置点光源),归一化后的结果和Monte Carlo结果较为接近,但是这样光源设置就无法模拟真实条件下的光源了。而且3D模型仿真耗时较长,理论上3D和2D axisymmetric结果应该一样,但是得到的结果却很不同,如下图,同样的附上mph文件和参考文献,希望懂得人或者相关领域的能一起交流,谢谢!Helmholtz方程仿真光传输问题(边界条件设置,3D简化2D)
Screen Shot 2015-03-31 at 10.38.52 AM.png
回复此楼

» 本帖附件资源列表

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

chentao807

铜虫 (小有名气)

引用回帖:
3楼: Originally posted by winterwac at 2015-04-01 20:59:43
结构就是一个简单的圆柱体,首先肯定是对称的,边界条件也很简单,就是所有的边界都一样的,考虑了一个边界反射,设置也都是一样的,但是结果却就是有很大的不同。不知您方便帮忙看下?非常感谢...

我的版本太旧打不开你的文件。我觉得你对称的那个界面上边界条件应该是错了。
4楼2015-04-01 21:47:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

chentao807

铜虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
2维对称和3维一样嘛?这可未必,首先你的结构必须有对称性,可以简化为对称结构。第二点更重要的是合理的施加边界条件。
如果两条都满足的话,可以负责任的说,结果基本是一样的。(对你研究的Helmholtz方程,应该不会例外,这个方程并不复杂)。
2楼2015-04-01 13:02:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

winterwac

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by chentao807 at 2015-04-01 13:02:24
2维对称和3维一样嘛?这可未必,首先你的结构必须有对称性,可以简化为对称结构。第二点更重要的是合理的施加边界条件。
如果两条都满足的话,可以负责任的说,结果基本是一样的。(对你研究的Helmholtz方程,应该 ...

结构就是一个简单的圆柱体,首先肯定是对称的,边界条件也很简单,就是所有的边界都一样的,考虑了一个边界反射,设置也都是一样的,但是结果却就是有很大的不同。不知您方便帮忙看下?非常感谢
3楼2015-04-01 20:59:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

winterwac

铁虫 (小有名气)

引用回帖:
4楼: Originally posted by chentao807 at 2015-04-01 21:47:21
我的版本太旧打不开你的文件。我觉得你对称的那个界面上边界条件应该是错了。...

是指对称轴上的么?2D中除了对称轴上没有设置边界条件,其他都是和3D中的边界条件一样的。
5楼2015-04-03 21:01:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

winterwac

铁虫 (小有名气)

引用回帖:
4楼: Originally posted by chentao807 at 2015-04-01 21:47:21
我的版本太旧打不开你的文件。我觉得你对称的那个界面上边界条件应该是错了。...

能否留个QQ或者邮箱进一步交流?我这个问题卡很久了。。。
6楼2015-04-03 21:10:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Dustin_q

铜虫 (初入文坛)

我也正在做光扩散的仿真,资料不太好找。你的这些工作对我很有帮助。谢谢楼主啦~
7楼2017-02-15 10:35:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Dustin_q

铜虫 (初入文坛)

【答案】应助回帖

兄弟,你仔细看看你3D的模型,你圆柱体之外的点源似乎没有用到,我认为你得到的结果是由于设置了内部的源项。2D的应该是正确的。
8楼2017-03-02 16:09:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Dustin_q

铜虫 (初入文坛)

【答案】应助回帖

兄弟,你仔细看看你3D的模型,你圆柱体之外的点源似乎没有用到,我认为你得到的结果是由于设置了内部的源项。2D的应该是正确的。
9楼2017-03-02 16:10:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Dustin_q

铜虫 (初入文坛)

【答案】应助回帖

兄弟,你仔细看看你3D的模型,你圆柱体之外的点源似乎没有用到,我认为你得到的结果是由于设置了内部的源项。2D的应该是正确的。
10楼2017-03-02 16:10:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 winterwac 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 建议基金发布提前给出明确的时间点 +11 kulium 2026-08-21 14/700 2026-08-24 02:12 by 流流伤
[基金申请] 2026国自然函评费到账 +16 羊腰板 2026-08-21 17/850 2026-08-23 23:02 by tianxiaochun
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +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 紫杉醇
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见