Vanishing viscosity limit for a coupled Navier-Stokes/Allen-Cahn system

In this paper, we study the vanishing viscosity limit for a coupled Navier-Stokes/Allen-Cahn system in a bounded domain. We first show the local existence of smooth solutions of the Euler/Allen-Cahn equations by modified Galerkin method. Then using the boundary layer function to deal with the mismatch of the boundary conditions between Navier-Stokes and Euler equations, and assuming that the energy dissipation for Navier-Stokes equation in the boundary layer goes to zero as the viscosity tends to zero, we prove that the solutions of the Navier-Stokes/Allen-Cahn system converge to that of the Euler/Allen-Cahn system in a proper small time interval. In addition, for strong solutions of the Navier-Stokes/Allen-Cahn system in 2D, the convergence rate is cν1/2.     

Boundary Layers Associated with a Coupled Navier-Stokes/Allem-Cahn System: the Non-characteristic Boundary Case

The goal of this article is to study the boundary layer of Navier-Stokes/Allen- Cahn system in a channel at small viscosity. We prove that there exists a boundary layer at the outlet (down-wind) of thickness n, where n is the kinematic viscosity. The convergence in L^2 of the solutions of the Navier-Stokes/Allen-Cahn equations to that of the Euler/Allen-Cahn equations at the vanishing viscosity was established. In two dimensional case we are able to derive the physically relevant uniform in space and time estimates, which is derived by the idea of better control on the tangential derivative and the use of an anisotropic Sobolve imbedding.     

非特征边界的MHD方程的边界层

考查了小粘性时非特征边界情况下MHD方程在边界附近的性质,说明速度在边界上不为零.源于之前非特征边界条件下不可压缩Navier-Stokes方程边界层的工作,证明了边界层的存在性,并得到了当粘性收敛于零时,MHD方程的解收敛于理想MHD方程的解.     

Zero dissipation limit of the 1D linearized Navier-Stokes equations for a compressible fluid

In this paper we study the asymptotic limiting behavior of the solutions to the initial boundary value problem for linearized one-dimensional compressible Navier-Stokes equations. We consider the characteristic boundary conditions, that is we assume that an eigenvalue of the associated inviscid Euler system vanishes uniformly on the boundary. The aim of this paper is to understand the evolution of the boundary layer, to construct the asymptotic ansatz which is uniformly valid up to the boundary, and to obtain rigorously the uniform convergence to the solution of the Euler equations without the weakness assumption on the boundary layer.     

Continuum models in problems of the hypersonic flow of a rarefied gas around blunt bodiest

All possible continuum (hydrodynamic) models in the case of two-dimensional problems of supersonic and hypersonic flows around blunt bodies in the two-layer model (a viscous shock layer and shock-wave structure) over the whole range of Reynolds numbers, Re, from low values (free molecular and transitional flow conditions) up to high values (flow conditions with a thin leading shock wave, a boundary layer and an external inviscid flow in the shock layer) are obtained from the Navier-Stokes equations using an asymptotic analysis. In the case of low Reynolds numbers, the shock layer is considered but the structure of the shock wave is ignored. Together with the well-known models (a boundary layer, a viscous shock layer, a thin viscous shock layer, parabolized Navier-Stokes equations (the single-layer model) for high, moderate and low Re numbers, respectively), a new hydrodynamic model, which follows from the Navier-Stokes equations and reduces to the solution of the simplified (“local”) Stokes equations in a shock layer with vanishing inertial and pressure forces and boundary conditions on the unspecified free boundary (the shock wave) is found at Reynolds numbers, and a density ratio, k, up to and immediately after the leading shock wave, which tend to zero subject to the condition that (k/Re)^1/2 → 0. Unlike in all the models which have been mentioned above, the solution of the problem of the flow around a body in this model gives the free molecular limit for the coefficients of friction, heat transfer and pressure. In particular, the Newtonian limit for the drag is thereby rigorously obtained from the Navier-Stokes equations. At the same time, the Knudsen number, which is governed by the thickness of the shock layer, which vanishes in this model, tends to zero, that is, the conditions for a continuum treatment are satisfied. The structure of the shock wave can be determined both using continuum as well as kinetic models after obtaining the solution in the viscous shock layer for the weak physicochemical processes in the shock wave structure itself. Otherwise, the problem of the shock wave structure and the equations of the viscous shock layer must be jointly solved. The equations for all the continuum models are written in Dorodnitsyn--Lees boundary layer variables, which enables one, prior to solving the problem, to obtain an approximate estimate of second-order effects in boundary-layer theory as a function of Re and the parameter k and to represent all the aerodynamic and thermal characteristic; in the form of a single dependence on Re over the whole range of its variation from zero to infinity.

