基于建筑风场数值模拟的网格层块加密法

在建筑风场的数值模拟中,当前普遍采用的离散网格多是计算前一次性布置的固定网格,通常很难适应实际流场变量的变化要求.为提高数值模拟的精度,基于结构化同位网格系统及控制容积离散微分方程的方法,将适应性网格局部加密(AMR)的思想引入到采用压力校正迭代算法的建筑风场模拟中,提出了一种半自适应的层块网格加密方法.该方法可结合误差分析对误差较大的区域网格实行自动判别并实施逐层块状加密.算例分析表明,该方法能在较高的效益下提高数值解的精度.     

自适应网格技术在流体力学界面处理中的应用

自适应网格技术是流体力学数值计算中的重要内容之一．不仅能够处理大变形问题，而且能够大大提高物质界面的计算精度。针对结构网格给出了网格细化结构、误差估计、数据管理、网格产生算法等。通过对强激波双马赫反射及Rayleigh-Taylor界面不稳定性现象进行数值模拟．获得符合精度要求的对激波阵面及物质界面自适应跟踪结果。     

并行有限元网格生成方法及其应用

大型工程数值仿真中,在前处理阶段需要生成千万甚至亿量级的网格,传统的串行网格生成方法由于内存和时间的限制,难以处理如此规模的网格。针对此问题,本文提出了一种大规模网格并行生成方法。首先基于推进波前法对几何模型进行初始体网格划分,接着利用图论理论进行区域分解,并通过表面单元恢复保持其几何精度,然后通过分裂法进行网格的并行生成。将所述方法应用到实际大型工程数值仿真前处理阶段,结果表明所述方法可以获得较好的并行效率,同时所产生的网格质量可以满足后续计算需要。     

非结构混合网格在气动热标模模拟中的应用

非结构混合网格消除了结构网格节点的结构性限制,可以较好地处理边界,同时兼顾了粘性边界层模拟的需求,具有灵活性大、对复杂外形适应能力强和生成耗时短等优点,在飞行器气动特性模拟中得到广泛应用.本文针对非结构混合网格的特点,把前期针对非结构混合网格气动力高精度模拟发展改进的梯度计算方法和Roe格式熵修正方法推广应用到气动热流的数值模拟.以典型钝锥标模外形的高超声速绕流为研究对象,开展了不同网格形式和第一层网格不同间距的影响研究.结果 表明,热流计算时,头部物面网格最好采用四边形或四边形交叉剖分得到的三角形网格,物面法向的网格雷诺数取20左右,为热流计算时非结构混合网格的生成提供了指导,同时验证了计算方法的有效性和可靠性.     

非结构动网格在多介质流体数值模拟中的应用

采用非结构动网格方法对含多介质的流场进行数值模拟.采用改进的弹簧方法来处理由于边界运动而产生的网格变形.采用基于格心的有限体积方法求解守恒型的ALE(Arbitrary Lagrangiall-Eulerian)方程,控制面通量的计算采用HLLC(Hartem,Lax,van Leer,Contact)方法,采用几何构造的方法使空间达到二阶精度,时间离散采用四阶Runge-Kutta方法.物质界面的处理采用虚拟流体方法.本文对含动边界的激波管、水下爆炸等流场进行数值模拟,取得较好的结果,不同时刻界面的位置和整个扩张过程被准确模拟.     

层裂问题的交错网格物质点法研究

交错网格物质点法(Staggered Grid Material Point Method, SGMP)有效地消除了物质点法(Material Point Method, MPM)的跨网格数值误差,是模拟冲击爆炸等涉及大变形与材料破坏问题的一种有效数值分析方法.为了研究层裂问题,论文将质点失效模型引入交错网格物质点法中,对Hopkinson杆碰撞、飞片撞击靶板碰撞和混凝土受爆炸载荷冲击等层裂问题进行了数值模拟研究,结果表明SGMP可以很好地模拟层裂问题.     

基于图像四叉树的改进型比例边界有限元法研究

为提高数值计算的精度, 断裂力学问题的数值模拟需要在裂纹扩展的局部区域采用较密的网格, 而远离裂纹扩展的区域可采用较疏的网格, 且对于裂纹扩展问题的数值模拟, 大多数数值方法又存在局部网格重剖分的问题. 论文提出了一种基于图像四叉树的改进型比例边界有限元法用于模拟裂纹扩展问题, 该方法可根据结构域几何外边界的图像全自动进行四叉树网格剖分, 无需任何人工干预, 网格剖分效率极高, 由于比例边界有限元法本身的优势, 四叉树网格的悬挂节点可以直接地视为新的节点, 无需任何特殊处理. 通过引入虚节点的思想, 将裂纹与四叉树单元边界交叉点作为虚节点, 虚节点的自由度作为附加自由度处理, 并采用水平集函数表征材料内部的裂纹面, 含不连续裂纹面的子域可通过节点水平集函数识别, 使得裂纹扩展时无需进行网格重剖分, 界面的几何特征通过比例边界有限元子域的附加自由度表征. 最后, 通过若干算例验证了该方法的性能, 建议的改进型比例边界有限元法在求解复合型应力强度因子和模拟材料内部裂纹扩展路径时均具有较高的精度.      

