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.

Slope Stability Analysis Using Anisotropic and Nonlinear Failure Criteria

The shear strength of soils or rocks developed in a landslide usually exhibits anisotropic and nonlinear behavior. The process of sedimentation and subsequent consolidation can cause anisotropy of sedimentary soils or rocks, for instance. Nonlinearity of failure envelope could be attributed to "interlocking" or "dilatancy" of the material, which is generally dependent upon the stress level. An analytical method considering both anisotropy and nonlinearity of the failure envelops of soil and rocks is presented in the paper. The nonlinearfailure envelopes can be determined from routine triaxial tests. A spreadsheet program, which uses the Janbu's Generalized Procedure of Slice and incorporates anisotropic, illustrates the implementation of the approach and nonlinearfailure envelops. In the analysis, an equivalent Mohr-Coulomb linear failure criterion is obtained by drawing a tangent to the nonlinear envelope of an anisotropic soil at an appropriate stress level. An illustrative example is presented to show the feasibility and numerical efficiency of the method.

COMPARISON OF THE CRACK-GROWTH METHOD AND THE CRACK STRESS WHITENING ZONE METHOD FOR THE FRACTURE-TOUGHNESS DETERMINATION OF PHENOLPHTHALEIN POLY(ETHER KETONE)

Fracture toughness values of phenolphthalein poly(ether ketone) (PEK-C) at 190 degrees C were determined by two different methods, i. e. the conventional crack growth method and the crack stress whitening zone method, which show consistent results. This indicates that the crack stress whitening zone method can be used to determine the crack initiation of some polymers for which the blunting line concept is unsuitable.

Numerical Simulation and Experimental Investigation about Internal and External Flows

In this paper, TASCflow3D is used to solve inner and outer 3D viscous incompressible turbulent flow (R-e = 5.6 X 10(6)) around axisymmetric body with duct. The governing equation is a RANS equation with standard k-epsilon turbulence model. The discrete method used is a finite volume method based on the finite element approach. In this method, the description of geometry is very flexible and at the same time important conservative properties are retained. The multi-block and algebraic multi-grid techniques are used for the convergence acceleration. Agreement between experimental results and calculation is good. It indicates that this novel approach can be used to simulate complex flow such as the interaction between rotor and stator or propulsion systems containing tip clearance and cavitation.

Numerical Simulation and Experimental Investigation about Internal and External Flows

In this paper, TASCflow3D is used to solve inner and outer 3D viscous incompressible turbulent flow (R-e = 5.6 X 10(6)) around axisymmetric body with duct. The governing equation is a RANS equation with standard k-epsilon turbulence model. The discrete method used is a finite volume method based on the finite element approach. In this method, the description of geometry is very flexible and at the same time important conservative properties are retained. The multi-block and algebraic multi-grid techniques are used for the convergence acceleration. Agreement between experimental results and calculation is good. It indicates that this novel approach can be used to simulate complex flow such as the interaction between rotor and stator or propulsion systems containing tip clearance and cavitation.

Multi-pipe string electromagnetic detection tool and its applications

The MID-K, a new kind of multi-pipe string detection tool is introduced. This tool provides a means of evaluating the condition of in-place pipe string, such as tubing and casino. It is capable of discriminating the defects of the inside and outside, and estimating the thickness of tubing and casing. It is accomplished by means of a low frequency eddy current to detect flaws on the inner surface and a magnetic flux leakage to inspect the full thickness. The measurement principle, the technology and applications are presented in this paper.

A new method of measuring the waveguide dispersion factor and the thermo-optic coefficient of long-period fiber gratings

Based on the Mach-Zehnder effect between the core mode and the cladding modes, the interference fringes are formed by a pair of cascaded long-period fiber gratings (CLPFGs). Theoretical analyses show that the spectral spacing and the wavelength of these fringes are functions of the waveguide dispersion factor gamma, which is a characterizing parameter to LPFG and with theoretical and applicational significance. By measuring the characteristics of the transmission spectra of CLPFGs, the absolute value of gamma can be obtained. At the same time, the thermo-optic coefficient of effective refractive index difference between core and cladding modes, p, can also be obtained by measured the temperature sensitivity of these fringes. In the experiments, \gamma\ and mu were measured by this method to be 0.874 and 4.08 x 10(-5) degreesC(-1), respectively, for LPFGs with period of 450 mum and with a HE14 resonant peak at 1554 nm. (C) 2004 Elsevier B.V. All rights reserved.

Origin of the double- and multi-pulse structure of echolocation signals in Yangtze finless porpoise (Neophocaena phocaenoides asiaeorientialis)

The signals of dolphins and porpoises often exhibit a multi-pulse structure. Here, echolocation signal recordings were made from four geometrically distinct positions of seven Yangtze finless porpoises temporarily housed in a relatively small, enclosed area. Some clicks demonstrated double-pulse, and others multi-pulse, structure. The interpulse intervals between the first and second pulse of the double- and multi-pulse clicks were significantly different among data from the four different positions (p < 0.01, one-way ANOVA). These results indicate that the interpulse interval and structure of the double- and multi-pulse echolocation signals depend on the hydrophone geometry of the animal, and that the double- and multi-pulse structure of echolocation signals in Yangtze finless porpoise is not caused by the phonating porpoise itself, but by the multipath propagation of the signal. Time delays in the 180 degrees phase-shifted surface reflection pulse and the nonphase-shifted bottom reflection pulse of the multi-pulse structures, relative to the direct signal, can be used to calculate the distance to a phonating animal. (c) 2005 Acoustical Society of America.

Study on the forecast effects of PI method to the north and southwest China

