
[Sponsors] 
August 1, 2017, 16:41 
Transport Equation in InterFoam using MULES

#1 
New Member
Thomas C. Sykes
Join Date: Jul 2017
Location: University of Leeds, UK
Posts: 11
Rep Power: 6 
Hi All,
OpenFOAM 4.1 I am adding a scalar transport equation to the InterFOAM solver $\frac{\partial c}{\partial t} + \nabla\cdot(c{\bf u})$ to track the progress of a scalar concentration field (c) which I define in setFields to be 0,1,2 initially in certain locations. I have successfully added a standard equation in the InterFoam.C file: Code:
fvScalarMatrix CEqn ( fvm::ddt(c) + fvm::div(phi, c) ); CEqn.solve(); I understand I need to use the MULES FCT scheme to try and solve this, as noted in Using MULES to solve transport equation with sharp interface but am unsure of how to do this. I've tried adding the following to the alphaEqn file: Code:
tmp<surfaceScalarField> phiCure ( fvc::flux ( phi, c, alphaScheme ) + fvc::flux ( fvc::flux(phir, alpha2, alpharScheme), c, alpharScheme ) ); Code:
MULES::explicitSolve(c, phiCN, phiCure, 2, 0); error: no matching function for call to ‘explicitSolve(Foam::volScalarField&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> >&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> >&, int, int)’ among many errors. Any help appreciated! Cheers. 

August 7, 2017, 16:18 

#2 
New Member
Thomas C. Sykes
Join Date: Jul 2017
Location: University of Leeds, UK
Posts: 11
Rep Power: 6 
Any help appreciated.


August 7, 2017, 18:48 

#3 
Senior Member
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 299
Rep Power: 15 
If you just advect a passive scalar without particular condition at the interfaces, why not just solve C from alpha? In your case you just have C=alpha.
Cheers, 

August 10, 2017, 22:12 

#4 
New Member
Thomas C. Sykes
Join Date: Jul 2017
Location: University of Leeds, UK
Posts: 11
Rep Power: 6 
Hi, I've got there regions (2 droplets coalescing) of flaid. Two have the same alpha value and the issue is I need to know where one droplet is going when it coalescea with the other. So alpha gets me the outline of the two droplets together but not of each separately if you see what I mean. Because the two droplets are (currently) the same fluid. Thanks for the suggestion though


August 15, 2017, 23:08 

#5 
Senior Member
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 299
Rep Power: 15 
When two droplets coalesce they become one... I am not sure your problem is well posed.


March 20, 2018, 11:57 

#6 
New Member
Marius Bünker
Join Date: Oct 2017
Posts: 13
Rep Power: 5 
I am stuck with a similar problem. Have you gotten any further?
Edit: Not a similar problem though. You were trying to compile two int, when the compiler wanted to see fields... 

Tags 
interfoam, mules, programming, transport equation 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Setting the height of the stream in the free channel  kevinmccartin  CFX  10  July 9, 2015 21:36 
Solving Transport equation for additional variable using CFX  creddy_trddc  CFX  11  April 2, 2015 21:07 
Calculation of the Governing Equations  Mihail  CFX  7  September 7, 2014 06:27 
Turbulent Heat Transfer Transport Equation  Flo.duck  Main CFD Forum  0  May 6, 2009 03:37 
Poisson equation combined with transport equation  SeHee  CFX  0  December 27, 2007 01:00 