Computational Analysis of Fluid Forces on an Obstacle in a Channel Driven Cavity: Viscoplastic Material Based Characteristics

In the current work, an investigation has been carried out for the Bingham fluid flow in a channel-driven cavity with a square obstacle installed near the inlet. A square cavity is placed in a channel to accomplish the desired results. The flow has been induced using a fully developed parabolic velocity at the inlet and Neumann condition at the outlet, with zero no-slip conditions given to the other boundaries. Three computational grids, C1, C2, and C3, are created by altering the position of an obstacle of square shape in the channel. Fundamental conservation and rheological law for viscoplastic Bingham fluids are enforced in mathematical modeling. Due to the complexity of the representative equations, an effective computing strategy based on the finite element approach is used. At an extra-fine level, a hybrid computational grid is created; a very refined level is used to obtain results with higher accuracy. The solution has been approximated using P2 − P1 elements based on the shape functions of the second and first-order polynomial polynomials. The parametric variables are ornamented against graphical trends. In addition, velocity, pressure plots, and line graphs have been provided for a better physical understanding of the situation Furthermore, the hydrodynamic benchmark quantities such as pressure drop, drag, and lift coefficients are assessed in a tabular manner around the external surface of the obstacle. The research predicts the effects of Bingham number (Bn) on the drag and lift coefficients on all three grids C1, C2, and C3, showing that the drag has lower values on the obstacle in the C2 grid compared with C1 and C3 for all values of Bn. Plug zone dominates in the channel downstream of the obstacle with augmentation in Bn, limiting the shear zone in the vicinity of the obstacle.

Аtmospheric wind flow around the mountain landscape in the vicinity of Danang airport and flight safety issues

Wind boundary layer flow over the mountain landscape and large structures located around runways (RWs) creates coherent vortex structures (CVSs) that can cross a glideslope and airspace in the vicinity of an airport. The aircraft, encountering a vortex structure, experiences significant changes of the aerodynamic forces and moments, what is especially hazardous due to proximity to terrain. From a mathematical point of view, the solution of this problem presents a challenge due to extremely large space – time scale of the phenomenon, the lack of relevant atmospheric models, as well as comprehensive initial – boundary conditions in numerical modeling. In this paper, a composite solution is constructed: the CVSs area generation is computed in sufficient details within the framework of the grid method. Based on the data obtained in the approximation of analytical functions, an initial vortex structure is formed, the evolution and stochastics of which are modeled within the potential approximation by means of Rankine vortices. The evaluation of the forces and moments increment from the impact of vortex structures on the aircraft was carried out by the panel method using the engineering approach. As an example, the CVSs, resulting from wind flow around the mountainous area of the Son Tra Peninsula, that is located short of RWs 35R-17L and 35L-17R of Da Nang airport, are investigated. To improve the computational grids quality and verify the method of solving the boundary value problem for the Reynolds-averaged Navier-Stokes equations, we used the criteria based on the principle of maximum pressure, requiring Q-parameter positivity property in the vortices cores and flow separation regions. A CVS related aviation event, involving a passenger aircraft MC-21, is studied. The aircraft, after takeoff from RW 35R-17L setting the course close to the direction of the vortex wind structure axis from the Son Tra Peninsula, encountered the mountainous area CVS.

Influence of gravity effect to the recovery rate at uranium in-situ leaching

In-Situ Leaching is a method of extracting minerals by selectively dissolving it with a leaching solution directly in the place of occurrence of the mineral. In practice, during the development of deposits with the In-Situ Leaching method, situations arise when the solution tends to go down below the active thickness of the stratum. This may be due to geological heterogeneity of the rock or gravitational sedimentation of the solution in the rock due to the difference in the densities of the solution and groundwater. As a result of the deposition of the solution along the height, there is a decrease in the recovery of the metal located in the upper part of the geological layers. This article examines the effect of gravity on the flow regime during the filtration of the solution in the rock. The influence of the gravitational effect on the flow of solution in the rock is studied for different ratios of the densities of the solution and groundwater without taking into account the interaction of the solution with the rock. The CUDA technology is used to improve the performance of calculations. The results show that the use of CUDA technology allows to increase the performance of calculations by 40-80 times compared to calculations on a central processing unit (CPU) for different computational grids.

Diffusive limits of 2D well-balanced schemes\\for kinetic models of neutron transport

