In this work, we establish the convergence of 2D, stationary Navier-Stokes flows with viscosity $\varepsilon> 0$, $(u^\varepsilon , v^\varepsilon )$ to the classical Prandtl boundary layer, $(\bar {u}_p, \bar {v}_p)$, posed on the domain $(0, \infty ) \times (0, \infty )$:

$$ \begin{align*} \| u^\varepsilon - \bar{u}_p \|_{L^\infty_y} \lesssim \sqrt{\varepsilon} \langle x \rangle^{- \frac 1 4 + \delta}, \qquad \| v^\varepsilon - \sqrt{\varepsilon} \bar{v}_p \|_{L^\infty_y} \lesssim \sqrt{\varepsilon} \langle x \rangle^{- \frac 1 2}. \end{align*} $$

This validates Prandtl’s boundary layer theory globally in the x-variable for a large class of boundary layers, including the entire one parameter family of the classical Blasius profiles, with sharp decay rates. The result demonstrates asymptotic stability in two senses simultaneously: (1) asymptotic as $\varepsilon \rightarrow 0$ and (2) asymptotic as $x \rightarrow \infty $. In particular, our result provides the first rigorous confirmation for the Navier-Stokes equations that the boundary layer cannot “separate” in these stable regimes, which is very important for physical and engineering applications.

It was proved in [11, J. Funct. Anal., 2020] that the Cauchy problem for some Oldroyd-B model is well-posed in $\dot{B}^{d/p-1}_{p,1}(\mathbb{R}^d) \times \dot{B}^{d/p}_{p,1}(\mathbb{R}^d)$ with $1\leq p \lt 2d$. In this paper, we prove that the Cauchy problem for the same Oldroyd-B model is ill-posed in $\dot{B}^{d/p-1}_{p,r}(\mathbb{R}^d) \times \dot{B}^{d/p}_{p,r}(\mathbb{R}^d)$ with $1\leq p\leq \infty$ and $1 \lt r\leq\infty$ due to the lack of continuous dependence of the solution.

This article is dedicated to investigating limit behaviours of invariant measures with respect to delay and system parameters of 3D Navier–Stokes–Voigt equations. Firstly, the well-posedness of such a system is obtained on arbitrary open sets that satisfy the Poincaré inequality, and then a unique minimal pullback attractor is attained by using the energy equation method and asymptotic compactness property. Furthermore, we construct a family of invariant Borel probability measures, which are supported on the pullback attractors. Specifically, when the external forcing terms are periodic in time, the periodic invariant measure can be obtained. Finally, as the delay approaches zero and system parameters tend to some numbers, the limit of the invariant measure sequences for this class of equations must be the invariant measure of the corresponding limit equations.

In this article, we give a comprehensive characterization of $L^1$-summability for the Navier-Stokes flows in the half space, which is a long-standing problem. The main difficulties are that $L^q-L^r$ estimates for the Stokes flow don't work in this end-point case: $q=r=1$; the projection operator $P: L^1\longrightarrow L^1_\sigma$ is not bounded any more; useful information on the pressure function is missing, which arises in the net force exerted by the fluid on the noncompact boundary. In order to achieve our aims, by making full use of the special structure of the half space, we decompose the pressure function into two parts. Then the knotty problem of handling the pressure term can be transformed into establishing a crucial and new weighted $L^1$-estimate, which plays a fundamental role. In addition, we overcome the unboundedness of the projection $P$ by solving an elliptic problem with homogeneous Neumann boundary condition.

The flow past an obstacle is a fundamental object in fluid mechanics. In 1967 Finn and Smith proved the unique existence of stationary solutions, called the physically reasonable solutions, to the Navier–Stokes equations in a two-dimensional exterior domain modeling this type of flows when the Reynolds number is sufficiently small. The asymptotic behavior of their solution at spatial infinity has been studied in detail and well understood by now, while its stability has remained open due to the difficulty specific to the two-dimensionality. In this paper, we prove that the physically reasonable solutions constructed by Finn and Smith are asymptotically stable with respect to small and well-localized initial perturbations.

The Navier-Stokes equations for viscous, incompressible fluids are studied in the three-dimensional periodic domains, with the body force having an asymptotic expansion, when time goes to infinity, in terms of power-decaying functions in a Sobolev-Gevrey space. Any Leray-Hopf weak solution is proved to have an asymptotic expansion of the same type in the same space, which is uniquely determined by the force, and independent of the individual solutions. In case the expansion is convergent, we show that the next asymptotic approximation for the solution must be an exponential decay. Furthermore, the convergence of the expansion and the range of its coefficients, as the force varies are investigated.

