Discontinuous cellular automaton method for crack growth analysis without remeshing

A numerical technical of discontinuous cellular automaton method for crack growth analysis without remeshing is developed. In this method, the level set method is employed to track the crack location and its growth path, where the level set functions and calculation grids are independent, so no explicit meshing for crack surface and no remeshing for crack growth are needed. Then, the discontinuous enrichment shape functions which are enriched by the Heaviside function and the exact near-tip asymptotic field functions are constructed to model the discontinuity of cracks. Finally, a discontinuous cellular automaton theory is proposed, which are composed of cell, neighborhood and updating rules for discontinuous case. There is an advantage that the calculation is only applied on local cell, so no assembled stiffness matrix but only cell stiffness is needed, which can overcome the stiffness matrix assembling difficulty caused by unequal degrees of nodal freedom for different cells, and much easier to consider the local properties of cells. Besides, the present method requires much less computer memory than that of XFEM because of it local property.     

二维Qausi-Geostrophic方程的几何约束与非爆炸性

研究二维Qausi-Geostrophic(QG)方程水平集的动态演化过程.在"涡线"上"涡度"的最大值与全局"涡度"最大值可比的假设下,得到若"涡线"的长度L(t)被O(ln ln ||▽~⊥θ(·,t)||_∞)控制(被O(ln ln ||▽~⊥θ(·,t)||_∞)控制),曲率最大值K(t)与L(t)的乘积被O(ln ln ln ||▽~⊥θ(·,t)||_∞)控制(被O(ln ln ||▽~⊥θ(·,t)||_∞)控制),则||▽~⊥θ(·,t)||_∞)的增长阶数最多可达四重指数阶,这样二维QG方程在有限时间内无爆炸发生.     

Modeling shear modulus distribution in magnetic resonance elastography with piecewise constant level sets

Magnetic resonance elastography (MRE) is designed for imaging the mechanical properties of soft tissues. However, the interpretation of shear modulus distribution is often confusing and cumbersome. For reliable evaluation, a common practice is to specify the regions of interest and consider regional elasticity. Such an experience-dependent protocol is susceptible to intrapersonal and interpersonal variability. In this study we propose to remodel shear modulus distribution with piecewise constant level sets by referring to the corresponding magnitude image. Optimal segmentation and registration are achieved by a new hybrid level set model comprised of alternating global and local region competitions. Experimental results on the simulated MRE data sets show that the mean error of elasticity reconstruction is 11.33% for local frequency estimation and 18.87% for algebraic inversion of differential equation. Piecewise constant level set modeling is effective to improve the quality of shear modulus distribution, and facilitates MRE analysis and interpretation.     

固定转发源级的目标回波强度测量法

本文首先介绍了两种传统的水中目标回波强度的测量方法,分析了其特点并指出其不足之处。然后提出了一种改进的回声转发测量法,我们称之为固定转发源级的目标回波强度测量法。固定转发源级的目标回波测量法,对传统的回声转发测量法进行了硬件结构和处理方法上的改进。通过在被测目标上加装声学数据采集记录装置,避免了变强度转发的需要。不但保留了其优点,而且克服了其主要缺点。明显降低了校准的工作量,提高了测量的精度和效率,拓展了系统的适用范围,具有一定的应用推广价值。     

Effect of a uniform magnetic field on dielectric two-phase bubbly flows using the level set method

In this study, the behavior of a single bubble in a dielectric viscous fluid under a uniform magnetic field has been simulated numerically using the Level Set method in two-phase bubbly flow. The two-phase bubbly flow was considered to be laminar and homogeneous. Deformation of the bubble was considered to be due to buoyancy and magnetic forces induced from the external applied magnetic field. A computer code was developed to solve the problem using the flow field, the interface of two phases, and the magnetic field. The Finite Volume method was applied using the SIMPLE algorithm to discretize the governing equations. Using this algorithm enables us to calculate the pressure parameter, which has been eliminated by previous researchers because of the complexity of the two-phase flow. The finite difference method was used to solve the magnetic field equation. The results outlined in the present study agree well with the existing experimental data and numerical results. These results show that the magnetic field affects and controls the shape, size, velocity, and location of the bubble.     

Variational piecewise constant level set methods for shape optimization of a two-density drum

