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

porousInterFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 30, 2013, 12:14
Default porousInterFoam
  #1
New Member
 
hesam
Join Date: Mar 2013
Posts: 13
Rep Power: 4
hesamgh is on a distinguished road
hi,

i am trying to compile the porousInterFoam solver, but whem i run the example, this error occurs ,

please help me to solve this problem

MULES: Solving for alpha1
Phase-1 volume fraction = 0.507435 Min(alpha1) = -0.321967 Max(alpha1) = 1.62518
MULES: Solving for alpha1
Phase-1 volume fraction = 0.507435 Min(alpha1) = -0.321915 Max(alpha1) = 1.62518
DICPCG: Solving for p_rgh, Initial residual = 0.924755, Final residual = 0.0281868, No Iterations 257
time step continuity errors : sum local = 2.13828e-05, global = 1.74228e-07, cumulative = -1.00629
DICPCG: Solving for p_rgh, Initial residual = 0.645343, Final residual = 0.030643, No Iterations 259
time step continuity errors : sum local = 0.000109586, global = 3.36113e-07, cumulative = -1.00629
DICPCG: Solving for p_rgh, Initial residual = 0.748674, Final residual = 9.91322e-08, No Iterations 446
time step continuity errors : sum local = 2.27742e-09, global = 3.95198e-11, cumulative = -1.00629
ExecutionTime = 3.25 s ClockTime = 3 s

Courant Number mean: 0.146919 max: 224.824
Interface Courant Number mean: 0.00106774 max: 6.39375
deltaT = 8.82449e-15
--> FOAM Warning :
From function Time:perator++()
in file db/Time/Time.C at line 1029
Increased the timePrecision from 10 to 11 to distinguish between timeNames at time 0.00459495
Time = 0.0045949473993

--> FOAM Warning :
From function Time:perator++()
in file db/Time/Time.C at line 1029
Increased the timePrecision from 11 to 12 to distinguish between timeNames at time 0.00459495
MULES: Solving for alpha1
Phase-1 volume fraction = 0.507435 Min(alpha1) = -0.321899 Max(alpha1) = 1.62518
--> FOAM Warning :
From function Time:perator++()
in file db/Time/Time.C at line 1029
Increased the timePrecision from 12 to 13 to distinguish between timeNames at time 0.00459495
MULES: Solving for alpha1
Phase-1 volume fraction = 0.507435 Min(alpha1) = -0.321884 Max(alpha1) = 1.62518
DICPCG: Solving for p_rgh, Initial residual = 0.997524, Final residual = 0.0482074, No Iterations 119
time step continuity errors : sum local = 1.98589e-05, global = -2.09955e-06, cumulative = -1.00629
DICPCG: Solving for p_rgh, Initial residual = 0.109669, Final residual = 0.00416443, No Iterations 291
time step continuity errors : sum local = 3.66002e-06, global = -1.02672e-08, cumulative = -1.00629
DICPCG: Solving for p_rgh, Initial residual = 0.0447736, Final residual = 9.61283e-08, No Iterations 531
time step continuity errors : sum local = 8.9133e-11, global = 3.26094e-14, cumulative = -1.00629
ExecutionTime = 3.91 s ClockTime = 4 s

Courant Number mean: 0.000801206 max: 1.18367
Interface Courant Number mean: 1.45436e-05 max: 0.066757
deltaT = 3.72759e-15
Time = 0.004594947399295

MULES: Solving for alpha1
Phase-1 volume fraction = 0.507435 Min(alpha1) = -0.321872 Max(alpha1) = 1.62518
MULES: Solving for alpha1
Phase-1 volume fraction = 0.507435 Min(alpha1) = -0.32186 Max(alpha1) = 1.62518
DICPCG: Solving for p_rgh, Initial residual = 0.751844, Final residual = 0.0271748, No Iterations 206
time step continuity errors : sum local = 1.02263e-05, global = -4.9812e-08, cumulative = -1.00629
DICPCG: Solving for p_rgh, Initial residual = 0.138542, Final residual = 0.00650515, No Iterations 255
time step continuity errors : sum local = 4.51487e-06, global = -1.06425e-08, cumulative = -1.00629
DICPCG: Solving for p_rgh, Initial residual = 0.0540249, Final residual = 9.98206e-08, No Iterations 410
time step continuity errors : sum local = 7.93589e-11, global = 3.34275e-14, cumulative = -1.00629
ExecutionTime = 4.54 s ClockTime = 5 s

