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.     

汽车液压阻尼器的四段线性化模型

通过对现有某车用液压阻尼器的实验研究和车辆悬架系统动力学仿真研究，发现这种阻尼器的动力学特性存在复杂的非线性行为。为准确描述这类阻尼器的动力学特性和进行汽车悬架系统的动力学仿真研究，本文根据阻尼器的实验数据，提出一种新型四段线性化阻尼器模型，并建立其相应的数学模型。通过汽车悬架系统的动力学仿真研究表明，采用四段线性化阻尼器模型，汽车悬架系统簧载质量在不同路面不平度函数作用下的加速度响应和均方根加速度响应随频率变化不会出现大的波动。与等效线性化阻尼器模型相比，四段线性化阻尼器模型，能较好降低汽车悬架系统簧载质量的加速度响应和均方根加速度响应，有利于提高汽车乘坐舒适性和行驶稳定性，并为新型阻尼器设计提供一个新的思路。     

Photoelastic model measurement with rotated principal axes by scattered-light photoelasticity

The method of measurement of a three-dimensional photoelastic model with rotated principal axes has not yet been fully experimentally established. It is known that a three-dimensional photoelastic model can be reduced to an optically equivalent model. In this paper, the optically equivalent model is realized from a stratified model consisting of two layers of the frozen stress model. The secondary principal stress direction and the relative phase retardation of the frozen disk model in the stratified frozen stress model are determined for the entire field from Stokes parameters obtained by scattered-light photoelasticity using unpolarized light. The accuracy of these values is confirmed by a comparison with results only from the frozen stress disk model.     

Principles of Integrated Flow Modeling

Integrated flow modeling is the combination of a traditional flow simulator with a petrophysical model. By combining a petrophysical model with a traditional flow model, it is possible to perform calculations that improve our ability to monitor fluid movement in porous media. This paper outlines the formulation of an integrated flow model IFLO and its multi-variable, Newton–Raphson IMPES solution procedure. The benefits of integrated flow modeling and the underlying principles involved in the integration of a flow model with a petrophysical model are presented. Results from the IFLO model are used to illustrate the principles.     

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.     

Solving parametric complex fluids models in rheometric flows

Inverse identification of complex fluid behaviors is a tricky task because sometimes there are several rheological parameters and the identification procedure itself is quite expensive from the computational time viewpoint. Standard inverse identification procedures solve the model for a choice of the model parameters, and then parameters are updated trying to minimize the gap between the model predictions and some available experimental measures. Thus, the model has to be evaluated for each trial set of the model parameters. When models involve a great number of degrees of freedom the identification procedure becomes a computationally expensive task. In this paper we propose a new procedure able to solve once the model for any value of the model parameters. For this purpose, all the model parameters are considered as extra-coordinates of the model. Thus, the model results finally defined in a multidimensional space including the physical space x, the time t and a number of extra-coordinates related to the model parameters. The solution of such model needs for circumventing the curse of dimensionality illness that suffer multidimensional models.     

Grid-averaged Lagrangian LES model for multiphase turbulent flow

A grid-averaged Lagrangian (GAL) model for dispersed particle motion in multiphase turbulent flow is presented to provide a large eddy simulation (LES) model for multiphase turbulent flow in which a quite large number of particles are involved. The GAL model is based on an averaging operation for a Lagrangian-type equation of motion of a particle over a computational grid volume and a procedure of reallocation of a dispersed particle cloud with its centroid movement to each grid. The model is therefore a mixed Eulerian–Lagrangian model which can effectively reduce computational time compared with existing Lagrangian-type models, without losing the advantage of Lagrangian-type models that they can properly describe the dynamical evolution of particles. Since the GAL model adopts the grid-volume averaging operation it can easily provide an effective SGS model for LES modeling of multiphase turbulent flow. The validity of the multiphase LES model developed, which is named the GAL-LES model, is confirmed through its application to a particle plume, in which the present model is found to simulate large-eddy motion usually observed in a jet and plume, and to give good agreements with experimental data.     

Triphasic FE Modeling of the Skin Water Barrier

A skin–air model is developed to model the water barrier function of skin. The skin model is a porous solid saturated with a monovalent salt solution. The air model is a vapor diffusion model in non-moving air. In vivo measurements of water loss from human skin under varying ambient conditions are used to validate the model.     

