Boundary knot method for 2D and 3D Helmholtz and convection-diffusion problems under complicated geometry

The boundary knot method (BKM) of very recent origin is an inherently meshless, integration-free, boundary-type, radial basis function collocation technique for the numerical discretization of general partial differential equation systems. Unlike the method of fundamental solutions, the use of non-singular general solution in the BKM avoids the unnecessary requirement of constructing a controversial artificial boundary outside the physical domain. The purpose of this paper is to extend the BKM to solve 2D Helmholtz and convection-diffusion problems under rather complicated irregular geometry. The method is also first applied to 3D problems. Numerical experiments validate that the BKM can produce highly accurate solutions using a relatively small number of knots. For inhomogeneous cases, some inner knots are found necessary to guarantee accuracy and stability. The stability and convergence of the BKM are numerically illustrated and the completeness issue is also discussed.

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.

Possibility of lifetime measurement in a flowing medium with phase fluorometry

The possibility of lifetime measurement in a flowing medium with phase fluorometry is investigated theoretically. A 3-D time dependent partial differential equation of the number density of atoms (or molecules) in the upper level of the fluorescence transition is solved analytically, taking flow, diffusion, optical excitation, decay, Doppler shift, and thickness of the excitation light sheet into account. An analytical expression of the intensity of the fluorescence signal in the flowing medium is deduced. Conditions are given, in which the principle of lifetime measurement with phase fluorometry in the static sample cell can be used in a flowing medium.

A Unified Tensor Level Set for Image Segmentation

This paper presents a new region-based unified tensor level set model for image segmentation. This model introduces a three-order tensor to comprehensively depict features of pixels, e.g., gray value and the local geometrical features, such as orientation and gradient, and then, by defining a weighted distance, we generalized the representative region-based level set method from scalar to tensor. The proposed model has four main advantages compared with the traditional representative method as follows. First, involving the Gaussian filter bank, the model is robust against noise, particularly the salt-and pepper-type noise. Second, considering the local geometrical features, e. g., orientation and gradient, the model pays more attention to boundaries and makes the evolving curve stop more easily at the boundary location. Third, due to the unified tensor pixel representation representing the pixels, the model segments images more accurately and naturally. Fourth, based on a weighted distance definition, the model possesses the capacity to cope with data varying from scalar to vector, then to high-order tensor. We apply the proposed method to synthetic, medical, and natural images, and the result suggests that the proposed method is superior to the available representative region-based level set method.

A Relay Level Set Method for Automatic Image Segmentation

This paper presents a new image segmentation method that applies an edge-based level set method in a relay fashion. The proposed method segments an image in a series of nested subregions that are automatically created by shrinking the stabilized curves in their previous subregions. The final result is obtained by combining all boundaries detected in these subregions. The proposed method has the following three advantages: 1) It can be automatically executed without human-computer interactions; 2) it applies the edge-based level set method with relay fashion to detect all boundaries; and 3) it automatically obtains a full segmentation without specifying the number of relays in advance. The comparison experiments illustrate that the proposed method performs better than the representative level set methods, and it can obtain similar or better results compared with other popular segmentation algorithms.

基于SVD和能量最小原则的图像自适应降噪算法

基于奇异值分解和能量最小原则,提出了一种自适应图像降噪算法,并给出了基于有界变差的能量降噪模型的代数形式。通过在矩阵范数意义下求能量最小,自适应确定去噪图像重构的奇异值个数。该算法的特点是将能量最小法则和奇异值分解结合起来,在代数空间中建立了一种自适应的图像降噪算法。与基于压缩比和奇异值分解的降噪方法相比,由于该算法避免了图像压缩比函数及其拐点的计算,因此具有快速去噪和简单可行的优点。实验结果证明,该算法是有效的。

Solitary magnon excitations in a one-dimensional antiferromagnet with Dzyaloshinsky-Moriya interactions

We investigate solitary excitations in a model of a one-dimensional antiferromagnet including a single-ion anisotropy and a Dzyaloshinsky-Moriya antisymmetric exchange interaction term. We employ the Holstein-Primakoff transformation, the coherent state ansatz and the time variational principle. We obtain two partial differential equations of motion by using the method of multiple scales and applying perturbation theory. By so doing, we show that the motion of the coherent amplitude must satisfy the nonlinear Schrodinger equation. We give the single-soliton solution.

A Study of Composite Beam with Shape Memory Alloy Arbitrarily Embedded under Thermal and Mechanical Loadings