Courant Number mean: 0.00291016 max: 2.35784
Interface Courant Number mean: 7.37114e-05 max: 0.327241
deltaT = 7.90465e-16
Time = 0.004594947399296

MULES: Solving for alpha1
Phase-1 volume fraction = 0.507435 Min(alpha1) = -0.321613 Max(alpha1) = 1.62518
--> FOAM Warning :
From function Time:perator++()
in file db/Time/Time.C at line 1029
Increased the timePrecision from 13 to 14 to distinguish between timeNames at time 0.00459495
MULES: Solving for alpha1
Phase-1 volume fraction = 0.507435 Min(alpha1) = -0.321366 Max(alpha1) = 1.62518
DICPCG: Solving for p_rgh, Initial residual = 0.510742, Final residual = 0.010461, No Iterations 272
time step continuity errors : sum local = 1.17433e-05, global = -5.60236e-09, cumulative = -1.00629
DICPCG: Solving for p_rgh, Initial residual = 0.379362, Final residual = 0.0155277, No Iterations 196
time step continuity errors : sum local = 3.53471e-05, global = -1.73532e-06, cumulative = -1.00629
DICPCG: Solving for p_rgh, Initial residual = 0.270411, Final residual = 5.46214e-08, No Iterations 468
time step continuity errors : sum local = 1.62082e-10, global = 2.19395e-13, cumulative = -1.00629
ExecutionTime = 5.2 s ClockTime = 5 s

