Relative energy approach to a diffuse interface model of a compressible two‐phase flow

We propose a simple model for a two‐phase flow with a diffuse interface. The model couples the compressible Navier‐Stokes system governing the evolution of the fluid density and the velocity field with the Allen‐Cahn equation for the order parameter. We show that the model is thermodynamically consistent, in particular, a variant of the relative energy inequality holds. As a corollary, we show the weak‐strong uniqueness principle, meaning any weak solution coincides with the strong solution emanating from the same initial data on the life span of the latter. Such a result plays a crucial role in the analysis of the associated numerical schemes. Finally, we perform the low Mach number limit obtaining the standard incompressible model.     

Development of a compressible multiphase cavitation approach for diesel spray modelling

The influence of in-nozzle phenomena including cavitation on the morphology of the spray from a diesel injector with a sharp nozzle inlet is investigated numerically. A compressible, multi-phase Volume of Fluid Large Eddy Simulation is implemented in the OpenFOAM environment. The volume fraction transport equations for liquid and gas phases are reformulated to include mass transfer source terms. These source terms are modelled with two cavitation models by Schnerr and Kunz, which are extended to eliminate non-physical mass transfer rates. Validation is carried out only for the Schnerr cavitation model due to its independence of empirical parameters. The numerical method is validated by comparing the simulated mass flow rates, pressure and velocity profiles at different cavitation conditions against published experimental data obtained using a slightly converging square channel. Favourable comparison between simulations and experiments is achieved with minor discrepancies attributable to uncertainties in fuel properties, experimental artefacts and assumptions made in numerical models. Application of the method to calculation of in-nozzle phenomena and primary breakup of a diesel spray reveals that in-nozzle flow separation, wall shear and cavitation contribute greatly to the fragmentation of the jet. Comparison of the two cavitation models shows that after the onset of complete flow detachment, the Kunz implementation predicts higher air inflow at the nozzle outlet than the Schnerr model.     

Model Predictive Control of two-phase flow using a diffuse interface approach

We propose a feedback control strategy for the stabilization of two-phase flow with variable densities, where we use an inexact variant of Model Predictive Control, called Instantaneous Control. The physical system is governed by the Cahn–Hilliard Navier–Stokes system with different densities proposed in [1]. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An evaluation of interface capturing methods in a VOF based model for multiphase flow of a non-Newtonian ceramic in tape casting

The aim of the present study is to evaluate the different interface capturing methods as well as to find the best approach for flow modeling of the ceramic slurry in the tape casting process. The conventional volume of fluid (VOF) method with three different interpolation methods for interface capturing, i.e. the Geometric Reconstruction Scheme (GRS), High Resolution Interface Capturing (HRIC) and Compressive Interface Capturing Scheme for Arbitrary Meshes (CICSAM), are investigated for the advection of the VOF, both for Newtonian and non-Newtonian cases. The main purpose is to find the best method for the free surface capturing during the flow of a ceramic slurry described by a constitutive power law equation in the tape casting process. First the developed model is tested against well-documented and relevant solutions from literature involving free surface tracking and subsequently it is used to investigate the flow of a La_0.85Sr_0.15MnO₃ (LSM) ceramic slurry modeled with the Ostwald de Waele power law. Results of the modeling are compared with corresponding experimental data and good agreement is found.     

A fully Eulerian multiphase model of windblown sand coupled with morphodynamic evolution: Erosion,transport, deposition,and avalanching

Modeling unsteady windblown sand dynamics requires not only treatment of the sand present in the air as a suspended constituent of a mixture but also consideration of erosion and sedimentation phenomena and consequently of the morphodynamic evolution of the sand-bed surface, including avalanching, especially in the presence of natural or human-built obstacles, artifacts, and infrastructures. With this aim in mind, we present a comprehensive multiphase model capable of accurately simulating all the physical phenomena mentioned above, producing satisfactory results, with reasonable computational effort. As test cases, two- and three-dimensional simulations of dune evolution are reported, as is windblown sand transport over a straight vertical wall. Examples of sand transport around other obstacles are given to show the flexibility of the model and its usefulness for such engineering applications.     

A gradient flow approach to an evolution problem arising in superconductivity

We study an evolution equation proposed by Chapman, Rubinstein, and Schatzman as a mean‐field model for the evolution of the vortex density in a superconductor. We treat the case of a bounded domain where vortices can exit or enter the domain. We show that the equation can be derived rigorously as the gradient flow of some specific energy for the Riemannian structure induced by the Wasserstein distance on probability measures. This leads us to some existence and uniqueness results and energy‐dissipation identities. We also exhibit some “entropies” that decrease through the flow and allow us to get regularity results (solutions starting in L^p, p > 1, remain in L^p). © 2007 Wiley Periodicals, Inc.     

