# InterFoam MULES solver

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 December 4, 2007, 12:40 Dear OpenFOAM uses and Develop #1 Senior Member   Join Date: Mar 2009 Posts: 248 Rep 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 --> 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 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

November 13, 2012, 05:36
#2
New Member

Join Date: May 2011
Posts: 15
Rep Power: 6
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

 November 13, 2012, 06:45 #3 Senior Member     Anton Kidess Join Date: May 2009 Location: Delft, Netherlands Posts: 919 Rep Power: 17 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

 November 13, 2012, 07:49 #4 New Member   Join Date: May 2011 Posts: 15 Rep Power: 6 yeah, I got it. it was because I should have written geometricOnesField() in place of OneField(). Thanks !

 November 21, 2012, 09:56 #5 Senior Member     Santiago Marquez Damian Join Date: Aug 2009 Location: Santa Fe, Santa Fe, Argentina Posts: 418 Rep Power: 14 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. 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

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post zou_mo OpenFOAM Running, Solving & CFD 127 May 25, 2011 16:30 feijooos OpenFOAM Running, Solving & CFD 4 December 8, 2008 12:01 m9819348 OpenFOAM Running, Solving & CFD 1 September 20, 2007 13:10 qiu OpenFOAM Running, Solving & CFD 0 May 6, 2007 22:48 mer OpenFOAM Running, Solving & CFD 5 May 31, 2006 12:22

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

 Contact Us - CFD Online - Top