A mixed-timescale SGS model for thermal field at various Prandtl numbers

A new subgrid-scale (SGS) model for the thermal field is proposed. The model is an extended version of the mixed-timescale (MTS) SGS model for velocity field by Inagaki et al. (2005), which has been confirmed to be a refined SGS model for velocity field suited to engineering-relevant practical large eddy simulation (LES). In the proposed model for the thermal field, a hybrid timescale between the timescales of the velocity and thermal fields is introduced in a manner similar to velocity-field modeling. Thus, the present model dispenses with an ambiguous SGS turbulent Prandtl number, like the dynamic SGS model. In addition, the wall-limiting behavior of turbulence is satisfied, which is not in the original MTS model, by incorporating the wall-damping function for LES based on the Kolmogorov velocity scale proposed by Inagaki et al. (2010). The model performance is tested in plane channel flows at various Prandtl numbers, and the results show that this model gives the ratio of the timescales between the velocity and thermal fields similar to that obtained using the dynamic Smagorinsky model with locally calculated model parameters. It is also shown that the proposed model predicts better mean and fluctuating temperature profiles in cooperation with the revised MTS model for the velocity field, than the Smagorinsky model and the dynamic Smagorinsky model. The present model is constructed with fixed model parameters, so that it does not suffer from computational instability with the dynamic model. Thus, it is expected to be a refined and versatile SGS model suited for practical LES of the thermal field.     

CFD ANALYSIS OF TURBULENT FLOW PAST SQUARE CYLINDER USING DYNAMIC LES

Computational fluid dynamics (CFD) analysis of turbulent flow past a square cylinder is conducted using large eddy simulation (LES). In the first part of the paper, the implementation of the three-dimensional (3-D) LES computation based on the conventional standard Smagorinsky model (S model) shows a good prediction capability compared with that of turbulence models based on Reynolds-averaged Navier-Stokes equations (RANS models). Next, the computation using the dynamic Smagorinsky model (DS model) is carried out, and the improvement by the use of DS model is described. Lastly, the Lagrangian dynamic Smagorinsky model (LDS model) is introduced to overcome the disadvantages of DS model, e.g., calculation instability. Results obtained from the various SGS models are compared with those from experiments.     

穿甲燃烧弹侵彻陶瓷复合装甲和玻璃复合装甲的FEM-SPH耦合计算模型

为了提高小口径穿甲燃烧弹侵彻陶瓷复合装甲和玻璃复合装甲（透明装甲）的仿真分析精度,本文将传统的FEM（finite element method）-SPH（smooth particle hydrodynamics）耦合计算模型中穿甲燃烧弹弹芯的有限元模型和JC（Johnson-Cook）材料模型分别替换为SPH模型和JH2（Johnson-Holmquist-ceramics）材料模型,提出了新型FEM-SPH耦合计算模型。研究表明,新型FEM-SPH耦合计算模型可以有效模拟弹芯碎裂现象,减少SPH粒子和有限元耦合计算量,进而显著提高仿真模型的计算精度和计算效率,并给出了新型FEM-SPH耦合计算模型的有限元/粒子尺度和建模尺寸的优选结果。     

中心逼近式灰色GM(1,1)模型在滑坡变形预测中的应用

黄龙西村滑坡位于甘肃天水,属黄土高势能滑坡,滑体体积3.9×10⁵m³,基底为花岗闪长岩。为了提高滑坡灰色GM(1,1)模型的预测精度,采用一种改变背景值的方法--中心逼近式灰色GM(1,1)模型。通过黄龙西村滑坡实例验证分析,结果表明中心逼近式灰色GM(1,1)模型的预测值与该滑坡实际监测值十分接近,且其残差平方和及平均误差百分比明显比传统灰色GM(1,1)模型的残差平方和及平均误差百分比小,具有较高的预测精度。同时,可通过调整模型中参数m的取值,使中心逼近式灰色GM(1,1)模型具有更高的预测精度。经计算,当m=6时,中心逼近式灰色GM(1,1)模型的预测精度比传统灰色GM(1,1)模型提高了5.34%。     

