Hi everyone
Help me please, I am at the end of the abyss. I work on: influence of the addition of surfactant on the displacement of the oil by the water-view of wettability changind and interface tracking. I made the simple model on comsol with functions, variables and equations all well defined.
When I run the simulation, it appears: "Failed to find a solution. Maximum number of Newton iterations reached. Returned solution is not converged."
The problem is probably in PDEs because the coefficients are also dependent on the main variables so incorrect Jacobian matrix I think.
To explain a little bit, I have 4 system of general equation governing:
1 - Naviers stokes equations to describe the flow in my channel (I choose Laminar flow-two phase flow)
Variable dependent is U (velocity)
2- Cahn hilliard equation to track interface for two immiscible fluid (choose phase field)
variable dependent is phipf
3- surfactant transport into water phase during injection (I choose coefficient form PDEs)
Here, main variable is F interface force or volume force related to equation of naviers stockes, but X is also function of phipf which is dependant in time and space
4- surfactant concentration evolution (I choose also coefficient form PDEs)
In the same case, main variable is f which traduced surfactant concentration, but we have delta which is function of "phipf".
Please do something for me. I am enclosing my model and The governing equations are also attached.
Thanks in advance,
Marc Felix