In earlier work, Maire developed a class of cell-centered Lagrangian schemes for solving Euler equations of compressible gas dynamics in cylindrical coordinates. These schemes use a node-based discretization of the numerical fluxes. The control volume version has several distinguished properties, including the conservation of mass, momentum and total energy and compatibility with the geometric conservation law (GCL). However it also has a limitation in that it cannot preserve spherical symmetry for one-dimensional spherical flow. An alternative is also given to use the first order area-weighted approach which can ensure spherical symmetry, at the price of sacrificing conservation of momentum. In this paper, we apply the methodology proposed in our recent work to the first order control volume scheme of Maire to obtain the spherical symmetry property. The modified scheme can preserve one-dimensional spherical symmetry in a two-dimensional cylindrical geometry when computed on an equal-angle-zoned initial grid, and meanwhile it maintains its original good properties such as conservation and GCL. Several two-dimensional numerical examples in cylindrical coordinates are presented to demonstrate the good performance of the scheme in terms of symmetry, non-oscillation and robustness properties.

A convex non-convex variational model is proposed for multiphase image segmentation. We consider a specially designed non-convex regularization term which adapts spatially to the image structures for better controlling of the segmentation and easy handling of the intensity inhomogeneities. The nonlinear optimization problem is effciently solved by an Alternating Directions Methods of Multipliers procedure. We provide a convergence analysis and perform numerical experiments on several images, showing the effectiveness of this procedure.

No abstract uploaded!

We propose a simple finite difference scheme for Navier--Stokes equations in primitive formulation on curvilinear domains. With proper boundary treatment and interplay between covariant and contravariant components, the spatial discretization admits exact Hodge decomposition and energy identity. As a result, the pressure can be decoupled from the momentum equation with explicit time stepping. No artificial pressure boundary condition is needed. In addition, it can be shown that this spatially compatible discretization leads to uniform inf-sup condition, which plays a crucial role in the pressure approximation of both dynamic and steady state calculations. Numerical experiments demonstrate the robustness and efficiency of our scheme.

The numerical solution of the one-dimensional nonlinear Klein-Gordon equation on an unbounded domain is studied in this paper. Split local absorbing boundary (SLAB) conditions are obtained by the operator splitting method, then the original problem is reduced to an initial boundary value problem on a bounded computational domain, which can be solved by the finite difference method. Several numerical examples are provided to show the advantages and effectiveness of the given method, and some interesting collision behaviors are also observed.