二阶非线性q-差分微分方程两点边值问题两个迭代正解的存在性*

研究了二阶非线性q-差分微分方程两点边值问题,给出了系统两个正解存在的充分条件. 与其他文献中使用的不动点定理不同,文章不仅证明了该系统正解的存在性, 而且还利用单调迭代技巧给出了逼近正解的迭代格式.     

Estimates for Divergence Velocities of Axially Moving Orthotropic Thin Plates

Some models for axially moving orthotropic thin plates are investigated analytically via methods of complex analysis to derive estimates for critical plate velocities. The linearized Kirchhoff plate theory is used, and the energy forms of steady-state models are considered with homogeneous and inhomogeneous tension profiles in the cross direction of the plate. With the help of the energy forms, some limits for the divergence velocity of the plate are found analytically. In numerical examples, the derived lower limits for the divergence velocity are analyzed for plates with small flexural rigidity.     

A stability analysis of periodic solutions to the steady forced Korteweg–de Vries–Burgers equation

The stability of periodic solutions to the steady forced Korteweg–de Vries–Burgers (fKdVB) equation is investigated here. This family of periodic solutions was identified by Hattam and Clarke (2015) using a multi-scale perturbation technique. Here, Floquet theory is applied to the governing equation. Consequently, two criteria are found that determine when the periodic solutions are stable. This analysis is then confirmed by a numerical study of the steady fKdVB equation.     

Displacement fields and self-energies of circular and polygonal dislocation loops in homogeneous and layered anisotropic solids

There are large classes of materials problems that involve the solutions of stress, displacement, and strain energy of dislocation loops in elastically anisotropic solids, including increasingly detailed investigations of the generation and evolution of irradiation induced defect clusters ranging in sizes from the micro- to meso-scopic length scales. Based on a two-dimensional Fourier transform and Stroh formalism that are ideal for homogeneous and layered anisotropic solids, we have developed robust and computationally efficient methods to calculate the displacement fields for circular and polygonal dislocation loops. Using the homogeneous nature of the Green tensor of order −1, we have shown that the displacement and stress fields of dislocation loops can be obtained by numerical quadrature of a line integral. In addition, it is shown that the sextuple integrals associated with the strain energy of loops can be represented by the product of a pre-factor containing elastic anisotropy effects and a universal term that is singular and equal to that for elastic isotropic case. Furthermore, we have found that the self-energy pre-factor of prismatic loops is identical to the effective modulus of normal contact, and the pre-factor of shear loops differs from the effective indentation modulus in shear by only a few percent. These results provide a convenient method for examining dislocation reaction energetic and efficient procedures for numerical computation of local displacements and stresses of dislocation loops, both of which play integral roles in quantitative defect analyses within combined experimental–theoretical investigations.     

土中爆炸成腔问题中确定炸点状态的一种计算反求方法

提出了一种通过给定的土中爆炸成腔毁伤效应确定炸点状态的计算反求方法。该方法将确定炸点状态的反问题转化为求解爆炸毁伤效应的计算值与给定值误差函数最小的优化问题。在反求过程中,采用基于误差减小比率技术的多项式近似模型代替土中爆炸数值分析模型,以便提高反求效率。采用Tikhonov正则化方法克服反求过程中出现的病态问题。在此基础上,引入信赖域管理策略判断当前近似模型与实际模型的逼近程度,以确定最优的反求向量。炸点状态反求结果与实验结果的对比分析表明,该方法能够有效且稳定地通过给定的毁伤效应实现炸点状态的反求,这可为炸点状态的设计提供参考。     

On the efficient application of weighted essentially nonoscillatory scheme

