Incompressible Navier-Stokes Equations Research Articles

Numerical simulation of droplet impact onto heated surfaces below the boiling point

Droplet impact onto heated surfaces is a widespread process in industrial applications, particularly in the context of spray cooling techniques. Therefore, it is essential to study the complex phenomenon of droplet spreading, heat removal and flow distributions during this process. The primary focus of this research centers on the wetting dynamics at a surface temperature below the saturation temperature of the liquid. This study employs a Computational fluid dynamics (CFD) framework to simulate the impact dynamics of a water droplet onto a heated surface. The impact dynamics, spreading, air entrapment and heat transfer in the process are discussed using simulation results. The Direct Numerical Simulation (DNS) tool Free Surface 3D (FS3D), an in-house code at the Institute of Aerospace Thermodynamics, University of Stuttgart, is utilized. The Finite-Volume method is used to solve the incompressible Navier-Stokes equations, and the interface is defined by the Volume of fluid (VOF) method. In this study, the temperature change above the entrapped air bubble formed due to droplet inertia is shown, which can influence the heat removal rate in the process. Therefore, conducting a comprehensive study on the temperature shift caused by the air bubble is crucial.

Maximal $L^{1}$-regularity and free boundary problems for the incompressible Navier–Stokes equations in critical spaces

Maximal $L^{1}$-regularity and free boundary problems for the incompressible Navier–Stokes equations in critical spaces

Arbitrary-order sensitivities of the incompressible base flow and its eigenproblem

First-order sensitivities and adjoint analysis are used widely to control the linear stability of unstable flows. Second-order sensitivities have recently helped to increase accuracy. In this paper, a method is presented to calculate arbitrary high-order sensitivities based on Taylor expansions of the incompressible base flow and its eigenproblem around a scalar parameter. For the incompressible Navier–Stokes equations, general expressions for the sensitivities are derived, into which parameter-specific information can be inserted. The computational costs are low since, for all orders, a linear equation system has to be solved, of which the left-hand-side matrix stays constant and thus its preconditioning can be exploited. Two flow scenarios are examined. First, the cylinder flow equations are expanded around the inverse of the Reynolds number, enabling the prediction of the two-dimensional cylinder base flow and its leading eigenvalue as a function of the Reynolds number. This approach computes accurately the base flow and eigenvalue even in the unstable regime, providing, when executed subsequently, a mean to calculate unstable base flows. This case gives a clear introduction into the method and allows us to discuss its constraints regarding convergence behaviour. Second, a small control cylinder is introduced into the domain of the cylinder flow for stabilization. Higher-order sensitivity maps are calculated by modelling the small cylinder with a steady forcing. These maps help to identify stabilizing areas of the flow field for Reynolds numbers within the laminar vortex shedding regime, with the required number of orders increasing as the Reynolds number rises. The results obtained through the proposed method align well with numerically calculated eigenvalues that incorporate the cylinder directly into the grid.

The Exact Limits and Improved Decay Estimates for All Order Derivatives of the Global Weak Solutions of n-Dimensional Incompressible Navier-Stokes Equations

The Exact Limits and Improved Decay Estimates for All Order Derivatives of the Global Weak Solutions of n-Dimensional Incompressible Navier-Stokes Equations

A robust computational framework for variational data assimilation of mean flows with sparse measurements corrupted by strong outliers

