|
[Sponsors] |
July 17, 2009, 07:59 |
Waves in simple channel
|
#1 |
New Member
Join Date: Jul 2009
Posts: 4
Rep Power: 17 |
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 |
|
|
|
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 |