Global response sensitivity analysis using probability distance measures and generalization of Sobol's analysis

The study introduces two new alternatives for global response sensitivity analysis based on the application of the L-2-norm and Hellinger's metric for measuring distance between two probabilistic models. Both the procedures are shown to be capable of treating dependent non-Gaussian random variable models for the input variables. The sensitivity indices obtained based on the L2-norm involve second order moments of the response, and, when applied for the case of independent and identically distributed sequence of input random variables, it is shown to be related to the classical Sobol's response sensitivity indices. The analysis based on Hellinger's metric addresses variability across entire range or segments of the response probability density function. The measure is shown to be conceptually a more satisfying alternative to the Kullback-Leibler divergence based analysis which has been reported in the existing literature. Other issues addressed in the study cover Monte Carlo simulation based methods for computing the sensitivity indices and sensitivity analysis with respect to grouped variables. Illustrative examples consist of studies on global sensitivity analysis of natural frequencies of a random multi-degree of freedom system, response of a nonlinear frame, and safety margin associated with a nonlinear performance function. (C) 2015 Elsevier Ltd. All rights reserved.

Ag-doped hydroxyapatite as efficient adsorbent for removal of Congo red dye from aqueous solution: Synthesis, kinetic and equilibrium adsorption isotherm analysis

In the present study a versatile and efficient adsorbent with high adsorption capacity for adsorption of Congo red dye in aqueous solution at ambient temperature without adjusting any pH is presented over the Ag modified calcium hydroxyapatite (CaHAp). CaHAp and Ag-doped CaHAp materials were synthesized using facile aqueous precipitation method. The physico-chemical properties of the materials were determined by powder X-ray diffraction (XRD), Fourier transform infrared (FT-IR) spectroscopy, Transmission electron microscopy (TEM), UV-Visible spectroscopy, N-2 physisorption and acidity was determined by n-butylamine titration and pyridine adsorption methods. XRD analysis confirmed all adsorbents exhibit hexagonal CaHAp structure with P6(3)/m space group. TEM analysis confirms the rod like morphology of the adsorbents and the average length of the rods were in the range of 40-45 nm. Pyridine adsorption results indicate increase in number of Lewis acid sites with Ag doping in CaHAp. Adsorption capacity of CaHAp was found increased with Ag content in the adsorbents. Ag (10): CaHAp adsorbent showed superior adsorption performance among all the adsorbents for various concentrations of Congo red (CR) dye in aqueous solutions. The amount of CR dye adsorbed on Ag (10): CaHAp was found to be 49.89-267.81 mg g(-1) for 50-300 ppm in aqueous solution. A good correlation between adsorption capacity and acidity of the adsorbents was observed. The adsorption kinetic data of adsorbents fitted well with pseudo second-order kinetic model with correlation coefficients ranged from 0.998 to 0.999. The equilibrium adsorption data was found to best fit to the Langmuir adsorption isotherm model. (C) 2015 Elsevier Inc. All rights reserved.

Prediction of electronically nonadiabatic decomposition mechanisms of isolated gas phase nitrogen-rich energetic salt: Guanidium-triazolate

Electronically nonadiabatic decomposition pathways of guanidium triazolate are explored theoretically. Nonadiabatically coupled potential energy surfaces are explored at the complete active space self-consistent field (CASSCF) level of theory. For better estimation of energies complete active space second order perturbation theories (CASPT2 and CASMP2) are also employed. Density functional theory (DFT) with B3LYP functional and MP2 level of theory are used to explore subsequent ground state decomposition pathways. In comparison with all possible stable decomposition products (such as, N-2, NH3, HNC, HCN, NH2CN and CH3NC), only NH3 (with NH2CN) and N-2 are predicted to be energetically most accessible initial decomposition products. Furthermore, different conical intersections between the S-1 and S-0 surfaces, which are computed at the CASSCF(14,10)/6-31G(d) level of theory, are found to play an essential role in the excited state deactivation process of guanidium triazolate. This is the first report on the electronically nonadiabatic decomposition mechanisms of isolated guanidium triazolate salt. (C) 2015 Elsevier B.V. All rights reserved.

Electromechanical analysis of piezoelectric bimorph actuator in static state considering the nonlinearity at high electric field

This article contains electromechanical analysis of a piezoelectric bimorph actuator at high electric field by incorporating second-order constitutive equations of piezoelectric material. Tip deflection, block force, block moment, block load, output strain energy, output energy density, input electrical energy, and energy efficiency are analytically derived for the actuator at high electric field. The analysis shows that output energy and energy density increase more rapidly at high electric field, compared to the prediction by the linear model. The analysis shows energy efficiency depends on electric field. Some analytical results are validated with the published experimental results.

