Hello All~
I conduct momentum transport and mass transports coupled simulation using matlab script~
but i have some problems~
when i simulated this topic, always failed to find a solutions~
====================================================
Warning:
Failed to find a solution:
No convergence, even when using the minimum damping factor.
Returned solution has not converged.
====================================================
[my script(geneal form formulation)] - momentum & mass transport (diffuction & convection)
fem.equ.ga={ { {'-p+2*eta*ux' 'eta*(uy+vx)'} {'eta*(uy+vx)' '(-p+2*eta*vy)'} {0 0} {'-D*cx-D*cy' '-D*cx-D*cy'} } };
fem.equ.f={ {'rho*(u*ux+v*uy)' 'rho*(u*vx+v*vy)' 'ux+vy' '-u*cx-v*cy'} };
fem.bnd.ind={[4] [1] [3] [2]}; % rectangular geometry inlet, wall, membrane wall, outlet
fem.bnd.r={ {'u+6*s*(1-s)*u0*nx' 'v+6*s*(1-s)*u0*ny' 0 '-c+c0'}...
{'u' 'v' 0 0} {'u' 'v-(deltap-fos*c)/Rm' 0 0} {0 0 'p' 0} };
fem.bnd.g={{0 0 0 0} {0 0 0 '-(-nx*u-ny*v)*c*rej'} {0 0 0 '-(-nx*u-ny*v)*c*rej'} {'nx*p' 'ny*p' 0 0}};
what is problems of my script or process?
thanks you~
I conduct momentum transport and mass transports coupled simulation using matlab script~
but i have some problems~
when i simulated this topic, always failed to find a solutions~
====================================================
Warning:
Failed to find a solution:
No convergence, even when using the minimum damping factor.
Returned solution has not converged.
====================================================
[my script(geneal form formulation)] - momentum & mass transport (diffuction & convection)
fem.equ.ga={ { {'-p+2*eta*ux' 'eta*(uy+vx)'} {'eta*(uy+vx)' '(-p+2*eta*vy)'} {0 0} {'-D*cx-D*cy' '-D*cx-D*cy'} } };
fem.equ.f={ {'rho*(u*ux+v*uy)' 'rho*(u*vx+v*vy)' 'ux+vy' '-u*cx-v*cy'} };
fem.bnd.ind={[4] [1] [3] [2]}; % rectangular geometry inlet, wall, membrane wall, outlet
fem.bnd.r={ {'u+6*s*(1-s)*u0*nx' 'v+6*s*(1-s)*u0*ny' 0 '-c+c0'}...
{'u' 'v' 0 0} {'u' 'v-(deltap-fos*c)/Rm' 0 0} {0 0 'p' 0} };
fem.bnd.g={{0 0 0 0} {0 0 0 '-(-nx*u-ny*v)*c*rej'} {0 0 0 '-(-nx*u-ny*v)*c*rej'} {'nx*p' 'ny*p' 0 0}};
what is problems of my script or process?
thanks you~