baran_foam |
August 31, 2015 06:02 |
1 Attachment(s)
Greeting all!!!
This is mesh quality of my geometry...I hope this is not made any problem...
Code:
Create polyMesh for time = 0
Enabling all (cell, face, edge, point) topology checks.
Enabling all geometry checks.
Time = 0
Mesh stats
points: 1116300
faces: 3270406
internal faces: 3191462
cells: 1076978
faces per cell: 6
boundary patches: 10
point zones: 0
face zones: 13
cell zones: 4
Overall number of cells of each type:
hexahedra: 1076978
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 0
Checking topology...
Boundary definition OK.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Topological cell zip-up check OK.
Face-face connectivity OK.
Number of regions: 1 (OK).
Checking patch topology for multiply connected surfaces...
Patch Faces Points Surface topology Bounding box
he_surface 13720 13736 ok (closed singly connected) (0.247 0.251 0.272) (0.917 0.941 1.202)
w_right_side_ref 5952 6111 ok (non-closed singly connected) (2.6999174e-24 0.15 0.212) (5.2351103e-16 1.036 1.352)
w_left_side_ref 5952 6111 ok (non-closed singly connected) (1.164 0.15 0.212) (1.164 1.036 1.352)
w_back_ref 12288 12513 ok (non-closed singly connected) (0 0.15 1.352) (1.164 1.036 1.502)
w_brick_bot 3256 3500 ok (non-closed singly connected) (0.217 -9.5154073e-16 0.212) (0.947 9.5707519e-16 1.352)
w_bot_ref 7104 7499 ok (non-closed singly connected) (0 -1.1571817e-15 0.212) (1.164 0.15 1.502)
w_door1 2588 2676 ok (non-closed singly connected) (0.332 0.356 -0.1) (0.832 0.85600003 -0.05)
w_door2 696 870 ok (non-closed singly connected) (0.332 0.356 -0.05) (0.832 0.85600003 0)
w_bricks 17028 17343 ok (non-closed singly connected) (0 -1.1571817e-15 0) (1.164 1.186 0.212)
w_top_ref 10360 10553 ok (non-closed singly connected) (0 1.036 0.212) (1.164 1.186 1.502)
Checking geometry...
Overall domain bounding box (0 -1.1571817e-15 -0.1) (1.164 1.186 1.502)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (3.2928793e-15 2.2620737e-15 7.0567305e-15) OK.
Max cell openness = 3.3154974e-16 OK.
Max aspect ratio = 13.162501 OK.
Minimum face area = 6.5776377e-06. Maximum face area = 0.00064648215. Face area magnitudes OK.
Min volume = 8.3282847e-08. Max volume = 7.9145335e-06. Total volume = 2.0928022. Cell volumes OK.
Mesh non-orthogonality Max: 69.14912 average: 8.4597446
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.66268921 OK.
Coupled point location match (average 0) OK.
Face tets OK.
Min/max edge length = 0.0016666666 0.032313618 OK.
All angles in faces OK.
Face flatness (1 = flat, 0 = butterfly) : min = 0.99388226 average = 0.99999605
All face flatness OK.
Cell determinant (wellposedness) : minimum: 0.028461416 average: 4.6877314
Cell determinant check OK.
Concave cell check OK.
Face interpolation weight : minimum: 0.10005566 average: 0.49007212
Face interpolation weight check OK.
Face volume ratio : minimum: 0.11131699 average: 0.9662018
Face volume ratio check OK.
Mesh OK.
Using the convection bc I have get my solution right though using radiation bc case is diverging as I told in my first place..Apart from this temperature variation with time is coming random with radiation bc....I have attached two point is air for this two case.. in picture(in total four...)
Can you suggest what might be the reason??
From the discussion of our previous case it is observed that case is diverging for radiation case...in bricks zone(thermal conductivity-0.15 J/kg K)
So I change the thermal conductivity of bricks from 0.15 to 10. Then run the same case with radiation bc....
Then case is diverging in other zone..fibrothal(thermal conductivity is around 0.15 J/kg K)...So is there is any problem with this radiation bc or radiation to handle high T with small thermal conductivity of insulation..
Please suggest something...mostly this problem is not involve with any mesh issue...??? I have attached the log file for the later case below...
Code:
Region: fluid Courant Number mean: 0.0198761 max: 0.29878319
Region: fibrothal_in Diffusion Number mean: 5.3692836e-08 max: 3.868369e-07
Region: solid_bricks Diffusion Number mean: 2.453461e-06 max: 1.317393e-05
Region: unifrax_out Diffusion Number mean: 1.0442085e-07 max: 3.2171136e-07
deltaT = 0.0011417394
Time = 132.83673
Solving for fluid region fluid
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 7.8428474e-05, Final residual = 2.3160953e-08, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 5.9872307e-05, Final residual = 2.1848024e-08, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 8.1347697e-05, Final residual = 1.7257565e-08, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 1.3754107e-05, Final residual = 1.8032295e-09, No Iterations 1
Min/max T:306.56532 939
GAMG: Solving for p_rgh, Initial residual = 0.004153507, Final residual = 8.7413984e-08, No Iterations 10
GAMG: Solving for p_rgh, Initial residual = 8.9189776e-06, Final residual = 8.0399224e-08, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 2.1315608e-06, Final residual = 7.7674832e-08, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (fluid): sum local = 6.4079763e-11, global = 8.4405831e-13, cumulative = 1.0816831e-07
GAMG: Solving for p_rgh, Initial residual = 7.9826171e-07, Final residual = 3.8468735e-08, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 2.8369056e-07, Final residual = 4.4607049e-08, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 1.209402e-07, Final residual = 2.6745305e-08, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (fluid): sum local = 2.2064599e-11, global = 8.3155456e-15, cumulative = 1.0816832e-07
Solving for solid region fibrothal_in
DICPCG: Solving for h, Initial residual = 4.8433956e-06, Final residual = 5.6555918e-12, No Iterations 1
DICPCG: Solving for h, Initial residual = 5.846885e-12, Final residual = 5.846885e-12, No Iterations 0
Min/max T:min(T) [0 0 0 1 0 0 0] 300 max(T) [0 0 0 1 0 0 0] 823.55181
Solving for solid region solid_bricks
DICPCG: Solving for h, Initial residual = 5.3780627e-06, Final residual = 1.8628724e-10, No Iterations 1
DICPCG: Solving for h, Initial residual = 1.9739328e-10, Final residual = 1.9739328e-10, No Iterations 0
Min/max T:min(T) [0 0 0 1 0 0 0] 300.00584 max(T) [0 0 0 1 0 0 0] 390.0248
Solving for solid region unifrax_out
DICPCG: Solving for h, Initial residual = 3.1120483e-06, Final residual = 3.4748384e-12, No Iterations 1
DICPCG: Solving for h, Initial residual = 3.5734963e-12, Final residual = 3.5734963e-12, No Iterations 0
Min/max T:min(T) [0 0 0 1 0 0 0] 300 max(T) [0 0 0 1 0 0 0] 317.42316
ExecutionTime = 48515.92 s ClockTime = 49698 s
Region: fluid Courant Number mean: 0.019876 max: 0.29877615
Region: fibrothal_in Diffusion Number mean: 5.3692908e-08 max: 3.8683897e-07
Region: solid_bricks Diffusion Number mean: 2.453461e-06 max: 1.317393e-05
Region: unifrax_out Diffusion Number mean: 1.0442085e-07 max: 3.2171136e-07
deltaT = 0.0011417394
Time = 132.83787
Solving for fluid region fluid
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 7.8490734e-05, Final residual = 2.314818e-08, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 5.9735226e-05, Final residual = 2.1923582e-08, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 8.1540717e-05, Final residual = 1.7182041e-08, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 1.3757046e-05, Final residual = 1.804325e-09, No Iterations 1
LU Back substitute C matrix..
Min/max T:306.56544 939
GAMG: Solving for p_rgh, Initial residual = 0.0041498473, Final residual = 9.054912e-08, No Iterations 9
GAMG: Solving for p_rgh, Initial residual = 8.3165628e-06, Final residual = 7.0594986e-08, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 2.032543e-06, Final residual = 6.9336629e-08, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (fluid): sum local = 5.7201248e-11, global = -9.3147623e-13, cumulative = 1.0816739e-07
GAMG: Solving for p_rgh, Initial residual = 7.6577489e-07, Final residual = 9.7339006e-08, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 2.8026848e-07, Final residual = 5.1413822e-08, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 1.298081e-07, Final residual = 2.9943948e-08, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (fluid): sum local = 2.4703296e-11, global = -6.8827074e-13, cumulative = 1.081667e-07
Solving for solid region fibrothal_in
DICPCG: Solving for h, Initial residual = 4.8426671e-06, Final residual = 5.6623743e-12, No Iterations 1
DICPCG: Solving for h, Initial residual = 4.4247984e-11, Final residual = 4.4247984e-11, No Iterations 0
Min/max T:min(T) [0 0 0 1 0 0 0] 168.02419 max(T) [0 0 0 1 0 0 0] 823.54449
Solving for solid region solid_bricks
DICPCG: Solving for h, Initial residual = 5.3882983e-06, Final residual = 1.8640842e-10, No Iterations 1
DICPCG: Solving for h, Initial residual = 3.0063655e-10, Final residual = 3.0063655e-10, No Iterations 0
Min/max T:min(T) [0 0 0 1 0 0 0] 300.00584 max(T) [0 0 0 1 0 0 0] 390.07683
Solving for solid region unifrax_out
DICPCG: Solving for h, Initial residual = 3.1120344e-06, Final residual = 3.4748206e-12, No Iterations 1
DICPCG: Solving for h, Initial residual = 3.5760929e-12, Final residual = 3.5760929e-12, No Iterations 0
Min/max T:min(T) [0 0 0 1 0 0 0] 300 max(T) [0 0 0 1 0 0 0] 317.42353
ExecutionTime = 48516.76 s ClockTime = 49698 s
Region: fluid Courant Number mean: 0.019875904 max: 0.29876907
Region: fibrothal_in Diffusion Number mean: 5.369298e-08 max: 3.8684106e-07
Region: solid_bricks Diffusion Number mean: 2.453461e-06 max: 1.317393e-05
Region: unifrax_out Diffusion Number mean: 1.0442085e-07 max: 3.2171136e-07
deltaT = 0.0011417394
Time = 132.83901
Solving for fluid region fluid
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 7.8363917e-05, Final residual = 2.3111003e-08, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 5.9808245e-05, Final residual = 2.181205e-08, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 8.1303226e-05, Final residual = 1.7190079e-08, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 1.3769818e-05, Final residual = 1.8025345e-09, No Iterations 1
Min/max T:-34.638404 939
GAMG: Solving for p_rgh, Initial residual = 0.0041635193, Final residual = 8.7220279e-08, No Iterations 18
GAMG: Solving for p_rgh, Initial residual = 3.0025266e-05, Final residual = 8.3696645e-08, No Iterations 7
GAMG: Solving for p_rgh, Initial residual = 6.9780425e-06, Final residual = 9.4629866e-08, No Iterations 3
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (fluid): sum local = 7.8066502e-11, global = -1.0715968e-11, cumulative = 1.0815598e-07
GAMG: Solving for p_rgh, Initial residual = 2.6392854e-06, Final residual = 8.7769725e-08, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 9.0381004e-07, Final residual = 4.9359257e-08, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 3.9857535e-07, Final residual = 5.3627197e-08, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (fluid): sum local = 4.4241536e-11, global = 2.0258428e-12, cumulative = 1.0815801e-07
Solving for solid region fibrothal_in
DICPCG: Solving for h, Initial residual = 4.8427064e-06, Final residual = 5.6629409e-12, No Iterations 1
DICPCG: Solving for h, Initial residual = 1.9117983e-11, Final residual = 1.9117983e-11, No Iterations 0
Min/max T:min(T) [0 0 0 1 0 0 0] -34.64088 max(T) [0 0 0 1 0 0 0] 823.55061
Solving for solid region solid_bricks
DICPCG: Solving for h, Initial residual = 5.3881957e-06, Final residual = 1.8640995e-10, No Iterations 1
DICPCG: Solving for h, Initial residual = 1.9768437e-10, Final residual = 1.9768437e-10, No Iterations 0
Min/max T:min(T) [0 0 0 1 0 0 0] 300.00584 max(T) [0 0 0 1 0 0 0] 390.07814
Solving for solid region unifrax_out
DICPCG: Solving for h, Initial residual = 3.1120205e-06, Final residual = 3.4748029e-12, No Iterations 1
DICPCG: Solving for h, Initial residual = 3.5713957e-12, Final residual = 3.5713957e-12, No Iterations 0
Min/max T:min(T) [0 0 0 1 0 0 0] 300 max(T) [0 0 0 1 0 0 0] 317.42391
ExecutionTime = 48517.36 s ClockTime = 49699 s
Region: fluid Courant Number mean: 0.019875807 max: 0.29875958
Region: fibrothal_in Diffusion Number mean: 5.3693051e-08 max: 3.8684315e-07
Region: solid_bricks Diffusion Number mean: 2.453461e-06 max: 1.317393e-05
Region: unifrax_out Diffusion Number mean: 1.0442085e-07 max: 3.2171136e-07
deltaT = 0.0011417394
Time = 132.84016
Solving for fluid region fluid
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 8.1108496e-05, Final residual = 2.5818077e-08, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 6.134127e-05, Final residual = 2.4355749e-08, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 8.5046122e-05, Final residual = 1.9803553e-08, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 1.3767232e-05, Final residual = 1.8025597e-09, No Iterations 1
Min/max T:168.04889 939
GAMG: Solving for p_rgh, Initial residual = 0.0041488398, Final residual = 9.7922006e-08, No Iterations 16
GAMG: Solving for p_rgh, Initial residual = 2.5548778e-05, Final residual = 8.2376033e-08, No Iterations 7
GAMG: Solving for p_rgh, Initial residual = 6.7243635e-06, Final residual = 5.251697e-08, No Iterations 4
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (fluid): sum local = 4.3325191e-11, global = 5.3072732e-12, cumulative = 1.0816332e-07
GAMG: Solving for p_rgh, Initial residual = 1.8993958e-06, Final residual = 6.9772195e-08, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 6.1908036e-07, Final residual = 8.0077444e-08, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 2.4772369e-07, Final residual = 4.5438216e-08, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (fluid): sum local = 3.7486014e-11, global = 1.2021921e-12, cumulative = 1.0816452e-07
Solving for solid region unifrax_out
DICPCG: Solving for h, Initial residual = 3.1120066e-06, Final residual = 3.4747852e-12, No Iterations 1
DICPCG: Solving for h, Initial residual = 3.5740577e-12, Final residual = 3.5740577e-12, No Iterations 0
Min/max T:min(T) [0 0 0 1 0 0 0] 300 max(T) [0 0 0 1 0 0 0] 317.42429
ExecutionTime = 48517.93 s ClockTime = 49700 s
Region: fluid Courant Number mean: 0.019875723 max: 0.29875153
Region: fibrothal_in Diffusion Number mean: 5.3693123e-08 max: 3.8684525e-07
Region: solid_bricks Diffusion Number mean: 2.453461e-06 max: 1.317393e-05
Region: unifrax_out Diffusion Number mean: 1.0442085e-07 max: 3.2171136e-07
deltaT = 0.0011417394
Time = 132.8413
Solving for fluid region fluid
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 7.8970139e-05, Final residual = 2.3571827e-08, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 6.0058436e-05, Final residual = 2.2336877e-08, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 8.2795778e-05, Final residual = 1.7421106e-08, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 1.3770425e-05, Final residual = 1.8043668e-09, No Iterations 1
Min/max T:-35.535871 939
GAMG: Solving for p_rgh, Initial residual = 0.0041525003, Final residual = 8.3978618e-08, No Iterations 16
GAMG: Solving for p_rgh, Initial residual = 2.0716891e-05, Final residual = 8.414588e-08, No Iterations 6
GAMG: Solving for p_rgh, Initial residual = 5.9691193e-06, Final residual = 6.361722e-08, No Iterations 4
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (fluid): sum local = 5.2482977e-11, global = -1.8655448e-12, cumulative = 1.0816265e-07
GAMG: Solving for p_rgh, Initial residual = 1.6250719e-06, Final residual = 7.813847e-08, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 4.7865432e-07, Final residual = 6.8283963e-08, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 1.8631846e-07, Final residual = 3.5284387e-08, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (fluid): sum local = 2.9109187e-11, global = -2.3103795e-13, cumulative = 1.0816242e-07
Solving for solid region fibrothal_in
DICPCG: Solving for h, Initial residual = 4.8424178e-06, Final residual = 5.6627448e-12, No Iterations 1
DICPCG: Solving for h, Initial residual = 7.5644313e-12, Final residual = 7.5644313e-12, No Iterations 0
Min/max T:min(T) [0 0 0 1 0 0 0] -35.538367 max(T) [0 0 0 1 0 0 0] 823.55579
Solving for solid region solid_bricks
DICPCG: Solving for h, Initial residual = 5.3879905e-06, Final residual = 1.864035e-10, No Iterations 1
DICPCG: Solving for h, Initial residual = 1.9767345e-10, Final residual = 1.9767345e-10, No Iterations 0
Min/max T:min(T) [0 0 0 1 0 0 0] 300.00584 max(T) [0 0 0 1 0 0 0] 390.08076
Solving for solid region unifrax_out
DICPCG: Solving for h, Initial residual = 3.1119927e-06, Final residual = 3.4747675e-12, No Iterations 1
DICPCG: Solving for h, Initial residual = 3.5731187e-12, Final residual = 3.5731187e-12, No Iterations 0
Min/max T:min(T) [0 0 0 1 0 0 0] 300 max(T) [0 0 0 1 0 0 0] 317.42467
ExecutionTime = 48518.51 s ClockTime = 49700 s
Region: fluid Courant Number mean: 0.019875633 max: 0.29874413
Region: fibrothal_in Diffusion Number mean: 5.3693195e-08 max: 3.8684734e-07
Region: solid_bricks Diffusion Number mean: 2.453461e-06 max: 1.317393e-05
Region: unifrax_out Diffusion Number mean: 1.0442085e-07 max: 3.2171136e-07
deltaT = 0.0011417394
Regards,
Baran
|