求助: 解有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取得足够小了。
