CFD Online Logo CFD Online URL
Home > Forums > CFX

Instable interface: capillary force driven imbibition of water in a porous domain

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   November 15, 2014, 05:03
Default Instable interface: capillary force driven imbibition of water in a porous domain
New Member
Join Date: Nov 2014
Posts: 7
Rep Power: 4
joha is on a distinguished road

I'm trying to model the imbibition of water in a porous domain. Therefore I set up a simple "1D" model.

My geometry has the following dimensions: 5x5x25 mm, the mesh consists of 1000 stacked cells with symmertries to left+right, front+back, an inlet at the bottom and an opening on the top.

The pressure drop in the porous domain is modelled by setting the appropriate permeabilty. In addition to that I configured a subdomain to model the capillary and the gravity force. Both are take into account as a bulk source by the following term in z-direction:

g*density -C [N/m^2]*WATER.vfs.grad z
C is a material constant (=f(porosity, surfacetension,.....))

The simulation works fine for a few thousand iterations. After that I get a verry blurred asymmetric interface (before the interface gets blurry the pressure field becomes asymmetric):

Detailed view with visible mesh.

I already tried various settings without success:
  • smaller timesteps
  • higher mesh resolution
  • switching inferface compression (0, 1, 2)
  • a different mesh (2x2x1000 cells, same geometry) with a free slip wall instead of a symmetrie boundary
  • serial run (CFX documentation mentions possbile instabilities if the interface crosses the partitition boundary)
The simulations ends with no errors if I remove the capillary source term and set the same pressure at the inlet.

How is it even possible to get this kind of asymmetry in a "1D" mesh?

I would appreciate any kind of help. Thanks!
joha is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Setting rotating frame of referece. RPFigueiredo CFX 3 October 28, 2014 05:59
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32
Modelling capillary forces in porous media s3s4m CFX 3 October 27, 2011 06:12
Force can not converge colopolo CFX 13 October 4, 2011 22:03
User-defined force transfer across FSI interface vivekcfd CFX 7 March 10, 2010 21:43

All times are GMT -4. The time now is 10:16.