# numerical scheme

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 LinkBack Thread Tools Display Modes
 October 10, 2000, 10:27 numerical scheme #1 ado Guest   Posts: n/a Dear friends, i am trying to use numerical methods in order to determine thermal diffusivity and thermal conductivity of unknown material ranged from highly diffused to lower one. The method used is well known and called "parameter estimation" which use the comparison between experimental and numerical temperature-time histories of the material under investigation. I am using a fully implicit numerical scheme. I got very good result in steady state (no variation of temperature with time). This allows to get the thermal conductivity of any material. However, I found something strange. A difference exists between the numerical and the experimental curves in transient state. This difference become more and more important with the "decrease " of thermal diffusivity value. Because of this, it is not possible at the present stage to get a satisfactory results for thermal diffusivity. Personaly, I guess this is due to the numerical scheme i am using which is well known to be unconditionaly stable but not very accurate..But on the other hand, how can this happen only for lower diffuse materials and not for the highers like metal liquids??? waiting for your answers and advise. ado

 October 10, 2000, 20:59 Re: numerical scheme #2 Sebastien Perron Guest   Posts: n/a Have you tried a second order schemes such as Gear? For this scheme, The approximation of the time derivative is (3u^(n+1)-4u^n+u^(n-1))/2dt for a constant time step dt. The other terms are evaluted at time t=t^(n+1) This scheme is also unconditialy stable and might gives better results.

 October 11, 2000, 11:09 Re: numerical scheme #3 Jim Park Guest   Posts: n/a Your accuracy will likely depend on the non-dimensional time step, which is delt(non) = delt(dim)*(alpha)/L^2, where L is a physical dimension, alpha is the thermal diffusivity, and delt(dim) is the dimensional time step. If you keep a constant dimensional time step and mesh, changing the material (diffusivity) changes the accuracy. You should experiment with various (and likely much smaller) time steps until you have a better 'feel' for the relationships between your numerical parameters. This will require running a lot of cases and carefully analyzing the differences in the solutions. Good luck!

 October 12, 2000, 08:20 Re: numerical scheme #4 K.S.Ravichandran Guest   Posts: n/a Ado, 1. I assume that your spatial discretization is 2nd order and time is 1st order. The steady state comparisons then will perhaps produce good results assuming that your linear model is a good approximation of the materails in question, and this is independent of the time step. To make sure that this is the case you must examine the modified equation actually approximated by the differential equation and verify that at steady state all delta_t terms vanish. As Perron has said, a second order in time computation may be a guide to the role played by time step in the accuracy of the computation in the transient stage. As for the error being more pronounced at lower diffusivity values, again the leading order truncation term must be examined to determine if for lower diffusivity, truncationerror magnitude can be held constant by a smaller mesh and time step. I must again empahasise that the mathematiacl model is linear and is good enough to approximate material behavior. Ravi

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post sonsiest Main CFD Forum 0 May 23, 2011 15:37 Furui Xiong FLUENT 0 July 29, 2010 01:30 sam CFX 1 November 18, 2004 17:25 libin Main CFD Forum 4 July 1, 2004 04:32 John S Main CFD Forum 2 March 14, 2004 08:52

All times are GMT -4. The time now is 02:24.

 Contact Us - CFD Online - Top