OneStopGate.Com
OnestopGate   OnestopGate
   Tuesday, December 3, 2024 Login  
OnestopGate
Home | Overview | Syllabus | Tutorials | FAQs | Downloads | Recommended Websites | Advertise | Payments | Contact Us | Forum
OneStopGate

GATE Resources
Gate Articles
Gate Books
Gate Colleges 
Gate Downloads 
Gate Faqs
Gate Jobs
Gate News 
Gate Sample Papers
Training Institutes

GATE Overview
Overview
GATE Eligibility
Structure Of GATE
GATE Coaching Centers
Colleges Providing M.Tech/M.E.
GATE Score
GATE Results
PG with Scholarships
Article On GATE
Admission Process For M.Tech/ MCP-PhD
GATE Topper 2012-13
GATE Forum




GATE 2025 Exclusive
Organizing Institute
Important Dates
How to Apply
Discipline Codes
GATE 2025 Exam Structure

GATE 2025 Syllabus
Aerospace Engg..
Agricultural Engg..
Architecture and Planning
Chemical Engg..
Chemistry
Civil Engg..
Computer Science / IT
Electronics & Communication Engg..
Electrical Engg..
Engineering Sciences
Geology and Geophysics
Instrumentation Engineering
Life Sciences
Mathematics
Mechanical Engg..
Metallurgical Engg..
Mining Engg..
Physics
Production & Industrial Engg..
Pharmaceutical Sciences
Textile Engineering and Fibre Science

GATE Study Material
Aerospace Engg..
Agricultural Engg..
Chemical Engg..
Chemistry
Civil Engg..
Computer Science / IT
Electronics & Communication Engg..
Electrical Engg..
Engineering Sciences
Instrumentation Engg..
Life Sciences
Mathematics
Mechanical Engg..
Physics
Pharmaceutical Sciences
Textile Engineering  and Fibre Science

GATE Preparation
GATE Pattern
GATE Tips N Tricks
Compare Evaluation
Sample Papers 
Gate Downloads 
Experts View

CEED 2013
CEED Exams
Eligibility
Application Forms
Important Dates
Contact Address
Examination Centres
CEED Sample Papers

Discuss GATE
GATE Forum
Exam Cities
Contact Details
Bank Details

Miscellaneous
Advertisment
Contact Us


Home » GATE Study Material » Engineering sciences » Truncation Errors

Truncation Errors

Looking for GATE Preparation Material? Join & Get here now!

** Gate 2013 Question Papers.. ** CEED 2013 Results.. ** Gate 2013 Question Papers With Solutions.. ** GATE 2013 CUT-OFFs.. ** GATE 2013 Results.. **

Truncation Errors

Truncation Errors

Applying a finite difference operator has the side effect of introducing truncation errors into the solution. In this section we discuss how truncation enters the system and how it affects it.  

  • Truncation and Dispersion3.5
  • Consequences and Interpretation

 

Truncation and Dispersion

Earlier we presented a generic equation describing the conservation laws. This allows us to extract information of the numeric properties of the system quite easily, as the analysis of the one generic equation can be propagated to all its special cases. We will not make any assumptions on the drains and sources in the system, and hence do not include any such in this presentation, instead we consider the simple conserved system:

 

\begin{displaymath}\frac{\partial \Gamma }{\partial t}=-\nabla \cdot \Gamma \mathbf{u}.\end{displaymath}
 

Now, to make the analysis easier we initially consider the one dimensional case with constant velocity u,

 

\begin{displaymath}\frac{\partial \Gamma }{\partial t}+u\frac{\partial \Gamma }{\partial x}=0.\end{displaymath}

 

In a previous section we discussed the type and order of the difference operators, in the following we will limit the analysis to a simple centered difference operator, but as we will see the analysis is straight-forward, and may be performed using operators of any order. In continuous time, discrete space, with \( \Gamma _{i} \)designating \( \Gamma \) evaluated at grid node i, we find

 

\begin{displaymath}\frac{\partial \Gamma }{\partial t}+u\frac{\Gamma _{i+1}-\Gamma _{i-1}}{2\Delta x}=0.\end{displaymath}

 

As we have noted earlier, this discretisation implies losing some higher-order terms. To see why this happens we will express the above approximation by the Taylor expansion of the terms involved, as we approximate \( \Gamma _{i+1} \)and \( \Gamma _{i-1} \) by respectively expanding forward and backward from the value \( \Gamma _{i} \) by the distance \( \Delta x \) between any two grid nodes:

 

