Incompressible Viscous Flow Research Articles

Simplified inverse distance weighting-immersed boundary method for simulation of fluid-structure interaction

When simulating fluid-structure interaction (FSI) problems involving moving objects, the implicit inverse distance weighting-immersed boundary method (IDW-IBM) developed by Du et al. [1] has to construct a large square correlation matrix and solve its inversion at each time step. In this work, a simplified inverse distance weighting-immersed boundary method (SIDW-IBM) is proposed to eliminate the intrinsic limitations in the original implicit IDW-IBM. Through error analysis using Taylor series expansion, a second order approximation can be derived, which allows us to approximate the large square correlation matrix into a diagonal matrix; thereby, we proposed the SIDW-IBM based on this second order approximation to explicitly evaluate the velocity corrections, where the needs to assemble the large correlation matrix and inverse it are circumvented. Owing to the fact that the inverse distance weighting interpolation removes the limitations in the Dirac delta function, the proposed SIDW-IBM has been successfully implemented on the non-uniform meshes to further improve the computational efficiency. The proposed SIDW-IBM is integrated with the reconstructed lattice Boltzmann flux solver (RLBFS) [2] to simulate some classic incompressible viscous flows, including flow past an in-line oscillating cylinder, flow past a heaving airfoil, and flow past a three-dimensional flexible plate. The good agreement between the present results and reference data demonstrates the capability and feasibility of the SIDW-IBM for simulating FSI problems with moving boundaries and large deformations.

An improved pressure gradient method for viscous incompressible flows

The pressure gradient method solves the viscous incompressible flow with the pressure gradients, rather than the pressure, as unknown variables. Two variants of the pressure gradient method have been developed in the past but have not received much attention due to their unsatisfactory performance or implementation complexity. Based on the artificial compressibility concept, this study proposes an improved pressure gradient method. One distinct feature of this method is that it requires no pressure/pressure gradient boundary condition or special treatment on wall boundaries. An auxiliary variable is introduced to represent the velocity dilatation, greatly simplifying the spatial discretization and computational procedure. The mathematical formulations are elaborated and compared with the previous pressure gradient methods, followed by discussions of compatibility relationships, boundary condition setup, and an extension to a pressure Poisson-like equation. Four validation examples are performed for various flow scenarios, and the solutions and domain solenoidity are examined for each case. The study also compares associated computational methods, different pressure boundary conditions, and flow characteristics, demonstrating the benefits of the present method.

Estimation of added effects and their frequency dependence in various fluid–structure interaction problems

This paper focuses on a fluid–structure interaction topic—the determination of added effects caused by fluid forces acting on a body, considering the standard linear equation of motion. We present various problems that assume small-displacement oscillations of single and multiple bodies in inviscid irrotational (potential) flow or viscous incompressible flow in both closed domain and external flow. For inviscid flow, effects of geometric parameters on the added effects were studied. The presented results extend results known from the literature. For viscous flow, frequency dependence of the added effects was studied for a wide range of frequency. The added effects were computed from data from numerical simulations of fluid flow, where the body oscillations were modeled using the dynamic mesh approach. Effects of the phase shift caused by the dynamic mesh were addressed. The added mass was compared with the corresponding value determined for inviscid flow where applicable. The results show strong dependence of the added effects on many parameters, making their proper computation challenging even for simplified cases.

Метод 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

Study on the integrated self-priming process of a vehicle-mounted self-priming pump

In order to explore the self-priming characteristics of the self-priming pump at the mobile pump truck, this paper established a complete three-dimensional circulatory piping system including the self-priming pump, tank, valves, inlet pipe and outlet pipe. The UDF(User Defined Functions) was used to realize the acceleration-constant speed operation process of the impeller, thus reflecting the actual changing state of the rotational speed. Based on the VOF(Volume Of Fluid) multiphase flow model and the Realizable k-ε turbulence model, a coupled numerical calculation of unsteady incompressible viscous flow was conducted for its self-priming process. The results show that the self-priming process of the pump can be roughly divided into four stages: the rapid suction stage, the shock exhaust stage, the rapid exhaust period and the pump residual gas discharge stage. The proportion of each stage in the total self-priming time showed an increasing trend. During the rapid suction stage, the water level in the vertical section of the inlet pipe showed a slow and then fast-rising pattern. During the shock exhaust stage, the average gas-phase volume fraction in the volute is lower than that of the impeller, and the gas content at the volute outlet is lower than that of the impeller inlet. The region at the inlet and outer edge of the impeller consistently experience significant energy losses.