Extreme wind pressure estimation based on the r largest order statistics mode

提出一种基于r个最大次序统计量(r largest order statistics, r-LOS)模型的极 值风压估算方法, 该模型包括广义极值(generalized extreme value distribution, GEV)联合分布形式及Gumbel联合分布形式. 提出了独立风压峰值 r-LOS序列构造方法、最优r值确定方法、r-LOS GEV模型和r-LOS Gumbel模型的优选方法. 将r-LOS模型方法应用于某低矮工业建筑刚性模型测压试验的极 值风压估算. 当采用多段风压时程估算极值风压时, r-LOS Gumbel模型优于r-LOS GEV模型和经典Gumbel模型. 当采用单段风压时程时, 与基于改进Hermite模型的峰值因子 法、Sadek-Simiu法相比, r-LOS Gumbel模型方法能更加精确地估算极值风压, 适用于非高斯风压的极值估算, 可以给出极值风压分位点的解析解. 分析表明, r-LOS Gumbel模型方法是一种多段时程和单段时程条件下均适用的极值风压估算方法.     

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

 
滑坡预报模型主要用来预报滑坡发生的时间,但是如何检验效果是滑坡预报者和决策者面临的一个重要课题。本文通过提出滑坡拟合效果指标（包括后验差指标、模型拟合效率指数和均方根误差）和试预报效果指标（包括试报效果指标和相关系数指标）,初步建立了滑坡的预报质量检验模型。依据建立的检验模型,针对链子崖危岩体监测资料分别采用灰色GM（1, 1）数学模型、三次指数平滑模型和时间序列模型进行预报。预报结果表明：三种模型中,综合指标GM（1, 1）数值最高,建议该滑坡采用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.     

基于石膏材料的小结构重力坝损伤本构模型研究

大坝模型试验的模型材料力学性能是试验成功与否的关键影响因素。基于石膏模型材料基本力学性质和混凝土材料基本相似的实际情况,考虑石膏模型材料的孔隙率,在能量损伤理论模型的基础上,结合工程规范,建立了石膏模型材料的受拉损伤本构模型。同时,分析了模型材料应变率对材料峰值应力的影响,引入动应力提高系数,并给出了与应变率相关的模型材料受拉损伤本构模型。利用该本构模型模拟了动荷载下石膏模型重力坝的破坏过程,并与模型试验结果进行对比。结果表明：数值模拟得到的重力坝开裂破坏的位置、形式与模型试验结果较为接近。     

从非线性动力学的视角认识细长压杆的稳定性

工程上大部分机构和结构都处于动载的作用下，受压细长杆的失稳是复杂的动力破坏事件．应用Lagrange描述法建立了两端角铰支受压细长杆的非线性动力学模型，通过对这种模型简化分别得到非线性静力学模型、线性动力学模型和含三次非线性项的动力学模型．利用谱截断方法，讨论了线性动力学模型的局部分岔．通过讨论平衡态存在性和稳定性，得到了含三次非线性项的动力学模型分岔条件．研究表明，受压细长杆的非线性动力学模型中存在叉形分岔．     

基于最小概率DWO的激光陀螺温度误差模型辨识

鉴于激光陀螺温度误差模型的非线性和时变性,从提高模型辨识质量的需要出发,运用基于最小概率DWO（直接加权优化）的非线性系统辨识方法进行激光陀螺温度误差模型辨识研究,提出了一种精度更高的激光陀螺温度误差模型。该模型以带宽作为唯一需要确定的模型参数,避免了以往温度误差模型研究中的结构与参数辨识等复杂问题,从而在保证模型辨识质量的基础上,也相应提高了模型的适应能力,并通过两种温度误差特性有显著区别的激光陀螺（四频差动激光陀螺和二频机抖激光陀螺）的温度实验数据验证了该模型的正确性和适应性。     

回采工作面周期来压步距的灰色预测

利用灰色系统理论, 将回采工作面周期来压视为灾变, 仿线性GM(1, 1)预测模型建模机理, 本文建立了包括线性模型在内的统一非线性GM(1, 1)预报模型。利用非线性模型特例Verhulst模型和线性模型分别对两个不同矿山条件老顶周期来压步距进行预测, 结果表明, 不同的地质条件将适合不同的预测模型, 不能简单一概用线性模型进行预测。对于地质条件发生灾变的情形可建立激励预测模型对周期来压步距进行预测。     

