CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

numerical scheme

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

LinkBack Thread Tools Display Modes
Old   October 10, 2000, 10:27
Default numerical scheme
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
  Reply With Quote

Old   October 10, 2000, 20:59
Default Re: numerical scheme
Sebastien Perron
Posts: n/a
Have you tried a second order schemes such as Gear?

For this scheme, The approximation of the time derivative is


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.
  Reply With Quote

Old   October 11, 2000, 11:09
Default Re: numerical scheme
Jim Park
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!
  Reply With Quote

Old   October 12, 2000, 08:20
Default Re: numerical scheme
Posts: n/a

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.

  Reply With Quote


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
Numerical viscosity due to the MUSCL and HLL coulpled scheme sonsiest Main CFD Forum 0 May 23, 2011 15:37
How to implement a numerical scheme by using UDF? Furui Xiong FLUENT 0 July 29, 2010 01:30
Numerical scheme in CFX ? sam CFX 1 November 18, 2004 17:25
problem about numerical scheme in LES. libin Main CFD Forum 4 July 1, 2004 04:32
the numerical scheme for LES. John S Main CFD Forum 2 March 14, 2004 08:52

All times are GMT -4. The time now is 04:34.