Note on the uniqueness of subsonic Euler flows in an axisymmetric nozzle

In this paper, we consider the uniqueness of globally subsonic compressible flows through an infinitely long axisymmetric nozzle. The flow is governed by the steady Euler equations and satisfies no-flow boundary conditions on the nozzle walls. We will show that for given mass flux and Bernoulli’s function in the upstream, the subsonic flow is unique in the class of all axisymmetric solutions, which possess the asymptotic behaviors at the far fields. This result extends the uniqueness of solutions in the previous paper Du and Duan (2011) [1].     

A multigrid method with unstructured adaptive grids for steady Euler equations

The multigrid method based on multi-stage Jacobi relaxation, earlier developed by the authors for structured grid calculations with Euler equations, is extended to unstructured grid applications. The meshes are generated with Delaunay triangulation algorithms and are adapted to the flow solution.     

延迟微分方程隐式Euler方法的收敛性

本文讨论了用隐式Euler方法求解一类延迟量满足Lipschitz条件且Lipschitz常数小于1的非线性变延迟微分方程初值问题的收敛性.获得了带线性插值的隐式Euler方法的收敛性结果.     

Multi-staging of Jacobi relaxation to improved smoothing properties of multigrid methods for steady Euler equations

Multi-stage versions of Jacobi relaxation are studied for use in multigrid methods for steady Euler equations. It is shown that these multi-stage versions allow much more general and much more efficient multigrid methods than possible with classic relaxation methods.     

Numerical algorithm for solving diffusion equations on the basis of multigrid methods

A new effective algorithm based on multigrid methods is proposed for solving parabolic equations. The algorithm preserves implicit-scheme advantages (such as stability, accuracy, and conservativeness) while it involves a considerably reduced amount of arithmetic operations at every time level. The absolute stability, conservativeness, and convergence of the algorithm is proved theoretically using one- and two-dimensional initial-boundary value model problems for the heat equation. The error of the solution is estimated. The good accuracy of the method is demonstrated using two-dimensional model problems, including ones with discontinuous coefficients.     

Algebraic multigrid preconditioning of the Hessian in optimization constrained by a partial differential equation

We construct an algebraic multigrid (AMG) based preconditioner for the reduced Hessian of a linear‐quadratic optimization problem constrained by an elliptic partial differential equation. While the preconditioner generalizes a geometric multigrid preconditioner introduced in earlier works, its construction relies entirely on a standard AMG infrastructure built for solving the forward elliptic equation, thus allowing for it to be implemented using a variety of AMG methods and standard packages. Our analysis establishes a clear connection between the quality of the preconditioner and the AMG method used. The proposed strategy has a broad and robust applicability to problems with unstructured grids, complex geometry, and varying coefficients. The method is implemented using the Hypre package and several numerical examples are presented.     

Global uniqueness of steady subsonic Euler flows in three-dimensional ducts

We show for certain boundary conditions, under suitable assumptions on the velocity field, the steady compressible Euler flows in three-dimensional straight ducts must be irrotational, and hence, prove global uniqueness of uniform subsonic flows in these ducts. The proof depends on careful analysis of transport equations of vorticity and theory of second-order elliptic equations in bounded or unbounded domains.     

On the Convergence Rate of Euler Scheme for SDE with Lipschitz Drift and Constant Diffusion

It is well known that the weak Euler approximation of a stochastic differential equation has order one, provided the coefficients of the equation are sufficiently smooth. We prove that the order of the approximation is still one in the case where the drift coefficient is a Lipschitz function and the diffusion coefficient is constant.     

Comments on algorithms for grid interfaces in simulating Euler flows

Some results concerning the algorithms for grid interfaces, which are crucial in simulating flows by zonal methods, are presented in this paper. It is indicated that the commonly used conservative interface scheme can ensure the discrete entropy condition, but it may be inconsistent and would bring a nonoverlapping solution on overlapping grids. A nonconservative interface matching obtained by interpolation can be monotonicity preserving, and it leads large conservation error when discontinuities are close to the interfaces. Methods for improvement of interface algorithms are also proposed.     

改进欧拉格式求解随机微分方程的收敛性

采用改进的欧拉格式求解随机微分方程,当方程的偏移系数和扩散系数均满足全局Lipschitz条件和线性增长条件时,证明改进格式的强收敛的阶是1/2.     

