The numerical solution of the problem of stratified fluid flow around a wedge using OpenFOAM

Based on the open source software 2D numerical simulations of incompressible stratified fluids flows have been performed. They are characterized by a wide range of values of internal scales that are not in a homogeneous liquid. Mathematical model is based on the fundamental set of differential equations of inhomogeneous multicomponent fluid mechanics. The problem is solved using the finite volume method in an open source package OpenFOAM. The method allows analyzing in a single formulation the dynamics and fine structure of flow patterns past obstacles in a wide range of flow parameters. A particular attention is focused at construction of a high quality computational grid which satisfies basic requirements for resolution of all the microscales of the problem in high-gradient regions of the flow. The calculations were performed in parallel regime on computational facilities of the web-laboratory UniHUB (www.unihub.ru). The same system of equations and a general numerical algorithm were used for the whole range of the parameters under consideration. The calculation results are in a qualitative agreement with the data from laboratory experiments. Transient flow patterns past obstacles are analyzed, and physical mechanisms are determined, which are responsible for formation of vortices in regions with high density gradients near the edges of an obstacle. For all the velocities of the body motion, the flow field is characterized by a complicated internal structure. In the flow pattern around motionless body dissipative gravity waves are manifested at the edges of the strip. Around the slowly moving body a group of attached waves, are formed in opposite phases at the edges of the wedge. Then, the main flow components become vortices, which are formed around the edge of the wedge and manifested downstream in the wake. With further increase in velocity of the body motion, the flow pattern becomes more non-stationary.

numerical simulation, open source computational packages, stratified flows

1. Chashechkin Yu. D. Structure and dynamics of natural flows: theoretical and laboratory simulation. In Proc.: Aktual'nye problemy mekhaniki. 50 let Institutu problem mekhaniki im. A. Ju. Ishlinskogo RAN [Actual Problems in Mechanics. 50 Years of the A. Yu. Ishlinskiy Institute for Problems in Mechanics of the RAS. Ed. F. L. Chernousko. Moskow: Nauka [Science], 2015, pp. 63-78. ISBN: 978-5-02-039181-9

2. Chashechkin Yu. D., Bardakov R. N., Zagumennyi Ia. V. Calculation and visualization of the fine structure of fields of two-dimensional attached internal waves. Physical Oceanography, 2010, no. 6, pp. 3-15. (in Russian)

3. Brouzet C., Dauxois T., Ermanyuk E., Joubaud S., Kraposhin M., Sibgatullin I. Direct numerical simulation of internal gravity wave attractor in trapezoidal domain with oscillating vertical wall. Trudy ISP RAN/Proc. ISP RAS, 2016, vol. 28, issue 1, pp. 117-141 (in Russian). DOI: 10.15514/ISPRAS-2014-26(5)-6

4. Phillips O. M. On flows induced by diffusion in a stably stratified fluid. Deep-Sea Res., 1970, vol. 17, pp. 435-443.

5. Shapiro A., Fedorovich E. A boundary-layer scaling for turbulent katabatic flow. Boundary-layer meteorology, 2014, vol. 153, issue 1, pp. 1-17. DOI: 10.1007/s10546-014-9933-3.

6. Mercier M. J., Ardekani F. M., Allshouse M. R., Doyle B., Peacock T. Self-propulsion of immersed object via natural convection. Physical Review Letters, 2014, vol. 112, pp. 204501(5). DOI: https://doi.org/10.1103/PhysRevLett.112.204501

7. Dimitrieva N. F., Zagumennyi Ia. V. Diffusion-driven flows on a wedge-shaped obstacle. Physica Scripta, vol. 91, no 8, 2016, pp. 084002. DOI: 10.1088/0031-8949/91/8/084002

8. Dimitrieva N.F., Chashechkin Yu.D. The structure of induced diffusion flows on a wedge with curved edges. Physical Oceanography, 2016, no 3, pp. 77-86. Available at: http://мгфж.рф/index.php/repository?id=120

9. Baydulov V.G., Chashechkin Yu.D. Comparative analysis of symmetries for the models of mechanics of nonuniform fluids. Doklady Physics, 2012, vol. 57, no. 5, pp. 192-196. doi: 10.1134/S1028335812050011

10. Dimitrieva N.F., Chashechkin Yu.D. High-performance numerical simulation of stratified flows around a wedge in OpenFOAM. Trudy ISP RAN/Proc. ISP RAS, 2016, vol. 28, issue 1, pp. 207-220 (in Russian). DOI: 10.15514/ISPRAS-2016-28(1)-12

11. Evstigneev N.M. A finite-volume TVD Riemann solver for the 2D shallow water equations. Vychislitel’nye Metody i Programmirovanie [Numerical Methods and Programming], 2006, vol. 7, no 1, pp. 108-112. Available at: http://num-meth.srcc.msu.ru/zhurnal/tom_2006/pdf/v7r113.pdf

12. Dimitrieva N.F., Chashechkin Yu.D. Numerical simulation of the dynamics and structure of a diffusion-driven flow on a wedge. Vychislitel'naya mekhanika sploshnykh sred [Computational continuum mechanics], 2015, vol. 8, no. 1, pp. 102- 110. (in Russian) doi: 10.7242/1999-6691/2015.8.1.9