Coupling of compressible and incompressible flow regions using the multiple pressure variables approach

In many cases, multiphase flows are simulated on the basis of the incompressible Navier–Stokes equations. This assumption is valid as long as the density changes in the gas phase can be neglected. Yet, for certain technical applications such as fuel injection, this is no longer the case, and at least the gaseous phase has to be treated as a compressible fluid. In this paper, we consider the coupling of a compressible flow region to an incompressible one based on a splitting of the pressure into a thermodynamic and a hydrodynamic part. The compressible Euler equations are then connected to the Mach number zero limit equations in the other region. These limit equations can be solved analytically in one space dimension that allows to couple them to the solution of a half‐Riemann problem on the compressible side with the help of velocity and pressure jump conditions across the interface. At the interface location, the flux terms for the compressible flow solver are provided by the coupling algorithms. The coupling is demonstrated in a one‐dimensional framework by use of a discontinuous Galerkin scheme for compressible two‐phase flow with a sharp interface tracking via a ghost‐fluid type method. The coupling schemes are applied to two generic test cases. The computational results are compared with those obtained with the fully compressible two‐phase flow solver, where the Mach number zero limit is approached by a weakly compressible fluid. For all cases, we obtain a very good agreement between the coupling approaches and the fully compressible solver. Copyright © 2014 John Wiley & Sons, Ltd.     

Properties of a solution to a nonstationary flow of a viscous weakly compressible gas

We study the dependence of integral norms of an exact solution to a nonstationary flow of a viscous weakly compressible gas on a parameter characterizing the compressibility and viscosity in the case of two spatial variables.     

Packing 4-Cycles in Eulerian and Bipartite Eulerian Tournaments with an Application to Distances in Interchange Graphs

We prove that every Eulerian orientation of K_m,n contains arc-disjoint directed 4-cycles, improving earlier lower bounds. Combined with a probabilistic argument, this result is used to prove that every regular tournament with n vertices contains arc-disjoint directed 4-cycles. The result is also used to provide an upper bound for the distance between two antipodal vertices in interchange graphs.Received February 6, 2004     

Asymptotic behaviour of solutions to the Navier-Stokes equations of a two-dimensional compressible flow

In this paper,we are concerned with the asymptotic behaviour of a weak solution to the Navier-Stokes equations for compressible barotropic flow in two space dimensions with the pressure function satisfying p(e) = a log d(e) for large .Here d > 2,a > 0.We introduce useful tools from the theory of Orlicz spaces and construct a suitable function which approximates the density for time going to infinity.Using properties of this function,we can prove the strong convergence of the density to its limit state.The behaviour of the velocity field and kinetic energy is also briefly discussed.     

On the well-posedness of a linearized plasma-vacuum interface problem in ideal compressible MHD

We study the initial-boundary value problem resulting from the linearization of the plasma-vacuum interface problem in ideal compressible magnetohydrodynamics (MHD). We suppose that the plasma and the vacuum regions are unbounded domains and the plasma density does not go to zero continuously, but jumps. For the basic state upon which we perform linearization we find two cases of well-posedness of the “frozen” coefficient problem: the “gas dynamical” case and the “purely MHD” case. In the “gas dynamical” case we assume that the jump of the normal derivative of the total pressure is always negative. In the “purely MHD” case this condition can be violated but the plasma and the vacuum magnetic fields are assumed to be non-zero and non-parallel to each other everywhere on the interface. For this case we prove a basic a priori estimate in the anisotropic weighted Sobolev space for the variable coefficient problem.     

Simulation of a compositional model for multiphase flow in porous media

In this article, we consider the simulation of a compositional model for three‐dimensional, three‐phase, multicomponent flow in a porous medium. This model consists of Darcy's law for volumetric flow velocities, mass conservation for hydrocarbon components, thermodynamic equilibrium for mass interchange between phases, and an equation of state for saturations. A discretization scheme based on the block‐centered finite difference method for pressures and compositions is developed. Numerical results are reported for the benchmark problem of the third comparative solution project (CSP) organized by the society of petroleum engineers (SPE). © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005.     

Atomic decomposition for the vorticity of a viscous flow in the whole space

We show that the vorticity of a viscous flow in ?³ admits an atomic decomposition of the form ω(x, t) = ω_k(x – x_k, t), with localized and oscillating building blocks ω_k, if such a property is satisfied at the beginning of the evolution. We also study the long time behavior of an isolated coherent structure and the special behavior of flows with highly oscillating vorticities. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Integrodifferential approaches to frequency analysis and control design for compressible fluid flow in a pipeline element

