Elliptic Object Detection Based on Shape Preserving and Active Contour

本文通过形状约束方程(组)与一般主动轮廓模型结合,将目标形状与主动轮廓模型融合到统一能量泛函模型中,提出了一种形状保持主动轮廓模型即曲线在演化过程中保持为某一类特定形状。模型通过参数化水平集函数的零水平集控制演化曲线形状,不仅达到了分割即目标的目的,而且能够给出特定目标的定量描述。根据形状保持主动轮廓模型,建立了一个用于椭圆状目标检测的统一能量泛函模型,导出了相应的Euler-Lagrange常微分方程并用水平集方法实现了椭圆状目标检测。此模型可以应用于眼底乳头分割,虹膜检测及相机标定。实验结果表明,此模型不仅能够准确的检测出给定图像中的椭圆状目标,而且有很强的抗噪、抗变形及遮挡性能。

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.

Some theorems in asymptotic analysis

Basing ourselves on the analysis of magnitude of order, we strictly prove fundamental lemmas for asymptotic integral, including the cases of infinite region. Then a general formula for asymptotic expansion of integrals is given. Finally, we derive a sufficient condition for an ordinary differential equation to possess a solution of the Frobenius series type at finite irregular singularities or branching points.

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.

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.

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.

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.

FORTRAN语言的非数值应用

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

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.