Navier Stokes Korteweg Research Articles

Time-Implicit High-Order Accurate Positivity-Preserving Discretizations for the Navier–Stokes and Navier–Stokes–Korteweg Equations

Time-Implicit High-Order Accurate Positivity-Preserving Discretizations for the Navier–Stokes and Navier–Stokes–Korteweg Equations

Large-time behavior of the solutions for an inflow problem of the full compressible Navier–Stokes–Korteweg equations

In this article, we are concerned with the large-time behavior of solutions to an inflow problem for the one-dimensional full compressible Navier–Stokes–Korteweg equations, which can be used to model the compressible viscous fluids with internal capillarity and heat conductivity. Under some suitable assumptions of the far fields and boundary values of the density, the velocity and the absolute temperature, the large-time behavior of the solutions for the initial boundary value problem is govern by the stationary solution and basic waves of the corresponding hyperbolic system. Then the asymptotic stability of the boundary layer, and the superposition of the boundary layer and the three-rarefaction wave are shown provided that the initial perturbation and the strength of the rarefaction wave and stationary wave are small. The proof is mainly based on L2-energy method, some time-decay estimates in Lp-norm for the smoothed rarefaction wave and the space variable decay estimates of the boundary layer. This can be viewed as the first result about the nonlinear stability of the combination of two different wave patterns for the inflow problem of the full compressible Navier–Stokes–Korteweg equations.

Stability of viscous contact wave for the full compressible Navier–Stokes–Korteweg equations with large perturbation

This study focuses on the stability of the viscous contact wave for the one-dimensional full Navier–Stokes–Korteweg equations with density-temperature dependent transport coefficients and large perturbation. Our findings demonstrate a Nishida–Smoller type result, indicating that the solution remains stable under large perturbation as long as γ−1 is sufficiently small. Notably, the smallness of the capillary coefficient is unnecessary. We then employ the initial layer analysis technique to investigate the asymptotic behaviour in the Wk,p norm. We show that the capillary term has a smoothing effect, which implies that the strong solution is indeed a smooth one. Our results represent an improvement over those previously reported in Chen and Sheng (2019 Nonlinearity 32 395–444). Furthermore, by applying the method in this study to the isothermal case, we can achieve a better outcome than Germain and LeFloch (2016 Commun. Pure Appl. Math. 69 3–61).

Global dynamics of large solution for the compressible Navier–Stokes–Korteweg equations

Global dynamics of large solution for the compressible Navier–Stokes–Korteweg equations

The existence of R$$ \mathcal{R} $$‐bounded solution operator for Navier–Stokes–Korteweg model with slip boundary conditions in half space

This paper proves the existence of ‐bounded solution operator families of the resolvent problem of Navier–Stokes–Korteweg model in half‐space with slip boundary condition. Especially we investigate the model for arbitrary constant viscosity and capillarity. We employ the ‐bounded solution operators of the model obtained from the whole space cases and partial Fourier transform techniques to analyze the model.

Vanishing capillarity–viscosity limit of the incompressible Navier–Stokes–Korteweg equations with slip boundary condition

In this paper, we investigate the vanishing capillarity–viscosity limit of the incompressible Navier–Stokes–Korteweg (NSK) equations in a three-dimensional horizontally periodic strip domain, in which the velocity of the fluid is supplemented with slip boundary condition and the gradient of density with Dirichlet boundary condition on the boundary. We prove that there exists an positive constant T0 independent on the capillarity and viscosity coefficients, such that the incompressible NSK equations have a unique strong solution on [0,T0] and the solution is uniformly bounded in H3. Based on the uniform estimates, we further give the convergence rate in H1 from the solutions of the incompressible NSK equations to the solution of the inhomogeneous incompressible Euler equations as the capillarity and viscosity coefficients go to zero simultaneously.

Investigation of the effect of curvature on the local mass flux of evaporating droplets using a phase field method

