Hi everyone,
I am a new COMSOL user and I tried to model a two phase flow (CO2 and Water) in porous medium. I did so in other software but I got recommendations to try to solved it in COMSOL . As you can see, the CO2 is injected to the brine-soil medium. However, when I try to analyse, I got the message below and I cannot find the problem.
Attempt to evaluate negative power of zero.
- Function: ^
Failed to evaluate temporary symbolic derivative variable.
- Variable: mod1.kr_w@VDN${abs@4}
- Defined as: (((((mod1.L*(mod1.Se_w^(mod1.L-1)))*(((mod1.Hc>0)*((mod1.N*(abs((mod1.Hc*mod1.alpha))^(mod1.N-1)))*((-mod1.M)*((1+(abs((mod1.Hc*mod1.alpha))^mod1.N))^((-mod1.M)-1)))))))*((1-((1-(mod1.Se_w^(1/mod1.M)))^mod1.M))^2))+((mod1.Se_w^mod1.L)*((-((-(((1/mod1.M)*(mod1.Se_w^((1/mod1.M)-1)))*(((mod1.Hc>0)*((mod1.N*(abs((mod1.Hc*mod1.alpha))^(mod1.N-1)))*((-mod1.M)*((1+(abs((mod1.Hc*mod1.alpha))^mod1.N))^((-mod1.M)-1))))))))*(((1-(mod1.Se_w^(1/mod1.M)))^(mod1.M-1))*mod1.M)))*((1-((1-(mod1.Se_w^(1/mod1.M)))^mod1.M))*2))))*(mod1.Hc>0))
Failed to evaluate temporary symbolic derivative variable.
- Variable: mod1.w.kappaxx@VDN${abs@4}
- Defined as: (((((((mod1.L*(mod1.Se_w^(mod1.L-1)))*(((mod1.Hc>0)*((mod1.N*(abs((mod1.Hc*mod1.alpha))^(mod1.N-1)))*((-mod1.M)*((1+(abs((mod1.Hc*mod1.alpha))^mod1.N))^((-mod1.M)-1)))))))*((1-((1-(mod1.Se_w^(1/mod1.M)))^mod1.M))^2))+((mod1.Se_w^mod1.L)*((-((-(((1/mod1.M)*(mod1.Se_w^((1/mod1.M)-1)))*(((mod1.Hc>0)*((mod1.N*(abs((mod1.Hc*mod1.alpha))^(mod1.N-1)))*((-mod1.M)*((1+(abs((mod1.Hc*mod1.alpha))^mod1.N))^((-mod1.M)-1))))))))*(((1-(mod1.Se_w^(1/mod1.M)))^(mod1.M-1))*mod1.M)))*((1-((1-(mod1.Se_w^(1/mod1.M)))^mod1.M))*2))))*(mod1.Hc>0)))*mod1.kap_s)
Failed to evaluate temporary symbolic derivative variable.
- Variable: mod1.w.u@VDN${abs@4}
- Defined as: ((((((((((mod1.L*(mod1.Se_w^(mod1.L-1)))*(((mod1.Hc>0)*((mod1.N*(abs((mod1.Hc*mod1.alpha))^(mod1.N-1)))*((-mod1.M)*((1+(abs((mod1.Hc*mod1.alpha))^mod1.N))^((-mod1.M)-1)))))))*((1-((1-(mod1.Se_w^(1/mod1.M)))^mod1.M))^2))+((mod1.Se_w^mod1.L)*((-((-(((1/mod1.M)*(mod1.Se_w^((1/mod1.M)-1)))*(((mod1.Hc>0)*((mod1.N*(abs((mod1.Hc*mod1.alpha))^(mod1.N-1)))*((-mod1.M)*((1+(abs((mod1.Hc*mod1.alpha))^mod1.N))^((-mod1.M)-1))))))))*(((1-(mod1.Se_w^(1/mod1.M)))^(mod1.M-1))*mod1.M)))*((1-((1-(mod1.Se_w^(1/mod1.M)))^mod1.M))*2))))*(mod1.Hc>0)))*mod1.kap_s))*((-mod1.p_wx)-((d(mod1.w.D,x)*mod1.w.rho)*mod1.w.g)))*(1/mod1.w.mu))
Failed to evaluate expression.
- Expression: d(d((mod1.w.rho*(-mod1.w.S*mod1.p_wt*test(mod1.p_w)+mod1.w.u*test(mod1.p_wx)+mod1.w.v*test(mod1.p_wy)))*(dvol),{test@2}),{abs@4})
I tried to find and fix these errors but I am not fully understand how to do it. Does anyone know what going wrong?
The model file is attached below
Thanks a lot for any help!
Raviv