Benchmark Problems for the Numerical Discretization of the Cahn–Hilliard Equation with a Source Term

In this paper, we present benchmark problems for the numerical discretization of the Cahn–Hilliard equation with a source term. If the source term includes an isotropic growth term, then initially circular and spherical shapes should grow with their original shapes. However, there is numerical anisotropic error and this error results in anisotropic evolutions. Therefore, it is essential to use isotropic space discretization in the simulation of growth phenomenon such as tumor growth. To test numerical discretization, we present two benchmark problems: one is the growth of a disk or a sphere and the other is the growth of a rotated ellipse or a rotated ellipsoid. The computational results show that the standard discrete Laplace operator has severe grid orientation dependence. However, the isotropic discrete Laplace operator generates good results.

Accurate numerical scheme for singularly perturbed parabolic delay differential equation

Objectives Numerical treatment of singularly perturbed parabolic delay differential equation is considered. Solution of the equation exhibits a boundary layer, which makes it difficult for numerical computation. Accurate numerical scheme is proposed using $$\theta$$ θ -method in time discretization and non-standard finite difference method in space discretization. Result Stability and uniform convergence of the proposed scheme is investigated. The scheme is uniformly convergent with linear order of convergence before Richardson extrapolation and second order convergent after Richardson extrapolation. Numerical examples are considered to validate the theoretical findings.

A Mixed Finite Element-Based Numerical Method for Elastodynamics Considering Adhesive Interface Damage for Dynamic Fracture

The numerical solution of the elastodynamic problem with kinematic boundary conditions is considered using mixed finite elements for the space discretization and a staggered leap-frog scheme for the discretization in time. The stability of the numerical scheme is shown under the usual CFL condition. Using the general form of Robin-type boundary conditions some cases of debonding and the resulting acoustic emission are studied. The methodology presented finds applications to geophysics such as in seismic waves simulation with dynamic rupture and energy release. In this paper, we focus on problems of fracture occurring at the interface of composite materials. Our goal is to study both the mechanism of crack initiation and propagation, as well as the acoustic emission of the released structure-borne energy which may be used in structural health monitoring and prognosis applications.

Second Order Fully Semi-Lagrangian Discretizations of Advection-Diffusion-Reaction Systems

We propose a second order, fully semi-Lagrangian method for the numerical solution of systems of advection-diffusion-reaction equations, which is based on a semi-Lagrangian approach to approximate in time both the advective and the diffusive terms. The proposed method allows to use large time steps, while avoiding the solution of large linear systems, which would be required by implicit time discretization techniques. Standard interpolation procedures are used for the space discretization on structured and unstructured meshes. A novel extrapolation technique is proposed to enforce second-order accurate Dirichlet boundary conditions. We include a theoretical analysis of the scheme, along with numerical experiments which demonstrate the effectiveness of the proposed approach and its superior efficiency with respect to more conventional explicit and implicit time discretizations.

Time and Space Discretization in the Digital Twin of the Airport Transport Network

The paper describes the concept of a Digital Twin (DT) creation, which could be used as part of a centralized ground traffic control system in the airport’s apron area. The work is focused on the automatic control systems expected in the future, which will completely replace the ground traffic controllers that solve the centralized control’s tasks in our days. It is assumed that the developed DT can be used both at the stage of training the control system and at the stage of its use. A special feature of the proposed DT concept is the clearly formulated principles of the observed process’ time discretization and the space in which the movement of ground vehicles and aircraft occurs. The experimental verification of the developed concept was carried out using a special simulation model designed to generate data streams on the airport’s transport network’s traffic participants’ location.

Adaptive Euler methods for stochastic systems with non-globally Lipschitz coefficients

We present strongly convergent explicit and semi-implicit adaptive numerical schemes for systems of semi-linear stochastic differential equations (SDEs) where both the drift and diffusion are not globally Lipschitz continuous. Numerical instability may arise either from the stiffness of the linear operator or from the perturbation of the nonlinear drift under discretization, or both. Typical applications arise from the space discretization of an SPDE, stochastic volatility models in finance, or certain ecological models. Under conditions that include montonicity, we prove that a timestepping strategy which adapts the stepsize based on the drift alone is sufficient to control growth and to obtain strong convergence with polynomial order. The order of strong convergence of our scheme is (1 − ε)/2, for ε ∈ (0,1), where ε becomes arbitrarily small as the number of finite moments available for solutions of the SDE increases. Numerically, we compare the adaptive semi-implicit method to a fully drift-implicit method and to three other explicit methods. Our numerical results show that overall the adaptive semi-implicit method is robust, efficient, and well suited as a general purpose solver.

Matrix Equation Techniques for Certain Evolutionary Partial Differential Equations

We show that the discrete operator stemming from time-space discretization of evolutionary partial differential equations can be represented in terms of a single Sylvester matrix equation. A novel solution strategy that combines projection techniques with the full exploitation of the entry-wise structure of the involved coefficient matrices is proposed. The resulting scheme is able to efficiently solve problems with a tremendous number of degrees of freedom while maintaining a low storage demand as illustrated in several numerical examples.

Linearly implicit energy-preserving Fourier pseudospectral schemes for the complex modified Korteweg-de Vries equation

We propose two linearly implicit energy-preserving schemes for the complex modified Korteweg–de Vries equation, based on the invariant energy quadratization method. First, a new variable is introduced and a new Hamiltonian system is constructed for this equation. Then the Fourier pseudospectral method is used for the space discretization and the Crank–Nicolson leap-frog schemes for the time discretization. The proposed schemes are linearly implicit, which is only needed to solve a linear system at each time step. The fully discrete schemes can be shown to conserve both mass and energy in the discrete setting. Some numerical examples are also presented to validate the effectiveness of the proposed schemes. doi:10.1017/S1446181120000218

A Fully Discrete SIPG Method for Solving Two Classes of Vortex Dominated Flows

To simulate incompressible Navier–Stokes equation, a temporal splitting scheme in time and high-order symmetric interior penalty Galerkin (SIPG) method in space discretization are employed, while the local Lax-Friedrichs flux is applied in the discretization of the nonlinear term. Under a constraint of the Courant–Friedrichs–Lewy (CFL) condition, two benchmark problems in 2D are simulated by the fully discrete SIPG method. One is a lid-driven cavity flow and the other is a circular cylinder flow. For the former, we compute velocity field, pressure contour and vorticity contour. In the latter, while the von Kármán vortex street appears with Reynolds number 50≤Re≤400, we simulate different dynamical behavior of circular cylinder flows, and numerically estimate the Strouhal numbers comparable to the existing experimental results. The calculations on vortex dominated flows are carried out to investigate the potential application of the SIPG method.