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

Buoyancydriven cavity flow

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 12, 2008, 11:40
Default Hi, I try to simulate a buo
  #1
New Member
 
Anita K.
Join Date: Mar 2009
Posts: 25
Rep Power: 8
anita is on a distinguished road
Hi,

I try to simulate a buoyancy-driven cavity flow like the example in "Computational Methods for Fluid Dynamics".
It is a square with adiabatic walls on the upper and lower side, a hot wall on the right side and a cold wall on left side. I define the temperature hot = 293 K, cold = 253 K and default 273 K.

The problem is that after starting the buoyantSimpleFoam solver there is an error.

Does somebody has a suggestion where the problem could be? I used the hotRoom-tutorial as basis.
I only changed turbulence model to laminar. Also with the k-epsilon modell there is the same problem.

Anita

Here is the log:

Exec : buoyantSimpleFoam . buoyancydriven32
Date : Mar 12 2008
Time : 16:26:08
Host : st31
PID : 25414
Root : /home/kienzani/OpenFOAM/kienzani-1.4.1/run/ownTests/ExamplesBook
Case : buoyancydriven32
Nprocs : 1
Create time

Create mesh for time = 0


Reading environmentalProperties
Reading thermophysical properties

Selecting thermodynamics package hThermo<puremixture<consttransport<speciethermo<hc onstthermo<perfectgas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model laminar
Calculating field g.h

Creating field pd


Starting time loop

Time = 1

DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 9.69948e-07, No Iterations 5
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 1.36774e-06, No Iterations 5
#0 Foam::error::printStack(Foam:stream&) in "/home/kienzani/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/kienzani/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xffffe420]
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/kienzani/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 Foam::tmp<foam::geometricfield<double,> > Foam::operator/<foam::fvspatchfield,>(Foam::GeometricField<double ,> const&, Foam::tmp<foam::geometricfield<double,> > const&) in "/home/kienzani/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/buoyantSi mpleFoam"
#5 main in "/home/kienzani/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/buoyantSi mpleFoam"
#6 __libc_start_main in "/lib/libc.so.6"
#7 Foam::regIOobject::readIfModified() in "/home/kienzani/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/buoyantSi mpleFoam"
anita is offline   Reply With Quote

Old   March 16, 2008, 09:54
Default Hi Anita, I have the same pro
  #2
New Member
 
aksouh mohamed
Join Date: Mar 2009
Location: Algiers
Posts: 2
Rep Power: 0
aksouh_med is on a distinguished road
Hi Anita,
I have the same problem, but for my problem is just 3d, and I have some results but I think that it didn't converge...so I had the errors message like :

Time = 139

Lookup gradScheme for grad(U)
Lookup divScheme for div((muEff*dev2(grad(U).T())))
Lookup laplacianScheme for laplacian(muEff,U)
Lookup divScheme for div(phi,U)
Lookup gradScheme for grad(rho)
Lookup gradScheme for grad(pd)
DILUPBiCG: Solving for Ux, Initial residual = 0.249057, Final residual = 1.70173e-06, No Iterations 5
DILUPBiCG: Solving for Uy, Initial residual = 0.0185925, Final residual = 4.28971e-06, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.104714, Final residual = 1.6076e-06, No Iterations 5
Lookup interpolationScheme for interpolate(rho)
Lookup interpolationScheme for interpolate(p)
Lookup interpolationScheme for interpolate(rho)
Lookup laplacianScheme for laplacian(alphaEff,h)
Lookup divScheme for div(phi,h)
DILUPBiCG: Solving for h, Initial residual = 0.00378344, Final residual = 2.23355e-06, No Iterations 3
Lookup interpolationScheme for interpolate((rho*U))
Lookup snGradScheme for snGrad(rho)
Lookup interpolationScheme for interpolate(((rho*gh)*(1|A(U))))
Lookup laplacianScheme for laplacian((rho*(1|A(U))),pd)
DICPCG: Solving for pd, Initial residual = 0.332361, Final residual = 9.81432e-09, No Iterations 429
Lookup fluxRequired for pd
time step continuity errors : sum local = 8.58238e-08, global = -1.04862e-18, cumulative = 3.09173e-18
Lookup gradScheme for grad(rho)
Lookup gradScheme for grad(pd)
rho max/min : 1.20886 1.20063
ExecutionTime = 11054.2 s ClockTime = 11061 s


Anita try to change yr laminar model to turbulence model, but for low Reynolds number, like "LaunderSharmaKE'

so for me my questions are :
1- I don't know if I should change the interpolation scheme, and if there is the case which the good scheme for my problem?
2- I know that my problem is the heat transfer, so for the density, I have the boussinesq density in order to create flow between the cold and hot wall, but I dont know where can I find the boussinesq equation fo density.

thank you very much for your help.
aksouh_med is offline   Reply With Quote

Old   March 16, 2008, 10:14
Default Sorry, but I forgot to sent y
  #3
New Member
 
aksouh mohamed
Join Date: Mar 2009
Location: Algiers
Posts: 2
Rep Power: 0
aksouh_med is on a distinguished road
Sorry,
but I forgot to sent you my fvScheme file :



ddtSchemes
{
default steadyState;
}

gradSchemes
{
default Gauss linear;
}

divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phi,h) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
div((muEff*dev2(grad(U).T()))) Gauss linear;
}

laplacianSchemes
{
default none;
laplacian(muEff,U) Gauss linear corrected;
laplacian((rho*(1|A(U))),pd) Gauss linear corrected;
laplacian(alphaEff,h) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
laplacian(DREff,R) Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
pd;
}



Thanks for all
Best regards
aksouh_med is offline   Reply With Quote

Old   March 26, 2008, 11:24
Default Hi, now the solver is worki
  #4
New Member
 
Anita K.
Join Date: Mar 2009
Posts: 25
Rep Power: 8
anita is on a distinguished road
Hi,

now the solver is working correct. I had initialized pressure with zero. After changing to 10000 calculation worked fine.

The new problem is that residual of U is falling the first 2500 slowly until below 0.001. Than it rises until 0.1. Afterwards it is nearly constant. Pressure residual does not even fall below 0.01.
anita 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
Buoyancydriven cavity results for comparision kar OpenFOAM 0 April 14, 2008 09:04
2D Cavity Flow ALI RAJABPOUR Main CFD Forum 0 July 27, 2005 01:55
Flow in a cavity Herry CFX 1 April 27, 2003 18:30
Cavity flow Shannag Main CFD Forum 5 September 7, 2001 14:03
Cavity flow Sergey Smirnov Main CFD Forum 8 October 26, 2000 11:50


All times are GMT -4. The time now is 13:42.