Shear flows of a new class of power-law fluids

We consider the flow of a class of incompressible fluids which are constitutively defined by the symmetric part of the velocity gradient being a function, which can be non-monotone, of the deviator of the stress tensor. These models are generalizations of the stress power-law models introduced and studied by J. Málek, V. Pr??a, K.R. Rajagopal: Generalizations of the Navier-Stokes fluid from a new perspective. Int. J. Eng. Sci. 48 (2010), 1907–1924. We discuss a potential application of the new models and then consider some simple boundary-value problems, namely steady planar Couette and Poiseuille flows with no-slip and slip boundary conditions. We show that these problems can have more than one solution and that the multiplicity of the solutions depends on the values of the model parameters as well as the choice of boundary conditions.     

Optimal control of time discrete two-phase flow governed by a diffuse interface model

We present results on optimal control of two-phase flows. The fluid is modeled by a thermodynamically consistent diffuse interface model and allows to treat fluids of different densities and viscosities. In earlier work we proposed an energy stable time discretization for this model that we now employ to derive existence of optimal controls for a time discrete optimal control problem. The control aim is to obtain a desired distribution of the two phases in the system. For this we investigate three control actions. We use tangential Dirichlet boundary control and distributed control. We further consider the inverse problem of finding an initial distribution such that the evolution over a given time horizon starting from this value is close to a desired distribution. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the development and generalizations of Cahn–Hilliard equations within a thermodynamic framework

We provide a thermodynamic basis for the development of models that are usually referred to as ??phase-field models?? for compressible, incompressible, and quasi-incompressible fluids. Using the theory of mixtures as a starting point, we develop a framework within which we can derive ??phase-field models?? both for mixtures of two constituents and for mixtures of arbitrarily many fluids. In order to obtain the constitutive equations, we appeal to the requirement that among all admissible constitutive relations that which is appropriate maximizes the rate of entropy production (see Rajagopal and Srinivasa in Proc R Soc Lond A 460:631?C651, 2004). The procedure has the advantage that the theory is based on prescribing the constitutive equations for only two scalars: the entropy and the entropy production. Unlike the assumption made in the case of the Navier?CStokes?CFourier fluids, we suppose that the entropy is not only a function of the internal energy and the density but also of gradients of the partial densities or the concentration gradients. The form for the rate of entropy production is the same as that for the Navier?CStokes?CFourier fluid. As observed earlier in Heida and Málek (Int J Eng Sci 48(11):1313?C1324, 2010), it turns out that the dependence of the rate of entropy production on the thermodynamical fluxes is crucial. The resulting equations are of the Cahn?CHilliard?CNavier?CStokes type and can be expressed both in terms of density gradients or concentration gradients. As particular cases, we will obtain the Cahn?CHilliard?CNavier?CStokes system as well as the Korteweg equation. Compared to earlier approaches, our methodology has the advantage that it directly takes into account the rate of entropy production and can take into consideration any constitutive assumption for the internal energy (or entropy).     

Further mathematical results concerning Burgers fluids and their generalizations

In this paper, we extend the earlier work by Quintanilla and Rajagopal (Math Methods Appl Sci 29: 2133?C2147, 2006) and establish qualitative new results for a proper generalization of Burgers?? original work that stems form a general thermodynamic framework. Such fluids have been used to describe the behavior of several geological materials such as asphalt and the earth??s mantle as well as polymeric fluids. We study questions concerning stability, uniqueness and continuous dependence on initial data for the solutions of the flows of these fluids. We show that if certain conditions are not satisfied by the material moduli, the solutions could be unstable. The spatial behavior of the solutions is also analyzed.     

On pressure boundary conditions for steady flows of incompressible fluids with pressure and shear rate dependent viscosities

We consider a class of incompressible fluids whose viscosities depend on the pressure and the shear rate. Suitable boundary conditions on the traction at the inflow/outflow part of boundary are given. As an advantage of this, the mean value of the pressure over the domain is no more a free parameter which would have to be prescribed otherwise. We prove the existence and uniqueness of weak solutions (the latter for small data) and discuss particular applications of the results.     

Examples of boundary layers associated with the incompressible Navier-Stokes equations

The author surveys a few examples of boundary layers for which the Prandtl boundary layer theory can be rigorously validated. All of them are associated with the incompressible Navier-Stokes equations for Newtonian fluids equipped with various Dirichlet boundary conditions (specified velocity). These examples include a family of (nonlinear 3D) plane parallel flows, a family of (nonlinear) parallel pipe flows, as well as flows with uniform injection and suction at the boundary. We also identify a key ingredient in establishing the validity of the Prandtl type theory, i.e., a spectral constraint on the approximate solution to the Navier-Stokes system constructed by combining the inviscid solution and the solution to the Prandtl type system. This is an additional difficulty besides the wellknown issue related to the well-posedness of the Prandtl type system. It seems that the main obstruction to the verification of the spectral constraint condition is the possible separation of boundary layers. A common theme of these examples is the inhibition of separation of boundary layers either via suppressing the velocity normal to the boundary or by injection and suction at the boundary so that the spectral constraint can be verified. A meta theorem is then presented which covers all the cases considered here.     

