|December 4, 2007, 12:40||
Dear OpenFOAM uses and Develop
Join Date: Mar 2009
Posts: 248Rep Power: 9
Dear OpenFOAM uses and Developers
I request the forum users to shed some light on the details of how MULES works. Please share the information if anybody is aware of it.
I looked into the code and this is what it contains:
1) overloaded function templates for
--> explicitSolve - bounds can be supplied as the function argument;
--> explicitSolve01 - for a scalar bounded between 0 and 1;
--> implicitSolve01 - for a scalar bounded between 0 and 1;
and a function template for
InterFoam uses the MULES solver to solve for the volScalarField gamma and calls
MULES::explicitSolve01(gamma, phi, phiGamma).
This eventually calls the function template
const surfaceScalarField& phi,
const SpType& Sp,
const SuType& Su,
const scalar psiMax,
const scalar psiMin
with Sp = zero
Su = zero
psiMax = 1.0
psiMin = 0.0
Could anybody please comment what does template arguments: class SpType, class SuType stand for?
Thanks in advance
|November 13, 2012, 05:36||
Join Date: May 2011
Posts: 15Rep Power: 6
I think they are representing eventual source terms in the advection equation. There is even this possibility :
void Foam::MULES::explicitSolve ( const RhoType& rho, volScalarField& psi, const surfaceScalarField& phi, surfaceScalarField& phiPsi, const SpType& Sp, const SuType& Su, const scalar psiMax, const scalar psiMin )
How can I call this special MULES function (I need to provide both source terms and rho) ?
I do :
|November 13, 2012, 06:45||
Join Date: May 2009
Location: Delft, Netherlands
Posts: 1,019Rep Power: 18
Did you look at interPhaseChangeFoam?
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
*Join the OpenFOAM stackexchange Q&A site: http://area51.stackexchange.com/prop...oHPxcPqde7HtA2
|November 21, 2012, 09:56||
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 420Rep Power: 15
Hi Jaswinder, how you doing? MULES is an Flux Corrected Transport explicit solver for hyperbolic equations. Its implementation mostly follow the ideas of Zalesak limiter but the 's are calculated iteratively. I'm going to give some details in my thesis soon but that's the basic idea.
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Research Center for Computational Mechanics (CIMEC) - CONICET/FICH-UNL
T.E.: 54-342-4511594 Ext. 1005
Güemes 3450 - (3000) Santa Fe
Santa Fe - Argentina
|Thread||Thread Starter||Forum||Replies||Last Post|
|About interFoam solver||zou_mo||OpenFOAM Running, Solving & CFD||127||May 25, 2011 16:30|
|Wmake problem interFoam solver||feijooos||OpenFOAM Running, Solving & CFD||4||December 8, 2008 12:01|
|DICPCG solver in interFoam||m9819348||OpenFOAM Running, Solving & CFD||1||September 20, 2007 13:10|
|About interfoam solver||qiu||OpenFOAM Running, Solving & CFD||0||May 6, 2007 22:48|
|Need documentation for interFOAM solver||mer||OpenFOAM Running, Solving & CFD||5||May 31, 2006 12:22|