A high‐order pure streamfunction method in general curvilinear coordinates for unsteady incompressible viscous flow with complex geometry

AbstractIn this paper, a high‐order compact finite difference method in general curvilinear coordinates is proposed for solving unsteady incompressible Navier‐Stokes equations. By constructing the fourth‐order spatial discretization schemes for all partial derivative terms of the pure streamfunction formulation in general curvilinear coordinates, especially for the fourth‐order mixed derivative terms, and applying a Crank‐Nicolson scheme for the second‐order temporal discretization, we extend the unsteady high‐order pure streamfunction algorithm to flow problems with more general non‐conformal grids. Furthermore, the stability of the newly proposed method for the linear model is validated by von‐Neumann linear stability analysis. Five numerical experiments are conducted to verify the accuracy and robustness of the proposed method. The results show that our method not only effectively solves problems with non‐conformal grids, but also allows grid generation and local refinement using commercial software. The solutions are in good agreement with the established numerical and experimental results.

A numerical model for the simulation of complex planar Newtonian interfaces

We present a numerical model for the simulation of complex planar interfaces at which moving solid objects can be immersed, reproducing a wide variety of experimental conditions. The mathematical model consists of the Navier-Stokes equations governing the incompressible viscous flow in the liquid subphase, the transport equation for the evolution of the surfactant concentration at the interface, and the interfacial stress balance equation. The equations are simplified by treating the problem as isothermal and the surfactant as insoluble. The bulk flow equations are discretized using a collocated finite volume method, while the interfacial flow equations are discretized using a finite area method. The Boussinesq-Scriven interface constitutive model and a variant form accounting for extensional viscosity are used to describe the extra surface stress tensor. The coupling between surfactant concentration, interfacial velocity, and bulk velocity is treated implicitly by solving the interfacial and bulk equations sequentially at each time step until a stopping criterion is satisfied. The motion of the solid is treated by an arbitrary Lagrangian-Eulerian method. The model has been implemented in the OpenFOAM framework and allows the incorporation of new interface models and solvers, making the developed new package a versatile and powerful tool in the field of computational rheology. Applications of the model include the numerical simulation of flow around objects, such as probes, immersed at a complex interface, reproducing given experimental conditions, and its use as a tool in the analysis and design of interfacial stress rheometers. Several test cases have been performed to validate the model by comparing the results obtained with analytical solutions and with numerical and experimental results available in the literature.

Hybrid ferrofluid flow on a stretching sheet with Stefan blowing and magnetic polarization effects in a porous medium

PurposeThis study investigates the behavior of viscous hybrid ferromagnetic fluids flowing through plain elastic sheets with the magnetic polarization effect. It examines flow in a porous medium using Stefan blowing and utilizes a versatile hybrid ferrofluid containing MnZnFe2O4 and Fe3O4 nanoparticles in the C2H2F4 base fluid, offering potential real-world applications. The study focuses on steady, laminar and viscous incompressible flow, analyzing heat and mass transfer aspects, including thermal radiation, Brownian motion, thermophoresis and viscous dissipation with convective boundary condition.Design/methodology/approachThe governing expression of the flow model is addressed with pertinent non-dimensional transformations, and the finite element method solves the obtained system of ordinary differential equations.FindingsThe variations in fluid velocity, temperature and concentration profiles against all the physical parameters are analyzed through their graphical view. The association of these parameters with local surface friction coefficient, Nusselt number and Sherwood number is examined with the numerical data in a table.Originality/valueThis work extends previous research on ferrofluid flow, investigating unexplored parameters and offering valuable insights with potential engineering, industrial and medical implications. It introduces a novel approach that uses mathematical simplification techniques and the finite element method for the solution.

