CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

rhoCentralFoam boundary issues with custom local time stepping

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

LinkBack Thread Tools Display Modes
Old   February 20, 2012, 11:15
Default rhoCentralFoam boundary issues with custom local time stepping
Laurens Van Dyck
Join Date: Jul 2011
Location: Netherlands/Germany
Posts: 34
Rep Power: 8
laurensvd is on a distinguished road
I have an axisymmetric problem which I am simulating with rhoCentralFoam and a k-Omega-SST turbulence model. At my outlet I have anomalies (see an example of the pressure distribution in the image below).

The reason for these anomalies is related to me using a local time stepping method which I got from this forum ( I dont have the anomalies in the original code without the local time stepping. Now, I would still like to have the local time stepping since I am only interested in the steady state solution and my calculations take way too long otherwise. Can anyone suggest some changes to get rid of these numerical errors?

My patches are defined as follows (not mentioning my model patches):
patch: frontAndBack
   all wedge

patch: symmetryLine
    all zeroGradient

patch: outlet
    scalar        k        inletOutlet
    scalar        alphat        calculated
    scalar        p        fixedValue
    scalar        T        fixedValue
    scalar        mut        calculated
    scalar        omega        inletOutlet
    vector        U        zeroGradient

patch: atmosphere
    scalar        k        fixedValue
    scalar        alphat        calculated
    scalar        p        slip
    scalar        T        zeroGradient
    scalar        mut        calculated
    scalar        omega        fixedValue
    vector        U        slip

patch: inlet
    scalar        k        fixedValue
    scalar        alphat        calculated
    scalar        p        zeroGradient
    scalar        T        zeroGradient
    scalar        mut        calculated
    scalar        omega        fixedValue
    vector        U        fixedValue
My fvSolution is:
        solver          diagonal;

        solver          smoothSolver;
        smoother        GaussSeidel;
        nSweeps         2;
        tolerance       1e-12;
        relTol          0;


        tolerance       1e-10;
        relTol          0;
        solver GAMG;
        tolerance 1e-08;     
        relTol 0.01;     
        smoother GaussSeidel;     
        cacheAgglomeration off;     
        nCellsInCoarsestLevel 20;     
        agglomerator faceAreaPair;     
         mergeLevels 1;

and my fvSchemes:
fluxScheme      Tadmor;

    default        localEuler rDeltaT;


    default cellLimited leastSquares 1;

    default none;
    div(tauMC) Gauss linear;
    div(phi,k) Gauss Gamma 1;
    div(phi,omega) Gauss Gamma 1;

    default Gauss linear corrected;

    default linear;
    reconstruct(rho) Gamma 1;
    reconstruct(U) GammaV 1;
    reconstruct(T) Gamma 1;

    default corrected;

    default no;    
laurensvd is offline   Reply With Quote


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
High Courant Number @ icoFoam Artex85 OpenFOAM Running, Solving & CFD 11 February 16, 2017 14:40
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32
Problem with FloatingObject Leech OpenFOAM Running, Solving & CFD 10 March 29, 2012 15:24
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 05:05
AMG versus ICCG msrinath80 OpenFOAM Running, Solving & CFD 2 November 7, 2006 16:15

All times are GMT -4. The time now is 14:05.