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

rhoPorousMRFPimpleFoam velocity

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 9, 2011, 05:50
Default rhoPorousMRFPimpleFoam velocity
  #1
New Member
 
Raph Raggatt
Join Date: Jul 2011
Posts: 9
Rep Power: 14
Rapha is on a distinguished road
Hi there,

I am using the rhoPorousMRFPimpleFoam solver for my base case however I am having problems with the output, getting increases in velocity and pressure when I should be getting the opposite. I have included the temperature equation into the solver.
The details are:
- Mesh imported from Salome with 5 meshes joined together in a compound mesh.
- Each boundary is specified, eg inlet, wall, regionCouple, etc
- Two porous regions either side of the inlet channel
- Inlet velocity of 30 m/s. Inlet temp of 700K
- the model should run by air flowing through the inlet at 700K, 30m/s, passing through and heating up the porous zones on either side of the inlet duct leaving the outlets on either side at a lower temperature.

The problems are several:
- The time scale I can reach is of the order 1e-4 s when really it should be of the order 1e3 s.
- The velocity increases rapidly to 267m/s at certain regions - why?
- The pressure increases from 1 bar to 1.5 bar
- The flow doesn't go from the inlet to the outlet through the porous zones.

If anyone could give any hints as to whether it's the fvschemes, solutions or solver I need to look at, I would be really grateful.

Regards,

Rapha
Rapha is offline   Reply With Quote

Old   August 9, 2011, 05:52
Default Error Code
  #2
New Member
 
Raph Raggatt
Join Date: Jul 2011
Posts: 9
Rep Power: 14
Rapha is on a distinguished road
This is the error code which I get for the last iteration:

Courant Number mean: 0.0166694 max: 0.605291
deltaT = 1.8885e-07
Time = 0.000333925

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.0115056, Final residual = 4.09179e-06, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.0112349, Final residual = 3.55671e-06, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 0.278506, Final residual = 1.42709e-07, No Iterations 2
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::hPsiThermo<Foam:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > >::calculate() in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libbasicThermophysicalModels.so"
#4 Foam::hPsiThermo<Foam:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > >::correct() in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libbasicThermophysicalModels.so"
#5
in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/HTTESS_Solver_P"
#6 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#7
in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/HTTESS_Solver_P"
Floating point exception
ra@ubuntu:~/OpenFOAM/ra-2.0.0/run/3_08/HTTESS$
Rapha is offline   Reply With Quote

Reply

Tags
rhoporousmrfpimplefoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Plotting Radial Velocity and Tangential Velocity in CFD Post ashtonJ CFX 5 July 13, 2015 03:49
Velocity in Porous medium : HELP! HELP! HELP! Kali Sanjay Phoenics 0 November 6, 2006 07:10
Neumann pressure BC and velocity field Antech Main CFD Forum 0 April 25, 2006 03:15
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 10:11
what the result is negatif pressure at inlet chong chee nan FLUENT 0 December 29, 2001 06:13


All times are GMT -4. The time now is 01:49.