
[Sponsors] 
increasing volume fraction in closed domain with compressibleInterFoam 

LinkBack  Thread Tools  Display Modes 
May 14, 2014, 12:58 
increasing volume fraction in closed domain with compressibleInterFoam

#1 
Member
VictorS. Eberhart
Join Date: Oct 2011
Posts: 33
Rep Power: 6 
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 timeacceleration tuples. As you can see in the excerpt of my accFoam.C, I only introduced the acceleration in zaxis direction and defined a vector called acc which also contains the gravitation: Code:
while (runTime.run()) { #include "readTimeControls.H" #include "CourantNo.H" #include "setDeltaT.H" runTime++; //  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 zdirection = " << 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; //  Pressurevelocity PIMPLE corrector loop while (pimple.loop()) Code:
fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(rhoPhi, U) + turbulence>divDevRhoReff(U)  rho*acc ); Code:
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) ) 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: Code:
thermoType { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState perfectFluid; specie specie; energy sensibleInternalEnergy; } mixture { specie { nMoles 1; molWeight 22.99; } equationOfState { R 5000; rho0 852; } thermodynamics { Cp 1277; Hf 0; } transport { mu 0.000264; Pr 0.00496; } } Thanks a lot, Victor 

Tags 
acceleration, compressibleinterfoam, multiphase, volume fraction 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
volume fraction in mixture model and VOF????  multiphaseflow  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 