A variational framework for the reconstruction of time-averaged mean flows using a sparse set of observations with large magnitudes of noise (referred to as outliers), is presented. The observations constitute a set of point-wise measurements of the mean flow with outliers at certain measurement locations and are incorporated into a numerical simulation governed by the two-dimensional, incompressible Reynolds-averaged Navier-Stokes (RANS) equations with an unknown momentum forcing. This forcing, which corresponds to the divergence of the Reynolds stress tensor, is calculated from a direct-adjoint optimization procedure to reduce the deviation between the measured and estimated mean velocities. ℓ2, ℓ1, Huber, and hybrid loss functions are used to represent the discrepancy in the mean velocity field between the measurements and the predictions. A variety of algorithms are considered to solve the optimization problem with these loss functions and a performance comparison in terms of the quality and physical features of the recovered mean flow is presented. The Huber loss function performed best as it remained robust to strong outliers in the measurements with its ℓ1 contribution and also ensured the uniqueness of the optimal solution with its ℓ2 contribution. Huber loss functions restrict the effect of outliers at the local measurement locations, thereby not affecting the quality of the high-dimensional reconstructed mean flow field. The hybrid loss function, a modified form of the continuous Huber loss function, also recovered the mean flow with high accuracy. We demonstrate the performance of the data assimilation framework for the case of two-dimensional laminar flow around a circular cylinder at Re=100. We then extend the analysis to the case of two-dimensional laminar flow over a backward-facing step at Re=500, to further assess the efficacy and robustness of the data assimilation framework.

Asymptotic behavior of the 3D incompressible Navier–Stokes equations with damping

In this paper, we consider the 3D incompressible Navier–Stokes equations with damping term |u|β−1u(β≥1). First, by using a different and simple method from Cai and Lei (2010), Jia et al. (2011), Jiang (2012) and Yu and Zheng (2019), for any β≥1, we prove that the weak solutions decay to zero in L2 as time tends to infinity; for any β≥3, we derive optimal decay rates of the L2-norm of the solutions. Second, we obtain the decay rate with some appropriate space weighted estimates, which is the first result on the 3D damped Navier–Stokes equations to our knowledge.

Global regularity and decay behavior for Leray equations with critical-dissipation and its application to self-similar solutions

In this paper, we show the global regularity and the optimal decay of weak solutions to the generalized Leray problem with critical dissipation. Our approach hinges on the maximal smoothing effect, L p L^{p} -type elliptic regularity of linearization, and the action of the heat semigroup generated by the fractional powers of Laplace operator on distributions with Fourier transforms supported in an annulus. As a by-product, we construct a self-similar solution to the three-dimensional incompressible Navier-Stokes equations. Most notably, we prove the global regularity and the optimal decay without the need for additional requirements found in existing literatures.

The Cut-Cell Method for the Conjugate Heat Transfer Topology Optimization of Turbulent Flows Using the “Think Discrete–Do Continuous” Adjoint

This paper presents a topology optimization (TopO) method for conjugate heat transfer (CHT), with turbulent flows. Topological changes are controlled by an artificial material distribution field (design variables), defined at the cells of a background grid and used to distinguish a fluid from a solid material. To effectively solve the CHT problem, it is crucial to impose exact boundary conditions at the computed fluid–solid interface (FSI); this is the purpose of introducing the cut-cell method. On the grid, including also cut cells, the incompressible Navier–Stokes equations, coupled with the Spalart–Allmaras turbulence model with wall functions, and the temperature equation are solved. The continuous adjoint method computes the derivatives of the objective function(s) and constraints with respect to the material distribution field, starting from the computation of derivatives with respect to the positions of nodes on the FSI and then applying the chain rule of differentiation. In this work, the continuous adjoint PDEs are discretized using schemes that are consistent with the primal discretization, and this will be referred to as the “Think Discrete–Do Continuous” (TDDC) adjoint. The accuracy of the gradient computed by the TDDC adjoint is verified and the proposed method is assessed in the optimization of two 2D cases, both in turbulent flow conditions. The performance of the TopO designs is investigated in terms of the number of required refinement steps per optimization cycle, the Reynolds number of the flow, and the maximum allowed power dissipation. To illustrate the benefits of the proposed method, the first case is also optimized using a density-based TopO that imposes Brinkman penalization terms in solid areas, and comparisons are made.

Layer separation of the 3D incompressible Navier–Stokes equation in a bounded domain

We provide an unconditional L 2 upper bound for the boundary layer separation of Leray–Hopf solutions in a smooth bounded domain. By layer separation, we mean the discrepancy between a (turbulent) low-viscosity Leray–Hopf solution u ν and a fixed (laminar) regular Euler solution u ¯ with similar initial conditions and body force. We show an asymptotic upper bound C | | u ¯ | | L ∞ 3 T on the layer separation, anomalous dissipation, and the work done by friction. This extends the previous result when the Euler solution is a regular shear in a finite channel. The key estimate is to control the boundary vorticity in a way that does not degenerate in the vanishing viscosity limit.

