Adaptive and anisotropic mesh strategy for thin shell problems. Case of inhibited parabolic shells

The aim of this paper is to show the reliability of an adaptive and anisotropic mesh procedure for thin shell problems. We consider singular perturbation problems only for parabolic shells whose behavior is described by the Koiter model. The corresponding system of equations, which depends on the relative thickness ε of the shell, is elliptic except at the limit for ε = 0 where it is parabolic. In a first part of this paper, we study theoretically the phenomena of internal layers appearing during the singular perturbation process, when the loading is somewhat singular. These layers have very different structures either they are along or across the asymptotic lines of the middle surface of the shell. In a second part, numerical computations are performed using a finite element software coupled with an adaptive anisotropic mesh generator. This technique enables to approach accurately the singularities and the layers predicted by the theory especially for very small values of the thickness. The efficiency of such a procedure in comparison with uniform meshes is put in a prominent position.     

Existence, number, and stability of limit cycles in weakly dissipative, strongly nonlinear oscillators

Oscillators control many functions of electronic devices, but are subject to uncontrollable perturbations induced by the environment. As a consequence, the influence of perturbations on oscillators is a question of both theoretical and practical importance. In this paper, a method based on Abelian integrals is applied to determine the emergence of limit cycles from centers, in strongly nonlinear oscillators subject to weak dissipative perturbations. It is shown how Abelian integrals can be used to determine which terms of the perturbation are influent. An upper bound to the number of limit cycles is given as a function of the degree of a polynomial perturbation, and the stability of the emerging limit cycles is discussed. Formulas to determine numerically the exact number of limit cycles, their stability, shape and position are given.     

Subcritical spatial transition of swept Hiemenz flow

It is known from experimental investigations that the leading-edge boundary layer of a swept wing exhibits transition to turbulence at subcritical Reynolds numbers, i.e. at Reynolds numbers which lie below the critical Reynolds number predicted by linear stability theory. In the present work, we investigate this subcritical transition process by direct numerical simulations of a swept Hiemenz flow in a spatial setting. The laminar base flow is perturbed upstream by a pair of stationary counter-rotating vortex-like disturbances. This perturbation generates high- and low-speed streaks by a non-modal growth mechanism. Further downstream, these streaky structures exhibit a strong instability to secondary perturbations which leads to a breakdown to turbulence.The observed transition mechanism has strong similarities to by-pass transition mechanisms found for two-dimensional boundary layers. It can be shown that transition strongly depends on the amplitude of the primary perturbation as well as on the frequency of the secondary perturbation.     

倾斜界面耦合弹性层中的渡越辐射能

物理学中,摄动源在非均匀介质中或非均匀介质附近匀速直线运动所产生的能量辐射现象称为渡越辐射.列车沿轨道运行,由轮轨接触产生的弹性波在非均匀轨道和基础中传播将发生渡越辐射,而轨道和基础的非均匀性集中体现在不同轨道基础之间的过渡段(如路桥过渡段、桥隧过渡段或有砟-无砟轨道过渡段).为研究车致弹性波在过渡段中引发的渡越辐射现象,本文以典型高速铁路路桥过渡段结构形式为依据,建立了二维平面应力渡越辐射能计算模型.其中,两个材料参数不同的半无限弹性层由一倾斜界面耦合,底端固定,上表面自由,一个集中载荷在自由表面上匀速运动.界面两侧弹性体中的波动方程均分解为本征场、自由场两个部分分别求解,其中自由场波动方程采用分离变量法数值求解.通过模型求解得到了不同载荷移动速度和界面倾斜角度条件下的渡越辐射能及界面附近应变能密度.结果表明,渡越辐射能的大小随载荷移动速度增大单调非线性增大,移动载荷速度达到刚度较大一侧介质表面波速的74%时产生的渡越辐射能就将超过载荷本身激发的本征场应变能;界面倾斜角度越大,即两侧介质刚度过渡距离越短,渡越辐射能与本征场应变能比值越大.     

