[atlan]

Hi,

is it possible to model porous zone with buoyantSimpleFoam or bouyantPimpleFoam? I have tried to add the porous zone into fvOptions but the solver failed (I think) due to messing thermodynamic properties of the porous zone. Can somebody advise me please how to implement this?


Thanks

[adambarfi]

hi atlan,

I tried it before by modifying the buoyantBoussinesqSimpleFoam but I gave up. and still I wanna do it, if my works let me do that!
what's the error? did u make the solver correctly?

Regards,
Mostafa

[atlan]

Hi Mostafa,
Thank you for the answer.

I have not made any solver, which is probably the problem. Is there any solver which is applicable for the calculation of heat transfer (wall conduction) including the porous wall?
Regards

[adambarfi]

AFAIK, nope, there isn't such a solver.
you should add energy equation for clear fluid and porous zone yourself.

[jherb]

I am not sure, that I understand what you want to do, but I used the following fvOptions successfully with buoyantPimpleFoam:

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

porosityBlockage
{
    type            explicitPorositySource;
    active          off;
    selectionMode   cellZone;
    cellZone        porous;

    explicitPorositySourceCoeffs
    {
        type            DarcyForchheimer;

        DarcyForchheimerCoeffs
        {
            d   d [0 -2 0 0 0] (1e9 1e9 0);
            f   f [0 -1 0 0 0] (1e9 1e9 0);

            coordinateSystem
            {
                e1  (1 0 0);
                e2  (0 1 0);
            }
        }
    }
}


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

(this example was used to force a certain velocity direction in the porous zone)

[atlan]

Thank you,

This really works. In fvOptions you can also add the porous zone temperature.

Atlan

[dewey]

I'm trying to run buoyantSimpleFoam with a porous media fvoptions like you said and then stop, giving the following message:




Could you help me, or give me an advice please

[dewey]

Hello., im trying to simulate a pipe with porous media and heat transfer using buoyant simplefoam with fvoptions porousmedia.

Do you know why the simulation stop?
I dont understand what mean that error.

thank u for your time


Time = 26

DILUPBiCG: Solving for Ux, Initial residual = 0.0623969, Final residual = 0.000334212, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.101444, Final residual = 0.000653273, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.149451, Final residual = 0.000506783, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 0.00689413, Final residual = 1.16579e-05, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 0.14369, Final residual = 0.00132583, No Iterations 4
time step continuity errors : sum local = 0.0131216, global = 0.00379873, cumulative = 0.00339527
rho max/min : 2.66735 0.225983
ExecutionTime = 336.84 s ClockTime = 389 s

Time = 27

DILUPBiCG: Solving for Ux, Initial residual = 0.0921637, Final residual = 0.000623586, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.12727, Final residual = 0.00119562, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.277155, Final residual = 0.00148346, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 0.00678693, Final residual = 6.33452e-06, No Iterations 2
[0] #0 Foam::error:rintStack(Foam::Ostream&) at ??:?
[0] #1 Foam::sigFpe::sigHandler(int) at ??:?
[0] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
[0] #4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
[0] #5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[0] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
[0] #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
[0] #8 Foam::fvMatrix<double>::solve() at ??:?
[0] #9 ? at ??:?
[0] #10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #11 ? at ??:?
[Guacamaya:09486] *** Process received signal ***
[Guacamaya:09486] Signal: Floating point exception (8)
[Guacamaya:09486] Signal code: (-6)
[Guacamaya:09486] Failing at address: 0x3e80000250e
[Guacamaya:09486] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f0ff0dce4b0]
[Guacamaya:09486] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f0ff0dce428]
[Guacamaya:09486] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f0ff0dce4b0]
[Guacamaya:09486] [ 3] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5scaleERNS_5Fi eldIdEES3_RKNS_9lduMatrixERKNS_10FieldFieldIS1_dEE RKNS_8UPtrListIKNS_17lduInterfaceFieldEEERKS2_h+0x ce)[0x7f0ff2039b4e]
[Guacamaya:09486] [ 4] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7 PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS 8_S9_S9_S9_S9_S9_RNS1_IS8_EESD_h+0x7b2)[0x7f0ff203ddb2]
[Guacamaya:09486] [ 5] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5Fi eldIdEERKS2_h+0x86d)[0x7f0ff204070d]
[Guacamaya:09486] [ 6] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegr egatedERKNS_10dictionaryE+0x15b)[0x7f0ff5a5f95b]
[Guacamaya:09486] [ 7] buoyantSimpleFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_ 10dictionaryE+0x191)[0x493971]
[Guacamaya:09486] [ 8] buoyantSimpleFoam(_ZN4Foam8fvMatrixIdE5solveEv+0x1 35)[0x493c25]
[Guacamaya:09486] [ 9] buoyantSimpleFoam[0x4306bc]
[Guacamaya:09486] [10] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f0ff0db9830]
[Guacamaya:09486] [11] buoyantSimpleFoam[0x433539]
[Guacamaya:09486] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 9486 on node Guacamaya exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------