A new compact scheme-based Lax–Wendroff method for high fidelity simulations

The classical Lax–Wendroff (LW) method employs explicit second order central difference schemes to solve partial differential equations. Recently, the authors in Sengupta et al. (2023) have identified numerical parameter ranges in performing large eddy simulation (LES) using this classical LW method. In a bid to improve these numerical parameter ranges for high fidelity computations, a new non-uniform grid based compact scheme for LW method (NUCLW) is developed here. This new method uses sixth order non-uniform compact scheme (NUC6) for the first and higher order derivatives. Reported global spectral analysis confirms substantial improvement in the resolution and accuracy of NUCLW scheme over the classical LW scheme. A significant advantage of the present scheme is its ability to compute solutions for non-uniform, structured grids, and that does not suffer degradation in accuracy. The potential of the NUCLW scheme for high accuracy computations is demonstrated by solving 2D incompressible Navier–Stokes equations for the benchmark unsteady flow inside a square lid driven cavity problem for a subcritical Reynolds number of 3200 and a post-critical Reynolds number of 10,000 and simulating the instability of the Taylor–Green vortex problem. These demonstrate the ability of the NUCLW scheme to solve the Navier–Stokes equations for high fidelity simulations such as LES/DNS with enhanced ranges of numerical parameters with respect to enhanced accuracy and faster computations.

Augmented Lagrangian Acceleration of Global-in-Time Pressure Schur Complement Solvers for Incompressible Oseen Equations

This work is focused on an accelerated global-in-time solution strategy for the Oseen equations, which highly exploits the augmented Lagrangian methodology to improve the convergence behavior of the Schur complement iteration. The main idea of the solution strategy is to block the individual linear systems of equations at each time step into a single all-at-once saddle point problem. By elimination of all velocity unknowns, the resulting implicitly defined equation can then be solved using a global-in-time pressure Schur complement (PSC) iteration. To accelerate the convergence behavior of this iterative scheme, the augmented Lagrangian approach is exploited by modifying the momentum equation for all time steps in a strongly consistent manner. While the introduced discrete grad-div stabilization does not modify the solution of the discretized Oseen equations, the quality of customized PSC preconditioners drastically improves and, hence, guarantees a rapid convergence. This strategy comes at the cost that the involved auxiliary problem for the velocity field becomes ill conditioned so that standard iterative solution strategies are no longer efficient. Therefore, a highly specialized multigrid solver based on modified intergrid transfer operators and an additive block preconditioner is extended to solution of the all-at-once problem. The potential of the proposed overall solution strategy is discussed in several numerical studies as they occur in commonly used linearization techniques for the incompressible Navier–Stokes equations.

Variational–hemivariational system for contaminant convection–reaction–diffusion model of recovered fracturing fluid

Abstract This work is devoted to study the convection–reaction–diffusion behavior of contaminant in the recovered fracturing fluid which flows in the wellbore from shale gas reservoir. First, we apply various constitutive laws for generalized non-Newtonian fluids, diffusion principles, and friction relations to formulate the recovered fracturing fluid model. The latter is a partial differential system composed of a nonlinear and nonsmooth stationary incompressible Navier-Stokes equation with a multivalued friction boundary condition, and a nonlinear convection–reaction–diffusion equation with mixed Neumann boundary conditions. Then, we provide the weak formulation of the fluid model which is a hemivariational inequality driven by a nonlinear variational equation. We establish existence of solutions to the recovered fracturing fluid model via a surjectivity theorem for multivalued operators combined with an alternative iterative method and elements of nonsmooth analysis.

Convergence analysis of high-order IMEX-BDF schemes for the incompressible Navier–Stokes equations