In this paper, the codes of Pattern Informatics (PI) method put forward by Rundle et al. have been worked out according to their algorithm published, and the retrospective forecast of PI method to North China (28.0 degrees-42.0 degrees N, 108.0 degrees-125.0 degrees E) and to Southwest China (22.0 degrees-28.3 degrees N, 98.0 degrees-106.0 degrees E) has been tested. The results show that the hit rates in different regions show a great difference. In Southwest China, 32 earthquakes with M(L)5.0 or larger have occurred during the predicted time period 2000-2007, and 26 out of the 32 earthquakes occurred in or near the hot spots. In North China, the total number of M(L)5.0 or larger was 12 during the predicted time period 2000-2007, and only 3 out of the 12 earthquakes occurred in or near the hot spots. From our results, we hold that if the PI method could be applied to all kinds of regions, the parameters associated with time points and time windows should be chosen carefully to obtain the higher hit rate. We also found that the aftershocks in a strong earthquake sequence affect the PI results obviously. Copyright (c) 2009 John Wiley & Sons, Ltd.

Study on the forecast effects of PI method to the north and southwest China

In this paper, the codes of Pattern Informatics (PI) method put forward by Rundle et al. have been worked out according to their algorithm published, and the retrospective forecast of PI method to North China (28.0 degrees-42.0 degrees N, 108.0 degrees-125.0 degrees E) and to Southwest China (22.0 degrees-28.3 degrees N, 98.0 degrees-106.0 degrees E) has been tested. The results show that the hit rates in different regions show a great difference. In Southwest China, 32 earthquakes with M(L)5.0 or larger have occurred during the predicted time period 2000-2007, and 26 out of the 32 earthquakes occurred in or near the hot spots. In North China, the total number of M(L)5.0 or larger was 12 during the predicted time period 2000-2007, and only 3 out of the 12 earthquakes occurred in or near the hot spots. From our results, we hold that if the PI method could be applied to all kinds of regions, the parameters associated with time points and time windows should be chosen carefully to obtain the higher hit rate. We also found that the aftershocks in a strong earthquake sequence affect the PI results obviously. Copyright (c) 2009 John Wiley & Sons, Ltd.

Photoionization of 1s electron and corresponding shake-up process in ground and excited lithium atoms

The cross sections of the 18 electron photoionization and corresponding shake-up processes for Li atoms in the ground state 1s(2)2s and excited states 1s(2)2p, 1s(2)3p, 1s(2)3p and 1s(2)3d are calculated using the multi-configuration Dirac-Fock method. The latest experimental photoelectron spectrum at hv = 100 eV [Cubaynes D et al. Phys. Rev. Lett. 99 (2007) 213004] has been reproduced by the present theoretical investigation excellently. The relative intensity of the shake-up satellites shows that the effects of correlation and relaxation become more important for the higher excited states of the lithium atom, which are explained very well by the spatial overlap of the initial and final state wavefunctions. In addition, strong dependence of the cross section on the atomic orbitals of the valence electrons are found, especially near the threshold.

Studies on the method of leaps and bounds regression

In this paper, we introduce the method of leaps and bounds regression which can be used to select variables quickly and obtain the best regression models. These models contain one variable, two variables, three variables and so on. The results obtained by using leaps and bounds regression were compared with those achieved by using stepwise regression to lead to the conclusion that leaps and bounds regression is an effective method.

Simultaneous enumeration of diatom, protozoa and meiobenthos from marine sediments using Ludox-QPS method

The Ludox-QPS method is a newly developed technique, which combines the Ludox HS 40 density centrifugation and quantitative protargol stain, to enumerate marine ciliates with good taxonomic resolution. We tested the method for simultaneous enumeration of diatoms, protozoa and meiobenthos and compared its extraction efficiency for meiobenthos with that of the routine Ludox-TM centrifugation and a modified protocol using Ludox HS 40. We conducted the evaluation with a sample size of 8.3 ml each from sandy, muddy-sand and muddy sediments collected from the intertidal area of the Yellow Sea in summer 2006 and spring 2007. The Ludox-QPS method not only produced high extraction efficiencies of 97 +/- 1.3% for diatoms and 97.6 +/- 0.8% for ciliates, indicating a reliable enumeration for eukaryotic microbenthos, but also produced excellent extraction efficiencies of on average 97.3% for total meiobenthos, 97.9% for nematodes and 97.8% for copepods from sands, muddy sands and mud. By contrast, the routine Ludox-TM centrifugation obtained only about 74% of total meiobenthos abundance with one extraction cycle, and the modified Ludox HS 40 centrifugation yielded on average 93% of total meiobenthos: 89.4 +/- 2.0% from sands, 93 +/- 4.1% from muddy sands and 97.1 +/- 3.0% from mud. Apart from the sediment type, sample volume was another important factor affecting the extraction efficiency for meiobenthos. The extraction rate was increased to about 96.4% when using the same modified Ludox centrifugation for a 4 ml sediment sample. Besides the excellent extraction efficiency, the Ludox-QPS method obtained higher abundances of meiobenthos, in particular nematodes, than the routine Ludox centrifugation, which frequently resulted in an uncertain loss of small meiobenthos during the sieving process. Statistical analyses demonstrated that there were no significant differences between the meiobenthos communities revealed by the Ludox-QPS method and the modified Ludox HS 40 centrifugation, showing the high efficiency of the Ludox-QPS method for simultaneous enumeration of diatom, protozoa and meiobenthos. Moreover, the comparatively high taxonomic resolution of the method, especially for diatoms and ciliates, makes it feasible to investigate microbial ecology at community level.