Global existence of weak solutions for a viscous two-phase model

The purpose of this paper is to explore a viscous two-phase liquid-gas model relevant for well and pipe flow. Our approach relies on applying suitable modifications of techniques previously used for studying the single-phase isothermal Navier-Stokes equations. A main issue is the introduction of a novel two-phase variant of the potential energy function needed for obtaining fundamental a priori estimates. We derive an existence result for weak solutions in a setting where transition to single-phase flow is guaranteed not to occur when the initial state is a true mixture of both phases. Some numerical examples are also included in order to demonstrate characteristic behavior of solutions. In particular, we illustrate how two-phase flow is genuinely different compared to single-phase flow concerning the behavior of an initial mass discontinuity.     

A blow-up criterion for a 2D viscous liquid-gas two-phase flow model

In this paper, we prove a blow-up criterion in terms of the upper bound of the liquid mass for the strong solution to the two-dimensional (2D) viscous liquid-gas two-phase flow model in a smooth bounded domain. The result also applies to three-dimensional (3D) case.     

Global weak solutions for a viscous liquid–gas model with transition to single-phase gas flow and vacuum

This work deals with a viscous two-phase liquid–gas model relevant to the flow in wells and pipelines. The liquid is treated as an incompressible fluid whereas the gas is assumed to be polytropic. The model is rewritten in terms of Lagrangian coordinates and is studied in a free boundary setting where the liquid and gas masses are of compact support initially, and continuous at the boundary. Consequently, the initial masses involve a transition to single-phase gas flow and vacuum at the boundary. An appropriate balance between pressure and viscous forces is identified which allows obtaining pointwise upper and lower estimates of masses. These estimates rely on the assumption of a certain relation between the rate of degeneracy of the viscosity coefficient and the rate that determines how fast the initial masses are vanishing at the boundary. By combining these estimates with basic energy type of estimates, higher order regularity estimates are obtained. The existence of global weak solutions is then proved by showing compactness for a class of semi-discrete approximations.     

Free boundary value problem for a viscous two-phase model with mass-dependent viscosity

