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

discretization of the 2nd wave equation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 6, 2012, 05:08
Default discretization of the 2nd wave equation
  #1
Senior Member
 
Join Date: Jun 2010
Posts: 101
Rep Power: 7
Hooman is on a distinguished road
Hi,
Can someone pls tell me how the 2nd order wave equation is discretized using the finite volume method? Any references would be appreciated.
Thanks,
Hooman is offline   Reply With Quote

Old   August 6, 2012, 05:17
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,645
Rep Power: 23
FMDenaro will become famous soon enough
Quote:
Originally Posted by Hooman View Post
Hi,
Can someone pls tell me how the 2nd order wave equation is discretized using the finite volume method? Any references would be appreciated.
Thanks,
In order to use FV method, you have to start form a conservative form of the governing equation. That is:

- differential conservative form ... Div (Flux) = ...
- Integral form ........................ Int [S] ( n . Flux) dS = ...

I suggest the book of Leveque on FV for hyperbolic equations
FMDenaro is offline   Reply With Quote

Old   August 6, 2012, 05:23
Default
  #3
Senior Member
 
Join Date: Jun 2010
Posts: 101
Rep Power: 7
Hooman is on a distinguished road
I'll check the book. Tanx.

I know the general rules of FVM just confused with being 2nd order in time and space and would appreciate examples and references.
Hooman is offline   Reply With Quote

Old   August 6, 2012, 05:34
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,645
Rep Power: 23
FMDenaro will become famous soon enough
Quote:
Originally Posted by Hooman View Post
I'll check the book. Tanx.

I know the general rules of FVM just confused with being 2nd order in time and space and would appreciate examples and references.
not all systems of hyperbolic equations can be solved with FVM ... for example, Euler equation in primitive variable [rho, u, p] have not a conservative form, you must go back to the original system in [rho, rho*u, rho*E]
FMDenaro is offline   Reply With Quote

Old   August 6, 2012, 05:42
Default
  #5
Senior Member
 
Join Date: Jun 2010
Posts: 101
Rep Power: 7
Hooman is on a distinguished road
So you're saying that the 2nd order wave equations cannot be solved using FVM?
Hooman is offline   Reply With Quote

Old   August 6, 2012, 05:50
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,645
Rep Power: 23
FMDenaro will become famous soon enough
Quote:
Originally Posted by Hooman View Post
So you're saying that the 2nd order wave equations cannot be solved using FVM?

Assuming the equation is in the form Phi_xx - Phi_yy=0,

you should define a divergence-like operator Div = (d/dx, d/dy) such that

Div F = 0

with F = ( u, v) where u = dPhi/dx, v = - dPhi/dy

and add an equation to close the problem

du/dy + dv/dt = 0
FMDenaro is offline   Reply With Quote

Old   August 6, 2012, 06:00
Default
  #7
Senior Member
 
Join Date: Aug 2011
Posts: 251
Rep Power: 7
leflix is on a distinguished road
Quote:
Originally Posted by Hooman View Post
So you're saying that the 2nd order wave equations cannot be solved using FVM?
the equation is U_tt= cē LAP(U) where U_tt is the second order time derivative and LAP is the laplacian operator.
As Filippo said it for spatial derivative integration is straightforward because LAP(U) = DIV(Grad(U))

For the time derivative discretization just use a forward finite difference formulae for U_tt multiplied by the volume of the cell.
leflix is offline   Reply With Quote

Old   August 6, 2012, 07:39
Default
  #8
Senior Member
 
Join Date: Jun 2010
Posts: 101
Rep Power: 7
Hooman is on a distinguished road
Thank you both. I think I more or less understand. So we need to make the spatial derivative to 1st order?
Hooman is offline   Reply With Quote

Old   August 6, 2012, 10:14
Default
  #9
Senior Member
 
Join Date: Aug 2011
Posts: 251
Rep Power: 7
leflix is on a distinguished road
Quote:
Originally Posted by Hooman View Post
Thank you both. I think I more or less understand. So we need to make the spatial derivative to 1st order?

The way you have to integrate the laplacian using finite volume may be found in any text book. check for Poisson equation using finite volume.
For the second order time derivative,you have to find a forward expression based on finite difference and you will have to multiplied it by the volume of the cell that's all.
leflix is offline   Reply With Quote

Old   August 6, 2012, 10:44
Default
  #10
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,645
Rep Power: 23
FMDenaro will become famous soon enough
Quote:
Originally Posted by leflix View Post
The way you have to integrate the laplacian using finite volume may be found in any text book. check for Poisson equation using finite volume.
For the second order time derivative,you have to find a forward expression based on finite difference and you will have to multiplied it by the volume of the cell that's all.

As leflix said, assuming your equation is Phi_tt = Div Grad (Phi), you integrate over a volume V of boundary S and apply the Gauss theorem to write

Int [V] Phi_tt dV = Int [S] (n . Grad (Phi) ) dS

This equation can be discretized at second order in time and space, for example using central derivatives. How course other issues as numerical stability, numerical oscillations, ecc have to be taken into account
FMDenaro is offline   Reply With Quote

Old   August 6, 2012, 11:04
Default
  #11
Senior Member
 
Join Date: Jun 2010
Posts: 101
Rep Power: 7
Hooman is on a distinguished road
Thank you bother very much.
Hooman 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
Discretization of the continuity equation for with low compressibility Hooman Main CFD Forum 0 April 12, 2012 07:17
Constant velocity of the material Sas CFX 15 July 13, 2010 08:56
Problem with a leapfrog scheme for wave equation Shiranui Main CFD Forum 0 June 22, 2010 09:19
Space and time discretization of Euler equation Hooman Main CFD Forum 2 June 6, 2010 08:30
discretization equation at a corner point (urgent) Sumeet Kumar Main CFD Forum 1 February 16, 2006 10:25


All times are GMT -4. The time now is 06:12.