On shock localization by digital image processing techniques (Q1088549)

A general approach is presented to a problem of localization of strong discontinuities in the spatial domain on the basis of numerical solution of multidimensional steady and unsteady problems by finite difference shock-capturing techniques. The essence of the aproach is in the use of ideas and methods of the digital image processing theory, in particular, of methods for the detection of edges in the image brightness. The algorithms of localization of strong discontinuities in the difference solutions developed on the above-mentioned basis do not require for their realization any a priori information on the orientation of shock surfaces with respect to the spatial co-ordinate axes. The questions of accuracy of the new localization algorithms and of increasing the accuracy of the numerical solutions with the aid of image restoration techniques are considered. The universality of the proposed algorithms is demonstrated in the examples of computations of a number of model two-dimensional unsteady problems by the Godunov scheme and by the MacCormack scheme.