Numerical solution of laminar incompressible generalized Newtonian fluids flow

This paper deals with the numerical solution of laminar viscous incompressible flows for generalized Newtonian fluids in the branching channel. The generalized Newtonian fluids contain Newtonian fluids, shear thickening and shear thinning non-Newtonian fluids. The mathematical model is the generalized system of Navier-Stokes equations. The finite volume method combined with an artificial compressibility method is used for spatial discretization. For time discretization the explicit multistage Runge-Kutta numerical scheme is considered. Steady state solution is achieved for t → ∞ using steady boundary conditions and followed by steady residual behavior. For unsteady solution a dual-time stepping method is considered. Numerical results for flows in two dimensional and three dimensional branching channel are presented.     

Numerical Solution of Newtonian and Non-Newtonian Flows in Bypass

This paper deals with a numerical solution of laminar incompressible steady flows of Newtonian and non–Newtonian fluids through bypass of a restricted vessel. Blood flow is considered to be Newtonian in the case of vessels of large diameters as aorta. On the other hand, with decreasing diameter of a vessel the non–Newtonian behavior of blood can play a significant role. One could describe these problems using Navier–Stokes equations and continuity equation as a model. In the case of Newtonian fluids one considers constant viscosity compared to non–Newtonian fluids where viscosity varies and can depend on the tensor of deformation. In order to find numerical solution, the system of equations is completed using an artificial compressibility method. The space derivatives are discretised using a cell centered finite volume method and arising system of ordinary differential equations is solved using an explicit multistage Runge–Kutta method with given steady boundary conditions. The steady solution is achieved for time t→∞ and steady boundary conditions. The results can be used in the field of cardiovascular research. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical solution of turbulent jet flows

The work deals with numerical modelling of several turbulent 3D jet flows: steady impinging jet, steady free jet in cross–flow, synthetic free jet (unsteady) and synthetic impinging jet (unsteady). The numerical method is based on artificial compressibility method with dual time extension for unsteady cases. Space discretization uses cell–centered finite volume method with third order accurate upwind approximation for convection, the time discretisations are implicit. Turbulence is modelled using two–equation eddy viscosity models and by explicit algebraic Reynolds stress model (EARSM by Wallin and Hellsten). The results of first three cases are compared with measurements. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical simulation of non-Newtonian blood flow in bypass models

The article presents the numerical investigation of non–Newtonian effects of steady blood flow in complete idealized 3–D bypass models, whose native artery is either coronary or femoral with average physiological parameters. Considering the blood to be a generalized Newtonian fluid, the shear–dependent viscosity is described by two well–known macroscopic non–Newtonian models (the Carreau–Yasuda model and the modified Cross model). The results were obtained by own developed computational software based on the pseudo–compressibility approach and on the cell–centred finite volume method defined on unstructured hexahedral grids. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical solution of steady and unsteady flow over a vibrating airfoil in a channel

The work deals with a numerical solution of 2D steady and unsteady inviscid incompressible flow over the profile NACA 0012 in a channel. The finite volume method (FVM) in a form of cell-centered explicit schemes at quadrilateral C-mesh is used. Governing system of equations is the system of incompressible Euler equations. The method of artificial compressibility and time dependent method is applied to steady computations. The small disturbance theory (SDT) applied to a numerical solution of flow over a rotated profile by a small angle only is mentioned. Brief introduction is given to the Arbitrary (Semi) Lagrangian-Eulerian (ALE) method used for unsteady computations. Some numerical results of unsteady flow over a vibrating profile achieved by both SDT and ALE method are presented. Unsteady flow is caused by prescribed oscillations of the profile (one degree of freedom) fixed in an elastic axis. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical Solution of Newtonian and Non-Newtonian Flows

This paper deals with the numerical solution of Newtonian and non-Newtonian flows. The flows are supposed to be laminar, viscous, incompressible and steady. The model used for non-Newtonian fluids is some variant of power-law. Governing equations in this model are incompressible Navier-Stokes equations. For numerical solution one could use artificial compressibility method with three stage Runge-Kutta finite volume method in cell centered formulation for discretization of space derivatives. Following cases of flows are solwed: flow through a bypass connected to main channel in 2D and 3D and non-Newtonian flow through branching channels in 2D. These results are presented for 2D and 3D case. This problem could have an application in the area of biomedicine. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Flow over a profile in a channel with dynamical effects

The work deals with a numerical solution of 2D inviscid incompressible flow over the profile NACA 0012 in a channel. The finite volume method in a form of cell‐centered scheme at quadrilateral C‐mesh is used. Governing system of equations is the system of Euler equations. Numerical results are partially compared with experimental data. Steady state solutions of the flow as well unsteady flows caused by prescribed oscillation of the profile were computed. The method of artificial compressibility and the time dependent method are used for computation of the steady state solution. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stationary Oberbeck–Boussinesq model of generalized Newtonian fluid governed by multivalued partial differential equations

