A multi-moment finite volume formulation for shallow water equations on unstructured mesh

A novel and accurate finite volume method has been presented to solve the shallow water equations on unstructured grid in plane geometry. In addition to the volume integrated average (VIA moment) for each mesh cell, the point values (PV moment) defined on cell boundary are also treated as the model variables. The volume integrated average is updated via a finite volume formulation, and thus is numerically conserved, while the point value is computed by a point-wise Riemann solver. The cell-wise local interpolation reconstruction is built based on both the VIA and the PV moments, which results in a scheme of almost third order accuracy. Efforts have also been made to formulate the source term of the bottom topography in a way to balance the numerical flux function to satisfy the so-called C-property. The proposed numerical model is validated by numerical tests in comparison with other methods reported in the literature. (C) 2010 Elsevier Inc. All rights reserved.

A global shallow water model using high order multi-moment constrained finite volume method and icosahedral grid

A novel accurate numerical model for shallow water equations on sphere have been developed by implementing the high order multi-moment constrained finite volume (MCV) method on the icosahedral geodesic grid. High order reconstructions are conducted cell-wisely by making use of the point values as the unknowns distributed within each triangular cell element. The time evolution equations to update the unknowns are derived from a set of constrained conditions for two types of moments, i.e. the point values on the cell boundary edges and the cell-integrated average. The numerical conservation is rigorously guaranteed. in the present model, all unknowns or computational variables are point values and no numerical quadrature is involved, which particularly benefits the computational accuracy and efficiency in handling the spherical geometry, such as coordinate transformation and curved surface. Numerical formulations of third and fourth order accuracy are presented in detail. The proposed numerical model has been validated by widely used benchmark tests and competitive results are obtained. The present numerical framework provides a promising and practical base for further development of atmospheric and oceanic general circulation models. (C) 2009 Elsevier Inc. All rights reserved.

耗散粒子动力学处理复杂固体壁面的一种有效方法

耗散粒子动力学(dissipative particle dynamics,DPD)作为一种介观尺度拉格朗日型粒子方法,已经成功地应用于微纳米流动和生化科技的研究中. 复杂固体壁面的处理和壁面边界条件的实施一直是DPD方法发展及应用的一个障碍. 提出了处理复杂固体壁面的一种新的方法. 复杂固体区域通过冻结随机分布并且达到平衡状态的DPD粒子代表;所冻结的DPD粒子位于临近流动区域的一个截距内;在靠近固体壁面的流动区域中设置流动反弹层,当流动DPD粒子进入此流动层后反弹回流动区域. 应用这种固体壁面处理方法.

Influence of Contact Geometry on Hardness Behavior in Nano-indentation

In this paper the influence of contact geometry, including the round tip of the indenter and the roughness of the specimen, on hardness behavior for elastic plastic materials is studied by means of finite element simulation. We idealize the actual indenter by an equivalent rigid conic indenter fitted smoothly with a spherical tip and examine the interaction of this indenter with both a flat surface and a rough surface. In the latter case the rough surface is represented by either a single spherical asperity or a dent (cavity). Indented solids include elastic perfectly plastic materials and strain hardening elastic-plastic materials, and the effects of the yield stress and strain hardening index are explored. Our results show that due to the finite curvature of the indenter tip the hardness versus indentation depth curve rises or drops (depending on the material properties of the indented solids) as the indentation depth decreases, in qualitative agreement with experimental results. Surface asperities and dents of curvature comparable to that of the indenter tip can appreciably modify the hardness value at small indentation depth. Their effects would appear as random variation in hardness.

Energy-level shifts of a stationary hydrogen atom in a static external gravitational field with Schwarzschild geometry

he first order perturbations of the energy levels of a stationary hydrogen atom in a static external gravitational field, with Schwarzschild metric, are investigated. The energy shifts are calculated for the relativistic 1S, 2S, 2P, 3S, 3P, 3D, 4S, 4P, 4D, and 4F levels. The results show that the energy-level shifts of the states with total angular momentum quantum number 1/2 are all zero, and the ratio of absolute energy shifts with total angular momentum quantum number 5/2 is 145. This feature can be used to help us to distinguish the gravitational effect from other effects.

The geometry of the NO2- anion: ab initio calculations and Franck-Condon analysis

Geometry optimization and harmonic vibrational frequency calculations have been performed on the (X) over bar (2)A(1) state of NO2 and (X) over bar (1)A(1) state of NO2-. Franck-Condon analyses and spectral simulations were carried out on the NO2((X) over bar (2)A(1))-NO2-((X) over bar (1)A(1)) photo detachment process. In addition, the equilibrium geometry parameters, r(NO)= 1.248 +/- 0.005 Angstrom and angle(ONO) 116.8 +/- 0.5degrees, of the (X) over bar (1)A(1) state of NO2-, are derived by employing an iterative Franck-Condon analysis procedure in the spectral simulation. Our conclusions regarding the anion geometry suggest a reinterpretation of the results of Woo et al. (C) 2004 Published by Elsevier B.V.

