CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

InterFoam MULES solver

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

Like Tree1Likes
  • 1 Post By santiagomarquezd

Reply
 
LinkBack Thread Tools Display Modes
Old   December 4, 2007, 12:40
Default Dear OpenFOAM uses and Develop
  #1
Senior Member
 
Join Date: Mar 2009
Posts: 248
Rep Power: 9
jaswi is on a distinguished road
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

--> limiter

InterFoam uses the MULES solver to solve for the volScalarField gamma and calls

MULES::explicitSolve01(gamma, phi, phiGamma).

This eventually calls the function template

template<class>
void Foam::MULES::explicitSolve
(
volScalarField& psi,
const surfaceScalarField& phi,
surfaceScalarField& phiPsi,
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

Jaswinder
jaswi is offline   Reply With Quote

Old   November 13, 2012, 05:36
Default
  #2
New Member
 
Join Date: May 2011
Posts: 15
Rep Power: 6
joris.hey is on a distinguished road
I think they are representing eventual source terms in the advection equation. There is even this possibility :
Code:
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  )
with field rho.

How can I call this special MULES function (I need to provide both source terms and rho) ?

I do :

Quote:
MULES::explicitSolve (oneField (), alpha1 , phi , phiAlpha , volScalarField ("Sp",fvc::div(voidfraction*U)),zeroField (), 1, 0);
But it doesn't compile... What's wrong ?

Cheers,
J
joris.hey is offline   Reply With Quote

Old   November 13, 2012, 06:45
Default
  #3
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 919
Rep Power: 17
akidess will become famous soon enough
Did you look at interPhaseChangeFoam?

See: http://foam.sourceforge.net/docs/cpp/a02688_source.html
__________________
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
*Help define the OpenFOAM stackexchange Q&A site: http://area51.stackexchange.com/prop...oam-technology
akidess is offline   Reply With Quote

Old   November 13, 2012, 07:49
Default
  #4
New Member
 
Join Date: May 2011
Posts: 15
Rep Power: 6
joris.hey is on a distinguished road
yeah, I got it. it was because I should have written geometricOnesField() in place of OneField().

Thanks !
joris.hey is offline   Reply With Quote

Old   November 21, 2012, 09:56
Default
  #5
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 418
Rep Power: 14
santiagomarquezd will become famous soon enough
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 \lambda's are calculated iteratively. I'm going to give some details in my thesis soon but that's the basic idea.

Regards.
mprinkey likes this.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Post-doctoral Fellow
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
http://www.cimec.org.ar
santiagomarquezd is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
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


All times are GMT -4. The time now is 15:30.