水下爆炸载荷作用下加筋板的毁伤模式

为研究加筋板结构在水下爆炸载荷作用下的破坏模式,进行了一系列模型试验。通过对试验结果的分析,将加筋板的毁伤模式总体分为塑性大变形（模式I）、拉伸破坏（模式II）和剪切破坏（模式III）等3大类。又根据载荷强弱和加强筋的相对强度,将每种破坏模式细分为3种子模式。提出了一个量纲一的损伤因子,作为预报毁伤模式的判据。     

杆系结构系统失效模式识别的几何法及其应用

在结构可靠性分析中,由于外载、材料性能参数等的随机性,在所有可能的失效模式中,必有少数几个失效模式其出现概率比其他失效模式要大得多,因此求取这为数不多的临界失效模式就是结构系统可靠性分析中一个主要课题.本文提出了求取杆系结构系统失效模式(即机构)的几何法.它是从几何方程出发,求得一机构,然后在此机构的基础上经过简单的矩阵相加.相乘就可得到其他机构.同时采用比例加载乘子λ这一物理意义明潦的参数来作为判别是否为临界机构的准则,效果很好.算例表明了此法的可行性,有效性.     

电场作用下流动聚焦的实验研究

通过在流动聚焦的同轴液-气射流区域施加电场, 开展了电场力和气动力共同作用下锥形以及带电射流的不稳定性特性实验研究. 实验在精密设计的流动聚焦装置上完成, 分析了外部电压、气体压力降和液体流量等主要控制参数对流动聚焦过程的影响, 获得了锥形的振动模式和稳定模式及其之间的转换, 得到了射流的滴模式、轴对称模式、共存模式和非轴对称模式及其转换并定量分析了电场对射流尺寸参数的影响. 结果表明, 相比于单一的流动聚焦, 该方法能够增强锥形的稳定性, 促进液体射流雾化, 减小颗粒的直径, 因此在科技领域和工程实际中具有重要的应用价值.      

发育生物学中模式形成的力学模型

对所有多细胞生物体而言,其个体发育均涉及由初始的单个细胞开始,经过在时间和空间上细胞有序地增殖、凋亡、分化、迁移等诸多过程,并最终生成物种特异的生物模式.对模式形成的机制的探讨一直是发育生物学的中心课题.目前已就这个问题积累了大量的分子生物学、生物化学、数学、力学等多学科的研究数据并提出了一些模式形成的理论,然而,迄今为止,模式形成的真正机制仍然很不清楚而需进一步的深入探索.本文简要介绍了生物发育过程中模式形成的一些典型的力学模型,重点在于力学模型的建立及其模型本身的介绍,而对其数值模拟和具体应用很少涉及.      

湍流模化的现状及发展趋势

本文讨论了湍流模化研究的发展现状,提出了今后湍流模式的可能改进途径。根据一组湍流封闭假定得到了各种2阶湍流模式,如Reynolds应力模式(RSM),代数应力模式(k-ε-A)及涡粘性模式(k-ε-E)。以自由剪切流、空腔流及通过偏心槽的流动为例作出了预报。虽然至今还不存在一个完备湍流模式,但2阶湍流模式已经具备了某些预报能力。湍流模化的不完备性可能要归咎于各向同性耗散和单湍流尺度等假定的不适当。利用多重湍流尺度概念,包括利用湍流涡的分维,可以改善湍流的预报。     

用自动矩阵方法枚举结构的主要失效模式

本文发展了一个考虑元件强度的自动矩阵方法，并基于这一方法和实际工程结构的特点，提出了一个枚举结构体系主要失效模式的方法，这一方法简单且能有效地枚举出结构体的主要失效模式。     

关于非传递零能模式的讨论

通常认为缩减积分引起的非传递零能模式在两个以上单元网格中是无害的.本文对此举出了一个反例.基于一约束网格的特征值和特征向量试验,提出了非传递零能模式在单元网格中仍可能以低能模式的形式出现,使得到的解不合理.     

激励小尺度模式在湍流圆管射流中的应用