Numerical modelling of aerodynamic response to gusts and gust effect mitigation

Unmanned aerial vehicle (UAV) flights in urban environments are challenging due to the complex flow structures and elevated turbulence around buildings. Consequently, research has shifted towards investigating the impact of gusts on the aerodynamic stability and control of UAVs. This study focuses on enhancing gust numerical modelling capabilities to understand the aerodynamic response, specifically exploring gust mitigation strategies for UAVs operating in turbulent urban environments. The split-velocity method, originally designed for two-dimensional compressible inviscid flows, where the velocity components were decomposed into a prescribed gust velocity and the remaining velocity components, is extended to three-dimensional incompressible viscous flows. To facilitate effective gust mitigation techniques, a radial basis function is applied to the modified split-velocity method to numerically model wings in pitching motions under gust encounters. A novel strategy is proposed to correct the discretized gust velocities and ensure gust flux conservation, showing effective improvement to the numerical predictions. The computed results agreed well with the experimental data available in the public domain, confirming that wing pitch motion can effectively mitigate the effects of gusts.

Heat transfer analysis of mixed convection boundary layer blood base nanofluids with the influence of viscous dissipation

In this research study, we delve into the analysis of two-dimensional boundary layer incompressible viscous flow of blood based nanofluids, incorporating a couple stress model. This investigation encompasses considerations of heat transmission and viscous dissipation effects. Specifically, focus on a particular type of nanoparticle, Magnesium oxide (MgO), dispersed in blood, the base fluid. By employing suitable similarity transformations, transform a set of partial differential equations (PDEs) into nonlinear ordinary differential equations (ODEs). To leverage the Homotopy Analysis Approach (HAM) to solve this system of equations. This study discerns the impacts of various factors on temperature and velocity distributions. Visualization through graphs depicting nanoparticle volume concentration, magnetic field strength, couple stress parameter, Eckert number, and Prandtl number aids in presenting our findings. These insights illuminate the intricate interplay among different parameters and offer insights for enhancing heat transfer in blood-based nanofluids. Additionally, scrutinize the thermal performance of the blood-based nanofluid, assessing the Nusselt number and local skin friction coefficient. Such analysis holds significance for ensuring patient safety and treatment efficacy, particularly in medical procedures involving blood circulation, where mixed convection plays a pivotal role. Understanding these dynamics is crucial for the design and optimization of various medical devices and therapies, ensuring optimal temperature maintenance during blood circulation-related operations.

A principle of maximum entropy for the Navier–Stokes equations

A principle of maximum entropy is proposed in the context of viscous incompressible flow in Eulerian coordinates. The relative entropy functional, defined over the space of L2 divergence-free velocity fields, is maximized relative to alternate measures supported over the energy–enstrophy surface. Since thermodynamic equilibrium distributions are characterized by maximum entropy, connections are drawn with stationary statistical solutions of the incompressible Navier–Stokes equations. Special emphasis is on the correspondence with the final statistics described by Kolmogorov’s theory of fully developed turbulence.

A face-centred finite volume method for laminar and turbulent incompressible flows

This work develops, for the first time, a face-centred finite volume (FCFV) solver for the simulation of laminar and turbulent viscous incompressible flows. The formulation relies on the Reynolds-averaged Navier–Stokes (RANS) equations coupled with the negative Spalart–Allmaras (SA) model and three novel convective stabilisations, inspired by Riemann solvers, are derived and compared numerically. The resulting method achieves first-order convergence of the velocity, the velocity-gradient tensor and the pressure. FCFV accurately predicts engineering quantities of interest, such as drag and lift, on unstructured meshes and, by avoiding gradient reconstruction, the method is less sensitive to mesh quality than other FV methods, even in the presence of highly distorted and stretched cells. A monolithic and a staggered solution strategies for the RANS-SA system are derived and compared numerically. Numerical benchmarks, involving laminar and turbulent, steady and transient cases are used to assess the performance, accuracy and robustness of the proposed FCFV method.

