CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

buoyantSimpleFoam + porous zone

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   July 14, 2014, 07:42
Default buoyantSimpleFoam + porous zone
  #1
New Member
 
Join Date: Jul 2014
Posts: 21
Rep Power: 5
atlan is on a distinguished road
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
atlan is offline   Reply With Quote

Old   July 14, 2014, 11:09
Default
  #2
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 319
Rep Power: 8
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to 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
adambarfi is offline   Reply With Quote

Old   July 16, 2014, 09:48
Default
  #3
New Member
 
Join Date: Jul 2014
Posts: 21
Rep Power: 5
atlan is on a distinguished road
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
atlan is offline   Reply With Quote

Old   July 16, 2014, 09:53
Default
  #4
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 319
Rep Power: 8
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
AFAIK, nope, there isn't such a solver.
you should add energy equation for clear fluid and porous zone yourself.
adambarfi is offline   Reply With Quote

Old   July 20, 2014, 09:03
Default
  #5
Senior Member
 
Joachim Herb
Join Date: Sep 2010
Posts: 425
Rep Power: 13
jherb is on a distinguished road
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)
jherb is offline   Reply With Quote

Old   July 24, 2014, 03:33
Default Porous zone + BuoyantSimpleFoam - SOLVED
  #6
New Member
 
Join Date: Jul 2014
Posts: 21
Rep Power: 5
atlan is on a distinguished road
Thank you,

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

Atlan
atlan is offline   Reply With Quote

Old   April 25, 2017, 13:05
Unhappy
  #7
Member
 
alberto
Join Date: Apr 2016
Location: Mexico
Posts: 76
Rep Power: 3
dewey is on a distinguished road
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
Attached Images
File Type: jpg Captura de pantalla de 2017-04-25 11-44-43.jpg (157.3 KB, 13 views)
dewey is offline   Reply With Quote

Old   May 3, 2017, 13:37
Default
  #8
Member
 
alberto
Join Date: Apr 2016
Location: Mexico
Posts: 76
Rep Power: 3
dewey is on a distinguished road
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).
--------------------------------------------------------------------------
dewey is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Modelling Combustion in Porous Zone tanjinjack FLUENT 2 September 26, 2016 04:10
Possible Bug in pimpleFoam (or createPatch) (or fluent3DMeshToFoam) cfdonline2mohsen OpenFOAM 3 October 21, 2013 09:28
fluentMeshToFoam multidomain mesh conversion problem Attesz OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 12 May 2, 2013 10:52
Problem in running ICEM grid in Openfoam Tarak OpenFOAM 6 September 9, 2011 17:51
Problem in IMPORT of ICEM input file in FLUENT csvirume FLUENT 2 September 9, 2009 01:08


All times are GMT -4. The time now is 08:54.