interFoam - validation for bubble/droplet flows in microfluidics
It is a bit ashamed but I would like to share one of our work on the validations of interFoam, titled "Benchmark numerical simulations of segmented two-phase flows in microchannels using the Volume of Fluid method". In this work, we presented three benchmark cases - the steady motion of bubbles in a straight two-dimensional channel, the formation of bubbles in two- and three-dimensional T-junctions, and the breakup of droplets in three-dimensional T-junctions - and provided guidelines to set appropriate numerical settings for a simulation of microbubble/microdroplet flows.
Please find here the link to the paper: http://dx.doi.org/10.1016/j.compfluid.2013.06.024
If you can not download the paper, I am more than happy to send you a copy of it.
With best regards,
Hi Dear Duong
Thanks for your good paper.
I am new in interFoam. I have two simple questions.
1) how can we apply the smoother function (Eq.9 in your paper) in our problem? I think it can be set: system/fvSolution dictionary, PISO loop subsection was set with m corrections (nCorrectors)
is it true?
2) where should be determined the value of adjustable coefficient "Cγ in Eq.7"?
Could you please share the case setup for openfoam? Best regards
I just got a bit of free time lately to reorganize my data. Here you can find the code and the test cases I have run for my validation. Please be aware that all these validations have been performed with OF-1.6-ext.
I'm not sure but I think that the function fvc::average already exists and do the same operation as your "smooth function".
In the OpenFOAM C++ doc, you can read:
volField = fvc::average(SurfaceField)
Area-weighted average a surfaceField creating a volField.
So you can do your smoothing by :
At the time I implemented that, I would like to have controls on the level of smoothing (averaging) and also to be able to implement different smoothers (which I did not find superior to that simple Laplacian). So that's why you have such a piece of code.
I just downloded your file (bubble_generation). I am able to blockMesh.
Number of internal faces : 44
Number of boundary faces : 56
Number of defined boundary faces : 56
Number of undefined boundary faces : 0
Checking patch -> block consistency
Creating block offsets
Creating merge list
Further when I run "interFoam" solver I get this error.
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.4.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
Build : 2.4.0-f0842aea0e77
Exec : interFoam
Date : Sep 30 2015
Time : 12:12:41
Host : "HP-WS3"
PID : 10778
Case : /home/paramesh/OpenFOAM/paramesh-2.4.0/run/tutorials/incompressible/pimpleFoam/bubble_generation
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create mesh for time = 0
--> FOAM FATAL IO ERROR:
keyword PIMPLE is undefined in dictionary "/home/paramesh/OpenFOAM/paramesh-2.4.0/run/tutorials/incompressible/pimpleFoam/bubble_generation/system/fvSolution"
file: /home/paramesh/OpenFOAM/paramesh-2.4.0/run/tutorials/incompressible/pimpleFoam/bubble_generation/system/fvSolution from line 55 to line 132.
From function dictionary::subDict(const word& keyword) const
in file db/dictionary/dictionary.C at line 648.
What is the problem?
Akidess. I notice the next error
PIMPLE: Operating solver in PISO mode
Reading field p_rgh
--> FOAM FATAL IO ERROR:
cannot find file
file: /home/paramesh/OpenFOAM/paramesh-3.0.0/run/tutorials/multiphase/interFoam/bubble_generation_trial/0/p_rgh at line 0.
From function regIOobject::readStream()
in file db/regIOobject/regIOobjectRead.C at line 73.
You are missing the file p_rgh. Please have a look at the interFoam tutorials.
hey.. do you still have these files?
I wanted to download them and i am unable to get the files from this link
Problem solved, able to download now.
|All times are GMT -4. The time now is 23:26.|