
[Sponsors] 
June 12, 2012, 06:05 

#121 
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 303
Rep Power: 9 
Hi anton,
i think i do not understand well your first question. what do you mean with "start filtering"?. I would first define grad(pc)(grad(pc) dot ns)ns at the cell centers, then i would interpolate it to face centers and dot product with "nf". We get a surfaceScalarField define at the face centers that can be multiplied with Cfc (which is a constant) and deltasf/(mag(deltasf+eps)) which is also a surfaceScalarField. Once fcf_filtered is calculated, fcf can be updated using (21). I would put the filtering just after the capillary pressure equation in interFoamSSF.C I think eps is used just to avoid division by zero. I would put some small numbers such as 1e6/1e8. Since neither fcf nor grad(pc) change in PISO loop i would put the filtering of capillary fluxes just before UEqn.H and just after the previous filtering. I would define the total capillary flux phi_c as it is defined in pag.9, then i would apply the filtering and i would use directly phi_c in Ueqn.H and pEqn.H, instead of using fcf and grad(pc) separeted. i am very sorry but i am preparing for a conference and i have no time to put my hands in the code now. i'll be back in 2 weeks best andrea 

June 12, 2012, 06:49 

#122  
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,188
Rep Power: 21 
Quote:
Quote:
Quote:
 Anton P.S. I'll post some code later today.
__________________
*On twitter @akidTwit *Spend as much time formulating your questions as you expect people to spend on their answer. 

June 12, 2012, 07:11 

#123 
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,188
Rep Power: 21 
Filtering code is now online: http://code.google.com/r/akidessinterfoamfsf/
__________________
*On twitter @akidTwit *Spend as much time formulating your questions as you expect people to spend on their answer. 

June 13, 2012, 03:20 

#124 
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,188
Rep Power: 21 
Results on the moving droplet case with filtering improved, but are nowhere as good as the reference :/
__________________
*On twitter @akidTwit *Spend as much time formulating your questions as you expect people to spend on their answer. 

June 13, 2012, 06:13 

#125 
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 303
Rep Power: 9 
just a small comment, at pag.9 they say that fcf_average is the average of capillary forces over all faces where they are nonzero. Does cmptAv account for that?
best andrea 

June 13, 2012, 09:10 

#126 
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,188
Rep Power: 21 
I believe so. Btw, using the filtering algorithm as published on the repo, even the static droplet testcase doesn't work any more. After deactivating the filtering of fcf (but keeping it on phic) the static droplet case works again, so there's a problem with the computation of fcf_filter.
__________________
*On twitter @akidTwit *Spend as much time formulating your questions as you expect people to spend on their answer. 

August 1, 2012, 10:51 
InterFoam: parasitic currents

#127 
New Member
M K Singh
Join Date: Sep 2009
Posts: 19
Rep Power: 9 
Dear Anton and Andrea,
Very nice discussion about parasitic currents. I am wondering if all the issues have been solved regarding spurious velocities at interfaces in interFoam. I am using OpenFoam 2.1.x to solve a 2D ladle draining from the bottom (water +air system). The initial condition is shown in Figure attached (t=0 s). I am getting very strange results, please see the attached pictures for t=100,120, 170 s. As you can see that after certain time (t>100s) the liquid level remains constatnt and a strange waves are appearing at the inteface. I am doubting that is it effect of spurious velocities generated at intefcae? My boundary conditions are as follows: ===================================== Valid fields: volVectorField U volScalarField p_rgh volScalarField alpha1 patch: SENOUTLET scalar p_rgh totalPressure scalar alpha1 inletOutlet vector U pressureInletOutletVelocity wall: SENWALL scalar p_rgh buoyantPressure scalar alpha1 zeroGradient vector U fixedValue empty: SYM4 scalar p_rgh empty scalar alpha1 empty vector U empty empty: SYM3 scalar p_rgh empty scalar alpha1 empty vector U empty empty: SYM2 scalar p_rgh empty scalar alpha1 empty vector U empty empty: SYM1 scalar p_rgh empty scalar alpha1 empty vector U empty patch: TW scalar p_rgh totalPressure scalar alpha1 inletOutlet vector U pressureInletOutletVelocity wall: BW scalar p_rgh buoyantPressure scalar alpha1 zeroGradient vector U fixedValue wall: SW scalar p_rgh buoyantPressure scalar alpha1 zeroGradient vector U fixedValue ====================================== Any suggestion in this regard? As suggested on forum, Now I have changed cAlpha=1 to see the effect of this on my result. Thanks in advance. M K Singh 

