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

Why decreasing time step makes diffusion more significant for diffusion equation?

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

Like Tree1Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 16, 2018, 12:34
Default
  #21
Senior Member
 
Join Date: Oct 2017
Location: United States
Posts: 233
Blog Entries: 1
Rep Power: 9
TurbJet is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
I checked in my old notes and I extracted the LTE expression, it appears clearly that, provided that the dx allows to get the Reynolds cell number <=2, you have a global positive diffusion coefficient.
Now the issue is that if you compare the solution at some time T, when you decreases the time step you need more iteration to reach that time and the error effects summ.
Thx very much for your file !

So to sum up: it appears to me now that the LTE will not increase by simply reducing the time step size; rather, the diffusion-like situation is due to the summation of errors for longer iteration. Am I correct?
TurbJet is offline   Reply With Quote

Old   February 16, 2018, 12:49
Default
  #22
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,813
Rep Power: 72
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by TurbJet View Post
Thx very much for your file !

So to sum up: it appears to me now that the LTE will not increase by simply reducing the time step size; rather, the diffusion-like situation is due to the summation of errors for longer iteration. Am I correct?
Well, that depends also on the value of the Re cell number you used, if it is greater than 2 then it introduces an anti-diffusion effect.
However, the evaluation of the time accuracy makes sense if you fix dx at a value quite low.
FMDenaro is offline   Reply With Quote

Old   February 16, 2018, 12:54
Default
  #23
Senior Member
 
Join Date: Oct 2017
Location: United States
Posts: 233
Blog Entries: 1
Rep Power: 9
TurbJet is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Well, that depends also on the value of the Re cell number you used, if it is greater than 2 then it introduces an anti-diffusion effect.
However, the evaluation of the time accuracy makes sense if you fix dx at a value quite low.
Ok, so anyway, to be simple enough: just make sure reducing dt and dx at the same time, right?
TurbJet is offline   Reply With Quote

Old   February 16, 2018, 13:09
Default
  #24
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,813
Rep Power: 72
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by TurbJet View Post
Ok, so anyway, to be simple enough: just make sure reducing dt and dx at the same time, right?
yes, right
FMDenaro is offline   Reply With Quote

Old   February 16, 2018, 13:09
Default
  #25
Senior Member
 
Join Date: Oct 2017
Location: United States
Posts: 233
Blog Entries: 1
Rep Power: 9
TurbJet is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
yes, right
appreciate it.
FMDenaro likes this.
TurbJet is offline   Reply With Quote

Old   February 17, 2018, 18:20
Default
  #26
Senior Member
 
Join Date: Oct 2017
Location: United States
Posts: 233
Blog Entries: 1
Rep Power: 9
TurbJet is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
yes, right
I re-ran my code with cosine wave instead of Dirac function, and with time step size decreased, the wave just diffuses faster; but eventually, the inf norms of the waves are almost of the same magnitude.

So I guess it is not accurate to say "more diffusion"; it should have be "diffuses slightly faster".
TurbJet is offline   Reply With Quote

Old   February 18, 2018, 01:14
Default
  #27
New Member
 
MOHAMMAD UMAIR
Join Date: Feb 2017
Location: ALIGARH, INDIA
Posts: 3
Rep Power: 9
mohd_umair is on a distinguished road
The initial condition should not have a sudden jump or discontinuities. Try to adjust time-step so that Pe becomes less than 2. Note that Pe<2 doesn't necessarily mean pure diffusion. If doesn't solve solve your problem then reduce dt even further.
mohd_umair is offline   Reply With Quote

Old   February 18, 2018, 03:01
Default
  #28
Senior Member
 
Join Date: Oct 2017
Location: United States
Posts: 233
Blog Entries: 1
Rep Power: 9
TurbJet is on a distinguished road
Quote:
Originally Posted by mohd_umair View Post
The initial condition should not have a sudden jump or discontinuities. Try to adjust time-step so that Pe becomes less than 2. Note that Pe<2 doesn't necessarily mean pure diffusion. If doesn't solve solve your problem then reduce dt even further.
I am a little bit confused by your notation: does your Pe stand for cell Re # or something else?
TurbJet is offline   Reply With Quote

Old   February 18, 2018, 03:05
Default
  #29
Senior Member
 
Join Date: Oct 2017
Location: United States
Posts: 233
Blog Entries: 1
Rep Power: 9
TurbJet is on a distinguished road
Quote:
Originally Posted by mohd_umair View Post
The initial condition should not have a sudden jump or discontinuities. Try to adjust time-step so that Pe becomes less than 2. Note that Pe<2 doesn't necessarily mean pure diffusion. If doesn't solve solve your problem then reduce dt even further.
If you mean by cell Re #, in my case (simple enough),

the velocity a = 1,
diffusivity = 1,
and the mesh size dx = 0.05,

so I guess the cell Re is much smaller than 2, and it should have not appear to have excessive diffusion for this reason.. I guess
TurbJet is offline   Reply With Quote

Old   February 18, 2018, 03:43
Default
  #30
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,813
Rep Power: 72
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by mohd_umair View Post
The initial condition should not have a sudden jump or discontinuities. Try to adjust time-step so that Pe becomes less than 2. Note that Pe<2 doesn't necessarily mean pure diffusion. If doesn't solve solve your problem then reduce dt even further.

The Pe_h<=2 is for the cell number, not for the global number. However, it does not depend on the time step but only on the mesh size.
FMDenaro is offline   Reply With Quote

Old   February 18, 2018, 03:45
Default
  #31
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,813
Rep Power: 72
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by TurbJet View Post
If you mean by cell Re #, in my case (simple enough),

the velocity a = 1,
diffusivity = 1,
and the mesh size dx = 0.05,

so I guess the cell Re is much smaller than 2, and it should have not appear to have excessive diffusion for this reason.. I guess
If you consider the physics you want to resolve, the global Reynolds (o Peclet) number is so small that you have a real physical diffusion that smoothes out your initial solution very fast.
FMDenaro is offline   Reply With Quote

Reply

Tags
diffusion equation, time step size

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
pressure in incompressible solvers e.g. simpleFoam chrizzl OpenFOAM Running, Solving & CFD 13 March 28, 2017 05:49
simpleFoam error - "Floating point exception" mbcx4jc2 OpenFOAM Running, Solving & CFD 12 August 4, 2015 02:20
Help for the small implementation in turbulence model shipman OpenFOAM Programming & Development 25 March 19, 2014 10:08
Micro Scale Pore, icoFoam gooya_kabir OpenFOAM Running, Solving & CFD 2 November 2, 2013 13:58
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33


All times are GMT -4. The time now is 14:57.