采用非涡黏性的激励小尺度（Stimulated Small Scale)模式对空间发展的轴对称湍流圆管射流进行了大涡模拟。以雷诺数为10000的流动为例,考证了激励小尺度模式在自由剪切流模拟中的可行性,描述了湍流强度、雷诺应力和湍流耗散量的变化,同时与标准的Smagorinsky涡黏性模式的计算结果进行了比较。数值结果显示,激励小尺度模式能够更为合理地描述湍流的耗散特性和能量传输特性,从而较为准确地展示出空间发展射流中由于流动不稳定而出现的旋涡产生、发展、破碎及合并等过程。     

黄土地层地表径流下潜模式与地质灾害

指出在黄土滑坡、地裂缝、黄土洞穴等地质灾害发生发展过程中下渗重力水所起的作用最为关键,对下渗重力水的基本特点进行了全面归纳总结。以黄土洞穴灾害为例,将黄土暗穴形成过程中黄土地层地表径流的下潜模式概化为渗透重力水漫渗型下潜模式、地表径流集中灌入下潜型模式、渗透重力水沿优势渗流通道集中下潜型模式以及混合灌渗集中下潜型模式等4种基本模式。     

用自动矩阵力法枚举结构的主要失效模式

本文发展了一个考虑元件强度的自动矩阵力法，并基于这一方法和实际工程结构的特点，提出了一个枚举结构体系主要失效模式的方法，这一方法简单且能有效地枚举出结构体系的主要失效模式．     

Turbulent particle dispersion in arbitrary wall-bounded geometries: A coupled CFD-Langevin-equation based approach

A Lagrangian continuous random walk (CRW) model is developed to predict turbulent particle dispersion in arbitrary wall-bounded flows with prevailing anisotropic, inhomogeneous turbulence. The particle tracking model uses 3D mean flow data obtained from the Fluent CFD code, as well as Eulerian statistics of instantaneous quantities computed from DNS databases. The turbulent fluid velocities at the current time step are related to those of the previous time step through a Markov chain based on the normalized Langevin equation which takes into account turbulence inhomogeneities. The model includes a drift velocity correction that considerably reduces unphysical features common in random walk models. It is shown that the model satisfies the well-mixed criterion such that tracer particles retain approximately uniform concentrations when introduced uniformly in the domain, while their deposition velocity is vanishingly small, as it should be. To handle arbitrary geometries, it is assumed that the velocity rms values in the boundary layer can locally be approximated by the DNS data of fully developed channel flows. Benchmarks of the model are performed against particle deposition data in turbulent pipe flows, 90° bends, as well as more complex 3D flows inside a mouth-throat geometry. Good agreement with the data is obtained across the range of particle inertia.     

A general one-equation turbulence model for free shear and wall-bounded flows

The purpose of thiswork is to introduce a complete and general one-equation model capable of correctly predicting a wide class of fundamental turbulent flows like boundary layer, wake, jet, and vortical flows. The starting point is the mature and validated two-equation k−ω turbulence model of Wilcox. The newly derived one-equation model has several advantages and yields better predictions than the Spalart-Allmaras model for jet and vortical flows while retaining the same efficiency and quality of the results for near-wall turbulent flows without using a wall distance. The derivation and validation of the new model using findings computed by the Spalart-Allmaras and the k−ω models are presented and discussed for several free shear and wall-bounded flows.     

A Dissipation Rate Equation for Low-Reynolds-Number and Near-Wall Turbulence

Two-equation turbulence models are usually formulated for specific flow types and are seldom validated against a variety of flows to account for near-wall and low-Reynolds-number effects simultaneously. In addition to low-Reynolds-number effects, near-wall flows also experience wall blocking, which is absent in free flows. Consequently, near-wall modifications to two-equation models could be quite different from low-Reynolds-number corrections. Besides, it is known that existing two-equation models perform poorly when used to calculate plane wall jets and two-dimensional backstep flows. These problems could be traced to the modeling of the dissipation rate equation. In this paper an attempt is made to improve the modeling of the dissipation rate equation so that it could successfully predict both free and wall-bounded shear flows including plane wall jets and backstep flows. The predictions are compared with experimental and direct numerical simulation data whenever available. Most of the data used are obtained at low Reynolds numbers. Good correlation with data is obtained. Therefore, for the first time, a model capable of correctly predicting free and wall-bounded shear flows, backstep flows, and plane wall jets is available. Received: 12 December 1995 and accepted 12 November 1996     

Direct numerical simulation of a freely decaying turbulent interfacial flow