Characteristics of singular entities of simple cell mappings

The method of cell-to-cell mapping has the potential to be a very effective and general method of global analysis for strongly non-linear systems. However, simple cell mappings being integer mappings, most of the classical methods of analysis based upon continuity and differentiability of the mapping are no longer applicable and, therefore, new notions need be introduced. In [6] the concept of singular multiplets is introduced for the cell functions associated with the cell mappings. In this paper we study the characteristics of these singular entities by examining the mapping properties of the cells in the singular entities and in their neighborhoods. The key tool used in classifying the mapping properties is the limit set of the mapping process of a cell. The work represents a continuing effort to develop the method of cell-to-cell mapping as a tool of global analysis and to provide the method with a sound and appropriate structure.     

Unsteady Flow to Wells Partially Penetrating in Two-Layered Aquifers

The unsteady flow of water to well in a layered aquifer with an interlayer flow is examined in this paper. The system studied comprises an aquifer consisting of two productive layers of finite thickness that are in a perfect hydraulic contact and a well which partially penetrates in one of the layers. Each layer is assumed to be homogeneous and isotropic and the water contained in the aquifer is of identical physical properties and small compressibility.The analytical solutions are derived for the case when the system is characterised by equality of hydraulic diffusivity of layers. These solutions give the results accurate enough for practical applications and allow to estimate the effects of partial penetration and contrast in parameters of formation on distribution of the flow potential both at large distances from wells and at the wells.The obtained solutions also provide a basis for predicting the interlayer flow arising from the performance of a pumping well as well as between a pair of wells which have an open interval located in different layers of the aquifer. Two cases have been analysed: (1) a pumping pair of wells which is used to prevent water inflow to the productive well, and (2) a pumping-injection pair of wells providing the demanded mixing of water from adjacent layers. Some examples of flow patterns and specially computed diagrams are given to illustrate the operation of such systems.     

Phantom ducks and models of excitability

We study a system which may serve as a prototype model for excitability. Its long-term behavior is governed by a globally attracting fixed point, but it has the interesting property that an appropriate small perturbation of the equilibrium situation elicits a large (transient) response. Using singular perturbation techniques, the threshold for excitation is determined. Asymptotic expansions for the threshold are given and their validity is proved. An associated system is studied, which yields further quantitative information on the position of the threshold.     

Dynamics of a quantized circle homeomorphism with quasiperiodic perturbation

We introduce the notion of a quantized circle homeomorphism that is a discontinuous mapping of an interval shift, which is widely used in modern digital radio electronics. For a two-dimensional dynamical system given by a triangular mapping, which is a quantized circle homeomorphism with quasiperiodic perturbation, we prove, under some assumptions, that there exist an invariant absorbing belt and a repulsive contour, study properties of these structures, and get estimates for their sizes. To make the exposition complete, we, first, study the corresponding problems for three less complicated systems, namely, a proper circle homeomorphism, a proper circle homeomorphism with quasiperiodic perturbation, and a quantized circle homeomorphism without perturbation.     

Non-unique turbulent boundary layer flows having a moderately large velocity defect: a rational extension of the classical asymptotic theory

The classical analysis of turbulent boundary layers in the limit of large Reynolds number Re is characterised by an asymptotically small velocity defect with respect to the external irrotational flow. As an extension of the classical theory, it is shown in the present work that the defect may become moderately large and, in the most general case, independent of Re but still remain small compared to the external streamwise velocity for non-zero pressure gradient boundary layers. That wake-type flow turns out to be characterised by large values of the Rotta–Clauser parameter, serving as an appropriate measure for the defect and hence as a second perturbation parameter besides Re. Most important, it is demonstrated that also this case can be addressed by rigorous asymptotic analysis, which is essentially independent of the choice of a specific Reynolds stress closure. As a salient result of this procedure, transition from the classical small defect to a pronounced wake flow is found to be accompanied by quasi-equilibrium flow, described by a distinguished limit that involves the wall shear stress. This situation is associated with double-valued solutions of the boundary layer equations and an unconventional weak Re-dependence of the external bulk flow—a phenomenon seen to agree well with previous semi-empirical studies and early experimental observations. Numerical computations of the boundary layer flow for various values of Re reproduce these analytical findings with satisfactory agreement.     