A peridynamic theory for linear elastic shells

A state-based peridynamic formulation for linear elastic shells is presented. The emphasis is on introducing, possibly for the first time, a general surface based peridynamic model to represent the deformation characteristics of structures that have one geometric dimension much smaller than the other two. A new notion of curved bonds is exploited to cater for force transfer between the peridynamic particles describing the shell. Starting with the three dimensional force and deformation states, appropriate surface based force, moment and several deformation states are arrived at. Upon application on the curved bonds, such states yield the necessary force and deformation vectors governing the motion of the shell. By incorporating a shear correction factor, the formulation also accommodates analysis of shells that have higher thickness. In order to attain this, a consistent second order approximation to the complementary energy density is considered and incorporated in peridynamics via constitutive correspondence. Unlike the uncoupled constitution for thin shells, a consequence of a first order approximation, constitutive relations for thick shells are fully coupled in that surface wryness influences the in-plane stress resultants and surface strain the moments. Our proposal on the peridynamic shell theory is numerically assessed against simulations on static deformation of spherical and cylindrical shells, that of flat plates and quasi-static fracture propagation in a cylindrical shell. (C) 2016 Elsevier Ltd. All rights reserved.

Perturbational Finite Volume Method For The Solution of 2-D Navier-Stokes Equations on Unstructured and Structured Colocated Meshes

Based on the first-order upwind and second-order central type of finite volume( UFV and CFV) scheme, upwind and central type of perturbation finite volume ( UPFV and CPFV) schemes of the Navier-Stokes equations were developed. In PFV method, the mass fluxes of across the cell faces of the control volume (CV) were expanded into power series of the grid spacing and the coefficients of the power series were determined by means of the conservation equation itself. The UPFV and CPFV scheme respectively uses the same nodes and expressions as those of the normal first-order upwind and second-order central scheme, which is apt to programming. The results of numerical experiments about the flow in a lid-driven cavity and the problem of transport of a scalar quantity in a known velocity field show that compared to the first-order UFV and second-order CFV schemes, upwind PFV scheme is higher accuracy and resolution, especially better robustness. The numerical computation to flow in a lid-driven cavity shows that the under-relaxation factor can be arbitrarily selected ranging from 0.3 to 0. 8 and convergence perform excellent with Reynolds number variation from 102 to 104.

双曲守恒型方程的二阶摄动有限差分格式

对双曲守恒型方程，将其一阶迎风格式空间差商的常系数摄动展开为时间步长和空间步长的幂级数，通过确定幂级数系数而获得二阶精度的摄动有限差分（PFD）格式。进而从双曲守恒型方程的通量分裂型一阶迎风格式出发，通过娄似的摄动展开方法，获得空间精度为二阶的通量分裂形式的摄动有限差分（FPFD）格式。这两类格式保留了一阶守恒迎风格式的简洁结构形式，使用三节点即可达到二阶精度，又避免了三点二阶格式的非物理数值振荡。并将这两类格式推广应用到双曲守恒型方程组，最后通过模型方程和一维激波管流动的数值算例验证了格式的高精度、高分辨率性质。

Thermocapillary Effect in a Thermal Rheological Jet

The process of die swell in polymer jets is an important feature within polymer processing and can be explained through a study of its rheological effects. The existence of a thermocapillary effect, driven by the gradient of its surface tension, should be considered when examining a thermal jet that has a non-uniform temperature distribution on its free surface, as in various polymer processings. Both the rheological effect and thermocapillary effect on die swell can be studied numerically through a finite element method as used on a two-dimensional and unsteady model, in which a Coleman-Noll second-order fluid model is employed. The results show that the expanding angle depends on both the rheological property of the fluid and the pressure at the vessel exit. Although both the thermocapillary and the rheological effects contribute to the cross-section expansion of the fluid jet, the latter is more important in determining the expansion.

Navier-Stokes方程二阶速度滑移边界条件的检验

对微尺度气体流动，Navier-Stokes方程和一阶速度滑移边界条件的结果与实验数据相比，在滑移区相互符合，在过渡领域则显著偏离，为改善Navier-Stokes方程在过渡领域的表现，有些研究者尝试引入二阶速度滑移边界条件，如Cercignani模型，Deissler模型和Beskok-Karniadakis模型．以微槽道气体流动为例，将Navier-Stokes方程在不同的二阶速度滑移模型下的结果与动理论的直接模拟Monte Carlo（DSMC）方法和信息保存（IP）方法以及实验数据进行比较．在所考察的3种具有代表性的二阶速度滑移模型中，Cercignani模型表现最好，其所给出的质量流率在Knudsen数为0．4时仍与DSMC和IP结果相符；然而，细致比较表明，Cercignani模型给出的物面滑移速度及其附近的速度分布在滑流区和过渡领域的分界处（Kn=0．1）已明显偏离DSMC和IP的结果．

