scholarly journals A Cartesian grid method for the three-dimensional numerical simulation of shock wave propagation in complex-shape domains with moving boundaries

Author(s):  
В.В. Елесин ◽  
Д.А. Сидоренко ◽  
П.С. Уткин

Статья посвящена разработке и количественной оценке свойств вычислительного алгоритма метода декартовых сеток для трехмерного математического моделирования распространения ударных волн в областях сложной изменяющейся формы. Представлено подробное описание вычислительного алгоритма, ключевым элементом которого является определение численного потока через грани, по которым внутренние, регулярные ячейки расчетной области соседствуют с внешними, пересекаемыми границами тел ячейками. Работоспособность алгоритма продемонстрирована в результате сравнения рассчитанных и экспериментальных данных в задачах о взаимодействии ударной волны с неподвижной сферой и подвижной частицей. This paper is devoted to the development and quantitative estimation of a numerical algorithm based on the Cartesian grid method for the threedimensional mathematical simulation of shock wave propagation in domains of complex varying shapes. A detailed description of the numerical algorithm is presented. Its key element is the specification of numerical fluxes through the edges that are common for the inner regular cells of the computational domain and the outer cells intersected by the boundaries of the bodies. The efficiency of the algorithm is shown by comparing the numerical and experimental data in the problems of interaction of a shock wave with a fixed sphere and a moving particle.

Author(s):  
Д.А. Сидоренко ◽  
П.С. Уткин

Статья посвящена разработке, программной реализации и количественной оценке свойств вычислительного алгоритма метода декартовых сеток для математического моделирования распространения ударных волн в областях сложной формы с криволинейными границами. Представлено подробное описание вычислительного алгоритма, основанного на методе "h-ячеек''. Работоспособность алгоритма продемонстрирована на задачах о регулярном и простом маховском отражении ударной волны от клина, а также о взаимодействии ударной волны с цилиндром. This paper is devoted to the development, software implementation, and quantitative estimation of a numerical algorithm based on the Cartesian grid method for the mathematical modeling of shock wave propagation in domains of complex shape with curvilinear boundaries. A detailed description of an algorithm based on the method of ``h-boxes'' is given. The efficiency of the algorithm is analyzed on the problems of regular and single Mach reflection of a shock wave from a wedge as well as on the problem of shock wave/cylinder interaction.


Author(s):  
V. V. Elesin ◽  
D. A. Sidorenko ◽  
P. S. Utkin

This paper is devoted to the development and quantitative evaluation of the properties of the numerical algorithm of the Cartesian grid method for three-dimensional (3D) simulation of shock-wave propagation in areas of varying shape. The detailed description of the algorithm is presented. The algorithm is relatively simple to implement and does not require solving the problem of determination of the shape of the body’s boundary intersection with regular computational cell. The accuracy of the algorithm is demonstrated by comparing the simulated and experimental data in the problems of the interaction of a shock wave (SW) with a nonmoving sphere and a moving particle.


2003 ◽  
Vol 67 (3) ◽  
Author(s):  
D. Samsonov ◽  
G. Morfill ◽  
H. Thomas ◽  
T. Hagl ◽  
H. Rothermel ◽  
...  

Author(s):  
Henrik Eklund ◽  
Sven Wedemeyer ◽  
Ben Snow ◽  
David B. Jess ◽  
Shahin Jafarzadeh ◽  
...  

Observations at millimetre wavelengths provide a valuable tool to study the small-scale dynamics in the solar chromosphere. We evaluate the physical conditions of the atmosphere in the presence of a propagating shock wave and link that to the observable signatures in mm-wavelength radiation, providing valuable insights into the underlying physics of mm-wavelength observations. A realistic numerical simulation from the three-dimensional radiative magnetohydrodynamic code Bifrost is used to interpret changes in the atmosphere caused by shock wave propagation. High-cadence (1 s) time series of brightness temperature ( T b ) maps are calculated with the Advanced Radiative Transfer code at the wavelengths 1.309 mm and 1.204 mm, which represents opposite sides of spectral band 6 of the Atacama Large Millimeter/submillimeter Array (ALMA). An example of shock wave propagation is presented. The brightness temperatures show a strong shock wave signature with large variation in formation height between approximately 0.7 and 1.4 Mm. The results demonstrate that millimetre brightness temperatures efficiently track upwardly propagating shock waves in the middle chromosphere. In addition, we show that the gradient of the brightness temperature between wavelengths within ALMA band 6 can potentially be used as a diagnostics tool in understanding the small-scale dynamics at the sampled layers. This article is part of the Theo Murphy meeting issue ‘High-resolution wave dynamics in the lower solar atmosphere’.


2015 ◽  
Vol 2015 ◽  
pp. 1-16 ◽  
Author(s):  
Ru-Chao Shi ◽  
Rui-Yuan Huang ◽  
Guang-Yong Wang ◽  
You-Kai Wang ◽  
Yong-Chi Li

This paper presents numerical simulation of underwater shock wave propagation nearby complex rigid wall. The Ghost Fluid Method (GFM) for the treatment of complex rigid wall is developed. The theoretical analysis on the GFM-based algorithm and relevant numerical tests demonstrate that the GFM-based algorithm is first-order accurate as applied to complex rigid wall. A large number of challenging numerical tests show that the GFM-based algorithm is robust and quite simple in various practical problems. The numerical results on shock wave propagation in the vicinity of rigid wall are verified by comparing to exact solution and the results by body-fitted-grid method.


Sign in / Sign up

Export Citation Format

Share Document