An efficient numerical method of global iterations, previously developed for solving viscous shock-layer equations, can be used to solve problems of supersonic and hypersonic flows around the windward side of blunt bodies using a single hydrodynamic model of a viscous shock layer for all Re numbers, subject to the condition that the limit (k/Re)^1/2 → 0 is satisfied in the case of small Re numbers. An aerodynamic and thermal calculation using different hydrodynamic models, corresponding to different ranges of variation Re (different types of flow) can thereby, in fact, be replaced by a single calculation using one model for the whole of the trajectory for the descent (entry) of space vehicles and natural cosmic bodies (meteoroids) into the atmosphere.     

Boundary Layers Associated with Incompressible Navier-Stokes Equations: The Noncharacteristic Boundary Case

The goal of this article is to study the boundary layer of wall bounded flows in a channel at small viscosity when the boundaries are uniformly noncharacteristic, i.e., there is injection and/or suction everywhere at the boundary. Following earlier work on the boundary layer for linearized Navier-Stokes equations in the case where the boundaries are characteristic (no-slip at the boundary and non-permeable), we consider here the case where the boundary is permeable and thus noncharacteristic. The form of the boundary layer and convergence results are derived in two cases: linearized equation and full nonlinear equations. We prove that there exists a boundary layer at the outlet (downwind) of the form e^−Uz/ε where U is the speed of injection/suction at the boundary, z is the distance to the outlet of the channel, and ε is the kinematic viscosity. We improve an earlier result of S. N. Alekseenko (1994, Siberian Math. J.35, No. 2, 209-230) where the convergence in L² of the solutions of the Navier-Stokes equations to that of the Euler equations at vanishing viscosity was established. In the two dimensional case we are able to derive the physically relevant uniform in space (L^∞ norm) estimates of the boundary layer. The uniform in space estimate is derived by properly developing our previous idea of better control on the tangential derivative and the use of an anisotropic Sobolev imbedding. To the best of our knowledge this is the first rigorously proved result concerning boundary layers for the full (nonlinear) Navier-Stokes equations for incompressible fluids.     

A Non-Conforming Least Squares Spectral Element Formulation for Oseen Equations with Applications to Navier-Stokes Equations

In this article, we propose a non-conforming exponentially accurate least-squares spectral element method for Oseen equations in primitive variable formulation that is applicable to both two- and three-dimensional domains. First-order reformulation is avoided, and the condition number is controlled by a suitable preconditioner for velocity components and pressure variable. A preconditioned conjugate gradient method is used to obtain the solution. Navier-Stokes equations in primitive variable formulation have been solved by solving a sequence of Oseen type iterations. For numerical test cases, similar order convergence has been achieved for all Reynolds number cases at the cost of higher iteration number for higher Reynolds number.     

Semi-implicit Runge-Kutta schemes for the Navier-Stokes equations

The stationary Navier-Stokes equations are solved in 2D with semi-implicit Runge-Kutta schemes, where explicit time-integration in the streamwise direction is combined with implicit integration in the body-normal direction. For model problems stability restrictions and convergence properties are studied. Numerical experiments for the flow over a flat plate show that the number of iterations for the semi-implicit schemes is almost independent of the Reynolds number.     

The transition between the Stokes equations and the Reynolds equation: A mathematical proof

