Stokes型积分微分方程的质量集中各向异性非协调有限元分析

本文将Crouzeix-Raviart型非协调三角形元应用到发展型Stokes积分微分方程,给出了其质量集中非协调有限元逼近格式.在各向异性网格下,导出了速度的L2模和能量模及压力的L2模的误差估计.     

非定常Navier-Stokes方程的质量集中各向异性非协调有限元逼近

该文将一个低阶Crouzeix-Raviart型非协调三角形元应用到非定常Navier-Stokes方程,给出了其质量集中有限元逼近格式.在不需要传统Ritz-Volterra投影下,通过引入两个辅助有限元空间对边界进行估计的技巧,在各向异性网格下导出了速度的L~2模和能量模及压力的L~2模的误差估计.     

双曲型方程的非协调变网格有限元方法

采用变网格的思想讨论了双曲型方程在各向异性网格下的Crouzeix-Raviart型非协调有限元逼近.在不需要引入传统分析中Riesz投影的情况下,得到了相应最优误差估计.     

变度量各向异性网格生成算法和各向异性椭圆偏微分方程匹配网格

对于各向异性问题而言,在进行数值求解时,一个相匹配的各向异性网格至关重要.本文针对二阶椭圆偏微分方程,当其系数矩阵为各向异性时,给出了系数矩阵的逆作为相匹配的各向异性度量,在此度量下生成的各向异性网格即为匹配网格.本文还提出了新的变度量各向异性网格生成算法,该算法通过在各向异性背景网格上结合各向异性Delaunay准则和基于力平衡的结点优化函数生成高质量的各向异性三角化网格.通过数值算例可知,在匹配的各向异性三角化网格上,不仅有条件数小的代数离散系统,还有高精度的数值解,更甚者,在匹配网格结点上的l2范数误差有超收敛现象.     

Lagrange中心型守恒格式

提出了Lagrange中心型守恒气体动力学格式．引入了当前时刻子网格密度与当前时刻网格声速产生的网格分片常数压力．初始网格密度乘以初始子网格体积得到子网格质量,这些子网格质量除以当前时刻子网格体积得到当前时刻子网格密度．应用网格分片常数压力,构造了满足动量守恒、总能量守恒的Lagrange中心型守恒气体动力学格式,格点速度以与网格面的数值通量相容的方式计算．对Saltzman活塞问题等进行了数值模拟,数值结果显示Lagrange中心型守恒气体动力学格式的有效性和精确性．     

Stokes问题在各向异性网格下的Bernadi-Raugel有限元逼近

在各向异性网格下,得到了Stokes问题著名的Bernadi-Raugel混合有限元格式的超逼近性质,而且通过构造插值后处理算子得到了关于速度的超收敛结果.     

特征值问题的Lagrange型各向异性有限元方法

在各向异性网格下首先研究了二阶椭圆特征值问题算子谱逼近的若干抽象结果.然后将这些结果具体应用于线性和双线性Lagrange型协调有限元,得到了与传统有限元网格剖分下相同的最优误差估计,从而拓宽了已有的成果.     

各向异性非内积型网格下旋转Q1元的收敛性

主要研究非协调旋转Q₁元在各向异性非内积型网格,即各向异性平行四边形网格下的收敛性．该元的插值误差在各向异性网格剖分时是发散的,但是我们证明了它对于二阶椭圆问题的各向异性收敛性．为了克服该元插值误差的这个缺陷,构造了一个适当的算子来证明它的逼近性质．利用一定的分析技巧,在各向异性平行四边形网格剖分下,得到了旋转Q₁元的最优逼近误差和相容误差．该结果可以引发人们继续分析一些不满足各向异性插值性质单元的收敛性．文章的最后用一些算例来验证了理论分析的正确性．     

带弱奇异核非线性积分微分方程的有限元分析

讨论了带弱奇异核的非线性抛物积分微分方程的Hermite型各向异性矩形元逼近.在各向异性网格下导出了关于Riesz投影的L~2和H~1模的误差估计.在半离散和向后欧拉全离散格式下,基于Riesz投影的性质并利用平均值技巧,分别得到了L~2模意义下的最优误差估计.     

广义神经传播方程的非协调变网格有限元方法

本文将Crouzeix-Raviart型各向异性非协调线性三角形元应用到广义神经传播方程,建立了其Crank-Nicolson变网格逼近格式.同时,直接利用插值技巧和单元的特殊性质给出了相应的收敛性分析和最优误差估计.     

Multi‐dimensional Lotka–Volterra systems for carcinogenesis mutations