The Impact of Installation Angle on the Wind Load of Solar Photovoltaic Panels

In order to explore the wind load characteristics acting on solar photovoltaic panels under extreme severe weather conditions, based on the Shear Stress Transport (SST) κ-ω turbulence model, numerical calculations of three-dimensional incompressible viscous steady flow were performed for four installation angles and two extreme wind directions of the solar photovoltaic panels. The wind load characteristics on both sides of the photovoltaic panels were obtained, and the vortex structure characteristics were analyzed using the Q criterion. The results indicate that, under different installation angles, the windward side pressure of the solar photovoltaic panel is generally higher than the leeward side. The leeward side is prone to forming larger vortices, increasing the fatigue and damage risk of the material, which significantly impacts the solar photovoltaic panel. As the installation angle increases, the windward side pressure of the solar photovoltaic panel also gradually increases. Therefore, optimal installation methods include installing the panel facing the wind at angles of 30° and 45°, or installing it facing away from the wind at a 60° angle, to minimize the impact of wind load on the solar photovoltaic panel.

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.

Modeling and simulation of pollution transport in the Mediterranean Sea using enriched finite element method

This paper presents a novel numerical method for simulating the transport and dispersion of pollutants in the Mediterranean sea. The governing mathematical equations consist of a barotropic ocean model with friction terms, bathymetric forces, Coriolis and wind stresses coupled to an advection–diffusion equation with anisotropic dispersion tensor and source terms. The proposed numerical solver uses a multilevel adaptive semi-Lagrangian finite element method that combines various techniques, including the modified method of characteristics, finite element discretization, coupled projection scheme based on a rotational pressure correction algorithm, and an adaptive L2-projection. The approach employs the gradient of the concentration as an error indicator for enrichment adaptations and increasing the number of quadrature points where needed without refining the mesh. The method is shown to provide accurate and efficient simulations for pollution transport in the Mediterranean sea. The proposed approach distinguishes itself from the well-established adaptive finite element methods for incompressible viscous flows by retaining the same structure and dimension of linear systems during the adaptation process.

A Fourier spectral immersed boundary method with exact translation invariance, improved boundary resolution, and a divergence-free velocity field

This paper introduces a new immersed boundary (IB) method for viscous incompressible flow, based on a Fourier spectral method for the fluid solver and on the nonuniform fast Fourier transform (NUFFT) algorithm for coupling the fluid with the immersed boundary. The new Fourier spectral immersed boundary (FSIB) method gives improved boundary resolution in comparison to the standard IB method. The interpolated velocity field, in which the boundary moves, is analytically divergence-free. The FSIB method is gridless and has the meritorious properties of volume conservation, exact translation invariance, conservation of momentum, and conservation of energy. We verify these advantages of the FSIB method numerically both for the Stokes equations and for the Navier-Stokes equations in both two and three space dimensions. The FSIB method converges faster than the IB method. In particular, we observe second-order convergence in various problems for the Navier-Stokes equations in three dimensions. The FSIB method is also computationally efficient with complexity of O(N3log⁡(N)) per time step for N3 Fourier modes in three dimensions.

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.

Numerical analysis of velocity slip and magnetic Soret effect on dissipative Maxwell liquid motion about a stretching wall with heat source/sink

