作者Lanjaja ()
看板Physics
标题[问题] 位势流的矛盾
时间Fri Dec 22 13:20:13 2017
【出处】(习题或问题的出处)
阅读流体力学课本後所产生的疑问。
【题目】(题目的文字叙述,如有图片亦可提供图片)
1.考虑流体力学中满足▽‧v=0的2维流体,可以旋度表示,即stream function,
这个地方没问题,
但是课文又提到如果进一步又加了无黏滞性的条件 => (▽^2)v = 0
利用向量恒等式可得LHS = ▽(▽‧v) - ▽ x (▽ x v) = RHS = 0
=>▽ x (▽ x v) = 0
=> 然後说所以▽ x v = 0。
一直想不清楚为什麽最後可以从▽x(▽xv) = 0 => ▽ x v = 0
2.另外一个相关问题问题
也有书本直接假设考虑流体力学中满足▽‧v=0的2维流体且无旋性▽ x v = 0
我的疑问是根据亥姆霍兹定理,一个向量v可以表示成-▽φ + ▽ x A
如果已经满足▽‧v=0,v将没有-▽φ的部分,
按照定理v = ▽ x A
但是因为无旋性=> v 又可表示成 -▽φ,
这不就上面说的满足▽‧v=0就没有-▽φ部份矛盾
这是否又表示无旋性和无散度不能够同时满足?否则v就只能是nontrivial的向量场?
【瓶颈】(解题瓶颈或思考脉络,请尽量详述以利回答者知道要从何处讲解指导)
(正确示范:我算出来的答案好像不太对,这是我的计算过程,哪里出问题?)
我已经将我的思路叙述在题目中,
恳请强者帮忙回答,
感谢!
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 61.56.10.112
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Physics/M.1513920015.A.91F.html
1F:推 Diluvius: 1.随便举都有反例 怀疑课本有其他前提 因为通常无黏滞不 12/22 15:20
2F:→ Diluvius: 会写Laplacian=0 应该有假设某个无因次量大小 12/22 15:22
纳维尔-史托克方程黏滞力项是ν▽^2 v
所以无黏滞力=> ▽^2 v=0
3F:→ Diluvius: 2.potential不唯一 12/22 15:23
4F:推 Diluvius: 1.通常令系数=0 现在想得到这样令的是Re#<<1 还是你说说 12/22 16:32
我知道你说的是化成无因次方程式,黏滞力项系数~1/Re#
书本是说ν▽^2 v = 0,不是ν=0造成的,而是▽^2 v = 0造成的,我也觉得好怪
书本进一步说▽^2 v = 0 = ▽(▽‧v) - ▽ x (▽ x v) = 0 - ▽ x (▽ x v)
=> ▽ x (▽ x v) = 0 及 ▽ x v = 0
^^^^^^^^^^
然後书本就引入位势流,
我想知道▽ x (▽ x v) = 0是靠什麽条件得到curl(v)=0?
这是一本叫做云物理的书,我怀疑作者不懂随便乱写。
其他流体力学如Fox是直接假设▽‧v = 0和人为限制非旋性流场的条件,
并没有讲到为何是非旋性
5F:→ Diluvius: 这是哪本书哪里 我找找 12/22 16:32
6F:→ Diluvius: 2.对啊 你给了两个potential啊? 12/22 16:33
我可能叙述没有很清楚,我的意思是
从v = -▽φ + ▽ x A为出发点,
因为▽‧v = 0 => v = ▽ x A,等同-▽φ=0
另一方面,无旋性流体的条件▽ x v = 0= > v = -▽φ,等同▽ x A = 0
^^^^^^^^因为前已满足▽‧v = 0
=>▽‧v = 0=>v = -▽φ = 0
所以v = 0?
两者不就互相抵触?
除非v = 0或常向量trivial的情况?
※ 编辑: Lanjaja (61.56.10.112), 12/22/2017 17:27:38
7F:推 KBmax: 你是不是搞错什麽了@@ 一个向量场的divergence等於零 不保 12/22 18:14
8F:→ KBmax: 证gradient potential等於零啊,是potential的laplacian等 12/22 18:14
请问K大如果▽‧v = 0的话,
如何保证一定存在一个gradient potentialφ使得v = -▽φ?
从v = -▽φ + ▽ x A为出发点,并没有保证存在v = -▽φ
如何能说有个东西叫gradient potential?
只有在v可以表示成-▽φ才称φ为gradient potential
我只知道▽‧v = 0 => 存在vector potential A使得v = ▽ x A
9F:→ KBmax: 於零,不然静电学一堆静电场都不存在了== 12/22 18:14
静电学是▽ x E = 0 => E = -▽V 是curl,不是divergence
KB是不是误会了
10F:推 KBmax: 不要乱背课本结论啊,自己把算子带进去算算看就知道自己该 12/22 18:23
11F:→ KBmax: 要注意什麽了,这些东西数学上的本质都一样 12/22 18:23
12F:→ silverywings: 如果引入 potential flow 那速度场就是一位能函数的 12/22 20:51
13F:→ silverywings: 梯度,一纯量取梯度其旋度一定会等於 0 12/22 20:51
没错阿,整篇脉络都是以
▽ x v = 0 <=> v = -▽φ
▽‧v = 0 <=> v = ▽ x A
为依归,并没有KB大指的搞错
14F:→ Lanjaja: 我的问题是从亥姆霍兹定理如果说明一个nontrivial向量场 12/22 22:01
15F:→ Lanjaja: 可以同时满足divergence=0和curl=0的条件,能够同时满足 12/22 22:01
16F:→ Lanjaja: 又不为trivial形式吗? 12/22 22:01
※ 编辑: Lanjaja (118.165.144.212), 12/22/2017 22:08:15
17F:→ Vulpix: ▽ x v = 0 <= v = -▽φ 是对的。但反过来不见得对,毕竟 12/22 22:37
18F:→ Vulpix: Helmholtz Thm. 需要 v 随着距离拉远而快速衰减才行。 12/22 22:38
19F:→ Vulpix: 会出现的差异项就是 harmonic field。所以其实应该说 12/22 22:44
20F:→ Vulpix: Helmholtz Thm. 并非一个完整的分解方式,他只能分解 12/22 22:45
21F:→ Vulpix: 1. v=0 on boundary (空间有界的时候) 2. v→0 "rapidly" 12/22 22:47
22F:→ Vulpix: for "far-away points" (空间各向无界) 3. 其实还有其他 12/22 22:49
23F:→ Vulpix: 空间的长相,不过用到的时候再谈就好。 12/22 22:49
24F:→ Vulpix: 但不管是哪一种状况都是为了让 harmonic term 消失。 12/22 22:53
谢谢V大对於亥姆霍兹定理的注解
其实我在一般书上没有看过真正严谨的限制和假设,所以用在什麽情况也不太确定。
但是文中说的同时满足▽ x v = 0 = ▽‧v的流场确实有v = (K/r) e_θ的例子,
就是您在推文中所说的harmonic field
所以我才会想把这个定理用来解释无散且无旋流体的问题上,因而产生这篇疑问。
25F:推 Diluvius: 虽然没讲作者 不过有的书应该都不是那本 暂时无法确定 12/22 23:07
Fox麦当劳流体力学课本和云物理学张宝贯着
26F:→ Diluvius: v=v1+v2 v1=无旋流 v2=无散流 要不v1=v v2=0 或者反过来 12/22 23:09
27F:→ Diluvius: 或者v1=v2=v/2 就有各种potential 12/22 23:10
没有错,我也曾经得到v1=v2(= v,这是流体力学课本的内容)这种奇怪结论,
可是如果v1=v2,除了常向量场之外,有可能还有其他场吗?
就像我上面的推论:
从v = -▽φ + ▽ x A --(1)出发,假设v满足前面V大推文中说的假设条件,
因为不可压缩性:▽‧v = 0,(1)代入 0 = -▽^2 φ
但是我又知道v = ▽ x A --(2),比较(1)後得-▽φ = 0 --(3)
v = ▽ x A,但是将(1)代入▽‧v = 0得-▽^2 φ = 0
另一方面,又因无旋性流体的条件▽ x v = 0,(1)代入 ▽ x(▽ x A) = 0
又因为▽ x v = 0,v可表示成-▽φ --(4),由(3)得v = -▽φ = 0 by(3)
可是v = c似乎也可以,只是仿上面的过程我得不出来,更别说其他例子,
可以请板上先近帮忙看一下推论的盲点,谢谢
※ 编辑: Lanjaja (118.165.144.212), 12/22/2017 23:53:58
28F:推 Diluvius: nontrivial的情况? 请看应数正交轨迹或柯西-黎曼方程 12/22 23:23
29F:推 ntust661: 单纯看你的文章,第一点是不成立的。▽ x v = 0 只是其 12/27 22:32
30F:→ ntust661: 中一个解而已。 12/27 22:32
31F:→ ntust661: 我只知道这些向量运算跟velocity potential一点关系都没 12/27 22:34
32F:→ ntust661: 有XD 12/27 22:34
33F:→ ntust661: stream function的定义没那麽麻烦~很简单的 12/27 22:36