[Juzzvy]

Good pm, I was running my case in OpenFOAM and now I am at the snappyHexMesh part. When I looked at my mesh, it didn't mesh anything, only the bottom part (see picture). What could be the problem? is it in the locationMesh or is my cells too low? Thank you

[aow]

You need to give more information to work from than just a picture... Your snappyHexMeshDict file and any other relevant meshing info like the blockMeshDict, .stl or other surface file you are using, etc... would probably help a lot.

[Juzzvy]

Good day. Thanks for the concern. I actually meshed it correctly, the locationInMesh seem to do the trick. Now I have another question, I am now using scalarTransport in my controlDict file in order to produce pollutant-like sources. I already tested the run and it is working. My question now is, how can I make the source emit pollutant at a specific time? Thank you

controlDict
application simpleFoam;

startFrom latestTime;

startTime 0;

stopAt endTime;

endTime 1000;

deltaT 1;

writeControl timeStep;

writeInterval 50;

purgeWrite 0;

writeFormat binary;

writePrecision 6;

writeCompression uncompressed;

timeFormat general;

timePrecision 6;

runTimeModifiable true;


libs
(
"libsolverFunctionObjects.so"
);
functions
{
s
{
type scalarTransport;
outputControl outputTime;
active true;
resetOnStartUp false;
autoSchemes true;
nCorr 0;
fvOptions
{
s
{
type scalarSemiImplicitSource;
active true;
timeStart 500;
duration 1000;

scalarSemiImplicitSourceCoeffs
{
selectionMode points;
points
(
(170 0 40)
(70 0 17)
(-30 0 -10)
(-60 0 -14)
(-90 0 -17)
(-120 0 -20)
(-170 0 -27)
(-170 0 -55)
(-110 0 -42)
(-80 0 -37)
(-50 0 -32)
(-20 0 -29)
(80 0 -17)
(180 0 11)
);
volumeMode absolute;
injectionRateSuSp
{
s (10 0);
}
}
}
};

}
}

[aow]

Hey Juzzvy,

How are you planning on injecting this pollutant into your domain? If you know exactly the value of the scalar you want at all times throughout your simulation, then I would recommend trying to find something like the timeVaryingMappedFixedValue boundary condition, which allows you to specify values for a field variable at any number of specific times. OpenFOAM does the rest of the work for you by interpolating between the points you specify to make sure the variable increases/decreases from one specified values to the next over the time intervals you also specified. Setting this up involves adding some stuff to the constant directory of the case directory, including time directories that contains the file specifying the variable values as well as a points file where you can specify the values on the face you're applying the timeVaryingMappedFixedValue boundary condition. I think I learned about it in the pitzDailyExptInlet tutorial case.

I hope this is the sort of thing you're looking for.

Best,

Andrew