Droplet evaporation, a fundamental process with wide-ranging applications, involves complex physical mechanisms that significantly influence mass transfer rates and droplet lifetimes. Despite its significance, the influence of droplet deformation induced by capillary forces on local saturation conditions and its subsequent impact on phase change have received limited attention. This study delves into the impact of capillary-induced droplet deformation on local saturation conditions and its subsequent effect on phase change within a quiescent vapor environment. Contrary to most analytical and numerical studies that impose temperature and heat flux conditions at the vapor–liquid interface, our investigation employs direct numerical simulations based on the thermal Navier–Stokes–Korteweg model, capturing the dynamics of the droplet and its interface without relying on ad-hoc interface conditions. Our findings demonstrate that droplet deformation induces free oscillatory motion, which, although not significantly altering droplet lifetime due to viscous damping, leads to non-uniform mass fluxes along the droplet interface. We observe that areas of higher curvature are associated with increased local pressure and reduced temperatures, a discrepancy primarily attributed to evaporative cooling. This effect is further elucidated through an analysis of pressure and energy fluxes, revealing that regions of high curvature experience significant heat consumption, driving evaporative cooling and creating temperature gradients along the interface. Conversely, areas with negative curvature exhibit a reverse trend, indicating the complex interplay of capillarity, deformation, and phase change in droplet evaporation. Our results highlight the intricate relationship between droplet deformation, local pressure, mass flux, and temperature, underscoring the necessity for improved modeling of these factors in future studies on droplet evaporation dynamics.

Interfacial vorticity dynamics for Navier–Stokes–Korteweg system: General theory and application to two-dimensional near-wall cavitation bubble

A theory of interfacial vorticity dynamics is developed for the single-component liquid–vapor flows described by the Navier–Stokes–Korteweg equations. The non-ideal equation of state and surface capillarity are collectively considered via the chemical potential gradient force. Explicit formulas of the interfacial enstrophy and chemical potential fluxes are derived to elucidate their generation mechanisms from the interface, which are expressed by using the fundamental surface physical quantities. Then, the transport equation of the Lamb dilatation (namely, the Lamb vector divergence) is rigorously derived, generalizing the previous results of the Lamb-vector dynamics for single-phase incompressible viscous flows. The derived results are evaluated in a near-wall cavitation bubble combined with the van der Waals equation of state. The roles of surface shear stress and vorticity, surface chemical potential gradient, surface acceleration, surface stretching-shrinking motion and dilatational waves are emphasized in generating the interfacial enstrophy and chemical potential fluxes, featured by the strong vorticity concentration, interface deformation and inward-propagating microjet due to the longitudinal pressure difference. These distinct dynamical features observed near the bubble interface are well captured by the spatial distribution of the Lamb dilatation. Analyzing the dynamics of the Lamb dilatation based on the transport equation demonstrates the significance of the high-curvature rotating surface region and its vicinity in view of their dominance to the unsteady diffusion effect of the Lamb dilatation. The present exposition could provide new physical insights into liquid–vapor flows from the perspective of interfacial vorticity dynamics.

On Rayleigh–Taylor instability in Navier–Stokes–Korteweg equations

This paper focuses on the Rayleigh–Taylor instability in the two-dimensional system of equations of nonhomogeneous incompressible viscous fluids with capillarity effects in a horizontal periodic domain with infinite height. First, we use the modified variational method to construct (linear) unstable solutions for the linearized capillary Rayleigh–Taylor problem. Then, motivated by the Grenier’s idea in (Grenier in Commun. Pure Appl. Math. 53(9):1067–1091, 2000), we further construct approximate solutions with higher-order growing modes to the capillary Rayleigh–Taylor problem and derive the error estimates between both the approximate solutions and nonlinear solutions of the capillary Rayleigh–Taylor problem. Finally, we prove the existence of escape points based on the bootstrap instability method of Hwang–Guo in (Hwang and Guo in Arch. Ration. Mech. Anal. 167(3):235–253, 2003), and thus obtain the nonlinear Rayleigh–Taylor instability result. Our instability result presents that the Rayleigh–Taylor instability can occur in the fluids with capillarity effects for any capillary coefficient kappa >0 if the critical capillary coefficient is infinite. In particular, it improves the previous Zhang’s result in (Zhang in J. Math. Fluid Mech. 24(3):70–23, 2022) with the assumption of smallness of the capillary coefficient.