In this paper, we consider developing high-order temporal integration schemes for the unsteady incompressible Navier–Stokes equations in bounded two-dimensional domain subjected to the periodic boundary conditions. Utilizing the k-step (k=3,4,5) backward differentiation formula (BDF) coupled with the implicit–explicit (IMEX) treatment of the nonlinear convective term in an anti-symmetry form, a class of IMEX-BDFk schemes up to fifth-order in time are constructed and analyzed. By imposing a zero-mean constrain on the finite-dimensional space for the pressure, the proposed numerical schemes are proven to be uniquely solvable. Based on the recent theoretical framework consisting of a class of discrete orthogonal convolution kernels, rigorous L2 norm error estimates for both the velocity and the pressure are established by using a novel divergence free projection system. The proposed schemes are then implemented in two benchmark experiments, including a Taylor–Green vortex problem and a double shear layer flow at various high Reynolds numbers. Numerical results demonstrate the expected solution accuracy and the computational effectiveness in simulating the realistic flow dynamics.

Unconditional H2$$ {H}^2 $$‐stability of the Euler implicit/explicit SAV‐based scheme for the 2D Navier–Stokes equations with smooth or nonsmooth initial data

AbstractIn this article, we propose ‐unconditional stable schemes for solving time‐dependent incompressible Navier–Stokes equations with smooth or nonsmooth initial data, , . The ‐stability analysis is established by leveraging the scalar auxiliary variable (SAV) approach. When dealing with nonsmooth initial data, we utilize a limited number of iteration of the semi‐implicit scheme followed by the SAV scheme. The overall efficiency is greatly enhanced due to the minimal computational cost of the semi‐implicit scheme and the explicit treatment of the nonlinear term within the SAV approach. The proposed schemes investigate two types of scalar auxiliary variables: the energy‐based variable and the exponential‐based variable. Rigorous proofs of the ‐unconditional stability of both schemes have been provided. Notice that both proposed numerical schemes enjoy unconditional long time stability for smooth and nonsmooth initial data when . Numerical experiments have been conducted to demonstrate the theoretical results.

Two-phase Porous Media Flow Model Based on the Incompressible Navier-Stokes Equation.

Two-phase porous media flow is important in many applications from drug delivery to groundwater diffusion and oil recovery and is of particular interest to biomedical diagnostic test developers using cellulose and nitrocellulose membranes with limited fluid sample volumes. This work presents a new two-phase porous media flow model based on the incompressible Navier-Stokes equation. The model aims to address the limitations of existing methods by incorporating a partial saturation distribution in porous media to account for limited fluid volumes. The basic parameters of the model are the pore size distribution and the contact angle. To validate the model, we solved five analytical solutions and compared them to corresponding experimental data. The experimentally measured penetration length data agreed with the model predictions, demonstrating model accuracy. Our findings suggest that this new two-phase porous media flow model can provide a valuable tool for researchers developing fluidic assays in paper and other porous media.

Free boundary problems of the incompressible Navier–Stokes equations with non-flat initial surface in the critical Besov space

Global well-posedness of the Navier–Stokes equations with a free boundary condition is considered in the scaling critical homogeneous Besov spaces B˙p,1-1+n/p(R+n)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\dot{B}}_{p,1}^{-1+n/p}({\\mathbb {R}}^n_+)$$\\end{document} with n-1<p<2n-1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$n-1< p< 2n-1$$\\end{document}. To show the global well-posedness, we establish end-point maximal L1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L^1$$\\end{document}-regularity for the initial-boundary value problem of the Stokes equations. Such an estimate is obtained via related estimate for the initial-boundary value problem of the heat equation with the inhomogeneous Neumann data as well as the pressure estimate in the critical Besov space framework. The proof heavily depends on the explicit expression of the fundamental integral kernel of the Lagrange transformed linearized Stokes equations and the almost orthogonal estimates with the space-time Littlewood–Paley dyadic decompositions. Our result here improves the initial space and boundary state than previous results by Danchin–Hieber–Mucha–Tolksdorf (Free boundary problems via Da Prato–Grisvard theory. arXiv:2011.07918v2) and ourselves (Ogawa and Shimizu in J Evol Equ 22(30):67, 2022; Ogawa and Shimizu in J Math Soc Jpn. arXiv:2211.06952v3).

