Question #683336 on Yade changed: https://answers.launchpad.net/yade/+question/683336
Anqi H posted a new comment: Hi Deepak, Thank you for the links. I've made deltaT 1/10 of the previous value, It now stalls at 0.01468s while both courant and max courant number is less than 0.5. Time = 0.01478 Courant Number mean: 0.00188364 max: 0.00654247 smoothSolver: Solving for Ux, Initial residual = 0.000895276, Final residual = 3.99666e-07, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.0012464, Final residual = 5.33346e-07, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.00136997, Final residual = 3.99407e-07, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.000356892, Final residual = 1.69725e-05, No Iterations 9 DICPCG: Solving for p, Initial residual = 1.69772e-05, Final residual = 7.49094e-07, No Iterations 10 DICPCG: Solving for p, Initial residual = 7.49094e-07, Final residual = 7.49094e-07, No Iterations 0 time step continuity errors : sum local = 6.37045e-12, global = -2.646e-14, cumulative = 1.25224e-08 DICPCG: Solving for p, Initial residual = 2.68288e-06, Final residual = 6.77363e-07, No Iterations 2 DICPCG: Solving for p, Initial residual = 6.77362e-07, Final residual = 6.77362e-07, No Iterations 0 DICPCG: Solving for p, Initial residual = 6.77362e-07, Final residual = 6.77362e-07, No Iterations 0 time step continuity errors : sum local = 5.76044e-12, global = -8.4556e-13, cumulative = 1.25216e-08 ExecutionTime = 128.46 s ClockTime = 129 s ControlDict: /*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: dev \\/ M anipulation | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object controlDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // application icoFoam; startFrom latestTime; startTime 0; stopAt endTime; endTime 0.05; deltaT 0.000005; writeControl timeStep; writeInterval 20; purgeWrite 0; writeFormat ascii; writePrecision 6; writeCompression off; timeFormat general; timePrecision 6; runTimeModifiable true; // ************************************************************************* // fv solution: /*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: dev \\/ M anipulation | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { p { solver PCG; preconditioner DIC; tolerance 1e-06; relTol 0.05; } pFinal { $p; relTol 0; } U { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-05; relTol 0; } } PISO { nCorrectors 2; nNonOrthogonalCorrectors 2; } // ************************************************************************* // fv scheme: /*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: dev \\/ M anipulation | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default Euler; } gradSchemes { default Gauss linear; } divSchemes { default none; div(phi,U) Gauss limitedLinearV 1; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } // ************************************************************************* // 0/p /*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: dev \\/ M anipulation | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -2 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type fixedValue; value uniform .2; } outlet { type fixedValue; value uniform 0; } walls { type zeroGradient; } } // ************************************************************************* // 0/U /*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: dev \\/ M anipulation | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volVectorField; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (1 0 0); boundaryField { inlet { type zeroGradient; } outlet { type zeroGradient; } walls { type fixedValue; value uniform (0 0 0); } } // ************************************************************************* // 0/uSource /*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 6 \\/ M anipulation | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volVectorField; object uSource; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -2 0 0 0 0]; internalField uniform (0 0 0); boundaryField { inlet { type calculated; value uniform (0 0 0); } outlet { type calculated; value uniform (0 0 0); } walls { type fixedValue; value uniform (0 0 0); } } // ************************************************************************* // blockMeshDict: /*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: dev \\/ M anipulation | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object blockMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // convertToMeters 0.001; vertices ( (-2 10 -2) (22 10 -2) (22 22 -2) (-2 22 -2) (-2 10 22) (22 10 22) (22 22 22) (-2 22 22) ); blocks ( hex (0 1 2 3 4 5 6 7) (20 4 20) simpleGrading (1 1 1) ); edges ( ); boundary ( inlet { type patch; faces ( (0 3 7 4) ); } outlet { type patch; faces ( (1 2 6 5) ); } walls { type wall; faces ( (0 1 2 3) (0 1 5 4) (4 5 6 7) (3 2 6 7) ); } ); // ************************************************************************* // -- You received this question notification because your team yade-users is an answer contact for Yade. _______________________________________________ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp