Assessment of 2D/3D numerical modeling for deep dynamic stall experiments

The results of computational fluid dynamics (CFD) simulations in two and three spatial dimensions are compared to pressure measurements and particle image velocimetry (PIV) flow surveys to assess the suitability of numerical models for the simulation of deep dynamic stall experiments carried out on a pitching NACA 23012 airfoil. A sinusoidal pitching motion with a 10° amplitude and a reduced frequency of 0.1 is imposed around two different mean angles of attack of 10° and 15°. The comparison of the airloads curves and of the pressure distribution over the airfoil surface shows that a three-dimensional numerical model can better reproduce the flow structures and the airfoil performance for the deep dynamic stall regime. Also, the vortical structures observed by PIV in the flow field are better captured by the three-dimensional model. This feature highlighted the relevance of three-dimensional effects on the flow field in deep dynamic stall.     

Study of linear response identification techniques and reduced‐order model generation for a 2D CFD scheme

Reduced‐Order Models (ROMs) have been the focus of research in various engineering situations, but it is only relatively recently that such techniques have begun to be introduced into the CFD field. The purpose of generating such models is to capture the dominant dynamics of the full set of CFD equations, but at much lower cost. One method that has been successfully implemented in the field of fluid flows is based on the calculation of the linear pulse responses of the CFD scheme coupled with an Eigensystem Realization algorithm (ERA), resulting in a compact aerodynamic model. The key to the models is the identification of the linear responses of the non‐linear CFD code. Two different methods have been developed and reported in literature for linear response identification; the first method linearizes the CFD code and the second method uses Volterra theory and the non‐linear code. As these methods were developed independently they have not previously been brought together and compared. This paper first explains the subtle, but fundamental differences between the two methods. In addition, a series of test cases are shown to demonstrate the performance and drawbacks of the ROMs derived from the different linear responses. The conclusions of this study provide useful guidance for the implementation of either of the two approaches to obtain the linear responses of an existing CFD code. Copyright © 2006 John Wiley & Sons, Ltd.     

Comparison of finite element techniques for 2D and 3D crack analysis under impact loading

For several technical applications the dynamic aspect in fracture mechanics cannot be neglected. When the reliability of components with macroscopic cracks has to be assessed, the consideration of dynamic effects may lead to much higher stress intensity factors than under static conditions. In this paper three different methods to calculate the dynamic stress intensity factor for the mode-I loading of stationary cracks are compared. Based on two- and three-dimensional finite element simulations, the dynamic stress intensity factor is computed with the dynamic J-integral, the modified crack closure integral and the displacement interpretation method. The theoretical fundamentals of all three methods are summarized in the paper and the numerical implementation is explained briefly. Results for different models are shown and compared to findings in the literature.     

基于连续-非连续单元法的三维脆性颗粒冲击破碎特性分析

通过在连续-非连续单元法(CDEM)中引入考虑应变率效应的断裂能本构以及能量统计算法,实现了球体冲击破碎过程中损伤破裂程度及能量演化的定量分析。计算结果表明,冲击破碎过程分为接触蓄能阶段、损伤破碎阶段和碎块飞散阶段。首先,颗粒的部分动能转化为单元弹性变形能,随后这部分变形能和动能迅速转化为摩擦消耗、阻尼消耗及弹簧断裂能,破碎基本完全后碎块继续飞散。不同冲击速度下,颗粒分别出现了反弹、开裂、破碎和粉碎的现象。随冲击速度的增加,D50的变化速率逐渐放缓,破碎块度逐渐趋于稳定;破裂度、损伤度以及平均损伤因子的变化速率先增加后放缓,颗粒破坏以拉伸破坏为主。以上结论可为脆性材料冲击破碎工艺的优化设计提供依据。     

基于细观力学的混凝土数值模型研究

详细分析了界面单元模型、接触面单元模型、CT图像三维重建模型三种细观混凝土数值模型的优缺点,并将数值计算结果与物理试验结果进行了对比.结果表明:混凝土接触面模型克服了界面单元模型中的界面层太厚、计算单元多、计算速度慢的弊端,且解决了原模型裂纹空洞区大的问题;而CT图像三维重建模型与实际物理模型能够一一对应,可以弥补传统随机混凝土骨料模型的不足.不仅可以使混凝土数值模型更接近于实际物理模型,而且可以较好地反映出了混凝土材料内部的力学特性.CT三维重建模型可为最终能够深入研究混凝土材料的细观破损机理开辟新的途径.     