Whereas Large Eddy Simulation (LES) of single-phase flows is already widely used in the CFD world, even for industrial applications, LES of two-phase interfacial flows, i.e. two-phase flows where an interface separates liquid and gas phases, still remains a challenging task. The main issue is the development of subgrid scale models well suited for two-phase interfacial flows. The aim of this work is to generate a detailed data base from direct numerical simulation (DNS) of two-phase interfacial flows in order to clearly understand interactions between small turbulent scales and the interface separating the two phases. This work is a first contribution in the study of the interface/turbulence interaction in the configuration where the interface is widely deformed and where both phases are resolved by DNS. To do this, the interaction between an initially plane interface and a freely decaying homogeneous isotropic turbulence (HIT) is studied. The densities and viscosities are the same for both phases in order to focus on the effect of the surface tension coefficient. Comparisons with existing theories built on wall-bounded or free-surface turbulence are carried out. To understand energy transfers between the interfacial energy and the turbulent one, PDFs of the droplet sizes distribution are calculated. An energy budget is carried out and turbulent statistics are performed including the distance to the interface as a parameter. A spectral analysis is achieved to highlight the energy transfer between turbulent scales of different sizes. The originality of this work is the study of the interface/turbulence interactions in the case of a widely deformed interface evolving in a turbulent flow.     

Compressibility effects due to turbulent fluctuations

This paper deals with intrinsic effects of compressibility, i.e. with dilatation fluctuations in response to pressure fluctuations. Three different types of turbulent flows are considered in more detail: homogeneous turbulent shear flow, wall-bounded turbulent shear flow and shock/turbulence interaction. A survey of the present knowledge in this field, mainly based on DNS data, is given. Using the linear inviscid perturbation equations a direct link between fluctuations of dilatation and of velocity in the direction of mean shear is presented for homogeneous shear flow. This relation might form the basis for a more universal pressure-dilatation model. It is conjectured that the insignificance of intrinsic compressibility effects in wall-bounded supersonic shear flow is mainly due to the impermeability constraint of the wall. To this end, a linear stability analysis of supersonic channel flow along cooled, but permeable walls has been performed based on Coleman et al.'s [5] mean flow data. It shows an increase in the moduli of eigenfunctions related to compressibility, like pressure, and in moduli of quantities derived from eigenfunctions such as ‘pressure dilatation’ and squared dilatation. Although these results do not prove our hypothesis they provide hints in this direction. Shock/turbulence interaction is viewed as a source of compressibility. Former DNS data of Hannappel and Friedrich [10] for shock/isotropic turbulence interaction showing the effect of compressibility on the amplification of fluctuations are interpreted based on linear perturbation equations.     

Simulation and Modelling of Turbulent Trailing-Edge Flow

Computations of turbulent trailing-edge flow have been carried out at a Reynolds number of 1000 (based on the free-stream quantities and the trailing-edge thickness) using an unsteady 3D Reynolds-Averaged Navier–Stokes (URANS) code, in which two-equation (k–ε) turbulence models with various low-Re near wall treatments were implemented. Results from a direct numerical simulation (DNS) of the same flow are available for comparison and assessment of the turbulence models used in the URANS code. Two-dimensional URANS calculations are carried out with turbulence mean properties from the DNS used at the inlet; the inflow boundary-layer thickness is 6.42 times the trailing-edge thickness, close to typical turbine blade flow applications. Many of the key flow features observed in DNS are also predicted by the modelling; the flow oscillates in a similar way to that found in bluff-body flow with a von Kármán vortex street produced downstream. The recirculation bubble predicted by unsteady RANS has a similar shape to DNS, but with a length only half that of the DNS. It is found that the unsteadiness plays an important role in the near wake, comparable to the modelled turbulence, but that far downstream the modelled turbulence dominates. A spectral analysis applied to the force coefficient in the wall normal direction shows that a Strouhal number based on the trailing-edge thickness is 0.23, approximately twice that observed in DNS. To assess the modelling approximations, an a priori analysis has been applied using DNS data for the key individual terms in the turbulence model equations. A possible refinement to account for pressure transport is discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Turbulent boundary-layer analysis using finite elements

