Exact solution and the multidimensional Godunov scheme for the acoustic equations

The acoustic equations derived as a linearization of the Euler equations are a valuable system for studies of multi-dimensional solutions. Additionally they possess a low Mach number limit analogous to that of the Euler equations. Aiming at understanding the behaviour of the multi-dimensional Godunov scheme in this limit, first the exact solution of the corresponding Cauchy problem in three spatial dimensions is derived. The appearance of logarithmic singularities in the exact solution of the 4-quadrant Riemann Problem in two dimensions is discussed. The solution formulae are then used to obtain the multidimensional Godunov finite volume scheme in two dimensions. It is shown to be superior to the dimensionally split upwind/Roe scheme concerning its domain of stability and ability to resolve multi-dimensional Riemann problems. It is shown experimentally and theoretically that despite taking into account multi-dimensional information it is, however, not able to resolve the low Mach number limit.

Quantity of the inverse problem data for the system of conservation laws

In this paper we study properties of the model, that describes the plane acoustic waves propagation. The model is based on the hyperboliv system of PDE, which is solved numerically by using the finite-volume method, based on Godunov scheme. After studying the direct problem we turn to the inverse one, where our goal is to recover the parameters of the system of PDE by using the initial data, measured in the receivers. We obtain the formula for the gradient of the misfits functional, which allows us to apply gradient-based optimization for recovering the density of the medium. We present the results of numerical experiments for different number of receivers, thus, studying the influence of the quantity of the data of inverse problem on the accuracy of the solution.

High Order Godunov Type Multimesh Method for 3d Impact Problems of Elastoplastic Media

A numerical method for calculating the three-dimensional processes of impact interaction of elastoplastic bodies under large displacements and deformations based on the multi mesh sharp interface method and modified Godunov scheme is presented. To integrate the equations of dynamics of an elastoplastic medium, the principle of splitting in space and in physical processes is used. The solutions of the Riemann problem for first and second order accuracy for compact stencil for an elastic medium in the case of an arbitrary stress state are obtained and presented, which are used at the “predictor” step of the Godunov scheme. A modification of the scheme is described that allows one to obtain solutions in smoothness domains with a second order of accuracy on a compact stencil for moving Eulerian-Lagrangian grids. Modification is performed by converging the areas of influence of the differential and difference problems for the Riemann’s solver. The “corrector” step remains unchanged for both the first and second order accuracy schemes. Three types of difference grids are used. The first – a moving surface grid – consists of a continuous set of triangles that limit and accompany the movement of bodies; the size and number of triangles in the process of deformation and movement of the body can change. The second – a regular fixed Eulerian grid – is limited to a surface grid; separately built for each body; integration of equations takes place on this grid; the number of cells in this grid can change as the body moves. The third grid is a set of local Eulerian-Lagrangian grids attached to each moving triangle of the surface from the side of the bodies and allowing obtain the parameters on the boundary and contact surfaces. The values of the underdetermined parameters in cell’s centers near the contact boundaries on all types of grids are interpolated. Comparison of the obtained solutions with the known solutions by the Eulerian-Lagrangian and Lagrangian methods, as well as with experimental data, shows the efficiency and sufficient accuracy of the presented three-dimensional methodology.

Efficient approach toward the application of the Godunov method to hydraulic transients

The proposed study investigated the applicability of the finite volume method (FVM) based on the Godunov scheme to transient water hammer with shock front simulation, in which intermediate fluxes were computed using either first-order or second-order Riemann solvers. Finite volume (FV) schemes are known to conserve mass and momentum and produce the efficient and accurate realization of shock waves. The second-order solution of the Godunov scheme requires an efficient slope or a flux limiter for error minimization and time optimization. The study examined a range of limiters and found that the MINMOD limiter is the best for modeling water hammer in terms of computational time and accuracy. The first- and second-order FVMs were compared with the method of characteristics (MOCs) and experimental water hammer measurements available in the literature. Both the FV methods accurately predicted the numerical and experimental results. Parallelization of the second-order FVM reduced the computational time similar to that of first-order. Thus, the study presented a faster and more accurate FVM which is comparable to that of MOC in terms of computational time and precision, therefore it is a good substitute for the MOC. The proposed study also investigated the implementation of a more complex convolution-based unsteady friction model in the FVM to capture real pressure dissipation. The comparison with experimental data proved that the first-order FV scheme with the convolution integral method is highly accurate for computing unsteady friction for sudden valve closures.

THE METHOD OF DECOMPOSITION OF GAPES IN THE THREE-DIMENSIONAL DYNAMICS OF ELASTOPLASTIC MEDIA

A numerical method for calculating the three-dimensional processes of impact interaction of elastoplastic bodies with large displacements and deformations based on the method of disintegration of discontinuities according to the Godunov scheme is presented. To integrate the equations of dynamics of an elastoplastic medium, the principle of splitting in space and in physical processes is used. The Riemann's solver for an elastic medium in the case of an arbitrary stress state are obtained and presented. A modification of the scheme is described that allows one to obtain solutions in smoothness domains with a second order of accuracy on a compact template for moving Eulerian – Lagrangian grids. Three types of difference grids are used. The first – a moving surface grid – consists of a continuous set of triangles that limit and accompany the movement of bodies; the size and number of triangles in the process of deformation and movement of the body can vary. The second – a regular fixed Eulerian grid – is limited to a surface grid; separately built for each body; integration of equations takes place on this grid; the number of cells in this grid can change as the body moves. The third grid is a set of local Eulerian – Lagrangian grids attached to each moving triangle of the surface from the side of the bodies and allowing to determine the parameters on the boundary and contact surfaces. The values of the underdetermined parameters near the contact boundaries on all types of grids are interpolated. Comparison of the obtained solutions with the known solutions and with the experimental data, shows the efficiency and sufficient accuracy of the presented three-dimensional methodology.