Assessment of different reconstruction techniques for implementing the NVSF schemes on unstructured meshes

Three new far‐upwind reconstruction techniques, New‐Technique 1, 2, and 3, are proposed in this paper, which localize the normalized variable and space formulation (NVSF) schemes and facilitate the implementation of standard bounded high‐resolution differencing schemes on arbitrary unstructured meshes. By theoretical analysis, it is concluded that the three new techniques overcome two inherent drawbacks of the original technique found in the literature. Eleven classic high‐resolution NVSF schemes developed in the past decades are selected to evaluate performances of the three new techniques relative to the original technique. Under the circumstances of arbitrary unstructured meshes, stretched meshes, and uniform triangular meshes, for each NVSF scheme, the accuracies and convergence properties, when implementing the four aforementioned far‐upwind reconstruction techniques respectively, are assessed by the pure convection of several scalar profiles. The numerical results clearly show that New‐Technique‐2 leads to a better performance in terms of overall accuracy and convergence behavior for the 11 NVSF schemes. Copyright © 2013 John Wiley & Sons, Ltd.     

Some recent results on a shell model for 3D turbulence

In this note we report a series of recently published results (by R. Benzi, M.H. Jensen, G. Paladin, G. Parisi, A. Vulpiani and myself), on a Shell Model for the three dimensional cascade of energy in Fully Developed Turbulence. We present numerical and analytical computations of the intermittent corrections to the K41 scaling of the velocity increments in the inertial range.     

基于转换模板的三维实体全六面体网格生成方法

以四面体一六面体基本转换模板为基础,提出一系列具有伸缩性的扩展转换模板．可根据需要选择不同模板及其组合,将四面体分解为不同数量、不同密度过渡形式的六面体单元。这样,初始四面体网格不需要划分得很细,生成的六面体单元数量也可以通过采用不同规格的扩展转换模板而得到控制。提出了基于CAD几何造型的边界节点坐标修正方法．从而使边界网格更好地拟合几何模型边界。     

岩石含水蠕变损伤模型的开发与应用

摘要：为研究岩石含水率对蠕变的影响,将含水开关与蠕变损伤阀值引入蠕变模型,实现了含水劣化与蠕变损伤在本构关系上的耦合,建立了岩石含水蠕变损伤模型。模型与某软岩的三轴压缩蠕变试验结果取得了较好的拟合效果,并由麦夸特法+通用全局优化（LM-UGO）算法反演出了模型参数随岩石含水率的变化规律。在C++语言和FLAC3D中Fish语言的环境下,通过推导模型本构方程的三维差分形式,并利用FLAC3D所提供的接口实现了模型的二次开发。通过对开发的岩石含水蠕变损伤模型进行算例验证和分析,验证了模型的可靠性,也得出了软岩巷道蠕变变形随围岩含水率增大而增加的结论,进一步反映出蠕变模拟计算中考虑岩石含水率的必要性。     

Stabilization of a time integrator for the 3D shallow water equations by smoothing techniques

A smoothing technique is applied to improve the stability of a semi-implicit time integrator for the three-dimensional shallow water equations. In this method the terms involving the vertical direction are treated implicitly. The stability condition on the time step depends only on the horizontal mesh sizes; therefore in the horizontal-direction a smoothing operator is added. Owing to the smoothing, the maximally stable time step increases considerably while the accuracy is hardly affected. Moreover, it turns out that the smoothing operator is efficient on vector and parallel computers.     

基于三维细观模型的全级配混凝土静态力学性能的数值模拟

根据混凝土材料的细观组成和力学特性,研究了骨料几何形状和空间分布规律,建立全级配混凝土三维凸多面体随机细观模型,引入了混凝土细观组份材料的本构模型,分别模拟了单轴、双轴和三轴状态下混凝土的静态力学性能,并建立混凝土梁的三维宏细观分析模型,研究了三点弯曲梁的变形及裂缝扩展情况。结果表明,本文建立的细观力学模型的计算结果与实验数据吻合较好,可以较好地模拟各种复杂应力条件下混凝土的静态力学性能和损伤破坏机理。     

