例子:
速度场:
v =
- x * sin(c * t) - y *(w / 2 - cos(c * t))
Y * sin(c * t)+ x *(w / 2 + cos(c * t))
0.
线性速度场可以写为:
x = [x,y,z];
a = [-sin(c * t),cos(c * t) - w / 2,0;cos(c * t)+ w / 2,sin(c * t),0; 0,0,0;];
v = a *转置(x);
转型(观察者改变):
在哪里
q = [cos(c / 2 * t),sin(c / 2 * t),0; -sin(c / 2 * t),cos(c / 2 * t),0; 0,0,1];
基于连续体力学,速度场变换为
。将此公式应用于:
VY = Diff(Q,T,1)*翻转(x)+ q * v;
我没有得到正确的答案。我相信这是因为v
和
还在旧框架里。如果我将它们更改为它的新帧。我使用的程序如下:
如果
也是线性的,那么它可以写作:
然后
B.
计算为:
qtr =转置(q);假设(QTR,'真实的')';
Q *(a * qtr - diff(qtr,t,1))
它有效。但是,现在问题是,我无法定义
一种
和
B.
对于非线性速度字段,对于不可能的数字速度字段根本不可能,所以我该怎么办?如何使用给定观察者更改转换速度字段?