Transition of unsteady flows of evaporation to steady state

We investigate the half-space problem of evaporation and condensation in the scope of discrete kinetic theory. Exact solutions are found to the boundary value problem and the initial boundary value problems of the flow in the half space for a discrete velocity model. The results are used to analyze the transition of the unsteady solutions towards steady states. To cite this article: A. d'Almeida, C. R. Mecanique 336 (2008).     

Discussion of Antiplane Singularities of Piezoelectric–Dielectric and Piezoelectric–Conductor Wedges

This technical note deals with two special topics from our previous paper (Chue and Chen in Arch Appl Mech 72 673–685, 2003) in Archive of Applied Mechanics: the effects of electrical conditions imposed on the edges and bonded interfaces of piezoelectric–dielectric and piezoelectric–conductor wedges on antiplane problems. After employing relatively realistic electrical conditions, we found that stress and electric displacement singularities are altered when boundary conditions and/or continuity conditions are changed, and we compared the results with those of previous studies.     

Temperature distribution in a Newtonian fluid injected between two semi-infinite plates

The analysis of the temperature distribution in time and place of a hot heat-conducting Newtonian fluid injected between two cooled parallel plates is presented. The 2-dimensional flow has a free flow front moving with constant velocity. The kernel of the fluid remains almost at the inlet temperature, but at the walls boundary layers occur with steeply descending temperature. The inner solutions inside these boundary layers are determined. To this end, the total region is divided into three distinct regions: the region G_I far behind the flow front, the flow front region G_II, and the intermediate region G_III between G_I and G_II. The asymptotics owing to each region are presented. The fundamental small parameter here is the thickness-to-length ratio of the 2-dimensional flow region. In most of the cases, similarity solutions are found. In the flow front region, for the formulation of the inner solution a Wiener–Hopf technique is used. Via matching procedures, the separate boundary layers are linked to each other to form one global boundary layer for the whole front region. All calculations in this paper are performed by analytical means, and all results are in analytical form. Comparison of our results with numerical solutions shows perfect agreement.     

Optimal Control of Turbulent Flow as Measure Solutions

Abstract

We use the concept of relaxed or measure-valued solutions for control problems of turbulent flow related to Navicr-Stokes equation. Sufficient conditions guaranteeing the existence of measure solutions arc presented. Results on existence of optimal controls for Blow up time and Bolza problems for such systems are also presented. New results on relaxed necessary conditions of optimality are proved. Further it is shown that the relaxed necessary conditions reduce to classical Pontryagin type necessary conditions if measure solutions degenerate into Dirac structure. The paper is concluded with an algorithm based on the new necessary conditions for computing optimal controls.     

Flow of viscoelastic fluids through a porous channel—I

The flow of viscoelastic fluids through a porous channel with one impermeable wall is computed. The flow is characterized by a boundary value problem in which the order of the differential equation exceeds the number of boundary conditions. Three solutions are developed: (i) an exact numerical solution, (ii) a perturbation solution for small R, the cross-flow Reynold's number and (iii) an asymptotic solution for large R. The results from exact numerical integration reveal that the solutions for a non-Newtonian fluid are possible only up to a critical value of the viscoelastic fluid parameter, which decreases with an increase in R. It is further demonstrated that the perturbation solution gives acceptable results only if the viscoelastic fluid parameter is also small. Two more related problems are considered: fluid dynamics of a long porous slider, and injection of fluid through one side of a long vertical porous channel. For both the problems, exact numerical and other solutions are derived and appropriate conclusions drawn.     

A consistent splitting scheme for unsteady incompressible viscous flows I. Dirichlet boundary condition and applications

