|
[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: 9 |
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: 9 |
Any help appreciated.
|
|
August 7, 2017, 18:48 |
|
#3 |
Senior Member
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 299
Rep Power: 18 |
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: 9 |
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: 18 |
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: 9 |
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 |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Setting the height of the stream in the free channel | kevinmccartin | CFX | 12 | October 13, 2022 21:43 |
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 | Se-Hee | CFX | 0 | December 27, 2007 01:00 |