CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   buoyantSimpleFoam + porous zone (https://www.cfd-online.com/Forums/openfoam-solving/138918-buoyantsimplefoam-porous-zone.html)

atlan July 14, 2014 07:42
buoyantSimpleFoam + porous zone
 
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 July 14, 2014 11:09
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 July 16, 2014 09:48
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 July 16, 2014 09:53
AFAIK, nope, there isn't such a solver.
you should add energy equation for clear fluid and porous zone yourself.

jherb July 20, 2014 09:03
I am not sure, that I understand what you want to do, but I used the following fvOptions successfully with buoyantPimpleFoam:
Code:
/*--------------------------------*- 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 July 24, 2014 03:33
Porous zone + BuoyantSimpleFoam - SOLVED
 
Thank you,

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

Atlan

dewey April 25, 2017 13:05
1 Attachment(s)
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 May 3, 2017 13:37
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::printStack(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).
--------------------------------------------------------------------------


All times are GMT -4. The time now is 04:17.