If the closed wire frame of a soap film having the shape of a Möbius strip is pulled apart and gradually deformed into a planar circle, the soap film transforms into a two-sided orientable surface. In the presence of a finite-time twist singularity, which changes the linking number of the film's Plateau border and the centreline of the wire, the topological transformation involves the collapse of the film toward the wire. In contrast to experimental studies of this process reported elsewhere, we use a numerical approach based on the immersed boundary method, which treats the soap film as a massless membrane in a Navier-Stokes fluid. In addition to known effects, we discover vibrating motions of the film arising after the topological change is completed, similar to the vibration of a circular membrane.

We present a moving-least-square immersed boundary method for solving viscous incompressible flow involving deformable and rigid boundaries on a uniform Cartesian grid. For rigid boundaries, noslip conditions at the rigid interfaces are enforced using the immersed-boundary direct-forcing method. We propose a reconstruction approach that utilizes moving least squares (MLS) method to reconstruct the velocity at the forcing points in the vicinity of the rigid boundaries. For deformable boundaries, MLS method is employed to construct the interpolation and distribution operators for the immersed boundary points in the vicinity of the rigid boundaries instead of using discrete delta functions. The MLS approach allows us to avoid distributing the Lagrangian forces into the solid domains as well as to avoid using the velocity of points inside the solid domains to compute the velocity of the deformable boundaries. The present numerical technique has been validated by several examples including a Poiseuille flow in a tube, deformations of elastic capsules in shear flow and dynamics of red-blood cell in microfluidic devices.

In this paper, we study an exponential time differencing method for solving the gauge system of incompressible viscous flows governed by Stokes or Navier-Stokes equations. The momentum equation is decoupled from the kinematic equation at a discrete level and is then solved by exponential time stepping multistep schemes in our approach. We analyze the stability of the proposed method and rigorously prove that the first order exponential time differencing scheme is unconditionally stable for the Stokes problem. We also present a compact representation of the algorithm for problems on rectangular domains, which makes FFT-based solvers available for the resulting fully discretized system. Various numerical experiments in two and three dimensional spaces are carried out to demonstrate the accuracy and stability of the proposed method.

This paper describes an adaptive preconditioner for numerical continuation of incompressible Navier–Stokes flows based on Stokes preconditioning [42] which has been used successfully in studies of pattern formation in convection. The preconditioner takes the form of the Helmholtz operator I–ΔtL which maps the identity (no preconditioner) for Δt≪1 to Laplacian preconditioning for Δt≫1. It is built on a first order Euler time-discretization scheme and is part of the family of matrix-free methods. The preconditioner is tested on two fluid configurations: three-dimensional doubly diffusive convection and a two-dimensional projection of a shear flow. In the former case, it is found that Stokes preconditioning is more efficient for , away from the values used in the literature. In the latter case, the simple use of the preconditioner is not sufficient and it is necessary to split the system of equations into two subsystems which are solved simultaneously using two different preconditioners, one of which is parameter dependent. Due to the nature of these applications and the flexibility of the approach described, this preconditioner is expected to help in a wide range of applications.

The lattice Boltzmann method is employed to simulate the steady flow in a two-dimensional lid-driven semi-elliptical cavity. Reynolds number (Re) and vertical-to-horizontal semi-axis ratio (D) are in the range of 500-5000 and 0.1-4, respectively. The effects of Re and D on the vortex structure and pressure field are investigated, and the evolutionary features of the vortex structure with Re and D are analyzed in detail. Simulation results show that the vortex structure and its evolutionary features significantly depend on Re and D. The steady flow is characterized by one vortex in the semi-elliptical cavity when both Re and D are small. As Re increases, the appearance of the vortex structure becomes more complex. When D is less than 1, increasing D makes the large vortexes more round, and the evolution of the vortexes with D becomes more complex with increasing Re. When D is greater than 1, the steady flow consists of a series of large vortexes which superimpose on each other. As Re and D increase, the number of the large vortexes increases. Additionally, a small vortex in the upper-left corner of the semi-elliptical cavity appears at a large Re and its size increases slowly as Re increases. The highest pressures appear in the upper-right corner and the pressure changes drastically in the upper-right region of the cavity. The total pressure differences in the semi-elliptical cavity with a fixed D decrease with increasing Re. In the region of themain vortex, the pressure contours nearly coincide with the streamlines, especially for the cavity flow with a large Re.

In this article, by applying the Stokes projection and an orthogonal projection with respect to curl and div operators, some new error estimates of finite element method (FEM) for the stationary incompressible magnetohydrodynamics (MHD) are obtained. To our knowledge, it is the first time to establish the error bounds which are explicitly dependent on the Reynolds numbers, coupling number and mesh size. On the other hand, The uniform stability and convergence of an Oseen type finite element iterative method for MHD with respect to high hydrodynamic Reynolds number Re and magnetic Reynolds number Rm , or small δ=1–σ with

(C 0, C 1 are constants depending only on Ω and F is related to the source terms of equations) are analyzed under the condition that . Finally, some numerical tests are presented to demonstrate the effectiveness of this algorithm.

