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

Waves in simple channel

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 17, 2009, 07:59
Default Waves in simple channel
  #1
New Member
 
Join Date: Jul 2009
Posts: 4
Rep Power: 17
taranov is on a distinguished road
Hello, forum!
I'm starting to learn the OpenFOAM and have a newbie question.
I would like to simulate the flow in very simple geometry:
a 2-phase parallelepiped with air-water interface. Without walls and bodies.
The OF version 1.5, the solver - interFoam. maxCo=0.3

The problem is as follows: at leading edge of free surface (at inlet)
raises the high value of velocity, which lead to raising of waves.
See attached pictures. Please advise how to avoid this problem.

Boundary conditions:

gamma:
in_air uniform 0;
in_water uniform 1;
out_water zeroGradient;
out_air zeroGradient;
top slip;
bottom slip;
side slip;

pd:
in_air zeroGradient;
in_water zeroGradient;
out_water zeroGradient;
out_air zeroGradient;
top slip;
bottom slip;
side slip;

U:
in_air uniform (1 0 0);
in_water uniform (1 0 0);
out_water zeroGradient;
out_air zeroGradient;
top slip;
bottom slip;
side slip;

================================================
fvSolvers:
solvers
{
pcorr PCG
{
preconditioner DIC;
tolerance 1e-10;
relTol 0;
};
pd PCG
{
preconditioner DIC;
tolerance 1e-7;
relTol 0.05;
};
pdFinal PCG
{
preconditioner DIC;
tolerance 1e-7;
relTol 0;
};
k PBiCG
{
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
};
omega PBiCG
{
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
};
U PBiCG
{
preconditioner DILU;
tolerance 1e-06;
relTol 0;
};
}

PISO
{
momentumPredictor no;
nCorrectors 3;
nNonOrthogonalCorrectors 0;
nGammaCorr 1;
nGammaSubCycles 2;
cGamma 1;

pdRefCell 1;
pRefCell 1;

pRefValue 1;
pdRefValue 0;
}
===============================================
fvSchemes:
ddtSchemes
{
default Euler;
}

gradSchemes
{
default Gauss linear;
grad(U) Gauss linear;
grad(gamma) Gauss linear;
}

divSchemes
{
div(rho*phi,U) Gauss limitedLinearV 1;
div(phi,gamma) Gauss vanLeer;
div(phirb,gamma) Gauss interfaceCompression;
div(phi,omega) Gauss upwind;
div(phi,k) Gauss upwind;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
pd;
pcorr;
gamma;
}

Thanks in advance
Andrey
Attached Images
File Type: png forum1.png (41.0 KB, 27 views)
File Type: png forum2.png (52.7 KB, 29 views)
taranov is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Why wallGradU result is different from dudy in the smooth channel bed xiuying OpenFOAM Running, Solving & CFD 12 March 2, 2009 08:58
simple airflow through a box channel Leo Lau FLUENT 4 October 14, 2007 06:23
SIMPLE Scheme for open channel flow ao Main CFD Forum 0 April 6, 2004 17:21
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 10:11
SIMPLE algorithm Jonathan Castro Main CFD Forum 3 December 10, 1999 05:59


All times are GMT -4. The time now is 06:18.