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.

     

Simulation of viscous water column collapse using adapting hierarchical grids

An adaptive hierarchical grid‐based method for predicting complex free surface flows is used to simulate collapse of a water column. Adapting quadtree grids are combined with a high‐resolution interface‐capturing approach and pressure‐based coupling of the Navier–Stokes equations. The Navier–Stokes flow solution scheme is verified for simulation of flow in a lid‐driven cavity at Re=1000. Two approaches to the coupling of the Navier–Stokes equations are investigated as are alternative face velocity and hanging node interpolations. Collapse of a water column as well as collapse of a water column and its subsequent interaction with an obstacle are simulated. The calculations are made on uniform and adapting quadtree grids, and the accuracy of the quadtree calculations is shown to be the same as those made on the equivalent uniform grids. Results are in excellent agreement with experimental and other numerical data. A sharp interface is maintained at the free surface. The new adapting quadtree‐based method achieves a considerable saving in the size of the computational grid and CPU time in comparison with calculations made on equivalent uniform grids. Copyright © 2005 John Wiley & Sons, Ltd.     

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.

     

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.     

Numerical simulation of thermally and chemically nonequilibrium flows and heat transfer in underexpanded induction plasmatron jets

The results of numerical simulation are presented for thermally and chemically nonequilibrium air plasma flows in a plasmatron discharge channel and underexpanded dissociated and partially ionized air jets flowing past a cylindrical model with a blunt leading edge and cooled copper surface under the experimental conditions realized in a VGU-4 100 kW induction plasmatron (Institute for Problems in Mechanics of the Russian Academy of Sciences) (see, for example, [1, 2]). The nonequilibrium excitation of the vibrational degrees of freedom of the molecules in the modal approximation and the difference between the electron and translational heavy-particle temperatures are taken into account in the calculations. The calculated data on the heat transfer and pressure at the stagnation point are compared with the results obtained within the framework of the thermally equilibrium model. Comparison with the experimental data obtained in the Institute for Problem in Mechanics of the Russian Academy of Sciences (Laboratory for interaction between plasma and radiation and materials) and kindly provided for comparison purposes gives satisfactory agreement.     

基于椭圆型方程的扭叶片造型方法的研究

本文以偏微分方程造型为基础,提出了一种基于椭圆型方程的扭叶片三维型面直接设计方法,详细推导了叶型曲面函数,给出了型面方程的求解及其前后缘修正。该方法具有设计叶型曲面自然光顺,设计参数少且各参数具有明显的几何意义,叶型曲面调整方便,利于采用非数值优化算法对其进行气动优化等优点。文中给出了设计实例,并通过数值实验分析了所设计叶片型面的流动特性。分析结果表明设计叶片具有良好的气动性能,同时也证明了本文提出的基于椭圆型方程的扭叶片三维型面设计方法的可靠性和实用性。     

Energy decay of solutions to Maxwell's equations with conductivity and polarization

In this paper, Maxwell's equations involving generally nonlinear polarization and field-dependent currents are studied. The main objective is the asymptotic behavior of the solution for t→∞ if no damping term occurs in the equation governing the polarization field.     

Besov-Morrey spaces: Function space theory and applications to non-linear PDE

This paper is devoted to the analysis of function spaces modeled on Besov spaces and their applications to non-linear partial differential equations, with emphasis on the incompressible, isotropic Navier-Stokes system and semi-linear heat equations. Specifically, we consider the class, introduced by Hideo Kozono and Masao Yamazaki, of Besov spaces based on Morrey spaces, which we call Besov-Morrey or BM spaces. We obtain equivalent representations in terms of the Weierstrass semigroup and wavelets, and various embeddings in classical spaces. We then establish pseudo-differential and para-differential estimates. Our results cover non-regular and exotic symbols. Although the heat semigroup is not strongly continuous on Morrey spaces, we show that its action defines an equivalent norm. In particular, homogeneous BM spaces belong to a larger class constructed by Grzegorz Karch to analyze scaling in parabolic equations. We compare Karch's results with those of Kozono and Yamazaki and generalize them by obtaining short-time existence and uniqueness of solutions for arbitrary data with subcritical regularity. We exploit pseudo-differential calculus to extend the analysis to compact, smooth, boundaryless, Riemannian manifolds. BM spaces are defined by means of partitions of unity and coordinate patches, and intrinsically in terms of functions of the Laplace operator.

     

Oscillation theorems for differential equations with a nonlinear damping term

In this paper, we are concerned with a class of nonlinear second-order differential equations with a nonlinear damping term. Passage to more general class of equations allows us to remove a restrictive condition usually imposed on the nonlinearity, and, as a consequence, our results apply to wider classes of nonlinear differential equations. Two illustrative examples are considered.     

使用非单调技术的不精确预条件牛顿类方法解非线性方程组

本文提供了预条件不精确牛顿型方法结合非单调技术解光滑的非线性方程组．在合理的条件下证明了算法的整体收敛性．进一步，基于预条件收敛的性质，获得了算法的局部收敛速率，并指出如何选择势序列保证预条件不精确牛顿型的算法局部超线性收敛速率．