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

TVD scheme Numerical Instability

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

LinkBack Thread Tools Display Modes
Old   June 29, 2011, 01:54
Default TVD scheme Numerical Instability
New Member
Join Date: Jun 2011
Posts: 1
Rep Power: 0
Tempa is on a distinguished road

I'm trying to calculate the following equation which is the derivative in 'x' of a distribution function:
d(dxF)/dt = d(Efield.(dvxF))/dx

dxF = Derivate of Distribution Function respect to "x"
Efield = Electric Field
dvxF = Derivative of Distribution Function respect to "velocity in x"

Im using kinetic model that's why I'm using Distributions functions and Phase-space. Anyway, that's the form of the equation I'm having trouble with

The problem comes because the right hand of the equation is solved by using central difference, but there is a zone where there is a discontinuity in the electric field (It suddenly increases and then goes down again, Its intended to be that way). The electric field is not exactly the same on each point (but more or less the same order of magnitude) but there is a zone where I have in increment so the difference between two adjacent grid is so big that creates an instability. Is there any way to smooth the numerical scheme?

= (E[i+1][j]. dvxF[i+1][j]-E[i-1][j].dvxF[i-1][j])/dx

This is more or less how it goes. Imagine there is a point where the difference between E[i+1][j] and E[i-1][j] is so big that, creating a big gradient in the dxF value.

I've been studying the MUSCL and TVD schemes but i don't quite understand well the procedure
Tempa is offline   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
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
numerical scheme ado Main CFD Forum 3 October 12, 2000 08:20
level set method with TVD scheme? zhanglei Main CFD Forum 1 April 25, 2000 14:01

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