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.

Complex potentials and singular integral equation for curve crack problem in antiplane elasticity

Four types of the fundamental complex potential in antiplane elasticity are introduced: (a) a point dislocation, (b) a concentrated force, (c) a dislocation doublet and (d) a concentrated force doublet. It is proven that if the axis of the concentrated force doublet is perpendicular to the direction of the dislocation doublet, the relevant complex potentials are equivalent. Using the obtained complex potentials, a singular integral equation for the curve crack problem is introduced. Some particular features of the obtained singular integral equation are discussed, and numerical solutions and examples are given.

Application of enthalpy method for moving boundary problem in laser surface melting

A numerical analysis was carried out to study the moving boundary problem in the physical process of pulsed Nd-YAG laser surface melting prior to vaporization. The enthalpy method was applied to solve this two-phase axisymmetrical melting problem Computational results of temperature fields were obtained, which provide useful information to practical laser treatment processing. The validity of enthalpy method in solving such problems is presented.

Analysis of strip electric saturation model of crack problem in piezoelectric materials

This paper presents a fully anisotropic analysis of strip electric saturation model proposed by Gao et al. (1997) (Gao, H.J., Zhang, T.Y., Tong, P., 1997. Local and global energy release rates for an electrically yielded crack in a piezoelectric ceramic. J. Mech. Phys. Solids, 45, 491-510) for piezoelectric materials. The relationship between the size of the strip saturation zone ahead of a crack tip and the applied electric displacement field is established. It is revealed that the critical fracture stresses for a crack perpendicular to the poling axis is linearly decreased with the increase of the positive applied electric field and increases linearly with the increase of the negative applied electric field. For a crack parallel to the poring axis, the failure stress is not effected by the parallel applied electric field. In order to analyse the existed experimental results, the stress fields ahead of the tip of an elliptic notch in an infinite piezoelectric solid are calculated. The critical maximum stress criterion is adopted for determining the fracture stresses under different remote electric displacement fields. The present analysis indicates that the crack initiation and propagation from the tip of a sharp elliptic notch could be aided or impeded by an electric displacement field depending on the field direction. The fracture stress predicted by the present analysis is consistent with the experimental data given by Park and Sun (1995) (Park, S., Sun, C.T., 1995. Fracture criteria for piezoelectric materials. J. Am. Ceram. Soc 78, 1475-1480).

Load-Unload Response Ratio and Accelerating Moment/Energy Release Critical Region Scaling and Earthquake Prediction

The main idea of the Load-Unload Response Ratio (LURR) is that when a system is stable, its response to loading corresponds to its response to unloading, whereas when the system is approaching an unstable state, the response to loading and unloading becomes quite different. High LURR values and observations of Accelerating Moment/Energy Release (AMR/AER) prior to large earthquakes have led different research groups to suggest intermediate-term earthquake prediction is possible and imply that the LURR and AMR/AER observations may have a similar physical origin. To study this possibility, we conducted a retrospective examination of several Australian and Chinese earthquakes with magnitudes ranging from 5.0 to 7.9, including Australia's deadly Newcastle earthquake and the devastating Tangshan earthquake. Both LURR values and best-fit power-law time-to-failure functions were computed using data within a range of distances from the epicenter. Like the best-fit power-law fits in AMR/AER, the LURR value was optimal using data within a certain epicentral distance implying a critical region for LURR. Furthermore, LURR critical region size scales with mainshock magnitude and is similar to the AMR/AER critical region size. These results suggest a common physical origin for both the AMR/AER and LURR observations. Further research may provide clues that yield an understanding of this mechanism and help lead to a solid foundation for intermediate-term earthquake prediction.

Dynamic function of damage and its implications

Dynamic function of damage is the key to the problem of damage evolution of solids. In order to understand it, one must understand its mesoscopic mechanisms and macroscopic formulation. In terms of evolution equation of microdamage and damage moment, a dynamic function of damage is strictly defined. The mesoscopic mechanism underlying self-closed damage evolution law is investigated by means of double damage moments. Numerical results of damage evolution demonstrate some common features for various microdamage dynamics. Then, the dynamic function of damage is applied to inhomogeneous damage field. In this case, damage evolution rate is no longer equal to the dynamic function of damage. It is found that the criterion for damage localization is closely related to compound damage. Finally, an inversion of damage evolution to the dynamic function of damage is proposed.

Moment Tensor Analysis of the Acoustic Emission Source in the Rock Damage Process

To further investigate the mechanism of acoustic emission (AE) in the rock fracture experiment, moment tensor analysis was carried out. The AE sources characterized by crack sizes, orientations and fracture modes, are represented by a time-dependent momen

A Lagrangian for von Karman Equations of Large Deflection Problem of Thin Circular Plate

By the semi-inverse method proposed by He, a Lagrangian is established for the large deflection problem of thin circular plate. Ritz method is used to obtain an approximate analytical solution of the problem. First order approximate solution is obtained, which is similar to those in open literature. By Mathematica a more accurate solution can be deduced.

Finite Boundary Effects in Problem of a Crack Perpendicular to and Terminating at a Bimaterial Interface

In this paper, the problem of a crack perpendicular to and terminating at an interface in bimaterial structure with finite boundaries is investigated. The dislocation simulation method and boundary collocation approach are used to derive and solve the basic equations. Two kinds of loading form are considered when the crack lies in a softer or a stiffer material, one is an ideal loading and the other one fits to the practical experiment loading. Complete solutions of the stress field including the T stress are obtained as well as the stress intensity factors. Influences of T stress on the stress field ahead of the crack tip are studied. Finite boundary effects on the stress intensity factors are emphasized. Comparisons with the problem presented by Chen et al. (Int. J. Solids and Structure, 2003, 40, 2731-2755) are discussed also.

On a Degenerate Parabolic Problem in Holder Spaces

In this paper, we study some degenerate parabolic equation with Cauchy-Dirichlet boundary conditions. This problem is considered in little Holder spaces. The optimal regularity of the solution v is obtained and is specified in terms of those of the second member when some conditions upon the Holder exponent with respect to the degeneracy are satisfied. The proofs mainly use the sum theory of linear operators with or without density of domains and the results of smoothness obtained in the study of some abstract linear differential equations of elliptic type.

Rotation Numbers Of Invariant Manifolds Around Unstable Periodic Orbits For The Diamagnetic Kepler Problem

In this paper, a method to construct topological template in terms of symbolic dynamics for the diamagnetic Kepler problem is proposed. To confirm the topological template, rotation numbers of invariant manifolds around unstable periodic orbits in a phase space are taken as an object of comparison. The rotation numbers are determined from the definition and connected with symbolic sequences encoding the periodic orbits in a reduced Poincare section. Only symbolic codes with inverse ordering in the forward mapping can contribute to the rotation of invariant manifolds around the periodic orbits. By using symbolic ordering, the reduced Poincare section is constricted along stable manifolds and a topological template, which preserves the ordering of forward sequences and can be used to extract the rotation numbers, is established. The rotation numbers computed from the topological template are the same as those computed from their original definition.

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.

A contact crack problem in an infinite plate

The problem of an infinite plate with crack of length 2a loaded by the remote tensile stress P and a pair of concentrated forces Q is discussed. The value of the force Q for the initial contact of crack face is investigated and the contact length elevated, while the Q force increases. The problem is solved assuming that the stress intensity factor vanishes at the end point of the contact portion. By the Fredholm integral equation for the multiple cracks, the reduction of stress intensity factor due to Q is found. (C) 1999 Elsevier Science Ltd. All rights reserved.