Flow Of Viscous Incompressible Fluid Research Articles

STRUCTURED PSEUDOSPECTRA IN PROBLEMS OF SPATIAL STABILITY OF BOUNDARY LAYERS

The paper is devoted to the numerical analysis of the sensitivity of the characteristics of spatial stability of boundary layers to the errors with which the main flow is specified. It is proposed to use structured pseudospectra for this purpose. It is shown that the obtained estimates are significantly more accurate than the estimates based on the unstructured pseudospectrum. The presentation is carried out using the example of a viscous incompressible fluid flow over a concave surface of small curvature with flow parameters favorable for the development of Goertler vortices and Tollmien–Schlichting waves.

Метод R-функцій для розв’язання задачі масообміну циліндричного тіла з твірною у вигляді кривої Ламе з в’язкою нестисливою рідиною

The problem of stationary mass exchange of a cylindrical body with viscous incompressible fluid flow is reduced to solving the equation for concentration with the corresponding boundary conditions on the surface of the body and at infinity. Applying the constructive apparatus of the R-functions theory allows us to accurately take into account the geometry of the domain, as well as the boundary conditions together with the condition at infinity. The R-function method was proposed by V. L. Rvachev, Ukraine National Academy of Sciences academician. For boundary value problems of mathematical physics, this method allows constructing the so-called structure of the solution of the boundary value problem, i.e. a bundle of functions that accurately takes into account all boundary conditions of the problem and depends on some uncertain components. The choice of these uncertain components is performed in such a way as to satisfy (in a certain sense) the equation of the problem. For this, one can use standard numerical methods of mathematical physics (the Ritz method, the Galerkin method, the least squares method, collocations method, etc.). It should be noted that the geometry of the region is taken into account exactly, i.e. in particular, there is no replacement of curvilinear sections of the boundary with polygonal lines inscribed in them, as it is done, for example, in the finite element method. The purpose of this article is to apply the R-functions method to the problem of mass exchange of a cylindrical body formed by the Lame curve, with a viscous incompressible fluid flowing around it. To achieve this goal, a complete structure of the solution of a linear boundary value problem for concentration is constructed by the R-functions theory methods, and a numerical algorithm based on the Galerkin method for approximating indefinite components in the problem structure is used. The degree of rigor and the conditions for using the considered equations of fluid motion are not considered in the work; these equations are considered as mathematical models subject to numerical algorithmization

Approximation of Classical Two-Phase Flows of Viscous Incompressible Fluids by a Navier–Stokes/Allen–Cahn System

We show convergence of the Navier–Stokes/Allen–Cahn system to a classical sharp interface model for the two-phase flow of two viscous incompressible fluids with same viscosities in a smooth bounded domain in two and three space dimensions as long as a smooth solution of the limit system exists. Moreover, we obtain error estimates with the aid of a relative entropy method. Our results hold provided that the mobility mε>0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_\\varepsilon >0$$\\end{document} in the Allen–Cahn equation tends to zero in a subcritical way, i.e., mε=m0εβ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_\\varepsilon = m_0 \\varepsilon ^\\beta $$\\end{document} for some β∈(0,2)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\beta \\in (0,2)$$\\end{document} and m0>0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_0>0$$\\end{document}. The proof proceeds by showing via a relative entropy argument that the solution to the Navier–Stokes/Allen–Cahn system remains close to the solution of a perturbed version of the two-phase flow problem, augmented by an extra mean curvature flow term mεHΓt\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$m_\\varepsilon H_{\\Gamma _t}$$\\end{document} in the interface motion. In a second step, it is easy to see that the solution to the perturbed problem is close to the original two-phase flow.

Perturbation method for laminar flow of viscous incompressible fluid between two porous boundary walls

The present study focuses on the investigation of the velocity profile of viscous fluid between two parallel porous walls while fluid is being ejected from both channel walls at the same rate. The distributions of velocity, pressure and shear stress were determined theoretically for incompressible, steady, fully developed laminar flow past a porous boundary wall. Due to the pressure difference across the porous wall, there is ejection or injection into the wall. The velocity specified as the boundary condition at the wall has been used to explore fluid flow phenomena in porous tubes and ducts. The variable wall velocity is based on the pressure differential across the wall, the fluid characteristics, the thickness of the wall, and the permeability of the structure. An appropriate wall condition must take these factors into account. This paper discusses a Reynolds number solution. We seek formulas for the velocity component. The perturbation approach is used to solve the governing differential equation and the general solution of the velocity is obtained for higher terms of Re. Fifth term of perturbation is used in the perturbation method. The velocity of the fluid is analyzed for different values of Re for a boundary condition. Finally, the result is analysis graphically for the range of Re from .1 to 50. In addition, another strategy is used to keep the Reynolds number constant when the diameter of the hole in the porous wall changes, i.e., changing the value of λ.

