非对称圆环拍的特性研究

建立微小非对称圆环点冲击下拍的理论模型,通过实验验证其有效性。微小非对称圆环是分析大钟拍现象的简化模型,微小非对称圆环的每一对径向振动模态的固有频率非常接近,这对模态互相干涉下出现拍现象。最后利用拍的理论模型绘制沿圆周拍的分布图,详细分析拍的分布特性,为校正钟的拍频声提供理论依据。     

Evaluation of Euler fluxes by a high-order CFD scheme: shock instability

The construction of Euler fluxes is an important step in shock-capturing/upwind schemes. It is well known that unsuitable fluxes are responsible for many shock anomalies, such as the carbuncle phenomenon. Three kinds of flux vector splittings (FVSs) as well as three kinds of flux difference splittings (FDSs) are evaluated for the shock instability by a fifth-order weighted compact nonlinear scheme. The three FVSs are Steger–Warming splitting, van Leer splitting and kinetic flux vector splitting (KFVS). The three FDSs are Roe's splitting, advection upstream splitting method (AUSM) type splitting and Harten–Lax–van Leer (HLL) type splitting. Numerical results indicate that FVSs and high dissipative FDSs undergo a relative lower risk on the shock instability than that of low dissipative FDSs. However, none of the fluxes evaluated in the present study can entirely avoid the shock instability. Generally, the shock instability may be caused by any of the following factors: low dissipation, high Mach number, unsuitable grid distribution, large grid aspect ratio, and the relative shock-internal flow state (or position) between upstream and downstream shock waves. It comes out that the most important factor is the relative shock-internal state. If the shock-internal state is closer to the downstream state, the computation is at higher susceptibility to the shock instability. Wall-normal grid distribution has a greater influence on the shock instability than wall-azimuthal grid distribution because wall-normal grids directly impact on the shock-internal position. High shock intensity poses a high risk on the shock instability, but its influence is not as much as the shock-internal state. Large grid aspect ratio is also a source of the shock instability. Some results of a second-order scheme and a first-order scheme are also given. The comparison between the high-order scheme and the two low-order schemes indicates that high-order schemes are at a higher risk of the shock instability. Adding an entropy fix is very helpful in suppressing the shock instability for the two low-order schemes. When the high-order scheme is used, the entropy fix still works well for Roe's flux, but its effect on the Steger–Warming flux is trivial and not much clear.     

Axisymmetric three-dimensional finger solutions in Mullins-Sekerka equation

In this paper, we study the existence problem of axisymmetricthree-dimensional finger solutions of Mullins–Sekerkaequation. The finger solutions are travelling wave solutionswhose finger-shaped interfaces are moving along a certain directionat a constant speed within a cylindrical domain. The existenceof at least one axisymmetric three-dimensional finger solutionis obtained through a fixed-point argument of the Hilbert transformation.     

A strategy for choosing Gegenbauer reconstruction parameters for numerical stability

The Gegenbauer reconstruction method was first proposed in 1992, but in early studies no attempts were made to optimize the relevant parameters of this method. These parameters were allowed to grow proportionally with the number of nodes which, in many cases, resulted in exponential convergence for a selected range of the proportionality constants. Early studies also made clear that very large error bounds could be expected if these key parameters were not chosen carefully. Subsequent studies then pointed out that, although unrelated to the method’s analytically predictable domains of poor accuracy, round-off errors could also sabotage the method’s accuracy. The challenge of successfully implementing a Gegenbauer reconstruction then rests on understanding the performance trade-offs we can expect when choosing the key parameters in accordance with different objectives.In this study, we propose a new strategy for choosing optimal parameters in the Chebyshev-Gegenbauer reconstruction method, specifically to achieve numerical stability. This strategy is based on asymptotic analysis as well as minimization problems in one and two dimensions. The effectiveness of our approach, which could also be applied to a wider selection of polynomials is then illustrated with results from numerical experiments.     

The cutoff phenomenon for randomized riffle shuffles

We study the cutoff phenomenon for generalized riffle shuffles where, at each step, the deck of cards is cut into a random number of packs of multinomial sizes which are then riffled together. © 2007 Wiley Periodicals, Inc. Random Struct. Alg., 2008     

C 1,α-solutions to non-autonomous anisotropic variational problems

We establish several smoothness results for local minimizers of non-autonomous variational integrals with anisotropic growth conditions. Mathematics Subject Classification (2000) 49N60, 49N99     

Digital total variation filtering as postprocessing for Chebyshev pseudospectral methods for conservation laws

Digital total variation filtering is analyzed as a fast, robust, post-processing method for accelerating the convergence of pseudospectral approximations that have been contaminated by Gibbs oscillations. The method, which originated in image processing, can be combined with spectral filters to quickly post-process large data sets with sharp resolution of discontinuities and with exponential accuracy away from the discontinuities.     

Regularity for obstacle problems without structure conditions

This paper deals with the Lipschitz regularity of minimizers for a class of variational obstacle problems with possible occurrence of the Lavrentiev phenomenon. In order to overcome this problem, the availment of the notions of relaxed functional and Lavrentiev gap is needed. The main tool used here is a crucial Lemma which reveals to be needed because it allows us to move from the variational obstacle problem to the relaxed-functional-related one. This is fundamental in order to find the solutions’ regularity that we intended to study. We assume the same Sobolev regularity both for the gradient of the obstacle and for the coefficients.     

野渡无人舟自横现象的数值研究

为了研究``野渡无人舟自横'这一自然现象,本文采用数值仿真技术模拟研究了小船在水流作用下的运动特性,得出：顺于水流以及横于水流都属于小船平衡位置。但顺于水流属于非稳定平衡位置,受扰动后将失去平衡,并再也恢复不到原平衡位置;而横于水流属于稳定平衡位置,此时即使小船受到扰动失去平衡,最终也能在自身力矩作用下重新恢复到原平衡位置。     

氧碘化学激光器中混合现象的数值模拟

给出了基于气体动力学、化学动力学和激光物理学的氧碘化学激光器理论模型,并对氧碘化学激光器中的混合现象进行了三维数值模拟,所得结果与文献上类似模拟结果一致。