August 1, 2012, 17:53 

#128 
New Member
Ali Q Raeini
Join Date: Feb 2010
Posts: 21
Rep Power: 9 
Dear Singh,
I was also used to get some nonphysical velocities near the wall boundaries, much higher than the usual spurious velocities. You can call them whatever you like! No, I think the filtering stuff are not still implemented in the interFoam code shipped with official openfoam, but you can check the links above. By the way, I am also going to release my code with a tutorial case as soon I find the time to clean it up, maybe in a few months time. Bests, Ali 

August 3, 2012, 04:38 

#129 
Senior Member
Eelco van Vliet
Join Date: Mar 2009
Location: The Netherlands
Posts: 123
Rep Power: 11 
Hi Ali
I just saw the paper of yours in J Comp Physics on this solver. Very nice work! Would be great if you would make your code public, I am sure a lot of us would like to use it. Are you going to post it on the forum? Would a 'prerelease' be possible? Regards Eelco 

August 3, 2012, 06:09 
Solution of my problem (emptying ladle)

#130  
New Member
M K Singh
Join Date: Sep 2009
Posts: 19
Rep Power: 9 
Quote:
I have found solution to my problem why my interface was getting stagnant at certain point. After analysing my simulation, I realized that my inteface was getting stagnant exactly at zero position (y=0, vertical direction, note my middle of domian is at (0,0,0)), see picture for t=10, 50,90 s. My guess is that during presure (rho*g*h) caluculation as interface reaches at y=0 the pressure is becoming zero which is leading to stagnant inteface!! However, this needs more explaination, may be some numerical experts can provide better argument. To solve the issue, I transformed my domain so that y=0 lies outside the domain, basically my all the domain lies in positive quadrant. This strategy has solved my problem, see attached pictures for t=10,60 s. If someone has more input in this regard, please share. Thanks. M K Singh PS: note that ladle_translation is picture where domain has been translated Last edited by mksingh; August 3, 2012 at 06:15. Reason: clarity 

February 20, 2013, 06:03 

#131 
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 303
Rep Power: 9 
Hi interFoamers,
Someone has advanced with the implementation of interFoamSSF? I finally have some time to spend on it but i'm still confused about filtering stuff. best andrea 

May 16, 2013, 03:02 

#132 
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,188
Rep Power: 21 
Ali Raeini has released some code now: http://www.cfdonline.com/Forums/ope...tml#post427813
__________________
*On twitter @akidTwit *Spend as much time formulating your questions as you expect people to spend on their answer. 

January 7, 2015, 16:45 
parasitic velocities in capillary flow

#133 
New Member
Christoph Gmelin
Join Date: Jan 2015
Posts: 1
Rep Power: 0 
Hi,
I'm trying to simulate a capillary twophase flow (water/Oil), where small droplets pass a junction with a sudden expansion (see first picture). U =0.05m/s, tube diameter = 0.5mm. My originalsettings are the same as in the dambreak tutorial and I'm still using interfoam. Everything is Ok, until the droplets reach the expansionzone. Then large Courant numbers occur in combination with high velocities at the interface of the first droplet (see red dot at picture 2). In my opinion the complete case is only almost stable (it's a very fragile equilibrium) an when the flow becomes a little more unstable due to its deceleration spurious currents occur. To resolve the situation I tried the following things (everything without success): reduce courant No. set cAlpha to 0 use finer grids at the wall finer grid in the expansion zone reduce time step more Alpha subcycles ... Surprisingly the unchanged dambreak settings still obtain the best results. In all other cases the parasitic velocities arise much earlier. Has anybody more Ideas how to get this problem fixed?? My last idea was to test the interfoamssfsolver (https://code.google.com/p/interfoamssf/) but I'm not very familiar with Cprogramming and compiling. My second question is: How can I compile the interfoamssf source code without destroying the rest of my FOAMinstallation? Thanks, Christoph 

Tags 
capillary flows, interfoam, parasitic currents 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
how to monitor free surface elevation vs time in OF?  ozgur  OpenFOAM PostProcessing  56  September 14, 2015 08:11 
parasitic currents  PeiYing Hsieh  Main CFD Forum  0  January 13, 2009 20:58 
Parasitic currents reduction  hsieh  OpenFOAM Running, Solving & CFD  0  January 13, 2009 16:44 
Parasitic currents reduction  hsieh  OpenFOAM Running, Solving & CFD  0  January 13, 2009 16:37 
Modelling ocean currents of the past Earth  pgm  Main CFD Forum  3  March 2, 2005 09:45 