Analysis of tensor product multigrid

We consider anisotropic second order elliptic boundary value problems in two dimensions, for which the anisotropy is exactly aligned with the coordinate axes. This includes cases where the operator features a singular perturbation in one coordinate direction, whereas its restriction to the other direction remains neatly elliptic. Most prominently, such a situation arises when polar coordinates are introduced.The common multigrid approach to such problems relies on line relaxation in the direction of the singular perturbation combined with semi-coarsening in the other direction. Taking the idea from classical Fourier analysis of multigrid, we employ eigenspace techniques to separate the coordinate directions. Thus, convergence of the multigrid method can be examined by looking at one-dimensional operators only. In a tensor product Galerkin setting, this makes it possible to confirm that the convergence rates of the multigrid V-cycle are bounded independently of the number of grid levels involved. In addition, the estimates reveal that convergence is also robust with respect to a singular perturbation in one coordinate direction.Finally, we supply numerical evidence that the algorithm performs satisfactorily in settings more general than those covered by the proof.     

Stability of Multidimensional Phase Transitions

In this paper, the author studies the multidimensional stability of subsonic phase transitions in a steady supersonic flow of van der Waals type. The viscosity capillarity criterion （in ＂Arch. Rat. Mech. Anal., 81（4）, 1983, 301-315＂） is used to seek physical admissible planar waves. By showing the Lopatinski determinant being non-zero, it is proved that subsonic phase transitions are uniformly stable in the sense of Majda （in ＂Mem. Amer. Math. Soc., 41（275）, 1983, 1-95＂） under both one dimensional and multidimensional perturbations.     

非线性变延迟微分方程隐式Euler方法的数值稳定性

在减弱对非线性刚性变延迟微分方程初值问题本身的约束条件的前提下 ,将已有的文献中隐式Euler方法数值稳定性的结论由常延迟的情形推广到了变延迟的情形 ,证明了隐式Euler方法是稳定的     

A large deviations principle for stochastic flows of viscous fluids

We study the well-posedness of a stochastic differential equation on the two dimensional torus ${\mathbb{T}}^2$, driven by an infinite dimensional Wiener process with drift in the Sobolev space $L^2(0,T;H^1\left({\mathbb{T}}^2\right))$. The solution corresponds to a stochastic Lagrangian flow in the sense of DiPerna Lions. By taking into account that the motion of a viscous incompressible fluid on the torus can be described through a suitable stochastic differential equation of the previous type, we study the inviscid limit. By establishing a large deviations principle, we show that, as the viscosity goes to zero, the Lagrangian stochastic Navier–Stokes flow approaches the Euler deterministic Lagrangian flow with an exponential rate function.     

V-cycle convergence of some multigrid methods for ill-posed problems

For ill-posed linear operator equations we consider some V-cycle multigrid approaches, that, in the framework of Bramble, Pasciak, Wang, and Xu (1991), we prove to yield level independent contraction factor estimates. Consequently, we can incorporate these multigrid operators in a full multigrid method, that, together with a discrepancy principle, is shown to act as an iterative regularization method for the underlying infinite-dimensional ill-posed problem. Numerical experiments illustrate the theoretical results.

     

带跳的随机比例微分方程的半隐式Euler数值解法(英文)

In this paper,we present the semi-implicit Euler（SIE）numerical solution for stochastic pantograph equations with jumps and prove that the SIE approximation solution converges to the exact solution in the mean-square sense under the Local Lipschitz condition.     

Relaxed Euler systems and convergence to Navier-Stokes equations

We consider the approximation of Navier-Stokes equations for a Newtonian fluid by Euler type systems with relaxation both in compressible and incompressible cases. This requires to decompose the second-order derivative terms of the velocity into first-order ones. Usual decompositions lead to approximate systems with tensor variables. We construct approximate systems with vector variables by using Hurwitz-Radon matrices. These systems are written in the form of balance laws and admit strictly convex entropies, so that they are symmetrizable hyperbolic. For smooth solutions, we prove the convergence of the approximate systems to the Navier-Stokes equations in uniform time intervals. Global-in-time convergence is also shown for the initial data near constant equilibrium states of the systems. These convergence results are established not only for the approximate systems with vector variables but also for those with tensor variables.