A generalization of the Oberbeck–Boussinesq model consisting of a system of steady state multivalued partial differential equations for incompressible, generalized Newtonian of the p-power type, viscous flow coupled with the nonlinear heat equation is studied in a bounded domain. The existence of a weak solution is proved by combining the surjectivity method for operator inclusions and a fixed point technique.     

Approximation of the incompressible non‐Newtonian fluid equations by the artificial compressibility method

This paper studies the approximation of the non‐Newtonian fluid equations by the artificial compressibility method. We first introduce a family of perturbed compressible non‐Newtonian fluid equations (depending on a positive parameter ε) that approximates the incompressible equations as ε → 0⁺. Then, we prove the unique existence and convergence of solutions for the compressible equations to the solutions of the incompressible equations. Copyright © 2012 John Wiley & Sons, Ltd.     

Bivariate spline method for numerical solution of steady state Navier–Stokes equations over polygons in stream function formulation

We use the bivariate spline finite elements to numerically solve the steady state Navier–Stokes equations. The bivariate spline finite element space we use in this article is the space of splines of smoothness r and degree 3r over triangulated quadrangulations. The stream function formulation for the steady state Navier–Stokes equations is employed. Galerkin's method is applied to the resulting nonlinear fourth‐order equation, and Newton's iterative method is then used to solve the resulting nonlinear system. We show the existence and uniqueness of the weak solution in H²(Ω) of the nonlinear fourth‐order problem and give an estimate of how fast the numerical solution converges to the weak solution. The Galerkin method with C¹ cubic splines is implemented in MATLAB. Our numerical experiments show that the method is effective and efficient. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 147–183, 2000     

Analytical solution of unsteady three-dimensional MHD boundary layer flow and heat transfer due to impulsively stretched plane surface

The unsteady magnetohydrodynamic viscous flow and heat transfer of Newtonian fluids induced by an impulsively stretched plane surface in two lateral directions are studied by using an analytic technique, namely, the homotopy method. The analytic series solution presented here is highly accurate and uniformly valid for all time in the entire region. The effects of the stretching ratio and the magnetic field on the surface shear stresses and heat transfer are studied. The surface shear stresses in x- and y-directions and the surface heat transfer are enchanced by increasing stretching ratio for a fixed value of the magnetic parameter. For a fixed stretching ratio, the surface shear stresses increase with the magnetic parameter, but the heat transfer decreases. The Nusselt number takes longer time to reach the steady state than the skin friction coefficients. There is a smooth transition from the initial unsteady state to the steady state.     

A finite element–finite volume discretization of convection‐diffusion‐reaction equations with nonhomogeneous mixedboundary conditions: Error estimates

We consider a time‐dependent and a steady linear convection‐diffusion‐reaction equation whose coefficients are nonconstant. Boundary conditions are mixed (Dirichlet and Robin–Neumann) and nonhomogeneous. Both the unsteady and the steady problem are approximately solved by a combined finite element–finite volume method: the diffusion term is discretized by Crouzeix–Raviart piecewise linear finite elements on a triangular grid, and the convection term by upwind barycentric finite volumes. In the unsteady case, the implicit Euler method is used as time discretization. The ‐ and the ‐error in the unsteady case and the H¹‐error in the steady one are estimated against the data, in such a way that no parameter enters exponentially into the constants involved. © 2016Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1591–1621, 2016     

Unsteady flow of a Maxwell fluid with fractional derivatives in a circular cylinder moving with a nonlinear velocity

Abstract

In this paper we determine the velocity field and the shear stress corresponding to the unsteady flow of a Maxwell fluid with fractional derivatives driven by an infinite circular cylinder that slides along its axes with a velocity At^a. The general solutions, obtained by means of integral transforms, satisfy all imposed initial and boundary conditions. They can be easily particularized to give the similar solutions for ordinary Maxwell and Newtonian fluids. Finally, the influence of the parameters α and β on the fluid motion as well as a comparison between models is underlined by graphical illustrations.     

High-order discrete-time orthogonal spline collocation methods for singularly perturbed 1D parabolic reaction–diffusion problems

Quasi-optimal error estimates are derived for the continuous-time orthogonal spline collocation (OSC) method and also two discrete-time OSC methods for approximating the solution of 1D parabolic singularly perturbed reaction–diffusion problems. OSC with C¹ splines of degree r ≥ 3 on a Shishkin mesh is employed for the spatial discretization while the Crank–Nicolson method and the BDF2 scheme are considered for the time-stepping. The results of numerical experiments validate the theoretical analysis and also exhibit additional quasi-optimal results, in particular, superconvergence phenomena.     

Global dynamics of approximate solutions to an age-structured epidemic model with diffusion