Energy equality for weak solutions of the compressible Navier–Stokes–Korteweg equations with general pressure law in a bounded domain

In this paper, we study the energy equality for weak solutions to the compressible Navier–Stokes–Korteweg equations with general pressure law in a bounded domain Ω⊂Rd (d=2,3). By exploring the special structure of the nonlinear terms, and using coarea formula, we formulate the sufficient conditions for the regularity of weak solutions to ensure that energy equality holds.

Asymptotic stability of a nonlinear wave for the compressible Navier–Stokes–Korteweg equations in the half space

Asymptotic stability of a nonlinear wave for the compressible Navier–Stokes–Korteweg equations in the half space

Stabilized formulation for phase-transforming flows with special emphasis on cavitation inception

Despite significant advancement in the field of cavitation, the inception of vapor bubbles from pure liquid is still not well understood. In this paper, we propose a numerical algorithm to solve the Navier–Stokes–Korteweg equations, which are phase-field equations that govern the dynamics of multiphase compressible flows. These equations can predict cavitation inception without phenomenological assumptions about mass transfer. We present a modification to the bulk free energy of the fluid that can help us simulate cavitation at much larger length scales than previously possible. We propose a numerical scheme based on a Taylor–Galerkin discretization, a residual-based discontinuity capturing operator and a spatial filtering scheme. Our numerical method can be used for high speed flows with large pressure gradients, which are distinctive features of flows where we observe cavitation inception. We present numerical examples that show the accuracy, stability and robustness of our scheme. Finally, we use the proposed numerical scheme to study benchmark problems in the field of cavitation inception such as flow past a wedge and flow past a bluff body. The proposed numerical scheme opens opportunities to understand different aspects of cavitation in general and inception in particular, in much more detail.

A Structure-Preserving Finite Volume Scheme for a Hyperbolic Reformulation of the Navier–Stokes–Korteweg Equations

In this paper, we present a new explicit second-order accurate structure-preserving finite volume scheme for the first-order hyperbolic reformulation of the Navier–Stokes–Korteweg equations. The model combines the unified Godunov-Peshkov-Romenski model of continuum mechanics with a recently proposed hyperbolic reformulation of the Euler–Korteweg system. The considered PDE system includes an evolution equation for a gradient field that is by construction endowed with a curl-free constraint. The new numerical scheme presented here relies on the use of vertex-based staggered grids and is proven to preserve the curl constraint exactly at the discrete level, up to machine precision. Besides a theoretical proof, we also show evidence of this property via a set of numerical tests, including a stationary droplet, non-condensing bubbles as well as non-stationary Ostwald ripening test cases with several bubbles. We present quantitative and qualitative comparisons of the numerical solution, both, when the new structure-preserving discretization is applied and when it is not. In particular for under-resolved simulations on coarse grids we show that some numerical solutions tend to blow up when the curl-free constraint is not respected.

Zero‐viscosity‐capillarity limit towards rarefaction wave for the full Navier–Stokes–Korteweg system of compressible fluids

We shall consider the one‐dimensional full Navier–Stokes–Korteweg (NSK) system which are used to model full compressible fluids with internal capillarity. Formally, the full NSK system converges, as the viscosity, heat‐conductivity, and capillary coefficient vanish, to the corresponding full Euler equation, and we do justify this for the case that the full Euler equation has a rarefaction wave. To prove this result, we first decompose the solution as a small perturbation of the smooth rarefaction wave, then make use of the rescaling technique and an elementary energy method on two time scales. Moreover, the uniform convergence rates are also obtained.

High friction limits of Euler–Navier–Stokes–Korteweg equations for multicomponent models

High friction limits of Euler–Navier–Stokes–Korteweg equations for multicomponent models

