
[Sponsors] 
August 6, 2012, 05:08 
discretization of the 2nd wave equation

#1 
Senior Member
Join Date: Jun 2010
Posts: 111
Rep Power: 9 
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, 

August 6, 2012, 05:17 

#2  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,417
Rep Power: 39 
Quote:
 differential conservative form ... Div (Flux) = ...  Integral form ........................ Int [S] ( n . Flux) dS = ... I suggest the book of Leveque on FV for hyperbolic equations 

August 6, 2012, 05:23 

#3 
Senior Member
Join Date: Jun 2010
Posts: 111
Rep Power: 9 
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. 

August 6, 2012, 05:34 

#4 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,417
Rep Power: 39 
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]


August 6, 2012, 05:42 

#5 
Senior Member
Join Date: Jun 2010
Posts: 111
Rep Power: 9 
So you're saying that the 2nd order wave equations cannot be solved using FVM?


August 6, 2012, 05:50 

#6  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,417
Rep Power: 39 
Quote:
Assuming the equation is in the form Phi_xx  Phi_yy=0, you should define a divergencelike 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 

August 6, 2012, 06:00 

#7  
Senior Member
Join Date: Aug 2011
Posts: 271
Rep Power: 9 
Quote:
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. 

August 6, 2012, 07:39 

#8 
Senior Member
Join Date: Jun 2010
Posts: 111
Rep Power: 9 
Thank you both. I think I more or less understand. So we need to make the spatial derivative to 1st order?


August 6, 2012, 10:14 

#9  
Senior Member
Join Date: Aug 2011
Posts: 271
Rep Power: 9 
Quote:
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. 

August 6, 2012, 10:44 

#10  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,417
Rep Power: 39 
Quote:
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 

August 6, 2012, 11:04 

#11 
Senior Member
Join Date: Jun 2010
Posts: 111
Rep Power: 9 
Thank you bother very much.


Thread Tools  
Display Modes  


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 