Physics 板


LINE

不好意思 我的表达能力太差 让大家看不懂我想问的东西...XD 其实我想问的 是下面这问题: =============================== 我们设定的最初浓度 在不同时间的时候 浓度应该也会随时间而衰减吧~~~ (譬如 我在0m处的起点 丢了一把浓度10的粉末 过了1分钟以後 0m处的浓度自然也会有衰减... 所以1分钟後的0m处之浓度 也不会是10了) 请问这个问题 大家是如何处理的 https://imgur.com/TxeHaFC 譬如这个一维浓度扩散方程式而言 是用迭代处理吗 颗颗 ..我打算直接用各种不同的T 手动带进去 ※ 引述《Ecampus (7.7)》之铭言: : : http://imgur.com/TxeHaFC : : 这是一维的污染物浓度扩散方程式 : : u是速度 c是浓度 : : E是扩散系数(延散系数) : : P是衰减系数 : : 我用台湾某条河川 做实地测验、与跑程式模拟两个结果 : : 实地测验得到的曲线 : : 就是高斯分布那一种型态(从0→到浓度最高点→到终点浓度衰退为0) : : 跑程式的话 : : 依照各参数输入 得到的结果却是(0→浓度最高点→到终点仍然浓度最高点 没衰退衰减过) : : → wohtp: 有图吗?你的初始条件和边界条件是什麽? 10/10 10:36 : : → wohtp: 终点是时间上的终点吗?多远? 10/10 10:39 : 您好 初始条件 我是假设开始的释放点(0公尺处) 浓度为100ppb : 边界条件 我假设终点站是N 又C(N-1)=C(N+1) 所以浓度会一直无法下降 : 後来我改成C(N)=0 : (距离是 0M~1600M) : (时间设定是0分钟~160分钟) : 但是 也是到後面快终点的时候 浓度才会降下来 : 我想得到的图形是类似高斯分布这样的 https://imgur.com/w11jVaI : 但我实际得到的结果是https://imgur.com/vZSSWeh : =.= : 我想做的是 时变的系统 浓度随时间变化的 : 不过现在卡关 我知道时不变系统怎麽做....不知道时变系统怎麽做 : 求高手帮忙解惑 薄酬2000P : CC(I)=CC(I)+E*DT*(CC(I+1)-2*CC(I)+CC(I-1))/DX/DX-& : U*DT*(CC(I+1)-CC(I-1))/2./DX-& : DE*DT*CC(I) : CC(N)=0 : 我的时不变系统是这样写 : 如果改成时变..我还在想要如何改 : : → Vulpix: 听起来除了没衰退有点诡异外,一段时间後浓度最高点位置不 10/10 20:20 : : → Vulpix: 变,这件事不是什麽大问题啊。当然视乎你的IC、BC、t。 10/10 20:22 : : → wohtp: 我本来想说是边界条件没设对,扩散到最旁边又被挡回来之类 10/11 10:21 : : → wohtp: 但是发现他还有用手放进去衰减 10/11 10:22 : : → wohtp: 最後没有变成零就是数值做错了,没话好说 10/11 10:23 : : → saltlake: 确定参数输入正确? 跑的程式谁写的? 10/11 11:17 --



※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 125.230.93.173
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Physics/M.1508409657.A.556.html ※ 编辑: Ecampus (125.230.93.173), 10/19/2017 18:44:36
1F:→ wohtp: 这个问题不就是你的微分方程式要描述的吗? 10/19 23:14
2F:→ wohtp: 怎麽处理,就...解方程式啊? 10/19 23:15
3F:→ yeahbo: 简单的说,切空间格点和时间格点,空间格点数值解微方边界 10/21 10:30
4F:→ yeahbo: 值条件,时间格点用叠代 10/21 10:31
5F:→ yeahbo: 你可以参考Finite-difference time-domain method的资料, 10/21 10:34
6F:→ yeahbo: 不过FDTD是处理电磁问题的微方 10/21 10:34
7F:→ bluecadence: 其实原po要的就是用数值分析解偏微分方程 建议他从 10/21 13:24
8F:→ bluecadence: Forward Euler, Backward Euler, Crank-Nicolson 的 10/21 13:25
9F:→ bluecadence: 方法学,运算法中包括楼上y大提到的格子点方法 10/21 13:27







like.gif 您可能会有兴趣的文章
icon.png[问题/行为] 猫晚上进房间会不会有憋尿问题
icon.pngRe: [闲聊] 选了错误的女孩成为魔法少女 XDDDDDDDDDD
icon.png[正妹] 瑞典 一张
icon.png[心得] EMS高领长版毛衣.墨小楼MC1002
icon.png[分享] 丹龙隔热纸GE55+33+22
icon.png[问题] 清洗洗衣机
icon.png[寻物] 窗台下的空间
icon.png[闲聊] 双极の女神1 木魔爵
icon.png[售车] 新竹 1997 march 1297cc 白色 四门
icon.png[讨论] 能从照片感受到摄影者心情吗
icon.png[狂贺] 贺贺贺贺 贺!岛村卯月!总选举NO.1
icon.png[难过] 羡慕白皮肤的女生
icon.png阅读文章
icon.png[黑特]
icon.png[问题] SBK S1安装於安全帽位置
icon.png[分享] 旧woo100绝版开箱!!
icon.pngRe: [无言] 关於小包卫生纸
icon.png[开箱] E5-2683V3 RX480Strix 快睿C1 简单测试
icon.png[心得] 苍の海贼龙 地狱 执行者16PT
icon.png[售车] 1999年Virage iO 1.8EXi
icon.png[心得] 挑战33 LV10 狮子座pt solo
icon.png[闲聊] 手把手教你不被桶之新手主购教学
icon.png[分享] Civic Type R 量产版官方照无预警流出
icon.png[售车] Golf 4 2.0 银色 自排
icon.png[出售] Graco提篮汽座(有底座)2000元诚可议
icon.png[问题] 请问补牙材质掉了还能再补吗?(台中半年内
icon.png[问题] 44th 单曲 生写竟然都给重复的啊啊!
icon.png[心得] 华南红卡/icash 核卡
icon.png[问题] 拔牙矫正这样正常吗
icon.png[赠送] 老莫高业 初业 102年版
icon.png[情报] 三大行动支付 本季掀战火
icon.png[宝宝] 博客来Amos水蜡笔5/1特价五折
icon.pngRe: [心得] 新鲜人一些面试分享
icon.png[心得] 苍の海贼龙 地狱 麒麟25PT
icon.pngRe: [闲聊] (君の名は。雷慎入) 君名二创漫画翻译
icon.pngRe: [闲聊] OGN中场影片:失踪人口局 (英文字幕)
icon.png[问题] 台湾大哥大4G讯号差
icon.png[出售] [全国]全新千寻侘草LED灯, 水草

请输入看板名称,例如:BuyTogether站内搜寻

TOP