Problem with JouleHeating simulation in a high aspectratio system
Hello,
I have a problem with a code which I have written to simulate JouleHeating 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 
Hi Farshad,
Rami 
Quote:
Thanks for your answer. I am using our inhouse code which is basically a FiniteDifference code. I think you are right. Because I try to shrink the aspect ratio of the grids in the highaspect ratio case and now it's giving me the right answer. Lz/Lx = 1. Thanks again, Farshad 
Quote:
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 
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 timesteps. The aspect ratio is implicitly dealt with when considering the characteristic length for stability.

All times are GMT 4. The time now is 12:59. 