CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

increasing volume fraction in closed domain with compressibleInterFoam

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

LinkBack Thread Tools Display Modes
Old   May 14, 2014, 12:58
Smile increasing volume fraction in closed domain with compressibleInterFoam
Victor-S. Eberhart
Join Date: Oct 2011
Posts: 35
Rep Power: 7
vigges is on a distinguished road
Hi Foamers!

I managed to extend the compressibleInterFoam solver by a time varying acceleration which, in my case, simulates a periodic motion and is provided by a file containing the time-acceleration tuples.

As you can see in the excerpt of my accFoam.C, I only introduced the acceleration in z-axis direction and defined a vector called acc which also contains the gravitation:
    while (
        #include "readTimeControls.H"
        #include "CourantNo.H"
        #include "setDeltaT.H"


        // --- time varying acceleration
        interpolationTable<scalar> lookup("acceleration.dat");

        scalar  acc_X = 0.0;
        scalar  acc_Y = 0.0;
        scalar  acc_Z = lookup(runTime.value());//timeSeriesAcceleration;
        Info<< "Acceleration in z-direction = " << acc_Z << nl << endl;

        const dimensionedScalar gunits("gunits", dimensionSet(0,1,-2,0,0,0,0), 1.0);
        dimensionedVector acc=g + vector(acc_X,acc_Y,acc_Z)*gunits;

        Info<< "Time = " << runTime.timeName() << nl << endl;

        // --- Pressure-velocity PIMPLE corrector loop
        while (pimple.loop())
In order to take the acceleration into account I added the acceleration multiplied by my mixture density to the UEqn.H file:
    fvVectorMatrix UEqn
        fvm::ddt(rho, U)
      + fvm::div(rhoPhi, U)
      + turbulence->divDevRhoReff(U)
      - rho*acc
Also, I extended the energy equation in the TEqn.H file as follows:
   fvScalarMatrix TEqn
        fvm::ddt(rho, T)
      + fvm::div(rhoPhi, T)
      - fvm::laplacian(twoPhaseProperties.alphaEff(turbulence->mut()), T)
      + (
            fvc::div(fvc::absolute(phi, U), p)
          + fvc::ddt(rho, K) + fvc::div(rhoPhi, K)
          - rho * (acc & U)
My test case, a closed test tube filled with 50% liquid sodium and 50% air "shaken" along the longitudinal axis according to a prescribed oscillating motion, yields good results for the liquid phase kinetic energy compared to a CFX case. However, during the simulation, the total liquid phase volume fraction defined by the initial filling level (0.5) is constantly increasing and reached a value of 0.58 so far. I have a suspicion that it has something to do with the compressibility of my fluids.

Speaking of which, I have doubts about the assigned value of R in my thermophysical properties for sodium. I had no idea how to calculate it, so I looked at the thermophysicalProperties.mercury from the compressibleMultiphaseInterFoam damBreak4Phase tutorial case and took a similar value, somy thermophysicalProperties.sodium looks like that:
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectFluid;
    specie          specie;
    energy          sensibleInternalEnergy;

        nMoles          1;
        molWeight       22.99;
        R               5000;
        rho0            852;
        Cp              1277;
        Hf              0;
        mu              0.000264;
        Pr              0.00496;
I scanned the forum but only found threads related to decreasing alpha values. Maybe someone of you has already run into that problem before and found a solution for it.

Thanks a lot,
vigges is offline   Reply With Quote


acceleration, compressibleinterfoam, multiphase, volume fraction

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
volume fraction in mixture model and VOF???? multiphase-flow FLUENT 4 August 7, 2014 10:35
dynamic Mesh is faster than MRF???? sharonyue OpenFOAM Running, Solving & CFD 14 August 26, 2013 07:47
Mass fraction and volume fraction eric weddle CFX 0 September 26, 2011 05:02
On the damBreak4phaseFine cases paean OpenFOAM Running, Solving & CFD 0 November 14, 2008 22:14
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 12:55

All times are GMT -4. The time now is 19:00.