|
|
发表于 2003-4-24 03:22:54
|
显示全部楼层
关于坐标变换后的差分格式问题
对与curvilinear finite volume, 将多出来的几项做为源项处理。如果physical domain 到computational domain的变型比较大,源项就比较大,稳定性就差了,这时FVM不好处理, 用FEM就有优势。
For example, laplace equation d^2T/dx^2+d^2T/dy^2=0, the final result after transformation can be written as:
x=x(k,e,t)
y=y(k,e,t)
A(d^2T/dk^2)+B(d^2T/de^2)+C(d^2T/dkde)+D(dT/dk)+E(dT/de)=0
A, B, C, D, E are coefficients functions of x,y, independent of T;
if x=x(k), y=y(e) only, then C=D=E=0
if 变型比较小, then A, B are large, and C, D, E are small, relatively stable.
if 变型比较大, then A, B are small relatively compared to C, D, E. Because C, D, E is 源项for FVM. =>relatively unstable.
----------
抛砖引玉。
|
|