Computational study on the addition of HCN to methanimine catalyzed by formamidine and formamide

In this paper, we have carried out a theoretical study on the addition of HCN to methanimine with formamidine or formamide using second-order Moller-Plesset perturbation (MP2) method with 6-31 + G(d,p) basis sets. At MP2 level. a high-energy, intermediate has been located for each pathway. The addition of HCN to methanimine with formamidine has the lowest free energy barrier according to the calculations at MP2 level. (C) 2004 Elsevier B.V. All rights reserved.

Stability shear stress and equilibrium cross-sectional geometry of sheltered tidal channels

This study relates tidal channel cross-sectional area (A) to peak spring discharge (Q) via a physical mechanism, namely the stability shear stress ( tau sub(S)) just necessary to maintain a zero gradient in net along-channel sediment transport. It is assumed that if bed shear stress ( tau ) is greater than tau sub(S), net erosion will occur, increasing A, and reducing tau similar to (Q/A) super(2) back toward tau sub(S). If tau < tau sub(S) there will be net deposition, reducing A and increasing tau toward tau sub(S). A survey of the literature allows estimates of Q and A at 242 sections in 26 separate sheltered tidal systems. Assuming a single value of tau sub(S) characterizes the entire length of a given tidal channel, it is predicted that along-channel geometry will follow the relation Ah sub(R) super(1) super(/) super(6) similar to Q. Along-channel regressions of the form Ah sub(R) super(1) super(/) super(6) similar to Q super( beta ) give a mean observed value for beta of 1.00 plus or minus 0.06, which is consistent with this concept. Results indicate that a lower bound on tau sub(S) (and an upper bound on A) for stable channels is provided by the critical shear stress ( tau sub(C)) just capable of initiating sediment motion. Observed tau sub(S) is found to vary among all systems as a function of spring tidal range (R sub(sp)) according to the relation tau sub(S) approximately 2.3 R sub(sp) super(0.79) tau sub(C). Observed deviations from uniform tau sub(S) along individual channels are associated with along-channel variation in the direction of maximum discharge (i.e., flood-versus ebb-dominance).

a computational proof of complexity of some restricted counting problems

South Central University

SAXS measurements of the interface in polyacrylate and epoxy interpenetrating networks with fractal geometry

The interface thickness in two-component interpenetrating polymer networks (IPN) system based on polyacrylate and epoxy were determined using small-angle X-ray scattering (SAXS) in terms of the theory proposed by Ruland. The thickness was found to be nonexistent for the samples at various compositions and synthesized at variable conditions-temperature and initiator concentration. By viewing the system as a two-phase system with a sharp boundary, the roughness of the interface was described by fractal dimension, D, which slightly varies with composition and synthesis condition. Length scales in which surface fractals are proved to be correct exist for each sample and range from 0.02 to 0.4 Angstrom(-1). The interface in the present IPN system was treated as fractal, which reasonably explained the differences between Pored's law and experimental data, and gained an insight into the interaction between different segments on the interface. (C) 1997 Elsevier Science Ltd.

Quantitative structure-property relationships for colour reagents and their colour reactions with cerium using computational neural networks

Quantitative structure-activity/property relationships (QSAR/QSPR) studies have been exploited extensively in the designs of drugs and pesticides, but few such studies have been applied to the design of colour reagents. In this work, the topological indices A(x1)-A(x3) suggested in this laboratory were applied to multivariate analysis in structure-property studies. The topological indices of 43 phosphone bisazo derivatives of chromotropic acid were calculated. The structure-property relationships between colour reagents and their colour reactions with cerium were studied using A(x1-Ax3) indices with satisfactory results. The purpose of this work was to establish whether QSAR can be used to predict the contrasts of colour reactions and in the longer term to be a helpful tool in colour reagent design.

Quantitative structure-property relationships for colour reagents and their colour reactions with ytterbium using regression analysis and computational neural networks

In this paper, the new topological indices A(x1)-A(x3) suggested in our laboratory and molecular connectivity indices have been applied to multivariate analysis in structure-property studies. The topological indices of twenty asymmetrical phosphono bisazo derivatives of chromotropic acid have been calculated. The structure-property relationships between colour reagents and their colour reactions with ytterbium have been studied by A(x1)-A(x3) indices and molecular connectivity indices with satisfactory results. Multiple regression analysis and neural networks were employed simultaneously in this study.

QUANTITATIVE STRUCTURE-ACTIVITY-RELATIONSHIPS FOR TOXICITY OF PHENOLS USING REGRESSION-ANALYSIS AND COMPUTATIONAL NEURAL NETWORKS

Quantitative structure-toxicity models were developed that directly link the molecular structures of a et of 50 alkYlated and/or halogenated phenols with their polar narcosis toxicity, expressed as the negative logarithm of the IGC50 (50% growth inhibitor