A new curvature technique calculation for surface tension contribution in PLIC-VOF method

The volume of fluid (VOF) methods have been used for numerous numerical simulations. Among these techniques used to define the moving interface, the piecewise linear interface reconstruction (PLIC-VOF) is one of the most accurate. A study of the superficial tension impact on two-phase flow with free surface is presented. A new method based on direct staggered grid is developped to include surface tension in PLIC-VOF. The new numerical curvature calculation method doesn't need smoothed colour function and leads to less “spurious current”. This technique is applied to the calculus of surface tension force in the case of the rise of air bubble in viscous liquid and the fall of liquid drop in the same liquid on free surface. Droplets, thin layer and capillarity waves are observed after the free surface rupture for different Bond number. The influence of surface tension calculus is then obvioused and when the drop hit the free surface, wavelets propagate toward the virtual boundaries imposed.     

LOW-CYCLE FATIGUE UNDER NON-PROPORTIONAL LOADING

A series of strain-controlled, low-cycle fatigue experiments have been conducted on 42CrMo steel under various loading paths including circular, square, cruciform, and rectangular paths. Present experiments have shown that there is additional hardening under non-proportional cyclic loading. Non-proportional cyclic additional hardening also results in a shorter life for multiaxial low cycle fatigue. A non-proportionality measure of strain path based on both a physical basis and macromechanical phenomena is proposed. The loading path effect on additional hardening is also described well. Low-cycle fatigue damage accumulation and the evolution process under non-proportional loading is analysed via the Continuum Damage Mechanics Model of Chaboche. A non-proportinality measure is introduced in the damage evolution equation and a modified Coffin-Manson type formula is derived. A novel fatigue life prediction approach based on the critical-plane concept of Brown and Miller is proposed.     

反向凝固母带拉速对结晶器流场影响的数值模拟

在连续统一模型的基础上建立了描述反向凝固结晶器中的流动与传热的数学模型。利用有限差分方法，计算出了反向凝固结晶器中母带厚度的变化规律，随着母带在钢液中停留时间的增加，母带厚度初期迅速增长，中期平衡相持，后期快速回熔；并进一步研究了反向凝固母带拉速对结晶器中流场的影响情况，拉速较小时，紧贴母带的钢液随母带-同运动，对整个流场的影响不大，当拉速较大或拉速变化较大时，部分甚至整个流场都发生了变化。计算结果表明采用数值模拟求解结果与试验数据是吻合的，证明假设合理，模型正确。     

柔性变形机翼后缘拓扑优化设计

为了实现机翼表面的自适应变形和结构轻量化,将柔件机构引入到机翼后缘形状变化结构设计中.应用连续体拓扑优化技术,以实际位移与目标位移之间的偏差为目标函数,材料用量为约束,建立SIMP(solid isotropic material with penalization)密度刚度插值模型.采用Matlab编程对柔性机构进行了优化设计,并对不同参数下的优化结果进行了讨论,最后进行机构的仿真分析.研究结果显示该柔性机构能够实现预期的形状变化,证明了方法的正确性,为柔性机翼设计提供理论基础.     

Structural topology optimization using ant colony optimization algorithm

The ant colony optimization (ACO) algorithm, a relatively recent bio-inspired approach to solve combinatorial optimization problems mimicking the behavior of real ant colonies, is applied to problems of continuum structural topology design. An overview of the ACO algorithm is first described. A discretized topology design representation and the method for mapping ant's trail into this representation are then detailed. Subsequently, a modified ACO algorithm with elitist ants, niche strategy and memory of multiple colonies is illustrated. Several well-studied examples from structural topology optimization problems of minimum weight and minimum compliance are used to demonstrate its efficiency and versatility. The results indicate the effectiveness of the proposed algorithm and its ability to find families of multi-modal optimal design.     

Set-valued Lyapunov functions for difference inclusions

The paper relates set-valued Lyapunov functions to pointwise asymptotic stability in systems described by a difference inclusion. Pointwise asymptotic stability of a set is a property which requires that each point of the set be Lyapunov stable and that every solution to the inclusion, from a neighborhood of the set, be convergent and have the limit in the set. Weak set-valued Lyapunov functions are shown, via an argument resembling an invariance principle, to imply this property. Strict set-valued Lyapunov functions are shown, in the spirit of converse Lyapunov results, to always exist for closed sets that are pointwise asymptotically stable.     

A model of cerebral cortex formation during fetal development using reaction-diffusion-convection equations with Turing space parameters

The cerebral cortex is a gray lamina formed by bodies of neurons covering the cerebral hemispheres, varying in thickness from 1.25 mm in the occipital lobe to 4 mm in the anterior lobe. The brain's surface is about 30 times greater that of the skull because of its many folds; such folds form the gyri, sulci and fissures and mark out areas having specific functions, divided into five lobes. Convolution formation may vary between individuals and is an important feature of brain formation; such patterns can be mathematically represented as Turing patterns. This article describes how a phenomenological model was developed by describing the formation pattern for the gyri occurring in the cerebral cortex by reaction diffusion equations with Turing space parameters. Numerical examples for simplified geometries of a brain were solved to study pattern formation. The finite element method was used for the numerical solution, in conjunction with the Newton–Raphson method. The numerical examples showed that the model can represent cerebral cortex fold formation and reproduce pathologies related to gyri formation, such as polymicrogyria and lissencephaly.     

Discrete choice in commuter space: Small area analysis of diurnal population change in the Tokyo Metropolitan Region

Considerations of metropolitan commuting have often related to urban space as a continuum, and to commuting itself as a process. This approach has led to escalating complexity in the measurement, modeling and analysis of commuting, often with no commensurate results. The present study extends an alternative approach that considers commuter space of a metropolis as a discrete set of small areas. Each small area throughout the metropolis is a nighttime location of some households and, as a possible discrete destination choice of individuals within households, each subarea is also a daytime location of some persons. Each subarea, therefore, may be thought of as attaining two values, corresponding to its daytime (e.g., noontime) and nighttime (e.g., midnight) populations. The diurnal demographic shift of subareas, as a set of binary events, or toggles, between daytime and nighttime populations, is associated with commuter access throughout the metropolis. Diurnal shift offers structural simplicity aimed at estimating small-area daytime populations across a metropolitan region that might suitably precede and complement conventional continuum considerations of commuting. Using data for the Tokyo Metropolitan Region, we identify average household size as the nighttime indicator of small-areas, and daytime population density as their daytime indicator. The potential for the estimation of daytime populations across metropolitan regions where daytime population data are unavailable, is underscored by the clear tendency towards inverse relation of these two indicators of diurnal shift.     

On the probabilistic-deterministic transition involved in a fragmentation process of brittle materials

Dynamic loadings produce high stress waves leading to the fragmentation of brittle materials such as ceramics, concrete, glass and rocks. The main mechanism used to explain the change of the number of fragments with stress rate is a shielding phenomenon. However, under quasi-static loading conditions, a weakest link hypothesis may be applicable. Therefore, depending on the local strain or stress rate, different fragmentation regimes are observed. One regime corresponds to single fragmentation for which a probabilistic approach is needed. Conversely, the multiple fragmentation regime may be described by a deterministic approach. The transition between the two fragmentation regimes is discussed for high performance concrete, glass and SiC ceramics.     

逆向工程中的光滑连接技术

介绍各类曲线段、曲面片光滑连接条件,总结各类曲面的特点和构造方法;给出曲线阶数的求解算法;探讨张量参数β对B样条曲线和曲面的影响,并给出实例和一些结论。