求助: 解有mixed derivative的PDE时遇到麻烦# Computation - 科学计算
y*i
1 楼
解一个2维PDE:
df/dt = D1*f_x + D11*f_xx + D2*f_y + D22*f_yy + D12*f_xy
f_x代表df/dx,f_xx代表d^2f/dx^2,f_xy是mixed derivative,代表d^2f/dxdy。
D11,D22都是正的,D11*D22>D12^2
在有些区域,D1, D2很大。
用explicit finite difference解,发现run了很长时间后有些区域f开始有negative值,
物理上f应该一直是正的。不知道为什么?用stability analsys,dt取得足够小了。
df/dt = D1*f_x + D11*f_xx + D2*f_y + D22*f_yy + D12*f_xy
f_x代表df/dx,f_xx代表d^2f/dx^2,f_xy是mixed derivative,代表d^2f/dxdy。
D11,D22都是正的,D11*D22>D12^2
在有些区域,D1, D2很大。
用explicit finite difference解,发现run了很长时间后有些区域f开始有negative值,
物理上f应该一直是正的。不知道为什么?用stability analsys,dt取得足够小了。