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.

A finite-element algorithm for modeling variably saturated flows

A general numerical algorithm in the context of finite element scheme is developed to solve Richards’ equation, in which a mass-conservative, modified head based scheme (MHB) is proposed to approximate the governing equation, and mass-lumping techniques are used to keep the numerical simulation stable. The MHB scheme is compared with the modified Picard iteration scheme (MPI) in a ponding infiltration example. Although the MHB scheme is a little inferior to the MPI scheme in respect of mass balance, it is superior in convergence character and simplicity. Fully implicit, explicit and geometric average conductivity methods are performed and compared, the first one is superior in simulation accuracy and can use large time-step size, but the others are superior in iteration efficiency. The algorithm works well over a wide variety of problems, such as infiltration fronts, steady-state and transient water tables, and transient seepage faces, as demonstrated by its performance against published experimental data. The algorithm is presented in sufficient detail to facilitate its implementation.

Design Study of a Heavily Loaded Ice Class Propeller Using an Advance Panel Method

A design and optimization procedure developed and used for a propeller installed on a twin-semitunnel-hull ship navigating in very shallow and icy water under heavy load conditions is presented. The base propeller for this vessel was first determined using classic design routines under open-water condition with existing model test data. In the optimization process, a panel method code (PROPELLA) was used to vary the pitch values and distributions and take into account the inflow wake distribution, tunnel gap, and cavitation effects. The optimized propeller was able to improve a ship speed of 0.02 knots higher than the desired speed and 0.06 knots higher than the classic B-series propeller. The analysis of the effect of inflow wake, hull tunnel, cavitation, and blade rake angle on propulsive performance is the focus of this paper.

Quantitative analysis analysis of element concentration in Al alloy by using laser-induced breakdown spectroscopy

Laser-induced breakdown plasma is produced by using Q-switched Nd: YAG laser operating at 532 nm, which interacts with the Al alloy sample target in air. The spectral lines in the 230-440 nm wavelength range have been identified, and based on the calibration-free method, the mass concentration of Al alloy are obtained, which is in good agreement with the standard value of the sample.

Enhancements of the simulation method on the edge effect in resin transfer molding processes

In resin transfer molding processes, small clearances exist between the fiber preform and the mold edges, which result in a preferential resin flow in the edge channel and then disrupt the flow patterns during the mold filling stage. A mathematical model including the effect of cavity thickness on resin flow was developed for flow behavior involving the interface between an edge channel and a porous medium. According to the mathematical analysis of momentum equations in a fully developed rectangular duct and formulations of the equivalent edge permeability, comparing with three-dimensional Navier-Stokes equations, the governing equations were modified in the edge channel. The volume of fluid (VOF) method was applied to track the flow front. A simple case is numerically simulated using the modified governing equations. The effects of edge channel width and cavity thickness on flow front and inlet pressure are analyzed, and the evolution characteristics of simulated results are in agreement with the experimental results. (c) 2007 Elsevier B.V. All rights reserved

Aptamer-based label-free method for hemin recognition and DNA assay by capillary electrophoresis with chemiluminescence detection

An aptamer-based label-free approach to hemin recognition and DNA assay using capillary electrophoresis with chemiluminescence detection is introduced here. Two guanine-rich DNA aptamers were used as the recognition element and target DNA, respectively. In the presence of potassium ions, the two aptamers folded into the G-quartet structures, binding hemin with high specificity and affinity. Based on the G-quartet-hemin interactions, the ligand molecule was specifically recognized with a K (d)approximate to 73 nM, and the target DNA could be detected at 0.1 mu M. In phosphate buffer of pH 11.0, hemin catalyzed the H2O2-mediated oxidation of luminol to generate strong chemiluminescence signal; thus the target molecule itself served as an indicator for the molecule-aptamer interaction, which made the labeling and/or modification of aptamers or target molecules unnecessary. This label-free method for molecular recognition and DNA detection is therefore simple, easy, and effective.

SPECTROPHOTOMETRIC DETERMINATION OF LANTHANUM AND SUM OF OTHER LIGHT RARE-EARTHS BY FLEXIBLE TOLERANCE SIMPLEX-METHOD

The absorptivities of color elements in a mixture can be obtained by using Gauas' elimination with selection of principal element in matrix to the standards. These values can be applied to flexible tolerance simplex method to give the composition of samples. In the exprimental design and data treatment, an effort was made to minimize the errors of results according the principal of optimization. When the difference of absorptivities of color material is significant to the exprimental error, the pr...

