Problem running buoyantFoam
Hello,
I try to get familiar was solving temperature Problems in Openfoam. at the Moment I'm still working with 1.4 until our admin will Install 1.6.. I was running the hotRoom tutorial which is usings Blockmesh. Because I will not use Ansa for further work not the Blockmesh in reconstructed the mesh using ansa. Ansa 13.xx is able to read and write OpenFoam decks (ok the writing has a bug).. I opened the OpenFoam deck with Ansa , saved it as fluentfile and used fluent3DMeshtofoam. There is now a cellZones file in the polymesh direktory, the boundary an initial conditions a took from the tutorial. After 3 steps the calculation crashes 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 kEpsilon Creating field dpdt Calculating field g.h Creating field pd Starting time loop Courant Number mean: 0 max: 0 Time = 0.25 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for Ux, Initial residual = 0.00508188, Final residual = 7.32216e06, No Iterations 20 DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 7.08488e06, No Iterations 36 DILUPBiCG: Solving for Uz, Initial residual = 0.00144893, Final residual = 9.13031e06, No Iterations 37 DILUPBiCG: Solving for h, Initial residual = 0.697128, Final residual = 7.01113e06, No Iterations 43 DICPCG: Solving for pd, Initial residual = 0.999995, Final residual = 6.62537e07, No Iterations 45 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 1.21113e05, global = 2.976e11, cumulative = 2.976e11 DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 3.92272e06, No Iterations 41 DICPCG: Solving for pd, Initial residual = 0.673709, Final residual = 7.38414e07, No Iterations 44 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 2.00443e05, global = 8.60306e11, cumulative = 5.62706e11 DILUPBiCG: Solving for epsilon, Initial residual = 0.126125, Final residual = 5.38938e06, No Iterations 11 DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 3.82471e06, No Iterations 33 ExecutionTime = 0.25 s ClockTime = 0 s Courant Number mean: 1.79445 max: 22.9819 Time = 0.5 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for Ux, Initial residual = 0.998631, Final residual = 7.49322e06, No Iterations 14 DILUPBiCG: Solving for Uy, Initial residual = 0.999779, Final residual = 3.38204e06, No Iterations 15 DILUPBiCG: Solving for Uz, Initial residual = 0.998631, Final residual = 7.49322e06, No Iterations 14 DILUPBiCG: Solving for h, Initial residual = 0.931951, Final residual = 7.36059e06, No Iterations 19 DICPCG: Solving for pd, Initial residual = 0.872791, Final residual = 9.98979e07, No Iterations 28 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.0984091, global = 0.000134958, cumulative = 0.000134958 DILUPBiCG: Solving for h, Initial residual = 0.997817, Final residual = 6.95462e06, No Iterations 48 DICPCG: Solving for pd, Initial residual = 0.699228, Final residual = 9.09491e07, No Iterations 27 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.0906934, global = 1.07524e06, cumulative = 0.000136033 DILUPBiCG: Solving for epsilon, Initial residual = 0.99387, Final residual = 6.26957e06, No Iterations 16 DILUPBiCG: Solving for k, Initial residual = 0.999733, Final residual = 3.80444e06, No Iterations 17 bounding k, min: 0.683897 max: 8.39597 average: 2.23756 ExecutionTime = 0.38 s ClockTime = 0 s Courant Number mean: 29703.4 max: 101397 Time = 0.75 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for Ux, Initial residual = 0.935259, Final residual = 7.06474e06, No Iterations 58 DILUPBiCG: Solving for Uy, Initial residual = 0.932357, Final residual = 8.07299e06, No Iterations 44 DILUPBiCG: Solving for Uz, Initial residual = 0.935259, Final residual = 7.06474e06, No Iterations 58 DILUPBiCG: Solving for h, Initial residual = 0.927174, Final residual = 3.32234e06, No Iterations 45 DICPCG: Solving for pd, Initial residual = 0.999669, Final residual = 5.44238e07, No Iterations 45 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.0249903, global = 1.14473e06, cumulative = 0.000137178 DILUPBiCG: Solving for h, Initial residual = 0.999947, Final residual = 3.96725e06, No Iterations 45 DICPCG: Solving for pd, Initial residual = 0.489969, Final residual = 8.87186e07, No Iterations 71 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.000912857, global = 1.0761e06, cumulative = 0.000138254 DILUPBiCG: Solving for epsilon, Initial residual = 0.999279, Final residual = 2.48876e06, No Iterations 6 bounding epsilon, min: 1.23011e+07 max: 8.26133e+09 average: 4.07851e+07 DILUPBiCG: Solving for k, Initial residual = 0.837395, Final residual = 6.2143e06, No Iterations 9 bounding k, min: 229.64 max: 426.695 average: 19.3091 ExecutionTime = 0.6 s ClockTime = 0 s Courant Number mean: 170.246 max: 13473.3 Time = 1 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for Ux, Initial residual = 0.736069, Final residual = 9.20155e06, No Iterations 61 DILUPBiCG: Solving for Uy, Initial residual = 0.726148, Final residual = 2.06373e06, No Iterations 54 DILUPBiCG: Solving for Uz, Initial residual = 0.736069, Final residual = 9.56633e06, No Iterations 61 DILUPBiCG: Solving for h, Initial residual = 0.893243, Final residual = 2.53873e06, No Iterations 52 > FOAM FATAL ERROR : Maximum number of iterations exceeded#0 Foam::error::printStack(Foam::Ostream&) in "/soft/OpenFOAM/linux64/OpenFOAM/OpenFOAM1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::error::abort() in "/soft/OpenFOAM/linux64/OpenFOAM/OpenFOAM1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so" #2 Foam::hThermo<Foam::pureMixture<Foam::constTranspo rt<Foam::specieThermo<Foam::hConstThermo<Foam::per fectGas> > > > >::calculate() in "/soft/OpenFOAM/linux64/OpenFOAM/OpenFOAM1.4.1/lib/linux64GccDPOpt/libbasicThermophysicalModels.so" #3 Foam::hThermo<Foam::pureMixture<Foam::constTranspo rt<Foam::specieThermo<Foam::hConstThermo<Foam::per fectGas> > > > >::correct() in "/soft/OpenFOAM/linux64/OpenFOAM/OpenFOAM1.4.1/lib/linux64GccDPOpt/libbasicThermophysicalModels.so" #4 main in "/soft/OpenFOAM/linux64/OpenFOAM/OpenFOAM1.4.1/applications/bin/linux64GccDPOpt/buoyantFoam" #5 __libc_start_main in "/lib64/libc.so.6" #6 Foam::regIOobject::readIfModified() in "/soft/OpenFOAM/linux64/OpenFOAM/OpenFOAM1.4.1/applications/bin/linux64GccDPOpt/buoyantFoam" From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const in file /home/dm2/henry/OpenFOAM/OpenFOAM1.4.1/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 83. Does somebody know the reason? Greetings Anke 
Problem solved,

Hi,
i'm not sure to have understood your question about ANSA and blockMesh but i think you have instability problems: the courant number increases untill the crash of the simulation. You should try to use smaller timesteps, maybe a finer mesh and check that your boundary conditions are correct. dd 
Quote:
Which ANSA version (v13.0.x, what is the x?) are you using and what kind of bug have you found in the output? 
1 Attachment(s)
Quote:
Can you kindly advise me how you resolved this issue. I have the same problem in rhoPimpleFoam. I have attached my error file. Can you kindly advise how I can solve it. Case is air flowing in pipe and assumed as compressible. I shall be grateful for your support. Regards 
All times are GMT 4. The time now is 16:40. 