Two-dimensional dissipative and isotropic kinetic models, like the ones used in neutron transport theory, are considered. Especially, steady-states are expressed for constant opacity and damping, allowing to derive a scattering $S$-matrix and corresponding ``truly 2D well-balanced'' numerical schemes. A first scheme is obtained by directly implementing truncated Fourier-Bessel series, whereas another proceeds by applying an exponential modulation to a former, conservative, one. Consistency with the asymptotic damped parabolic approximation is checked for both algorithms. These findings are confirmed by means of practical benchmarks carried out on coarse Cartesian computational grids.

Non-isothermal vortex flow in the T-junction channel

In this work we present the numerical simulation of coolant mixing modes in the T-junction. We shows that the RANS approach is beneficial for a qualitative flow analysis to obtain relatively agreed averaged velocity and temperature. Moreover, traditionally, the RANS approach calculates only the averaged temperature distribution. It should also be emphasized that unlike the LES approach, the steady RANS approach cannot express a local flow structure in intense mixing zones. Nevertheless, apparently the used RANS approach should be used for assessing the quality of computational grids, boundary conditions in order to use the LES approach for further numerical simulation.

An improved third-order HWCNS for compressible flow simulation on curvilinear grids

Due to the very high requirements on the quality of computational grids, stability property and computational efficiency, the application of high-order schemes to complex flow simulation is greatly constrained. In order to solve these problems, the third-order hybrid cell-edge and cell-node weighted compact nonlinear scheme (HWCNS3) is improved by introducing a new nonlinear weighting mechanism. The new scheme uses only the central stencil to reconstruct the cell boundary value, which makes the convergence of the scheme more stable. The application of the scheme to Euler equations on curvilinear grids is also discussed. Numerical results show that the new HWCNS3 achieves the expected order in smooth regions, captures discontinuities sharply without obvious oscillation, has higher resolution than the original one and preserves freestream and vortex on curvilinear grids.

Numerical simulation of fire in road tunnel. Selection of calculation grid

Рассматривается вопрос выбора расчетной сетки при моделировании пожара в тоннеле с помощью полевого метода и проводится оценка возможного влияния размеров ячеек сетки, а также граничного условия постоянства давления на результаты расчета. Выполнено моделирование пожара для четырех размеров расчетной сетки. Обоснована возможность применения наиболее грубой из используемых сеток с точки зрения инженерных расчетов, в том числе с оговоркой относительно постановки граничного условия. The issue of fire safety of road tunnels is currently an urgent task. Road tunnels are usually not standard typical facilities, but the unique structures. Therefore, it is necessary to study the influence of various parameters on the development of fire in order to take into account the characteristics of a particular object and make decisions on its effective fire protection. Implementation of field tests in this case is expensive and time-consuming. In this regard, numerical modeling is one of the most effective methods of such research. Field models are the most common and currently used for numerical calculations. These models are based on the numerical solution of the system of conservation equations for small control volumes of the calculation grid. This paper examines the issues of selection the calculation grid when modeling a fire in tunnels using the field method is considered and the possible influence of the size of the grid cells is estimated. The mathematical model used in this work is based on a set of differential equations of hydrodynamics, heat transfer, as well as the equation of conservation of the masses of components. Four computational grids were selected for a horizontal (without slope) model tunnel to determine the optimal cell size. As a result of conducted calculations it was established the following: the size of calculated grid is not fundamental for the initial stage of the fire; the use of smaller grid may be preferable at further development of fire, accompanied by increase of combustion capacity to the maximum; the maximum temperature values, especially in the far sections, are obtained on the coarsest grid. The use of such a grid for estimated engineering calculations can be allowed.

Application of Video-to-Video Translation Networks to Computational Fluid Dynamics

In recent years, the evolution of artificial intelligence, especially deep learning, has been remarkable, and its application to various fields has been growing rapidly. In this paper, I report the results of the application of generative adversarial networks (GANs), specifically video-to-video translation networks, to computational fluid dynamics (CFD) simulations. The purpose of this research is to reduce the computational cost of CFD simulations with GANs. The architecture of GANs in this research is a combination of the image-to-image translation networks (the so-called “pix2pix”) and Long Short-Term Memory (LSTM). It is shown that the results of high-cost and high-accuracy simulations (with high-resolution computational grids) can be estimated from those of low-cost and low-accuracy simulations (with low-resolution grids). In particular, the time evolution of density distributions in the cases of a high-resolution grid is reproduced from that in the cases of a low-resolution grid through GANs, and the density inhomogeneity estimated from the image generated by GANs recovers the ground truth with good accuracy. Qualitative and quantitative comparisons of the results of the proposed method with those of several super-resolution algorithms are also presented.