Some Exact Solutions of Unsteady Viscous Incompressible Fluid Flows with Topology Conservation of Vorticity Vector Field by Ansatz Method

Some Exact Solutions of Unsteady Viscous Incompressible Fluid Flows with Topology Conservation of Vorticity Vector Field by Ansatz Method

On well-posedness of Navier–Stokes variational inequalities

In this paper, well-posedness is carried out on variational inequalities in a general form for the stationary Navier–Stokes equations modeling viscous incompressible fluid flows with a nonsmooth slip boundary condition of friction type. For a special form of the variational inequality commonly studied in existing literature, our proof technique leads to a substantial improvement on the well-posedness condition.

Numerical study of transitions in lid-driven flows in shallow cavities

In this article, three dimensional (3D) lid-driven flow in shallow cavities with a unit square base are studied. The numerical solution of the Navier–Stokes equations modeling incompressible viscous fluid flow in a cavity is obtained via a methodology combining a first order accurate operator-splitting scheme, a L 2 -projection Stokes solver, a wave-like equation treatment of the advection and finite element space approximations. Numerical results of a lid-driven flow in a cubic cavity show a good agreement with those reported in literature. The critical Reynolds numbers (Re cr ) for having flow with increasing of oscillating amplitude (a Hopf bifurcation) in different shallow cavities are obtained and associated oscillating modes are studied.

Exact solutions to vorticity of the fractional nonuniform Poiseuille flows

Abstract Closed-form expressions for the dimensionless velocity, shear stresses, and the flow vorticity fields corresponding to the isothermal unsteady Poiseuille flows of a fractional incompressible viscous fluid over an infinite flat plate are established. The fluid motion induced by a pressure gradient in the flow direction is also influenced by the flat plate that oscillates in its plane. The vorticity field is dependent on two spatial coordinate and time, and it is an arbitrary trigonometric polynomial in the horizontal coordinate. The exact solutions, obtained by generalized separation of variables and Laplace transform technique, are presented in terms of the Wright function and complementary error function of Gauss. Their advantage consists in the fact that the values of the fractional parameter can be chosen so that the predicted material properties by them to be in agreement with the corresponding experimental results. In addition, they describe motions for which the nontrivial shear stresses are influenced by history of the shear rates. It is found that the flow vorticity is stronger near the plate, but it could be attenuated in the case of fractional model.

Optimization of mass and heat flux of MHD viscous fluid flow with constant proportional Caputo derivative by using response surface methodology: Sensitivity analysis

The study of viscous dissipation in the context of unsteady incompressible magnetohydrodynamic viscous fluid flow over a moving plate is a complex and important topic in fluid dynamics and heat transfer. This scenario combines several physical phenomena and has practical applications in various engineering and scientific fields. The fractional mathematical model is proposed by considering the effect of viscous dissipation force of the flowing fluid. For this purpose, fluid is assumed to be lying over a vertical plate that is moving in its own plane. The governing equations are transformed into the dimensionless form and developed fractional model with a new operator Constant Proportional Caputo Derivative. We used the Laplace transform technique to find the solution of the dimensionless governing equation analytically. The transformed solutions for energy and momentum balances developed in terms of series. The validation of the present model comparison graphs is plotted for temperature. The sensitivity of different input parameters on output response are analyzed and displayed the sensitivity results in tabular and graphical form and found that the permeable parameter is most sensitive to skin friction coefficient, and the Eckert number is most sensitive parameter among others to Nusselt Number. The novelty of present work is to develop a correlation between input parameters and output responses by using Response Surface Methodology (RSM). As no such correlation has been developed for optimization of mass and heat flux of MHD viscous fluid flow with Constant Proportional Caputo derivative by using RSM. By using this correlation, we have performed sensitivity analysis for output responses.

A Hermite-type collocation mesh-free approach for simulating incompressible viscous fluid flows

