体育场环状悬挑屋盖脉动风压数值模拟

本文提出了预测悬挑屋盖上脉动风压的实用方法,用数值模拟中稳态计算得到的平均压力、平均速度和湍动能的分布计算大跨屋盖的脉动风压。以昆山体育场为例,通过数值模拟的平均压力系数、湍动能分布及脉动风压与风洞实验结果的比较,验证了用数值模拟方法预测环状悬挑屋盖脉动风压的有效性,并建议了需选用的合适湍流模型和参考点的风速。     

飞机风挡结构抗鸟撞一体化设计技术研究

结合大型非线性有限元分析软件,实践了飞机风挡结构由初步设计、数值仿真以及实验验证的一体化设计过程,建立了适于风挡结构材料的非线性黏弹性本构计算接口程序,对合理选取单元类型、材料模型、边界条件影响以及试验设计等工作细节进行了深入研究,通过实验与仿真计算结果的对比,验证了数值仿真计算模型的精度．为飞机风挡的抗鸟撞结构设计-数值仿真-实验验证一体化技术提供了支持论证．     

混凝土板爆破作用数值分析

混凝土板是爆破拆除经常遇到的爆破对象,是建筑物的基本构件．爆炸应力波峰值随距离衰减是工程爆破参数设计的理论依据之一．通过数值计算,研究爆炸应力波峰值随距离的衰减规律．研究结果表明,与无限混凝土介质相比,爆炸应力波峰值在混凝土板中衰减较快．应力波峰值的衰减指数随装药量增大而减小,衰减指数的变化率也随装药量的增大而减小．     

一种岩石概率型本构及其应用

在岩石材料非均匀特性的基础上,建立了弹塑性概率材料本构．认为材料单元的切线模量与破坏强度都是服从二参数Weibull分布的随机量,通过用户材料子程序将该本构导入到LS—DYNA3D中,推导了程序流程,分析了导入过程中的关键技术．算例验证表明,应用该本构模拟得出的现象与实验观测完全一致,这是确定性本构无法实现的．     

Particle behavior in homogeneous isotropic turbulence

Direct numerical simulations were conducted to investigate the behavior of heavy particles in homogeneous isotropic turbulence. The present study focused on the effect of particle inertia and drift on the autocorrelations of the particle velocity and the fluid seen by particles and the dispersion characteristics of particles. The Lagrangian integral time scale of particles monotonically increased as the magnitude of the particle response time increased, while that of the fluid seen by particles remained relatively constant; it reached a maximum when the particle response time was close to the Kolmolgorov time scale of the flow. Particle dispersion increased as the particle inertia increased for small particles, while for larger particles, it decreased as particle inertia increased; particle eddy diffusion coefficient was maximal, and greater than that of the fluid by about 30%, at the preferential concentration. The concentration field of the particles with τp/τk≈1.0 showed that particles tend to collect in regions of low vorticity (high strain) due to preferential concentration. As the drift velocity of a particle is increased it crosses the paths of fluid elements more rapidly and will tend to lose correlation with its previous velocity faster than a fluid element will. And the correlation of particle velocities along the drift direction is more persistent than that perpendicular to the direction of drift. Simulations also showed that the continuity effect and the crossing-trajectory effect are weakened for particles with infinite inertia.     

A multigrid method with higher‐order discretization schemes

The implementation of the multigrid method into the SIMPLE algorithm presents interesting aspects concerning the mass fluxes conservation on coarser grids, the k–ε turbulence model and the higher‐order discretization schemes. Higher‐order discretization schemes for the convection terms are increasingly used in order to guarantee accuracy in demanding engineering applications. However, when used in single‐grid algorithms, their convergence is considerably slower compared with the first‐order schemes. Unbounded higher‐order schemes offer maximum accuracy, but quite often they do not converge due to their oscillatory behaviour. This paper demonstrates the dual function of the multigrid method: reduction of CPU time and stabilization of the iterating procedure, making it possible to perform computations with the third‐order accurate QUICK scheme in all cases. The method is applied to the calculation of two‐ and three‐dimensional flows with or without turbulence modelling. The results show that the convergence rate of the present algorithm does not deteriorate when QUICK is used and that, if applied on complex engineering cases, large gains in computational time can be achieved. Copyright © 2001 John Wiley & Sons, Ltd.     

非线性热传导方程MLPG/RBF-FD无网格数值模拟

结合无网格局部彼得洛夫-伽辽金(MLPG)方法和径向基函数有限差分(RBF-FD)无网格方法求解非线性热传导问题.MLPG方法属于弱式无网格方法,具有处理边界条件方便的优点,然而因其要做大量的插值、积分运算而计算效率偏低;RBF-FD无网格方法属于强式无网格方法,直接对微分算子进行数值离散,计算效率高,然而其边界条件的...     

Radiative Heat Transfer in a Resistance Heated Floating Zone Furnace: A Numerical Study with FIDAPTM

This paper presents a numerical study of radiative heat transfer in a floating zone (FZ) furnace which was performed by using the commercial finite element program FIDAP^TM. This resistance furnace should provide a temperature higher than the melting temperature of silicon (i.e. T_max ≈ 1500 °C) and a variable temperature gradient at the liquid/solid interface (≥ 25 K/cm). Due to the high working temperatures, heat radiation plays the dominant role for the heat transfer in the furnace. For this reason, the quality of view factors used in the wall‐to‐wall model was carefully inspected with energy‐balance checks. A numerical model with two control parameters is applied to study the influence of material and geometrical parameters on the temperature field. In addition, this model allows us to estimate the internal thermal conditions which were used as thermal boundary conditions for partial 3D simulations. The influences of an optical lens system on the radial symmetry of the temperature field were examined with these partial 3D simulations. Furthermore, we used the inverse modeling method to achieve maximum possible temperature gradients at the liquid/solid interface according to the limitation of maximum available power and the maximum stable height of a melt zone.     

Numerical stability study and error estimation for two implicit schemes in a moving boundary problem

Two algorithms are described [Ferris D. H. (fixed time‐step method) and Gupta and Kumar (variable time‐step method)] that solve a mathematical model for the study of the one‐dimensional moving boundary problem with implicit boundary conditions. Landau's transformation is used, in order to work with a fixed number of nodes at each time‐step. The p.d.e. is discretized using an implicit finite difference scheme. The mathematical model describes the oxygen diffusion in absorbing tissues. An important application is the estimation of time‐variant radiation treatments of cancerous tumors. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 42–61, 2000