塑料压力管道热板焊接接头瞬时应力有限元分析

基于热粘弹性积分型本构关系,考虑材料性能依赖于温度变化及相变潜热的影响,利用有限元软件ANSYS热-力耦合及载荷步功能模拟结晶型高密度聚乙烯塑料压力管道热板焊接过程。并对焊接接头的应力分布进行有限元分析,得到了环向、轴向及径向瞬态应力分布规律。采用盲孔法和锯切法测量焊后残余应力,实测结果与数值分析基本吻合。

高密度聚乙烯塑料压力管道热板焊接头应力分布有限元分析

基于热粘弹性积分型本构关系,考虑材料性能依赖于温度变化及相变潜热的影响,利用AN SYS热-力耦合及载荷步功能模拟结晶型高密度聚乙烯(HDPE)塑料压力管道热板焊接过程,并对焊接接头的应力分布进行有限元分析,得到了环向、轴向以及径向瞬态应力分布的基本规律。采用盲孔法和锯切法测量焊后残余应力,实测结果与数值分析基本符合

基于窗口精细结构描述的岩体质量评价及力学参数综合优选

This thesis bases on horizontal research project “The research about the fine structure and mechanical parameters of abutment jointed rock mass of high arch dam on Jinping Ⅰ Hydropower Station, Yalong River” and “The research about the fine structure and mechanical parameters of the columnar basalt rock mass on Baihetan Hydropower Station, Jinsha River”. A rounded system about the fine structure description and rock mass classification is established. This research mainly contains six aspects as follow: (1) Methods about fine structure description of the window rock mass; (2) The window rock mass classification about the fine structure; (3) Model test study of intermittent joints; (4) Window rock mass strength theory; (5) Numerical experimentations about window rock mass; (6) The multi-source fusion of mechanical parameters based on Bayes principle. Variation of intact rock strength and joint conditions with the weathering and relaxation degree is studied through the description of window rock mass. And four principal parameters: intact rock point load strength, integration degree of window rock mass, joint conditions, and groundwater condition is selected to assess the window rock mass. Window rock mass is classified into three types using the results of window rock mass fine structure description combined with joints develop model. Scores about intact rock strength, integrality condition, divisional plane condition and groundwater conditions are given based on window rock mass fine structure description. Then quality evaluation about two different types of rock mass: general joint structure and columnar jointing structure are carried out to use this window rock mass classification system. Application results show that the window rock mass classification system is effective and applicable. Aimed at structural features of window structure of “the rock mass damaged by recessive fracture”, model tests and numerical models are designed about intermittent joints. By conducting model tests we get shear strength under different normal stress in integrated samples, through samples and intermittent joints samples. Also, the changing trends of shear strength in various connectivity rates are analyzed. We numerically simulate the entire process of direct shear tests by using PFC2D. In order to tally the stress-strain curve of numerical simulation with experimental tests about both integrated samples and through samples, we adjust mechanical factors between particles. Through adopting the same particle geometric parameter, the numerical sample of intermittent joints in different connective condition is re-built. At the same time, we endow the rock bridges and joints in testing samples with the fixed particle contacting parameters, and conduct a series of direct shear tests. Then the destructive process and mechanical parameters in both micro-prospective and macro-prospective are obtained. By synthesizing the results of numerical and sample tests and analyzing the evolutionary changes of stress and strain on intermittent joints plane, we conclude that the centralization of compressive stress on rock bridges increase the shear strength of it. We discuss the destructive mechanics of intermittent joints rock under direct shear condition, meanwhile, divide the whole shear process into five phases, which are elasticity phase, fracture initiation phase, peak value phase, after-peak phase and residual phase. In development of strength theory, the shear strength mechanisms of joint and rock bridge are analyzed respectively. In order to apply the deducted formulation conveniently in the real projects, a relationship between these formulations and Mohr-Coulomb hypothesis is built up. Some sets of numerical simulation methods, i.e. the distinct element method (UDEC) based on in-situ geology mapping are developed and introduced. The working methods about determining mechanical parameters of intact rock and joints in numerical model are studied. The operation process and analysis results are demonstrated detailed from the research on parameters of rock mass based on numerical test in the Jinping Ⅰ Hydropower Station and Baihetan Hydropower Station. By comparison，the advantages and disadvantages are discussed. Results about numerical simulation study show that we can get the shear strength mechanical parameters by changing the load conditions. The multi-source rock mass mechanical parameters can be fused by the Bayes theory, which are test value, empirical value and theoretical value. Then the value range and its confidence probability of different rock mass grade are induced and these data supports the reliability design.

