[samiam1000]

Thank you so much Tobias.

This really helped.

[samiam1000]

Hi Tobias,

can I ask you to share the fvSolution and fvSchemes files that you used for the solid regions. I think there's something wrong in mine, since I notice a strange behvior in domain2 region (look at the strange values of temperature) and the simulation is taking a lot of time to solve domain3.

Time = 1


Solving for fluid region fluid
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 0.00545485, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 0.0211023, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 0.00876899, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 1, Final residual = 0.0972692, No Iterations 11
Min/max T:270.103 305.567
DICPCG:  Solving for p_rgh, Initial residual = 0.61935, Final residual = 0.00590249, No Iterations 316
time step continuity errors : sum local = 0.0921942, global = 0.000144678, cumulative = 0.000144678
Min/max rho:1.13988 1.28954
--> FOAM Warning : 
    From function gaussConvectionScheme
    in file finiteVolume/convectionSchemes/gaussConvectionScheme/gaussConvectionScheme.H at line 123
    Reading "/home/zampini/Documenti/personali/Epta/PCMPacks/steady/system/fluid/fvSchemes.divSchemes.div(phi,epsilon)" at line 35
    Unbounded 'Gauss' div scheme used in steady-state solver, use 'bounded Gauss' to ensure boundedness.
    To remove this warning switch off 'boundedGauss' in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/etc/controlDict"
DILUPBiCG:  Solving for epsilon, Initial residual = 0.210492, Final residual = 0.00306992, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 1, Final residual = 0.0467604, No Iterations 1

Solving for solid region domain1
DICPCG:  Solving for h, Initial residual = 1, Final residual = 0.0253284, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 297.999 max(T) [0 0 0 1 0 0 0] 298.003

Solving for solid region domain2
DICPCG:  Solving for h, Initial residual = 1, Final residual = 0.0693808, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 275.845 max(T) [0 0 0 1 0 0 0] 316.145

Solving for solid region domain3

Any idea?

Thanks a lot,
Samuele

[Tobi]

Originally Posted by samiam1000

Hi Tobias,

can I ask you to share the fvSolution and fvSchemes files that you used for the solid regions. I think there's something wrong in mine, since I notice a strange behvior in domain2 region (look at the strange values of temperature) and the simulation is taking a lot of time to solve domain3.

Time = 1


Solving for fluid region fluid
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 0.00545485, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 0.0211023, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 0.00876899, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 1, Final residual = 0.0972692, No Iterations 11
Min/max T:270.103 305.567
DICPCG:  Solving for p_rgh, Initial residual = 0.61935, Final residual = 0.00590249, No Iterations 316
time step continuity errors : sum local = 0.0921942, global = 0.000144678, cumulative = 0.000144678
Min/max rho:1.13988 1.28954
--> FOAM Warning : 
    From function gaussConvectionScheme
    in file finiteVolume/convectionSchemes/gaussConvectionScheme/gaussConvectionScheme.H at line 123
    Reading "/home/zampini/Documenti/personali/Epta/PCMPacks/steady/system/fluid/fvSchemes.divSchemes.div(phi,epsilon)" at line 35
    Unbounded 'Gauss' div scheme used in steady-state solver, use 'bounded Gauss' to ensure boundedness.
    To remove this warning switch off 'boundedGauss' in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/etc/controlDict"
DILUPBiCG:  Solving for epsilon, Initial residual = 0.210492, Final residual = 0.00306992, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 1, Final residual = 0.0467604, No Iterations 1

Solving for solid region domain1
DICPCG:  Solving for h, Initial residual = 1, Final residual = 0.0253284, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 297.999 max(T) [0 0 0 1 0 0 0] 298.003

Solving for solid region domain2
DICPCG:  Solving for h, Initial residual = 1, Final residual = 0.0693808, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 275.845 max(T) [0 0 0 1 0 0 0] 316.145

Solving for solid region domain3

Any idea?

Thanks a lot,
Samuele

Hi samuele,

didn't you notice the things I wrote you in skype. There is a link for the modified case there you find the fvSolutions.

1. your Warning belong to the none used "bound" scheme. Its a numerical trick I always use for stabilisation.

2. You Temperature of your fluid seems to be bad too.
3. Why you are not using only hexaeder meshes? It would be very nice.
4. For pressure equations you should / could use GAMG solvers

Regards
Tobi

[samiam1000]

Dear Tobias,

thanks a lot for answering. Actually, I think that I have not received anything via Skype. Maybe I was online or something went wrong.

Let's try again.
Sorry.

Samuele

[Sherlock_1812]

Hi Samuelle,

Can you tell me how you got rid of this error?

 --> FOAM FATAL IO ERROR: 
keyword type is undefined in dictionary ".T"

file: .T from line 66 to line 75.

I have tried changing the name of the case but it persists.

[samiam1000]

Please,

attach the dictionary!

[esujby]

Hello samiam,

have you resolved ur issue?

best regards