In the paper we consider three classes of models describing carcinogenesis mutations. Every considered model is described by the system of (n+1) equations, and in each class three models are studied: the first is expressed as a system of ordinary differential equations (ODEs), the second—as a system of reaction–diffusion equations (RDEs) with the same kinetics as the first one and with the Neumann boundary conditions, while the third is also described by the system of RDEs but with the Dirichlet boundary conditions. The models are formulated on the basis of the Lotka–Volterra systems (food chains and competition systems) and in the case of RDEs the linear diffusion is considered. The differences between studied classes of models are expressed by the kinetic functions, namely by the form of kinetic function for the last variable, which reflects the dynamics of malignant cells (that is the last stage of mutations). In the first class the models are described by the typical food chain with favourable unbounded environment for the last stage, in the second one—the last equation expresses competition between the pre‐malignant and malignant cells and the environment is also unbounded, while for the third one—it is expressed by predation term but the environment is unfavourable. The properties of the systems in each class are studied and compared. It occurs that the behaviour of solutions to the systems of ODEs and RDEs with the Neumann boundary conditions is similar in each class; i.e. it does not depend on diffusion coefficients, but strongly depends on the class of models. On the other hand, in the case of the Dirichlet boundary conditions this behaviour is related to the magnitude of diffusion coefficients. For sufficiently large diffusion coefficients it is similar independently of the class of models, i.e. the trivial solution that is unstable for zero diffusion gains stability. Copyright © 2009 John Wiley & Sons, Ltd.     

Modelling and simulation of micro-well formation

Physico-chemical processes on the micro-scale require new modelling concepts because some effects become dominating that are negligible for macroscopic systems. This is illustrated by a new method for the production of micro-wells based on the placement of a small drop of toluene on a plate of polystyrene. After droplet evaporation, a micro-well is left. A mathematical model has been developed to understand the elementary processes of the micro-well formation. The model accounts for: (1) growth of the drop on the substrate, (2) evaporation process of the solvent, (3) dissolution of the substrate, (4) flow rate in the evaporating drop caused by the pinning effect, including the vertical velocity profile, and (5) increase in the concentration of dissolved material followed by precipitation. In the modelling and simulation process, it could be shown that the method of drop production also has a significant influence on the shape of the micro-wells.     

Forces of rails for electromagnetic railguns

In this paper, using electromagnetic theory, the equations for the total magnetic force on the rails are deduced. Besides this, the equations for the changes of the magnetic force along with time and the running position of the armature under jointed pulse operating current are presented. Then, the equations for changes of the bending stress and shear stress in the rail along with time and the running position of the armature are given. Changes of the magnetic force distribution between the two rails along with the running position of the armature and other parameters are investigated. The total electromagnetic forces on the rail under jointed pulse operating current are analyzed. Distribution and changes of the bending stresses in the rail are studied. A number of results are obtained. The results are useful for design and application of the railgun system.     

基于韧性度的低轨卫星通信网络抗毁性度量及优化

低轨卫星通信网络的抗毁性是描述网络安全可靠的有效工具,在网络体系结构设计和路由策略等领域得到了广泛的应用。根据低轨卫星通信网络中卫星在轨道平面内移动,需要不断进行切换的特点,从建立抗毁性测度模型以及网络抗毁性优化两个角度来评估和提高网络抗毁性,提出一种基于韧性度的低轨卫星通信网络抗毁性度量方法。通过对移动模型以及切换模型的结构分析,对每种结构以一定概率出现的低轨卫星通信网络,应用韧性度函数,求得网络在某个时刻及某一段时间段内的抗毁性,并针对切换模型的不足之处进行优化,用赋权韧性度来体现优化的效果,得到了优化后的网络抗毁性。以铱星系统为应用实例进行仿真,结果表明:任意时刻网络的抗毁性跟拓扑结构的韧性度值有关,并且是一种线性关系,即随着韧性度的增加,其抗毁性也增加。通过对铱星通信系统切换模型的优化,网络的抗毁性与平均抗毁性都得到了提升,说明本文所构建模型的有效性和实用性。     

Piecewise parametric polynomial fuzzy sets

