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

Why 3D solid-pore geometry showing diverged solution?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 3, 2012, 15:45
Default Why 3D solid-pore geometry showing diverged solution?
  #1
Member
 
Sangeeta
Join Date: Jul 2012
Location: Kingston, Canada
Posts: 70
Rep Power: 4
Sargam05 is on a distinguished road
Hello everyone,

I am using solidDisplacement solver for 3D solid-pore cube geometry. Following are displacement D boundary conditions:

boundaryField
{
electron_minX
{
type tractionDisplacement;
traction uniform ( 0 0 0 );
pressure uniform 0;
value uniform (0 0 0);


}

electron_maxX
{
type tractionDisplacement;
traction uniform ( 0 0 0 );
pressure uniform 0;
value uniform (0 0 0);

}
electron_minY
{
type tractionDisplacement;
traction uniform ( 0 10 0 );
pressure uniform 0;
value uniform (0 0 0);



}
electron_maxY
{
type tractionDisplacement;
traction uniform ( 0 -10 0 );
pressure uniform 0;
value uniform (0 0 0);


}
electron_minZ
{
type tractionDisplacement;
traction uniform ( 0 0 0);
pressure uniform 0;
value uniform (0 0 0);


}
electron_maxZ
{
type tractionDisplacement;
traction uniform ( 0 0 0 );
pressure uniform 0;
value uniform (0 0 0);

}
electron_to_pore
{
type tractionDisplacement;
traction uniform ( 0 0 0 );
pressure uniform 0;
value uniform (0 0 0);

}
electron_to_npore
{
type zeroGradient;
}



}


OpenFOAM simulation is giving diverge solution. Following is part of log file:
GAMG: Solving for Dx, Initial residual = 0.0407966, Final residual = 0.00378547, No Iterations 1
GAMG: Solving for Dy, Initial residual = 0.0521521, Final residual = 0.00500495, No Iterations 2
GAMG: Solving for Dz, Initial residual = 0.0856927, Final residual = 0.00255286, No Iterations 2
Min/Max/Ave DEqn-7.00897e+65 -2.76576e+77 -1.2423e+65) (7.19343e+65 -2.76576e+77 1.17129e+65) (1.36111e+64 -2.76576e+77 -2.05129e+64)
GAMG: Solving for Dx, Initial residual = 0.030862, Final residual = 0.0010546, No Iterations 2
GAMG: Solving for Dy, Initial residual = 0.0446989, Final residual = 0.00435236, No Iterations 3
GAMG: Solving for Dz, Initial residual = 0.069533, Final residual = 0.00225699, No Iterations 2
Min/Max/Ave DEqn-6.97373e+65 -2.76576e+77 -1.2305e+65) (7.17125e+65 -2.76576e+77 1.07291e+65) (1.47816e+64 -2.76576e+77 -2.09031e+64)
GAMG: Solving for Dx, Initial residual = 0.0282657, Final residual = 0.000732531, No Iterations 2
GAMG: Solving for Dy, Initial residual = 0.0351312, Final residual = nan, No Iterations 1000
GAMG: Solving for Dz, Initial residual = 0.0653158, Final residual = 0.00194052, No Iterations 2
Min/Max/Ave DEqn-6.92125e+65 nan -1.19559e+65) (7.11257e+65 nan 1.03401e+65) (1.54778e+64 nan -2.09798e+64)


I have done similar 3D cube geometry problem without pore in which got converge solution. But I do not know why similar boundary conditions are not giving converge solution for solid-pore 3D cube geometry. I have also tried OpenFOAM-ext but similar problem is coming. When I am using three symmetry planes I get converge solution but when I am taking all surfaces as solution domain I get diverge solution. How I can get converge solution with considering all surfaces of the 3D cube (solid-pore phase) geometry?

Does someone has any idea why this problem is coming? How I can solve this problem?

Best regards,
Sangeeta
Sargam05 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
CFD in ProE Lev Main CFD Forum 19 September 26, 2001 07:29
CFL Condition Matt Umbel Main CFD Forum 14 January 12, 2001 14:34
Adaptive optimization of geometry Jan Rusås CFX 1 September 16, 2000 02:38
Mesh for 3 dim Geometry Phil FLUENT 9 July 12, 2000 04:39
Wall functions Abhijit Tilak Main CFD Forum 6 February 5, 1999 01:16


All times are GMT -4. The time now is 03:25.