The Reynolds equation is used to calculate the pressure distribution in a thin layer of lubricant film between two surfaces. Using the asymptotic expansion in the Stokes equations, we show the existence of singular perturbation phenomena whenever the two surfaces are in relative motion. We prove that the Reynolds equation is an approximation of the Stokes equations and that the kind of convergence is strongly related with the boundary conditions on the velocity field.     

Zero Dissipation Limit to Rarefaction Waves for the 1-D Compressible Navier-Stokes Equations

The zero dissipation limit for the one-dimensional Navier-Stokes equations of compressible,isentropic gases in the case that the corresponding Euler equations have rarefaction wave solutions is investi...     

DNS and Lie group analysis of turbulent boundarly layer flow

Lie group approach is used to derive new scaling laws for zero-pressure gradient turbulent boundary layer flow. A direct numerical simulation of the flow at Reynolds number Re_θ = 2240 was performed for the verify theoretical results. Navier-Stokes equations were numerically solved using spectral method with up to 160 million grid points. The numerical simulation shows validity of the theoretical results. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

椭圆-抛物偏微分方程奇异摄动问题的差分解法

本文讨论了含有小参数在高阶导数项的椭圆型方程奇异摄动问题的差分解法.当ε=0时椭圆型方程退化为抛物型方程.作者根据此问题解的边界层性质,构造了特殊的差分格式:研究了它的收敛性和解的渐近性态.最后给出一个数值例题.     

Saturation-free numerical scheme for computing the flow past a lattice of airfoils and the determination of separation points in a viscous fluid

A numerical method for computing the potential flow past a lattice of airfoils is described. The problem is reduced to a linear integrodifferential equation on the lattice contour, which is then approximated by a linear system of equations with the help of specially derived quadrature formulas. The quadrature formulas exhibit exponential convergence in the number of points on an airfoil and have a simple analytical form. Due to its fast convergence and high accuracy, the method can be used to directly optimize the airfoils as based on any given integral characteristics. The shear stress distribution and the separation points are determined from the velocity distribution at the airfoil boundary calculated by solving the boundary layer equations. The method proposed is free of laborious grid generation procedures and does not involve difficulties associated with numerical viscosity at high Reynolds numbers.     

Remarks on Vanishing Viscosity Limits for the 3D Navier-Stokes Equations with a Slip Boundary Condition

The authors study vanishing viscosity limits of solutions to the 3-dimensional incompressible Navier-Stokes system in general smooth domains with curved boundaries for a class of slip boundary conditions. In contrast to the case of flat boundaries, where the uniform convergence in super-norm can be obtained, the asymptotic behavior of viscous solutions for small viscosity depends on the curvature of the boundary in general. It is shown, in particular, that the viscous solution converges to that of the ideal Euler equations in C([0, T];H1(Ω)) provided that the initial vorticity vanishes on the boundary of the domain.     

Wall effects on the slow steady motion of a particle in a viscous incompressible fluid

In this paper, asymptotic expansions with respect to small Reynolds numbers are proved for the slow steady motion of an arbitrary particle in a viscous, incompressible fluid past a single plane wall. The flow problem is modelled by a certain boundary value problem for the stationary, nonlinear Navier-Stokes equations. The coefficients of these expansions are the solutions of various, linear Stokes problems which can be constructed by single layer potentials and corresponding boundary integral equations on the boundary surface of the particle. Furthermore, some asymptotic estimates at small Reynolds numbers are presented for the slow steady motion of an arbitrary particle in a viscous, incompressible fluid between two parallel, plane walls and in an infinitely long, rectilinear cylinder of arbitrary cross section. In the case of the flow problem with a single plane wall, the paper is based on a thesis, which the author has written under the guidance of Professor Dr. Wolfgang L. Wendland.     

Boundary Shape Control of the Navier-Stokes Equations and Applications

