HYDRODYNAMICS AND ACOUSTICS

2021 ◊ Volume 2 (92) ◊ Issue 2 ◊ p. 126-148

A. O. Borysyuk ^*

^* Institute of Hydromechanics of NAS of Ukraine, Kyiv, Ukraine

Studying the flow in a straight flat rigid channel with two axisymmetric rigid-wall narrowings of a rectangular shape. Part 1. Theory

Gidrodin. akust. 2021, 2(2):126-148

https://doi.org/10.15407/jha2021.02.126

TEXT LANGUAGE: Ukrainian

ABSTRACT

The paper deals with developing a numerical method for solving the problem of fluid movement in a straight flat rigid channel with two axisymmetric rigid-wall narrowings of a rectangular shape. The method has the second order of accuracy in terms of spatial coordinates and time. At the same time, the Navier-Stokes and continuity equations are solved in the velocity-pressure variables using integration over elementary volumes into which the computational domain is divided. The resulting integral equations are subject to spatiotemporal discretization, which leads to the need to solve nonlinear algebraic equations. The base of the temporal discretization part is an implicit three-point asymmetric scheme with backward differences, and the base of the spatial one is the TVD scheme with the appropriate discretization scheme of spatial derivatives. The solving of specified algebraic equations includes three stages. First, rewrite the discrete momentum equation as an equation for velocity. Then, based on the discrete continuity equation, derive the pressure equation. After that, apply the procedure of finding and matching successive approximations of these quantities to the resulting coupled nonlinear algebraic equations for velocity and pressure. When calculating the first approximation in the velocity equation, substitute the unknown pressure and velocity (in the expression for the flow) by their values found at the previous time moment. When finding the following approximations, replace these values with already-known previous approximations. The given solution accuracy determines the approximation number. It makes it possible to move from the solution of coupled systems of nonlinear algebraic equations for velocity and pressure to the corresponding independent linear equations. To solve systems of linear algebraic equations, apply an iterative algorithm using methods of delayed correction and conjugate gradients, as well as ICCG (for symmetric matrices) and Bi-CGSTAB (for asymmetric matrices) solvers.

KEY WORDS

flow, channel, narrowing, the Navier-Stokes equations, iterative method

REFERENCES

[1] S. A. Berger and L.-D. Jou, "Flows in stenotic vessels," Annual Review of Fuid Mechanics Research, vol. 32, pp. 347-382, 2000.
https://doi.org/10.1146/annurev.fluid.32.1.347

[2] A. O. Borisyuk, "Experimental study of wall pressure fluctuations in rigid and elastic pipes behind an axisymmetric narrowing," Journal of Fluids and Structures, vol. 26, no. 4, pp. 658-674, 2010.
https://doi.org/10.1016/j.jfluidstructs.2010.03.005

[3] A. O. Borysyuk and Y. A. Borysyuk, "Wall pressure fluctuations behind a pipe narrowing," Science-Based Technologies, vol. 34, no. 2, pp. 162-170, 2017.
https://doi.org/10.18372/2310-5461.34.11615

[4] D. F. Young, "Fluid mechanics of arterial stenoses," Journal of Biomechanical Engineering, vol. 101, no. 3, pp. 157-175, 1979.
https://doi.org/10.1115/1.3426241

[5] A. Aldrovandi, F. Cademartiri, D. Arduini, D. Lina, F. Ugo, E. Maffei, A. Menozzi, C. Martini, A. Palumbo, F. Bontardelli, T. Gherli, L. Ruffini, and D. Ardissino, "Computed tomography coronary angiography in patients with acute myocardial infarction without significant coronary stenosis," Circulation, vol. 126, no. 25, pp. 3000- 3007, 2012.
https://doi.org/10.1161/CIRCULATIONAHA.112.117598

[6] G. Pedrizzetti, "Unsteady tube flow over an expansion," Journal of Fluid Mechanics, vol. 310, no. 1, pp. 89-111, 1996.
https://doi.org/10.1017/S0022112096001747

[7] A. O. Borisyuk, "Study of the flow and acoustic fields in a rigid-walled channel of circular cross-section with a local axisymmetric narrowing. part 1. a theory," International Journal of Fluid Mechanics Research, vol. 34, no. 2, pp. 99-114, 2007.
https://doi.org/10.1615/InterJFluidMechRes.v34.i2.10

[8] A. O. Borisyuk, "Study of the flow and acoustic fields in a rigid-walled channel of circular cross-section with a local axisymmetric narrowing. part 2. numerical results," International Journal of Fluid Mechanics Research, vol. 34, no. 3, pp. 191-209, 2007.
https://doi.org/10.1615/InterJFluidMechRes.v34.i3.10

