CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   Pressure convergence problems with directMapped (

Taltan March 14, 2012 08:02

Pressure convergence problems with directMapped
1 Attachment(s)
I am trying to simulate an infinite corrugated pipe using a single corrugation and a wedge. The solver is rhoPimpleFoam with LES Smagorinsky. The mesh is attached. checkMesh reports no errors with

Mesh non-orthogonality Max: 39.808437694 average: 4.9037417816
Max skewness = 0.5954216663 OK.

I used the pitzDailyDirectMapped setup but changed the sampleMode to nearestPatchFace. My entry in the boundary file is

        type            directMappedPatch;
        nFaces          37;
        startFace      8051;
        sampleMode        nearestPatchFace;
        samplePatch        outlet;
        offset                (0.0029 0 0);
        //sampleMode      nearestCell;
        //sampleRegion    region0;
        //samplePatch    none;
        //offsetMode      uniform;

Boundary conditions:

        type            directMapped;
        value          uniform (10 0 0);
        interpolationScheme cell;
        setAverage        true;
        average                (10 0 0);
        type            inletOutlet;
        inletValue      uniform (0 0 0);
        value          uniform (0 0 0);


        type            zeroGradient;
        type            waveTransmissive;
        field          p;
        phi            phi;
        rho            rho;
        psi            psi;
        gamma          1.3;
        fieldInf        1e5;
        lInf            1;
        value          uniform 1e5;

I found that waveTransmissive performs better than fixedValue. The problem is that the pressure never converges properly. When running the outer correctors of the PIMPLE algorithm the pressure converges either at 10 loops or never (50 max). Mostly it does not converge. I think that the pressure and the flow field across the mapped patch is the problem, but I do not know how to fix this.

Any help would be greatly appreciated.

Taltan March 15, 2012 05:44

If anyone is interested in running my case for themselves, I have uploaded my whole case with mesh.

Taltan March 16, 2012 05:56

Did you get that error with my case setup? I am running OF 2.0. The solver I am using for the pressure is GAMG, so that should not give an error.

EDIT: responded to a post here which has since been removed.

klio April 27, 2012 03:59


Im trying to do a multi Region sim and now I have to set Bc. It is not very clear for me what compressible::turbulentHeatFluxTemperature does and how I can change my BC so that OF don't want so cast my bc and how i can define an energy heatsource on my solid patch.

Can you help me with my problem?

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