A well‐recognized approach for handling the incompressibility constraint by operating directly on the discretized Navier–Stokes equations is used to obtain the decoupling of the pressure from the velocity field. By following the current developments by Guermond and Shen, the possibilities of obtaining accurate pressure and reducing boundary‐layer effect for the pressure are analysed. The present study mainly reports the numerical solutions of an unsteady Navier–Stokes problem based on the so‐called consistent splitting scheme (J. Comput. Phys. 2003; 192 :262–276). At the same time the Dirichlet boundary value conditions are considered. The accuracy of the method is carefully examined against the exact solution for an unsteady flow physics problem in a simply connected domain. The effectiveness is illustrated viz. several computations of 2D double lid‐driven cavity problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Finite element implementation of boundary conditions for the pressure Poisson equation of incompressible flow

In this paper we address the problem of the implementation of boundary conditions for the derived pressure Poisson equation of incompressible flow. It is shown that the direct Galerkin finite element formulation of the pressure Poisson equation automatically satisfies the inhomogeneous Neumann boundary conditions, thus avoiding the difficulty in specifying boundary conditions for pressure. This ensures that only physically meaningful pressure boundary conditions consistent with the Navier-Stokes equations are imposed. Since second derivatives appear in this formulation, the conforming finite element method requires C¹ continuity. However, for many problems of practical interest (i.e. high Reynolds numbers) the second derivatives need not be included, thus allowing the use of more conventional C⁰ elements. Numerical results using this approach for a wall-driven contained flow within a square cavity verify the validity of the approach. Although the results were obtained for a two-dimensional problem using the p-version of the finite element method, the approach presented here is general and remains valid for the conventional h-version as well as three-dimensional problems.     

Empirical test of a microscopic three-phase traffic theory

A review of dynamic nonlinear features of spatiotemporal congested patterns in freeway traffic is presented. The basis of the review is a comparison of theoretical features of the congested patterns that are shown by a microscopic traffic flow model in the context of the Kerner's three-phase traffic theory and empirical microscopic and macroscopic pattern characteristics measured on different freeways over various days and years. In this test of the microscopic three-phase traffic flow theory, a model of an "open" road is applied: Empirical time-dependence of traffic demand and drivers' destinations are used at the upstream model boundaries. At downstream model boundary conditions for vehicle freely leaving a modeling freeway section(s) are given. Spatiotemporal congested patterns emerge, develop, and dissolve in this open freeway model with the same types of bottlenecks as those in empirical observations. It is found that microscopic three-phase traffic models can explain all known macroscopic and microscopic empirical congested pattern features (e.g., probabilistic breakdown phenomenon as a first-order phase transition from free flow to synchronized flow, moving jam emergence in synchronized flow rather than in free flow, spatiotemporal features of synchronized flow and general congested patterns at freeway bottlenecks, intensification of downstream congestion due to upstream congestion at adjacent bottlenecks). It turns out that microscopic optimal velocity (OV) functions and time headway distributions are not necessarily qualitatively different, even if local congested traffic behavior is qualitatively different. Model performance with respect to spatiotemporal pattern emergence and evolution cannot be tested using these traffic characteristics. The reason for this is that important spatiotemporal features of congested traffic patterns are lost in these and many other macroscopic and microscopic traffic characteristics, which are widely used as the empirical basis for a test of traffic flow models. PACS: 89.40. + k, 47.54. + r, 64.60.Cn, 64.60.Lx     

Corner Singularities and Regularity Results for the Reissner/Mindlin Plate Model

The theory for elliptic boundary value problems for general elliptic systems is used in order to investigate systematically corner singularities and regularity for weak solutions to a broad class of boundary value problems for the Reissner/Mindlin plate model in polygonal domains. The regularity results for the deflection of the midplane and for the rotation of fibers normal to the midplane are formulated in Sobolev spaces H ^s , where s>1 is a real number. The number s depends on the geometry, the material parameters and the boundary conditions in general and is related to a decomposition of the fields in a singular and a regular part. The leading singular terms are calculated for a wide class of boundary conditions (36 combinations). The results are critically compared with those known from a stress potential approach.     

ON THE WELL POSEDNESS OF INITIAL VALUE PROBLEM FOR EULER EQUATIONS OF INCOMPRESSIBLE INVISCID FLUID(Ⅱ)