We apply the piecewise constant level set method to a class of eigenvalue related two-phase shape optimization problems. Based on the augmented Lagrangian method and the Lagrange multiplier approach, we propose three effective variational methods for the constrained optimization problem. The corresponding gradient-type algorithms are detailed. The first Uzawa-type algorithm having applied to shape optimization in the literature is proven to be effective for our model, but it lacks stability and accuracy in satisfying the geometry constraint during the iteration. The two other novel algorithms we propose can overcome this limitation and satisfy the geometry constraint very accurately at each iteration. Moreover, they are both highly initial independent and more robust than the first algorithm. Without penalty parameters, the last projection Lagrangian algorithm has less severe restriction on the time step than the first two algorithms. Numerical results for various instances are presented and compared with those obtained by level set methods. The comparisons show effectiveness, efficiency and robustness of our methods. We expect our promising algorithms to be applied to other shape optimization and multiphase problems.     

History matching of facies distribution with the EnKF and level set parameterization

In this work, we develop a methodology to combine the Ensemble Kalman filter (EnKF) and the level set parameterization for history matching of facies distribution. With given prior knowledge about the facies of the reservoir geology, initial realizations are generated by commonly used software as the prior guesses of the unknown field. Furthermore, level set functions are used to reparameterize these initial realizations. In the reparameterization process, a representing node system is set up, on which the values of level set functions are assigned using Gaussian random numbers. The mean and the standard deviation of the Gaussian random numbers are designed according to the facies proportion, and the sign of the random numbers depends on the facies type at the representing nodes. The values of the level set functions at the other grid nodes are obtained by linear interpolation. The level set functions on the representing nodes are the model parameters of the EnKF state vector and are updated in the data assimilation process. On the basis of our numerical examples for two-dimensional reservoirs with two or three facies, the proposed method is demonstrated to be able to capture the main features of the reference facies distributions.     

The moving boundary node method: A level set-based,finite volume algorithm with applications to cell motility

Eukaryotic cell crawling is a highly complex biophysical and biochemical process, where deformation and motion of a cell are driven by internal, biochemical regulation of a poroelastic cytoskeleton. One challenge to built quantitative models that describe crawling cells is solving the reaction–diffusion–advection dynamics for the biochemical and cytoskeletal components of the cell inside its moving and deforming geometry. Here we develop an algorithm that uses the level set method to move the cell boundary and uses information stored in the distance map to construct a finite volume representation of the cell. Our method preserves Cartesian connectivity of nodes in the finite volume representation while resolving the distorted cell geometry. Derivatives approximated using a Taylor series expansion at finite volume interfaces lead to second order accuracy even on highly distorted quadrilateral elements. A modified, Laplacian-based interpolation scheme is developed that conserves mass while interpolating values onto nodes that join the cell interior as the boundary moves. An implicit time stepping algorithm is used to maintain stability. We use the algorithm to simulate two simple models for cellular crawling. The first model uses depolymerization of the cytoskeleton to drive cell motility and suggests that the shape of a steady crawling cell is strongly dependent on the adhesion between the cell and the substrate. In the second model, we use a model for chemical signalling during chemotaxis to determine the shape of a crawling cell in a constant gradient and to show cellular response upon gradient reversal.     

A level set method for the semiclassical limit of the Schrödinger equation with discontinuous potentials

We propose a level set method for the semiclassical limit of the Schrödinger equation with discontinuous potentials. The discontinuities in the potential corresponds to potential barriers, at which incoming waves can be partially transmitted and reflected. Previously such a problem was handled by Jin and Wen using the Liouville equation – which arises as the semiclassical limit of the Schrödinger equation – with an interface condition to account for partial transmissions and reflections (S. Jin, X. Wen, SIAM J. Num. Anal. 44 (2006) 1801–1828). However, the initial data are Dirac-delta functions which are difficult to approximate numerically with a high accuracy. In this paper, we extend the level set method introduced in (S. Jin, H. Liu, S. Osher, R. Tsai, J. Comp. Phys. 210 (2005) 497–518) for this problem. Instead of directly discretizing the Delta functions, our proposed method decomposes the initial data into finite sums of smooth functions that remain smooth in finite time along the phase flow, and hence can be solved much more easily using conventional high order discretization schemes.     

Z/(2e)上本原序列不同压缩映射的导出序列

设 f( x)是 Z/ ( 2 e)上 n次强本原多项式 ,对形如 xe- 1 +η( x0 ,… ,xe- 2 )的二个 e元布尔函数 Φ( x0 ,… ,xe- 1 )和 Ψ( x0 ,… ,xe- 1 )及二条序列 a,b∈G( f( x) ) e,若Φ( a0 ,… ,ae- 1 ) =Ψ ( b0 ,… ,be- 1 ) ,给出了函数Φ ( x0 ,… ,xe- 1 )和Ψ ( x0 ,… ,xe- 1 )之间的关系与序列 a和 b之间的关系 .所给出的结论进一步说明了导出的二元序列具有良好的密码性质