☆─────────────────────────────────────☆
laser2000 (掌上智能版都是大厦挂) 于 (Tue Aug 12 00:45:49 2014) 提到:
我知道有无粘线性稳定性，层流线性稳定性
但是对于不存在稳定解的问题，怎么选择basic flow?
某些流动有很强的周期变化，选择任何状态作为basic flow都会与小扰动矛盾，怎么
搞??
※ 来源:·水木社区
http://newsmth.net·[FROM: 61.150.43.*]
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Tue Aug 12 10:41:55 2014) 提到:
比如剪切射流？
设个假的理想解……
然后稳定性能算出来那个周期变化的尺度……
什么？你说湍流里面的相干结构？ 我建议你选择平均流……
理论上大尺度的失稳结构即便流动是非稳定的 是可以用平均解解析这个非稳结构的
因为这个时候的增长率都普遍较大
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: 我知道有无粘线性稳定性，层流线性稳定性
: 但是对于不存在稳定解的问题，怎么选择basic flow?
: 某些流动有很强的周期变化，选择任何状态作为basic flow都会与小扰动矛盾，怎么
: ...................
☆─────────────────────────────────────☆
LaPi (坎切拉皮) 于 (Tue Aug 12 11:29:17 2014) 提到:
我猜如果周期比较长的话
可不可以每个时刻的线性稳定性都算一下？采用该时刻的瞬时流动做平均流
如果周期短的话这个方法应该是不太行的
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: 我知道有无粘线性稳定性，层流线性稳定性
: 但是对于不存在稳定解的问题，怎么选择basic flow?
: 某些流动有很强的周期变化，选择任何状态作为basic flow都会与小扰动矛盾，怎么
: ...................
☆─────────────────────────────────────☆
laser2000 (掌上智能版都是大厦挂) 于 (Tue Aug 12 15:24:52 2014) 提到:
3x
那有些复杂流动，频率成分很多，没法确定明显的周期，那在足够长的时间内取平均流可
以么？
此外，如果采用DES之类来计算平均流，可否用来做线性稳定性分析呢？
【 在 shookware (@viocy) 的大作中提到: 】
: 比如剪切射流？
: 设个假的理想解……
: 然后稳定性能算出来那个周期变化的尺度……
: ...................
☆─────────────────────────────────────☆
LaPi (坎切拉皮) 于 (Tue Aug 12 15:55:12 2014) 提到:
圡问一个：
都已经是湍流了，再做线性稳定性分析还有意义不。。。。。
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: 3x
: 那有些复杂流动，频率成分很多，没法确定明显的周期，那在足够长的时间内取平均流可
: 以么？
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Tue Aug 12 16:11:02 2014) 提到:
大尺度的也许可以
DES算的平均流可以不可以不知道 我倾向于需要看情况
原则上只要失稳的时间尺度空间尺度相对于基本流足够小 就可以分析
虽然不一定准 但是不会差太多
不过可能做湍流中的失稳模态分析需要做global instability的分析
这样对于不搞这个方向的来说 会有一些难度 （这些难度本身并不是由于GI的分析工具
有多难写 而是对GI的物理本质的认识）
换句话说 看你分析的对象到底是什么了 如果是尾迹流或者是射流 我倾向于认为可以用
平均流+线性稳定性分析试试 看看结果怎么样…… 我知道有人这么搞的 而且理论上由
于这种剪切流动的增长率很大 所以不会差太多。。
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: 3x
: 那有些复杂流动，频率成分很多，没法确定明显的周期，那在足够长的时间内取平均流可
: 以么？
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Tue Aug 12 16:11:46 2014) 提到:
不过如果你用来做壁面湍流分析 我觉得这问题就需要仔细考察了
恐怕会比较难…… 虽然我觉得这题目还是很有趣就是了。。。
呵呵
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: 3x
: 那有些复杂流动，频率成分很多，没法确定明显的周期，那在足够长的时间内取平均流可
: 以么？
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Tue Aug 12 16:14:53 2014) 提到:
可以
stocks流动就是这种流动的一个例子
既有global的稳定性问题 又有一个local的稳定性问题
总之复杂又有趣就是了
【 在 LaPi (坎切拉皮) 的大作中提到: 】
: 我猜如果周期比较长的话
: 可不可以每个时刻的线性稳定性都算一下？采用该时刻的瞬时流动做平均流
: 如果周期短的话这个方法应该是不太行的
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Tue Aug 12 16:16:21 2014) 提到:
有
我没看到 但是据说最近有比较振奋的消息是 用线性稳定性算出来的湍流的相干结构与
实验观测到的尺度相似……
另外 T大的许春晓老师的办公室不是也用这个来分析湍流么……
【 在 LaPi (坎切拉皮) 的大作中提到: 】
: 圡问一个：
: 都已经是湍流了，再做线性稳定性分析还有意义不。。。。。
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Tue Aug 12 16:17:06 2014) 提到:
当然 做当地稳定性分析的前提是 基本流动的变化尺度远远大于失稳扰动的尺度……
【 在 LaPi (坎切拉皮) 的大作中提到: 】
: 我猜如果周期比较长的话
: 可不可以每个时刻的线性稳定性都算一下？采用该时刻的瞬时流动做平均流
: 如果周期短的话这个方法应该是不太行的
: ...................
☆─────────────────────────────────────☆
laser2000 (掌上智能版都是大厦挂) 于 (Tue Aug 12 16:17:57 2014) 提到:
……
我一个搞CFD算法的，赶鸭子上架做流体物理分析，确实是对稳定性分析的物理理解不熟
悉
我要搞的是高Re数剪切层，非自由的，受到几何空间干扰的那类
我的理解是，希望用稳定性分析方法得到不稳定模态及其频率、增长率等性质，因为这些
东西只从模拟是没法得到的。然后希望用这些性质通过参数辨识来搞个低维的模型。
【 在 shookware (@viocy) 的大作中提到: 】
: 大尺度的也许可以
: DES算的平均流可以不可以不知道 我倾向于需要看情况
: 原则上只要失稳的时间尺度空间尺度相对于基本流足够小 就可以分析
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Tue Aug 12 16:22:12 2014) 提到:
非自由剪切层的话 我个人认为只能用稳定性算算看了。。
不知道壁面边界对它的影响到底有多大
如果有兴趣并且可能的话 倒是可以合作做这个看看 我对这玩意倒是确实挺感兴趣的 但
是苦于不大会算湍流 囧~~~（ 手残网格画起来是个二把刀 湍流模式啥的也没真正的做
这一块 DNS啥的虽然算但是仅仅是验证理论模型的工具这么用）
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: ……
: 我一个搞CFD算法的，赶鸭子上架做流体物理分析，确实是对稳定性分析的物理理解不熟
: 悉
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Tue Aug 12 16:22:53 2014) 提到:
对了
你这个被分析对象是不可压缩流动还是可压缩流动？
如果可压缩流动的话 Ma大概多少？
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: ……
: 我一个搞CFD算法的，赶鸭子上架做流体物理分析，确实是对稳定性分析的物理理解不熟
: 悉
: ...................
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Tue Aug 12 17:08:52 2014) 提到:
Floquet Operator + Arnoldi method
在周期解的情况下可以做进一步的稳定性分析
比如Von Karman flow 向三维A模态和B模态的过渡
参看 Barkley & Henderson JFM 1995
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: ……
: 我一个搞CFD算法的，赶鸭子上架做流体物理分析，确实是对稳定性分析的物理理解不熟
: 悉
: ...................
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Tue Aug 12 17:10:37 2014) 提到:
长周期就是WKBJ的方法了
【 在 LaPi (坎切拉皮) 的大作中提到: 】
: 我猜如果周期比较长的话
: 可不可以每个时刻的线性稳定性都算一下？采用该时刻的瞬时流动做平均流
: 如果周期短的话这个方法应该是不太行的
: ...................
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Tue Aug 12 17:15:16 2014) 提到:
好多人都用时间平均的流场作为基解
个人认为不对
可以考虑使用selective frequency damping 做基解
其实这个方法也是有些问题的 不过比时间平均的流场好一点 毕竟可以找到驻点解
再就是用Newton迭代法直接找驻点解，作为基解，这个我个人认为比较靠普点
比如Von Karman流动 Re > 47， 可以用Newton迭代法做continuation 找到一个不稳定的驻点解，然后以此为基解，做稳定性分析
【 在 shookware (@viocy) 的大作中提到: 】
: 比如剪切射流？
: 设个假的理想解……
: 然后稳定性能算出来那个周期变化的尺度……
: ...................
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Tue Aug 12 17:15:40 2014) 提到:
re
【 在 LaPi (坎切拉皮) 的大作中提到: 】
: 圡问一个：
: 都已经是湍流了，再做线性稳定性分析还有意义不。。。。。
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Tue Aug 12 17:17:40 2014) 提到:
【 在 shookware (@viocy) 的大作中提到: 】
: 有
: 我没看到 但是据说最近有比较振奋的消息是 用线性稳定性算出来的湍流的相干结构与
: 实验观测到的尺度相似……
这个RB对流里好多年了，但是问题是物理机制是否真的一样
最后还是POD和DMD了
: ...................
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Tue Aug 12 17:19:47 2014) 提到:
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: ……
: 我一个搞CFD算法的，赶鸭子上架做流体物理分析，确实是对稳定性分析的物理理解不熟
: 悉
: 我要搞的是高Re数剪切层，非自由的，受到几何空间干扰的那类
高Re ，识别频率和增长率
可以考虑使用DMD
不过Koopman Operator的数学背景我还没有完全理解
不保证完全可行
: 我的理解是，希望用稳定性分析方法得到不稳定模态及其频率、增长率等性质，因为这些
: 东西只从模拟是没法得到的。然后希望用这些性质通过参数辨识来搞个低维的模型。
☆─────────────────────────────────────☆
laser2000 (掌上智能版都是大厦挂) 于 (Tue Aug 12 17:29:31 2014) 提到:
你是说用DMD得到离散系统
x[n+1] = A x[n]
的系统矩阵A的降阶形式，再做特征值分析么？
【 在 yihongyuan (怡红之主) 的大作中提到: 】
: 高Re ，识别频率和增长率
: 可以考虑使用DMD
: 不过Koopman Operator的数学背景我还没有完全理解
: ...................
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Tue Aug 12 17:37:19 2014) 提到:
对的
看Peter Schmid 2010 JFM就好了 网上有
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: 你是说用DMD得到离散系统
: x[n+1] = A x[n]
: 的系统矩阵A的降阶形式，再做特征值分析么？
: ...................
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Tue Aug 12 17:38:04 2014) 提到:
至于频率和增长率，就复杂一点了
看Rowley等人的文章吧
【 在 yihongyuan (怡红之主) 的大作中提到: 】
: 对的
: 看Peter Schmid 2010 JFM就好了 网上有
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Tue Aug 12 20:05:38 2014) 提到:
他这个问题不是周期边界的
所以不能用floquet展开。。
【 在 yihongyuan (怡红之主) 的大作中提到: 】
: Floquet Operator + Arnoldi method
: 在周期解的情况下可以做进一步的稳定性分析
: 比如Von Karman flow 向三维A模态和B模态的过渡
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Tue Aug 12 20:07:46 2014) 提到:
你这说的都是纯理论模型问题。。
他这个是个实际的可压缩流动 对特定频率做滤波的特定频率没办法选的。。
算出来什么就是什么的那种。。
【 在 yihongyuan (怡红之主) 的大作中提到: 】
: 好多人都用时间平均的流场作为基解
: 个人认为不对
: 可以考虑使用selective frequency damping 做基解
: ...................
☆─────────────────────────────────────☆
LaPi (坎切拉皮) 于 (Tue Aug 12 21:56:52 2014) 提到:
哦哦哦，算出来不稳定，然后可以用最大增长率或者其他成分分析一下如何演化？
或者说线性假设下能演化成啥样子？
如果用线性稳定性算出来的和实验观测到有相似的东西
是不是也只是某些特例，不具有推广性？
【 在 shookware (@viocy) 的大作中提到: 】
: 有
: 我没看到 但是据说最近有比较振奋的消息是 用线性稳定性算出来的湍流的相干结构与
: 实验观测到的尺度相似……
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Tue Aug 12 23:29:47 2014) 提到:
结论是 不知道……
【 在 LaPi (坎切拉皮) 的大作中提到: 】
: 哦哦哦，算出来不稳定，然后可以用最大增长率或者其他成分分析一下如何演化？
: 或者说线性假设下能演化成啥样子？
: 如果用线性稳定性算出来的和实验观测到有相似的东西
: ...................
☆─────────────────────────────────────☆
wealth (阳光普照) 于 (Wed Aug 13 00:01:20 2014) 提到:
for the stability of a periodic flow, Floquet analysis is needed.
perturbation can be written as to be proportional to exp(sigma*t+n*f*t) where f is the frequency of base flow, n is an integer and sigma is the so-called Floquet frequency (the eigenvalue of the linear problem).
【 在 laser2000 的大作中提到: 】
: 我知道有无粘线性稳定性，层流线性稳定性
: 但是对于不存在稳定解的问题，怎么选择basic flow?
: 某些流动有很强的周期变化，选择任何状态作为basic flow都会与小扰动矛盾，怎么
: ...................
☆─────────────────────────────────────☆
rich1983 (红叶疯了) 于 (Wed Aug 13 10:38:18 2014) 提到:
请教。看到过一篇文章是在流动方向取若干截面做local 稳定性，得到每个截面的增长率之类的，这种做法能理解为global稳定性吗？
【 在 shookware 的大作中提到: 】
: 可以
: stocks流动就是这种流动的一个例子
: 既有global的稳定性问题 又有一个local的稳定性问题
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Wed Aug 13 10:42:38 2014) 提到:
不能
【 在 rich1983 (红叶疯了) 的大作中提到: 】
: 请教。看到过一篇文章是在流动方向取若干截面做local 稳定性，得到每个截面的增长率之类的，这种做法能理解为global稳定性吗？
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Wed Aug 13 10:43:34 2014) 提到:
全局稳定性是说 在流动里面全局的失稳扰动都有同样的一个增长率
要增长一起涨 要衰减一起衰
【 在 rich1983 (红叶疯了) 的大作中提到: 】
: 请教。看到过一篇文章是在流动方向取若干截面做local 稳定性，得到每个截面的增长率之类的，这种做法能理解为global稳定性吗？
☆─────────────────────────────────────☆
rich1983 (红叶疯了) 于 (Wed Aug 13 11:06:02 2014) 提到:
什么方法做？比如每个截面速度型等都不一样，那基本流怎么给？
【 在 shookware 的大作中提到: 】
: 全局稳定性是说 在流动里面全局的失稳扰动都有同样的一个增长率
: 要增长一起涨 要衰减一起衰
:
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Wed Aug 13 11:25:31 2014) 提到:
整个流动都是基本流啊~
直接L*v=-i*omega*v 这里面omega就是特征值啊~
L就是那个流场的控制方程的空间算子
【 在 rich1983 (红叶疯了) 的大作中提到: 】
: 什么方法做？比如每个截面速度型等都不一样，那基本流怎么给？
☆─────────────────────────────────────☆
laser2000 (掌上智能版都是大厦挂) 于 (Wed Aug 13 22:52:30 2014) 提到:
global (in)stability 不是考虑三维扰动情况下的（不）稳定性么？
【 在 shookware (@viocy) 的大作中提到: 】
: 全局稳定性是说 在流动里面全局的失稳扰动都有同样的一个增长率
: 要增长一起涨 要衰减一起衰
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Wed Aug 13 22:54:46 2014) 提到:
主要是基本流动的算子不是只与法向有关的。。
或者说没有办法在其它方向上引入波假设。。
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: global (in)stability 不是考虑三维扰动情况下的（不）稳定性么？
☆─────────────────────────────────────☆
laser2000 (掌上智能版都是大厦挂) 于 (Wed Aug 13 22:59:44 2014) 提到:
3x
最近看了一些文献，有许多不解，请各位指点：
1 对于多个频率成分，且并非倍频关系，Floquet方法没法搞了吧
2 关于平均流作为基本流场，似乎看得到不同看法，有的说可以，有的说不好。
3 见到文献用DMD方法建立线性讲解模型，再做EVP分析，效果不错。但是前提是，流动快照处于小扰动线性区，如果流动本身饱和了，就不能这么搞。
【 在 wealth (阳光普照) 的大作中提到: 】
: for the stability of a periodic flow, Floquet analysis is needed.
: perturbation can be written as to be proportional to exp(sigma*t+n*f*t) where f is the frequency of base flow, n is an integer and sigma is the so-called Floquet frequency (the eigenvalue of the linear problem).
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Wed Aug 13 23:04:02 2014) 提到:
流动饱和的话 得引入非线性模型了 那这问题就复杂多了。。
但是我觉得尺度上用线性模型来考虑应该不会差的太严重
但是增长率必然差到姥姥家了。。
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: 3x
: 最近看了一些文献，有许多不解，请各位指点：
: 1 对于多个频率成分，且并非倍频关系，Floquet方法没法搞了吧
: ...................
☆─────────────────────────────────────☆
wealth (阳光普照) 于 (Thu Aug 14 00:01:02 2014) 提到:
: 1 对于多个频率成分，且并非倍频关系，Floquet方法没法搞了吧
my guess: maybe perturbation can be written to be proportional to exp(sigma*t+n1*f1*t+n2*f2*t+...). You need to check if this is correct. Even if this works, the matrix is huge. So it is better to use time-stepping rather than eigenvalue calculation.
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Thu Aug 14 06:33:24 2014) 提到:
这办法不行我猜……
你要是同时有多个频率在里面 这俩不同的频率之间还有非线性作用呢。。
那恐怕至少还要再去用伴随方程去解弱非线性才行
【 在 wealth (阳光普照) 的大作中提到: 】
: my guess: maybe perturbation can be written to be proportional to exp(sigma*t+n1*f1*t+n2*f2*t+...). You need to check if this is correct. Even if this works, the matrix is huge. So it is better to use time-stepping rather than eigenvalue calculation.
☆─────────────────────────────────────☆
wealth (阳光普照) 于 (Thu Aug 14 08:25:32 2014) 提到:
this linear stability analysis will give rise to the triad resonance of modes of base flow and modes of instability but cannot give rise to the interaction of modes of base flow.
【 在 shookware 的大作中提到: 】
: 这办法不行我猜……
: 你要是同时有多个频率在里面 这俩不同的频率之间还有非线性作用呢。。
: 那恐怕至少还要再去用伴随方程去解弱非线性才行
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Thu Aug 14 09:43:41 2014) 提到:
你那个主流里面的周期解一样也可以认为是一种具有一定尺度的扰动的。能做的办法是如果认为是个线性分析，那么做全局稳定性分析就好了。或者失稳扰动尺度远小于基本流扰动尺度的话用局部解做。但是我觉得这问题没这简单，得仔细做再验证才知道。而且也用不着对基本流做展开直接对算子求特征值就行了。做展开的目的就是矩阵规模小一些。而且他这个问题，基本流应该也不是真的纯周期的的，实际上只是用des算出来的一个湍流场。
【 在 wealth 的大作中提到: 】
: this linear stability analysis will give rise to the triad resonance of modes of base flow and modes of instability but cannot give rise to the interaction of modes of base flow.
:
☆─────────────────────────────────────☆
wealth (阳光普照) 于 (Thu Aug 14 09:52:40 2014) 提到:
Yes, you are right. Linear (local and global) stability analysis is in terms of laminar flow. If the base flow is already unstable or even turbulent then it is better to do time-stepping calculation.
【 在 shookware 的大作中提到: 】
: 你那个主流里面的周期解一样也可以认为是一种具有一定尺度的扰动的。能做的办法是如果认为是个线性分析，那么做全局稳定性分析就好了。或者失稳扰动尺度远小于基本流扰动尺度的话用局部解做。但是我觉得这问题没这简单，得仔细做再验证才知道。而且也用不着对基本流做展开直接对算子求特征值就行了。做展开的目的就是矩阵规模小一些。而且他这个问题，基本流应该也不是真的纯周期的的，实际上只是用des算出来的一个湍流场。
: :
☆─────────────────────────────────────☆
LaPi (坎切拉皮) 于 (Thu Aug 14 11:24:24 2014) 提到:
湍流应该是一种线性小扰动充分发展以后的现象了。
如果说基本流是怎么样，从道理上讲对应的是流动变成湍流之前的层流，
或者说失稳之前的流动。
用这个时刻的流动去分析该状态对应的稳定性情况没问题
但是拿这个来分析充分发展的非线性流动，显然有问题。
而拿湍流本身当做基本流直接做线性分析，这个根基上一样有点不稳的感觉。。。
【 在 shookware (@viocy) 的大作中提到: 】
: 你那个主流里面的周期解一样也可以认为是一种具有一定尺度的扰动的。能做的办法是如果认为是个线性分析，那么做全局稳定性分析就好了。或者失稳扰动尺度远小于基本流扰动尺度的话用局部解做。但是我觉得这问题没这简单，得仔细做再验证才知道。而且也用不着对基本流做展
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Thu Aug 14 11:32:34 2014) 提到:
Floquet要求基解时间上是周期的就好啦
空间边界条件可以任意
【 在 shookware (@viocy) 的大作中提到: 】
: 他这个问题不是周期边界的
: 所以不能用floquet展开。。
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Thu Aug 14 11:38:01 2014) 提到:
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: 3x
: 最近看了一些文献，有许多不解，请各位指点：
: 1 对于多个频率成分，且并非倍频关系，Floquet方法没法搞了吧
多频搞不了了
: 2 关于平均流作为基本流场，似乎看得到不同看法，有的说可以，有的说不好。
如果是避免流动，比如channel flow
Gibson的group 基本上是用Newton法得到的驻点解 fixed point 作为基解
: 3 见到文献用DMD方法建立线性讲解模型，再做EVP分析，效果不错。但是前提是，流动快照处于小扰动线性区，如果流动本身饱和了，就不能这么搞。
DMD的实质是Koopman算子，是一个线性算子，可以作用在非线性系统上，所以饱和的系统，湍流情况下是可以使用的。
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Thu Aug 14 11:45:43 2014) 提到:
re
可以看看Lyapunov spectrum 和 特征值有一定的关系 但是不完全是特征值
还有length scale of chaotic attractor
但是湍流里面还是存在线性失稳的问题
一般是通过滤波，去掉所有的高频扰动
然后对比低Re 的线性失稳流动，基本可以确定是不是一个机制
其实这个方法还是回到了POD或者DMD里面的第一模态
【 在 wealth (阳光普照) 的大作中提到: 】
: Yes, you are right. Linear (local and global) stability analysis is in terms of laminar flow. If the base flow is already unstable or even turbulent then it is better to do time-stepping calculation.
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Thu Aug 14 11:47:24 2014) 提到:
FSD的方法，实际流场，完全可以使用
Newton迭代法也可以用于湍流
【 在 shookware (@viocy) 的大作中提到: 】
: 你这说的都是纯理论模型问题。。
: 他这个是个实际的可压缩流动 对特定频率做滤波的特定频率没办法选的。。
: 算出来什么就是什么的那种。。
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Thu Aug 14 12:45:58 2014) 提到:
就是不稳的
所以基本上对应于一个强剪切流动 你用湍流的平均流场来分析 还可以
例如射流或者混合层这种的 因为在这类流动状态中 流动的失稳扰动的增长率都非常大
而尺度差异也很明显 所以没啥问题……
如果不是 那麻烦就比较大了 我觉得不一定可以 需要试试才知道尺度是否存在一致性。
。 但是增长率肯定是不对的
【 在 LaPi (坎切拉皮) 的大作中提到: 】
: 湍流应该是一种线性小扰动充分发展以后的现象了。
: 如果说基本流是怎么样，从道理上讲对应的是流动变成湍流之前的层流，
: 或者说失稳之前的流动。
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Thu Aug 14 12:49:54 2014) 提到:
求介绍FSD相关的处理办法的文献……
-_-
另外 如果有POD的经典的 最好也帮忙推荐一下~~
【 在 yihongyuan (怡红之主) 的大作中提到: 】
: FSD的方法，实际流场，完全可以使用
: Newton迭代法也可以用于湍流
☆─────────────────────────────────────☆
mcx (娶柠檬妻) 于 (Sat Aug 23 10:19:30 2014) 提到:
大神好，这个你都会破？那你会破空化和空蚀不?
【 在 shookware (@viocy) 的大作中提到: 】
: 标 题: Re: help:流动稳定性分析
: 发信站: 水木社区 (Tue Aug 12 10:41:55 2014), 站内
:
: 比如剪切射流？
: 设个假的理想解……
: 然后稳定性能算出来那个周期变化的尺度……
:
: 什么？你说湍流里面的相干结构？ 我建议你选择平均流……
:
: 理论上大尺度的失稳结构即便流动是非稳定的 是可以用平均解解析这个非稳结构的
: 因为这个时候的增长率都普遍较大
:
:
: 【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: : 我知道有无粘线性稳定性，层流线性稳定性
: : 但是对于不存在稳定解的问题，怎么选择basic flow?
: : 某些流动有很强的周期变化，选择任何状态作为basic flow都会与小扰动矛盾，怎么
: : ...................
:
: --
:
: ※ 来源:·水木社区 newsmth.net·[FROM: 60.25.235.*]
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Sat Aug 23 11:24:04 2014) 提到:
不会……
【 在 mcx (娶柠檬妻) 的大作中提到: 】
: 大神好，这个你都会破？那你会破空化和空蚀不?
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Sat Aug 23 15:15:05 2014) 提到:
你要做空化和空蚀的什么方面？
数值模拟？ 可压还是不可压？
【 在 mcx (娶柠檬妻) 的大作中提到: 】
: 大神好，这个你都会破？那你会破空化和空蚀不?
☆─────────────────────────────────────☆
mcx (娶柠檬妻) 于 (Sat Aug 23 17:39:58 2014) 提到:
没有，大神抬举了。我只是感兴趣，随便问下，我要做的还是板壳力学方面。
其实我现在对于指定截面形状和长度的管道里头的无滑移边界的局部水头损失很感兴趣，典型问题是一个直径5~500um的孔，长度在2~400um，里头流体是空气，流动的典型速度是1um/s。也不知道有没有办法计算，因为自己根本就不懂。不知道大神方便不方便指导我一下。
【 在 yihongyuan (怡红之主) 的大作中提到: 】
: 标 题: Re: help:流动稳定性分析
: 发信站: 水木社区 (Sat Aug 23 15:15:05 2014), 站内
:
: 你要做空化和空蚀的什么方面？
: 数值模拟？ 可压还是不可压？
: 【 在 mcx (娶柠檬妻) 的大作中提到: 】
: : 大神好，这个你都会破？那你会破空化和空蚀不?
:
:
: --
: 开辟鸿蒙，
: 谁为情种？
: 都只为风月情浓。
: 趁着这奈何天、
: 伤怀日、
: 寂寥时，
:
:
: ※ 来源:·水木社区 newsmth.net·[FROM: 58.49.207.*]
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Sat Aug 23 17:43:25 2014) 提到:
这个不是Poiseuille 流动吗？
Fluent就能搞定啊
【 在 mcx (娶柠檬妻) 的大作中提到: 】
: 没有，大神抬举了。我只是感兴趣，随便问下，我要做的还是板壳力学方面。
: 其实我现在对于指定截面形状和长度的管道里头的无滑移边界的局部水头损失很感兴趣，典型问题是一个直径5~500um的孔，长度在2~400um，里头流体是空气，流动的典型速度是1um/s。也不知道有没有办法计算，因为自己根本就不懂。不知道大神方便不方便指导我一下。
☆─────────────────────────────────────☆
mcx (娶柠檬妻) 于 (Sat Aug 23 17:46:17 2014) 提到:
泊肃叶流动只能算里头啊，我现在发现算的不准，想考察一下进出口处的局部水头损失。fluent也可以算咩？我之前是用ansys flotran算的，里头泊肃叶流动也算出来了，水头损失不会算...
【 在 yihongyuan (怡红之主) 的大作中提到: 】
: 标 题: Re: help:流动稳定性分析
: 发信站: 水木社区 (Sat Aug 23 17:43:25 2014), 站内
:
: 这个不是Poiseuille 流动吗？
: Fluent就能搞定啊
: 【 在 mcx (娶柠檬妻) 的大作中提到: 】
: : 没有，大神抬举了。我只是感兴趣，随便问下，我要做的还是板壳力学方面。
: : 其实我现在对于指定截面形状和长度的管道里头的无滑移边界的局部水头损失很感兴趣，典型问题是一个直径5~500um的孔，长度在2~400um，里头流体是空气，流动的典型速度是1um/s。也不知道有没有办法计算，因为自己根本就不懂。不知道大神方便不方便指导我一下。
:
:
: --
: 开辟鸿蒙，
: 谁为情种？
: 都只为风月情浓。
: 趁着这奈何天、
: 伤怀日、
: 寂寥时，
:
:
: ※ 来源:·水木社区 newsmth.net·[FROM: 58.49.207.*]
☆─────────────────────────────────────☆
logarithm (踢死你) 于 (Sun Aug 24 00:00:12 2014) 提到:
我觉得这个得用分子动力学模拟，这么小尺度空气不能看成连续介质的吧，当然可以试下，误差估计不小。
【 在 yihongyuan (怡红之主) 的大作中提到: 】
: 这个不是Poiseuille 流动吗？
: Fluent就能搞定啊
☆─────────────────────────────────────☆
logarithm (踢死你) 于 (Sun Aug 24 00:03:21 2014) 提到:
看过一篇文章说小尺度里的压力损失主要就在进出口。
【 在 mcx (娶柠檬妻) 的大作中提到: 】
: 泊肃叶流动只能算里头啊，我现在发现算的不准，想考察一下进出口处的局部水头损失。fluent也可以算咩？我之前是用ansys flotran算的，里头泊肃叶流动也算出来了，水头损失不会算...
: ※ 修改:·mcx 于 Aug 23 17:46:55 2014 修改本文·[FROM: 211.99.222.*]
☆─────────────────────────────────────☆
mcx (娶柠檬妻) 于 (Sun Aug 24 08:05:46 2014) 提到:
是啊，我仿真不出来啊，是不是很丢人？
【 在 logarithm (踢死你) 的大作中提到: 】
: 标 题: Re: help:流动稳定性分析
: 发信站: 水木社区 (Sun Aug 24 00:03:21 2014), 站内
:
: 看过一篇文章说小尺度里的压力损失主要就在进出口。
: 【 在 mcx (娶柠檬妻) 的大作中提到: 】
: : 泊肃叶流动只能算里头啊，我现在发现算的不准，想考察一下进出口处的局部水头损失。fluent也可以算咩？我之前是用ansys flotran算的，里头泊肃叶流动也算出来了，水头损失不会算...
: : ※ 修改:·mcx 于 Aug 23 17:46:55 2014 修改本文·[FROM: 211.99.222.*]
:
:
: --
: 水木是这样一个地方，你可以说：‘’spreading the wealth around“，但你不能
: 说：‘’均贫富“
:
:
: ※ 来源:·水木社区 newsmth.net·[FROM: 193.48.172.*]
☆─────────────────────────────────────☆
mcx (娶柠檬妻) 于 (Sun Aug 24 08:06:09 2014) 提到:
分子动力学？不懂...
【 在 logarithm (踢死你) 的大作中提到: 】
: 标 题: Re: help:流动稳定性分析
: 发信站: 水木社区 (Sun Aug 24 00:00:12 2014), 站内
:
: 我觉得这个得用分子动力学模拟，这么小尺度空气不能看成连续介质的吧，当然可以试下，误差估计不小。
: 【 在 yihongyuan (怡红之主) 的大作中提到: 】
: : 这个不是Poiseuille 流动吗？
: : Fluent就能搞定啊
:
:
: --
: 水木是这样一个地方，你可以说：‘’spreading the wealth around“，但你不能
: 说：‘’均贫富“
:
:
: ※ 来源:·水木社区 newsmth.net·[FROM: 193.48.172.*]
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Sun Aug 24 15:32:48 2014) 提到:
这主要是微流体的问题
需要考虑液体表面张力的问题
算法还是有一些的，见过有人拿FLUENT算过类似的问题
准不准就不知道了
【 在 logarithm (踢死你) 的大作中提到: 】
: 我觉得这个得用分子动力学模拟，这么小尺度空气不能看成连续介质的吧，当然可以试下，误差估计不小。
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Sun Aug 24 15:33:40 2014) 提到:
液体表面张力吧？
【 在 logarithm (踢死你) 的大作中提到: 】
: 看过一篇文章说小尺度里的压力损失主要就在进出口。
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Sun Aug 24 16:31:19 2014) 提到:
不仅仅是
如果克努森数很小 连边界条件都要重新搞 即便你认为NS方程是适用的
张力肯定要考虑 但是微流体这个问题的难点不是有张力 而是N-S方程是否有效的问题
【 在 yihongyuan (怡红之主) 的大作中提到: 】
: 这主要是微流体的问题
: 需要考虑液体表面张力的问题
: 算法还是有一些的，见过有人拿FLUENT算过类似的问题
: ...................
☆─────────────────────────────────────☆
logarithm (踢死你) 于 (Sun Aug 24 19:31:39 2014) 提到:
在这个尺度提表面张力，我觉得不太合适了，表面张力也是微观分子作用在宏观的表现，还是要按非连续性介质考虑。
【 在 yihongyuan (怡红之主) 的大作中提到: 】
: 液体表面张力吧？
☆─────────────────────────────────────☆
logarithm (踢死你) 于 (Sun Aug 24 19:35:48 2014) 提到:
这个有不少开源软件，计算量比较大。简单说分子势能加牛顿公式来模拟分子运动，然后用统计的方法得到宏观量。
【 在 mcx (娶柠檬妻) 的大作中提到: 】
: 分子动力学？不懂...
☆─────────────────────────────────────☆
logarithm (踢死你) 于 (Sun Aug 24 19:39:07 2014) 提到:
对，要么ns在一定程度上可以继续凑合，边界基本不能用无滑移边界，得提新的假设。要么干脆直接上分子动力学模拟，不考虑ns适用不适用问题。
【 在 shookware (@viocy) 的大作中提到: 】
: 不仅仅是
: 如果克努森数很小 连边界条件都要重新搞 即便你认为NS方程是适用的
: 张力肯定要考虑 但是微流体这个问题的难点不是有张力 而是N-S方程是否有效的问题
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Sun Aug 24 21:53:08 2014) 提到:
关键是算法问题 算起来比N-S方程慢不少吧？
【 在 logarithm (踢死你) 的大作中提到: 】
: 这个有不少开源软件，计算量比较大。简单说分子势能加牛顿公式来模拟分子运动，然后用统计的方法得到宏观量。
☆─────────────────────────────────────☆
logarithm (踢死你) 于 (Mon Aug 25 01:55:52 2014) 提到:
那慢的不是一星半点，有号称算了一亿个分子而洋洋得意的，你可以想象下有多大计算量
【 在 shookware (@viocy) 的大作中提到: 】
: 关键是算法问题 算起来比N-S方程慢不少吧？
☆─────────────────────────────────────☆
laser2000 (掌上智能版都是大厦挂) 于 (Sun Jan 4 12:00:10 2015) 提到:
最近用burgs方程为模型，试了试DMD
DMD的重复那是非常好，搞得我兴奋了一阵
但是DMD的伴随解求得的敏感度和burgs的伴随解求的敏感度完全不是一码事，差了几个
数量级……白兴奋了
看来，虽然DMD能够把场的演化重复的非常好，但是拟合的线性模型和原始非线性模型完
全不是一回事。说白了，原始模型是线性的或者生成的数据是线性近似的，DMD才有效
果。
【 在 shookware (@viocy) 的大作中提到: 】
: 流动饱和的话 得引入非线性模型了 那这问题就复杂多了。。
: 但是我觉得尺度上用线性模型来考虑应该不会差的太严重
: 但是增长率必然差到姥姥家了。。
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Sun Jan 4 12:06:37 2015) 提到:
肯定啊
DMD推导的时候就是假设系统是线性的啊……
-_-
这个不需要验证也可以推测出来吧？
不过DMD这个对于做局部的稳定性线性分析里面还是有点用的…… 可以直接用DMD提取数
值解数据中的模态和理论解比 省事……
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: 最近用burgs方程为模型，试了试DMD
: DMD的重复那是非常好，搞得我兴奋了一阵
: 但是DMD的伴随解求得的敏感度和burgs的伴随解求的敏感度完全不是一码事，差了几个
: ...................
☆─────────────────────────────────────☆
laser2000 (掌上智能版都是大厦挂) 于 (Sun Jan 4 12:29:04 2015) 提到:
外行，理解不到位
【 在 shookware (@viocy) 的大作中提到: 】
: 肯定啊
: DMD推导的时候就是假设系统是线性的啊……
: -_-
: ...................
-来自 Windows Phone 8 水木社区客户端
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Sun Jan 4 13:45:21 2015) 提到:
那给讲讲……
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: 外行，理解不到位
: -来自 Windows Phone 8 水木社区客户端
☆─────────────────────────────────────☆
laser2000 (掌上智能版都是大厦挂) 于 (Sun Jan 4 14:24:19 2015) 提到:
你误会了
我说我自己外行
【 在 shookware (@viocy) 的大作中提到: 】
: 那给讲讲……
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Sun Jan 4 14:24:52 2015) 提到:
………………
-_-
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: 你误会了
: 我说我自己外行
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Sun Jan 4 16:41:53 2015) 提到:
DMD基不是正交的
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: 最近用burgs方程为模型，试了试DMD
: DMD的重复那是非常好，搞得我兴奋了一阵
: 但是DMD的伴随解求得的敏感度和burgs的伴随解求的敏感度完全不是一码事，差了几个
: ...................
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Sun Jan 4 16:43:18 2015) 提到:
DMD是用的Koopman算符，是线性算符作用在非线性系统上
并没有假设系统是线性的
【 在 shookware (@viocy) 的大作中提到: 】
: 肯定啊
: DMD推导的时候就是假设系统是线性的啊……
: -_-
: ...................
☆─────────────────────────────────────☆
laser2000 (掌上智能版都是大厦挂) 于 (Sun Jan 4 17:20:07 2015) 提到:
1 不是用QR或者SVD计算的么？怎么不正交?
2 我的理解是，既然假设了 x[n+1] = A x[n] 这种形式，难道不是假设系统是线性
的了??
【 在 yihongyuan (怡红之主) 的大作中提到: 】
: DMD是用的Koopman算符，是线性算符作用在非线性系统上
: 并没有假设系统是线性的
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Sun Jan 4 17:23:06 2015) 提到:
没吧？
我看文章 它的推导过程也是假设整个系统得是个线性的 那些玩意才能成立啊……
【 在 yihongyuan (怡红之主) 的大作中提到: 】
: DMD是用的Koopman算符，是线性算符作用在非线性系统上
: 并没有假设系统是线性的
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Sun Jan 4 17:24:16 2015) 提到:
我也这么认为……
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: 1 不是用QR或者SVD计算的么？怎么不正交?
: 2 我的理解是，既然假设了 x[n+1] = A x[n] 这种形式，难道不是假设系统是线性
: 的了??
: ...................
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Sun Jan 4 17:26:38 2015) 提到:
看Rowley还有mezic的文章吧
一句话说不清楚
不能只看Peter Schmid的文章，数学背景都没有交代清楚
他上课也是如此，大家都觉得好清楚啊，回来一看，后面好多东西的
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: 1 不是用QR或者SVD计算的么？怎么不正交?
不是QR 是SVD但是对应的不是线性方程的正交基，如果你用的是非线性系统
: 2 我的理解是，既然假设了 x[n+1] = A x[n] 这种形式，难道不是假设系统是线性
: 的了??
不是假设系统是线性的，只是假设recurrence可以写成这种形式
: ...................
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Sun Jan 4 17:27:31 2015) 提到:
你说的是Schmid的JFM那篇吧？
【 在 shookware (@viocy) 的大作中提到: 】
: 没吧？
: 我看文章 它的推导过程也是假设整个系统得是个线性的 那些玩意才能成立啊……
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Sun Jan 4 17:41:59 2015) 提到:
是啊……
【 在 yihongyuan (怡红之主) 的大作中提到: 】
: 你说的是Schmid的JFM那篇吧？
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Sun Jan 4 17:43:29 2015) 提到:
关键问题在于 如果你的系统都是非线性的
那么按照DMD的求解方法 那么你求出来的东西是什么？ 是非线性算子对应的线性算子化的
什么玩意 对吧？
但是这个东西的物理意义是啥？
【 在 yihongyuan (怡红之主) 的大作中提到: 】
: 看Rowley还有mezic的文章吧
: 一句话说不清楚
: 不能只看Peter Schmid的文章，数学背景都没有交代清楚
: ...................
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Sun Jan 4 17:48:13 2015) 提到:
【 在 shookware (@viocy) 的大作中提到: 】
: 关键问题在于 如果你的系统都是非线性的
: 那么按照DMD的求解方法 那么你求出来的东西是什么？ 是非线性算子对应的线性算子化的
: 什么玩意 对吧？
还是非线性系统的几个动力学模态
只是通过一个线性算子做了一个projection
这是我的理解，可能不完全对
mezic的文章，我没有100%理解
不过非常建议看看他的那篇ann jfm
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Sun Jan 4 17:51:07 2015) 提到:
回头看看
但是总得有个可以供比较的对象 就是怎么证明你这个分析是对的的参照物……
反正给我的印象是 如果对于小扰动 本身的算子就是个线性的 那么DMD可以认为是一种
基于数据分析的求解线性特征值的方法…… -_- 如果是非线性的算子 根据DMD给出来的
特征值的虚部的物理意义与DNS给出的对应模态的增长过程是否是一回事？
【 在 yihongyuan (怡红之主) 的大作中提到: 】
: 还是非线性系统的几个动力学模态
: 只是通过一个线性算子做了一个projection
: 这是我的理解，可能不完全对
: ...................
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Sun Jan 4 18:42:06 2015) 提到:
【 在 shookware (@viocy) 的大作中提到: 】
: 回头看看
: 但是总得有个可以供比较的对象 就是怎么证明你这个分析是对的的参照物……
: 反正给我的印象是 如果对于小扰动 本身的算子就是个线性的 那么DMD可以认为是一种
: 基于数据分析的求解线性特征值的方法…… -_-
不仅仅是
：如果是非线性的算子 根据DMD给出来的
: 特征值的虚部的物理意义与DNS给出的对应模态的增长过程是否是一回事？
是，也不仅仅是，看情况
☆─────────────────────────────────────☆
willianhe (溜达) 于 (Mon Jan 5 20:08:16 2015) 提到:
http://www.google.es/url?sa=t&rct=j&q=d.heningson%20spectral%20nek5000&source=web&cd=1&ved=0CCYQFjAA&url=http%3A%2F%2Fwww.springer.com%2Fcda%2Fcontent%2Fdocument%2Fcda_downloaddocument%2F9783319016009-c1.pdf%3FSGWID%3D0-0-45-1429314-p175348417&ei=eX6qVJPTG8LlUsTagKgL&usg=AFQjCNG-dMIUavUQyltHmxjD2MeIVsCyjQ&sig2=qnt6Icj3D0C_rhBhNc5SVA&bvm=bv.82001339,d.d24
不知道是否符合你的要求？
【 在 shookware 的大作中提到: 】
: 求介绍FSD相关的处理办法的文献……
: -_-
: 另外 如果有POD的经典的 最好也帮忙推荐一下~~
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Mon Jan 5 22:42:00 2015) 提到:
文章名？中国大陆没办法访问google
【 在 willianhe 的大作中提到: 】
:
http://www.google.es/url?sa=t&rct=j&q=d.heningson%20spectral%20nek5000&source=web&cd=1&ved=0CCYQFjAA&url=http%3A%2F%2Fwww.springer.com%2Fcda%2Fcontent%2Fdocument%2Fcda_downloaddocument%2F9783319016009-c1.pdf%3FSGWID%3D0-0-45-1429314-p175348417&ei=eX6qVJPTG8LlUsTagKgL&usg=AFQjCNG-dMIUavUQyltHmxjD2MeIVsCyjQ&sig2=qnt6Icj3D0C_rhBhNc5SVA&bvm=bv.82001339,d.d24
: 不知道是否符合你的要求？
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Tue Jan 6 01:21:38 2015) 提到:
Stability Tools for the Spectral-Element Code
Nek5000: Application to Jet-in-Crossflow
【 在 shookware (@viocy) 的大作中提到: 】
: 文章名？中国大陆没办法访问google
☆─────────────────────────────────────☆
willianhe (溜达) 于 (Tue Jan 6 07:14:50 2015) 提到:
Stability Tools for the Spectral-Element Code
Nek5000: Application to Jet-in-Crossflow
A. Peplinski, P. Schlatter, P.F. Fischer, and D.S. Henningson
Abstract
【 在 shookware 的大作中提到: 】
: 文章名？中国大陆没办法访问google
☆─────────────────────────────────────☆
bullet01 ( ) 于 (Tue Jan 6 09:28:29 2015) 提到:
直接贴这儿吧
【 在 yihongyuan (怡红之主) 的大作中提到: 】
: 标 题: Re: help:流动稳定性分析
: 发信站: 水木社区 (Tue Jan 6 01:21:38 2015), 转信
:
: Stability Tools for the Spectral-Element Code
: Nek5000: Application to Jet-in-Crossflow
: 【 在 shookware (@viocy) 的大作中提到: 】
: : 文章名？中国大陆没办法访问google
:
:
: --
: 试遣愚衷。
: 因此上，
: 演出这怀金悼玉的《红楼梦》。
:
:
: ※ 来源:·水木社区 newsmth.net·[FROM: 82.241.173.*]
☆─────────────────────────────────────☆
laser2000 (掌上智能版都是大厦挂) 于 (Tue Jan 6 09:29:59 2015) 提到:
可以访问啊
难道各个学校不一样?
【 在 shookware (@viocy) 的大作中提到: 】
: 文章名？中国大陆没办法访问google
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Tue Jan 6 09:52:57 2015) 提到:
确实不能访问啊……
【 在 laser2000 (掌上智能版都是大厦挂) 的大作中提到: 】
: 可以访问啊
: 难道各个学校不一样?
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Tue Jan 6 09:53:01 2015) 提到:
：）
【 在 yihongyuan (怡红之主) 的大作中提到: 】
: Stability Tools for the Spectral-Element Code
: Nek5000: Application to Jet-in-Crossflow
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Wed Jan 7 01:16:07 2015) 提到:
不过话说回来,这个稳定性的module
只有他们KTH用,网上搜不到
所以只能找他们要,估计完了还得属他们的名
【 在 shookware (@viocy) 的大作中提到: 】
: ：）
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Wed Jan 7 10:12:24 2015) 提到:
暂时不需要。。
这个代码只能做低速问题……
不过我们做低速也有自己的代码…… 无所谓。。
【 在 yihongyuan (怡红之主) 的大作中提到: 】
: 不过话说回来,这个稳定性的module
: 只有他们KTH用,网上搜不到
: 所以只能找他们要,估计完了还得属他们的名
: ...................
☆─────────────────────────────────────☆
feng666 (小枫) 于 (Wed Jan 7 22:09:37 2015) 提到:
请问这个有文章吗？可以共享下？
【 在 shookware 的大作中提到: 】
: 有
: 我没看到 但是据说最近有比较振奋的消息是 用线性稳定性算出来的湍流的相干结构与
: 实验观测到的尺度相似……
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Thu Jan 8 13:55:10 2015) 提到:
可能是近期JFM或者POF上的文章 我当时是听同事老师说的
【 在 feng666 (小枫) 的大作中提到: 】
: 请问这个有文章吗？可以共享下？
☆─────────────────────────────────────☆
willianhe (溜达) 于 (Tue Jan 13 03:15:17 2015) 提到:
代码都是自己写还是基于原来的自己再开发？
【 在 shookware 的大作中提到: 】
: 暂时不需要。。
: 这个代码只能做低速问题……
: 不过我们做低速也有自己的代码…… 无所谓。。
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Tue Jan 13 13:33:21 2015) 提到:
自己写吧
一代一代的传下来的 然后再修改 课题组资产
就是方法上都是simple系 和 涡流函数法 比较原始 呵呵
【 在 willianhe (溜达) 的大作中提到: 】
: 代码都是自己写还是基于原来的自己再开发？
☆─────────────────────────────────────☆
willianhe (溜达) 于 (Tue Jan 13 18:17:24 2015) 提到:
我们组没有自己的code，本来想自己试着写，发现这完全远远超出我的的能力了，还是用开源的了，
紧跟着的问题也来了，到底不清楚人家是具体怎么实现的，怕被人诘问啊
【 在 shookware 的大作中提到: 】
: 自己写吧
: 一代一代的传下来的 然后再修改 课题组资产
: 就是方法上都是simple系 和 涡流函数法 比较原始 呵呵
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Tue Jan 13 21:26:02 2015) 提到:
当年念硕士的时候写了一年的不可压缩流code…… 所以其实很多都烂熟了……
-_-
不过多少年不做不可压缩了 囧~
不过做稳定性的code其实还好 谱方法或者直接高阶格式其实都还差不多~
当然 最好谱方法 省不少网格量
【 在 willianhe (溜达) 的大作中提到: 】
: 我们组没有自己的code，本来想自己试着写，发现这完全远远超出我的的能力了，还是用开源的了，
: 紧跟着的问题也来了，到底不清楚人家是具体怎么实现的，怕被人诘问啊
☆─────────────────────────────────────☆
willianhe (溜达) 于 (Tue Jan 13 21:47:39 2015) 提到:
现在用nektar++，的确比openfoam节省不少计算量
等静下心来 还得写自己的code啊
【 在 shookware 的大作中提到: 】
: 当年念硕士的时候写了一年的不可压缩流code…… 所以其实很多都烂熟了……
: -_-
: 不过多少年不做不可压缩了 囧~
: ...................
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Tue Jan 13 21:49:16 2015) 提到:
能用现成的就用现成的吧
满足需要的前提下 没必要非要用自己的
除非你本身就是做方法的…… 做物理问题的犯不着搞这种费力不讨好的事
【 在 willianhe (溜达) 的大作中提到: 】
: 现在用nektar++，的确比openfoam节省不少计算量
: 等静下心来 还得写自己的code啊
☆─────────────────────────────────────☆
Z5boy (二五仔) 于 (Tue Jan 13 22:53:41 2015) 提到:
外行纯好奇问一下
你们自己写code处理的模型形状是不是都是比较简单的？
譬如就是个矩形或长方体之类的
【 在 shookware (@viocy) 的大作中提到: 】
: 当年念硕士的时候写了一年的不可压缩流code…… 所以其实很多都烂熟了……
: -_-
: 不过多少年不做不可压缩了 囧~
: ...................
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Wed Jan 14 02:43:05 2015) 提到:
nektar++ 貌似还有些bug
而且不支持热对流等 耦合问题
自己编的花，难度不小
doc也不是特别多，除了那本书
感觉不是他们课题组的，只能当黑箱使用了
【 在 willianhe (溜达) 的大作中提到: 】
: 现在用nektar++，的确比openfoam节省不少计算量
: 等静下心来 还得写自己的code啊
☆─────────────────────────────────────☆
willianhe (溜达) 于 (Wed Jan 14 04:18:50 2015) 提到:
的确有些bug，只是因为这个用起来比较简单，如您所说，黑箱子，给个输入拿输出，对错自己心里有时候也打鼓。热对流什么的我暂时还用不到，不过他们发布了个最新版本的，加了一些东西......
不知道Spencer为什么舍弃了老版本的code转而开发新的，莫非是为了增加并行能力？
【 在 yihongyuan 的大作中提到: 】
: nektar++ 貌似还有些bug
: 而且不支持热对流等 耦合问题
: 自己编的花，难度不小
: ...................
☆─────────────────────────────────────☆
yihongyuan (怡红之主) 于 (Wed Jan 14 04:37:00 2015) 提到:
C++重新写的，结构全部重新优化过了
这样不仅可以算流体，还能做磁场电场什么的
主要是求解pde的一个工具箱
但是对于非C++熟练的人员，有些难度，包括我自己
而且他们用到了好多好多的库，boost什么的，程序达到编译都得1个小时
我曾经试图读懂他们的程序，最后还是放弃了
程序太大了，特别是又不在他们课题组，太难自己编了
【 在 willianhe (溜达) 的大作中提到: 】
: 的确有些bug，只是因为这个用起来比较简单，如您所说，黑箱子，给个输入拿输出，对错自己心里有时候也打鼓。热对流什么的我暂时还用不到，不过他们发布了个最新版本的，加了一些东西......
: 不知道Spencer为什么舍弃了老版本的code转而开发新的，莫非是为了增加并行能力？
☆─────────────────────────────────────☆
shookware (@viocy) 于 (Wed Jan 14 09:43:28 2015) 提到:
也可以比较复杂……
无非就是有个度量系数……
当然 这都是结构网格 要是非结构网格 要稍微麻烦一点 但是主要麻烦也是数据结构。
【 在 Z5boy (二五仔) 的大作中提到: 】
: 外行纯好奇问一下
: 你们自己写code处理的模型形状是不是都是比较简单的？
: 譬如就是个矩形或长方体之类的
: ...................
修改:mcx FROM 211.99.222.*
FROM 82.241.173.*