# TVD Schemes

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

 January 24, 2000, 23:09 TVD Schemes #1 Matt Umbel Guest   Posts: n/a I have been experimenting with using TVD schemes to solve the scalar convection equation given by dphi/dt + U*dphi/dx = 0. I am using 2nd Order Upwind differencing as higher order convection sheme. I am using explicit differencing for temporal discretisation. Eventually I plan to also test Crank-Nicholson and implict. So far I have made the following observations and would appreciate comments/answers. (1) For initial condition that contain sharp discontinuities given by a square wave, the TVD schemes perform, if I may use a subjective description, quite well. (2) For initial condition that is relatively smooth, the TVD schemes tend to sharpen the solution until it almost looks like a square wave, though it does not create any new mins or maxes. It seems for this case that it sharpens the profile, in contrast to first order upwinds drastic smearing. (3) Observation (2) shows a courant number dependancy. For lower courant numbers the TVD scheme does not sharpen the profile as much. (4)Are TVD schemes not valid for problems that do not have sharp discontinuities? (5)If so does anyone have a suggestion for a general-purpose 2nd Order scheme that has flux limiting capabilities? Thanks

 January 25, 2000, 05:21 Re: TVD Schemes #2 Dr. Hrvoje Jasak Guest   Posts: n/a Hi, You have to be a bit careful here! It is true that some TVD schemes sharpen the profile of the scalar if they introduce a certain level of Downwinding. The examples here would be SUPERBEE and van Leer. However, that is NOT the effect you're seeing. If you want to prove this to yourself, you can pick a 2-D test case (like a convection of a profile across a square domain - there is a nice selection of profiles in Leonard's paper). This test case is steady and you will see exactly what one can expect from TVD schemes. What you are seeing is an effect of temporal discretisation. A while back I have managed to derive that explicit time-scheme introduces a negative numerical diffusion (sharpening of the profile) similar and opposite in effect to numerical diffusion of Upwind Differencing (UD). In fact, in 1-D with Co=1 the two effects exactly cancel out and you get the right solution with two first-order schemes! You'll get the opposite effect with implicit temporal discretisation (positive numerical diffusion = smearing of the profile), whereas Crank-Nicholson is diffusion-neutral, but also unbounded (depending on exactly which TVD scheme you choose). If you want the gory details (e.g. derivation of all this), please send me an E-mail and I'll pass it on to you. Hrv

 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 joern OpenFOAM Programming & Development 5 July 7, 2010 15:34 Shuo Main CFD Forum 5 September 19, 2008 21:47 Shuo Main CFD Forum 0 August 19, 2008 21:51 Shuo Main CFD Forum 3 July 3, 2007 22:04 Sudeep FLUENT 0 October 14, 2000 21:31

All times are GMT -4. The time now is 06:45.