Filters

2010 - 201715

More filters

2019 - 2019142020 - 20221
Reset filters

Settings

End date: 2019 × Subject: Navier-Stokes equations ×

Search Results

Now showing 1 - 10 of 15
  • Thumbnail Image
    Item
    Strong solutions to nonlocal 2D Cahn-Hilliard-Navier-Stokes systems with nonconstant viscosity, degenerate mobility and singular potential
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2016) Frigeri, Sergio; Gal, Ciprian G.; Grasselli, Maurizio; Sprekels, Jürgen
    We consider a nonlinear system which consists of the incompressible Navier-Stokes equations coupled with a convective nonlocal Cahn-Hilliard equation. This is a diffuse interface model which describes the motion of an incompressible isothermal mixture of two (partially) immiscible fluids having the same density. We suppose that the viscosity depends smoothly on the order parameter as well as the mobility. Moreover, we assume that the mobility is degenerate at the pure phases and that the potential is singular (e.g. of logarithmic type). This system is endowed with no-slip boundary condition for the (average) velocity and homogeneous Neumann boundary condition for the chemical potential. Thus the total mass is conserved. In the two-dimensional case, this problem was already analyzed in some joint papers of the first three authors. However, in the present general case, only the existence of a global weak solution, the (conditional) weak-strong uniqueness and the existence of the global attractor were proven. Here we are able to establish the existence of a (unique) strong solution through an approximation procedure based on time discretization. As a consequence, we can prove suitable uniform estimates which allow us to show some smoothness of the global attractor. Finally, we discuss the existence of strong solutions for the convective nonlocal Cahn-Hilliard equation, with a given velocity field, in the three dimensional case as well.
  • Thumbnail Image
    Item
    On a reduced sparsity stabilization of grad-div type for incompressible flow problems
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2013) Linke, Alexander; Rebholz, Leo
    We introduce a new operator for stabilizing error that arises from the weak enforcement of mass conservation in finite element simulations of incompressible flow problems. We show this new operator has a similar positive effect on velocity error as the well-known and very successful grad-div stabilization operator, but the new operator is more attractive from an implementation standpoint because it yields a sparser block structure matrix. That is, while grad-div produces fully coupled block matrices (i.e. block-full), the matrices arising from the new operator are block-upper triangular in two dimensions, and in three dimensions the 2,1 and 3,1 blocks are empty. Moreover, the diagonal blocks of the new operators matrices are identical to those of grad-div. We provide error estimates and numerical examples for finite element simulations with the new operator, which reveals the significant improvement in accuracy it can provide. Solutions found using the new operator are also compared to those using usual grad-div stabilization, and in all cases, solutions are found to be very similar.
  • Thumbnail Image
    Item
    Development of a numerical workflow based on μ-CT imaging for the determination of capillary pressure–saturation-specific interfacial area relationship in 2-phase flow pore-scale porous-media systems: a case study on Heletz sandstone
    (Göttingen : Copernicus Publ., 2016) Peche, Aaron; Halisch, Matthias; Bogdan Tatomir, Alexandru; Sauter, Martin
    In this case study, we present the implementation of a finite element method (FEM)-based numerical pore-scale model that is able to track and quantify the propagating fluid–fluid interfacial area on highly complex micro-computed tomography (μ-CT)-obtained geometries. Special focus is drawn to the relationship between reservoir-specific capillary pressure (pc), wetting phase saturation (Sw) and interfacial area (awn). The basis of this approach is high-resolution μ-CT images representing the geometrical characteristics of a georeservoir sample. The successfully validated 2-phase flow model is based on the Navier–Stokes equations, including the surface tension force, in order to consider capillary effects for the computation of flow and the phase-field method for the emulation of a sharp fluid–fluid interface. In combination with specialized software packages, a complex high-resolution modelling domain can be obtained. A numerical workflow based on representative elementary volume (REV)-scale pore-size distributions is introduced. This workflow aims at the successive modification of model and model set-up for simulating, such as a type of 2-phase problem on asymmetric μ-CT-based model domains. The geometrical complexity is gradually increased, starting from idealized pore geometries until complex μ-CT-based pore network domains, whereas all domains represent geostatistics of the REV-scale core sample pore-size distribution. Finally, the model can be applied to a complex μ-CT-based model domain and the pc–Sw–awn relationship can be computed.
  • Thumbnail Image
    Item
    Analysis of improved Nernst-Planck-Poisson models of compressible isothermal electrolytes. Part III: Compactness and convergence
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2017) Dreyer, Wolfgang; Druet, Pierre-Étienne; Gajewski, Paul; Guhlke, Clemens
    We consider an improved Nernst-Planck-Poisson model first proposed by Dreyer et al. in 2013 for compressible isothermal electrolytes in non equilibrium. The model takes into account the elastic deformation of the medium that induces an inherent coupling of mass and momentum transport. The model consists of convectiondiffusionreaction equations for the constituents of the mixture, of the Navier-Stokes equation for the barycentric velocity, and of the Poisson equation for the electrical potential. Due to the principle of mass conservation, crossdiffusion phenomena must occur and the mobility matrix (Onsager matrix) has a kernel. In this paper, which continues the investigations of [DDGG17a, DDGG17b], we prove the compactness of the solution vector, and existence and convergence for the approximation schemes. We point at simple structural PDE arguments as an adequate substitute to the AubinLions compactness Lemma and its generalisations: These familiar techniques attain their limit in the context of our model in which the relationship between time derivatives (transport) and diffusion gradients is highly non linear.
  • Thumbnail Image
    Item
    Analysis of improved Nernst-Planck-Poisson models of compressible isothermal electrolytes. Part II: Approximation and a priori estimates
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2017) Dreyer, Wolfgang; Druet, Pierre-Étienne; Gajewski, Paul; Guhlke, Clemens
    We consider an improved NernstPlanckPoisson model first proposed by Dreyer et al. in 2013 for compressible isothermal electrolytes in non equilibrium. The model takes into account the elastic deformation of the medium that induces an inherent coupling of mass and momentum transport. The model consists of convectiondiffusionreaction equations for the constituents of the mixture, of the Navier-Stokes equation for the barycentric velocity, and of the Poisson equation for the electrical potential. Due to the principle of mass conservation, crossdiffusion phenomena must occur and the mobility matrix (Onsager matrix) has a kernel. In this paper, which continues the investigation of [DDGG17a], we derive for thermodynamically consistent approximation schemes the natural uniform estimates associated with the dissipations. Our results essentially improve our former study [DDGG16], in particular the a priori estimates concerning the relative chemical potentials.
  • Thumbnail Image
    Item
    Analysis of improved Nernst-Planck-Poisson models of isothermal compressible electrolytes subject to chemical reactions: The case of a degenerate mobility matrix
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2016) Druet, Pierre-Etienne
    We continue our investigations of the improved NernstPlanckPoisson model introduced in [DGM13]. In the paper [DDGG16] the analysis relies on the hypothesis that the mobility matrix has maximal rank under the constraint of mass conservation (rank N-1 for a mixture of N species). In this paper we allow for the case that the positive eigenvalues of the mobility matrix tend to zero along with the partial mass densities of certain species. In this approach the mobility matrix has a variable rank between zero and N-1 according to the number of locally available species. We set up a concept of weak solution able to deal with this scenario, showing in particular how to extend the fundamental notion of differences of chemical potentials that supports the modelling and the analysis in [DDGG16]. We prove the global-in-time existence in this solution class.
  • Thumbnail Image
    Item
    Analysis of improved Nernst-Planck-Poisson models of compressible isothermal electrolytes. Part I: Derivation of the model and survey of the results
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2017) Dreyer, Wolfgang; Druet, Pierre-Étienne; Gajewski, Paul; Guhlke, Clemens
    We consider an improved NernstPlanckPoisson model first proposed by Dreyer et al. in 2013 for compressible isothermal electrolytes in non equilibrium. The model takes into account the elastic deformation of the medium that induces an inherent coupling of mass and momentum transport. The model consists of convectiondiffusionreaction equations for the constituents of the mixture, of the Navier-Stokes equation for the barycentric velocity, and of the Poisson equation for the electrical potential. Due to the principle of mass conservation, crossdiffusion phenomena must occur and the mobility matrix (Onsager matrix) has a kernel. In this paper we establish the existence of a globalintime weak solution for the full model, allowing for a general structure of the mobility tensor and for chemical reactions with highly non linear rates in the bulk and on the active boundary. We characterise the singular states of the system, showing that the chemical species can vanish only globally in space, and that this phenomenon must be concentrated in a compact set of measure zero in time. With respect to our former study [DDGG16], we also essentially improve the a priori estimates, in particular concerning the relative chemical potentials.
  • Thumbnail Image
    Item
    Existence of weak solutions for improved Nernst-Planck-Poisson models of compressible reacting electrolytes
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2016) Dreyer, Wolfgang; Druet, Pierre-Étienne; Gajewski, Paul; Guhlke, Clemens
    We consider an improved Nernst-Planck-Poisson model for compressible electrolytes first proposed by Dreyer et al. in 2013. The model takes into account the elastic deformation of the medium. In particular, large pressure contributions near electrochemical interfaces induce an inherent coupling of mass and momentum transport. The model consists of convection-diffusion-reaction equations for the constituents of the mixture, of the Navier-Stokes equation for the barycentric velocity and the Poisson equation for the electrical potential. Cross-diffusion phenomena occur due to the principle of mass conservation. Moreover, the diffusion matrix (mobility matrix) has a zero eigenvalue, meaning that the system is degenerate parabolic. In this paper we establish the existence of a global-in-time weak solution for the full model, allowing for cross-diffusion and an arbitrary number of chemical reactions in the bulk and on the active boundary.
  • Thumbnail Image
    Item
    On nonlocal Cahn-Hilliard-Navier-Stokes systems in two dimensions
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2014) Frigeri, Sergio; Gal, Cipian G.; Grasselli, Maurizio
    We consider a diffuse interface model which describes the motion of an incompressible isothermal mixture of two immiscible fluids. This model consists of the Navier-Stokes equations coupled with a convective nonlocal Cahn-Hilliard equation. Several results were already proven by two of the present authors. However, in the two-dimensional case, the uniqueness of weak solutions was still open. Here we establish such a result even in the case of degenerate mobility and singular potential. Moreover, we show the strong-weak uniqueness in the case of viscosity depending on the order parameter, provided that the mobility is constant and the potential is regular. In the case of constant viscosity, on account of the uniqueness results we can deduce the connectedness of the global attractor whose existence was obtained in a previous paper. The uniqueness technique can be adapted to show the validity of a smoothing property for the difference of two trajectories which is crucial to establish the existence of an exponential attractor.
  • Thumbnail Image
    Item
    On the convergence rate of grad-div stabilized Taylor-Hood to Scott-Vogelius solutions for incompressible flow problems
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2011) Linke, Alexander; Rebholz, Leo G.; Wilson, Nicholas E.
    It was recently proven that, under mild restrictions, grad-div stabilized Taylor-Hood solutions of Navier-Stokes problems converge to the Scott-Vogelius solution of that same problem. However, even though the analytical rate was only shown to be gamma^-frac12 (where gamma is the stabilization parameter), the computational results suggest the rate may be improvable gamma^-1. We prove herein the analytical rate is indeed gamma^-1, and extend the result to other incompressible flow problems including Leray-alpha and MHD. Numerical results are given that verify the theory.