大型反倾岩体变形破坏机制及稳定性分析———以金沙江宽谷坝址龙蟠边坡为例

A large number of catastrophic accidents were aroused by the instability and destruction of anti-dip rock masses in the worldwide engineering projects, such as hydropower station, mine, railways and so on. Problems in relation to deformation and failure about anti-dip rock slopes are significant for engineering geology research. This dissertation takes the Longpan slope in the Jinsha River as a case to study the deformation mechanism of large-scale anti-dip rock masses and the slope stability analysis method. The primary conclusions are as follows. The Dale Reach of Jinsha River, from Longpan to the debouchment of Chongjiang tributary, is located in the southeastern margin of the Qinghai-Tibet Plateau. Longpan slope is the right embankment of Dale dam, it is only 26 km to the Shigu and 18 km to Tiger Leaping Gorge. The areal geology tectonic structures here area are complicated and blurry. Base on the information of geophysical exploration (CSAMT and seismology) and engineering geological investigation, the perdue tectonic pattern of Dale Reach is put forward for the first time in this paper. Due to the reverse slip of Longpan fault and normal left-rotation of Baihanchang fault, the old faulted valley came into being. The thick riverbed sediments have layered characters of different components and corresponding causes, which attribute to the sedimentary environments according with the new tectonic movements such as periodic mountain uplifting in middle Pleistocene. Longpan slope consists of anti-dip alternate sandstone and slate stratums, and the deformable volume is 6.5×107m3 approximately. It was taken for an ancient landslide or toppling failure in the past so that Dale dam became a vexed question. Through the latest field surveying, displacement monitoring and rock masses deforming characters analyses, the geological mechanism is actually a deep-seated gravitational bending deformation. And then the discrete element method is used to simulate the deforming evolution process, the conclusion accords very well with the geo-mechanical patterns analyses. In addition strength reduction method based on DEM is introduced to evaluate the factor of safety of anti-dip rock slope, and in accordance with the expansion way of the shear yielding zones, the progressive shear failure mechanism of large-scale anti-dip rock masses is proposed for the first time. As an embankment or a close reservoir bank to the lower dam, the stability of Longpan slope especially whether or not resulting in sliding with high velocity and activating water waves is a key question for engineering design. In fact it is difficult to decide the unified slip surface of anti-dip rock slope for traditional methods. The author takes the shear yielding zones acquired form the discrete element strength reduction calculation as the potential sliding surface and then evaluates the change of excess pore pressure and factor of stability of the slope generated by rapid drawdown of ponded water. At the same time the dynamic response of the slope under seismic loading is simulated through DEM numerical modeling, the following results are obtained. Firstly the effective effect of seismic inertia force is resulting in accumulation of shear stresses. Secondly the discontinuous structures are crucial to wave transmission. Thirdly the ultimate dynamic response of slope system takes place at the initial period of seismic loading. Lastly but essentially the effect of earthquake load to bringing on deformation and failure of rock slope is the coupling effect of shear stresses and excess pore water pressure accumulation. In view of limitations in searching the critical slip surface of rock slope of the existing domestic and international software for limit equilibrium slope stability analyses, this article proposes a new method named GA-Sarma Algorithm for rock slope stability analyses. Just as its name implies, GA-Sarma Algorithm bases on Genetic Algorithm and Sarma method. GA-Sarma Algorithm assumes the morphology of slip surface to be a broken line with traceability to extend along the discontinuous surface structures, and the slice boundaries is consistent with rock mass discontinuities such as rock layers, faults, cracks, and so on. GA-Sarma Algorithm is revolutionary method that is suitable for global optimization of the critical slip surface for rock slopes. The topics and contents including in this dissertation are closely related to the difficulties in practice, the main conclusions have been authorized by the engineering design institute. The research work is very meaningful and useful for the engineering construction of Longpan hydropower station.

龙游石窟洞室群稳定性分析及最佳旅游路线选择

