Domain decomposition hybrid method for numerical simulation of bluff body flows:theoretical model and application

The discrete vortex method is not capable of precisely predicting the bluff body flow separation and the fine structure of flow field in the vicinity of the body surface. In order to make a theoretical improvement over the method and to reduce the difficulty in finite-difference solution of N-S equations at high Reynolds number, in the present paper, we suggest a new numerical simulation model and a theoretical method for domain decomposition hybrid combination of finite-difference method and vortex method. Specifically, the full flow. field is decomposed into two domains. In the region of O(R) near the body surface (R is the characteristic dimension of body), we use the finite-difference method to solve the N-S equations and in the exterior domain, we take the Lagrange-Euler vortex method. The connection and coupling conditions for flow in the two domains are established. The specific numerical scheme of this theoretical model is given. As a preliminary application, some numerical simulations for flows at Re=100 and Re-1000 about a circular cylinder are made, and compared with the finite-difference solution of N-S equations for full flow field and experimental results, and the stability of the solution against the change of the interface between the two domains is examined. The results show that the method of the present paper has the advantage of finite-difference solution for N-S equations in precisely predicting the fine structure of flow field, as well as the advantage of vortex method in efficiently computing the global characteristics of the separated flow. It saves computer time and reduces the amount of computation, as compared with pure N-S equation solution. The present method can be used for numerical simulation of bluff body flow at high Reynolds number and would exhibit even greater merit in that case.

Some theoretical aspects of the simplified Navier-Stokes(SNS) equations

This study deals with the formulation, mathematical property and physical meaning of the simplified Navier-Stokes (SNS) equations. The tensorial SNS equations proposed is the simplest in form and is applicable to flow fields with arbitrary body boundaries. The zones of influence and dependence of the SNS equations, which are of primary importance to numerical solutions, are expounded for the first time from the viewpoint of subcharacteristics. Besides, a detailed analysis of the diffusion process in flow fields shows that the diffusion effect has an influence zone globally windward and an upwind propagation greatly depressed by convection. The maximum upwind influential distance of the viscous effect and the relative importance of the viscous effect in the flow direction to that in the direction normal to the flow are represented by the Reynolds number, which illustrates the conversion of the complete Navier-Stokes (NS) equations to the SNS equations for flows with large Reynolds number.

Prediction of mixed-mode cracking direction in random, short-fibre composite-materials

The prediction of cracking direction in composite materials is of significance to the design of composite structures. This paper presents several methods for predicting the cracking direction in the double grooved tension-shear specimen which gives mixed-mode cracking. Five different criteria are used in this analysis: two of them have been used by other investigators and the others are proposed by the present authors. The strain energy density criterion proposed by G.C. Sih is modified to take account of the influence of the anisotropy of the strength on the direction of crack. The two failure criteria of Tsai-Hill and Norris are extended to predict the crack orientation. The stress distributions in the near-notch zone are calculated by using the 8-node quadrilateral isoparametric finite element method. The predictions of all the criteria except one are in good agreement with the experimental measurement. In addition, on the basis of the FEM results, the size of the zone in which the singular term is dominant is estimated.

Application of the strain energy density criterion to ductile fracture

The strain energy density criterion due to Sih is used to predict fracture loads of two thin plates subjected to large elastic-plastic deformation. The prediction is achieved with a finite element analysis which is based on Hill's variational principle for incremental deformations capable of solving gross yielding problems involving arbitrary amounts of deformation. The computed results are in excellent agreement with those obtained in Sih's earlier analysis and with an experimental observation.

EXPERIMENTAL AND THEORETICAL STUDY OF BINARY DROPLET HEAD-ON COLLISIONS IN MEMS

The experimental and theoretical investigations into the head-on collision between a landing droplet with another one resting on the PDMS substrate were addressed in this talk. The colliding process of the two droplets was recorded with highspeed camera. Four different responses after collision were observed in our experiments: complete rebound, coalescence, partial rebound with conglutination, and coalescence accompanied by conglutination. The contact time between the two colliding droplets was found to be in the range of 10-20 milliseconds. For the complete bouncing case, Hertz contact model was applied to estimate the contact time of the binary head-on colliding droplets with both the droplets considered as elastic bodies. The estimated contact time was in good agreement with the experimental result.

Prediction of particle distribution in isotropic turbulence by large-eddy simulation

