← Back to team overview

yade-users team mailing list archive

Re: [Question #683336]: Yade-OpenFoam-coupling delete the particles outside the fluid cells

 

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 <question683336@xxxxxxxxxxxxxxxxxxxxx>
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@xxxxxxxxxxxxxxxxxxx
> 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.