Multidimensional transonic shocks and free boundary problems for nonlinear equations of mixed type

We establish the existence and stability of multidimensional transonic shocks for the Euler equations for steady potential compressible fluids. The Euler equations, consisting of the conservation law of mass and the Bernoulli law for the velocity, can be written as a second-order, nonlinear equation of mixed elliptic-hyperbolic type for the velocity potential. The transonic shock problem can be formulated into the following free boundary problem: The free boundary is the location of the transonic shock which divides the two regions of smooth flow, and the equation is hyperbolic in the upstream region where the smooth perturbed flow is supersonic. We develop a nonlinear approach to deal with such a free boundary problem in order to solve the transonic shock problem. Our results indicate that there exists a unique solution of the free boundary problem such that the equation is always elliptic in the downstream region and the free boundary is smooth, provided that the hyperbolic phase is close to a uniform flow. We prove that the free boundary is stable under the steady perturbation of the hyperbolic phase. We also establish the existence and stability of multidimensional transonic shocks near spherical or circular transonic shocks.

     

Convergence analysis of a modified inexact implicit method for general mixed monotone variational inequalities

We consider a useful modification of the inexact implicit method with a variable parameter in Wang et al. J Optim Theory 111: 431–443 (2001) for generalized mixed monotone variational inequalities. One of the contributions of the proposed method in this paper is that the restrictions imposed on the variable parameter are weaker than the ones in Wang et al. J Optim Theory 111: 431–443 (2001). Another contribution is that we establish a sufficient and necessary condition for the convergence of the proposed method to a solution of the general mixed monotone variational inequality.     

Substructuring preconditioners for saddle-point problems arising from Maxwell's equations in three dimensions

This paper is concerned with the saddle-point problems arising from edge element discretizations of Maxwell's equations in a general three dimensional nonconvex polyhedral domain. A new augmented technique is first introduced to transform the problems into equivalent augmented saddle-point systems so that they can be solved by some existing preconditioned iterative methods. Then some substructuring preconditioners are proposed, with very simple coarse solvers, for the augmented saddle-point systems. With the preconditioners, the condition numbers of the preconditioned systems are nearly optimal; namely, they grow only as the logarithm of the ratio between the subdomain diameter and the finite element mesh size.

     

Variational Derivative for Buckling Loads of Beams and Frames

The finite-difference method is a numerical technique for obtaining approximate solutions to differential equations. The main objective of the present study is to give a new aspect to the finite-difference method by using a variational derivative. By applying this formulation, accurate values of the buckling loads of beams and frames with various end supports are obtained. The performance of this formulation is verified by comparison with numerical examples in the literature __________ Published in Prikladnaya Mekhanika, Vol. 41, No. 7, pp. 139–144, July 2005.     

A WKB analysis of the buckling condition for a cylindrical shell of arbitrary thickness subjected to an external pressure

We consider the plane-strain buckling of a cylindrical shellof arbitrary thickness which is made of a Varga material andis subjected to an external hydrostatic pressure on its outersurface. The WKB method is used to solve the eigenvalue problemthat results from the linear bifurcation analysis. We show thatthe circular cross-section buckles into a non-circular shapeat a value of µ₁ which depends on A₁/A₂ and a mode number,where A₁ and A₂ are the undeformed inner and outer radii, andµ₁ is the ratio of the deformed inner radius to A₁. Inthe large mode number limit, we find that the dependence ofµ₁ on A₁/A₂ has a boundary layer structure: it is constantover almost the entire region of 0 < A₁/A₂ < 1 and decreasessharply from this constant value to unity as A₁/A₂ tends tounity. Our asymptotic results for A₁ – 1 = O(1) and A₁– 1 = O(1/n) are shown to agree with the numerical resultsobtained by using the compound matrix method.     

Absorbing boundary conditions for reaction-diffusion equations

We treat here of the question of absorbing boundary conditionsfor nonlinear diffusion equations. We use the conditions designedfor the linear equation, we prove them to be well posed forthe nonlinear problem, and through numerical experiments thatthey are well suited for reaction–diffusion equations.     

CAP-子群与p-幂零性

设G是有限群。H是G的一个正规子群。P是|H|的一个素因子。P是H的一个Sylow p-子群。若下列条件之一满足。则H是p-幂零：1)P的极大子群都是G的CAP-子群且(|H|，P-1)=1；2)P的二次极大子群都是G的CAP-子群且(|H|，P^2-1)=1。特别地，若进一步假设G／H为P-幂零且(|G|，P-1)=1或(|G|，P^2-1)=1。则G自身亦为P-幂零。     

CFD在燃煤细粒子凝聚过程中的应用

燃煤产生的细粒子富集了大量的有毒痕量元素,对大气环境和人类健康造成严重危害。本文针对燃煤细粒子的形成过程,在计算流体力学(CFD)的基础上,结合气溶胶动力学理论,模拟了烟气中细粒子在圆柱形流场中由于碰撞而凝聚的过程,并通过计算结果分析细粒子颗粒特征参量(速度、质量、直径和颗粒数目)随流场的变化趋势,为深入研究燃煤细粒子的形成演化机理提供理论基础。     

On computation of the greatest common divisor of several polynomials over a finite field

We propose a probabilistic algorithm to reduce computing the greatest common divisor of m polynomials over a finite field (which requires computing m−1 pairwise greatest common divisors) to computing the greatest common divisor of two polynomials over the same field.     

On the degree of convergence of lemniscates in finite connected domains

For an arbitrary bounded closed set E in the complex plane with complement Ω of finite connectivity, we study the degree of convergence of the lemniscates in Ω.