[ChukerSweet]

Hi guys!
I need to impose a velocity source term, that simulates a pump, in a closed loop. The fluid is turbulent and incompressible, interested both in stationary and transitory condition.
Could someone indicate me a similar case study so I can extrapolate how to impose such condition?
THX in advance.

[cnsidero]

Originally Posted by ChukerSweet

Hi guys!
I need to impose a velocity source term, that simulates a pump, in a closed loop. The fluid is turbulent and incompressible, interested both in stationary and transitory condition.
Could someone indicate me a similar case study so I can extrapolate how to impose such condition?
THX in advance.

In OF >= 2.2, there is a pressureGradientExplicitSource source available through fvOptions. This allows one to specify a target velocity and the source term will adjust the pressure gradient to achieve the target velocity. See $FOAM_SRC/fvOptions/sources/derived/pressureGradientExplicitSource for all the details. It can even be used for specific cell zones. This however, is a steady-state value.

If you need a time-varying source, there is - to my knowledge - none available. Your options would be to write a new fvOptions source for a time varying velocity or use swak4Foam which gives the ability to code in generic source term into your solver and one specifies the form of the source term in a dictionary.

The nice thing about the swak4Foam approach is once the generic source term is hooked into the solver, the mathematical form is specified through the dictionary. Therefore, if one wants to modify either the form or the parameters, no compilation is necessary. So one can easily experiment with the source term and the parameters. In fact, if one was so inclined, a parametric study of the source would be straightforward with a script.

[ChukerSweet]

Thank you very much for the reply. The fvOption will work for now.
Eventually I will need to use a variable pump condition but for the time being, and benchamrking porpoises a static one will do just fine.
Again THX allot

[ChukerSweet]

Just one little more. I can't seem able to "activate" the option...

One could write, from the header suggestion:
pressureGradientExplicitSourceCoeffs
{
fieldNames (U); // name of velocity field
Ubar (10.0 0 0); // desired average velocity
gradPini gradPini [0 2 -2 0 0] 0; // initial pressure gradient
flowDir (1 0 0); // flow direction
}

Or from the PimpleFoam/Channel395 example (found it ! )

momentumSource
{
type pressureGradientExplicitSource;
active on; //on/off switch
selectionMode all; //cellSet // points //cellZone

pressureGradientExplicitSourceCoeffs
{
fieldNames (U);
Ubar ( 0.1335 0 0 );
}
}
In which the pressure is applied all over the place, I think. Anyhow whter I put one or the other file in, let's say, the Cavity example (in the systemo folder) nothing happens. This is logical, the control files need to know that I'm trying to impose a velocity somewere, but how do I do this?

Sorry, and THX

[cnsidero]

Originally Posted by ChukerSweet

In which the pressure is applied all over the place, I think.
Sorry, and THX

Yes but as I mentioned and as indicated in the sample fvOptions from the channel395 tutorial, one can apply to specific regions.

momentumSource
{
    type            pressureGradientExplicitSource;
    active          on;        // on/off switch
    selectionMode   all;       // cellSet // points //cellZone <= HERE

    pressureGradientExplicitSourceCoeffs
    {
        fieldNames  (U);
        Ubar        ( 0.1335 0 0 );
    }
}

You would replace 'all' with one of cellSet/points/cellZone and then specify the name of your set, e.g.

selectionMode cellSet;
cellZone rotor_zone;

and you would have to had previously created a cell zone called 'rotor_zone'. If your mesh generator doesn't support creating zones, use the setSet utility. The details of how fvOptions works on sub-regions is found in $FOAM_SRC/fvOptions/fvOptions/fvOption.C

[ChukerSweet]

THX allot for the info!
Yes, the fvOption should be accompanied with the topoSet function. Looking into it. The multiphase-> throttle exampe seams to implement it, though not combined with the fvOption