Minimal surfaces: a geometric three dimensional segmentation approach

Summary. A novel geometric approach for three dimensional object segmentation is presented. The scheme is based on geometric deformable surfaces moving towards the objects to be detected. We show that this model is related to the computation of surfaces of minimal area (local minimal surfaces). The space where these surfaces are computed is induced from the three dimensional image in which the objects are to be detected. The general approach also shows the relation between classical deformable surfaces obtained via energy minimization and geometric ones derived from curvature flows in the surface evolution framework. The scheme is stable, robust, and automatically handles changes in the surface topology during the deformation. Results related to existence, uniqueness, stability, and correctness of the solution to this geometric deformable model are presented as well. Based on an efficient numerical algorithm for surface evolution, we present a number of examples of object detection in real and synthetic images. Received January 4, 1996 / Revised version received August 2, 1996     

An efficient linear second order unconditionally stable direct discretization method for the phase-field crystal equation on surfaces

We develop an unconditionally stable direct discretization scheme for solving the phase-field crystal equation on surfaces. The surface is discretized by using an unstructured triangular mesh. Gradient, divergence, and Laplacian operators are defined on triangular meshes. The proposed numerical method is second-order accurate in space and time. At each time step, the proposed computational scheme results in linear elliptic equations to be solved, thus it is easy to implement the algorithm. It is proved that the proposed scheme satisfies a discrete energy-dissipation law. Therefore, it is unconditionally stable. A fast and efficient biconjugate gradients stabilized solver is used to solve the resulting discrete system. Numerical experiments are conducted to demonstrate the performance of the proposed algorithm.     

Fully discrete energy stable scheme for a phase-field moving contact line model with variable densities and viscosities

In this study, we propose a fully discrete energy stable scheme for the phase-field moving contact line model with variable densities and viscosities. The mathematical model comprises a Cahn–Hilliard equation, Navier–Stokes equation, and the generalized Navier boundary condition for the moving contact line. A scalar auxiliary variable is employed to transform the governing system into an equivalent form, thereby allowing the double well potential to be treated semi-explicitly. A stabilization term is added to balance the explicit nonlinear term originating from the surface energy at the fluid–solid interface. A pressure stabilization method is used to decouple the velocity and pressure computations. Some subtle implicit–explicit treatments are employed to deal with convention and stress terms. We establish a rigorous proof of the energy stability for the proposed time-marching scheme. A finite difference method based on staggered grids is then used to spatially discretize the constructed time-marching scheme. We also prove that the fully discrete scheme satisfies the discrete energy dissipation law. Our numerical results demonstrate the accuracy and energy stability of the proposed scheme. Using our numerical scheme, we analyze the contact line dynamics based on a shear flow-driven droplet sliding case. Three-dimensional droplet spreading is also investigated based on a chemically patterned surface. Our numerical simulation accurately predicts the expected energy evolution and it successfully reproduces the expected phenomena where an oil droplet contracts inward on a hydrophobic zone and then spreads outward rapidly on a hydrophilic zone.     

Decoupled energy‐law preserving numerical schemes for the Cahn–Hilliard–Darcy system

We study two novel decoupled energy‐law preserving and mass‐conservative numerical schemes for solving the Cahn‐Hilliard‐Darcy system which models two‐phase flow in porous medium or in a Hele–Shaw cell. In the first scheme, the velocity in the Cahn–Hilliard equation is treated explicitly so that the Darcy equation is completely decoupled from the Cahn–Hilliard equation. In the second scheme, an intermediate velocity is used in the Cahn–Hilliard equation which allows for the decoupling. We show that the first scheme preserves a discrete energy law with a time‐step constraint, while the second scheme satisfies an energy law without any constraint and is unconditionally stable. Ample numerical experiments are performed to gauge the efficiency and robustness of our scheme. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 936–954, 2016     

Regularization of a discrete scheme governing the three-dimensional evolution of a fluid-fluid interface

