Picture of person typing on laptop with programming code visible on the laptop screen

World class computing and information science research at Strathclyde...

The Strathprints institutional repository is a digital archive of University of Strathclyde's Open Access research outputs. Strathprints provides access to thousands of Open Access research papers by University of Strathclyde researchers, including by researchers from the Department of Computer & Information Sciences involved in mathematically structured programming, similarity and metric search, computer security, software systems, combinatronics and digital health.

The Department also includes the iSchool Research Group, which performs leading research into socio-technical phenomena and topics such as information retrieval and information seeking behaviour.


Stability analysis of second- and fourth-order finite-difference modelling of wave propagation in orthotropic media

Veres, Istvan (2010) Stability analysis of second- and fourth-order finite-difference modelling of wave propagation in orthotropic media. Ultrasonics, 50 (3). pp. 431-438. ISSN 0041-624X

Full text not available in this repository. Request a copy from the Strathclyde author


The stability of the finite-difference approximation of elastic wave propagation in orthotropic homogeneous media in the three-dimensional case is discussed. The model applies second- and fourth-order finite-difference approaches with staggered grid and stress-free boundary conditions in the space domain and second-order finite-difference approach in the time domain. The numerical integration of the wave equation by central differences is conditionally stable and the corresponding stability criterion for the time domain discretisation has been deduced as a function of the material properties and the geometrical discretization. The problem is discussed by applying the method of VonNeumann. Solutions and the calculation of the critical time steps is presented for orthotropic material in both the second- and fourth-order case. The criterion is verified for the special case of isotropy and results in the well-known formula from the literature. In the case of orthotropy the method was verified by long time simulations and by calculating the total energy of the system.