My Math Forum  

Go Back   My Math Forum > Science Forums > Computer Science

Computer Science Computer Science Forum


Reply
 
LinkBack Thread Tools Display Modes
March 20th, 2011, 03:31 AM   #1
Newbie
 
Joined: Mar 2011

Posts: 6
Thanks: 0

predictor-corrector scheme

Hi!
I'm asking for help to understand the realization of the iterative predictor-corrector scheme (I'm not a specialist in this field) for PDE system . The description of a method attached in article.
The question concerns algorithm realization under the formula (A5), resulted in article
(A5)
where - function in the right part of the equations, and - variables.
, as I have understood, these are predicted value of all variables from, for example, Euler's method. Then it is necessary to calculate variables and to compare them with predicted, repeating algorithm before achievement of necessary accuracy. But how to calculate if in the right part there are the same , for me not clearly? In any way I don't understand the dependence $ of two variables means? How to use this formula in calculations? I will be glad to any concrete algorithms or the references to it.
Wait for your suggestions.
Attached Images
File Type: tif FDTD_PC.tif (155.5 KB, 1089 views)
-nad- is offline  
 
March 20th, 2011, 08:53 PM   #2
Member
 
Joined: Nov 2007

Posts: 73
Thanks: 0

Re: predictor-corrector scheme

Without more references (I can't see the attachment), I guess the corrector method is implicit. Roughly speaking, you have both sides. You can solve the equation for using any method for zeros, like Newton's method.

Be happy!
ecam is offline  
March 21st, 2011, 07:49 AM   #3
Newbie
 
Joined: Mar 2011

Posts: 6
Thanks: 0

Re: predictor-corrector scheme

Thank you for answer
I have partial differential equations system so the using of zero methods is more complex, some help I have found in the article where the iterative PC method was used. Sorry for attachment, here it is http://zalil.ru/30711340
Attached Images
File Type: tif FDTD_PC.tif (137.1 KB, 1059 views)
-nad- is offline  
March 21st, 2011, 07:59 AM   #4
Member
 
Joined: Nov 2007

Posts: 73
Thanks: 0

Re: predictor-corrector scheme

Yes, I know it is a PDE system. But when you have the discretized equation, this one on the form u^{new}... is not longer a PDE. That's the point in almost all numerical methods. Think for instance in an implicit Runge-Kutta method for ODE's.
ecam is offline  
March 21st, 2011, 08:27 AM   #5
Newbie
 
Joined: Mar 2011

Posts: 6
Thanks: 0

Re: predictor-corrector scheme

The discretization of ODE and using predictor corrector method is more simple. For me not clear how I must use corrector cycle: for every equation separately or for all system in one cycle. For example the second equation discretization will be does it right to use such scheme where is calculate from other equation on previous step?
-nad- is offline  
March 22nd, 2011, 04:36 AM   #6
Member
 
Joined: Nov 2007

Posts: 73
Thanks: 0

Re: predictor-corrector scheme

No, I think your suggestion is not correct. If you use the U^{predict} instead of you're not going to obtain a correction as n increases.
You have a system for the correction, and this system involves all components of the unknown.
Let's see. Your unknown is . Using the predictor you obtain . Let's take .
The correction procedure is an iteration in the form:
More precisely:





The point is in the right hand side of every equation you have ALL the unknowns. This forces to solve ALL the system at the same time. One method to compute U^{n+1} is to write:






This is a problem about find a zero of a function . You can use the Newton method (or any method for zeros).

Be happy!
ecam is offline  
March 22nd, 2011, 08:44 AM   #7
Newbie
 
Joined: Mar 2011

Posts: 6
Thanks: 0

Re: predictor-corrector scheme

Thank you for developed answer
Please explain what you mean when use . For example for equation (A4c) . What the dependence on two parametrs means?
-nad- is offline  
March 22nd, 2011, 01:52 PM   #8
Member
 
Joined: Nov 2007

Posts: 73
Thanks: 0

Re: predictor-corrector scheme

I am not sure about what you meant.
The F dependence on two parameters just meant (in this case) that F depends on the n and the n+1 approximated solutions.
ecam is offline  
March 22nd, 2011, 08:37 PM   #9
Newbie
 
Joined: Mar 2011

Posts: 6
Thanks: 0

Re: predictor-corrector scheme

I don't understand how the expression for will look exactly (n - time step or number of iteration in time-spatial step ? ). If and
so if we use corrector scheme does it right to write and from where we get on the first corrector iteration for example ?
-nad- is offline  
March 23rd, 2011, 05:07 AM   #10
Member
 
Joined: Nov 2007

Posts: 73
Thanks: 0

Re: predictor-corrector scheme

That's the point. The method it's not explicit. It's an implicit method. That means you need to solve an equation (normal equation, not differential one) to compute every approximation U^{n+1} using U^n. Check this wikipedia article: http://en.wikipedia.org/wiki/Explicit_a ... it_methods
I think it can helps you. Note in the wikipedia example the implicit formula is easy to solve, but in your case probably you need to solve it numerically too (Newton method or any other method for zeros).
ecam is offline  
Reply

  My Math Forum > Science Forums > Computer Science

Tags
predictorcorrector, scheme



Search tags for this page
Click on a term to search for related topics.
Thread Tools
Display Modes


Similar Threads
Thread Thread Starter Forum Replies Last Post
Predictor-corrector RoseNeedsHelp Computer Science 1 December 4th, 2013 12:36 PM
Show the scheme is inconsistent lu5t Applied Math 0 June 10th, 2012 06:18 PM
Stability of implicit spatial backward scheme for 1D wave eq lionthink Applied Math 0 April 7th, 2012 09:56 AM
ODE + predictor and corrector praethorian Pre-Calculus 0 November 2nd, 2008 01:53 PM
Horner's scheme and extrema in matlab hobz Applied Math 0 October 6th, 2007 04:35 AM





Copyright © 2019 My Math Forum. All rights reserved.