The large ancient underground rock caverns in Longyou is an important component of grotto cultural. Current task facing the long-term preservation of these unmovable cultural relics is arduous and challenging. The deformation failure of the caverns' surrounding rock is deteriorating. The weathering velocity of these caverns is accelerating. With the strength of caverns' surrounding rock worsening, critical rocks were generated in local regions of the caverns' vault and posing a threat to the security of people passing by. Selection of a maximum-security route and construction a aisle in the caverns might be an efficient way to ensure the security of tourists and reach the target of long-term preservation. The deformation and destruction of the ancient underground caverns is primarily dominated by geological conditions and the special structure of caverns. Based on field investigation, several fundamental conditions for deformation and failure are recognized, and nine deformation and fracture patterns of the Longyou grotto are proposed. In order to judge the stability of caverns’ surrounding rock, the element safety coefficient method is presented. An explicit explanation for the meaning of the method is deduced using Mohr-Coulomb strength criterion. Numerical analyses are carried out in the dissertation through FLAC3D code. Through numerical analysis, the stress distribution regularities of the caverns’ roofs, piles and public side wall are analysed, and the stability properties of caverns’ surrounding rock are also assessed. At the same time, the element safety coefficient method is introduced to contrast the stability degree of different regions in caverns. The above analyses are bases for choosing the optimal tourism routes in the caverns of Longyou grotto. The impact of surface load on the stability of shallow buried cavities in Longyou grotto is evaluated, the results show that building load has significant influence on the stability of the No.1 cavern’s roof, pile and public side wall between the No.1 cavern and the No.2 cavern, pedestrian load has less impact on the stability of surrounding rock than building load. The principles for choosing the optimal tourism routes in the caverns are discussed. With these principles, the dissertation makes a systematic research on the geological analytic method, numerical analytic method and meeting tourism requirements method, which are used in selecting the optimal tourism routes in the caverns. In order to achieve the best effect in the process of tourism routes selection, the above three method are integrated through Theory of Engineering Geomechanics Meta-system(EGMS). According to field investigations, numerical analyses, tourism requirements and expert experiences, the optimal tourism routes through No.1 to No.5 cavern are determined preliminarily. The obtained results from the research work are useful for the security aisle's construction, they also have reference value to other projects in practice.

岩石类材料静态、动态拉伸弹性模量测试方法研究——以巴西圆盘试验和霍普金森压杆、效应为手段

As we all know, rock-like materials will absolutely show very different mechanical properties under the compressive stress and tensile stress respectively. Similarly, under the dynamic compressive stress or dynamic tensile stress, the characteristics of the dynamics showed by the rock-like materials also have great differences from the mechanical behavior under static force. Studying their similarities and differences in rock mechanics theory and practical engineering will be of great significance. Generally, there are compression modulus of elasticity and tensile modulus of elasticity corresponding to compressive stress state and the tensile stress state in the rock. Both the two kinds of elastic modulus play an extremely important role in calculation of engineering mechanics. Their reliability directly affects the accuracy and reliability of the calculation results of internal stress field and displacement field of engineering rock mass. At present, it is easy to obtain the compression modulus of elasticity in laboratory; but it is very difficult to determine the tensile modulus of elasticity with direct tensile test due to that direct tensile test is difficult to perform in laboratory in general. In order to solve this problem, this thesis invents and develops several indirect test methods to determine the static or dynamic tensile modulus of elasticity of rock-type materials with high reliability and good interoperability. For the static tensile modulus of elasticity, the analytical stress field solution has been given out for the Brazilian disc under the radial and linear concentration load with Airy stress function method. At the same time, the stress field has been modeled for the Brazilian disc test by using the finite element software of ANSYS and ADINA. The analytical stress field solution is verified to be right by comparatively researching the analytical stress field solution and the numerical stress field solution. Based on the analytical stress field solution, this thesis proposes that a strain gauge is pasted at the Brazilian disc center along the direction perpendicular to the applied force to indirectly determine the static tensile modulus of elasticity, and related measurement theory also has been developed. The method proposed here has good feasibility and high accuracy verified by the experimental results. For the dynamic tensile modulus of elasticity, two measuring methods and theories are invented here. The first one is that the Split Hopkinson Pressure Bar is used to attract the Brazilian disc to generate the dynamic load, make the dynamic tensile stress is formed at the Brazilian disc center; and also a strain gauge is pasted at the Brazilian disc center to record the deformation. The second is that, in the Hopkinson effect phenomenon, the reflection tensile stress wave is formed when the shock wave propagates to the free end of cylindrical rock bar and reflect, which can make the rock bar is under dynamic tensile stress state; and some strain gauges are pasted at the appropriate place on the rock bar to record the strain coursed by the tensile or compressive stress wave. At last, the dynamic tensile modulus of elasticity can be determined by the recorded strain and the dynamic tensile stress which can be determined by related theories developed in this thesis.