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.

The plane problem of collinear rigid lines under arbitrary loads

The elastic plane problem of collinear rigid lines under arbitrary loads is dealt with. Applying the Riemann-Schwarz symmetry principle integrated with the analysis of the singularity of complex stress functions, the general formulation is presented, and the closed-form solutions to several problems of practical importance are given, which include some published results as the special cases. Lastly the stress distribution in the immediate vicinity of the rigid line end is examined.

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).

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 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.

A method to find unstable periodic orbits for the diamagnetic Kepler problem

A method to determine the admissibility of symbolic sequences and to find the unstable periodic orbits corresponding to allowed symbolic sequences for the diamagnetic Kepler problem is proposed by using the ordering of stable and unstable manifolds. By investigating the unstable periodic orbits up to length 6, a one to one correspondence between the unstable periodic orbits and their corresponding symbolic sequences is shown under the system symmetry decomposition.

A rigorous analytical method for doubly periodic cylindrical inclusions under longitudinal shear and its application

An infinite elastic solid containing a doubly periodic parallelogrammic array of cylindrical inclusions under longitudinal shear is studied. A rigorous and effective analytical method for exact solution is developed by using Eshelby's equivalent inclusion concept integrated with the new results from the doubly quasi-periodic Riemann boundary value problems. Numerical results show the dependence of the stress concentrations in such heterogeneous materials on the periodic microstructure parameters. The overall longitudinal shear modulus of composites with periodic distributed fibers is also studied. Several problems of practical importance, such as those of doubly periodic holes or rigid inclusions, singly periodic inclusions and single inclusion, are solved or resolved as special cases. The present method can provide benchmark results for other numerical and approximate methods. (C) 2003 Elsevier Ltd. All rights reserved.

Stress intensity factors calculation in anti-plane fracture problem by orthogonal integral extraction method based on femol

For an anti-plane problem, the differential operator is self-adjoint and the corresponding eigenfunctions belong to the Hilbert space. The orthogonal property between eigenfunctions (or between the derivatives of eigenfunctions) of anti-plane problem is exploited. We developed for the first time two sets of radius-independent orthogonal integrals for extraction of stress intensity factors (SIFs), so any order SIF can be extracted based on a certain known solution of displacement (an analytic result or a numerical result). Many numerical examples based on the finite element method of lines (FEMOL) show that the present method is very powerful and efficient.