Computation of stress intensity factors by the sub-region mixed finite element method of lines

Based on the sub-region generalized variational principle, a sub-region mixed version of the newly-developed semi-analytical 'finite element method of lines' (FEMOL) is proposed in this paper for accurate and efficient computation of stress intensity factors (SIFs) of two-dimensional notches/cracks. The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used, with the sought SIFs being among the unknown coefficients. The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements. A mixed system of ordinary differential equations (ODEs) and algebraic equations is derived via the sub-region generalized variational principle. A singularity removal technique that eliminates the stress parameters from the mixed equation system eventually yields a standard FEMOL ODE system, the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver. A number of numerical examples, including bi-material notches/cracks in anti-plane and plane elasticity, are given to show the generally excellent performance of the proposed method.

Magnetohydrodynamics equilibrium of a self-confined elliptic plasma ball

A variational principle is applied to the problem of magnetohydrodynamics (MHD) equilibrium of a self-contained elliptical plasma ball, such as elliptical ball lightning. The principle is appropriate for an approximate solution of partial differential equations with arbitrary boundary shape. The method reduces the partial differential equation to a series of ordinary differential equations and is especially valuable for treating boundaries with nonlinear deformations. The calculations conclude that the pressure distribution and the poloidal current are more uniform in an oblate self-confined plasma ball than that of an elongated plasma ball. The ellipticity of the plasma ball is obviously restricted by its internal pressure, magnetic field, and ambient pressure. Qualitative evidence is presented for the absence of sighting of elongated ball lightning.

The thermal and mechanical structure of a two-dimensional plume in the earths mantle

On the condition that the distribution of velocity and temperature at the mid-plane of a mantle plume has been obtained (pages 213–218, this issue), the problem of determining the lateral structure of the plume at a given depth is reduced to solving an eigenvalue problem of a set of ordinary differential equations with five unknown functions, with an eigenvalue being related to the thermal thickness of the plume at this depth. The lateral profiles of upward velocity, temperature and viscosity in the plume and the thickness of the plume at various depths are calculated for two sets of Newtonian rheological parameters. The calculations show that the precondition for the existence of the plume, δT/L 1 (L = the height of the plume, δT = lateral distance from the mid-plane), can be satisfied, except for the starting region of the plume or near the base of the lithosphere. At the lateral distance, δT, the upward velocity decreases to 0.1 – 50% of its maximum value at different depths. It is believed that this model may provide an approach for a quantitative description of the detailed structure of a mantle plume.

Hybrid numerical methods for exponential models of growth

A series of novel numerical methods for the exponential models of growth are proposed. Based on these methods, hybrid predictor-corrector methods are constructed. The hybrid numerical methods can increase the accuracy and the computing speed obviously, as well as enlarge the stability domain greatly. (c) 2005 Published by Elsevier Inc.

On the soft wall guiding potentials in realistic quantum waveguides

A transfer matrix method is presented for the study of electron conduction in a quantum waveguide with soft wall lateral confinement. By transforming the two-dimensional Schrodinger equation into a set of second order ordinary differential equations, the total transfer matrix is obtained and the scattering probability amplitudes are calculated. The proposed method is applied to the evaluation of the electron transmission in two types of cavity structure with finite-height square-well confinement. The results obtained by our method, which are found to be in excellent agreement with those from another transfer matrix method, suggest that the infinite square-well potential is a good approximation to finite-height square-well confinement for electrons propagating in the ground transverse mode, but softening of the walls has an obvious effect on the electron transmission and mode-mixing for propagating in the excited transverse mode. (C) 1996 American Institute of Physics.

Analysis for temperature and pressure fields in process of hydrate dissociation by depressurization

An improved axisymmetric mathematic modeling is proposed for the process of hydrate dissociation by depressurization around vertical well. To reckon in the effect of latent heat of gas hydrate at the decomposition front, the energy balance equation is employed. The semi-analytic solutions for temperature and pressure fields are obtained by using Boltzmann-transformation. The location of decomposition front is determined by solving initial value problem for system of ordinary differential equations. The distributions of pressure and temperature along horizontal radiate in the reservoir are calculated. The numeric results indicate that the moving speed of decomposition front is sensitively dependent on the well pressure and the sediment permeability. Copyright (C) 2010 John Wiley & Sons, Ltd.

NONLINEAR GRAVITY-WAVES IN DEEP-WATER

In consideration of the problem on the boundary condition of nonlinear free water wave, coordinate transform is used to handle the free boundary. Supposing the solution form be the traveling wave, the ordinary differential equations of the one-order autonomous system with two variables are caused, then expanding the nonlinear terms at the equilibrium point with the Taylor expansion, we obtained the solution to traveling wave. The linear approximate equation near the equilibrium point is the small amplitude wave. A new nonlinear periodic traveling wave and nonlinear dispersion relation are shown when expanding to the second-order terms. A conclusion that the expansion of dispersion relation does not contain any odd-power terms of wave steepness and because of the nonlinear effort an oscillate structure is produced in the vertical direction is drawn.