The use of finite element methods for turbulent boundary-layer flow is relatively recent and of limited extent.¹ In the present study, we extend the group variable approach of Fletcher and Fleer^2,3 to treat turbulent boundary layer flows with heat transfer using a two-equation turbulence model. The main concepts in the formulations include a Dorodnitsyn-type transformation which uses a velocity component as the transverse variable, a ‘variational’ formulation for the transformed equations using special test functions and development of a two-equation turbulence model in terms of the turbulent kinetic energy and turbulence dissipation rate as additional field variables. Several numerical test cases have been examined comparing the results with finite difference calculations and comparing the two-equation turbulence model with an algebraic turbulence model.     

Nonlinear turbulence models for predicting strong curvature effects

Prediction of the characteristics of turbulent flows with strong streamline curvature, such as flows in turbomachines, curved channel flows, flows around airfoils and buildings, is of great importance in engineering applications and poses a very practical challenge for turbulence modeling. In this paper, we analyze qualitatively the curvature effects on the structure of turbulence and conduct numerical simulations of a turbulent Uduct flow with a number of turbulence models in order to assess their overall performance. The models evaluated in this work are some typical linear eddy viscosity turbulence models, nonlinear eddy viscosity turbulence models （NLEVM） （quadratic and cubic）, a quadratic explicit algebraic stress model （EASM） and a Reynolds stress model （RSM） developed based on the second-moment closure. Our numerical results show that a cubic NLEVM that performs considerably well in other benchmark turbulent flows, such as the Craft, Launder and Suga model and the Huang and Ma model, is able to capture the major features of the highly curved turbulent U-duct flow, including the damping of turbulence near the convex wall, the enhancement of turbulence near the concave wall, and the subsequent turbulent flow separation. The predictions of the cubic models are quite close to that of the RSM, in relatively good agreement with the experimental data, which suggests that these models may be employed to simulate the turbulent curved flows in engineering applications.     

LES-based filter-matrix lattice Boltzmann model for simulating fully developed turbulent channel flow

In this paper, a three-dimensional filter-matrix lattice Boltzmann (FMLB) model based on large eddy simulation (LES) was verified for simulating wall-bounded turbulent flows. The Vreman subgrid-scale model was employed in the present FMLB–LES framework, which had been proved to be capable of predicting turbulent near-wall region accurately. The fully developed turbulent channel flows were performed at a friction Reynolds number Re_τ of 180. The turbulence statistics computed from the present FMLB–LES simulations, including mean stream velocity profile, Reynolds stress profile and root-mean-square velocity fluctuations greed well with the LES results of multiple-relaxation-time (MRT) LB model, and some discrepancies in comparison with those direct numerical simulation (DNS) data of Kim et al. was also observed due to the relatively low grid resolution. Moreover, to investigate the influence of grid resolution on the present LES simulation, a DNS simulation on a finer gird was also implemented by present FMLB–D3Q19 model. Comparisons of detailed computed various turbulence statistics with available benchmark data of DNS showed quite well agreement.     

Direct numerical simulation of fluid flow in a 5x5 square rod bundle

Characterization of parallel flow through rod bundles is of key importance in assessing the performance and safety of several engineering systems, including a majority of nuclear reactor concepts. Inhomogeneities in the bundle cross-section can present complex flow phenomena, including varying local conditions of turbulence. With the ever-increasing capabilities of high-performance computing, Direct Numerical Simulation (DNS) of turbulent flows is becoming more feasible. Through resolving all scales of turbulence, DNS can serve as a “numerical experiment,” and can provide substantial insight into flow physics, but at considerable computational cost. Thus to date, the DNS in open literature for rod bundle flows is relatively scarce, and largely limited to unit-cell domains. Since wall effects are important in rod bundle flows, a multiple-pin DNS study can expand understanding of rod bundle flows while providing valuable reference data for evaluating reduced-resolution techniques. In this work, DNS of a 5x5 square bare rod bundle representative of typical light water reactor fuel dimensions was performed using the spectral element code Nek5000. Turbulent microscales based on an advanced Reynolds-Averaged Navier–Stokes model were used to establish the required DNS resolution. Velocity and Reynolds stress fields are analyzed in detail, and invariant analysis is used for further investigation into flow physics. The results show stark changes in the structure of turbulence in the edge gaps, suggesting the presence of gap vortices in these regions. In addition, turbulent kinetic energy budgets are presented to more fully illustrate the various turbulent processes. These data can prove useful for rigorous evaluation of lower-fidelity turbulence modeling approaches.