滑坡预报模型的质量检验研究

 
滑坡预报模型主要用来预报滑坡发生的时间,但是如何检验效果是滑坡预报者和决策者面临的一个重要课题。本文通过提出滑坡拟合效果指标（包括后验差指标、模型拟合效率指数和均方根误差）和试预报效果指标（包括试报效果指标和相关系数指标）,初步建立了滑坡的预报质量检验模型。依据建立的检验模型,针对链子崖危岩体监测资料分别采用灰色GM（1, 1）数学模型、三次指数平滑模型和时间序列模型进行预报。预报结果表明：三种模型中,综合指标GM（1, 1）数值最高,建议该滑坡采用GM（1, 1）进行预报。通过模型的综合分析评判,不仅分析了模拟效果,而且直接比较了模型间的拟合效果,并为建立模拟效果好的模型提供了充分的依据,从而表明滑坡预报质量检验模型是一种有效、实用的方法。


     

抗差自适应模型预测滤波及其在组合导航中的应用

为了提高捷联惯导(SINS)/天文导航(CNS)/合成孔径雷达(SAR)组合导航系统的定位精度,在吸收模型预测滤波和抗差自适应滤波算法优点的基础上,提出了一种新的抗差自适应模型预测滤波算法。该算法首先利用模型预测滤波估计出系统模型误差,并对其进行实时修正,以抑制系统模型误差对导航解算精度的影响;然后利用抗差自适应因子控制观测异常,抑制观测噪声对导航解算精度的影响。将提出的算法应用于SINS/CNS/SAR组合导航系统进行仿真验证,并与抗差自适应滤波进行比较,结果表明,提出的算法得到的姿态误差、速度误差和位置误差分别在[0.2,0.2]、[0.3m/s,0.3m/s]和[6 m,6 m]以内,滤波性能明显优于抗差自适应滤波算法,说明该算法能有效抑制系统模型误差及观测异常对导航解的影响,提高组合导航的解算精度。     

Fourier neural operator approach to large eddy simulation of three-dimensional turbulence

Fourier neural operator (FNO) model is developed for large eddy simulation (LES) of three-dimensional (3D) turbulence. Velocity fields of isotropic turbulence generated by direct numerical simulation (DNS) are used for training the FNO model to predict the filtered velocity field at a given time. The input of the FNO model is the filtered velocity fields at the previous several time-nodes with large time lag. In the a posteriori study of LES, the FNO model performs better than the dynamic Smagorinsky model (DSM) and the dynamic mixed model (DMM) in the prediction of the velocity spectrum, probability density functions (PDFs) of vorticity and velocity increments, and the instantaneous flow structures. Moreover, the proposed model can significantly reduce the computational cost, and can be well generalized to LES of turbulence at higher Taylor-Reynolds numbers.     

Integrated numerical model for vegetated surface and saturated subsurface flow interaction

The construction of an integrated numerical model is presented in this paper to deal with the interactions between vegetated surface and saturated subsurface flows. A numerical model is built by integrating the previously developed quasi-three-dimensional (Q3D) vegetated surface flow model with a two-dimensional (2D) saturated groundwater flow model. The vegetated surface flow model is constructed by coupling the explicit finite volume solution of 2D shallow water equations (SWEs) with the implicit finite difference solution of Navier-Stokes equations (NSEs) for vertical velocity distribution. The subsurface model is based on the explicit finite volume solution of 2D saturated groundwater flow equations (SGFEs). The ground and vegetated surface water interaction is achieved by introducing source-sink terms into the continuity equations. Two solutions are tightly coupled in a single code. The integrated model is applied to four test cases, and the results are satisfactory.     

CFD modeling using heterogeneous reaction kinetics for catalytic dehydrogenation syngas reactions in a fixed-bed reactor

A comprehensive 2D computational fluid dynamics （CFD） model was developed to simulate the flow behavior and catalytic dehydrogenation reaction of syngas in a heterogenous fixed-bed reactor （FBR）. The model combined the porous medium CFD model with a reaction kinetics model. To acquire an accu- rate reaction kinetics model, a comprehensive reaction mechanism was studied for the heterogeneous catalytic dehydrogenation reaction ofsyngas over a supported metal catalyst. Based on the reaction mech- anism and a statistical test, a reliable kinetics model was proposed. The CFD model combined with the above kinetics model was validated with one set of experimental data. The CFD model was also used to predict key reaction variable distributions such as the temperature and the component concentrations in the reactor.     

混凝土黏聚开裂模型若干进展

黏聚模型是用来描述混凝土断裂行为的基本模型, 首先介绍了混凝土的黏聚开裂模型的基本概念,总结了确定黏聚区的本构方程的各种方法,即直接单轴拉伸测试、J积分方法、R曲线法、柔度法和逆推法.然后介绍了黏聚模型在I型和复合型裂纹问题、疲劳断裂问题中的应用以及黏聚模型与混凝土尺寸效应的关系.最后对黏聚开裂模型与桥联模型、带状裂缝模型进行了比较和总结, 指出了该模型存在的问题, 并对其以后的发展方向提出了建议.      

