Higher-Order Numerical Methods for Transient Wave Equations

2003 ◽  
Vol 114 (1) ◽  
pp. 21-21 ◽  
Author(s):  
Gary C. Cohen ◽  
Qing Huo Liu
Author(s):  
Jairo Rodríguez-Padilla ◽  
Daniel Olmos-Liceaga

The implementation of numerical methods to solve and study equations for cardiac wave propagation in realistic geometries is very costly, in terms of computational resources. The aim of this work is to show the improvement that can be obtained with Chebyshev polynomials-based methods over the classical finite difference schemes to obtain numerical solutions of cardiac models. To this end, we present a Chebyshev multidomain (CMD) Pseudospectral method to solve a simple two variable cardiac models on three-dimensional anisotropic media and we show the usefulness of the method over the traditional finite differences scheme widely used in the literature.


Author(s):  
Manzoor Hussain ◽  
Sirajul Haq

In this paper, meshless spectral interpolation technique using implicit time stepping scheme is proposed for the numerical simulations of time-fractional higher-order diffusion wave equations (TFHODWEs) of variable coefficients. Meshless shape functions, obtained from radial basis functions (RBFs) and point interpolation method (PIM), are used for spatial approximation. Central differences coupled with quadrature rule of [Formula: see text] are employed for fractional temporal approximation. For advancement of solution, an implicit time stepping scheme is then invoked. Simulations performed for different benchmark test problems feature good agreement with exact solutions. Stability analysis of the proposed method is theoretically discussed and computationally validated to support the analysis. Accuracy and efficiency of the proposed method are assessed via [Formula: see text], [Formula: see text] and [Formula: see text] error norms as well as number of nodes [Formula: see text] and time step-size [Formula: see text].


Author(s):  
Yanbin He ◽  
Tianning Chen ◽  
Jinghuai Gao

Abstract The perfectly matched layer (PML) has been demonstrated to be an efficient absorbing boundary for near-field wave simulation. For heterogeneous media, the property of the PML needs to be carefully specified to avoid numerical instability and artificial reflection because part of it lies at the discontinuous interface. Coupled acoustic-poroelastic (A-P) media or coupled elastic-poroelastic (E-P) media often arise in the field of geophysics. However, PMLs that appropriately terminate these heterogeneous poroelastic media are still lacking. The main purpose of this paper is to explore the application of unsplit PMLs for transient wave modeling in infinite, heterogeneous, coupled A-P media or coupled E-P media. To this end, a consistent derivation of memory-efficient PML formulations for the second-order Biot's equations, elastic wave equations and acoustic wave equations is performed based on complex coordinate transformation using auxiliary differential equations. Furthermore, the interface boundary conditions inside the absorbing layer are rigorously derived for the considered A-P and E-P cases. Finally, the weak form of PML formulations for coupled poroelastic problems is presented. The finite element method is used to validate the proposed PML based on several two-dimensional benchmarks. The accuracy and stability of weak PML formulations are investigated. In particular, for coupled acoustic-poroelastic PML, two extreme (open-pore and sealed-pore) interface conditions are considered and PML results are compared with known analytical solutions. This study demonstrates the ability of the PML to effectively eliminate outgoing bulk waves and surface waves in coupled poroelastic media.


2012 ◽  
Vol 14 (05) ◽  
pp. 1250032 ◽  
Author(s):  
BRUNO FRANCHI ◽  
MARIA CARLA TESI

In this paper we define Maxwell's equations in the setting of the intrinsic complex of differential forms in Carnot groups introduced by M. Rumin. It turns out that these equations are higher-order equations in the horizontal derivatives. In addition, when looking for a vector potential, we have to deal with a new class of higher-order evolution equations that replace usual wave equations of the Euclidean setting and that are no more hyperbolic. We prove equivalence of these equations with the "geometric equations" defined in the intrinsic complex, as well as existence and properties of solutions.


2010 ◽  
Vol 80 (8) ◽  
pp. 1734-1745
Author(s):  
K. Kozel ◽  
P. Louda ◽  
J. Příhoda

Sign in / Sign up

Export Citation Format

Share Document