Question #683336 on Yade changed: https://answers.launchpad.net/yade/+question/683336
Deepak proposed the following answer: hi, can you send the yade script as well? On Fri, Sep 13, 2019 at 3:57 PM Anqi H <question683...@answers.launchpad.net> wrote: > 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 > -- 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