一个基于非均匀网格上的二阶基本无振荡差分格式

在基本无振荡格式的构造中，将通常的对流通量f的逼近方式推广到对通量导数的逼近，这一构造方法可以有效地应用到非均匀或非结构网格。直接基于非均匀网格上，构造了一个二阶的基本无振荡(ENO)差分格式。该格式具有形式简单，对网格的划分灵活，与传统格式相比不增加计算量等优点，儿个数值算例证明了格式的有效性。

Simulation of Coal Combustion by AUSM Turbulence-Chemistry Char Combustion Model and a Full Two-Fluid Model

An algebraic unified second-order moment (AUSM) turbulence-chemistry model of char combustion is introduced in this paper, to calculate the effect of particle temperature fluctuation on char combustion. The AUSM model is used to simulate gas-particle flows, in coal combustion in a pulverized coal combustor, together with a full two-fluid model for reacting gas-particle flows and coal combustion, including the sub-models as the k-epsilon-k(p) two-phase turbulence niodel, the EBU-Arrhenius volatile and CO combustion model, and the six-flux radiation model. A new method for calculating particle mass flow rate is also used in this model to correct particle outflow rate and mass flow rate for inside sections, which can obey the principle of mass conservation for the particle phase and can also speed up the iterating convergence of the computation procedure effectively. The simulation results indicate that, the AUSM char combustion model is more preferable to the old char combustion model, since the later totally eliminate the influence of particle temperature fluctuation on char combustion rate.

Studies of the Effect of a Coal Concentrator on NO Formation in Swirling Coal Combustion

To develop low-pollution burners, the effect of a coal concentrator on NO formation in swirling coal combustion is studied using both numerical simulation and experiments. The isothermal gas-particle two-phase velocities and particle concentration in a cold model of swirl burners with and without coal concentrators were measured using the phase Doppler particle anemometer (PDPA). A full two-fluid model of reacting gas-particle flows and coal combustion with an algebraic unified second-order moment (AUSM) turbulence-chemistry model for the turbulent reaction rate of NO formation are used to simulate swirling coal combustion and NO formation with different coal concentrators. The results give the turbulent kinetic energy, particle concentration, temperature and NO concentration in cases of with and without coal concentrators. The predicted results for cold two-phase flows are in good agreement with the PDPA measurement results, showing that the coal concentrator increases the turbulence and particle concentration in the recirculation zone. The combustion modeling results indicate that although the coal concentrator increases the turbulence and combustion temperature, but still can remarkably reduce the NO formation due to creating high coal concentration in the recirculation zone.

An Improved Ce/Se Scheme And Its Application To Detonation Propagation

An improved two-dimensional space-time conservation element and solution element ( CE/ SE) method with second-order accuracy is proposed, examined and extended to simulate the detonation propagations using detailed chemical reaction models. The numerical results of planar and cellular detonation are compared with corresponding results by the Chapman-Jouguet theory and experiments, and prove that the method is a new reliable way for numerical simulations of detonation propagation.

Quasi-Dammann Grating With Proportional Intensity Array Spots

A quasi-Dammann grating is proposed to generate array spots with proportional-intensity orders in the far field. To describe the performance of the grating, the uniformities of the array spots are redefined. A two-dimensional even-sampling encode scheme is adopted to design the quasi-Dammann grating. Numerical solutions of the binary-phase quasi-Dammann grating with proportional-intensity orders are given. The experimental results with a third-order quasi-Dammann grating, which has an intensity proportion of 3:2:1 from zero order to second order, are presented. (C) 2008 Optical Society of America

Hydrodynamic interaction between two vertical cylinders in water waves

The hydrodynamic interaction between two vertical cylinders in water waves is investigated based on the linearized potential flow theory. One of the two cylinders is fixed at the bottom while the other is articulated at the bottom and oscillates with small amplitudes in the direction of the incident wave. Both the diffracted wave and the radiation wave are studied in the present paper. A simple analytical expression for the velocity potential on the surface of each cylinder is obtained by means of Graf's addition theorem. The wave-excited forces and moments on the cylinders, the added masses and the radiation damping coefficients of the oscillating cylinder are all expressed explicitly in series form. The coefficients of the series are determined by solving algebraic equations. Several numerical examples are given to illustrate the effects of various parameters, such as the separation distance, the relative size of the cylinders, and the incident angle, on the first-order and steady second-order forces, the added masses and radiation-damping coefficients as well as the response of the oscillating cylinder.