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

How to solve this alpha Eqn using MULES?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 29, 2013, 14:44
Default How to solve this alpha Eqn using MULES?
  #1
Member
 
David Long
Join Date: May 2012
Location: Germany
Posts: 54
Rep Power: 5
keepfit is on a distinguished road
Hi Foamers,

I want to modify the interFoam solver where the alpha Eqn is like:

Eqn 1:


The fluid Cell contains fluid phase 1 and phase 2, along with the solid phase, here \epsilon_f stands for the volume fraction occupied by fluid (V_phase1 + V_phase2).

\epsilon_f filed is also defined in createFields.H
Code:
    volScalarField voidfractionFluid
    (
        IOobject
        (
            "voidfractionFluid",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE
        ),
        mesh
    );
The original VOF method in interFoam is:

Eqn 2:

the implementation code of alpha equation:

Code:
 for (int aCorr=0; aCorr<nAlphaCorr; aCorr++)
    {
        surfaceScalarField phiAlpha
        (
            fvc::flux
            (
                phi,
                alpha1,
                alphaScheme
            )
          + fvc::flux
            (
                -fvc::flux(-phir, scalar(1) - alpha1, alpharScheme),
                alpha1,
                alpharScheme
            )
        );

        MULES::explicitSolve(alpha1, phi, phiAlpha, 1, 0);

        rhoPhi = phiAlpha*(rho1 - rho2) + phi*rho2;
So I wonder how to add the volScalarField voidfractionFluid to the MULES solver similar like above? Any hints will be highly appreciated!

Cheers,

David Long

Last edited by keepfit; October 2, 2013 at 12:21.
keepfit is offline   Reply With Quote

Old   October 11, 2013, 13:36
Default
  #2
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, I guess it could be done using the source term included in MULES.

Regards.
__________________
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

Old   October 23, 2013, 04:06
Default
  #3
New Member
 
Mahdi S
Join Date: Jul 2012
Posts: 12
Rep Power: 4
Mahdi2010 is on a distinguished road
Hi Santiago
I followed your PhD thesis where you explain how MULE works. Concerning this problem
how the source terms should look like? because MULE needs Su and Sp source terms
to solve it implicitly or explicitly, but the question is how to derive Su and Sp from this source term on the right hand side.
Mahdi2010 is offline   Reply With Quote

Reply

Tags
alphaeqn, mules, source term

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
Using MULES to solve transport equation with sharp interface derkermit OpenFOAM Programming & Development 3 February 27, 2014 07:37
interFoam: Solve gives better results than MULES hrushi.397 OpenFOAM Running, Solving & CFD 6 September 10, 2013 12:57
How to solve energy eqn with enthalpy terms? sabitha FLUENT 0 April 6, 2006 04:58
Alpha to be phased out Heiko Gerhauser Main CFD Forum 17 July 25, 2001 14:40
Determining alpha and beta for porous baffle Liaqat Khan CD-adapco 1 October 27, 2000 04:44


All times are GMT -4. The time now is 05:02.