Search Results

Now showing 1 - 10 of 14
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • Item
    A Stokes-consistent backflow stabilization for physiological flows
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2014) Bertoglio, Cristobal; Caiazzo, Alfonso
    In computational fluid dynamics incoming flow at open boundaries, or emphbackflow, often yields to unphysical instabilities for high Reynolds numbers. It is widely accepted that this is due to the incoming energy arising from the convection term, which cannot be empha priori controlled when the velocity field is unknown at the boundary. In order to improve the robustness of the numerical simulations, we propose a stabilized formulation based on a penalization of the residual of a weak Stokes problem on the open boundary, whose viscous part controls the incoming convective energy, while the inertial term contributes to the kinetic energy. We also present different strategies for the approximation of the boundary pressure gradient, which is needed for defining the stabilization term. The method has the advantage that it does not require neither artificial modifications or extensions of the computational domain. Moreover, it is consistent with the Womersley solution. We illustrate our approach on numerical examples