We consider an age-dependent s-i-s epidemic model with diffusion whose mortality is unbounded. We approximate the solution using Galerkin methods in the space variable combined with backward Euler along the characteristic direction in the age and time variables. It is proven that the scheme is stable and convergent in optimal rate in l ^∞,2 (L ²) norm. To investigate the global behavior of the discrete solution resulting from the algorithm, we reformulate the resulting system into a monotone form. Positivity of the nonlocal birth process is proved using the positivity of the first eigenvalue of the resulting matrix system and using the fact that the positivity is preserved along the characteristics. The difference equation of the steady state coupled with nonlocal birth process is solved by developing monotone iterative schemes. The stability of the discrete solution of the steady state is then analyzed by constructing suitable positive subsolutions. Mathematics subject classifications (2000) 65M12, 65M25, 65M60, 92D25 M.-Y. Kim: This work was supported by Korea Research Foundation Grant (KRF-2001-041-D00037).     

Two‐Grid method for nonlinear parabolic equations by expanded mixed finite element methods

In this article, we develop a two‐grid algorithm for nonlinear reaction diffusion equation (with nonlinear compressibility coefficient) discretized by expanded mixed finite element method. The key point is to use two‐grid scheme to linearize the nonlinear term in the equations. The main procedure of the algorithm is solving a small‐scaled nonlinear equations on the coarse grid and dealing with a linearized system on the fine space using the Newton iteration with the coarse grid solution. Error estimation to the expanded mixed finite element solution is analyzed in detail. We also show that two‐grid solution achieves the same accuracy as long as the mesh sizes satisfy H = O(h^1/2). Two numerical experiments are given to verify the effectiveness of the algorithm. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Viscous Flows in a Half Space Caused by Tangential Vibrations on Its Boundary

The paper is devoted to the studies of viscous flows caused by a vibrating boundary. The fluid domain is a half‐space, its boundary is a nondeformable plane that exhibits purely tangential vibrations. Such a simple geometrical setting allows us to study general boundary velocity fields and to obtain general results. From a practical viewpoint, such boundary conditions may be seen as the tangential vibrations of the material points of a stretchable plane membrane. In contrast to the classical boundary layer theory, we aim to build a global solution. To achieve this goal we employ the Vishik–Lyusternik approach, combined with two‐timing and averaging methods. Our main result is: we obtain a uniformly valid in the whole fluid domain approximation to the global solutions. This solution corresponds to general boundary conditions and to three different settings of the main small parameter. Our solution always include the inner part and outer part that both contain oscillating and non‐oscillating components. It is shown that the nonoscillating outer part of the solution is governed either by the full Navier–Stokes equations or the Stokes equations (both with the unit viscosity) and can be interpreted as a steady or unsteady streaming. In contrast to the existing theories of a steady streaming, our solutions do not contain any secular (infinitely growing with the inner normal coordinate) terms. The examples of the spatially periodic vibrations of the boundary and the angular torsional vibrations of an infinite rigid disc are considered. These examples are still brief and illustrative, while the core of the paper is devoted to the adaptation of the Vishik–Lyusternik method to the development of the general theory of vibrational boundary layers.     

Continuous interior penalty finite element method for the time-dependent Navier–Stokes equations: space discretization and convergence

This paper focuses on the numerical analysis of a finite element method with stabilization for the unsteady incompressible Navier–Stokes equations. Incompressibility and convective effects are both stabilized adding an interior penalty term giving L ²-control of the jump of the gradient of the approximate solution over the internal faces. Using continuous equal-order finite elements for both velocities and pressures, in a space semi-discretized formulation, we prove convergence of the approximate solution. The error estimates hold irrespective of the Reynolds number, and hence also for the incompressible Euler equations, provided the exact solution is smooth.     

The existence,uniqueness, and regularity for an incompressible Newtonian flow with intrinsic degree of freedom

We consider the regularity and uniqueness of solution to the Cauchy problem of a mathematical model for an incompressible, homogeneous, Newtonian fluid, taking into account internal degree of freedom. We first show there exist uniquely a local strong solution. Then we show this solution can be extend to the whole interval [0,T] if the velocity u, or its gradient ? u, or the pressure p belongs to some function class, which are similar with that of incompressible Navier–Stokes equations. Our result shows that the solution is unique in these classes, and that velocity field plays a more prominent role in the existence theory of strong solution than the angular velocity field. Finally, if the L^{3 ∕ 2}‐norm of the initial angular velocity vector and some homogeneous Besov norm of initial velocity field are small, then there exists uniquely a global strong solution. Copyright © 2012 John Wiley & Sons, Ltd.     

An Existence Proof for the Gravitating BPS Monopole

We prove the existence of the gravitating BPS monopole in Einstein-Yang-Mills-Higgs (EYMH) theory. Existence is established using a Newtonian perturbation argument which shows that a Yang-Mills-Higgs BPS monopole solution can be be continued analytically in powers of 1/c² to an EYMH solution. Communicated by Sergiu Klainerman submitted 2/04/04, accepted 29/08/05