
[Sponsors] 
August 8, 2005, 14:12 
Higher order Hyperbolic systems

#1 
Guest
Posts: n/a

Does anyone know how I can solve (using finite differences) higher order hyperbolic systems (e.g. wave equations in elasticity)? The system I want to solve looks like:
u = u(x,t) u_tt + a(x)*u_t + b(x)*u_x + c(x)*u_xx + d(x)*u_xxx + e(x)*u_xxxx + f(x)*u = g(x,t) _t and _x is partial differentiation with respect to t and x. 

August 8, 2005, 15:13 
Re: Higher order Hyperbolic systems

#2 
Guest
Posts: n/a

The best method is Level Set Method by Stanley Osher and James Sethian
Vinod Dhiman 

August 8, 2005, 16:03 
Re: Higher order Hyperbolic systems

#3 
Guest
Posts: n/a

Why would level set methods be of any help? I have not started thinking about boundary conditions as yet ...


August 8, 2005, 16:22 
Re: Higher order Hyperbolic systems

#4 
Guest
Posts: n/a

Is there a way to solve an equation numerically without specifying boundary conditions. Well, it could be, I really don't know then.
Vinod Dhiman 

August 8, 2005, 16:34 
Re: Higher order Hyperbolic systems

#5 
Guest
Posts: n/a

I would suggest using DRP schemes or compact schemes.They were derived mainly for aeroacoustics but should work fine with the problem you have in hand provided since the schemes are dispersion preserving over a wide range of wave numbers. I would suggest you to look up for papers relating to them using scholar.google.com.
H 

August 9, 2005, 02:35 
Re: Higher order Hyperbolic systems

#6 
Guest
Posts: n/a

Such equations arise in structural analysis (e.g., beam, plate and shell transient response). Actually my M.Sc. thesis was about developing a FEM for shells, where the equations may be degenerated to yours if most of the nonlinearities are omitted. However, I suggest you start from a simpler approach. If you use FEM or FVM, the BC may be simple for treatment if they are of usual types (Dirichlet or Neuman). For FDM, you may need to add ghost cells to take care of this, or otherwise use onesided differences at the boundaries, but ones which are of the same order as your scheme away from the boundaries (e.g., if your scheme is 2nd order, do NOT use the usual 1st order forward / backward schemes for the BC).


August 9, 2005, 09:37 
Re: Higher order Hyperbolic systems

#7 
Guest
Posts: n/a

Rami: Can you tell me which scheme to use? I can't find a way to fit this into the standard hyperbolic conservation form, and so no scheme really comes to mind.


August 9, 2005, 10:14 
Re: Higher order Hyperbolic systems

#8 
Guest
Posts: n/a

Have you checked that your equation is strictly hyperbolic? and more specifically is it wellposed for your choice of boundary conidtions and coefficients (for numerical problems a wellposed problem, in the sense of Cauchy, can appear to be illposed because of the growth of gridscale noise).
As an observation if e(x)>0 ( =1 say) then the dominant balance in your equation is u_tt + u_xxxx = 0 which has solutions of the form exp(+/ i.t.k^2) which oscillate rapidly with increasing t (k is the wavenumber)  the inverse transform will be a Fresnel integral. If e<0 your problem will be illposed since the solution grows like exp(t.k^2). The simplest way to solve this equation is to use explicit 2nd order central differences and march in time with a small timestep. 

August 9, 2005, 10:46 
Re: Higher order Hyperbolic systems

#9 
Guest
Posts: n/a

>>u = u(x,t)
>>u_tt + a(x)*u_t + b(x)*u_x + c(x)*u_xx + d(x)*u_xxx + e(x)*u_xxxx + f(x)*u = g(x,t) Do any of the coefficients of the u_nx terms ever become negative? If so, you will have to be careful of both 'static stability' (xrelated) & 'dynamic stability' (trelated)issues... Smallenough time steps will probably be enough to settle the dynamic stability issues, but static stability issues will put a bound on maximum dx step... diaw... 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Fifth order anti diffusive scheme for 1D hyperbolic heat transfer  areffallah  Main CFD Forum  0  November 3, 2010 05:48 
divergence with higher order scheme  shekharc  Main CFD Forum  1  July 23, 2009 13:53 
Higher order terms in the equations  Argyle  Main CFD Forum  7  October 31, 2008 08:41 
CFL condition for higher order schemes  Shyam  Main CFD Forum  2  February 14, 2008 15:24 
Higher order discretization on staggered grid  Chandra Shekhar  Main CFD Forum  9  January 27, 2005 17:31 