Note: This discussion is about an older version of the COMSOL Multiphysics® software. The information provided may be out of date.

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.

Transport through diluted species-Time dependent problem

Please login with a confirmed email address before reporting spam

Hello,
I am performing a time dependent analysis for transport through diluted species based on the solution I obtain by performing transport through free and porous region in the same geometry. Basically first I solve Naviers stokes and Brinkmann equation in stationary mode and use the velocity profile obtained in the Convection -diffusion equations . I am using Direct Paradiso solver. Using BDF solver, the solution never converges. When I use the generalisd alpha method, it converges for small time interval upto 500 second and gives a good solution.

However when I use a high value for time 10000 seconds, the solution starts rapidly diverging after 1077 seconds. can you please suggest where could be possible error.

Regards
Pranay Ghosh


3 Replies Last Post Apr 24, 2015, 7:01 a.m. EDT
COMSOL Moderator

Hello Pranay Ghosh

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.


Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Jan 17, 2011, 9:40 p.m. EST
Pranay, I also have this problem. I too am modelling a diluted species transport model with a modified navier stokes for flow through a porous media. The time point at which is starts diverging seems random, however does change when timestep and numerical damping is varied (amplification factor for high frequency) in the generalized alpha time stepping method. Sometimes its at 16sec, others it is at 160 second.

If you have managed to find a solution, please update forum- Any help would be greatly appreciated.

Caroline
Pranay, I also have this problem. I too am modelling a diluted species transport model with a modified navier stokes for flow through a porous media. The time point at which is starts diverging seems random, however does change when timestep and numerical damping is varied (amplification factor for high frequency) in the generalized alpha time stepping method. Sometimes its at 16sec, others it is at 160 second. If you have managed to find a solution, please update forum- Any help would be greatly appreciated. Caroline

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Jan 18, 2011, 6:48 a.m. EST

Pranay, I also have this problem. I too am modelling a diluted species transport model with a modified navier stokes for flow through a porous media. The time point at which is starts diverging seems random, however does change when timestep and numerical damping is varied (amplification factor for high frequency) in the generalized alpha time stepping method. Sometimes its at 16sec, others it is at 160 second.

If you have managed to find a solution, please update forum- Any help would be greatly appreciated.

Caroline


Do you have a high Pe and Sc? if so, it seems you have the common problem of divergence in C-D. You can try increasing the mesh where possible, and use artificial diffusion. But yet again, I had to invest a year to get a feeling of my model :(


--
Comsol 4.1
Ubuntu 10.04.1
[QUOTE] Pranay, I also have this problem. I too am modelling a diluted species transport model with a modified navier stokes for flow through a porous media. The time point at which is starts diverging seems random, however does change when timestep and numerical damping is varied (amplification factor for high frequency) in the generalized alpha time stepping method. Sometimes its at 16sec, others it is at 160 second. If you have managed to find a solution, please update forum- Any help would be greatly appreciated. Caroline [/QUOTE] Do you have a high Pe and Sc? if so, it seems you have the common problem of divergence in C-D. You can try increasing the mesh where possible, and use artificial diffusion. But yet again, I had to invest a year to get a feeling of my model :( -- Comsol 4.1 Ubuntu 10.04.1

Please login with a confirmed email address before reporting spam

Posted: 10 years ago Apr 24, 2015, 7:01 a.m. EDT
Hi,

Since this topic is quite old, I re-open it because we currently have the same king of problem with the diluted species solver. Our case is the following:
- the case: an axi-symetric geometry (a simple cylinder)
- first study with bubbly flow: computation of an airlift system to get the fluid velocity field inside the cylinder
- second study with diluted species: a patch of concentration is added at the top of the cylinder, and a time dependent computation is performed. The boundary conditions are no flux everywhere.

The concentration first correctly convects/diffuses inside the cylinder, but after a while the concentration begins to grow exponentially. Looking at the solutions shows that there is actually several mesh points on the walls (on different boundaries) where the concentration diverge and then diffuse in the whole domain.
If we change the boundary condition from no flux to a constant concentration at these locations, no divergence happens.
If we desactivate convection (ie diffusion only), no divergence too.

Any idea or news about this problem ? Is the mesh refinement the only solution to reduce the divergence ? I could provided our case if necessary. Thank you !
Hi, Since this topic is quite old, I re-open it because we currently have the same king of problem with the diluted species solver. Our case is the following: - the case: an axi-symetric geometry (a simple cylinder) - first study with bubbly flow: computation of an airlift system to get the fluid velocity field inside the cylinder - second study with diluted species: a patch of concentration is added at the top of the cylinder, and a time dependent computation is performed. The boundary conditions are no flux everywhere. The concentration first correctly convects/diffuses inside the cylinder, but after a while the concentration begins to grow exponentially. Looking at the solutions shows that there is actually several mesh points on the walls (on different boundaries) where the concentration diverge and then diffuse in the whole domain. If we change the boundary condition from no flux to a constant concentration at these locations, no divergence happens. If we desactivate convection (ie diffusion only), no divergence too. Any idea or news about this problem ? Is the mesh refinement the only solution to reduce the divergence ? I could provided our case if necessary. Thank you !

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.