The solvability of the Euler equations about incompressible inviscid fluid based on the stratification theory is discussed. And the conditions for the existence of formal solutions and the methods are presented for calculating all kinds of ill-posed initial value problems. Two examples are given as the evidences that the initial problems at the hyper surface does not exist any unique solution. Foundation item: the National Natural Science Foundation of China (19971054) Biography: Shen Zhen (1977−)     

Convection regime flow in a vertical slot: continuum of solutions from capped to open ends

In a recent publication Bühler (Heat Mass Transfer 39:631–638, 2003) reported new results for conduction regime flow between vertical differentially-heated walls that provide a continuum of solutions between capped and open ends. In this paper we extend Bühler’s work to realize a continuum of solutions of convection regime flow using empirical results for the vertical temperature gradient that develops in tall aspect ratio geometries. The mass flux is determined analytically for this three-parameter family of solutions. Identical viscous and thermal boundary layers exist at the opposing walls when the cavity is capped. However, as the flow evolves to one with open ends, there is an intensification (attenuation) of the boundary layers near the hot (cold) walls. In the limit corresponding to an open-ended cavity, the boundary layer at the cold wall vanishes altogether.     

Uniqueness of the self-similar solutions of three-dimensional laminar boundary-layer equations

The equations of the three-dimensional laminar boundary layer on lines of flow outflow and inflow are studied for conical outer flow under the assumption that the Prandtl number and the productρμ are constant. It is shown that in the case of a positive velocity gradient of the secondary flow (α₁>0) the additional conditions which result from the physical flow pattern determine a unique solution of the system of boundary-layer equations. For a negative velocity gradient of the secondary flow (α₁≤0) these conditions are satisfied by two solutions. An approximate solution is obtained for the boundary layer equations which is in rather good agreement with the numerical integration results. Compressible gas flow in a three-dimensional laminar boundary layer is described by a system of nonlinear differential equations whose solution is not unique for given boundary conditions. Therefore additional conditions resulting from the physical pattern of the gas flow are imposed on the resulting solution. In the solution of problems with a negative pressure gradient these additional conditions are sufficient for a unique selection of the solution of the boundary-layer equations. However, in the case of a positive pressure gradient the solution of the boundary-layer equations satisfying the boundary and additional conditions may not be unique. In particular, in [1] in a study of a three-dimensional laminar boundary layer in the vicinity of the stagnation point it was shown that for $$c = {{\frac{{\partial v_e }}{{\partial y}}} \mathord{\left/ {\vphantom {{\frac{{\partial v_e }}{{\partial y}}} {\frac{{\partial u_e }}{{\partial x}}}}} \right. \kern-\nulldelimiterspace} {\frac{{\partial u_e }}{{\partial x}}}} > 0$$ the solution is unique, while for c<0 there are two solutions. In the present paper we study the question of the uniqueness of the self-similar solution of the three-dimensional laminar boundary-layer equations on lines of flow outflow and inflow for a conical outer flow.     

Effects of Fracture Boundary Conditions on Matrix-fracture Transfer Shape Factor

The matrix-fracture transfer shape factor is one of the important parameters in modeling naturally fractured reservoirs. Four decades after Warren and Root (1963, SPEJ, 245–255.) introduced the double porosity concept and suggested a relation for it, this parameter is still not completely understood. Even for a single-phase flow problem, investigators report different shape factors. This study shows that for a single-phase flow in a particular matrix block, the shape factor that Warren and Root defined is not unique and depends on the pressure in the fracture and how it changes with time. We use the Laplace domain analytical solutions of the diffusivity equation for different geometries and different boundary conditions to show that the shape factor depends on the fracture pressure change with time. In particular, by imposing a constant fracture pressure as it is typically done, one obtains the shape factor that Lim and Aziz (1995, J. Petrolean Sci. Eng. 13, 169.) calculated. However, other shape factors, similar to those reported in other studies are obtained, when other boundary conditions are chosen. Although, the time variability of the boundary conditions can be accounted for by the Duhamel’s theorem, in practice using large time-steps in numerical simulations can potentially introduce large errors in simulation results. However, numerical simulation models make use of a stepwise approximation of this theorem. It is shown in this paper that this approximation could lead to large errors in matrix-fracture transfer rate if large time-steps are chosen.     

