两相可混溶驱动问题中近似压力方程的一种新方法

1　介　　绍ΩR2为凸多边形区域,Ω上的两相可混溶驱动问题可由以下微分方程系统来描述a)-.[a(x,c)(p-r(c)]=.u=q,b)φ(x)ct+u.c-(Dc)=(c-c)q=g(c),(1.1)其中a(x,c)=-k(x)μ(c),k(x)为介质的渗透率,μ(c)为流体的粘度,p为流体的压力,φ(x)为介质的孔隙度,c为一相流体的体积浓度,q为外部流体的体积流速,且满足相容性条件∫Ωqdx=0.D是2×2阶矩阵,D=φ(x)[dmI+|u|(dlE(u)+dtE⊥(u))],E(u)=(uiuj/|u|2)2×2,dm为分子扩散系数,dl,dt分别为横向、纵向弥散系数.系统的边界条件、初始条件:n为边界单位外法向a)u.n=0,(x,t)∈Ω×Jb)2i,j=1Dij(…