In this paper, the problem of magnetohydrodynamics (MHD) boundary layer flow of nanofluid with heat and mass transfer through a porous media in the presence of thermal radiation, viscous dissipation and chemical reaction is studied. Three types of nanofluids, namely Copper (Cu)-water, Alumina (Al 2 O 3)-water and Titanium Oxide (TiO 2)-water are considered. The governing set of partial differential equations of the problem is reduced into the coupled nonlinear system of ordinary differential equations (ODEs) by means of similarity transformations. Finite element solution of the resulting system of nonlinear differential equations is obtained using continuous Galerkin-Petrov discretization together with the well-known shooting technique. The obtained results are validated using MATLAB “bvp4c” function and with the existing results in the literature. Numerical results for the dimensionless velocity, temperature and concentration profiles are obtained and the impact of various physical parameters such as the magnetic parameter M, solid volume fraction of nanoparticles 𝜙 and type of nanofluid on the flow is discussed. The results obtained in this study confirm the idea that the finite element method (FEM) is a powerful mathematical technique which can be applied to a large class of linear and nonlinear problems arising in different fields of science and engineering.

This paper presents the extension of a well-established Immersed Boundary (IB)/cut-cell method, the LS-STAG method (Y. Cheny & O. Botella, J. Comput. Phys. Vol. 229, 1043-1076, 2010), to viscoelastic flow computations in complex geometries. We recall that for Newtonian flows, the LS-STAG method is based on the finite-volume method on staggered grids, where the IB boundary is represented by its level-set function. The discretization in the cut-cells is achieved by requiring that global conservation properties equations be satisfied at the discrete level, resulting in a stable and accurate method and, thanks to the level-set representation of the IB boundary, at low computational costs.

In the present work, we consider a general viscoelastic tensorial equation whose particular cases recover well-known constitutive laws such as the Oldroyd-B, White-Metzner and Giesekus models. Based on the LS-STAG discretization of the Newtonian stresses in the cut-cells, we have achieved a compatible velocity-pressure-stress discretization that prevents spurious oscillations of the stress tensor. Applications to popular benchmarks for viscoelastic fluids are presented: the four-to-one abrupt planar contraction flows with sharp and rounded re-entrant corners, for which experimental and numerical results are available. The results show that the LS-STAG method demonstrates an accuracy and robustness comparable to body-fitted methods.

Removing geometric details from the computational domain can significantly reduce the complexity of downstream task of meshing and simulation computation, and increase their stability. Proper estimation of the sensitivity analysis error induced by removing such domain details, called defeaturing errors, can ensure that the sensitivity analysis fidelity can still be met after simplification. In this paper, estimation of impacts of removing arbitrarily constrained domain details to the analysis of incompressible fluid flows is studied with applications to fast analysis of incompressible fluid flows in complex environments. The derived error estimator is applicable to geometric details constrained by either Dirichlet or Neumann boundary conditions, and has no special requirements on the outer boundary conditions. Extensive numerical examples were presented to demonstrate the effectiveness and efficiency of the proposed error estimator.

A genuine finite volume method based on the lattice Boltzmann equation (LBE) for nearly incompressible flows is developed. The proposed finite volume lattice Boltzmann method (FV-LBM) is grid-transparent, i.e., it requires no knowledge of cell topology, thus it can be implemented on arbitrary unstructured meshes for effective and efficient treatment of complex geometries. Due to the linear advection term in the LBE, it is easy to construct multi-dimensional schemes. In addition, inviscid and viscous fluxes are computed in one step in the LBE, as opposed to in two separate steps for the traditional finite-volume discretization of the Navier-Stokes equations. Because of its conservation constraints, the collision term of the kinetic equation can be treated implicitly without linearization or any other approximation, thus the computational efficiency is enhanced. The collision with multiple-relaxation-time (MRT) model is used in the LBE. The developed FV-LBM is of second-order convergence. The proposed FV-LBM is validated with three test cases in two-dimensions: (a) the Poiseuille flow driven by a constant body force; (b) the Blasius boundary layer; and (c) the steady flow past a cylinder at the Reynolds numbers Re=10, 20, and 40. The results verify the designed accuracy and efficacy of the proposed FV-LBM.

In the first part, we study the convergence of discrete solutions to splitting schemes for two-phase flow with different mass densities suggested in [Guillen-Gonzalez, Tierra, J.Comput.Math. (6)2014]. They have been formulated for the diffuse interface model in [Abels, Garcke, Grün, M3AS, 2012, DOI:10.1142/S0218202511500138] which is consistent with thermodynamics. Our technique covers various discretization methods for phase-field energies, ranging from convex-concave splitting to difference quotient approaches for the double-well potential. In the second part of the paper, numerical experiments are presented in two space dimensions to identify discretizations of Cahn-Hilliard energies which are ϕ-stable and which do not reduce the acceleration of falling droplets. Finally, 3d simulations in axial symmetric geometries are shown to underline even more the full practicality of the approach.