Convergence of a Penalty-Finite Element Approximation for an Obstacle Problem.
An equilibrium triangular block-element, proposed by Watwood and Hartz, is subjected to an analysis and its approximability property is proved. If the solution is regular enough, a quasi-optimal error estimate follows for the dual approximation to the mixed boundary value problem of elasticity (based on Castigliano's principle). The convergence is proved even in a general case, when the solution is not regular.
We formulate a finite element method for the computation of solutions to an anisotropic phase-field model for a binary alloy. Convergence is proved in the -norm. The convergence result holds for anisotropy below a certain threshold value. We present some numerical experiments verifying the theoretical results. For anisotropy below the threshold value we observe optimal order convergence, whereas in the case where the anisotropy is strong the numerical solution to the phase-field equation does not...
The problem of modeling acoustic waves scattered by an object with Neumann boundary condition is considered. The boundary condition is taken into account by means of the fictitious domain method, yielding a first order in time mixed variational formulation for the problem. The resulting system is discretized with two families of mixed finite elements that are compatible with mass lumping. We present numerical results illustrating that the Neumann boundary condition on the object is not always correctly...
The problem of modeling acoustic waves scattered by an object with Neumann boundary condition is considered. The boundary condition is taken into account by means of the fictitious domain method, yielding a first order in time mixed variational formulation for the problem. The resulting system is discretized with two families of mixed finite elements that are compatible with mass lumping. We present numerical results illustrating that the Neumann boundary condition on the object is not always...
Discrete-velocity approximations represent a popular way for computing the Boltzmann collision operator. The direct numerical evaluation of such methods involve a prohibitive cost, typically O(N2d + 1) where d is the dimension of the velocity space. In this paper, following the ideas introduced in [C. Mouhot and L. Pareschi, C. R. Acad. Sci. Paris Sér. I Math. 339 (2004) 71–76, C. Mouhot and L. Pareschi, Math. Comput. 75 (2006) 1833–1852], we derive fast summation techniques for the evaluation of...
We develop a method for counting number of cells and extraction of approximate cell centers in 2D and 3D images of early stages of the zebra-fish embryogenesis. The approximate cell centers give us the starting points for the subjective surface based cell segmentation. We move in the inner normal direction all level sets of nuclei and membranes images by a constant speed with slight regularization of this flow by the (mean) curvature. Such multi- scale evolutionary process is represented by a geometrical...
Elastic two-layer curved composite beam with partial shear interaction is considered. It is assumed that each curved layer separately follows the Euler-Bernoulli hypothesis and the load slip relation for the flexible shear connection is a linear relationship. The curved composite beam at one of the end cross sections is fixed and the other end cross section is subjected by a concentrated radial load. Two cases are considered. In the first case the loaded end cross section is closed by a rigid plate...