\begin{eqnarray*}\Gamma _{i+1} & = & \Gamma _{i}+\frac{\Delta x}{1!}\frac{\parti...
...}\frac{\partial ^{3}\Gamma _{i}}{\partial x^{3}}+O(\Delta x^{5})
\end{eqnarray*}

By inserting into the continuous time, discrete space approximation we find the equation we are actually solving

 

\begin{displaymath}\frac{\partial \Gamma }{\partial t}+u\left( \frac{\partial \G...
...ial ^{3}\Gamma _{i}}{\partial x^{3}}+O(\Delta x^{4})\right) =0.\end{displaymath}
 

Consequences and Interpretation

To give an impression of what this means for the solution, we now let \( \Gamma \)represent a simple sinusoidal wave in time and space. We then isolate the temporal frequency and calculate the wave speed, from which we can gain some insight into how the system evolves. The wave can be described as \( \Gamma =e^{ikx-i\omega t} \)which enables us to extend the argument, at least to periodic solutions as we simply approximate the shape through Fourier series. By simple insertion we find

 

\begin{eqnarray*}& \frac{\partial \Gamma }{\partial t}+u\left( \frac{\partial \G...
...^{2}}{3!}k^{2}+...+(-1)^{n}\frac{\Delta x^{n}}{n!}k^{n}+...). &
\end{eqnarray*}

 


With spatial wave length \( l=\frac{2\pi }{k} \) and time frequency \( T=\frac{2\pi }{\omega } \)we find the wave speed
 

 

\begin{displaymath}s=\frac{l}{T}=\frac{\omega }{k}=u(1-\frac{(k\Delta x)^{2}}{3!}...+(-1)^{n}\frac{(k\Delta x)^{n}}{n!}+...),\end{displaymath}

 

 

which converges as the faculty operator is stronger than the power. We see that the speed at which a solution propagates is dependent on the spatial frequency, and that the most significant term enters with a negative coefficient, causing a slow down of frequencies. At constant grid spacing \( \Delta x \), and as a low value of k models a high spatial frequency wave, we can now see high frequency modes will move slower than low frequency ones. This means that over time a solution consisting of a broad range of frequencies will change shape, as high mode frequencies will be left behind. This is especially evident in the solution of the square wave, which has non-vanishing amplitude in all frequency modes, and hence all terms in the Taylor expansion are of importance. This is the origin of the dispersion phenomena mentioned in an earlier section. In our case this is especially worrying, as we have already seen, shocks may evolve naturally in the system, and hence the near discontinuities in the shock fronts result in high order frequencies dominating the system. This can be seen in figure 3.5
 
   Figure 3.5: Constant velocity propagation of a steep wave using a naive implementation of the difference operator.
\resizebox*{0.8\textwidth}{0.2\textheight}{\includegraphics{nodiff.eps}}

 

 


, where we see how such a system evolves. We see that the truncation results in dispersion as high frequency modes are left behind by the bulk of the wave causing numerical instability due to the steep gradients flooding the solution. At roughly iteration number 130 the valley immediately following the hill became negative, crashing the simulation. As a rule of thumb truncation is of no immediate concern if the shock is resolved on the grid, i.e. the wave length of the highest mode frequency in the shock is a number of times longer than the grid spacing \( \Delta x \), as is also hinted by the term \( k\Delta x \) entering under the power operator.

Another view at the problem is noticing that in the difference operation of order n we truncated
 

\begin{displaymath}\frac{\Delta x^{n}}{n!}\frac{\partial ^{n}\Gamma _{i}}{\partial x^{n}}\end{displaymath}

 

terms from the Taylor series from the order n and up, terms for which the coefficients \( \Delta x^{n}/n! \) was in some sense small, under the assumption that the differential operator \( \partial ^{n}\Gamma _{i}/\partial x^{n} \)is in some sense limited, making all the truncated terms insignificant for the solution. But at near discontinuities this is not the case, as \( \partial ^{n}\Gamma _{i}/\partial x^{n} \)becomes of significance, breaking the assumption and infecting the solution.