Biased compensation recursive least squares-based threshold algorithm for time-delay rational models via redundant rule

This paper develops a biased compensation recursive least squares-based threshold algorithm for a time-delay rational model. The time-delay rational model is transformed into an augmented model by using the redundant rule, and then, a recursive least squares algorithm is proposed to estimate the parameters of the augmented model. Since the output of the augmented model is correlated with the noise, a biased compensation method is derived to eliminate the bias of the parameter estimates. Furthermore, based on the structures of the augmented model parameter vector and the rational model parameter vector, the unknown time delay can be computed by using a threshold given in prior. A simulated example is used to illustrate the efficiency of the proposed algorithm.     

Gravitational forces in dual-porosity systems: II. Computational validation of the homogenized model

Three models are considered for single component, single phase flow in naturally fractured porous media. The microscopic model holds on the Darcy scale, and it is considered to govern the system. The macroscopic, dual-porosity model was derived in Part I of this work from the microscopic model by two-scale mathematical homogenization. In this paper, we show that the dual-porosity model predicts well the behavior of the microscopic model by comparing their computed solutions in certain reasonable test cases. Homogenization gives a complex formula for a key parameter in the dual-porosity model; herein a simple approximation to this formula is presented. The third model considered is a single-porosity model with averaged parameters. It is shown that this type of model cannot predict the behavior of the microscopic flow.This work was supported in part by the National Science Foundation and by the State of Texas Governor's Energy Office.     

Numerical simulation of pollutant transport acted by wave for a shallow water sea bay

A numerical model of pollutant transport acted by water waves on a shallow‐water mild‐slope beach is established in this study. The numerical model is combined with a wave propagation model, a multiple wave‐breaking model, a wave‐induced current model and a pollutant convection–dispersion model. The wave propagation model is based on the higher‐order approximation of parabolic mild‐slope equation which can be used to simulate the wave refraction, diffraction and breaking in a large area of near‐shore zone combined with the wave‐breaking model. The wave‐induced current model is established using the concept of the radiation stress and considering the effect of bottom resistance caused by waves. The numerical model is verified by laboratory experiment results of regular and irregular waves over two mild beaches with different slopes. The numerical results agree well with experimental results. The numerical model has been applied in the near‐shore zone of Bohai bay in China. It is concluded that pollutant transport parallel to the shoreline due to the action of waves, which will induce serious pollution on the beach. Copyright © 2005 John Wiley & Sons, Ltd.     

复合材料层板基于剩余刚度比的剩余强度模型

本文从唯裂 学的观点提出了一个复合材料疲劳剩余刚度衰退模型,根据该模型导出了基于疲劳乘余刚度比的复合材料层合板疲劳剩余强度概率分布模型,给出了确定该模型中各参数的计算方法,根据该模型可预报复合材料层合板在某一应力水平的疲劳载荷作用下疲劳剩余刚度比为某一确定值时的疲劳剩余强度概率分布。为验证该模型的合理性,用典型复合材料层合板做了几组疲劳实验,并用实验数据对模型中的参数进行了估计,计算结果证明,该模型预报的疲劳剩余强度概率理论分布与实验结果符合较好。     

Aeroelastic prediction of the limit cycle oscillations of a cropped delta wing

The flutter and limit cycle oscillation (LCO) behavior of a cropped delta wing are investigated using a newly developed computational aeroelastic solver. This computational model includes a well-validated Euler finite difference solver coupled to a high-fidelity finite element structural solver. The nonlinear structural model includes geometric nonlinearities which are modelled using a co-rotational formulation. The LCOs of the cropped delta wing are computed and the results are compared to previous computations and to experiment. Over the range of dynamic pressures for which experimental results are reported, the LCO magnitudes computed using the current model are comparable to those from a previous computation which used a lower-order von Karman structural model. However, for larger dynamic pressures, the current computational model and the model which used the von Karman theory start to differ significantly, with the current model predicting larger deflections for a given dynamic pressure. This results in a LCO curve which is in better qualitative agreement with experiment. Flow features which were present in the previous computational model such as a leading edge vortex and a shock wave are enhanced in the current model due to the prediction of larger deflections and rotations at the higher dynamic pressures.