In this paper, the efficient application of high‐order weighted essentially nonoscillatory (WENO) reconstruction to the subsonic and transonic engineering problems is studied. On the basis of the physical considerations, two techniques are proposed to enhance the accuracy and efficiency of the WENO reconstruction. First, it is observed that the WENO scheme using characteristic variable has better accuracy and convergence speed than the scheme using primitive variable. For engineering problems with shock of moderate amplitude, on the basis of the Rankine–Hugoniot conditions, a simplified characteristic‐variable‐based WENO is developed. The simplified version significantly reduces the cost overhead without sacrificing the shock‐capturing capability. Second, in this work, it is found for viscous case that it is better to include the viscous effect. On the basis of a simple analysis, the viscous correction to the parameter ε in the WENO reconstruction is proposed. Numerical results indicate, with the proposed simplified characteristic‐variable‐based reconstruction and the viscous correction, that the nonlinear WENO interpolation is sharply activated in the region of shock jump, whereas in the shockless area, the WENO interpolation weights are tuned towards the designed optimal value for better accuracy. Compared with the original characteristic‐variable‐based WENO, the current implementation has similar accuracy and reduced cost. At the same time, compared with the primitive variable‐based WENO, better accuracy and convergence speed are obtained at marginal cost overhead. Several practical cases are calculated to demonstrate the accuracy and efficiency of the current methodology. Copyright © 2012 John Wiley & Sons, Ltd.     

Simulation of two‐phase flow–body interaction problems using direct forcing/fictitious domain–level set method

In the present paper, a direct forcing/fictitious domain (DF/FD)–level set method is proposed to simulate the twophase flow–body interaction. The DF/FD does not sacrifice accuracy and robustness by employing a discrete δ (Dirac delta) function to transfer quantities between the Eulerian nodes and Lagrangian points explicitly as the immersed boundary method. The advantages of this approach are the simple concept, the easy implementation and the utilization of original governing equation without modification. The main idea is to combine DF/FD method with the level set method in the Cartesian coordinates. We present the results of a number of test cases to illustrate the effectiveness of the proposed method for single‐phase flow–body interaction problem and the two‐phase flows with a stationary body. Eventually, the simulations of various water entry problems have been conducted to validate the capability and the accuracy of the present method on solving the twophase flow–body interaction. Consequently, the present results are found to be in good agreement with those of previous studies. Copyright © 2013 John Wiley & Sons, Ltd.     

A generic framework for time-stepping partial differential equations (PDEs): general linear methods,object-oriented implementation and application to fluid problems

Time-stepping algorithms and their implementations are a critical component within the solution of time-dependent partial differential equations (PDEs). In this article, we present a generic framework – both in terms of algorithms and implementations – that allows an almost seamless switch between various explicit, implicit and implicit–explicit (IMEX) time-stepping methods. We put particular emphasis on how to incorporate time-dependent boundary conditions, an issue that goes beyond classical ODE theory but which plays an important role in the time-stepping of the PDEs arising in computational fluid dynamics. Our algorithm is based upon J.C. Butcher's unifying concept of general linear methods that we have extended to accommodate the family of IMEX schemes that are often used in engineering practice. In the article, we discuss design considerations and present an object-oriented implementation. Finally, we illustrate the use of the framework by applications to a model problem as well as to more complex fluid problems.     

Application of a line implicit method to fully coupled system of equations for turbulent flow problems

For unstructured finite volume methods, we present a line implicit Runge–Kutta method applied as smoother in an agglomerated multigrid algorithm to significantly improve the reliability and convergence rate to approximate steady-state solutions of the Reynolds-averaged Navier–Stokes equations. To describe turbulence, we consider a one-equation Spalart–Allmaras turbulence model. The line implicit Runge–Kutta method extends a basic explicit Runge–Kutta method by a preconditioner given by an approximate derivative of the residual function. The approximate derivative is only constructed along predetermined lines which resolve anisotropies in the given grid. Therefore, the method is a canonical generalisation of point implicit methods. Numerical examples demonstrate the improvements of the line implicit Runge–Kutta when compared with explicit Runge–Kutta methods accelerated with local time stepping.     

Non-linear global dynamics of an axially moving plate

In the present study, the geometrically non-linear dynamics of an axially moving plate is examined by constructing the bifurcation diagrams of Poincaré maps for the system in the sub and supercritical regimes. The von Kármán plate theory is employed to model the system by retaining in-plane displacements and inertia. The governing equations of motion of this gyroscopic system are obtained based on an energy method by means of the Lagrange equations which yields a set of second-order non-linear ordinary differential equations with coupled terms. A change of variables is employed to transform this set into a set of first-order non-linear ordinary differential equations. The resulting equations are solved using direct time integration, yielding time-varying generalized coordinates for the in-plane and out-of-plane motions. From these time histories, the bifurcation diagrams of Poincaré maps, phase-plane portraits, and Poincaré sections are constructed at points of interest in the parameter space for both the axial speed regimes.