温盐双扩散系统的多平衡稳定性

采用数值试验的方法研究了温盐双扩散系统中的多平衡稳定性在系统中引入非线性状态方程、温通量边界条件和变强迫因子在数值模拟中引入边界拟会坐标系,并按照温盐浮力比λ大于1,等于1和小于1将流场分为温度控制场、温盐均势场和盐度控制场得到了四种稳定平衡状态,特别是,大强迫因子和强非线性影响下的四胞环流对称状态．在对称的定常强迫边条件作用下的流场、温度场和盐度场均为对称分布引入温通量边条件后,对称流场在小干扰下失去对称性,发展为非对称分布,并因为干扰位置的不同而造成不同的非对称性,出现了从对称状态到非对称状态的音叉分叉,得到了a-Ras的稳定曲线及分叉图,并解释了分叉的本质     

Laminar flow through a channel with uniformly porous walls of different permeability

Summary The steady laminar flow of a viscous incompressible fluid through a two-dimensional channel, having fluid sucked or injected with different velocities through its uniformly porous parallel walls is considered. A solution for small suction Reynolds number has been given by the authors in a previous paper. The purpose of this paper is to present a solution valid for large Reynolds numbers for the cases of (i) suction at both walls, and (ii) suction at one wall and injection at the other. A technique of matching outer and inner expansions is used to obtain an asymptotic solution for both of these cases. Further a perturbation solution for the case of suction at one wall and injection at the other is obtained by choosing the difference between two wall velocities as the perturbation parameter. Both asymptotic and perturbation solutions are confirmed by exact numerical solutions. As expected, the resulting solutions show the presence of the usual suction boundary layers in both types of flow considered in this paper.     

Bifurcation of limit cycles from a Liénard system with a heteroclinic loop connecting two nilpotent saddles

In this article, we study the limit cycles bifurcated from a Liénard system with a heteroclinic loop connecting two nilpotent saddles. We apply expansion theory of a first-order Melnikov function to investigate the number of limit cycles near the heteroclinic loop and the center, and by some perturbation theory we find 3 limit cycles with 7 different distributions. Last, the least upper bound of the number of limit cycles bifurcated from the annulus is given by an algebraic criterion developed in J. Differ. Equ. 251, 1656–1669 (2011).     

Dissipativity of a class of cellular neural networks with proportional delays

In this paper, the problem of dissipativity is investigated for cellular neural networks with proportional delays. Without assuming monotonicity, differentiability, and boundedness of activation functions, two new delay-independent criteria for checking the dissipativity of the addressed neural networks are established by using inner product properties and matrix theory. Two examples and their simulation results are given to show the effectiveness and less conservatism of the proposed criteria.     

A METHOD TO CORRECT CATALOG ORBIT DETERMINATION VELOCITY 1)

对于大批量空间目标,监测资源有限,测轨数据稀疏,导致编目定轨结果误差较大。本文分析了不同轨道类型的编目轨道预报误差演化特性,分析结果表明,轨道预报误差主要分布在沿迹方向,且主要是由于大气阻力摄动模型误差和初始径向速度误差引起的。进一步的理论推导显示,在忽略初始位置误差的假设条件下,轨道初值误差引起的预报位置误差前后具有近似对称特性,利用仿真数据,验证了近似对称特性的正确性。基于上述分析,本文提出了一种校准编目定轨初始速度的方法,即通过减小向前预报的位置与已知位置的偏差来迭代地校准定轨结果的速度,从而提高目标向后预报的轨道精度。利用实际轨道数据的试验结果表明,该方法对于采用稀疏数据的近圆轨道目标定轨结果修正效果明显,可以有效改进自主编目定轨结果的精度,对提高大批量空间目标的编目管理能力具有重要价值。     