In this study, modelling, frequency analysis, and optimization of control processes are considered for the fluid flow in pipeline systems. A mathematical model of controlled pipeline elements with distributed parameters is proposed to describe the dynamical behaviour of compressible fluid which is transported in a long rigid tube. By exploiting specific functions representing cross-sectional forces and effective displacements as well as linear approximations of fluidic resistances, the original problem with non-uniform parameters is reduced to a partial differential equation (PDE) system with constant coefficients and homogeneous initial and boundary conditions. Three numerical approaches are applied to an efficient analysis of natural vibrations and reliable control-oriented modelling of pipeline elements. The conventional Galerkin method is compared with the method of integrodifferential relations based on a weak formulation of the constitutive laws. In the latter approach, the original initial-boundary value problem is reduced to the minimization of an error functional which provides explicit energy estimates of the solution quality. A novel projection approach is implemented on the basis of the Petrov–Galerkin method combined with the method of integrodifferential relations. This technique benefits from the advantages of the above-mentioned projection and variational approaches, namely sufficient numerical stability, a lower differential order, and an explicit quality estimation. Numerical optimization procedures, making use of a modified finite element technique, are proposed to obtain a feedforward control strategy for changing the pressure and mass flow inside the pipeline system to a desired operating state. At this given finite point of time, residual elastic oscillations inside the pipeline are minimized. Numerical results, obtained for ideal as well as viscous fluid models, are analysed and discussed.     

The numerical simulation of multiphase flow through a porous medium and its application to reservoir engineering

The analysis of multiphase flow in porous media is of considerable significance in the field of petroleum reservoir simulation, where accurate predictions of fluid flow are important in assessing the performance of oil and gas fields. The specific case of two-phase immiscible flow is considered by first deriving the governing nonlinear partial differential equations. The space discretization is then carried out making use of the additional versatility of the finite element method compared with originally used finite differences. By using a pair of dependent variables P and R, the bandwidth of the discrete space-continuous time equations may be reduced to increase significantly the speed of the algorithm. A discussion of time-stepping methods is followed by an application of the technique to a five-spot,an extraction pattern used in the field. The boundary conditions used to simulate this flow pattern are also discussed.     

A numerical model of the reacting multiphase flow in a pulp digester

This paper summarises the development, implementation and typical results from a new full-scale three-dimensional numerical model of the multiphase chemically reacting flow in a continuous digester. The model is based on continuity and momentum equations for wood chips and free liquid. A new sub-model describing the tangential stresses for the solid phase has been developed and incorporated into the digester model. The solution procedure that takes into account the high inter-phase friction terms and the high values of the solid pressure has been introduced. In order to illustrate the model’s performance, simulations have been carried on for the industrial digester. The results are in quantitative agreement with available field measurements and in qualitative agreement with operating personnel observations. The model utilises curvilinear body fitted coordinates and can be used to simulate the operation of various pulp digesters or any other chemical reactors working in a similar regime.     

The large deformation of nonlinearly elastic rings in a two-dimensional compressible flow

The nonlinear nonlocal system of the equilibrium equations ofan elastic ring under the action of an external two-dimensionaluniformly subsonic potential barotropic steady-state gas flowis considered. The configurations of the elastic ring are identifiedby a pair of functions (, ). The simple curve represents theshape of the ring and the real-valued function identifies theorientation of the material sections of the ring. The pressurefield on the ring depends nonlocally on , and on two parametersU and P which represent the pressure and the velocity at infinity.The system is shown to be equivalent to a fixed-point problem,which is then treated with continuation methods. It is shownthat the solution branch ensuing from certain equilibrium states((₀, ₀), 0, P₀) in the solution-parameter space of ((₀, ₀),0, P₀) either approaches the boundary of the admissible ((,), U,p)'s in a well-defined sense, or is unbounded, or is homotopicallynontrivial in the sense that there exists a continuous map from the branch to a two-dimensional sphere which is not homotopicin the sphere to a constant, while restricted to the branchminus ((₀, ₀), 0, P₀) is homotopic to a constant in the sphere.Furthermore, by fixing the pressure parameter at P₀ and by consideringthe one-parameter problem in ((, ), U), the following holds.Every hyperplane in the solution-parameter space of the ((,), U)'s which contains the equilibrium state ((₀, ₀), 0) anddoes not include a welldetermined one-dimensional subspace intersectsthe solution branch above at a point different from ((₀, ₀),0).     

Turbulent compressible fluid flow computation in a male rotor-housinggap of screw compressors

This contribution is devoted to a turbulent compressible fluid flow computation through a 2-D model of the male rotor-housing gap in screw compressors. For the numerical solution of the nonlinear conservative system of the Favre-averaged Navier-Stokes equations the cell-centred finite volume formulation of the explicit two-step TVD MacCormack scheme proposed by Causon on a structured quadrilateral grid is used. The turbulent viscosity is calculated by using the algebraic Baldwin-Lomax turbulence model which is implemented into the own developed numerical code. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)