[Kojote]

Hi,

i am struggeling to add some material via fvOptions. Did anyone here manage to add mass e.g. over a cellZone?

Thanks for helping me on that one.

Br

Christian

[Fauster]

Not sure to understand what you are trying to accomplish. You would like to simulate phase change ?

[Kojote]

Originally Posted by Fauster

Not sure to understand what you are trying to accomplish. You would like to simulate phase change ?

No i want to have a meshless inlet on a certain point.

[Fauster]

You want to have a "boundary condition" on a certain point on your mesh for defining velocity, alpha, temperature etc.. ?

I dont think it's possible to do that on OpenFOAM. fvOptions class allows you to have source terms. For Momentum equation you will be able to define an acceleration source but this is not what you want to do. For energy it will be power source and so on.

[Kojote]

Hi

thanks for the replay. I found code like this -->

Mass source sounds like what i want to do.

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

injector1
{
timeStart 0.1;
duration 5;
selectionMode points;
points
(
(0.075 0.2 0.05)
);
}

options
{
massSource1
{
type scalarSemiImplicitSource;

$injector1;

volumeMode absolute;
injectionRateSuSp
{
thermo:rho.air (1e-3 0); // kg/s
}
}

momentumSource1
{
type vectorSemiImplicitSource;

$injector1;

volumeMode absolute;
injectionRateSuSp
{
U.air ((0 -1e-2 0) 0); // kg*m/s^2
}
}

energySource1
{
type scalarSemiImplicitSource;

$injector1;

volumeMode absolute;
injectionRateSuSp
{
e.air (500 0); // kg*m^2/s^3
}
}
}


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

[BIRAJ]

Hello, for this you have to mention the region in the topoSet and introduce the source term in fvOption mentioning the same region you mentioned in the topoSet.

[dewey]

Originally Posted by Kojote

Hi

thanks for the replay. I found code like this -->

Mass source sounds like what i want to do.

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

injector1
{
timeStart 0.1;
duration 5;
selectionMode points;
points
(
(0.075 0.2 0.05)
);
}

options
{
massSource1
{
type scalarSemiImplicitSource;

$injector1;

volumeMode absolute;
injectionRateSuSp
{
thermo:rho.air (1e-3 0); // kg/s
}
}

momentumSource1
{
type vectorSemiImplicitSource;

$injector1;

volumeMode absolute;
injectionRateSuSp
{
U.air ((0 -1e-2 0) 0); // kg*m/s^2
}
}

energySource1
{
type scalarSemiImplicitSource;

$injector1;

volumeMode absolute;
injectionRateSuSp
{
e.air (500 0); // kg*m^2/s^3
}
}
}


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

Hello, I have a similar problem. Did you find a solution?

[NBad]

Looking at the code of compressibleMultiphaseInterFoam it looks like fvOptions is not implemented into the equations. For example the equation for U is:

    
    fvVectorMatrix UEqn
    (
        fvm::ddt(rho, U)
      + fvm::div(mixture.rhoPhi(), U)
      + turbulence->divDevRhoReff(U)
    );

What you can do is to look at the code of some other solver and implement fvOptions. For example multiphaseInterFoam has fvOptions included:

    fvVectorMatrix UEqn
    (
        fvm::ddt(rho, U) + fvm::div(rhoPhi, U)
      - fvm::Sp(contErr, U)
      + MRF.DDt(rho, U)
      + turbulence.divDevRhoReff(U)
     ==
        fvOptions(rho, U)
    );