Courant Number mean: 0.00519716 max: 4.91789
Interface Courant Number mean: 0.000167649 max: 0.556441
deltaT = 8.03662e-17
Time = 0.0045949473992963

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4 void Foam::MULES::limit<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::zeroField const&, Foam::zeroField const&, double, double, int, bool) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5 Foam::MULES::explicitSolve(Foam::GeometricField<do uble, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#6
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/interFoam"
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/interFoam"
Floating point exception (core dumped)

hesamgh is offline   Reply With Quote

Old   October 30, 2013, 14:22
Default
  #2
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,619
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Hallo Hesam,

Your time steps are in the order of machine precision, so something went wrong long time before. I suspect that your boundary conditions or the definition of your porous zone(s) are not correct.

Kind regards

Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is offline   Reply With Quote

Old   November 15, 2013, 04:05
Default porous
  #3
New Member
 
hesam
Join Date: Mar 2013
Posts: 13
Rep Power: 4
hesamgh is on a distinguished road
hi dear friend
excuse me for delay

do you know how can i identify the porous medium in openfoam when import the shape from gambit to openfoam?
i make the porous media in gambit from this address:
operation,zones,specify continuum types

thanks alot for your help
hesamgh is offline   Reply With Quote

Old   November 15, 2013, 07:04
Default
  #4
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,123
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
you can use topoSet to create your zone look following tutorials to find how to work with topoSet:
Quote:
./run/tutorials/lagrangian/reactingParcelFilmFoam/splashPanel/system/wallFilmRegion.topoSet
./run/tutorials/lagrangian/reactingParcelFilmFoam/splashPanel/system/wallFilmRegion/createWallFilmRegionPatches.topoSet
./run/tutorials/lagrangian/reactingParcelFilmFoam/hotBoxes/system/topoSetDict
./run/tutorials/lagrangian/reactingParcelFilmFoam/cylinder/system/topoSetDict
./run/tutorials/lagrangian/reactingParcelFilmFoam/rivuletPanel/system/wallFilmRegion.topoSet
./run/tutorials/lagrangian/coalChemistryFoam/simplifiedSiwek/system/topoSetDict
./run/tutorials/lagrangian/reactingParcelFoam/filter/system/topoSetDict
./run/tutorials/mesh/moveDynamicMesh/simpleHarmonicMotion/system/topoSetDict
./run/tutorials/multiphase/MRFMultiphaseInterFoam/mixerVessel2D/system/topoSetDict
./run/tutorials/multiphase/twoPhaseEulerFoam/mixerVessel2D/system/topoSetDict
./run/tutorials/multiphase/interDyMFoam/ras/damBreakWithObstacle/log.topoSet
./run/tutorials/multiphase/interDyMFoam/ras/damBreakWithObstacle/system/topoSetDict
./run/tutorials/multiphase/interDyMFoam/ras/floatingObject/system/topoSetDict
./run/tutorials/multiphase/MRFInterFoam/mixerVessel2D/system/topoSetDict
./run/tutorials/multiphase/interFoam/ras/damBreakPorousBaffle/system/topoSetDict
./run/tutorials/multiphase/interFoam/les/nozzleFlow2D/system/topoSetDict.2
./run/tutorials/multiphase/interFoam/les/nozzleFlow2D/system/topoSetDict.1
./run/tutorials/multiphase/compressibleTwoPhaseEulerFoam/mixerVessel2D/system/topoSetDict
./run/tutorials/multiphase/cavitatingFoam/ras/throttle/system/topoSetDict.3
./run/tutorials/multiphase/cavitatingFoam/ras/throttle/system/topoSetDict.2
./run/tutorials/multiphase/cavitatingFoam/ras/throttle/system/topoSetDict.1
./run/tutorials/multiphase/cavitatingFoam/les/throttle3D/system/topoSetDict.3
./run/tutorials/multiphase/cavitatingFoam/les/throttle3D/system/topoSetDict.2
./run/tutorials/multiphase/cavitatingFoam/les/throttle3D/system/topoSetDict.1
./run/tutorials/multiphase/cavitatingFoam/les/throttle/system/topoSetDict.3
./run/tutorials/multiphase/cavitatingFoam/les/throttle/system/topoSetDict.2
./run/tutorials/multiphase/cavitatingFoam/les/throttle/system/topoSetDict.1
./run/tutorials/multiphase/multiphaseEulerFoam/mixerVessel2D/system/topoSetDict
./run/tutorials/combustion/PDRFoam/flamePropagationWithObstacles/system/topoSetDict
./run/tutorials/combustion/fireFoam/les/smallPoolFire3D/system/topoSetDict
./run/tutorials/combustion/fireFoam/les/oppositeBurningPanels/f.topoSetDict
./run/tutorials/combustion/fireFoam/les/oppositeBurningPanels/fBurner.topoSetDict
./run/tutorials/combustion/fireFoam/les/oppositeBurningPanels/cRefine.topoSetDict
./run/tutorials/combustion/fireFoam/les/smallPoolFire2D/system/topoSetDict
./run/tutorials/incompressible/pimpleDyMFoam/mixerVesselAMI2D/system/topoSetDict
./run/tutorials/incompressible/pimpleDyMFoam/propeller/system/createAMIFaces.topoSetDict~
./run/tutorials/incompressible/pimpleDyMFoam/propeller/system/createInletOutletSets.topoSetDict
./run/tutorials/incompressible/pimpleDyMFoam/propeller/system/removeRedundantZones.topoSetDict
./run/tutorials/incompressible/pimpleDyMFoam/propeller/system/createAMIFaces.topoSetDict
./run/tutorials/incompressible/simpleFoam/turbineSiting/system/topoSetDict
./run/tutorials/incompressible/simpleFoam/pipeCyclic/system/topoSetDict
./run/tutorials/incompressible/pimpleFoam/TJunctionFan/system/topoSetDict
./run/tutorials/incompressible/pimpleFoam/elipsekkLOmega/system/topoSetDict
./run/tutorials/incompressible/potentialFreeSurfaceFoam/oscillatingBox/system/topoSetDict
./run/tutorials/heatTransfer/chtMultiRegionSimpleFoam/heatExchanger/log.topoSet.air.2
./run/tutorials/heatTransfer/chtMultiRegionSimpleFoam/heatExchanger/log.topoSet.air.1
./run/tutorials/heatTransfer/chtMultiRegionSimpleFoam/heatExchanger/system/air/topoSetDict.2
./run/tutorials/heatTransfer/chtMultiRegionSimpleFoam/heatExchanger/system/air/topoSetDict.1
./run/tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeater/system/topoSetDict
./run/tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeaterRadiation/system/topoSetDict
./run/tutorials/heatTransfer/buoyantSimpleFoam/circuitBoardCooling/system/topoSetDict
./run/tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/topoSetDict
./run/tutorials/heatTransfer/chtMultiRegionFoam/multiRegionLiquidHeater/system/topoSetDict
__________________
Training Course on OpenFOAM at (http://www.isme.ir/)
My Weblog (http://openfoam.blogfa.com/)
nimasam 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
Porosity and permeability in PorousInterFoam jasouza1974 OpenFOAM 14 August 12, 2015 04:16
help adding the energy equation to porousinterfoam using the enthalpy formulation nadine OpenFOAM Programming & Development 18 June 17, 2014 08:39
porousInterFoam sachinlb OpenFOAM Running, Solving & CFD 0 December 27, 2012 07:30
porousinterFoam mcdonalds OpenFOAM Running, Solving & CFD 2 July 17, 2012 15:44
porousInterFoam JF.JTT OpenFOAM 2 May 1, 2011 14:05


All times are GMT -4. The time now is 09:56.