Different effects of economic and structural performance indexes on model construction of structural topology optimization

The objective and constraint functions related to structural optimization designs are classified into economic and performance indexes in this paper.The influences of their different roles in model construction of structural topology optimization are also discussed.Furthermore,two structural topology optimization models,optimizing a performance index under the limitation of an economic index,represented by the minimum compliance with a volume constraint(MCVC)model,and optimizing an economic index under the limitation of a performance index,represented by the minimum weight with a displacement constraint(MWDC)model,are presented.Based on a comparison of numerical example results,the conclusions can be summarized as follows:(1)under the same external loading and displacement performance conditions,the results of the MWDC model are almost equal to those of the MCVC model;(2)the MWDC model overcomes the difficulties and shortcomings of the MCVC model;this makes the MWDC model more feasible in model construction;(3)constructing a model of minimizing an economic index under the limitations of performance indexes is better at meeting the needs of practical engineering problems and completely satisfies safety and economic requirements in mechanical engineering,which have remained unchanged since the early days of mechanical engineering.     

Multi-layer analytic solution for k-ω model equations via a symmetry approach

Despite being one of the oldest and most widely-used turbulence models in engineering computational fluid dynamics (CFD), the k-ω model has not been fully understood theoretically because of its high nonlinearity and complex model parameter setting. Here, a multi-layer analytic expression is postulated for two lengths (stress and kinetic energy lengths), yielding an analytic solution for the k-ω model equations in pipe flow. Approximate local balance equations are analyzed to determine the key parameters in the solution, which are shown to be rather close to the empirically-measured values from the numerical solution of the Wilcox k-ω model, and hence the analytic construction is fully validated. The results provide clear evidence that the k-ω model sets in it a multilayer structure, which is similar to but different, in some insignificant details, from the Navier-Stokes (N-S) turbulence. This finding explains why the k-ω model is so popular, especially in computing the near-wall flow. Finally, the analysis is extended to a newly-refined k-ω model called the structural ensemble dynamics (SED) k-ω model, showing that the SED k-ω model has improved the multi-layer structure in the outer flow but preserved the setting of the k-ω model in the inner region.     

A rheology model of high damping rubber bearings for seismic analysis: Identification of nonlinear viscosity

The mechanical behavior of high damping rubber bearings (HDRBs) is investigated under horizontal cyclic shear deformation with a constant vertical compressive load. On the basis of experimental observations, an elasto-viscoplastic rheology model of HDRBs for seismic analysis is developed. In this model, the Maxwell model is extended by adding a nonlinear elastic spring and an elasto-plastic model (spring-slider) in parallel. In order to identify constitutive relations of each element in the rheology model, an experimental scheme comprised of three types of tests, namely a cyclic shear (CS) test, a multi-step relaxation (MSR) test, and a simple relaxation (SR) test, are carried out at room temperature. HDRB specimens with the standard ISO geometry and three different high damping rubber materials are employed in these tests. A nonlinear viscosity law of the dashpot in the Maxwell model is deduced from the experimental scheme, and incorporated into the rheology model to reproduce the nonlinear rate dependent behavior of HDRBs. Finally, numerical simulation results for sinusoidal loading are presented to illustrate capability of the proposed rheology model in reproducing the mechanical behavior of HDRBs.     

Spatial artificial neural network model for subgrid-scale stress and heat flux of compressible turbulence

The subgrid-scale(SGS) stress and SGS heat flux are modeled by using an artificial neural network(ANN) for large eddy simulation(LES) of compressible turbulence. The input features of ANN model are based on the first-order and second-order derivatives of filtered velocity and temperature at different spatial locations. The proposed spatial artificial neural network(SANN)model gives much larger correlation coefficients and much smaller relative errors than the gradient model in an a priori analysis. In an a posteriori analysis, the SANN model performs better than the dynamic mixed model(DMM) in the prediction of spectra and statistical properties of velocity and temperature, and the instantaneous flow structures.     

Modeling of the mode dynamics generated by Madison Symmetric Torus machine utilizing a modified sine-Gordon equation

