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

Question about FDM and supersonic flows

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 20, 2013, 14:44
Default Question about FDM and supersonic flows
  #1
New Member
 
Guillaume Fournier
Join Date: Jun 2010
Posts: 3
Rep Power: 7
Guillaume_Fournier is on a distinguished road
Hello everyone,

I have a naive question about Finite Differences Method in the case of supersonic flows.

Correct me if I'm mistaken but when deriving the differential form of the NS equations from the integral form, the use of the divergence theorem implies that the velocity field need to be continuous and differentiable. Therefore, it seems to me that FDM, that come directly from the differential form of the NS equation, are irrelevant to simulate compressible flows with shocks that are, by nature, discontinuous. Yet there are many studies based on FDM for compressible flows.

Am I missing something here ? Is there a numerical artifact I am not aware of that can compensate this mathematical irrelevance ? Is it really mathematically irrelevant ?

Thanks for your help

Cheers,

Guillaume
Guillaume_Fournier is offline   Reply With Quote

Old   January 20, 2013, 18:50
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,659
Rep Power: 23
FMDenaro will become famous soon enough
Quote:
Originally Posted by Guillaume_Fournier View Post
Hello everyone,

I have a naive question about Finite Differences Method in the case of supersonic flows.

Correct me if I'm mistaken but when deriving the differential form of the NS equations from the integral form, the use of the divergence theorem implies that the velocity field need to be continuous and differentiable. Therefore, it seems to me that FDM, that come directly from the differential form of the NS equation, are irrelevant to simulate compressible flows with shocks that are, by nature, discontinuous. Yet there are many studies based on FDM for compressible flows.

Am I missing something here ? Is there a numerical artifact I am not aware of that can compensate this mathematical irrelevance ? Is it really mathematically irrelevant ?

Thanks for your help

Cheers,

Guillaume

Hello,
your doubts are right and I can address some issues, hoping they can be useful. Sorry if I will be lenghty in the response...

As you correctly said, the integral (weak) form of conservation law is the only one for which discontinuous solution can be managed. The differential form has no meaning across a discontinuity.
However this is what one can say for the continuous form of the equations and for the Euler equation for which a singularity can appear also from regular initial condition.

Before discussing the discrete form, let me first use the 1D Burgers equation as example:

Differential (strong) form
du/dt + d/dx(u^2/2)=0 in a periodic domain
u(x,0) = u0(x)= C*sin(x) initial condition

the solution is regular until the time tb = -1/min(u0'(x)), then, you must consider that weak form. The local averaged form can be demonstrated to be a weak form for a particular choice of the test function.
Then, if you integrate the differential Burgers equation over a local volume of measure h centred around x, use the Gauss theorem you get

Integral (weak) form
d u_bar/dt + [u^2(x+h/2,t) - u^2(x-h/2,t)]/2h=0
with u_bar = (1/h) int [x-h/2, x+h/2] u(x',t) dx'.

The issue of having the closing relation u=u[u_bar] is called flux reconstruction from averaged cells. The most simple reconstruction can be to assume u = u_bar, that is accurate to O(h^2) for centred cell.

Now the problem is the "confusion" one can find in several book that use the Finite Difference (FD) or Finite Volume (FV) approximation.
FD discretizations are only for the differential form of the equations and FV for the integral ones. But...

1) Traditionally second order central discretization were used, then if you see the discrete equation:

dU/dt (i,n) + [U^2(i+1,n) - U^2(i-1,n)]/2h =0

it looks like a FD discretization... but if you start from the integral formulation to develop a FV method and use a linear reconstruction for the fluxes you will get the same discrete equation! Thet are a consistent discretization to the integral form of the conservation laws. For discretizations at higher order, FD and FV form can no longer be confused.

2) Sometime the use of FD or FV appears misleading in reading some textbooks. To clarify the point, you can write a differential form for the pointwise field u(x,t) but you can also write a differential form governing the averaged (weak) field u_bar(x,t). Owing to the averaging operator, the field is regularized. This is practically the same thing that leads to write the LES equations: differential equations for a filtered field.

In conclusion, if you read FD in cases where inviscid flows with singularity are treated, you must not think of the discretization of the strong form of the equation but is related to the weak one.

My opinion is that is better to refer to FV formulation coming from the integral form of the equations but I know that the literature is not univocal about this point.

You can find some reading in the book of LeVeque on FV methods as well as in the book of Peric and Ferziger
FMDenaro is online now   Reply With Quote

Old   January 20, 2013, 23:41
Default
  #3
Senior Member
 
Martin Hegedus
Join Date: Feb 2011
Posts: 467
Rep Power: 9
Martin Hegedus is on a distinguished road
Is this what you're trying to say? http://hdl.handle.net/2060/19880020687 From ntrs.nasa.gov.
Martin Hegedus 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
FVM (vs) FDM for compressible flows Rajani Kumar Akula Main CFD Forum 6 April 12, 2002 01:56


All times are GMT -4. The time now is 05:25.