Using digital holographic microscopy for simultaneous measurements of 3D near wall velocity and wall shear stress in a turbulent boundary layer

A digital holographic microscope is used to simultaneously measure the instantaneous 3D flow structure in the inner part of a turbulent boundary layer over a smooth wall, and the spatial distribution of wall shear stresses. The measurements are performed in a fully developed turbulent channel flow within square duct, at a moderately high Reynolds number. The sample volume size is 90 × 145 × 90 wall units, and the spatial resolution of the measurements is 3–8 wall units in streamwise and spanwise directions and one wall unit in the wall-normal direction. The paper describes the data acquisition and analysis procedures, including the particle tracking method and associated method for matching of particle pairs. The uncertainty in velocity is estimated to be better than 1 mm/s, less than 0.05% of the free stream velocity, by comparing the statistics of the normalized velocity divergence to divergence obtained by randomly adding an error of 1 mm/s to the data. Spatial distributions of wall shear stresses are approximated with the least square fit of velocity measurements in the viscous sublayer. Mean flow profiles and statistics of velocity fluctuations agree very well with expectations. Joint probability density distributions of instantaneous spanwise and streamwise wall shear stresses demonstrate the significance of near-wall coherent structures. The near wall 3D flow structures are classified into three groups, the first containing a pair of counter-rotating, quasi streamwise vortices and high streak-like shear stresses; the second group is characterized by multiple streamwise vortices and little variations in wall stress; and the third group has no buffer layer structures.     

基于表面力模型的液滴冲击弹性基底SPH模拟

采用光滑粒子动力学SPH方法建立液滴冲击弹性基底的流固耦合数值模型,给出描述粘性流体和弹性固体运动的SPH离散方程和数值处理格式,引入人工耗散项来抑制标准SPH方法的数值震荡。为模拟液滴的表面张力效应,通过精确检测边界粒子,采用拉格朗日插值方法计算表面法向量和曲率,结合界面理论中的连续表面力CSF方法,建立了适用于自由表面液滴的表面力模型,方形液滴变形的模拟结果与拉普拉斯理论解吻合较好。随后,采用SPH流固耦合模型模拟1.0 mm直径水滴以不同速度(0.2 m/s～3.0 m/s)冲击两种薄板型基底,分析了基底弹性变形对液滴铺展、收缩以及回弹行为的影响。     

Phase transitions induced by extension in a slender SMA cylinder: Analytical solutions for the hysteresis loop based on a quasi-3D continuum model

In this paper, we propose a quasi-3D continuum model to study the rate-independent hysteresis phenomenon in phase transitions of a slender shape memory alloy (SMA) cylinder subject to the uniaxial tension. Based on the three-dimensional field equations and the traction-free boundary conditions, by using a coupled series-asymptotic expansion method, we manage to express the total elastic potential energy of the cylinder in terms of the leading order term of the axial strain. We further consider the rate-independent dissipation effect in a purely one-dimensional setting. The mechanical dissipation functions are also expressed in terms of the axial strain. The equilibrium configuration of the cylinder is then determined by using the principle of maximizing the total energy dissipation. An illustrative example with some special chosen material constants is further considered. Free end boundary conditions are proposed at the two ends of the cylinder. By conducting a phase plane analysis and through some calculations, we obtain the analytical solutions of the equilibrium equation. We find that the engineering stress–strain curves corresponding to the obtained solutions can capture some important features of the experimental results. It appears that the analytical results obtained in this paper reveal the multiple solutions nature of the problem and shed certain light on the instability phenomena during the phase transition process.     

An evaluation of a 3D free‐energy‐based lattice Boltzmann model for multiphase flows with large density ratio