In this work, a novel Hermite-type mesh-free model for incompressible viscous fluid flows is proposed, as a way, to overcome the numerical difficulties experienced by solving under weak or strong formulation the pure stream-function Navier–Stokes equation. In this model a dimensionless virtual domain is utilized, in order to make an accurate dimensionless shape-functions of the recently developed Hermite-type collocation Weighted Least Square approximation, by using dimensionless scale parameters. Then, a Jacobian transformation matrix is built to express the shape-functions spatial derivatives from the virtual to the considered domain. The resolution of the obtained nonlinear algebraic system is performed by the High Order Continuation Method. Numerical tests were investigated on incompressible viscous isothermal fluid flows in lid-driven cavity, backward-facing step, sudden expansions and around various obstacle shapes. They show the advantages of the proposed model compared to the reference Hermite-type collocation models used for pure stream-function equation and the weak Galerkin procedures as Finite Element Method and Element Free Galerkin method, in term of computation accuracy, number of degree of freedom, implementation simplicity and computation time. The numerical and experimental results of literature confirm the validity and capability of the proposed dimensionless model to simulate fluid flows in a expended-lengths inflow–outflow geometry and around obstacles of circular, inclined elliptical and Joukowski airfoil shapes.

Effect of the Non-Linear Radiative Unsteady Mixed Convective Flow over a Curved Stretching Surface with Soret and Dufour Effects: A Numerical Study

The subject of unsteady convective flow with non-linear thermal radiation has become an important issue of research, due to its implications in advanced energy conversion systems operating at high temperature, solar energy technology and chemical process at high operation temperature. Due to the importance of this issue, a time dependent incompressible viscous fluid flow, heat and mass transfer over a curved stretching surface has been numerically analysed by taking into account the heat flux due to concentration gradient and mass flux due to temperature gradient. Together with this the Rosseland approximation is being employed for the nonlinear thermal radiation impact in presence of thermal slip. With the aid of non-dimensional variables and the corresponding physical boundary conditions, the leading nonlinear momentum, energy, and species equations are converted into a set of coupled ordinary differential equations. These equations are then resolved using the MATLAB bvp4c solver. The stability of the numerical technique has been verified and compared with available literatures. The resultant parameters of engineering interest and the boundary layer flow field parameters and have been presented using tables and graphically plots. The study concludes that for lesser curvature parameter (0.5≤K≤0.7) the surface drag force, heat and mass transfer rates can improve by about 9.59%, 2.87% and 1.67% each respectively. The presence of the temperature ratio parameter and the non-linear thermal radiation are found to greatly influence the temperature profile and the heat transfer rate of the system. Results show that the heat transfer rate improves by about 24.39% and 16.66% for varying non-linear thermal radiation (1≤Rd≤1.5) and temperature ratio parameter (1.2≤θw≤1.4) respectively. Results obtained also show that improving the thermal slip parameter (0.4≤L≤0.6) can reduce heat transfer rate by about 13.62% and reduce the surface temperature profile.

Asymptotic behavior of a viscous incompressible fluid flow in a fractal network of branching tubes

<p>We considered a viscous incompressible fluid flow in a varying bounded domain consisting of branching thin cylindrical tubes whose axes are line segments that form a network of pre-fractal curves constituting an approximation of the Sierpinski gasket. We supposed that the fluid flow is driven by volumic forces and governed by Stokes equations with boundary conditions for the velocity and the pressure on the wall of the tubes and inner continuity conditions for the normal velocity on the interfaces between the junction zones and the rest of the pipes. We constructed local perturbations, related to boundary layers in the junction zones, from solutions of Leray problems in semi-infinite cylinders representing the rescaled junctions. Using $ \Gamma $-convergence methods, we studied the asymptotic behavior of the fluid as the radius of the tubes tends to zero and the sequence of the pre-fractal curves converges in the Hausdorff metric to the Sierpinski gasket. Based on the constructed local perturbations, we derived, according to a critical parameter related to a typical Reynolds number of the flow in the junction zones, three effective flow models in the Sierpinski gasket, consisting of a singular Brinkman flow, a singular Darcy flow, and a flow with constant velocity.</p>

Dissipation in nonequilibrium thermodynamics and its connection to the Rayleighian functional

We examine quantitatively the role of dissipation in nonequilibrium thermodynamics and its connection to variational principles and the Rayleighian functional. The extremum of the Rayleighian is sometimes used to describe the inertialess (dissipation-dominated) dynamics of continuum systems, and it has been applied recently for the modeling of soft matter dynamics. We discuss how dissipation is considered within one of the modern complete descriptions of nonequilibrium thermodynamics, namely the single generator bracket formalism. Within this formalism, dissipation is introduced through the use of the dissipation bracket, describing irreversible dynamics, which is added to a Poisson bracket that describes the reversible dynamics of the system. A possible connection with the Rayleighian functional is then demonstrated that in all cases considered herein, the Rayleighian is equal to minus one half of the effective dissipation rate of the Lagrangian functional. The effective dissipation rate is obtained starting with an inertial (i.e., flux-based or velocity-based) system description, involving the Poisson bracket and the primitive part (i.e., without the entropy correction term) of the dissipative bracket. Several examples are discussed in detail, ranging from an algebraic model (damped oscillator) to continuum ones: modeling of fluid flow in porous particle media, viscous Newtonian compressible and incompressible fluid flows, and more interestingly, flow of a nematic liquid-crystalline material.

