scholarly journals A Fourth Order Finite Difference Method for the Good Boussinesq Equation

2014 ◽  
Vol 2014 ◽  
pp. 1-10 ◽  
Author(s):  
M. S. Ismail ◽  
Farida Mosally

The “good” Boussinesq equation is transformed into a first order differential system. A fourth order finite difference scheme is derived for this system. The resulting scheme is analyzed for accuracy and stability. Newton’s method and linearization techniques are used to solve the resulting nonlinear system. The exact solution and the conserved quantity are used to assess the accuracy and the efficiency of the derived method. Head-on and overtaking interactions of two solitons are also considered. The numerical results reveal the good performance of the derived method.

1989 ◽  
Vol 79 (4) ◽  
pp. 1210-1230
Author(s):  
C. R. Daudt ◽  
L. W. Braile ◽  
R. L. Nowack ◽  
C. S. Chiang

Abstract The Fourier method, the second-order finite-difference method, and a fourth-order implicit finite-difference method have been tested using analytical phase and group velocity calculations, homogeneous velocity model calculations for disperson analysis, two-dimensional layered-interface calculations, comparisons with the Cagniard-de Hoop method, and calculations for a laterally heterogeneous model. Group velocity rather than phase velocity dispersion calculations are shown to be a more useful aid in predicting the frequency-dependent travel-time errors resulting from grid dispersion, and in establishing criteria for estimating equivalent accuracy between discrete grid methods. Comparison of the Fourier method with the Cagniard-de Hoop method showed that the Fourier method produced accurate seismic traces for a planar interface model even when a relatively coarse grid calculation was used. Computations using an IBM 3083 showed that Fourier method calculations using fourth-order time derivatives can be performed using as little as one-fourth the CPU time of an equivalent second-order finite-difference calculation. The Fourier method required a factor of 20 less computer storage than the equivalent second-order finite-difference calculation. The fourth-order finite-difference method required two-thirds the CPU time and a factor of 4 less computer storage than the second-order calculation. For comparison purposes, equivalent runs were determined by allowing a group velocity error tolerance of 2.5 per cent numerical dispersion for the maximum seismic frequency in each calculation. The Fourier method was also applied to a laterally heterogeneous model consisting of random velocity variations in the lower half-space. Seismograms for the random velocity model resulted in anticipated variations in amplitude with distance, particularly for refracted phases.


2014 ◽  
Vol 580-583 ◽  
pp. 2924-2927 ◽  
Author(s):  
Pavel A. Akimov ◽  
Mojtaba Aslami

This paper is devoted to correct and efficient method of local static analysis of Bernoulli beam on elastic foundation. First of all, problem discretized by finite difference method, and then transformed to a localized one by using the Haar wavelets. Finally, imposing an optimal reduction in wavelet coefficients, the localized, reduced results can be obtained. It becomes clear after comparison with analytical solutions, that the localization of the problem by multiresolution wavelet approach gives exact solution in desired regions of beam even in high level of reduction in wavelet coefficients. This localization can be applied to any arbitrary region of the beam by choosing optimum reduction matrix and obtaining exact solutions with an acceptable reduced size of the problem.


Sign in / Sign up

Export Citation Format

Share Document