无网格MPM法三维前处理系统设计

无网格法在解决爆炸冲击、裂纹扩展、高速碰撞以及大变形等问题时,由于与网格无关,能够避免有限元法因网格畸变计算精度降低或计算终止的矛盾,是目前求解固体力学问题的有效方法之一,也是近年来研究的热点。前处理对于无网格法数值计算是非常重要的环节,基于面向对象程序设计思想,应用C 和OpenGL开发了三维图形交互方式的MPM法前处理系统,提出了应用递归四叉树法、扫描表示法和CSG法划分任意形状非等间距的MPM法质点单元,为建立通用的无网格前处理平台提出了有效的解决方案。     

陡峭山坡风场数值模拟方法研究

地形起伏是造成山地风场复杂多变的主要原因,其风场特性与基于均匀粗糙平面的风场有很大的区别.为准确模拟山坡地形风场,以坡角45°的简化陡峭山坡为研究对象,采用CFD数值模拟方法进行了流场分析,通过与风洞试验和各国规范对比,详细分析了网格分辨率、湍流模型和坡顶局部光滑处理等因素对数值模拟风场精度的影响.结果 表明,采用基本网格尺度布置以及Realizable k-ε湍流模型,在坡顶位置的风压系数值与风洞试验存在一定偏差;在坡顶分离点处采用具有二阶连续性的曲线进行局部光滑,可使得数值模拟所得风速比和风压系数与试验结果更好地吻合,且光滑曲线的过渡段水平距离越短,模拟效果相对越好,坡顶位置的地形加速效应模拟结果与文献试验、中国以及澳大利亚规范具有更好的一致性.     

基于节点分离Lagrange有限元方法的超高速碰撞碎片云数值模拟

提出了一种基于节点分离概念的Lagrange有限元冲击破碎分析的算法,采用节点分离技术实现 了网格断裂,并用畸变侵蚀技术处理严重畸变单元。利用C++编程实现了节点分离计算模型的创建和畸变 侵蚀,结合LS-dyna的重启动分析功能,编程调用LS-dyna求解器以及畸变侵蚀程序实现了破碎分析。分别 用节点分离的Lagrange有限元算法、SPH 无网格算法对超高速碰撞问题进行了数值模拟,并与文献中的实 验结果和Euler有限元数值模拟结果进行了对比。结果表明,节点分离算法具有计算速度快、数值稳定、边界 明确等优点,能准确有效模拟超高速碰撞问题。     

Design for a high energy density Kelvin–Helmholtz experiment

While many high energy density physics (HEDP) Rayleigh–Taylor and Richtmyer–Meshkov instability experiments have been fielded as part of basic HEDP and astrophysics studies, not one HEDP Kelvin–Helmholtz (KH) experiment has been successfully performed. Herein, a design for a novel HEDP X-ray driven KH experiment is presented along with supporting radiation-hydrodynamic simulation and theory.     

Accelerating piston action and plasma heating in high-energy density laser plasma interactions

In the field of high-energy density physics (HEDP), lasers in both the nanosecond and picosecond regimes can drive conditions in the laboratory relevant to a broad range of astrophysical phenomena, including gamma-ray burst afterglows and supernova remnants. In the short-pulse regime, the strong light pressure (>Gbar) associated ultraintense lasers of intensity I > 10¹⁸ W/cm² plays a central role in many HEDP applications. Yet, the behavior of this nonlinear pressure mechanism is not well-understood at late time in the laser–plasma interaction. In this paper, a more realistic treatment of the laser pressure ‘hole boring’ process is developed through analytical modeling and particle-in-cell simulations. A simple Liouville code capturing the phase space evolution of ponderomotively-driven ions is employed to distill effects related to plasma heating and ion bulk acceleration. Taking into account these effects, our results show that the evolution of the laser-target system encompasses ponderomotive expansion, equipartition, and quasi-isothermal expansion epochs. These results have implications for light piston-driven ion acceleration scenarios, and astrophysical applications where the efficiencies of converting incident Poynting flux into bulk plasma flow and plasma heat are key unknown parameters.     

Stress and displacement fields in soft cylindrical multilayers

Multilayered gels play an important role in biomedical engineering as drug delivery vehicles, replacement tissues and bio-mimetic substrates for cell cultures. It has been established that the gel elasticity strongly influences the intended functionalities. In view of this, second-order elastic solutions for the stresses and displacements in cylindrical multilayered hydrogels subjected to various dilatation profiles are developed in this paper. The results emphasize the importance of nonlinearity in gel mechanics, and suggest the possibility of a rational selection of layer elasticities, layer thicknesses and dilatation profiles for improved mechanical responses such as maximum/minimum swelling and multiaxial stress states.     

Fully computerized evaluation of interferograms from fluid flow investigations