Analytical Approach to Nonlinear Oscillatory Flow Through a Semi-Infinite Vertical Permeable Wall in a Porous Medium with MHD Mixed Convection

The present work illustrates mixed MHD convection flow of an unstable two-dimensional viscous incompressible fluid flow in a porous material under the following conditions: semi-infinite vertical permeable wall support (i) There is a constant mean oscillation in the suction velocity normal to the wall throughout time. (ii) A constant mean velocity of the free stream; (iii) a constant wall temperature; and (iv) a somewhat high temperature differential between the wall and the free stream, which generates free convection currents. Conservation of mass, energy and momentum are used to obtain the governing equations. These equations are non-linear in nature and connected with one another. The nonlinear differential equations are non-dimensionalized employing non-dimensionalized parameters such as Prandtl number Pr, Grashof number Gr, Eckert number Ek, Hartmann number Mh, and Viscous ratio. An analytical method is used to solve approximate solutions for nonlinear PDEs utilizing the double regular perturbation technique. The numerical results for temperature, skin friction, heat flow, velocity, and skin friction are computed for various parameters and are demonstrated to be in good agreement with previous findings.

Зарождение вихреобразования на границе кругового цилиндра при его обтекании потоком вязкой жидкости

The paper considers the problem of the generation of vortex formation on the surface of a circular cylinder placed in a flow of incompressible viscous fluid. To study this process, a new time-iterative approach is formu-lated for solving the 2D problem of flow around the circle based on the BIE method with Green's functions for the outer part of the cylinder with the Dirichlet boundary condition on its surface. A new approach, based on the transient Navier-Stokes equation in the variables stream function - vortex, using classical methods of potential theory and the boundary element method, allows us to obtain a complete system of three integral equations that determine the fields of stream functions and vorticity at a new time step, if these fields are known from the previous step. The assumption about the initial unexcited state of the flow (the vorticity function in the entire flow region is zero) allows us to obtain an analytical solution at the first time iteration step. The analysis of the obtained solution depending on the time step allows us to formulate a hypothesis about the nature of the generation of vortex formation on the surface of the circular cylinder in a flow of incompressible viscous fluid. The results obtained describe the appearance on the surface of the cylinder of two generalized vortices, having infinitely high intensity and infinitely small localization near the surface of the cylinder, which “instantly” break away from the surface of the cylinder, with a significant decrease in intensity, but with a large coverage of the flow around the cylinder.

Analysis of incompressible viscous fluid flow in convergent and divergent channels with a hybrid meta-heuristic optimization techniques in ANN: An intelligent approach

Analysis of incompressible viscous fluid flow in convergent and divergent channels with a hybrid meta-heuristic optimization techniques in ANN: An intelligent approach

Exact solutions for the description of nonuniform unidirectional flows of magnetic fluids in the Lin–Sidorov–Aristov class

The paper considers the exact integration of magnetic hydrodynamic equations for describing nonuniform unidirectional flows of viscous incompressible fluids. The construction of an exact solution is based on the well-known representation of hydrodynamic fields as the Lin–Sidorov–Aristov class. The 3d magnetic field is described by linear forms with respect to two spatial coordinates (longitudinal, or horizontal). The coefficients of the linear forms depend on the third coordinate and time. In view of the incompressibility condition, the 1D velocity field depends on two coordinates and time. The pressure is shown to be determined by a quadratic form with constant coefficients. These coefficients are determined by pressure distribution on the known (free) boundary. The exact solution is illustrated by the integration of non-1D hydrodynamic fields in the case of the steady motion of a conducting viscous incompressible fluid. This solution is polynomial, and it will be useful for the formulation of new problems of hydrodynamic stability.

Weakly viscous two-dimensional incompressible fluid flows using efficient isogeometric finite element method