A regularized discrete scheme is developed that describes the three-dimensional evolution of the interface between fluids with different viscosities and densities in the Leibenzon-Muskat model. The regularization is achieved by smoothing the kernel of the singular integral involved in the differential equation governing the moving interface. The discrete scheme is tested by solving the problem of a drop of one fluid evolving in a translational flow of another.     

一种求解运动曲面上对流扩散方程的三维水平集方法

给出了一种求解运动曲面上对流扩散方程的三维水平集算法. 水平集函数被用来表示曲面.曲面上的微分方程及其解通过水平集方法被延拓到包含曲面的一个小邻域中. 一种半隐式的Crank-Nicholson 格式被用来做时间推进, 中心差分和三阶加权实质无振荡(WENO) 格式被分别用来离散方程中的扩散项和对流项. 分析证明了它在标准的Courant-Friedrichs-Lewy (CFL) 条件下的稳定性. 数值算例显示了它能取得二阶精度.     

Error analysis of a fully discrete projection method for magnetohydrodynamic system

In this paper, we develop and analyze a finite element projection method for magnetohydrodynamics equations in Lipschitz domain. A fully discrete scheme based on Euler semi-implicit method is proposed, in which continuous elements are used to approximate the Navier–Stokes equations and H ( curl ) conforming Nédélec edge elements are used to approximate the magnetic equation. One key point of the projection method is to be compatible with two different spaces for calculating velocity, which leads one to obtain the pressure by solving a Poisson equation. The results show that the proposed projection scheme meets a discrete energy stability. In addition, with the help of a proper regularity hypothesis for the exact solution, this paper provides a rigorous optimal error analysis of velocity, pressure and magnetic induction. Finally, several numerical examples are performed to demonstrate both accuracy and efficiency of our proposed scheme.     

Error Estimate of a New Conservative Finite Difference Scheme for the Klein-Gordon-Dirac System

In this paper, we derive and analyze a conservative Crank-Nicolson-type finite difference scheme for the Klein-Gordon-Dirac (KGD) system. Differing from the derivation of the existing numerical methods given in literature where the numerical schemes are proposed by directly discretizing the KGD system, we translate the KGD equations into an equivalent system by introducing an auxiliary function, then derive a nonlinear Crank-Nicolson-type finite difference scheme for solving the equivalent system. The scheme perfectly inherits the mass and energy conservative properties possessed by the KGD, while the energy preserved by the existing conservative numerical schemes expressed by two-level's solution at each time step. By using energy method together with the 'cut-off' function technique, we establish the optimal error estimate of the numerical solution, and the convergence rate is $\mathcal{O}(τ^2 + h^2)$ in $l^∞$-norm with time step $τ$ and mesh size $h.$ Numerical experiments are carried out to support our theoretical conclusions.     

A Vertex-Centered and Positivity-Preserving Finite Volume Scheme for Two-Dimensional Three-Temperature Radiation Diffusion Equations on General Polygonal Meshes

Two-dimensional three-temperature (2-D 3-T) radiation diffusion equations are widely used to approximately describe the evolution of radiation energy within a multimaterial system and explain the exchange of energy among electrons, ions and photons. In this paper, we suggest a new positivity-preserving finite volume scheme for 2-D 3-T radiation diffusion equations on general polygonal meshes. The vertex unknowns are treated as primary ones for which the finite volume equations are constructed. The edge-midpoint and cell-centered unknowns are used as auxiliary ones and interpolated by the primary unknowns, which makes the final scheme a pure vertex-centered one. By comparison, most existing positivity-preserving finite volume schemes are cell-centered and based on the convex decomposition of the co-normal. Here, the co-normal decomposition is not convex in general, leading to a fixed stencil of the flux approximation and avoiding a certain search algorithm on complex grids. Moreover, the new scheme effectively alleviates the numerical heat-barrier issue suffered by most existing cell-centered or hybrid schemes in solving strongly nonlinear radiation diffusion equations. Numerical experiments demonstrate the second-order accuracy and the positivity of the solution on various distorted grids. For the problem without analytic solution, the contours of the numerical solutions obtained by our scheme on distorted meshes accord with those on smooth quadrilateral meshes.     