Connection between Volterra series and perturbation method in nonlinear systems analyses

This paper is concerned with the connection between the Volterra series and the regular perturbation method in nonlinear systems analyses. It is revealed for the first time that, for a forced polynomial nonlinear system, if its derived linear system is a damped dissipative system, the steady response obtained through the regular perturbation method is exactly identical to the response given by the Volterra series. On the other hand, if the derived linear system is an undamped conservative system, then the Volterra series is incapable of modeling the forced polynomial nonlinear system. Numerical examples are further presented to illustrate these points. The results provide a new criterion for quickly judging whether the Volterra series is applicable for modeling a given polynomial nonlinear system.     

On the persistence of regular reflection

It is well accepted that the persistence of regular reflection of a shock wave off a wedge beyond the ideal theoretical prediction is due to viscous and thermal boundary layers induced behind the reflection point. Experiments have been done by reflecting two shock waves of equal strength off each other so that the plane of symmetry between them becomes an ideal inviscid and adiabatic reflection plane thereby experimentally mimicking the assumptions of the theory. There is one definitive experiment done at a wall angle of 40° using a bifurcated shock tube that indicates that the actual transition angle is the theoretical detachment condition. This paper extends these results to two cases near limiting conditions; one at a very low incidence shock Mach number and one at a wall angle very close to the theoretical transition limit. The first confirms the reasons for the von Neumann Paradox but cannot discriminate between sonic and detachment conditions, but is within about 0.5% of them, and the second shows transition much closer to the sonic than the detachment condition but with both within the experimental error bounds. In both cases, the results are notably different from transition conditions off a wedge and confirm the effects of transport properties being the cause of persistence of regular reflection.     

Shooting method in singular perturbation problem of ordinary differential equations with boundary layers

By using the adaptive steplength integration scheme with a shooting technique, a rather difficult singular perturbation problem of ordinary differential equations with boundary layers can be calculated effectively. Computing examples are given in this paper which show the convergence within one iteration of the method in the case of a linear problem, the efficiency of the method for many boundary layers and turning points, especially the convenience in calculating multiple solutions. A comparison with traditional difference method is given at the end of this paper.     

可压缩自由剪切流混合转捩大涡模拟

针对湍流气动光学效应与冲压发动机气体混合机理问题,开展了可压缩混合层流动空间模式大涡模拟和时间模式直接数值模拟研究.通过对流场(包含亚/亚混合、超/亚混合两种情况)失稳、转捩直至完全湍流的空间发展过程的研究表明,对流Mach数0.4状态下流场失稳以二维最不稳定扰动为主;非线性发展中,基频涡对并及展向涡撕裂主控流动转捩,流场发生混合转捩;转捩后脉动流场基本达到局部各向同性,此时,湍流Mach数低于0.3,流动压缩性可近似忽略.     

Design Sensitivity Analysis in Structural Mechanics. I. Static Response Variations

ABSTRACT

The dependence of the solution of boundary-value problems of structural mechanics on design variables that specify material properties and distribution is characterized. Prototype problems treated include beams, plates, and plane elastic solids. Symmetry and positive definiteness properties of the elliptic differential operators that govern system response are used to show that their inverses, hence the displacement fields, are Frechet differentiable with respect to design variables. Formulas for the derivatives are given and used to obtain computable formulas for design sensitivity coefficients (first variation) of integrals that arise in optimal design formulations. The results establish an extension of the concept of “well-posed” problems of structural mechanics to include continuity (in fact, differentiability” of static structural response with respect to distributed design variables and design parameters     

Stability characteristics of transition layer solutions

The multiple existences and their stability properties of stationary solutions with a single transition layer in some scalar reaction-diffusion equation are shown. Each solution is constructed by using classical singular perturbation methods and its stability property is determined by a simple algebraic quantity, say index, appearing in the construction of a solution.