CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM (
-   -   bubbleInterTrackFoam and Mass transport (

kel85uk July 12, 2011 00:10

bubbleInterTrackFoam and Mass transport
2 Attachment(s)
Hi all,

I am trying to model mass transfer across a bubble interface using bubbleInterTrackFoam, but when I added the scalar transport equation into the solver, it does not solve anything.

The following is an excerpt of the solver. Any help and suggestions on how to make this work? Thank you



#              include "continuityErrs.H"

                // Momentum corrector
                U -= fvc::grad(p)/AU;

//#          include "solveBulkSurfactant.H"
#          include "species_transfer.H"


#          include "freeSurfaceContinuityErrs.H"

#      include "updateMovingReferenceFrame.H"

#      include "volContinuity.H"


Info << "Solving the species field" << endl;

DT = fluidIndicator*( DTA - DTB ) + DTB;

  + fvm::div(phi - fvc::meshPhi(rho,U), C_O2)
  //+ fvm::div(phi, C_O2)
  == fvm::laplacian(DT, C_O2)

Bernhard July 12, 2011 03:47

I have the impression that you are now solving the convection-diffusion equation in two seperate domains (you can easily check that I think with a simple case).
For the mass transfer across the interface, you should somehow communicate between the patches freeSurface(Shadow).

hfsf December 24, 2012 08:25

bubbleInterTrackFoam + Mass Transport

I would like to know if you could share your bubbleInterTrackFam + Mass transport implementation with me. I'm working with the same topíc, a falling drop with the conjugated solute transport to the exterior, but I'm having trouble even setting the case. Any help is appreciated :D.

Best regards,

Hanniel Freitas

Eli1993 May 11, 2016 07:48

Hello everyone
I am using bubbleInterTrackFoam to simulate a bubble rising in a liquid column.I got some error like bellow

Create time

Create dynamic mesh for time = 0

Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: laplace
Selecting motion diffusivity: uniform

Reading field p

Reading field U

Reading/calculating face flux field phi

Found free surface patch. ID: 1
Found free surface shadow patch. ID: 0

Error in B-to-A face patchToPatchInterpolation.

From function freeSurface::makeInterpolators()
in file makeFreeSurfaceData.C at line 105.

FOAM aborting

Can anyone pls tell me what to do:(

All times are GMT -4. The time now is 05:04.