Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Problem incompressible Navier Stokes, using Flow in porous media (Brinkman Equation)

Please login with a confirmed email address before reporting spam

Hi guys,

I have a problem when a use the incompressible Navier Stokes with the activation of Flow in porous media (Brinkman Equation).

I get the next error:

Error:
Failed to find consistent initial values.
Matrix_has_zero_on_diagonal

What can I do?, I have tried everithing but nothing happen.

Regards

Pablo

12 Replies Last Post 2010年10月27日 GMT-4 00:44

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010年6月25日 GMT-4 05:29
Hi Pablo,

This is because of the solver, as the first guess. Switch to a Direct solver if possible to see if the problem is the solver or you have missed some physics. Also, since you enabled another physic, your stiffness matrix changes character and without fine tuning of pre-conditioners for iterative solvers (venka) it will give that error. I'm sure more experienced people can give your more detailed answer.
Hi Pablo, This is because of the solver, as the first guess. Switch to a Direct solver if possible to see if the problem is the solver or you have missed some physics. Also, since you enabled another physic, your stiffness matrix changes character and without fine tuning of pre-conditioners for iterative solvers (venka) it will give that error. I'm sure more experienced people can give your more detailed answer.

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010年6月25日 GMT-4 05:34
Hi, Thank for your help

What you mean with "Direct solver"? Is a option or something in the program?



Hi, Thank for your help What you mean with "Direct solver"? Is a option or something in the program?

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010年6月25日 GMT-4 05:35
Iterative solvers vs. direct solvers. In the solver box, there are several direct solvers (UMFPACK, PARDISO) and some iterative ones. Direct solvers are usually much easier to use than iterative ones.
Iterative solvers vs. direct solvers. In the solver box, there are several direct solvers (UMFPACK, PARDISO) and some iterative ones. Direct solvers are usually much easier to use than iterative ones.

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010年6月25日 GMT-4 05:51
Hi,

Thank a lot i didn't know this part.
I will try with this now and I will tell you the results

regards

Pablo
Hi, Thank a lot i didn't know this part. I will try with this now and I will tell you the results regards Pablo

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010年6月25日 GMT-4 06:13
Hi,

now is running with the "Flow in porous media", thanks a lot for your help

Other question: There are a lot of different between the direct solve and the iteration solve? or is a goos solution too?

regards
Hi, now is running with the "Flow in porous media", thanks a lot for your help Other question: There are a lot of different between the direct solve and the iteration solve? or is a goos solution too? regards

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010年6月25日 GMT-4 07:04
It depends on your degree of freedom (DOF), complexity of your physics, 2D-3D, and plenty more. I have used Direct solvers for 700k DOF, but usually it works the best for <300k. They have less tuning and easier to setup. But for special cases e.g. 3D direct solvers are problematic; both speed and memory usage.
It depends on your degree of freedom (DOF), complexity of your physics, 2D-3D, and plenty more. I have used Direct solvers for 700k DOF, but usually it works the best for

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010年6月29日 GMT-4 07:02
Hi,
I would like to know what solver and what preconditionner that you used for your problem.
Regard
Hi, I would like to know what solver and what preconditionner that you used for your problem. Regard

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010年6月29日 GMT-4 09:12
Hi,

I was using BigStab, but was not working, after i have change to PARDISO and is work.
But i have also the problem that when i try to mix with the Convection and diffusion equation (using the ouput of the navir stokes u,v,w ) only work with small values of the Initial Velocity in the Navier and this value is like 10000 time more small that the real value that i nedd.

What you mean with precondition. which values you need?

Hi, I was using BigStab, but was not working, after i have change to PARDISO and is work. But i have also the problem that when i try to mix with the Convection and diffusion equation (using the ouput of the navir stokes u,v,w ) only work with small values of the Initial Velocity in the Navier and this value is like 10000 time more small that the real value that i nedd. What you mean with precondition. which values you need?

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010年6月29日 GMT-4 10:02
Thanks for your answer.
we have preconditioner when we use indirects solvers like GMRES or FGMRES
Thanks for your answer. we have preconditioner when we use indirects solvers like GMRES or FGMRES

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010年6月29日 GMT-4 10:22
ok, so i don't have in this case.

You have some ideo why is not working with big values in the Initial Velocity and only for small values?

regards
ok, so i don't have in this case. You have some ideo why is not working with big values in the Initial Velocity and only for small values? regards

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010年10月15日 GMT-4 19:31
Which solver do you think suitable incase of solving porous media problem?
Which solver do you think suitable incase of solving porous media problem?

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010年10月27日 GMT-4 00:44
what about the default setting of solver provided by COMSOL40a? is there a downside to using these defaults? it's just unclear what rule of thumb must be applied per problem base. thanks for your suggestion.

--
merci
what about the default setting of solver provided by COMSOL40a? is there a downside to using these defaults? it's just unclear what rule of thumb must be applied per problem base. thanks for your suggestion. -- merci

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.