We present a scheme for tractable parametric representation of fuzzy set membership functions based on the use of a recursive monotonic hierarchy that yields different polynomial functions with different orders. Polynomials of the first order were found to be simple bivalent sets, while the second order polynomials represent the typical saw shape triangles. Higher order polynomials present more diverse membership shapes. The approach demonstrates an enhanced method to manage and fit the profile of membership functions through the access to the polynomials order, the number and the multiplicity of anchor points as wells as the uniformity and periodicity features used in the approach. These parameters provide an interesting means to assist in fitting a fuzzy controller according to system requirements. Besides, the polynomial fuzzy sets have tractable characteristics concerning the continuity and differentiability that depend on the order of the polynomials. Higher order polynomials can be differentiated as many times as the order of the polynomial less the multiplicity of the anchor points. An algorithmic optimization approach using the steepest descent method is introduced for fuzzy controller tuning. It was shown that the controller can be optimized to model a certain output within small number of iterations and very small error margins. The mathematics generated by the approach is consistent and can be simply generalized to standard applications. The recursive propagation was noticed for its clarity and ease of calculations. Further, the degree of association between the sets is not limited to the neighbors as in traditional applications; instead, it may extend beyond.Such approach can be useful in dynamic fuzzy sets for adaptive modeling in view of the fact that the shape parameters can be easily altered to get different profiles while keeping the math unchanged. Hypothetically, any shape of membership functions under the partition of unity constraint can be produced. The significance of the mentioned characteristics of such modeling can be observed in the field of combinatorial and continuous parameter optimization, automated tuning, optimal fuzzy control, fuzzy-neural control, membership function fitting, adaptive modeling, and many other fields that require customized as well as standard fuzzy membership functions. Experimental work of different scenarios with diverse fuzzy rules and polynomial sets has been conducted to verify and validate our results.     

Solution of two-dimensional boundary-value problems of the theory of thin composite shells

Discrete analogues of the boundary-value problems of a two-dimensional refined theory of anisotropic shells taking into account the transverse shear deformation are presented. The systems of resolving equations in the general form are obtained for arbitrary nonshallow shells of variable curvature whose coordinate lines of the reduction surface may not coincide with the lines of principal curvatures. The algebraic problems of determining the stress-strain state in shells made of composite materials with stress concentrators under various kinds of loads are obtained as particular cases of the schemes presented. The results of calculating the stress concentration near a nonsmall circular hole in a transversely isotropic nonshallow spherical shell under internal pressure are presented. The dependences of stress concentration factors on the hole dimension and on a change in the shear stiffness of the shells are studied. A comparison between the calculation results obtained within the framework of the theories of shallow and nonshallow shells is given.Presented at the 11th International Conference on the Mechanics of Composite Materials (Riga, June 11–15, 2000).Timoshenko Institute of Mechanics, Ukranian National Academy of Sciences, Kiev, Ukraine. Translated from Mekhanika Kompozitnykh Materialov, Vol. 36, No. 4, pp. 465–472, July–August, 2000.     

The efficiency of remanufacturing in a dynamic input–output model

For industrialised economy of ourdays, remanufacturing represents perhaps the largest unexploited resource and opportunity for realising a greater growth of the economy in an environmental-conscious manner. The aim of this paper is to investigate of the impact of remanufacturing in the economy from an economic-efficiency point of view. In static context this phenomenon was analysed in the literature. We use the multi-sector input–output framework in a dynamic context to study intra-period relationships of the sectors of economy. We extend the classical dynamic input–output model taking into consideration the activity of remanufacturing .We try to answer the question, whether the remanufacturing/reuse increases the growth possibility of an economy. We expose a sufficient condition concerning the effectivity of an economy with remanufacturing. By this evaluation we analyse a possible sustainable development of the economy on the basis of the product recovery management of industries.     

Investigation of the Stress State of Pressure Vessels of Inhomogeneous Structure Made of Composite Materials

The problem on the stress-strain state of layered cylindrical shells with bottoms of intricate shape under the action of internal pressure is considered. The elastic system examined is formed by spiral-circular winding. Two variants of the shell bottom structure are investigated. In the first variant, one spiral layer is installed, which leads to great variations in the bottom thickness along the meridian. In the second one, the bottoms are formed according to the zone-winding scheme. The stress state of the shell constructions of the classes considered is determined by solving boundary-value problems for systems of ordinary differential equations. The solution results for cylindrical shells with elliptic bottoms for the two types of winding are given. It is shown that the zone winding leads to smaller deflections and stresses than the conventional ways of reinforcing shell bottoms.     

ML Estimation of the MultivariatetDistribution and the EM Algorithm

Maximum likelihood estimation of the multivariatetdistribution, especially with unknown degrees of freedom, has been an interesting topic in the development of the EM algorithm. After a brief review of the EM algorithm and its application to finding the maximum likelihood estimates of the parameters of thetdistribution, this paper provides new versions of the ECME algorithm for maximum likelihood estimation of the multivariatetdistribution from data with possibly missing values. The results show that the new versions of the ECME algorithm converge faster than the previous procedures. Most important, the idea of this new implementation is quite general and useful for the development of the EM algorithm. Comparisons of different methods based on two datasets are presented.