The net result is that the truncation error will be of relevance in shock regions, and that we cannot approximate the solution correctly in an extreme shock-region, no matter the order of the difference operator. If ignored this will in many cases cause the discrete approximations to explode, and in return drown the correct solution in numerical noise.

But what if we could avoid the discontinuities altogether?

 

 
 



Discussion Center

Discuss/
Query

Papers/
Syllabus

Feedback/
Suggestion

Yahoo
Groups

Sirfdosti
Groups

Contact
Us

MEMBERS LOGIN
  
Email ID:
Password:

  Forgot Password?
 New User? Register!

INTERVIEW EBOOK
Get 9,000+ Interview Questions & Answers in an eBook. Interview Question & Answer Guide
  • 9,000+ Interview Questions
  • All Questions Answered
  • 5 FREE Bonuses
  • Free Upgrades
GATE RESOURCES
 
  • Gate Books
  • Training Institutes
  • Gate FAQs
  • GATE BOOKS
     
  • Mechanical Engineeering Books
  • Robotics Automations Engineering Books
  • Civil Engineering Books
  • Chemical Engineering Books
  • Environmental Engineering Books
  • Electrical Engineering Books
  • Electronics Engineering Books
  • Information Technology Books
  • Software Engineering Books
  • GATE Preparation Books
  • Exciting Offers



    GATE Exam, Gate 2009, Gate Papers, Gate Preparation & Related Pages


    GATE Overview | GATE Eligibility | Structure Of GATE | GATE Training Institutes | Colleges Providing M.Tech/M.E. | GATE Score | GATE Results | PG with Scholarships | Article On GATE | GATE Forum | GATE 2009 Exclusive | GATE 2009 Syllabus | GATE Organizing Institute | Important Dates for GATE Exam | How to Apply for GATE | Discipline / Branch Codes | GATE Syllabus for Aerospace Engineering | GATE Syllabus for Agricultural Engineering | GATE Syllabus for Architecture and Planning | GATE Syllabus for Chemical Engineering | GATE Syllabus for Chemistry | GATE Syllabus for Civil Engineering | GATE Syllabus for Computer Science / IT | GATE Syllabus for Electronics and Communication Engineering | GATE Syllabus for Engineering Sciences | GATE Syllabus for Geology and Geophysics | GATE Syllabus for Instrumentation Engineering | GATE Syllabus for Life Sciences | GATE Syllabus for Mathematics | GATE Syllabus for Mechanical Engineering | GATE Syllabus for Metallurgical Engineering | GATE Syllabus for Mining Engineering | GATE Syllabus for Physics | GATE Syllabus for Production and Industrial Engineering | GATE Syllabus for Pharmaceutical Sciences | GATE Syllabus for Textile Engineering and Fibre Science | GATE Preparation | GATE Pattern | GATE Tips & Tricks | GATE Compare Evaluation | GATE Sample Papers | GATE Downloads | Experts View on GATE | CEED 2009 | CEED 2009 Exam | Eligibility for CEED Exam | Application forms of CEED Exam | Important Dates of CEED Exam | Contact Address for CEED Exam | CEED Examination Centres | CEED Sample Papers | Discuss GATE | GATE Forum of OneStopGATE.com | GATE Exam Cities | Contact Details for GATE | Bank Details for GATE | GATE Miscellaneous Info | GATE FAQs | Advertisement on GATE | Contact Us on OneStopGATE |
    Copyright © 2024. One Stop Gate.com. All rights reserved Testimonials |Link To Us |Sitemap |Privacy Policy | Terms and Conditions|About Us
    Our Portals : Academic Tutorials | Best eBooksworld | Beyond Stats | City Details | Interview Questions | India Job Forum | Excellent Mobiles | Free Bangalore | Give Me The Code | Gog Logo | Free Classifieds | Jobs Assist | Interview Questions | One Stop FAQs | One Stop GATE | One Stop GRE | One Stop IAS | One Stop MBA | One Stop SAP | One Stop Testing | Web Hosting | Quick Site Kit | Sirf Dosti | Source Codes World | Tasty Food | Tech Archive | Software Testing Interview Questions | Free Online Exams | The Galz | Top Masala | Vyom | Vyom eBooks | Vyom International | Vyom Links | Vyoms | Vyom World
    C Interview Questions | C++ Interview Questions | Send Free SMS | Placement Papers | SMS Jokes | Cool Forwards | Romantic Shayari