New conservative difference schemes for a coupled nonlinear Schrödinger system

In this paper, two conservative difference schemes for solving a coupled nonlinear Schrödinger (CNLS) system are numerically analyzed. Firstly, a nonlinear implicit two-level finite difference scheme for CNLS system is studied, then a linear three-level difference scheme for CNLS system is presented. An induction argument and the discrete energy method are used to prove the second-order convergence and unconditional stability of the linear scheme. Numerical examples show the efficiency of the new scheme.     

An extended tube model for thermo-viscoelasticity of rubberlike materials: Parameter identification and examples

The extended tube-model was presented by KALISKE & HEINRICH (RCT 72, 602-632) in 1999 as a novel approach for isothermal hyperelasticity of rubberlike materials. This contribution is dedicated to its further development to finite non-linear thermo-viscoelasticity. A non-linear evolution law and a thermo-mechanical coupled free energy formulation are the kernel of the phenomenological approach where the elastic material response is inspired by statistical-mechanical theory. The representation of viscoelasticity is based on a multiplicative decomposition of the deformation gradient. The Helmholtz free energy of the material is formulated in terms of isothermal free energy functions multiplicatively coupled with non-linear temperature evolution functions. The non-linear evolution law for the viscous material branch is solved by applying a predictor-corrector algorithm with an exponential mapping scheme. In today's literature, several sophisticated thermo-mechanical material models are available. However, they are built upon a considerable number of material parameters governing the mechanical and thermal material response which need to be identified for practical application. Therefore, particular emphasis is given to an appropriate parameter identification technique for the thermal field. For the latter, a uniaxial extension test is carried out where the recorded data of the temperature field of the rubber specimen under cyclic loading is used for parameter identification. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A new fourth‐order numerical algorithm for a class of three‐dimensional nonlinear evolution equations

In this article, a new compact alternating direction implicit finite difference scheme is derived for solving a class of 3‐D nonlinear evolution equations. By the discrete energy method, it is shown that the new difference scheme has good stability and can attain second‐order accuracy in time and fourth‐order accuracy in space with respect to the discrete H¹ ‐norm. A Richardson extrapolation algorithm is applied to achieve fourth‐order accuracy in temporal dimension. Numerical experiments illustrate the accuracy and efficiency of the extrapolation algorithm. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

An Eulerian–Lagrangian method for optimization problems governed by multidimensional nonlinear hyperbolic PDEs

We present a numerical method for solving tracking-type optimal control problems subject to scalar nonlinear hyperbolic balance laws in one and two space dimensions. Our approach is based on the formal optimality system and requires numerical solutions of the hyperbolic balance law forward in time and its nonconservative adjoint equation backward in time. To this end, we develop a hybrid method, which utilizes advantages of both the Eulerian finite-volume central-upwind scheme (for solving the balance law) and the Lagrangian discrete characteristics method (for solving the adjoint transport equation). Experimental convergence rates as well as numerical results for optimization problems with both linear and nonlinear constraints and a duct design problem are presented.     

Finite Difference Scheme for Barotropic Gas Equations

An implicit finite difference scheme approximating the equations of barotropic gas flow is proposed. This scheme ensures the positivity of density and the validity of an energy inequality and the mass conservation law. The continuity equation is approximated implicitly. It is proved that the resulting system of nonlinear equations has a solution for any time and space stepsizes. An iterative method for solving the system of nonlinear equations at each time step is proposed.     

Analysis of signal optimization problems for dielectric media

The authors have recently presented the results of the numerical implementation of a scheme for designing an incident electromagnetic plane wave which will deposit a maximum amount of energy in any selected portion of an inhomogeneous dielectric slab. This paper presents the mathematical analysis underlying this scheme for solving the signal design problem. The analysis indicates that the maximum energy is achieved by choosing the incident signal to be an appropriately normalized eigenfunction of a compact Fredholm integral operator. A numerical scheme for solving the optimization problem is discussed briefly.This research was supported by Ames Laboratory-USDOE under Contract No. W-7405-Eng-82.     