The recent application of large-eddy simulation (LES) to particle-laden turbulence requires that the LES with a subgrid scale (SGS) model could accurately predict particle distributions. Usually, a SGS particle model is used to recover the small-scale structures of velocity fields. In this study, we propose a rescaling technique to recover the effects of small-scale motions on the preferential concentration of inertial particles. The technique is used to simulate particle distribution in isotropic turbulence by LES and produce consistent results with direct numerical simulation (DNS). Key words: particle distribution, particle-laden turbulence, large-eddy simulation, subgrid scale model.

中国力学学会史 = A history of Chinese society of theoretical and applied mechanics

《中国力学学会史》是《中国学会史丛书》之一。是一部全面系统记述中国力学学会建立与发展历程的专著。 《中国力学学会史》全书30万字，书中不但重点对学会的初创情况、发展过程、组织建设、学术交流、分支机构等进行了专门介绍，还特别收录了记述学会重大活动情况的大事记、名人与学会发展的丰富资料和一些极有史料价值的历史照片，旨在反映学会在不同时期的活动概况及其在中国力学界中发挥的桥梁与纽带作用。 中国力学学会是中国科协的组成部分，也是我国著名的学术团体之一，仅以此书的编著出版，纪念中国科协成立50周年和中国力学学会成立50多周年。本书可供力学界和科技界有关部门及工作者、各学会相关人员、大专院校师生参阅，也可作为组织和开展国内外学术交流研究的参考资料。

Prediction Of Shear-Band Thickness In Metallic Glasses

We derive an explicit expression for predicting the thicknesses of shear bands in metallic glasses. The model demonstrates that the shear-band thickness is mainly dominated by the activation size of the shear transformation zone (STZ) and its activation free volume concentration. The predicted thicknesses agree well with the results of measurements and simulations. The underlying physics is attributed to the local topological instability of the activated STZ. The result is of significance in understanding the origin of inhomogeneous flow in metallic glasses. (C) 2009 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Prediction of Effective Moduli of Carbon Nanotube-Reinforced Composites with Waviness and Debonding

Carbon nanotubes (CNTs) have been regarded as ideal reinforcements of high-performance composites with enormous applications. However, the waviness of the CNTs and the interfacial bonding condition between them and the matrix are two key factors that influence the reinforcing efficiency. In this paper, the effects of the waviness of the CNTs and the interfacial debonding between them and the matrix on the effective moduli of CNT-reinforced composites are studied. A simple analytical model is presented to investigate the influence of the waviness on the effective moduli. Then, two methods are proposed to examine the influence of the debonding. It is shown that both the waviness and debonding can significantly reduce the stiffening effect of the CNTs. The effective moduli are very sensitive to the waviness when the latter is small, and this sensitivity decreases with the increase of the waviness. (C) 2008 Elsevier Ltd. All rights reserved.

Switch Region for Pathogenic Structural Change in Conformational Disease and Its Prediction

Many diseases are believed to be related to abnormal protein folding. In the first step of such pathogenic structural changes, misfolding occurs in regions important for the stability of the native structure. This destabilizes the normal protein conformation, while exposing the previously hidden aggregation-prone regions, leading to subsequent errors in the folding pathway. Sites involved in this first stage can be deemed switch regions of the protein, and can represent perfect binding targets for drugs to block the abnormal folding pathway and prevent pathogenic conformational changes. In this study, a prediction algorithm for the switch regions responsible for the start of pathogenic structural changes is introduced. With an accuracy of 94%, this algorithm can successfully find short segments covering sites significant in triggering conformational diseases (CDs) and is the first that can predict switch regions for various CDs. To illustrate its effectiveness in dealing with urgent public health problems, the reason of the increased pathogenicity of H5N1 influenza virus is analyzed; the mechanisms of the pandemic swine-origin 2009 A(H1N1) influenza virus in overcoming species barriers and in infecting large number of potential patients are also suggested. It is shown that the algorithm is a potential tool useful in the study of the pathology of CDs because: (1) it can identify the origin of pathogenic structural conversion with high sensitivity and specificity, and (2) it provides an ideal target for clinical treatment.

Numerical investigation on a transition prediction model

A new transition prediction model is introduced, which couples the intermittency effect into the turbulence transport equations and takes the characteristics of fluid transition into consideration to mimic the exact process of transition. Test cases include a two-dimensional incompressible plate and a two-dimensional NACA0012 airfoil. Performance of this transition model for incompressible flows is studied, with numerical results consistent to experimental data. The requirement of grid resolution for this transition model is also studied.