In this paper, a new dynamic model is presented for the experimental data generated by the Madison Symmetric Torus (MST) machine. The model is based on a modified sine-Gordon (SG) dynamic equation. The modified sine-Gordon equation model effectively captures the behavior of the slinky mode in reversed-field pinch experiments. In addition, this paper demonstrates how the derived model accurately describes the behavior of the localized magnetohydrodynamic mode (slinky mode) that appears in reversed-field pinch toroidal magnetic confinement systems. The modified SG equation model is solved analytically by using the perturbation method. The resulting model is fit to match a variety of experimental results in the MST reversed-field pinch experiment. The efficacy of the newly developed model in effectively representing the slinky mode is verified by comparing obtained analytical solution to experimentally measured data.     

Numerical Modeling of Hydraulic Hysteresis in Unsaturated Soils

This article presents the implementation of the constitutive model of Wheeler (Geotechnique 53(1):41–54, 2003) for coupling of hydraulic hysteresis and mechanical behavior of unsaturated soils in a fully coupled transient hydro-mechanical finite element (FE) model (computer code UNSATEX) developed by the authors. The constitutive model considers the effects of irreversible changes of degree of saturation on stress–strain behavior and the influence of plastic volumetric strains on the water retention behavior. The mathematical framework and the numerical implementation of the constitutive model are presented and discussed. The FE model is verified and validated against analytical predictions [obtained using the model of Wheeler (Geotechnique 53(1):41–54, 2003] as well as experimental results from the literature involving unsaturated soils undergoing various combinations of drying, wetting, loading, unloading, and reloading paths. Comparison of the results shows that the developed FE model can be used to predict various aspects of the behavior of unsaturated soils under drying and wetting as well as loading and unloading paths. The merits and limitations of the FE model are highlighted.     

基于最优加权组合模型及高斯－牛顿法的滑坡变形预测研究

在对最优加权组合理论和高斯－牛顿法优化非线性模型参数的方法研究的基础上,依托于洒勒山滑坡的实际变形监测资料,建立了该滑坡变形预测的3个非线性预测模型：指数模型、Verhulst模型和灰色GM（1,1）模型;利用最优加权组合理论建立了洒勒山滑坡的最优加权组合预测模型,并运用高斯－牛顿法对各单一模型和组合模型的参数进行了优化。通过对比分析得出：组合模型的预测精度高于任何单一模型的预测精度;参数优化后各单一模型的预测精度都有不同程度的提高;参数优化后的组合模型预测精度是最高的。因此,综合运用最优组合理论和高斯－牛顿法处理滑坡预测预报模型,是提高滑坡预测预报精度的行之有效的方法。     

A crack-mechanics based model for damage and plasticity of brittle materials under dynamic loading

A rate-dependent model for damage and plastic deformation of brittle materials under dynamic loading is presented. The model improves upon a recently developed micromechanical damage model (Zuo et al., 2006) by incorporating plastic deformation of the material. The distribution of the microcracks in the material is assumed to remain isotropic, and the damage evolution is through the growth of the average crack size. Plasticity is considered through an additive decomposition of the total strain rate, and a rate-independent, von Mises model is used. The model was applied to simulate the response of a model material (SiC) under uniaxial strain loading. To further examine the behavior of the model, cyclic loading and large-strain compressive loading were considered. Numerical results of the model predictions are presented, and comparisons with those from a previous model are provided.     

A new car-following model with consideration of inter-vehicle communication

In this paper, we construct a new car-following model with inter-vehicle communication (IVC) to study the driving behavior under an accident. The numerical results show that the proposed model can qualitatively describe the effects of IVC on each vehicle’s speed, acceleration, movement trail, and headway under an accident and that the new model can overcome the full velocity difference (FVD) model’s shortcoming that collisions occur under an accident, which illustrates that the new model can better describe the driving behavior under an accident than the FVD model.     

双材料界面断裂力学模型与实验方法

纤维增强聚合物(FRP)质轻、高强, 可提高结构的刚度、强度、抗震性能和耐久性, 近年来在结构加固及工程改造中得到广泛应用. FRP与传统复合材料之间形成双材料黏结界面, 界面断裂特性是决定双材料结构性能的关键因素. 对双材料界面裂纹尖端应力场理论、界面裂纹模型、黏结界面I型、II型及混合型断裂试验及理论研究现状进行综合评述和分析. 界面模型主要有经典梁/板理论和刚性节点模型、考虑剪切变形的双亚层理论和半刚性节点模型、基于双亚层理论的柔性节点模型、考虑剪切变形的多层亚层理论和多亚层柔性节点模型、弹性地基梁模型以及黏聚模型. 还介绍了双材料界面断裂力学在FRP-混凝土研究中的应用.