The main aim of present numerical analysis is to explore the effects of viscous dissipation and velocity slip on the two-dimensional steady-state viscous incompressible flow of Maxwell fluid over a stretching sheet under the infleunce of magnetic field. Novel collective effects of heat source/sink and Soret impacts have been included in the governing equations to explore the salient features of heat and mass transport mechanism. Further, the velocity slip is introduced into the boundary conditions to depict the momentum transport behavior in the flow region. A typical Maxwell fluid model is deployed to separate the non-Newtonian flow behavior from that of classical Newtonian fluids. However, the visco-elastic Maxwell fluid finds its abundant applications in various fields of engineering and medicine such as plastic sheet drawing, metallurgy, polymer sheet extrusion, chemical engineering, electronic cooling, injection molding and nuclear cooling and many others. Inspired by these applications, authors have made an attempt to disclose the magneto-thermo salient features of slip flow of Maxwell fluid over stretching wall under Soret effect. The current physical situation results the highly nonlinear coupled two-dimensional steady-state partial differential equations and which are not amenable to any of the direct techniques. Due to this, the suitable similarity transformations are deployed to reduce the dimensional complexity of the constructed partial differential equations. A robust Matlab-based bvp4c scheme is utilized as a basic tool to produce the similarity solutions of the governing equations. The computer generated numerical data is presented in view of flow, temperature and concentration profiles including physical quantities of interest like skin-friction coefficient, heat and mass transport rates for the various values of physical parameters range such as, 0 ≤ M ≤ 1 , 0 ≤ β ≤ 1 , − 0.3 ≤ λ ≤ 0.6 , − 0.2 ≤ Q ≤ 0.2 , 0 ≤ Ec ≤ 2.2 , 1.2 ≤ Pr ≤ 2.2 , 0 ≤ Sr ≤ 0.5 and 1.2 ≤ Sc ≤ 2.2 . Accordingly, it is noticed that, enhancing magnetic parameter decreases the Maxwell fluid flow and amplifies the temperature and concentration fields. Enlarging Maxwell fluid parameter decreases the velocity field and increases the temperature and concentration profiles. Amplified slip parameter diminished the Maxwell flow inside the flow regime. Increasing heat source/sink parameter amplifies the temperature field. Further, magnifying Soret number amplifies the concentration field. Magnitude of skin-friction coeffiecient enlarged with amplified magnetic and Maxwell fluid parameters. Heat and mass transport rates decayed with enhancement in heat source/sink and Soret parameters. Finally, it is described that, the present similarity solutions showing good agreement with the previously reported solutions in the literature and this fact confirms the accurateness of the used numerical method and the generated similarity solutions.

A boundary condition-enhanced direct-forcing immersed boundary method for simulations of three-dimensional phoretic particles in incompressible flows

In this paper we propose an improved three-dimensional immersed boundary method coupled with a finite-difference code to simulate self-propelled phoretic particles in viscous incompressible flows. We focus on the phenomenon of diffusiophoresis which, using the driving of a concentration gradient, can generate a slip velocity on a surface. In such a system, both the Dirichlet and Neumann boundary conditions are involved. In order to enforce the boundary conditions, we propose two improvements to the basic direct-forcing immersed boundary method. The main idea is that the immersed boundary terms are corrected by adding the force of the previous time step, in contrast to the traditional method which relies only on the instantaneous forces in each time step. For the Neumann boundary condition, we add two auxiliary layers inside the body to precisely implement the desired concentration gradient. To verify the accuracy of the improved method, we present problems of different complexity: The first is the pure diffusion around a sphere with Dirichlet and Neumann boundary conditions. Then we show the flow past a fixed sphere. In addition, the motion of a self-propelled Janus particle in the bulk and the spontaneously symmetry breaking of an isotropic phoretic particle are reported. The results are in very good agreements with the data that are reported in previously published literature.

Numerical heat and mass transfer in magnetized Williamson liquid motion over a sensor wall engulfed in a micro-cantilever system: An application to thin-film fluidic cells

Central focus of this numerical analysis is to investigate the influence of magnetic body force and thermal radiation on the time-dependent electrically conducting two-dimensional viscous incompressible boundary layer flow of non-Newtonian Williamson fluid over a microcantilever sensor surface suspended in a squeezed regime numerically. A novel heat source/sink and Soret effects are included in the respective governing equations. The constructed time-dependent nonlinear coupled partial differential equations are solved by using a bobust RK-4 technique. Increasing magnetic number increases the velocity distribution and decreases the temperature and concentration profiles in the flow regime.Enhancing Soret number increases the mass diffusion profile. Increasing magnetic number suppress the skin-friction coefficient. Heat transport rate booted with rising radiation number and Sherwood number decreased with enlarged Schmidt number in the channel regime. Finally, the numerical accuracy of the present similarity solutions and used numerical method is validated with existing results and observed a remarkable agreement.