An interactive computerized system capable of evaluating interferograms, obtained for example in gas flow studies, has been developed. The program takes full advantage of the facilities offered by the Apple 2 microcomputer and is designed so that the user needs no previous knowledge of computers and/or software. The program is constructed of several main programs invoking other sub-programs; each program can be active at each step of the process and is menu driven. Thus, the user does not have to remember either the questions or their order of appearance.     

Comment on “Studying astrophysical collisionless shocks with counterstreaming plasmas from high power lasers”

It is noted that the evaluation of the mean free path of ions by Park et al. [HEDP 8 (2012), 38–45] is incomplete without taking into account ion–electron collisions for colliding streams of high-speed laser-produced plasmas. As a result, the calculated scale of their interaction, ～few mm, can be explained by Coulomb collisions rather than as claimed in the paper, by the formation of a laboratory collisionless shock that simulates astrophysical shocks mediated by a Weibel-like instability. Advantages of CO₂-laser's irradiation and plane targets for such kind of experiments are discussed.     

User defined finite element for modeling and analysis of active piezoelectric shell structures

This paper presents the user defined nine-node piezoelectric shell element and its implementation within a user element subroutine for modeling and analysis of thin-walled active composite structures. The element has been implemented within commercial finite element software ABAQUS. Application of the element covers modeling of arbitrary thin-walled structures also with complex geometries, whereby the automated mesh generation has been accomplished by developing a Python based interface for meshing procedure. In order to be able to perform the post-processing, a special adaptation of the user element had to be performed for visualization purposes. The implemented element regards the piezoelectric thin layers polarized in the thickness direction and it is based on the \(e_{31}\) piezoelectric effect. It has been also shown that this biquadratic nine-node element based on degenerated shell approach is less prone to locking effects and more suitable for implementation with curved structures. Through several examples the accuracy of the implemented user defined shell element as well as of the Python-based mesh extension has been demonstrated, along with the possibilities for post-processing. Meshing the structures with the nine-node user element is not third party software dependent.     

Fluid flow and heat transfer in microchannels with rectangular cross section

Forced convective heat transfer coefficients and friction factors for flow of water in microchannels with a rectangular cross section were measured. An integrated microsystem consisting of five microchannels on one side and a localized heater and seven polysilicon temperature sensors along the selected channels on the other side was fabricated using a double-polished-prime silicon wafer. For the microchannels tested, the friction factor constant obtained are values between 53.7 and 60.4, which are close to the theoretical value from a correlation for macroscopic dimension, 56.9 for D _h  = 100 μm. The heat transfer coefficients obtained by measuring the wall temperature along the micro channels were linearly dependent on the wall temperature, in turn, the heat transfer mechanism is strongly dependent on the fluid properties such as viscosity. The measured Nusselt number in the laminar flow regime tested could be correlated by which is quite different from the constant value obtained in macrochannels.     

复杂用户网络在移动通信方面的应用研究

在传统的移动蜂窝系统研究中,所有用户都假设具有相同的行为,即他们以相等的概率打出/接听电话,以及他们以相同的移动性在网络内移动. 还有,基础用户网络假设为一个全互联的网络, 这暗示着所有的用户都互相认识,且以同一概率进行相互通讯. 但在现实的环境中,每个用户有各自的熟人(包括家人、亲戚、朋友和同事)通讯录,其通讯录的大小也视乎个别用户而有所不同. 此外,由于用户的移动性取决于各种各样的因素(如工作性质等),不同用户的移动程度也不尽相同. 为了估计移动系统更接近实际的性能,我们使用无尺度网络作为潜在的用户网络模型,即一个用户的熟人数目服从幂律分布,并且研究在3种移动性(即零移动性、相同移动性和幂律移动性)下的移动系统.最后, 通过计算机仿真来估计系统的性能指标(如承载话务量和阻塞率等).结果表明, 用户网络结构以及用户移动性都会对移动系统的性能产生影响.      

Stress–strain and fracture behaviour of 0°/90° and plain weave ceramic matrix composites from tow multi-axial properties

A computationally economic finite-element-based multi-linear elastic orthotropic materials approach has been developed to predict the stress–strain and fracture behaviour of ceramic matrix composites with strain-induced damage. The finite element analysis utilises a solid element to represent a homogenised orthotropic medium of a heterogeneous uni-directional tow. The non-linear multi-axial stress–strain behaviour has been discretised to multi-linear elastic curves, which have been implemented by a user defined subroutine or UMAT in the commercial finite element package, ABAQUS. The model has been used to study the performance of two CMC composites: a SiC (Nicalon) fibre/calcium aluminosilicate (CAS) matrix 0°/90° cross-ply laminate Nicalon/CAS; and, a carbon fibre/carbon matrix–SiC matrix (C/C–SiC) plain weave laminate DLR-XT. The global stress–strain curves with catastrophic fracture behaviour and effects of fibre waviness have been predicted. Comparisons have been made between the predictions and experimental data for both materials. The predicted results when fibre waviness is taken into account compare well with the experimental data.