Strong solutions and attractor dimension for 2D NSE with dynamic boundary conditions

We consider incompressible Navier–Stokes equations in a bounded 2D domain, complete with the so-called dynamic slip boundary conditions. Assuming that the data are regular, we show that weak solutions are strong. As an application, we provide an explicit upper bound of the fractal dimension of the global attractor in terms of the physical parameters. These estimates comply with analogous results in the case of Dirichlet boundary condition.

Monge–Ampère geometry and vortices

We introduce a new approach to Monge–Ampère geometry based on techniques from higher symplectic geometry. Our work is motivated by the application of Monge–Ampère geometry to the Poisson equation for the pressure that arises for incompressible Navier–Stokes flows. Whilst this equation constitutes an elliptic problem for the pressure, it can also be viewed as a non-linear partial differential equation connecting the pressure, the vorticity, and the rate-of-strain. As such, it is a key diagnostic relation in the quest to understand the formation of vortices in turbulent flows. We study this equation via an associated (higher) Lagrangian submanifold in the cotangent bundle to the configuration space of the fluid. Using our definition of a (higher) Monge–Ampère structure, we study an associated metric on the cotangent bundle together with its pull-back to the (higher) Lagrangian submanifold. The signatures of these metrics are dictated by the relationship between vorticity and rate-of-strain, and their scalar curvatures can be interpreted in a physical context in terms of the accumulation of vorticity, strain, and their gradients. We show explicity, in the case of two-dimensional flows, how topological information can be derived from the Monge–Ampère geometry of the Lagrangian submanifold. We also demonstrate how certain solutions to the three-dimensional incompressible Navier–Stokes equations, such as Hill’s spherical vortex and an integrable case of Arnol’d–Beltrami–Childress flow, have symmetries that facilitate a formulation of these solutions from the perspective of (higher) symplectic reduction.

A centered limited finite volume approximation of the momentum convection operator for low‐order nonconforming face‐centered discretizations

AbstractWe propose in this article a discretization of the momentum convection operator for fluid flow simulations on quadrangular or generalized hexahedral meshes. The space discretization is performed by the low‐order nonconforming Rannacher–Turek finite element: the scalar unknowns are associated with the cells of the mesh while the velocities unknowns are associated with the edges or faces. The momentum convection operator is of finite volume type, and its expression is derived, as in MUSCL schemes, by a two‐step technique: computation of a tentative flux, here, with a centered approximation of the velocity, and limitation of this flux using monotonicity arguments. The limitation procedure is of algebraic type, in the sense that its does not invoke any slope reconstruction, and is independent from the geometry of the cells. The derived discrete convection operator applies both to constant or variable density flows and may thus be implemented in a scheme for incompressible or compressible flows. To achieve this goal, we derive a discrete analogue of the computation (with the velocity, one of its component, the density, and assuming that the mass balance holds) and discuss two applications of this result: first, we obtain stability results for a semi‐implicit in time scheme for incompressible and barotropic compressible flows; second, we build a consistent, semi‐implicit in time scheme that is based on the discretization of the internal energy balance rather than the total energy. The performance of the proposed discrete convection operator is assessed by numerical tests on the incompressible Navier–Stokes equations, the barotropic and the full compressible Navier–Stokes equations and the compressible Euler equations.

Analysis of quasi-variational–hemivariational inequalities with applications to Bingham-type fluids

In this paper we study the sensitivity analysis of elliptic quasi-variational–hemivariational inequalities with constraint. The upper semicontinuity property of the solution map with respect to a parameter is established. An application to the steady-state incompressible Navier–Stokes equation with mixed boundary conditions in a model for a generalized Newtonian fluid of Bingham-type is provided. The boundary conditions represent a generalization of the no leak condition, and a multivalued and nonmonotone version of a nonlinear Navier–Fujita frictional slip condition. Furthermore, a sensitivity result is proved for the weak formulation of the model when all the data are subjected to perturbations. Finally, for the Bingham-type fluids, an optimal control problem is studied.