Transport Equation in InterFoam using MULES

 Register Blogs Members List Search Today's Posts Mark Forums Read

 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(); However, this is affected by a large amount of numerical diffusion. 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 phiCure ( fvc::flux ( phi, c, alphaScheme ) + fvc::flux ( -fvc::flux(-phir, alpha2, alpharScheme), c, alpharScheme ) ); together with Code: MULES::explicitSolve(c, phiCN, phiCure, 2, 0); but this does not compile with wmake: error: no matching function for call to ‘explicitSolve(Foam::volScalarField&, Foam::tmp >&, Foam::tmp >&, 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
Quote:
 Originally Posted by Cyp 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,
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