[hangs]

Hello,

I want to simulate the flow in a shock tube by using rhoCentralFoam solver with a wedge geometry (axisymmetric, 2D). In order to set the different pressure values for the driven and driver sections, I edited the file "setFieldsDict" and run "setFields" function before I run "rhoCentralFoam" solver.

In "setFieldsDict" file, I used "boxToCell" to identify the cells in driver section (although I do not think it is a best way), so I draw a sufficiently big box (at least I think it was sufficient) to cover all the driver section so that I can set a high pressure. However, I did not work well, in the test results, only the pressure at the side wall of the driver section was the changed, all others were not. So I am writing to look for some tips about how to specify the nonuniform initial conditions.

By the way, I generated the mesh with Gmesh where I specified the physical volume "driver_section" and "driven_section", and I can see them after I run "gmshToFoam" in "constant/polyMesh/cellZones.c" with groups of numbers. So is any way that I can take advantage of these "zone names" for specify initial conditions?

Also, please let me know if I missed any important step. This is my first time using OpenFoam.

Thanks so much for your help!

[usv001]

Hi there,

You probably can use 'zoneToCell' in setFieldsDict instead of 'boxToCell'. Your setFieldsDict should look something like this:

defaultFieldValues
(
    ... // set"driven_section" values
);

regions
(
    zoneToCell
    {
        name    "driver_section";
        fieldValues
        (
            ... // set "driver_section" values
        );
    }
);

USV