Problem with Joule-Heating simulation in a high aspect-ratio system

farshadsm · November 28, 2012, 15:27

Hello,

I have a problem with a code which I have written to simulate Joule-Heating through a 3D cubic system made of silver. When the aspect ratio between l_x/l_y (l_x and l_y are the length of the system in x and y direction) is less than 100, I get good results from the code independent of the number of grid points generated in the system. However, when the l_x/l_y is larger than 100, the temperature distribution predicted by the model shows deviations from a true symmetric distribution for even 30000 grid points which are not really large. What could be the problem? I have used a resistive networks to solve electrical and heat conduction equations.

I appreciate for any help.

Farshad

Rami · November 29, 2012, 03:28

Hi Farshad,

Are you using a CFD package, and if so - which?
With a large aspect ratio, did you also use high aspect ration cells? If so - it can explain the difficultiy you are experiencing.
You did not mention the Lz / Lx ratio. Please specify.

Rami

farshadsm · November 29, 2012, 06:57

Quote:

Originally Posted by Rami

Hi Farshad,

Are you using a CFD package, and if so - which?
With a large aspect ratio, did you also use high aspect ration cells? If so - it can explain the difficultiy you are experiencing.
You did not mention the Lz / Lx ratio. Please specify.

Rami

Hi Rami,

Thanks for your answer. I am using our in-house code which is basically a Finite-Difference code. I think you are right. Because I try to shrink the aspect ratio of the grids in the high-aspect ratio case and now it's giving me the right answer. Lz/Lx = 1.

Thanks again,

Farshad

farshadsm · November 30, 2012, 12:00

Quote:

Originally Posted by farshadsm

Hi Rami,

Thanks for your answer. I am using our in-house code which is basically a Finite-Difference code. I think you are right. Because I try to shrink the aspect ratio of the grids in the high-aspect ratio case and now it's giving me the right answer. Lz/Lx = 1.

Thanks again,

Farshad

Rami,

Is there any mathematical rule or basis for choosing the correct value for the time step as well as grid boxes for large geometrical aspect ratio?

Farshad

Rami · December 4, 2012, 03:09

The time step should be selected by considering stability and accuracy. The stability depends on the scheme you use, and - if your boundary conditions, initial conditions and sources are extremely time- or space- dependent (e.g., have discontinuities of high derivatives) - they may also influence stability. Accuracy is usually improved by smaller time-steps. The aspect ratio is implicitly dealt with when considering the characteristic length for stability.

November 28, 2012, 15:27	Problem with Joule-Heating simulation in a high aspect-ratio system	#1
farshadsm New Member Join Date: Nov 2012 Posts: 3 Rep Power: 13	Hello, I have a problem with a code which I have written to simulate Joule-Heating through a 3D cubic system made of silver. When the aspect ratio between l_x/l_y (l_x and l_y are the length of the system in x and y direction) is less than 100, I get good results from the code independent of the number of grid points generated in the system. However, when the l_x/l_y is larger than 100, the temperature distribution predicted by the model shows deviations from a true symmetric distribution for even 30000 grid points which are not really large. What could be the problem? I have used a resistive networks to solve electrical and heat conduction equations. I appreciate for any help. Farshad

November 29, 2012, 03:28		#2
Rami Senior Member Rami Ben-Zvi Join Date: Mar 2009 Posts: 155 Rep Power: 17	Hi Farshad, Are you using a CFD package, and if so - which? With a large aspect ratio, did you also use high aspect ration cells? If so - it can explain the difficultiy you are experiencing. You did not mention the Lz / Lx ratio. Please specify. Rami

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[Commercial meshers] [fluent3DMeshToFoam] High aspect ratio cells found (mesh from ICEM)	alquimista	OpenFOAM Meshing & Mesh Conversion	10	January 29, 2013 15:28
use of high aspect ratio cells near walls	Usman	FLUENT	0	January 23, 2008 07:16
Problem on high density ratio in Level Set method	Kai Yan	Main CFD Forum	10	December 25, 2007 06:12
mesh for a high aspect ratio geometry	matthias	CFX	3	October 20, 2006 02:55
high aspect ratio elements	Olivier	CFX	5	April 3, 2003 06:55

December 4, 2012, 03:09		#5
Rami Senior Member Rami Ben-Zvi Join Date: Mar 2009 Posts: 155 Rep Power: 17	The time step should be selected by considering stability and accuracy. The stability depends on the scheme you use, and - if your boundary conditions, initial conditions and sources are extremely time- or space- dependent (e.g., have discontinuities of high derivatives) - they may also influence stability. Accuracy is usually improved by smaller time-steps. The aspect ratio is implicitly dealt with when considering the characteristic length for stability.