The constitutive relations and kinematic assumptions on the composite beam with shape memory alloy (SMA) arbitrarily embedded are discussed and the results related to the different kinematic assumptions are compared. As the approach of mechanics of materials is to study the composite beam with the SMA layer embedded, the kinematic assumption is vital. In this paper, we systematically study the kinematic assumptions influence on the composite beam deflection and vibration characteristics. Based on the different kinematic assumptions, the equations of equilibrium/motion are different. Here three widely used kinematic assumptions are presented and the equations of equilibrium/motion are derived accordingly. As the three kinematic assumptions change from the simple to the complex one, the governing equations evolve from the linear to the nonlinear ones. For the nonlinear equations of equilibrium, the numerical solution is obtained by using Galerkin discretization method and Newton-Rhapson iteration method. The analysis on the numerical difficulty of using Galerkin method on the post-buckling analysis is presented. For the post-buckling analysis, finite element method is applied to avoid the difficulty due to the singularity occurred in Galerkin method. The natural frequencies of the composite beam with the nonlinear governing equation, which are obtained by directly linearizing the equations and locally linearizing the equations around each equilibrium, are compared. The influences of the SMA layer thickness and the shift from neutral axis on the deflection, buckling and post-buckling are also investigated. This paper presents a very general way to treat thermo-mechanical properties of the composite beam with SMA arbitrarily embedded. The governing equations for each kinematic assumption consist of a third order and a fourth order differential equation with a total of seven boundary conditions. Some previous studies on the SMA layer either ignore the thermal constraint effect or implicitly assume that the SMA is symmetrically embedded. The composite beam with the SMA layer asymmetrically embedded is studied here, in which symmetric embedding is a special case. Based on the different kinematic assumptions, the results are different depending on the deflection magnitude because of the nonlinear hardening effect due to the (large) deflection. And this difference is systematically compared for both the deflection and the natural frequencies. For simple kinematic assumption, the governing equations are linear and analytical solution is available. But as the deflection increases to the large magnitude, the simple kinematic assumption does not really reflect the structural deflection and the complex one must be used. During the systematic comparison of computational results due to the different kinematic assumptions, the application range of the simple kinematic assumption is also evaluated. Besides the equilibrium study of the composite laminate with SMA embedded, the buckling, post-buckling, free and forced vibrations of the composite beam with the different configurations are also studied and compared.

Statistical and stochastic description of microvoid evolution observed in dynamic failure of ductile materials

A brief review is presented of statistical approaches on microdamage evolution. An experimental study of statistical microdamage evolution in two ductile materials under dynamic loading is carried out. The observation indicates that there are large differences in size and distribution of microvoids between these two materials. With this phenomenon in mind, kinetic equations governing the nucleation and growth of microvoids in nonlinear rate-dependent materials are combined with the balance law of void number to establish statistical differential equations that describe the evolution of microvoids' number density. The theoretical solution provides a reasonable explanation of the experimentally observed phenomenon. The effects of stochastic fluctuation which is influenced by the inhomogeneous microscopic structure of materials are subsequently examined (i.e. stochastic growth model). Based on the stochastic differential equation, a Fokker-Planck equation which governs the evolution of the transition probability is derived. The analytical solution for the transition probability is then obtained and the effects of stochastic fluctuation is discussed. The statistical and stochastic analyses may provide effective approaches to reveal the physics of damage evolution and dynamic failure process in ductile materials.

Tensionless Contact Of A Finite Beam Resting On Reissner Foundation

A more generalized model of a beam resting on a tensionless Reissner foundation is presented. Compared with the Winkler foundation model, the Reissner foundation model is a much improved one. In the Winkler foundation model, there is no shear stress inside the foundation layer and the foundation is assumed to consist of closely spaced, independent springs. The presence of shear stress inside Reissner foundation makes the springs no longer independent and the foundation to deform as a whole. Mathematically, the governing equation of a beam on Reissner foundation is sixth order differential equation compared with fourth order of Winkler one. Because of this order change of the governing equation, new boundary conditions are needed and related discussion is presented. The presence of the shear stress inside the tensionless Reissner foundation together with the unknown feature of contact area/length makes the problem much more difficult than that of Winkler foundation. In the model presented here, the effects of beam dimension, gap distance, loading asymmetry and foundation shear stress on the contact length are all incorporated and studied. As the beam length increases, the results of a finite beam with zero gap distance converge asymptotically to those obtained by the previous model for an infinitely long beam. (C) 2008 Elsevier Ltd. All rights reserved.

A theoretical analysis on stochastic behaviour of short cracks and fatigue damage of metals

In this paper, the closed form of solution to the stochastic differential equation for a fatigue crack evolution system is derived. and the relationship between metal fatigue damage and crack stochastic behaviour is investigated. It is found that the damage extent of metals is independent of crack stochastic behaviour ii the stochastic deviation of the crack growth rate is directly proportional to its mean value. The evolution of stochastic deviation of metal fatigue damage in the stage close to the transition point between short and long crack regimes is also discussed.

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.

FORTRAN语言的非数值应用

<正> 1.系统设计 利用高速计算机,进行微分方程公式的推导是数学研究工作的一个新的尝试,它将影响今后微分方程研究工作的现代化。 微分方程程序系统(Differential Equation Program System,简称DEPS)的主要用途是