On uniqueness and time regularity of flows of power‐law like non‐Newtonian fluids

Large class of non‐Newtonian fluids can be characterized by index p, which gives the growth of the constitutively determined part of the Cauchy stress tensor. In this paper, the uniqueness and the time regularity of flows of these fluids in an open bounded three‐dimensional domain is established for subcritical ps, i.e. for p>11/5. Our method works for ‘all’ physically relevant boundary conditions, the Cauchy stress need not be potential and it may depend explicitly on spatial and time variable. As a simple consequence of time regularity, pressure can be introduced as an integrable function even for Dirichlet boundary conditions. Moreover, these results allow us to define a dynamical system corresponding to the problem and to establish the existence of an exponential attractor. Copyright © 2010 John Wiley & Sons, Ltd.     

AN ARTIFICIAL BOUNDARY CONDITION FOR THE INCOMPRESSIBLE VISCOUS FLOWS IN A NO-SLIP CHANNEL

1.IntroductionWhencomputingthenumericals0luti0nsofviscousfluidfl0wproblemsinallun-boundedd0main,0neoftenintroducesartificialboundaries,andsetsupanartificialbopundarycondition0nthem;thenthe0riginalproblemisreducedtoaproblemonab0undedc0mputationald0main.InordertoIimitthecomputatio11alcost,theseboundariesmustnotbet00farfromthedomainofinterest.Theref0re,theartificialboundaryc0nditi0nsmustbegoodapprotimationt0the"exact"boundaryconditions(sothatthes0lutionoftheproblemintheboundeddonlainisequaltothes…     

Approximate controllability for linearized compressible barotropic Navier–Stokes system in one and two dimensions

In this paper we consider compressible barotropic Navier–Stokes equations in one and two dimensions linearized around a constant steady state with Dirichlet boundary conditions. We explore the controllability of this linearized system using a control only for the velocity equation. We prove that the system with homogeneous Dirichlet boundary conditions, is approximately controllable by a localized interior control when time is sufficiently large.     

Numerical study of spherical Couette flows for certain zenith-angle-dependent rotations of boundary spheres at low Reynolds numbers

The numerical method with splitting of boundary conditions developed previously by the first and third authors for solving the stationary Dirichlet boundary value problem for the Navier-Stokes equations in spherical layers in the axisymmetric case at low Reynolds numbers and a corresponding software package were used to study viscous incompressible steady flows between two con-centric spheres. Flow regimes depending on the zenith angle ?? of coaxially rotating boundary spheres (admitting discontinuities in their angular velocities) were investigated. The orders of accuracy with respect to the mesh size of the numerical solutions (for velocity, pressure, and stream function in a meridional plane) in the max and L ₂ norms were studied in the case when the velocity boundary data have jump discontinuities and when some procedures are used to smooth the latter. The capabilities of the Richardson extrapolation procedure used to improve the order of accuracy of the method were investigated. Error estimates were obtained. Due to the high accuracy of the numerical solutions, flow features were carefully analyzed that were not studied previously. A number of interesting phenomena in viscous incompressible flows were discovered in the cases under study.     

General Nonlocal Model Describing the Laminar and Turbulent Flows of Viscous and Nonlinear Viscous Fluids and Its Investigation

A general nonlocal model describing the flows of viscous and nonlinear viscous fluids for both laminar and turbulent flows is introduced and studied. For this model, the viscosity of the fluid depends on the second invariant of the rate of the strain tensor and on a nonlocal (integral) characteristic of the flow. This characteristic is a vector that, in the simplest case, is an analog of the Reynolds number. For slow flows, the model turns into the Navier–Stokes equations or into the equations of a nonlinear viscous fluid. Problems on steady and nonsteady flows with mixed boundary conditions when velocities and surface forces are prescribed on different parts of the boundary are studied. Existence results without restrictions on the smallness of data and on the length of the interval of time are proved.     

二维不可压缩磁流体方程边界层分离的条件

该文研究了二维不可压缩磁流体方程的解,其中要求磁流体的速度满足Dirichlet边界条件、磁场在边界上的值与时间无关. 利用Taylor展开式和不可压缩流的结构分歧理论, 得到了磁流体方程发生边界层分离的条件, 它取决于外力、初值和磁场在边界上的取值, 并且该条件可以预测磁流体边界层分离发生的时间与地点.     

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 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.     

On non-Newtonian p-fluids. The pseudo-plastic case

In the following we study a class of stationary Navier-Stokes equations with shear dependent viscosity, under the non-slip (Dirichlet) boundary condition. We consider pseudo-plastic fluids. A fluid is said pseudo-plastic, or shear thinning, if in Eq. (1.1) below one has p<2. We are interested in global (i.e., up to the boundary) regularity results, in dimension n=3, for the second order derivatives of the velocity and the first order derivatives of the pressure. We consider a cubic domain Ω and impose the non-slip boundary condition only on two opposite faces. On the other faces we assume periodicity, as a device to avoid effective boundary conditions. This choice is made so that we work in a bounded domain Ω and simultaneously with a flat boundary.     

Thermally coupled,stationary, incompressible MHD flow; existence,uniqueness, and finite element approximation

This article addresses the questions of existence, uniqueness, and finite element approximation (including some computational aspects) of solutions to the equations of steady-state magnetohy-drodynamic (MHD) when buoyancy effects due to temperature differences in the flow cannot be neglected. We couple the MHD equations to the heat equation and employ the well-known Boussinesq approximation. We consider the equations posed on a bounded three-dimensional domain. The boundary conditions for the velocity are of Dirichlet type; the boundary conditions for the temperature are mixed (of Dirichlet type and of Neumann type); we also specify the normal component of the magnetic field and tangential component of the electric field on the boundary. We point out that these problems are relevant to many physical phenomena such as the cooling of nuclear reactors by electrically conducting fluids, continuous metal casting, crystal growth, and semi-conductor manufacture. © 1995 John Wiley & Sons, Inc.     

Stationary solutions to the equations of dynamics of mixtures of heat-conductive compressible viscous fluids

We consider the equations describing the three-dimensional steady motions of binary mixtures of heat-conductive compressible viscous fluids. An existence theorem for the boundary value problem that corresponds to flows in a bounded domain is proved in the class of weak generalized solutions.     

Exponential attractors for a class of reaction-diffusion problems with time delays

We consider a reaction-diffusion system subject to homogeneous Neumann boundary conditions on a given bounded domain. The reaction term depends on the population densities as well as on their past histories in a very general way. This class of systems is widely used in population dynamics modelling. Due to its generality, the longtime behavior of the solutions can display a certain complexity. Here we prove a qualitative result which can be considered as a common denominator of a large family of specific models. More precisely, we demonstrate the existence of an exponential attractor, provided that a bounded invariant region exists and the past history decays exponentially fast. This result will be achieved by means of a suitable adaptation of the l-trajectory method coming back to the seminal paper of Málek and Nečas. The first author was partially supported by the Italian PRIN 2006 research project Problemi a frontiera libera, transizioni di fase e modelli di isteresi. The second author was supported by the research project MŠM 0021620839 and by the project LC06052 (Jindřich Nečas Center for Mathematical Modeling).     

Nonlinear Eigenvalue Problems in the Stability Analysis of Morphogen Gradients

This paper is concerned with several eigenvalue problems in the linear stability analysis of steady state morphogen gradients for several models of Drosophila wing imaginal discs including one not previously considered. These problems share several common difficulties including the following: (a) The steady state solution which appears in the coefficients of the relevant differential equations of the stability analysis is only known qualitatively and numerically. (b) Though the governing differential equations are linear, the eigenvalue parameter appears nonlinearly after reduction to a problem for one unknown. (c) The eigenvalues are determined not only as solutions of a homogeneous boundary value problem with homogeneous Dirichlet boundary conditions, but also by an alternative auxiliary condition to one of the Dirichlet conditions allowed by a boundary condition of the original problem. Regarding the stability of the steady state morphogen gradients, we prove that the eigenvalues must all be positive and hence the steady state morphogen gradients are asymptotically stable. The other principal finding is a novel result pertaining to the smallest (positive) eigenvalue that determines the slowest decay rate of transients and the time needed to reach steady state. Here we prove that the smallest eigenvalue does not come from the nonlinear Dirichlet eigenvalue problem but from the complementary auxiliary condition requiring only to find the smallest zero of a rational function. Keeping in mind that even the steady state solution needed for the stability analysis is only known numerically, not having to solve the nonlinear Dirichlet eigenvalue problem is both an attractive theoretical outcome and a significant computational simplification.     

Acceleration of the Process of Entering Stationary Mode for Solutions of a Linearized System of Viscous Gas Dynamics. II

For finite-difference approximation of the linearized system of differential equations of viscous gas dynamics, the governing Dirichlet boundary conditions are constructed to guarantee the acceleration of the process of reaching the steady state solution. Necessary estimates are presented for the rate of convergence in the case of zero boundary conditions as well as the calculation results for stabilization in the case of initial conditions with jumps of pressure arid/or density.