Analytic regularity for Navier–Stokes–Korteweg model on pseudo-measure spaces

Analytic regularity for Navier–Stokes–Korteweg model on pseudo-measure spaces

A relaxation model for the non-isothermal Navier-Stokes-Korteweg equations in confined domains

The Navier-Stokes-Korteweg (NSK) system is a classical diffuse interface model which is based on van der Waals' theory of capillarity. Diffuse interface methods have gained much interest to model two-phase flow in porous media. However, for the numerical solution of the NSK equations two major challenges have to be faced. First, an extended numerical stencil is required due to a third-order term in the linear momentum and the total energy equations. In addition, the dispersive contribution in the linear momentum equations prevents the straightforward use of contact angle boundary conditions. Secondly, any real gas equation of state is based on a non-convex Helmholtz free energy potential which may cause the eigenvalues of the Jacobian of the first-order fluxes to become imaginary numbers inside the spinodal region.In this work, a thermodynamically consistent relaxation model is presented which is used to approximate the NSK equations. The model is complimented by thermodynamically consistent non-equilibrium boundary conditions which take contact angle effects into account. Due to the relaxation approach, the contribution of the Korteweg tensor in the linear momentum and total energy equations can be reduced to first-order terms which enables a straightforward implementation of contact angle boundary conditions in a numerical scheme. Moreover, the definition of a modified pressure function enables to formulate first-order fluxes which remain strictly hyperbolic in the entire spinodal region. The present work is a generalization of a previously presented parabolic relaxation model for the isothermal NSK equations.A high-order discontinuous Galerkin spectral element method which supports curved elements and hanging nodes is employed to discretize the system. The relaxation model and its corresponding boundary conditions are validated using solutions of the original NSK model and analytical results for one-, two- and three-dimensional test cases. The simulation of a spinodal decomposition in a three-dimensional porous structure underlines the capability of the presented approach.

Existence of global strong solution for the Navier–Stokes–Korteweg system in one dimension for strongly degenerate viscosity coefficients

Existence of global strong solution for the Navier–Stokes–Korteweg system in one dimension for strongly degenerate viscosity coefficients

Towards a consistent lattice Boltzmann model for two-phase fluids

We revisit the construction of discrete kinetic models for single-component isothermal two-phase flows. Starting from a kinetic model for a non-ideal fluid, we show that, under conventional scaling, the Navier–Stokes equations with a non-ideal equation of state are recovered in the hydrodynamic limit. A scaling based on the smallness of velocity increments is then introduced, which recovers the full Navier–Stokes–Korteweg equations. The proposed model is realized on a standard lattice and validated on a variety of benchmarks. Through a detailed study of thermodynamic properties including co-existence densities, surface tension, Tolman length and sound speed, we show thermodynamic consistency, well-posedness and convergence of the proposed model. Furthermore, hydrodynamic consistency is demonstrated by verification of Galilean invariance of the dissipation rate of shear and normal modes and the study of visco-capillary coupling effects. Finally, the model is validated on dynamic test cases in three dimensions with complex geometries and large density ratios such as drop impact on textured surfaces and mercury drops coalescence.

Boundary vorticity dynamics of two-phase viscous flow

From the Navier–Stokes–Korteweg equations, the exact relations between the fundamental surface physical quantities for the two-phase viscous flow with the diffuse interface are derived, including density gradient, shear stress, vorticity, pressure, enstrophy flux, and surface curvature. These theoretical results provide a solid foundation of the boundary/interfacial vorticity dynamics and a new tool for the analysis of complex interfacial phenomena in two-phase viscous flows. To demonstrate the application of the developed results, simulation of a droplet impacting and spreading on a solid wall is conducted by using a recently developed well-balanced discrete unified gas kinetic scheme, focusing on the spreading process when the separation bubbles form inside the droplet. The distributions of shear stress, pressure, and enstrophy flux at the interface and the wall are analyzed, particularly near the moving contact points and other characteristic points. This example gives an unique perspective to the physics of droplet impingement on a wall.