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

Method of characterestics (MOC)

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 4, 2013, 11:56
Default Method of characterestics (MOC)
  #1
New Member
 
CFDLearner
Join Date: Jul 2013
Posts: 16
Rep Power: 5
Farouk is on a distinguished road
Hi, I am currently trying to understand the position and state diagrams used in the Method Of Characteristics (used to solve PDE for CFD)...I can't find much useful literature about the topic..so if any one has suggestions, that would be of great help.
Thanks in advance.

Last edited by Farouk; July 7, 2013 at 12:27.
Farouk is offline   Reply With Quote

Old   July 4, 2013, 12:34
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,368
Rep Power: 29
FMDenaro will become famous soon enoughFMDenaro will become famous soon enough
Quote:
Originally Posted by Farouk View Post
Hi, I am curently trying to understand the position and state diagrams used in the Method Of Characterestics (used to solve PDE for CFD)...I can't find much usefull literature about the topic..so if any one has suggestions, that would be of great help.
Thanks in advance.

see the classical book of Zucrow
FMDenaro is offline   Reply With Quote

Old   July 7, 2013, 12:25
Default
  #3
New Member
 
CFDLearner
Join Date: Jul 2013
Posts: 16
Rep Power: 5
Farouk is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
see the classical book of Zucrow
Hi FMDenaro and thanks for the reply.
The text you suggested was very helpful, but still, it didn't answer all my questions, mainly, the use of left going \lambda and right going \betacharacteristics.
For those who are interested, you can see below the two diagrams for a simple wave.
Does anyone happen to have a practical explanation?
Thanks
Attached Images
File Type: jpg Position and State Diagram.jpg (40.5 KB, 20 views)
Farouk is offline   Reply With Quote

Old   July 7, 2013, 12:38
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,368
Rep Power: 29
FMDenaro will become famous soon enoughFMDenaro will become famous soon enough
have a look to the simple wave analysis, homoentropic flows have two Riemann invariants, one along each characteristic curve that allow you to get an exact solution. Generally, the application is reported as a piston moving in a tube without dissipation. The book of Zucrow illustrated such case. Furthermore, the more modern book of LeVeque can be also useful
FMDenaro is offline   Reply With Quote

Old   July 8, 2013, 16:57
Default
  #5
Member
 
Jonny
Join Date: Aug 2009
Posts: 61
Rep Power: 8
Jonny6001 is on a distinguished road
Hello.

Generally speaking the method of characteristics is used to change a PDE to an ODE. It does this by making one of the variables a function of the other variable.
I like to think of the method by considering a second coordinate system that moves along the second variable. It's position along the second variable is a function of the first variable. Then say that the derivative of the function with respect to the new coordinate system is 0. Therefore the characteristic must stay constant along the curves.

This is a very simplified analogy but here goes!

On the 'x-t' diagram with t=0 each characteristic curve has a value. The gradient of the curves emerging from the t=0 represents the speed of information propagation. Think of it as several cars starting from different positions at the same time. If all cars move parallel to each other (have the same velocity) they will never hit each other. However if the car behind moves quicker (has a different gradient to the characteristic curve) it will eventually hit the car in front. This is the point of a shockwave.

The method can be used for the exact solution of equations such as linear advection. A non-linear equation (Burgers) results in curved characteristic curves.

Hope that's not confused you further.
Jonny6001 is offline   Reply With Quote

Old   July 9, 2013, 03:16
Default
  #6
New Member
 
CFDLearner
Join Date: Jul 2013
Posts: 16
Rep Power: 5
Farouk is on a distinguished road
Quote:
Originally Posted by Jonny6001 View Post
Hello.

Generally speaking the method of characteristics is used to change a PDE to an ODE. It does this by making one of the variables a function of the other variable.
I like to think of the method by considering a second coordinate system that moves along the second variable. It's position along the second variable is a function of the first variable. Then say that the derivative of the function with respect to the new coordinate system is 0. Therefore the characteristic must stay constant along the curves.

This is a very simplified analogy but here goes!

On the 'x-t' diagram with t=0 each characteristic curve has a value. The gradient of the curves emerging from the t=0 represents the speed of information propagation. Think of it as several cars starting from different positions at the same time. If all cars move parallel to each other (have the same velocity) they will never hit each other. However if the car behind moves quicker (has a different gradient to the characteristic curve) it will eventually hit the car in front. This is the point of a shockwave.

The method can be used for the exact solution of equations such as linear advection. A non-linear equation (Burgers) results in curved characteristic curves.

Hope that's not confused you further.
Hi Jonny6001 and thanks for the answer. That's exactly the type of answers I was looking for (no mathematics, but logic)...I will try to work some theory and examples, hopefully, I will manage to implement something simple, and i will post back to this thread.

Cheers.
Farouk is offline   Reply With Quote

Old   July 9, 2013, 09:16
Default
  #7
New Member
 
CFDLearner
Join Date: Jul 2013
Posts: 16
Rep Power: 5
Farouk is on a distinguished road
Hi there, does any one happen to have an implemented MOC code? possibly in Matlab? That would be a good start for me.

Regards
Farouk is offline   Reply With Quote

Old   January 15, 2015, 17:23
Default
  #8
New Member
 
yang Shao
Join Date: Jan 2015
Posts: 3
Rep Power: 3
Yashao is on a distinguished road
I'm studying the 3d MOC in reacting flow recently, I read several papers and books and now I understand the characteristic surfaces, and I can get the the compatibility equation, but the following numerical part really makes me headache. How to solve the compatibility equation in finite difference is very hard for me.
I'd appreciate it if anyone can give me some suggestion!
Yashao is offline   Reply With Quote

Reply

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
discretizer - gmshToFoam Andyjoe Open Source Meshers: Gmsh, Netgen, CGNS, ... 13 March 14, 2012 05:35
Fluent 6.3.26 vs 12.1 and partition method Anorky FLUENT 0 April 27, 2010 10:55
Comparison: Finite Volume Method vs. Analytic Method m-fry Main CFD Forum 1 April 20, 2010 14:40
Code for most powerfull FDV Method D.S.Nasan Main CFD Forum 6 September 4, 2008 02:08
Help about the preconditioner in gmres method Dan Gao Main CFD Forum 3 July 5, 2008 01:18


All times are GMT -4. The time now is 03:17.