Shallow Water Model On Cubed-Sphere By Multi-Moment Finite Volume Method

A global numerical model for shallow water flows on the cubed-sphere grid is proposed in this paper. The model is constructed by using the constrained interpolation profile/multi-moment finite volume method (CIP/MM FVM). Two kinds of moments, i.e. the point value (PV) and the volume-integrated average (VIA) are defined and independently updated in the present model by different numerical formulations. The Lax-Friedrichs upwind splitting is used to update the PV moment in terms of a derivative Riemann problem, and a finite volume formulation derived by integrating the governing equations over each mesh element is used to predict the VIA moment. The cubed-sphere grid is applied to get around the polar singularity and to obtain uniform grid spacing for a spherical geometry. Highly localized reconstruction in CIP/MM FVM is well suited for the cubed-sphere grid, especially in dealing with the discontinuity in the coordinates between different patches. The mass conservation is completely achieved over the whole globe. The numerical model has been verified by Williamson's standard test set for shallow water equation model on sphere. The results reveal that the present model is competitive to most existing ones. (C) 2008 Elsevier Inc. All rights reserved.

A Cip/Multi-Moment Finite Volume Method For Shallow Water Equations With Source Terms

A novel finite volume method has been presented to solve the shallow water equations. In addition to the volume-integrated average (VIA) for each mesh cell, the surface-integrated average (SIA) is also treated as the model variable and is independently predicted. The numerical reconstruction is conducted based on both the VIA and the SIA. Different approaches are used to update VIA and SIA separately. The SIA is updated by a semi-Lagrangian scheme in terms of the Riemann invariants of the shallow water equations, while the VIA is computed by a flux-based finite volume formulation and is thus exactly conserved. Numerical oscillation can be effectively avoided through the use of a non-oscillatory interpolation function. The numerical formulations for both SIA and VIA moments maintain exactly the balance between the fluxes and the source terms. 1D and 2D numerical formulations are validated with numerical experiments. Copyright (c) 2007 John Wiley & Sons, Ltd.

Estudio del desempeño de OLSR en una red mallada inalámbrica en escenario real

Resumen: El objetivo de este trabajo es analizar el comportamiento del protocolo OLSR sobre una red mallada configurada con firmware OpenWrt utilizando distintos equipos de hardware. Se presentan los resultados empíricos de varias pruebas utilizando el mismo escenario. El escenario que se presenta es una red de mundo real (no de laboratorio) con pruebas reales y no simulaciones. OpenWrt es un software perfectamente válido que puede ser utilizado en una gran variedad de dispositivos y su configuración para utilizarlo con protocolo OLSR es sencilla de realizar y no presenta problemas de funcionamiento con dicho protocolo.

Experimental study of gaseous detonation propagation over acoustically absorbing walls

The paper presents results front an experimental investigation of the propagation of gaseous detonation waves over tube sections lined with acoustically absorbent materials. The measurements were compared with results from control tests in a smooth wall section. The results show the increasing effectiveness of a perforated steel plate, wire mesh and steel wool in attenuating detonation.

Novel experimental phenomena of fine-particle fluidized beds

The dynamic response of bed height and concentration waves in liquid-solid fluidized beds to a step change in the fluidization velocity is considered. We experimentally study the liquid-solid fluidized beds, spherical beadings, with sizes ranging from 230 to 270 mesh and the inner diameter of columns made from glass is 2.4 mm. Experimental results find that under certain conditions, fine particles with large Richardson-Zaki exponent n display different dynamic behavior from usual particles with smaller n during expansion and collapse of the fluidized state. (c) 2007 Elsevier Inc. All rights reserved.

A new method for evaluation of stress intensities for interface cracks

A new method is presented for calculating the values of K-I and K-II in the elasticity solution at the tip of an interface crack. The method is based on an evaluation of the J-integral by the virtual crack extension method. Expressions for calculating K-I and K-II by using the displacements and the stiffness derivative of the finite element solution and asymptotic crack tip displacements are derived. The method is shown to produce very accurate solutions even with coarse element mesh.

Collapsed isoparametric element as a singular element for a crack normal to the bi-material interface

It is shown that the variable power singularity of the strain field at the crack tip can be obtained by the simple technique of collapsing quadrilateral isoparametric elements into triangular elements around the crack tip and adequately shifting the side-nodes adjacent to this crack tip. The collapsed isoparametric elements have the desired singularity at crack tip along any ray. The strain expressions for a single element have been derived and in addition to the desired power singularity, additional singularities are revealed. Numerical examples have shown that triangular elements formed by collapsing one side lead to excellent results.

Local stresses and strains in axially cracked cylindrical shells

A perturbation solution is obtained for the local stress-strain fields in an axially cracked cylindrical shell. The tenth-order differential equations are used that take into account the transverse shear deformation. The perturbation of a curvature parameter, λ, is employed, where . The stress intensity factors for finite size cylindrical shells subjected to bending and internal pressure are evaluated. Sufficient accuracy can be obtained without using fine mesh sizes in regions near the crack tip. Also analyzed are the influence of cylinder diameter and shearing stiffness on bulging.

Prediction of mixed-mode cracking direction in random, short-fibre composite-materials

The prediction of cracking direction in composite materials is of significance to the design of composite structures. This paper presents several methods for predicting the cracking direction in the double grooved tension-shear specimen which gives mixed-mode cracking. Five different criteria are used in this analysis: two of them have been used by other investigators and the others are proposed by the present authors. The strain energy density criterion proposed by G.C. Sih is modified to take account of the influence of the anisotropy of the strength on the direction of crack. The two failure criteria of Tsai-Hill and Norris are extended to predict the crack orientation. The stress distributions in the near-notch zone are calculated by using the 8-node quadrilateral isoparametric finite element method. The predictions of all the criteria except one are in good agreement with the experimental measurement. In addition, on the basis of the FEM results, the size of the zone in which the singular term is dominant is estimated.

