Hi,
I was wondering if anyone had an example of using the twoLiquidMixingFoam solver that they would be willing to share? The user and programer's reference don't seem to have one. I am interested in modeling some microfluidic mixers.

Hi Hrv, (and others)
I went to run the example you sent and got the following error:

--> FOAM FATAL ERROR : gradientInternalCoeffs cannot be called for a defaultFvPatchField.
(Actual type fixedMeanValue)
You are probably trying to solve for a field with a default boundary conditions.

From function defaultFvPatchField<type>::gradientInternalCoeffs( ) const
in file fields/fvPatchFields/basicFvPatchFields/default/defaultFvPatchField.C at line 421.

FOAM exiting

I noticed that when defining the boundary conditions you used an auxilary file outside of the blockMesh dictionary. The file boundary in the polyMesh/ directory. I'm a little new to OpenFOAM and am not entirely sure how this works. If somone has the time to clue me in I would greatly appreciate it. All the examples in the programing guide appear to have b.c.'s self contained in the blockMeshDictonary. Putting them in a seperate file has obvious advantages. I suspect my problem is that I have not figured out how to tell interFoam to look for the boundary file.

Hi,
I still a little confused. Is defaultFvPatchField a boundary condition you wrote? If so where is it referenced that I need to replace with a fixedValue? I didn't see it in the .OpenFOAM-1.3/Fields directory. Or do I need to replace some of the conditions in the blockMeshDictionary with fixedValue? Sorry for all the newbie questions.

Hi,
I was wondering if their is a simple way to make the viscosity in twoLiquidMixing foam a function of gamma. I am interested in mixing two newtonian liquids where the viscosity is a function of the fraction gamma. It seems like this should be fairly straightforward to code, but I wanted to ask the various c++ and OpenFoam masters before I invest the time. Also, if you could point me to the part of the code I need to tinker with that would be great.

Aha, that's one of my new boundary conditions: change it to fixedValue and all will be well.

Hrv

Hi Hrv,

is there a version of fixedMeanValue for Openfoam 1.6 yet ?
I would be really glad to get hold of it, if possible, since i have been dealing with a problem which might be solved by your boundary condition. But I kind of dislike the thought of going back do OF 1.5...or earlier.

Does anyone of you know where I can find some informations about model being implemented in TwoLiquidMixingFoam, especially descriptions of the extra terms in NS equations ?
1. (fvc::grad(U) & fvc::grad(muEff))
2. what function reconstruct does ?

They use flux reconstruction to introduce the pressure gradient and the buoyancy term (combined now with the introduction of p_rgh) in the momentum predictor. Due to this, the pressure gradient introduced as the contribution it gives to the flux, as argument of reconstruct().

__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.OpenQBMM - An open-source implementation of quadrature-based moment methods