Characteristic Dimensionless Numbers in Multi-Scale and Rate-Dependent Processes

For this sake, the macroscopic equations of mechanics and the kinetic equations of the microstructural transformations should form a unified set that be solved simultaneously. As a case study of coupling length and time scales, the trans-scale formulation

求解Euler方程的区域分解方法与并行算法

将复杂形状区域划分成多块子区域，研究发展了一种多块区域之间迎风守恒型的内边界耦合方法，实现相邻子区域解的光滑过渡，使多区耦合得到总体流场的数值解。对二维翼型跨音速流动和圆板形隆起物超音速流动等进行了分区数值计算，并将计算结果与单区计算结果和实验结果作了比较。并行分区计算引入“先进先出”的同步控制等待机制，实现了高效率并行计算，还分析了影响并行效率的主要因素。

A Lagrangian lattice Boltzmann method for Euler equations

A Lagrangian lattice Boltzmann method for solving Euler equations is proposed. The key step in formulating this method is the introduction of the displacement distribution function. The equilibrium distribution function consists of macroscopic Lagrangian variables at time steps n and n + 1. It is different from the standard lattice Boltzmann method. In this method the element, instead of each particle, is required to satisfy the basic law. The element is considered as one large particle, which results in simpler version than the corresponding Eulerian one, because the advection term disappears here. Our numerical examples successfully reproduce the classical results.

Effects of the Spanwise Characteristic Length on the Transition Feature in the Wake of a Cylinder

The transition features of the wake behind a uniform circular cylinder at Re = 200, which is just beyond the critical Reynolds number of 3-D transition, are investigated in detail by direct numerical simulations of 3-D incompressible Navier-Stokes equations. The spanwise characteris-tic length determines the transition features and global properties of the wake.

A Study on the Microstructure Characteristic of the Cold-Rolled Deformed Nanocrystalline Nickel

In this paper the microstructure characteristic of the cold-rolled deformed nanocrystalline Nickel metal has been studied by transmission electron microscopy (TEM). The results show that there were step structures near by grain boundary (GB), and the contrast of stress field in front of the step corresponds to the step in the shape. It indicates that the interaction between twins and dislocations is not a necessary condition to realizing the deformation. In the later stage of the deformation when the grain size became about 100 nm, the deformation occurs only depend upon the moving of the boundary of the stack faults (SFs) which result from the imperfection dislocations emitted from GBs. In the other word, the movement of the boundary dislocations of SFs results to growing-up of the size of the SFs, therefore realizes deformation. However, when the size of stack faults grows up, the local internal stress which is in front of the step gradually becomes higher. When this stress reach a critical value stopping the gliding of the partial dislocations, the SFs will stop growing up and leave a step structure behind.

一种易于并行求解Euler方程组的分区技术

在具有复杂边界的计算区域内,求解偏微分方程组时,经常需要分区和并行计算,分区方法直接关系到数值计算的并行化程度,本文在应用时间算子分裂方法求解Euler方程组的过程中,提出了一种非常容易实现并行化计算的分区技术.

二维Euler方程的非结构网格再生成自适应方法

给出了提高二维Euler方程定常解质量的非结构网格自适应方法和模拟结果。计算了无粘激波在固壁上的反射、NAC0012翼型跨声速绕流和马赫数为3的前台阶绕流，自适应效果较好。

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.

Euler方程的双分布函数格子Boltzmann Godunov方法

应用双分布函数系统，通过Ｇｏｄｕｎｏｖ分解，构造了一维Ｅｕｌｅｒ方程的格子Ｂｏｌｔｚｍａｎｎ算法。解决了传统格子气固有的ＧＣ问题与能量方程之间的矛盾，实现了分布函数与宏观物理量之间的一一对应。

基于CE/SE方法的二维Euler型多物质流体弹塑性问题计算

将CE/SE方法推广到二维固体流体弹塑性问题的数值计算,同时结合杂交粒子水平集方法追踪物质界面和合适的边界条件,提出一套完整的二维Euler型流体弹塑性计算方案.通过长钨杆侵彻装甲钢实验的数值模拟,对方法的精度和有效性进行验证.

全机绕流Euler方程多重网格分区计算方法

全机三维复杂形状绕流数值求解只能采用分区求解的方法，本文采用可压缩Euler方程有限体积方法以及多重网格分区方法对流场进行分区计算。数值方法采用改进的van Leer迎风型矢通量分裂格式和MUSCL方法，基于有限体积方法和迎风型矢通量分裂方法，建立一套处理子区域内分界面的耦合条件。各个子区域之间采用显式耦合条件，区域内部采用隐式格式和局部时间步长等，以加快收敛速度。计算结果飞机表面压力分布等气动力特性与实验值进行了比较，二者基本吻合。计算结果表明采用分析“V”型多重网格方法，能提高计算效率，加快收敛速度达到接近一个量级。根据全机数值计算结果和可视化结果讨论了流场背风区域旋涡的形成过程。

A Multiscale-Domain Decomposition Method for High Reynolds Number Flows

The high Reynolds number flow contains a wide range of length and time scales, and the flow domain can be divided into several sub-domains with different characteristic scales. In some sub-domains, the viscosity dissipation scale can only be considered in a certain direction; in some sub-domains, the viscosity dissipation scales need to be considered in all directions; in some sub-domains, the viscosity dissipation scales are unnecessary to be considered at all. For laminar boundary layer region, the characteristic length scales in the streamwise and normal directions are L and L Re-1/ 2 , respectively. The characteristic length scale and the velocity scale in the outer region of the boundary layer are L and U, respectively. In the neighborhood region of the separated point, the length scale l<Euler equations, NS- and diffusion parabolized (DP) NS equations. In this paper, a new multiscale-domain decomposition method is developed for the high Reynolds number flow. First, the whole domain is decomposed to different sub-domains with the different characteristic scales. Then the different dominant equation of all sub-domains is defined according to the diffusion parabolized (DP) theory of viscous flow. Finally these different equations are solved simultaneously in whole computational region. For numerical tests of high Reynolds numerical flows, two-dimensional supersonic flows over rearward and frontward steps as well as an interaction flow between shock wave and boundary layer were solved numerically. The pressure distributions and local coefficients of skin friction on the wall are given. The numerical results obtained by the multiscale-domain decomposition algorithm are well agreement with those by NS equations. Comparing with the usual method of solving the Navier-Stokes equations in the whole flow, under the same numerical accuracy, the present multiscale domain decomposition method decreases CPU consuming about 20% and reflects the physical mechanism of practical flow more accurately.