[anthonygay1812]

Hello Everyone

I'm running into some errors with dsmcFoam, specifically with properly implementing BC's. I'm working on running a shocktube case but I've created a more simplified example to show my issue.

I'm running a 1D simulation with a zeroGradient BC's with just particles sitting there, but for some reason I seem to be losing particles at the boundaries. I've double checked my implementations of the BC's and just can't figure out whats going on.

I've tried to attach my case but the file is too large, I'm currently working with WSL and git so I think files may not be compressing properly. I went ahead and attached a link to my git hub and a picture of whats going on.

Link to Git: https://github.com/Anthony-Gay/HPC_workspace
Specific Case: https://github.com/Anthony-Gay/HPC_w...dDsmcShockTube

Thank you!

blockMesh File:

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 1;

vertices // places vertices in coordinate plane
(
    (-5 -1 -1)
    (5 -1 -1)
    (5 1 -1)
    (-5 1 -1)
    (-5 -1 1)
    (5 -1 1)
    (5 1 1)
    (-5 1 1)
);

blocks
(
    hex (0 1 2 3 4 5 6 7) (1000 1 1) simpleGrading (1 1 1) // gives labels to vertices and #cells in each dim and there sizes?
);

boundary
(
    sides
    {
        type patch; // patch is generic type for inlet or outlet NOTE: the sides tags correspond to the boundaries specified in the 0 dict
        faces
        (
            (1 2 6 5)
            (0 4 7 3)
        );
    }
    empty
    {
        type empty; // Which faces do not require a solution, compresses to 1D
        faces
        (
            (0 1 5 4)
            (5 6 7 4)
            (3 7 6 2)
            (0 3 2 1)
        );
    }
);

mergePatchPairs
(
);


// ************************************************************************* //

Example 0 File (rhoN):

dimensions      [0 -3 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    sides
    {
        type            zeroGradient;
    }

    empty
    {
        type            empty;
    }
}

// ************************************************************************* //

[anthonygay1812]

I managed to figure out that dsmcFoam really only has an inflow model for patches and that applying a zerogradient to most of the parameters doesnt mean anything. To get the Dirichlet BC's I needed, I switched to a symmetry BC.

If youre running into issues this thread helped me out a lot. 0 Folder for a dsmcFoam simulation