In this paper, we study a free boundary value problem for two-phase liquid-gas model with mass-dependent viscosity coefficient when both the initial liquid and gas masses connect to vacuum with a discontinuity. This is an extension of the paper [S. Evje, K.H. Karlsen, Global weak solutions for a viscous liquid-gas model with singular pressure law, http://www.irisresearch.no/docsent/emp.nsf/wvAnsatte/SEV]. Just as in [S. Evje, K.H. Karlsen, Global weak solutions for a viscous liquid-gas model with singular pressure law, http://www.irisresearch.no/docsent/emp.nsf/wvAnsatte/SEV], the gas is assumed to be polytropic whereas the liquid is treated as an incompressible fluid. We give the proof of the global existence and uniqueness of weak solutions when β∈(0,1], which have improved the previous result of Evje and Karlsen, and get the asymptotic behavior result, also we obtain the regularity of the solutions by energy method.     

Existence of weak solutions for a non-classical sharp interface model for a two-phase flow of viscous,incompressible fluids

We introduce a new sharp interface model for the flow of two immiscible, viscous, incompressible fluids. In contrast to classical models for two-phase flows we prescribe an evolution law for the interfaces that takes diffusional effects into account. This leads to a coupled system of Navier–Stokes and Mullins–Sekerka type parts that coincides with the asymptotic limit of a diffuse interface model. We prove the long-time existence of weak solutions, which is an open problem for the classical two-phase model. We show that the phase interfaces have in almost all points a generalized mean curvature.     

Application of finite element method in aeroelasticity

The subject of this paper is the numerical simulation of the interaction of two-dimensional incompressible viscous flow and a vibrating airfoil, which can rotate around the elastic axis and oscillate in the vertical direction. The numerical simulation consists of the finite element approximation of the Navier–Stokes equations coupled with the system of ordinary differential equations describing the airfoil motion. The arbitrary Lagrangian–Eulerian (ALE) formulation of the Navier–Stokes equations, stabilization the finite element discretization and coupling of both models is discussed. Moreover, the Reynolds averaged Navier–Stokes (RANS) system of equations together with the Spallart–Almaras turbulence model is also discussed. The computational results of aeroelastic calculations are presented and compared with the NASTRAN code solutions.     

Numerical simulation of the performance of a snow fence with airfoil snow plates by FVM

The aim of this work is to analyze the efficiency of a snow fence with airfoil snow plates to avoid the snowdrift formation, to improve visibility and to prevent blowing snow disasters on highways and railways. In order to attain this objective, it is necessary to solve particle transport equations along with the turbulent fluid flow equations since there are two phases: solid phase (snow particles) and fluid phase (air). In the first place, the turbulent flow is modelled by solving the Reynolds-averaged Navier-Stokes (RANS) equations for incompressible viscous flows through the finite volume method (FVM) and then, once the flow velocity field has been determined, representative particles are tracked using the Lagrangian approach. Within the particle transport models, we have used a particle transport model termed as Lagrangian particle tracking model, where particulates are tracked through the flow in a Lagrangian way. The full particulate phase is modelled by just a sample of about 15,000 individual particles. The tracking is carried out by forming a set of ordinary differential equations in time for each particle, consisting of equations for position and velocity. These equations are then integrated using a simple integration method to calculate the behaviour of the particles as they traverse the flow domain. Finally, the conclusions of this work are exposed.     

Global weak solutions for a compressible gas-liquid model with well-formation interaction

The objective of this work is to explore a compressible gas-liquid model designed for modeling of well flow processes. We build into the model well-reservoir interaction by allowing flow of gas between well and formation (surrounding reservoir). Inflow of gas and subsequent expansion of gas as it ascends towards the top of the well (a so-called gas kick) represents a major concern for various well operations in the context of petroleum engineering. We obtain a global existence result under suitable assumptions on the regularity of initial data and the rate function that controls the flow of gas between well and formation. Uniqueness is also obtained by imposing more regularity on the initial data. The key estimates are to obtain appropriate lower and upper bounds on the gas and liquid masses. For that purpose we introduce a transformed version of the original model that is highly convenient for analysis of the original model. In particular, in the analysis of the transformed model additional terms, representing well-formation interaction, can be treated by natural extensions of arguments that previously have been employed for the single-phase Navier-Stokes model. The analysis ensures that transition to single-phase regions do not appear when the initial state is a true gas-liquid mixture.     

Modeling of two-phase flows with surface tension by finite pointset method (FPM)

A meshfree method for two-phase immiscible incompressible flows including surface tension is presented. The continuum surface force (CSF) model is used to include the surface tension force. The incompressible Navier–Stokes equation is considered as the mathematical model. Application of implicit projection method results in linear second-order partial differential equations for velocities and pressure. These equations are then solved by the finite pointset method (FPM), which is a meshfree and Lagrangian method. The fluid is represented as finite number of particles and the immiscible fluids are distinguished by the color of each particle. The interface is tracked automatically by advecting the color functions for each particle. Two test cases, Laplace's law and the Rayleigh–Taylor instability in 2D have been presented. The results are found to be consistent with the theoretical results.     

On subsonic compressible flows in a two-dimensional duct

This paper concerns subsonic flows passing a two-dimensional duct for the steady compressible Euler system. If the Bernoulli constant is uniform in the flow field, the density at the entry and both the pressures at the entrance and the exit are given, we show that the problem is generally ill-posed; but if we give the pressure at the exit with a constant difference, then under the same other conditions as above we establish the existence of subsonic flows.     

Systems of partial differential equations of mixed hyperbolic-parabolic type

We prove the global existence and uniqueness of solutions of certain mixed hyperbolic-parabolic systems of partial differential equations in one space dimension with initial data that is assumed to be pointwise bounded with possibly large oscillation and with small total energy. The systems we consider are general enough to include the Navier-Stokes equations of compressible flow, the equations of compressible MHD, models of chemical combustion, and others. In particular, the application of our results to the MHD system gives an existence result which is new.     

Crocco variable formulation for uniform shear flow over a stretching surface with transpiration: Multiple solutions and stability

The simultaneous effects of transpiration through and tangential movement of a semi-infinite flat plate on the self-similar boundary layer flow driven by uniform shear in the far field is considered. Difficulties with standard shooting techniques are overcome using Crocco variables which also serve to better elucidate the solution structure. The stabilities of dual, triple and even quadruple steady flow solutions encountered in different ranges of plate stretching and wall stress are determined using a linear temporal stability analysis for the self-similar flow.      

Statistics of numerically generated turbulence

This paper considers numerically generated turbulence obtained by integrating the complete time-dependent three-dimensional Navier-Stokes equations. The simulated unidirectional turbulent flow, bounded by two parallel planes, is strongly inhomogeneous in the direction normal to the planes but homogeneous in the parallel directions. The resulting flow field, which is considered a numerical realization of fully developed turbulent channel flow, contains detailed information on spatial coherent flow structures as well as on the time-dependency and statistics of the three-dimensional velocity and pressure fields. Focussing here on the statistics of the numerically generated turbulence, second-moments and higher-moments are presented and compared with the most recent PTV and LDV laboratory measurements. It is concluded that direct numerical simulations are an invaluable approach to turbulence which complements field studies and laboratory investigations. Numerical experiments are now becoming a principal source of detailed and reliable information, which play a key role in the deepening of our understanding of turbulent flow phenomena.     

Mixed convection boundary layer flow over a permeable vertical surface with prescribed wall heat flux

The effects of suction/injection on the laminar mixed convection boundary-layer flow on a vertical wall with a prescribed heat flux are considered. The conditions which allow the equations to be reduced to similarity form are derived and numerical solutions of the resulting equations are obtained for a range of values of the suction/injection and buoyancy parameters. Two specific cases, corresponding to a stagnation point flow and uniform wall heat flux, are treated in detail. Results are presented in terms of the skin friction and wall temperature with a selection of velocity and temperature profiles also being given. Dual solutions are found to exist for assisting flow, these are an addition to what has been reported previously for opposing flows. Solutions for some limiting values of the parameters are also derived.      

Analysis and finite element approximation of a Ladyzhenskaya model for viscous flow in streamfunction form

In this paper we consider a model for the motion of incompressible viscous flows proposed by Ladyzhenskaya. The Ladyzhenskaya model is written in terms of the velocity and pressure while the studied model is written in terms of the streamfunction only. We derived the streamfunction equation of the Ladyzhenskaya model and present a weak formulation and show that this formulation is equivalent to the velocity–pressure formulation. We also present some existence and uniqueness results for the model. Finite element approximation procedures are presented. The discrete problem is proposed to be well posed and stable. Some error estimates are derived. We consider the 2D driven cavity flow problem and provide graphs which illustrate differences between the approximation procedure presented here and the approximation for the streamfunction form of the Navier–Stokes equations. Streamfunction contours are also displayed showing the main features of the flow.     

宾汉流体稠密两相湍流流动中的二阶矩模型

建立的Bingham流体稠密两相流动的二阶矩-颗粒动力论湍流模型（USM-theta模型）既体现了两相的作用,又体现了屈服应力所引起的附加项,并提出了USM-theta模型下考虑浓度修正值影响的两相湍流流动的算法．利用该模型对圆管内Bingham流体的单相湍流流动、稠密液固两相的湍流流动进行了计算,并和五方程湍流模型进行了比较,结果表明该模型的预测效果更好．利用USM-theta模型对含颗粒的Bingham流体的两相湍流流动进行了模拟,随着屈服应力的增加,Bingham流体相与颗粒相在管道中心附近的主流速度减小．液固两相湍流和Bingham流体两相湍流的计算结果表明屈服应力引起的附加项对流动有很重要的影响．     

Low-order spectral models of the atmospheric circulation: A survey

A quasi-geostrophic potential vorticity equation is derived from the Navier-Stokes equations for atmospheric motions. It describes the evolution of a quasi-horizontally flow on time scales of a few days and more. The associated boundary-value problem is analyzed by projection of the equation onto orthonormal eigenfunctions (modes) of a Sturm-Liouville operator. The result is a spectral model, consisting of an infinite number of nonlinear ordinary differential equations for the evolution of the mode amplitudes. Low-order spectral models, in which only a few modes are resolved, appear to have properties which agree with observations of the atmospheric circulation. However, little justification is available for truncating the spectral expansion at low resolution numbers. It is argued that stochastic forcing terms should be added to the equations, but it is not a priori clear how they should be specified.A derivation is presented of a specific low-order spectral model of the quasi-geostrophic potential vorticity equation. Some of its subsystems are analyzed for their physical and mathematical properties. It appears that topography can act as a triggering mechanism to generate multiple equilibria. The corresponding flow patterns resemble preference states of the atmospheric circulation. The systems can vacillate between three characteristic regimes with transitions provided either by external or internal mechanisms. A discussion is presented on the validity of stochastically forced spectral models and deterministic chaotic models for the atmospheric circulation.Present affiliation: Institute of Meteorology and Oceanography, Princetonplein 5, 3584 CC, Utrecht, The Netherlands.     

A blow-up criterion for 3D compressible viscoelastic flow with large initial data

In this paper, we consider a Cauchy problem for the three-dimensional compressible viscoelastic flow with large initial data. We establish a blow-up criterion for the strong solutions in terms of the gradient of velocity only, which is similar to the Beale-Kato-Majda criterion for ideal incompressible flow (cf. Beale et al. (1984) [20]) and the blow-up criterion for the compressible Navier-Stokes equations (cf. Huang et al. (2011) [21]).     

Particle trajectories beneath small amplitude shallow water waves in constant vorticity flows

We investigate the particle trajectories in a constant vorticity shallow water flow over a flat bed as periodic waves propagate on the water’s free surface. Within the framework of small amplitude waves, we find the solutions of the nonlinear differential equations system which describes the particle motion in the considered case, and we describe the possible particle trajectories. Depending on the relation between the initial data and the constant vorticity, some particle trajectories are undulating curves to the right, or to the left, others are loops with forward drift, or with backward drift, others can follow some peculiar shapes.