[Alasir]

Hi.

Im trying to use directionalPressureGradientExplicitSource as a momentumsource to get my channel going.

A problem with this is that it requires a Facezone set. Ive aquired one using the topoSet application, but I seem to have encountered a problem.

I did the orgiinal mesh using gmsh, with gave me a "Internal" cell list.

toposetDict:

actions
(
    {
        name    fs;
        type    faceSet;
        action  new;
        source  cellToFace;
        sourceInfo
        {
        set Internal; option both;
        }
    }
    
    {
        name    fz;
        type    faceZoneSet;
        action  new;
        source  setToFaceZone;

        sourceInfo {faceSet fs;}
    }
);

fvOptions:

Creating finite volume options from "constant/fvOptions"

Selecting finite volume options model type directionalPressureGradientExplicitSource
    Source: fluidDeflection
    - selecting all cells
    - selected 249411 cell(s) with volume 1.01957e-10
    Initial pressure gradient = 249411{(0 0 0)}


Starting time loop

forces forceCoeffs:
    Not including porosity effects

forceCoeffs forceCoeffs:
    Not including porosity effects

Time = 1e-14

Courant Number mean: 0 max: 0
smoothSolver:  Solving for Ux, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Uz, Initial residual = 0, Final residual = 0, No Iterations 0
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib64/libc.so.6"
#3  Foam::fv::directionalPressureGradientExplicitSource::correct(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) at ??:?
#4  ? at ??:?
#5  ? at ??:?
#6  __libc_start_main in "/lib64/libc.so.6"
#7  ? at ??:?
Floating point exception
[ofuser@default MasterRe1000003x4kEq]$

Note that pisoFoam runs if I dont use fvOptions.

Any help on how to use the correct facezones would be appreciated.