Menu

Discussion Forum

New Discussion
Filter

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.

Schiller Neuman SLIP Velocity Mixture Model V3.5a

Posted Jun 15, 2011, 2:39 p.m. EDT Fluid & Heat, Computational Fluid Dynamics (CFD), Studies & Solvers Version 4.0a, Version 5.2 6 Replies

Rui Silva
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Hi folks.

I have trying to use the Schiller Neuman Slip velocity model for a few days but no success.

I am successful in using the Hadamard-Rybczynski but I need to use the Schiller Model for comparison
purposes.

However I get the error 7248 - Undefined value found.

But I can not understand why...

The reason I am using V3.5a is because Helpdesk have said to me that the Schiller model will only work
in version V4.2.

Hope someone can help, because I am really stuck.

The file is attached.

Cheers.

Attachments:

6 Replies Last Post Apr 5, 2017, 10:58 a.m. EDT
COMSOL Moderator

Hello Rui Silva

Your Discussion has gone 30 days without a reply. If you still need help with COMSOL and have an on-subscription license, please visit our Support Center for help.

If you do not hold an on-subscription license, you may find an answer in another Discussion or in the Knowledge Base.

Cathal Cummins
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Aug 12, 2011, 10:56 a.m. EDT
Hi there,


I am running a mixture model in comsol 4.0a.

* laminar
* low gas concentration
* no slip on all walls except top gas vent, which has slip
* gravity

the particle Reynolds number is about 3 and the sedimenting dispersed phase are about 100 microns diameter. Hence, while the Hadamard-Rybczynski model gives reasonable results, it is far from ideal. I had the model working fine with this drag model however, I feel that a Schiller-Naumann model would suit better. This means solving an implicit equation for u_slip which is okay but comsol tells me that it can't evaluate Jacobian:

Failed to evaluate variable Jacobian.
- Variable: mod1.mm.udz
- Geometry: 1
- Boundary: 3

Boundary 3 is the top of the geometry. However, it works fine for other drag models.

The initial conditions are 5% void fraction, zero velocity everywhere both gas and liquid. My question is, when you choose SN drag, is there anything else you need to specify to help it solve the implicit equation?
Hi there, I am running a mixture model in comsol 4.0a. * laminar * low gas concentration * no slip on all walls except top gas vent, which has slip * gravity the particle Reynolds number is about 3 and the sedimenting dispersed phase are about 100 microns diameter. Hence, while the Hadamard-Rybczynski model gives reasonable results, it is far from ideal. I had the model working fine with this drag model however, I feel that a Schiller-Naumann model would suit better. This means solving an implicit equation for u_slip which is okay but comsol tells me that it can't evaluate Jacobian: Failed to evaluate variable Jacobian. - Variable: mod1.mm.udz - Geometry: 1 - Boundary: 3 Boundary 3 is the top of the geometry. However, it works fine for other drag models. The initial conditions are 5% void fraction, zero velocity everywhere both gas and liquid. My question is, when you choose SN drag, is there anything else you need to specify to help it solve the implicit equation?

Attachments:

Cathal Cummins
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Aug 12, 2011, 5:40 p.m. EDT
According to support there's a bug in the slip velocity model in 4.0a?

Could this be the reason why the Schiller-Naumann model wasn't working?

Cheers,
According to support there's a bug in the slip velocity model in 4.0a? Could this be the reason why the Schiller-Naumann model wasn't working? Cheers,
Cathal Cummins
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Aug 15, 2011, 2:17 p.m. EDT
Hi,

Just in case anyone has the same problem, my problem was simply a badly defined initial pressure distribution. I figured a zero pressure ic would suffice but it didn't. Specifying the hydrostatic pressure p(z)=rho g (H-z) avoided the singularity in the slip velocity which happens when p is identically zero everywhere.

Hope this helps!
Hi, Just in case anyone has the same problem, my problem was simply a badly defined initial pressure distribution. I figured a zero pressure ic would suffice but it didn't. Specifying the hydrostatic pressure p(z)=rho g (H-z) avoided the singularity in the slip velocity which happens when p is identically zero everywhere. Hope this helps!
Rui Silva
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Aug 15, 2011, 7:27 p.m. EDT
Hi Cathal,

sorry it took me so long to reply but I have been away on vacation.

To me the problem was solved using version 4.2 and using an estimative for the slip
velocity using the terminal velocity expression.

I don't have the references with me but If you need them I will gladly provide them.

Hope this helps.

Cheers!
Hi Cathal, sorry it took me so long to reply but I have been away on vacation. To me the problem was solved using version 4.2 and using an estimative for the slip velocity using the terminal velocity expression. I don't have the references with me but If you need them I will gladly provide them. Hope this helps. Cheers!
Cathal Cummins
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Aug 16, 2011, 2:39 a.m. EDT
Thank you for your help Rui. As per the last post, it was actually the initial condition for the pressure (also the slip squared velocity) which was causing me the problems.

Initially I ran the H-R model which, since the slip velocity vector is directly proportional to the pressure gradient, could handle p_{t=0}=0 everywhere this doesn't carry over to the S-N. In any case, as you say there is a bug in the older versions of comsol.

Thanks again!
Thank you for your help Rui. As per the last post, it was actually the initial condition for the pressure (also the slip squared velocity) which was causing me the problems. Initially I ran the H-R model which, since the slip velocity vector is directly proportional to the pressure gradient, could handle p_{t=0}=0 everywhere this doesn't carry over to the S-N. In any case, as you say there is a bug in the older versions of comsol. Thanks again!
Andres Arevalo
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Posted: 7 years ago Apr 5, 2017, 10:58 a.m. EDT
Hi,
I´m using the machinery rotating, mixture model and I have a problem with Schiller-Naumann slip model, When I use Schiller-Naumann slip velocity in the model, the program produce a error.

"Failed to find consistent initial values"

What happen? I had changed the initial value for pressure and squared slip velocity and nothing! How can I fix this problem?

Thanks for help me!
Hi, I´m using the machinery rotating, mixture model and I have a problem with Schiller-Naumann slip model, When I use Schiller-Naumann slip velocity in the model, the program produce a error. "Failed to find consistent initial values" What happen? I had changed the initial value for pressure and squared slip velocity and nothing! How can I fix this problem? Thanks for help me!

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.

Suggested Content