Dynamic-response of a laminated plate with friction damping

快速Delaunay逐点插入网格生成算法

对插入形心的Delaunay逐点插入算法,提出按单元可插度分组的双向链表组数据结构,避免了对最大可插度单元的搜索。采用了邻接单元搜索、双向链表存储、随机方向搜索、邻接旋转、几何量继承等技术,使算法的计算时间与生成单元数近似呈线性关系,时间复杂度达到O(N1.05),N为生成单元数。算例表明,在一台AMD Athlon3200+(主频2.0 GHz)PC上,该算法的四面体单元生成速度达50 000个/s以上。

Water quality impacts of mechanical shredding of aquatic macrophytes

We examined the impacts of mechanical shredding (i.e.. shredding plants and leaving biomass in the system) of the water chestnut (Trapa natans) on water quality and nutrient mobilization in a control and experimental site in Lake Champlain (Vermont-New York). A 1-ha plot was mechanically shredded within 1 h on 26 July, 1999. Broken plant material was initially concentrated on the lake surface of the experimental station after shredding, and was noticeable on the lake surface for 19 d. Over a two week period after shredding. concentrations of total nitrogen (N) and phosphorus (P), and soluble reactive P increased in the lower water column of the experimental station, coinciding with decomposition of water chestnut. Sediments in the control and experimental stations exhibited vet-v low rates of N and P release and could not account for increases in nutrient concentrations in the water column after mechanical shredding. Shredded plant material deployed in mesh bags at the experimental station lost similar to 70% of their total mass, and 42%, N and 70% P within 14 d. indicating Substantial nutrient mobilization via autolysis and decomposition. Chlorophyll a concentrations increased to 35 g/L at the experimental station on day 7 after shredding, compared to a concentration of 4 g/L at the control station. suggesting uptake of mobilized nutrients by phytoplankton. Disruption Of the Surface canopy of water chestnut by shredding was associated with marked increases in turbidity and dissolved oxygen, suggesting increased mixing at the experimental site.

Influence of Dilute Acetic Acid Treatments on Survival of Monoecious Hydrilla Tubers in the Oregon House Canal, California

Hydrilla (Hydrilla verticillata (L.f.)Royle), a serious aquatic weed, reproduces through formation of underground tubers. To date, attacking this life-cycle stage has been problematic. The purpose of this study was to measure the impact of exposure to dilute acetic acid on monoecious hydrilla tubers under field conditions. In this field experiment, treatments were acetic acid concentration (0, 2.5, or 5%) and sediment condition (perforated or not perforated). Each of 60, 1x1 m plots (in the Oregon House Canal) were randomly assigned to one treatment. Two weeks after treatment, we collected three samples from each plot. One was washed over 2 mm wire mesh screens to separate tubers from sediment. Relative electrolyte leakage was measured for one tuber from each plot. Five additional tubers from each plot were placed in a growth chamber and sprouting monitored for four weeks. A second sample from each plot was placed in a plastic tub and placed in an outdoor tank, filled with water. These samples were monitored for tuber sprouting. Relative electrolyte leakage increased significantly for tubers exposed to 2.5% or 5% acetic acid. Effects on tubers in perforated sediment were reduced. Exposure to acetic acid inhibited tuber sprouting by 80 to 100%, in both chamber and outdoor tests. These results confirm findings from earlier laboratory/greenhouse experiments, and suggest that this approach may be useful in the management of hydrilla tuber banks in habitats where the water level can be lowered to expose the sediments.

Numerical Simulation of R-M Instability

In order to capture shock waves and contact discontinuities in the field and easy to program with parallel computation a new algorithm is developed to solve the N-S equations for simulation of R-M instability problems. The method with group velocity control is used to suppress numerical oscillations, and an adaptive non-uniform mesh is used to get fine resolution. Numerical results for cylindrical shock-cylindrical interface interaction with a shock Mach number Ms=1.2 and Atwood number A=0.818, 0.961, 0.980 (the interior density of the interface/outer density p(1)/p(2) = 10, 50, 100, respectively), and for the planar shock-spherical interface interaction with Ms=1.2 and p(1)/p(2) = 14.28are presented. The effect of Atwood number and multi-mode initial perturbation on the R-M instability are studied. Multi-collisions of the reflected shock with the interface is a main reason of nonlinear development of the interface instability and formation of the spike-bubble structures In simulation with double mode perturbation vortex merging and second instability are found. After second instability the small vortex structures near the interface produced. It is important factor for turbulent mixing.

Catastrophic Rupture Induced Damage Coalescence in Heterogeneous Brittle Media

In heterogeneous brittle media, the evolution of damage is strongly influenced by the multiscale coupling effect. To better understand this effect, we perform a detailed investigation of the damage evolution, with particular attention focused on the catastrophe transition. We use an adaptive multiscale finite-element model (MFEM) to simulate the damage evolution and the catastrophic failure of heterogeneous brittle media. Both plane stress and plane strain cases are investigated for a heterogeneous medium whose initial shear strength follows the Weibull distribution. Damage is induced through the application of the Coulomb failure criterion to each element, and the element mesh is refined where the failure criterion is met. We found that as damage accumulates, there is a stronger and stronger nonlinear increase in stress and the stress redistribution distance. The coupling of the dynamic stress redistribution and the heterogeneity at different scales result in an inverse cascade of damage cluster size, which represents rapid coalescence of damage at the catastrophe transition.