This paper presents a straightforward and efficient numerical simulation method for solving the Navier–Stokes equations for weakly viscous incompressible fluids describing steady flow. Our approach utilizes isogeometric finite elements to handle higher-order partial differential operators associated with weakly viscous incompressible flow problems. Specifically, our numerical formulation employs a principle of virtual power (PVP)-based weak formulation that utilizes a stream-function field, which distinguishes it from the more commonly used bi-harmonic type formulations. The usage of a stream-function field ensures a pointwise divergence-free velocity condition, making the present method suitable for low to moderately high Reynolds number flow problems. In contrast to the bi-harmonic formulation, which is typically used for describing internal flow and requires special treatment of outflow boundary conditions, the PVP-based formulation is more general and does not require special treatment at the outflow boundary. It is also demonstrated that both bi-harmonic and PVP-based weak formulations yield identical results for internal flow problems. Our method employs non-uniform rational B-spline basis functions, and we present a simple stitching technique for imposing no-slip Dirichlet boundary conditions. Finally, we solve Poisson's equation to recover the pressure field. Furthermore, we use an appropriate Gaussian quadrature that is exact for splines to speed up the computation of various element matrices, especially for high polynomial degrees. The proposed formulation is evaluated by solving a set of numerical problems, including internal flow and channel flow problems.

Analysis of Fluid Flows in Bounded Domain with Particular Shape of a Cavity using Lattice Boltzmann Method

Background: The present work numerically models the incompressible, continuous phase, viscous flow of Newtonian fluid flow in a bounded domain of two-dimensional cavity that is driven by walls and contains grooves in the shape of squares on the lower wall. With the help of the mesoscopic lattice Boltzmann method (LBM) and D2Q9 square lattice model, simulation results are found stable and reliable. The flow physics of the problem by varying Reynolds number, the height and quantity of lower wall grooves, and other fluid flow characteristics within the bounded domain are studied in detail. It is seen that the effects of the groove heights and wavelengths on the fluid flow are structured within the bounded domain. The study is performed from low Re = 100 to high Re = 3200, with minimum two and maximum four-wavelength grooves evaluated on the bottom surface, each having a height of low 0.25 and high 0.75. Additionally, a thorough discussion of complicated vortex dynamics is provided regarding the input parameters and geometry. Objective: The current study aims to use mesoscopic LBM to analyze incompressible viscous fluid flows on complex geometries other than rectangular shapes. Methods: Mesoscopic approach of kinetic theory-based Lattice Boltzmann method (LBM) is implemented in the current work. The popular Single Relaxation Time Lattice Boltzmann method with D2Q9 square lattice model and second-order accurate boundary condition is adopted for the current study. Results: The numerical approach of LBM is used to simulate fluid flows in a 2D bounded domain with grooved bottom surfaces. The influence of different factors, such as the height of bottom-wall surface grooves, flow Reynolds number, and wavelength of these grooves on flow patterns, is then investigated. Conclusion: The numerical study of the bounded domain is considered, and the Reynolds number is varied from 100 to 3200, with two and four-wavelength grooves evaluated on the bottom surface, each having a height of 0.25 and 0.75. The impacts on the flow pattern both within and slightly above the grooves of the computational findings for different Reynolds numbers, groove heights, and groove wavelengths are evaluated. As the Reynolds number rises, the mixing phenomenon of fluid is shown to flow more quickly in the wall-driven enclosures.

On Exact Solutions to Problems of Steady Flow of Incompressible Viscous Fluid in a Channel and in a Circular Pipe

Turning to the well-known Couette flow, we again, but now in a new formulation, consider the steady motion of an incompressible viscous fluid in a channel with two parallel flat walls, one of which is stationary while the other moves at a constant speed in its plane. In the classical formulation of this problem on a layered (laminar) fluid flow in a channel, the wall velocity and the constant pressure drop along the channel are independent quantities and are considered predetermined. However, when a solid body moves along the fluid surface, the fluid particles are set in motion. This occurs due to the arising pressure generated in the fluid caused by the movement of the solid body. Therefore, here the pressure or pressure drop is not predetermined but is determined as a result of solving the corresponding two-dimensional boundary value problem of a steady flow of a viscous fluid in a channel. For this, the corresponding boundary value problem of a fluid flow in the channel is formulated based on the linearized Navier-Stokes equations obtained from the general nonlinear Navier-Stokes equations for a steady flow of a viscous fluid neglecting the convective terms, which is true for small Reynolds numbers. Using the Fourier integral transform method, an exact (closed) solution to this boundary value problem is constructed; the velocity and pressure components are determined. Then it is shown that the resulting solution, different from the known Couette flow solution, also satisfies the original nonlinear Navier-Stokes equations, the continuity equation, and boundary conditions. A comparative analysis of the new solution with the Couette flow solution is carried out. The well-known Hagen-Poiseuille problem is also discussed from the point of view of establishing the conditions for the laminar axisymmetric and steady flow of an incompressible viscous fluid in a straight circular pipe According to the study's results, to implement laminar flow under the indicated conditions, the pressure drop along the pipe must be constant. The study demonstrates that this condition is simultaneously sufficient, and therefore, the necessary and sufficient condition is established for a laminar flow in a straight circular pipe.