Perturbation method for solving the nonlinear eigenvalue problem arising from fatigue crack growth problem in a damaged medium

An analytical solution of the nonlinear eigenvalue problem arising from the fatigue crack growth problem in a damaged medium in coupled formulation is obtained. The perturbation technique for solving the nonlinear eigenvalue problem is used. The method allows to find the analytical formula expressing the eigenvalue as the function of parameters of the damage evolution law. It is shown that the eigenvalues of the nonlinear eigenvalue problem are fully determined by the exponents of the damage evolution law. In the paper the third-order (four-term) asymptotic expansions of the angular functions determining the stress and continuity fields in the neighborhood of the crack tip are given. The asymptotic expansions of the angular functions permit to find the closed-form solution for the problem considered.     

A conservative extension of the method of characteristics for 1-D shallow flows

The method of characteristics (MOC) has been used for a long time in open channels and pipes flows. It is based on non-conservative equations, and hence it cannot be used directly for solving discontinuous shallow flows. In this paper we develop a conservative version of the MOC scheme for 1-D shallow flows by imposing the conservation law at the interpolation step. The conservation property of the scheme ensures the production of an accurate shock modeling and enables the MOC scheme to simulate dam-break type flows. By using a proper interpolation function, the proposed method can also produce quite accurate low-oscillatory results. A number of challenging test cases show considerable improvement compared to the traditional non-conservative MOC scheme in the case of dam-break type and trans-critical flow simulations.     

The Isoenergetic KAM-Type Theorem at Resonant Case for Nearly Integrable Hamiltonian Systems

In this paper, we study the persistence of resonant invariant tori on energy surfaces for nearly integrable Hamiltonian systems under the usual R$\ddot{u}$ssmann nondegenerate condition. By a quasilinear iterative scheme, we prove the following things: (1) The majority of resonant tori on a given energy surface will be persisted under R$\ddot{u}$ssmann nondegenerate condition. (2) The maximal number about the preserved frequency components on a perturbed torus is characterized by the smaller of the maximal rank of the Hessian matrices of the unperturbed system and the nondegeneracy of resonance. (3) If unperturbed systems admit subisoenergetic nondegeneracy on an energy surface, then the majority of the unperturbed resonant tori on the energy surface will be persisted and give rise to a family of perturbed tori with the same energy, whose frequency ratios among respective ''nondegenerate'' components are preserved.     

An analysis of the Perona‐Malik scheme

We investigate how the Perona‐Malik scheme evolves piecewise smooth initial data in one dimension. By scaling a natural parameter that appears in the scheme in an appropriate way with respect to the grid size, we obtain a meaningful continuum limit. The resulting evolution can be seen as the gradient flow for an energy, just as the discrete evolutions are gradient flows for discrete energies. It involves, except at special isolated times, solving a system of heat equations coupled to each other through nonlinear boundary conditions. At the special times, the solutions experience gradient blowup; nevertheless, there is a natural continuation for the solutions beyond these singular times. © 2001 John Wiley & Sons, Inc.     

Convergence Analysis of a Numerical Scheme for the Porous Medium Equation by an Energetic Variational Approach

The porous medium equation (PME)is a typical nonlinear degenerate parabolic equation. We have studied numerical methods for PME by an energetic variational approach in [C. Duan et al., J. Comput. Phys., 385 (2019), pp. 13–32], where the trajectory equation can be obtained and two numerical schemes have been developed based on different dissipative energy laws. It is also proved that the nonlinear scheme, based on $f$log $f$ as the total energy form of the dissipative law, is uniquely solvable on an admissible convex set and preserves the corresponding discrete dissipation law. Moreover, under certain smoothness assumption, we have also obtained the second order convergence in space and the first order convergence in time for the scheme. In this paper, we provide a rigorous proof of the error estimate by a careful higher order asymptotic expansion and two step error estimates. The latter technique contains a rough estimate to control the highly nonlinear term in a discrete $W$^1,∞norm and a refined estimate is applied to derive the optimal error order.