[9] S. G. Mirolyubov, "Hydrodynamics of a stenosis," Modern Problems of Biomechanics, vol. 1, pp. 73-136, 1983.

[10] I. V. Vovk, V. T. Grinchenko, and V. S. Malyuga, "Special features of fluid motion in a duct with stenoses," Applied Hydromecanics, vol. 11 (83), no. 4, pp. 17-30, 2009.

[11] V. S. Malyuga, "Numerical investigation of flow in a duct with two serial stenosis. an algorithm of the solution," Applied Hydromecanics, vol. 12 (84), no. 4, pp. 45-62, 2010.

[12] A. V. Shaldenko and A. A. Gourzhii, "Analysis of the heat-transfer processes in a straight channel with insertions at low Reynolds numbers," Applied Hydromecanics, vol. 17 (89), no. 3, pp. 55-66, 2017.

[13] E. V. Brujatskiy, A. G. Kostin, and E. I. Nikiforovich, "Numerical investigation of the velocity and pressure fields in a flat duct having a rectangular ledge on the wall," Applied Hydromecanics, vol. 13 (85), no. 3, pp. 33-47, 2011.

[14] C. Hirsch, Numerical computation of internal and external flows. Butterworth-Heinemann, 2007.

[15] J. H. Ferziger and M. Peri'c, Computational methods for fluid dynamics. Springer, 2002.
https://doi.org/10.1007/978-3-642-56026-2

[16] B. van Leer, "Towards the ultimate conservative difference scheme. II. monotonicity and conservation combined in a second-order scheme," Journal of Computational Physics, vol. 14, no. 4, pp. 361-370, 1974.
https://doi.org/10.1016/0021-9991(74)90019-9

[17] S. Chakravarthy and S. Osher, "High resolution applications of the Osher upwind scheme for the Euler equations," in 6th Computational Fluid Dynamics Conference. American Institute of Aeronautics and Astronautics, 1983.
https://doi.org/10.2514/6.1983-1943

[18] P. K. Sweby, "High resolution schemes using flux limiters for hyperbolic conservation laws," SIAM Journal on Numerical Analysis, vol. 21, no. 5, pp. 995-1011, 1984.
https://doi.org/10.1137/0721062

[19] N. P. Waterson and H. Deconinck, "Design principles for bounded higher-order convection schemes - a unified approach," Journal of Computational Physics, vol. 224, no. 1, pp. 182-207, 2007.
https://doi.org/10.1016/j.jcp.2007.01.021

[20] R. I. Issa, "Solution of the implicitly discretised fluid flow equations by operator-splitting," Journal of Computational Physics, vol. 62, no. 1, pp. 40-65, 1986.
https://doi.org/10.1016/0021-9991(86)90099-9

[21] H. K. Versteeg and W. Malalasekera, An introduction to computational fluid dynamics. Pearson Education Ltd, 2007.

[22] J. D. j. Anderson, Computational fluid dynamics. The basics with applications. McGraw-Hill, 1995.

[23] P. K. Khosla and S. G. Rubin, "A diagonally dominant second-order accurate implicit scheme," Computers & Fluids, vol. 2, no. 2, pp. 207-209, 1974.
https://doi.org/10.1016/0045-7930(74)90014-0

[24] M. R. Hestenes and E. L. Stiefel, "Method of conjugate gradients for solving linear systems," Journal of Research of the National Bureau of Standards, vol. 49, no. 6, pp. 409-436, 1952.
https://doi.org/10.6028/jres.049.044

[25] D. A. H. Jacobs, Preconditioned conjugate gradient methods for solving systems of algebraic equations, 1980, p. 31.

[26] R. Barrett, M. Berry, T. F. Chan, J. Demmel, J. Donato, J. Dongarra, V. Eijkhout, R. Pozo, C. Romine, and H. Van der Vorst, Templates for the solution of linear systems: Building blocks for iterative methods. Philadelphia, PA: SIAM, 1994.
https://doi.org/10.1137/1.9781611971538

[27] H. A. van der Vorst, "Bi-CGSTAB: A fast and smoothly converging variant of bi-CG for the solution of nonsymmetric linear systems," SIAM Journal on Scientific and Statistical Computing, vol. 13, no. 2, pp. 631-644, 1992.
https://doi.org/10.1137/0913035

[28] Netlib repository at utk and ornl.

[29] E. V. Brujatskiy and A. G. Kostin, "Numerical investigation of fluid motion in a closed rectangular cavity with movable upper lid," Applied Hydromecanics, vol. 11 (83), no. 1, pp. 3-15, 2009.

[30] A. A. Prikhodko and D. A. Redchidz, "Numerical simulation of unsteady flow in a wake behind a cylinder on the basis of the Navier-Stokes equations," Applied Hydromecanics, vol. 7 (79), no. 1, pp. 56-71, 2005.