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

multiphaseEulerFoam loosing fluid with zeroGradient for p on walls???

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 28, 2012, 07:33
Default multiphaseEulerFoam loosing fluid with zeroGradient for p on walls???
  #1
Member
 
Petr Furmanek
Join Date: Jan 2012
Location: Prague, Czech Republic
Posts: 60
Rep Power: 5
petr.f. is on a distinguished road
Dear fellow foamers,

I've run into another strange problem with multiphaseEulerFoam (mpEF). In order to perform scalability tests with mpEF, I've prepared a simple 3d dam break case (only 2 phases - air and water). I've used standard boundary conditions and solver settings taken from tutorial files (see below). When I use the zeroGradient BC for pressure on solid walls all the water somehow leaks through the walls (with no effect on Courant number!). fixedFluxPressure works just well... Any idea where this behaviour comes from?

--
0.org/alphaair

boundaryField
{
walls
{
type zeroGradient;
}
atmosphere
{
type inletOutlet;
phi phiair;
inletValue uniform 1;
value uniform 1;
}
}

--
alphawater is almost the same, only on atmosphere the inletValue and value are uniform 0
--
0.org/p

boundaryField
{
walls
{
type zeroGradient;
}
atmosphere
{
type totalPressure;
p0 uniform 0;
U Uair;
phi phiair;
rho rho;
psi none;
gamma 1;
value uniform 0;
}

}

--
0.0rg/Uair

boundaryField
{
walls
{
type fixedValue;
value uniform (0 0 0);
}
atmosphere
{
type fluxCorrectedVelocity;
value uniform (0 0 0);
}
}
--
... again, Uwater is the same
--
system/fvSolution

solvers
{
p
{
solver GAMG;
tolerance 1e-7;
relTol 0.05;
smoother DIC;//GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration on;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}

pFinal
{
solver PCG;
preconditioner
{
preconditioner GAMG;
tolerance 1e-7;
relTol 0;
nVcycles 2;
smoother DIC;//GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration on;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}
tolerance 1e-7;
relTol 0;
maxIter 20;
}

pcorr
{
$pFinal;
tolerance 1e-5;
relTol 0;
}

U
{
solver smoothSolver;
smoother DIC;//GaussSeidel;
tolerance 1e-8;
relTol 0.1;
nSweeps 1;
}

UFinal
{
$U;
tolerance 1e-7;
relTol 0;
}
}

PIMPLE
{
// nOuterCorrectors 1;
nCorrectors 3;
nNonOrthogonalCorrectors 0;
nAlphaSubCycles 3;
}

relaxationFactors
{

"U.*" 1;
p 1;

}
--
As for fvSchemes, I use vanLeer and limitedLinear as in damBreak tutorial.

Thanks, Petr.
Attached Images
File Type: gif mpEF_damBreak3D.gif (96.8 KB, 12 views)

Last edited by petr.f.; November 28, 2012 at 09:18.
petr.f. 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
Fluid Structure Interaction Apollo Main CFD Forum 5 July 4, 2011 16:15
What is the total energy for incompressible fluid? Harry Dong Main CFD Forum 12 February 4, 2006 01:55
Intl Conf Computational Methods in Fluid Power Jacek Stecki Main CFD Forum 0 November 10, 2002 06:49
My Revised "Time Vs Energy" Article For Review Abhi Main CFD Forum 2 July 9, 2002 09:08
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 09:11


All times are GMT -4. The time now is 01:34.