On the linking up between Bingham fluid and plugged flow

When Bingham fluid is in motion, plugged flow often occurs at places far from the boundary walls. As there is not a decisive formula of constitutive relation for plugged flow, in some problems the solutions obtained may be indefinite. In this paper, annular flow and pipe flow are discussed, and unique solution is obtained in each case by utilizing the analytic property of shear stress. The solutions are identical in form with the commonly used formula for the pressure drop of mud flow in petroleum engineering.     

The analysis of the channel flow sensitivity to the parameters of the k–ε method

This paper deals with the problem of using sensitivity analysis for fluid mechanics solutions to the constants of the standard k–ε method for 2D, incompressible and steady flows. The problem is described and analysed on the basis of a channel flow. Sensitivity coefficients of the following properties were determined: a pressure, two components of a velocity, a turbulence kinetic energy, a dissipation rate of turbulence kinetic energy and a turbulence dynamic viscosity. The calculated property values depend on five model constants that are parameters of the sensitivity analysis in this paper. Sensitivity coefficients are derivatives of the above properties, for individual parameters. In this paper these coefficients are determined using a finite difference approximation to the sensitivities coefficients. The author of this paper compares three models of the boundary layer with regard to the sensitivity of properties to the parameters. Irrespective of the boundary layer model used here, the analysis of sensitivity coefficients for the channel flow properties shows that the most sensitive property is the turbulence dissipation rate. Next properties of consequence, although of significantly smaller values of sensitivity coefficients, are the turbulence viscosity and the turbulence kinetic energy. All flow properties are mostly sensitive to the C_µ parameter. One of the final conclusions in this paper is that the analysis of sensitivity coefficient fields allows the reliable checking of results and indicates those areas most prone to calculation difficulties. Copyright © 2008 John Wiley & Sons, Ltd.     

A Refined Shear Deformation Theory of Bending of Isotropic Plates*

ABSTRACT

A nonlinear, in-plane displacement assumption is proposed, based on an undetermined variation df/dz of transverse shear strains through the plate thickness. A second-order ordinary differential equation for f(z) and two surface conditions, as well as a set of eighth-order partial differential equations and four associated boundary conditions, are derived from the principle of minimum potential energy. Coupling exists between the partial and ordinary differential equations. In the homogeneous solutions for the former, in addition to an interior solution contribution, there exist two edge-zone solution contributions, one of which induces self-equilibrated (in the thickness direction) boundary stresses. Three examples are calculated using the present theory. The last gives the stress couple and maximum-stress concentration factors at the free edge of a circular hole in a large bent plate. Numerical results for the examples are compared with those given by three-dimensional elasticity theory and several two-dimensional theories. It is found that the present theory can accurately predict nonlinear variations of in-plane stresses through the thickness of a plate.     

An unsymmetric constitutive equation for anisotropic viscoelastic fluid

A continuum constitutive theory of corotational derivative type is developed for the anisotropic viscoelastic fluid–liquid crystalline (LC) polymers. A concept of anisotropic viscoelastic simple fluid is introduced. The stress tensor instead of the velocity gradient tensor D in the classic Leslie–Ericksen theory is described by the first Rivlin–Ericksen tensor A and a spin tensor W measured with respect to a co-rotational coordinate system. A model LCP-H on this theory is proposed and the characteristic unsymmetric behaviour of the shear stress is predicted for LC polymer liquids. Two shear stresses thereby in shear flow of LC polymer liquids lead to internal vortex flow and rotational flow. The conclusion could be of theoretical meaning for the modern liquid crystalline display technology. By using the equation, extrusion–extensional flows of the fluid are studied for fiber spinning of LC polymer melts, the elongational viscosity vs. extension rate with variation of shear rate is given in figures. A considerable increase of elongational viscosity and bifurcation behaviour are observed when the orientational motion of the director vector is considered. The contraction of extrudate of LC polymer melts is caused by the high elongational viscosity. For anisotropic viscoelastic fluids, an important advance has been made in the investigation on the constitutive equation on the basis of which a series of new anisotropic non-Newtonian fluid problems can be addressed. The project supported by the National Natural Science Foundation of China (10372100, 19832050) (Key project). The English text was polished by Yunming Chen.     