Governing Equations and Numerical Solutions of Tension Leg Platform with Finite Amplitude Motion

It is demonstrated that when tension leg platform (TLP) moves with finite amplitude in waves, the inertia force, the drag force and the buoyancy acting on the platform are nonlinear functions of the response of TLP. The tensions of the tethers are also nonlinear functions of the displacement of TLP. Then the displacement, the velocity and the acceleration of TLP should be taken into account when loads are calculated. In addition, equations of motions should be set up on the instantaneous position. A theoretical model for analyzing the nonlinear behavior of a TLP with finite displacement is developed, in which multifold nonlinearities are taken into account, i.e., finite displacement, coupling of the six degrees of freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force. Based on the theoretical model, the comprehensive nonlinear differential equations are deduced. Then the nonlinear dynamic analysis of ISSC TLP in regular waves is performed in the time domain. The degenerative linear solution of the proposed nonlinear model is verified with existing published one. Furthermore, numerical results are presented, which illustrate that nonlinearities exert a significant influence on the dynamic responses of the TLP.

First-passage failure of quasi-integrable Hamiltonian systems

The first-passage failure of quasi-integrable Hamiltonian si-stems (multidegree-of-freedom integrable Hamiltonian systems subject to light dampings and weakly random excitations) is investigated. The motion equations of such a system are first reduced to a set of averaged Ito stochastic differential equations by using the stochastic averaging method for quasi-integrable Hamiltonian systems. Then, a backward Kolmogorov equation governing the conditional reliability function and a set of generalized Pontryagin equations governing the conditional moments of first-passage time are established. Finally, the conditional reliability function, and the conditional probability density and moments of first-passage time are obtained by solving these equations with suitable initial and boundary conditions. Two examples are given to illustrate the proposed procedure and the results from digital simulation are obtained to verify the effectiveness of the procedure.

First-passage time of Duffing oscillator under combined harmonic and white-noise excitations

The first-passage time of Duffing oscillator under combined harmonic and white-noise excitations is studied. The equation of motion of the system is first reduced to a set of averaged Ito stochastic differential equations by using the stochastic averaging method. Then, a backward Kolmogorov equation governing the conditional reliability function and a set of generalized Pontryagin equations governing the conditional moments of first-passage time are established. Finally, the conditional reliability function, and the conditional probability density and moments of first-passage time are obtained by solving the backward Kolmogorov equation and generalized Pontryagin equations with suitable initial and boundary conditions. Numerical results for two resonant cases with several sets of parameter values are obtained and the analytical results are verified by using those from digital simulation.

An asymptotic solution of hydrodynamic equations at the mid-plane of a plume in the earths mantle

From the partial differential equations of hydrodynamics governing the movements in the Earth's mantle of a Newtonian fluid with a pressure- and temperature-dependent viscosity, considering the bilateral symmetry of velocity and temperature distributions at the mid-plane of the plume, an analytical solution of the governing equations near the mid-plane of the plume was found by the method of asymptotic analysis. The vertical distribution of the upward velocity, viscosity and temperature at the mid-plane, and the temperature excess at the centre of the plume above the ambient mantle temperature were then calculated for two sets of Newtonian rheological parameters. The results obtained show that the temperature at the mid-plane and the temperature excess are nearly independent of the rheological parameters. The upward velocity at the mid-plane, however, is strongly dependent on the rheological parameters.

Topology of toroidal helical fields in non-circular cross-sectional tokamaks

The ordinary differential magnetic field line equations are solved numerically; the tokamak magnetic structure is studied on Hefei Tokamak-7 Upgrade (HT-7U) when the equilibrium field with a monotonic q-profile is perturbed by a helical magnetic field. We find that a single mode (m, n) helical perturbation can cause the formation of islands on rational surfaces with q = m/n and q = (m +/- 1, +/- 2, +/- 3,...)/n due to the toroidicity and plasma shape (i.e. elongation and triangularity), while there are many undestroyed magnetic surfaces called Kolmogorov-Arnold-Moser (KAM) barriers on irrational surfaces. The islands on the same rational surface do not have the same size. When the ratio between the perturbing magnetic field B-r(r) and the toroidal magnetic field amplitude B(phi)0 is large enough, the magnetic island chains on different rational surfaces will overlap and chaotic orbits appear in the overlapping area, and the magnetic field becomes stochastic. It is remarkable that the stochastic layer appears first in the plasma edge region.

Static shape control of repetitive structures integrated with piezoelectric actuators

In this paper, several simplification methods are presented for shape control of repetitive structures such as symmetrical, rotational periodic, linear periodic, chain and axisymmetrical structures. Some special features in the differential equations governing these repetitive structures are examined by considering the whole structures. Based on the special properties of the governing equations, several methods are presented for simplifying their solution process. Finally, the static shape control of a cantilever symmetrical plate with piezoelectric actuator patches is demonstrated using the present simplification method. The result shows that present methods can effectively be used to find the optimal control voltage for shape control.

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.