In this paper, the geometrical design for the blade's surface in an impeller or for the profile of an aircraft, is modeled from the mathematical point of view by a boundary shape control problem for the Navier-Stokes equations. The objective function is the sum of a global dissipative function and the power of the fluid. The control variables are the geometry of the boundary and the state equations are the Navier-Stokes equations. The Euler-Lagrange equations of the optimal control problem are derived, which are an elliptic boundary value system of fourth order, coupled with the Navier-Stokes equations. The authors also prove the existence of the solution of the optimal control problem, the existence of the solution of the Navier-Stokes equations with mixed boundary conditions, the weak continuity of the solution of the Navier-Stokes equations with respect to the geometry shape of the blade's surface and the existence of solutions of the equations for the Gateaux derivative of the solution of the Navier-Stokes equations with respect to the geometry of the boundary.     

An application of a boundary element method to natural convection

The direct boundary element method is applied to the numerical modelling of thermal fluid flow in a transient state. The Navier-Stokes equations are considered under the Boussinesq approximation and the viscous thermal flow equations are expressed in terms of stream function, vorticity, and temperature in two dimensions. Boundary integral equations are derived using logarithmic potential and time-dependent heat potential as fundamental solutions. Boundary unknowns are discretized by linear boundary elements and flow domains are divided into a series of triangular cells. Charged points are translated upstream in the numerical evaluation of convective terms. Unknown stream function, vorticity, and temperature are staggered in the computational scheme.

Simple iteration is found to converge to the quasi steady-state flow. Boundary solutions for two-dimensional examples at a Reynolds number 100 and Grashoff number 10⁷ are obtained.     

A High-Order Kernel-Free Boundary Integral Method for Incompressible Flow Equations in Two Space Dimensions

This paper presents a fourth-order kernel-free boundary integral method for the time-dependent, incompressible Stokes and Navier-Stokes equations defined on irregular bounded domains. By the stream function-vorticity formulation, the incompressible flow equations are interpreted as vorticity evolution equations. Time discretization methods for the evolution equations lead to a modified Helmholtz equation for the vorticity, or alternatively, a modified biharmonic equation for the stream function with two clamped boundary conditions. The resulting fourth-order elliptic boundary value problem is solved by a fourth-order kernel-free boundary integral method, with which integrals in the reformulated boundary integral equation are evaluated by solving corresponding equivalent interface problems, regardless of the exact expression of the involved Green's function. To solve the unsteady Stokes equations, a four-stage composite backward differential formula of the same order accuracy is employed for time integration. For the Navier-Stokes equations, a three-stage third-order semi-implicit Runge-Kutta method is utilized to guarantee the global numerical solution has at least third-order convergence rate. Numerical results for the unsteady Stokes equations and the Navier-Stokes equations are presented to validate efficiency and accuracy of the proposed method.     

Numerical study of spherical Couette flows for certain zenith-angle-dependent rotations of boundary spheres at low Reynolds numbers

The numerical method with splitting of boundary conditions developed previously by the first and third authors for solving the stationary Dirichlet boundary value problem for the Navier-Stokes equations in spherical layers in the axisymmetric case at low Reynolds numbers and a corresponding software package were used to study viscous incompressible steady flows between two con-centric spheres. Flow regimes depending on the zenith angle ?? of coaxially rotating boundary spheres (admitting discontinuities in their angular velocities) were investigated. The orders of accuracy with respect to the mesh size of the numerical solutions (for velocity, pressure, and stream function in a meridional plane) in the max and L ₂ norms were studied in the case when the velocity boundary data have jump discontinuities and when some procedures are used to smooth the latter. The capabilities of the Richardson extrapolation procedure used to improve the order of accuracy of the method were investigated. Error estimates were obtained. Due to the high accuracy of the numerical solutions, flow features were carefully analyzed that were not studied previously. A number of interesting phenomena in viscous incompressible flows were discovered in the cases under study.     

Solving the Dirichlet problem for Navier–Stokes equations by probabilistic approach

We construct a number of layer methods for Navier-Stokes equations (NSEs) with no-slip boundary conditions. The methods are obtained using probabilistic representations of solutions to NSEs and exploiting ideas of the weak sense numerical integration of stochastic differential equations. Despite their probabilistic nature, the proposed methods are nevertheless deterministic.