[dasith0001]

Hi Formers,

I have this question for some time and I could not yet figure out a way to do it. Here is my question.

1. There is two fluid domains with same thermophysical properties. Each fluid domain can have a inlet and a outlet. (fluid domains should have at lease one external boundary)

2. Two domains are physically connected by a narrow channel.

3. I want to define a pressure gradient between the internal surfaces of the narrow channel. ( so that the narrow channel basically acts as a relief valve).

Could you suggest me an any direction setting up such a case? Or is it even possible to do this in OpenFoam ? Your help is greatly appreciated.

Thank you.
Dasith

[Bloerb]

I do not quite understand what you mean by two domains...do you mean two separate solvers? Or just two physically separate regions, but with a channel between them. Hence fully connected.

You can set a pressure gradient inside your domain using fvOptions...e.g semiImplicitSource or directionalPressureGradientExplicitSource

You might want to look into baffles (zero thickness walls) if you are talking about two domains, with a wall like interface between them. If those domains are somehow separate, but the same solver. Otherwise something like the multiple world coupling allows coupling several solvers (i think this was introduced in v2012)


So do you want to set this pressure gradient as a boundary condition or as a volumetric source term?

[dasith0001]

Hi Bloerb,

Thank you very much for your response.

sorry for the confusion, it is just two physically separated regions but with a channel between them. So you are right, just one fluid region.

I have been abled to follow on the ''directionalPressureGradientExplicitSource'' with fvOptions. created a faceZone to accommodate the pressure gradient.

No, I am not interested in the velocity so I set it up as a constant pressureDrop.

Please see the fvOptions file, the model just stops running at the first time step without giving an error. Please see the following fvOption

PHP Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5                                     |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object fvOptions;
}


reliefValve1
{
    type                 directionalPressureGradientExplicitSource;
    active            on;

    directionalPressureGradientExplicitSourceCoeffs

    {
    selectionMode          all; //cellZone; //all;
    faceZone          valveR;
    fields              (U);
    flowDir              (0 0 1);
    relaxationFactor      0.3;
    model               constant; //DarcyForchheimer;
    
    //constant model
    pressureDrop          1216;    //;

    }
}

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


Is something wrong with the fvOptions file?

I highly appreciate your help, thank you.

Best Regards,
Dasith

[Bloerb]

You should probably not use all but apply the pressure gradient only in the small region. Hence define a cellZone where the pressure gradient should be added.

[dasith0001]

Hi Bloerb,

Thank you very much for the quick reply.
I am still relatively new to OpenFoam and your directions on this highly appreciated.

I tried to set the a cellZone with cellSet but I have not been successful so far.

Please correct me if I am wrong.

I have already define fluid region in topoSet (say regionA).

I want to define a sub region within regionA, say region 'channel'.

In the toposet file I have tried it with 'subset' but I could not define a cellzone named 'channel' with the regionA.

It is much appreciated if you could direct me to a example or for any suggestions. my subset looks like this

PHP Code:

   {
        name    channelCellSet; // chamberTop Section top
        type    cellSet;
        action  new;
        source  cylinderToCell;
        p1      (0 0 0.736); //
        p2      (0 0 0.730); // bottom centre 0.686
        radius   0.028; //   =0.028 - 0.001
        innerRadius 0.001;
    }

    {
        name    channel;
        type    cellZoneSet;
        action  new;
        source  setToCellZone;
        set     channelCellSet;
    }

    {
        name channel;
        type cellSet;;
        action subset;
        source cellToCell;
        sourceInfo
        {
            name "regionA";
        }
    } 


I could not find much literature on this. I tried with couple of other combinations but either topoSet or splitMeshRegions or both don't like it.

Thanks,
Dasith