A comparative study of poiseuille flow of a polar fluid under various boundary conditions with applications to blood flow

In this paper, Poiseuille flow of a polar fluid (model of a red blood cell suspension) under various boundary conditions at the wall, viz., slip or no-slip in the axial velocity and couple stresses zero or non-zero at the boundary, is considered from the point of view of its applications to blood flow. Analytic expressions for axial and rotational velocities, flow rate, effective viscosity and stresses are obtained. The magnitudes of the length ratioL and the coupling number N are determined in accordance with concentration and tube radius (in the existing literature, values ofL andN are chosen arbitrarily). Velocity profiles (both axial and rotational) and the variation of the effective viscosity with concentration, tube radius and for various values of the boundary condition parameters are shown graphically. The analytic results obtained are compared with experimental results (for blood flow). It is found that they are in a reasonably good agreement. The effective viscosity exhibits the Inverse Fahraeus-Lindquist Effect in all the cases (including the slip or no-slip in the velocity fields). A method is given for determining the non-zero couple stress boundary condition for a given concentration. Applications of this theory to blood flow are briefly discussed.     

Viscoelastic stresses in the stagnation flow of a dilute polymer solution

In this paper, we consider viscoelastic stresses T₁₁, T₁₂ and T₂₂ arising in the stagnation flow of a dilute polymer solution; in particular, we consider an upper convected Maxwell (UCM) fluid. We present exact solutions to the coupled partial differential equations describing the viscoelastic stresses and deduce the results for the stress T₂₂ of Becherer et al. [P. Becherer, A.N. Morozov, W. van Saarloos, Scaling of singular structures in extensional flow of dilute polymer solutions, J. Non-Newtonian Fluid Mech. 153 (2008) 183–190]. As we considered the viscoelastic stresses over two spatial variables, we are able to study the effect of variable boundary data at the inflow. As such, our results are applicable to a wider range of fluid flow problems.     

Existence and Stability of Multidimensional Transonic Flows through an Infinite Nozzle of Arbitrary Cross-Sections

We establish the existence and stability of multidimensional steady transonic flows with transonic shocks through an infinite nozzle of arbitrary cross-sections, including a slowly varying de Laval nozzle. The transonic flow is governed by the inviscid potential flow equation with supersonic upstream flow at the entrance, uniform subsonic downstream flow at the exit at infinity, and the slip boundary condition on the nozzle boundary. Our results indicate that, if the supersonic upstream flow at the entrance is sufficiently close to a uniform flow, there exists a solution that consists of a C ^1,α subsonic flow in the unbounded downstream region, converging to a uniform velocity state at infinity, and a C ^1,α multidimensional transonic shock separating the subsonic flow from the supersonic upstream flow; the uniform velocity state at the exit at infinity in the downstream direction is uniquely determined by the supersonic upstream flow; and the shock is orthogonal to the nozzle boundary at every point of their intersection. In order to construct such a transonic flow, we reformulate the multidimensional transonic nozzle problem into a free boundary problem for the subsonic phase, in which the equation is elliptic and the free boundary is a transonic shock. The free boundary conditions are determined by the Rankine–Hugoniot conditions along the shock. We further develop a nonlinear iteration approach and employ its advantages to deal with such a free boundary problem in the unbounded domain. We also prove that the transonic flow with a transonic shock is unique and stable with respect to the nozzle boundary and the smooth supersonic upstream flow at the entrance.