In this paper, the 3D Navier–Stokes (N–S) equation and Cahn–Hilliard (C–H) equations were solved using a free‐energy‐based lattice Boltzmann (LB) model. In this model, a LB equation with a D3Q19 velocity model is used to recover continuity and N–S equations while another LB equation with D3Q7 velocity model for solving C–H equation (Int. J. Numer. Meth. Fluids, 2008; 56 :1653–1671) is applied to solve the 3D C–H equation. To avoid the excessive use of computational resources, a moving reference frame is adopted to allow long‐time simulation of a bubble rising. How to handle the inlet/outlet and moving‐wall boundary conditions are suggested. These boundary conditions are simple and easy for implementation. This model's performance on two‐phase flows was investigated and the mass conservation of this model was evaluated. The model is validated by its application to simulate the 3D air bubble rising in viscous liquid (density ratio is 1000). Good agreement was obtained between the present numerical results and experimental results when Re is small. However, for high‐Re cases, the mass conservation seems not so good as the low‐Re case. Copyright © 2009 John Wiley & Sons, Ltd.     

Multiscale evaluation method of the drag effect on shallow water flow through coastal forests based on 3D numerical simulations

This study presents a method for determining the drag parameter in the 2D shallow water (SW) equation for flows through a coastal forest by conducting a series of 3D numerical simulations (3D NSs). Following the theory of multiscale modeling, an evaluation method procedure is proposed. We first prepare a local test domain that contains a sufficient number of trees to constitute part of a coastal forest. Then, 3D NSs are conducted in this test domain with various inflow conditions. Based on the corresponding results, the momentum losses over the test domain are converted into the drag parameter of the global SW equation. A response surface of the drag parameter is constructed as a function of the flow conditions. The stabilized finite element method is employed for both the local and the global NSs, and the phase-field method is utilized to represent 3D free surfaces. Comparisons between the 2D SW calculation results and the 3D NS results are also performed to verify the validity of the proposed method.     

Experimental and numerical analysis of fretting crack formation based on 3D X-FEM frictional contact fatigue crack model

Nowadays, numerical simulation of 3D fatigue crack growth is easily handled using the eXtended Finite Element Method coupled with level set techniques. The finite element mesh does not need to conform to the crack geometry. Most difficulties associated to complex mesh generation around the crack and the re-meshing steps during the possible propagation are hence avoided. A 3D two-scale frictional contact fatigue crack model developed within the X-FEM framework is presented in this article. It allows the use of a refined discretization of the crack interface independent from the underlying finite element mesh and adapted to the frictional contact crack scale. A stabilized three-field weak formulation is also proposed to avoid possible oscillations in the local solution linked to the LBB condition when tangential slip is occurring. Two basic three-dimensional numerical examples are presented. They aim at illustrating the capacities and the high level of accuracy of the proposed X-FEM model. Stress intensity factors are computed along the crack front. Finally an experimental 3D ball/plate fretting fatigue test with running conditions inducing crack nucleation and propagation is modeled. 3D crack shapes defined from actual experimental ones and fretting loading cycle are considered. This latter numerical simulation demonstrates the model ability to deal with challenging actual complex problems and the possibility to achieve tribological fatigue prediction at a design stage based on the fatigue crack modeling.     

基于三维弹性杆模型的DNA拉伸响应的研究

通过Young-Laplace方程将界面张力引入Kirchhoff方程,并结合Gibbs与Langmuir吸附方程建立了受溶液浓度影响的三维DNA弹性杆模型.基于此模型,引入DNA端部的边界条件,运用打靶法来模拟计算溶液中的DNA链段受溶液界面张力和端部拉力共同作用下的几何构型.进一步分析了在界面能与弹性应变能的耦合作用下,DNA链段平衡构型的变化规律.并在三维DNA弹性杆模型的基础上拟合了DNA链段的FR曲线,其变化趋势与蠕虫链段模型(WLC)的计算结果相一致.     

Development of an efficient numerical model for hail impact simulation based on experimental data obtained from pressure sensitive film

Aircraft are subjects to a number of unpredictable loadings that can seriously affect their performance. In the spirit of ever increasing the safety of passengers, hail impact has been studied. This paper shows the progress that has been made using pressure sensitive film to measure the hail impact event. Moreover, the smooth particle hydrodynamics (SPH) method in LS-DYNA is used to create a numerical model in order to validate the numerical hail model so that it can be used in future advanced simulations of hail impact on components of aircraft. Results show that the SPH method can be effectively used to create a numerical hail model and that pressure sensitive film is a simple and inexpensive tool to capture the experimental data.