A modified predictor–corrector scheme for the two-dimensional sine-Gordon equation

A three-time level finite-difference scheme based on a fourth order in time and second order in space approximation has been proposed for the numerical solution of the nonlinear two-dimensional sine-Gordon equation. The method, which is analysed for local truncation error and stability, leads to the solution of a nonlinear system. To avoid solving it, a predictor–corrector scheme using as predictor a second-order explicit scheme is proposed. The procedure of the corrector has been modified by considering as known the already evaluated corrected values instead of the predictor ones. This modified scheme has been tested on the line and circular ring soliton and the numerical experiments have proved that there is an improvement in the accuracy over the standard predictor–corrector implementation. This research was co-funded by E.U. (75%) and by the Greek Government (25%).     

An Efficient Numerical Model for Planar Spiral Inductor On-Chip

In this paper, a very useful numerical technique has been developed for analyzing the transient characteristics of a planar-spiral inductor on-chip. A locally conformal technique and an alternating-direction implicit scheme are applied to the finite-difference time-domain method. A formulation for solving three dimensional Maxwell’s equations is proposed. Using the proposed method, various parameters of the planar-spiral inductors have been analyzed and an equivalent circuit, which includes frequency-independent circuit elements, has been introduced. Highly computational efficiency is implemented. Numerical results show excellent agreement with the measured data over a wide frequency range.     

The Convergence and Stability of Splitting Finite-Difference Schemes for Nonlinear Evolutionary Equations

We consider a splitting finite-difference scheme for an initial-boundary value problem for a two-dimensional nonlinear evolutionary equation. The problem is split into nonlinear and linear parts. The linear part is also split into locally one-dimensional equations. We prove the convergence and stability of the scheme in L ₂ and C norms. Printed in Lietuvos Matematikos Rinkinys, Vol. 45, No. 3, pp. 413–434, July–September, 2005.     

Parametric Quintic-Spline Approach to the Solution of a System of Fourth-Order Boundary-Value Problems

In this paper, we use parametric quintic splines to derive some consistency relations which are then used to develop a numerical method for computing the solution of a system of fourth-order boundary-value problems associated with obstacle, unilateral, and contact problems. It is known that a class of variational inequalities related to contact problems in elastostatics can be characterized by a sequence of variational inequations, which are solved using some numerical method. Numerical evidence is presented to show the applicability and superiority of the new method over other collocation, finite difference, and spline methods.     

Optimized prefactored compact schemes

The numerical simulation of aeroacoustic phenomena requires high-order accurate numerical schemes with low dispersion and dissipation errors. In this paper we describe a strategy for developing high-order accurate prefactored compact schemes, requiring very small stencil support. These schemes require fewer boundary stencils and offer simpler boundary condition implementation than existing compact schemes. The prefactorization strategy splits the central implicit schemes into forward and backward biased operators. Using Fourier analysis, we show it is possible to select the coefficients of the biased operators such that their dispersion characteristics match those of the original central compact scheme and their numerical wavenumbers have equal and opposite imaginary components. This ensures that when the forward and backward stencils are added, the original central compact scheme is recovered. To extend the resolution characteristic of the schemes, an optimization strategy is employed in which formal order of accuracy is sacrificed in preference to enhanced resolution characteristics across the range of wavenumbers realizable on a given mesh. The resulting optimized schemes yield improved dispersion characteristics compared to the standard sixth- and eighth-order compact schemes making them more suitable for high-resolution numerical simulations in gas dynamics and computational aeroacoustics. The efficiency, accuracy and convergence characteristics of the new optimized prefactored compact schemes are demonstrated by their application to several test problems.     

Case of a Boltzmann gas leading to the Smoluchowski coagulation equation

A special model of a rarefied hard-sphere gas is considered. The hard-sphere particles undergo absolutely elastic collisions. It is assumed that particles can collide only if their nonzero velocities are orthogonal to each other. The model makes it possible to proceed from the Boltzmann equation to the Smoluchowski coagulation equation, where coagulation means that the kinetic energies of the colliding particles are added. A Monte Carlo scheme for simulation of the phenomenon is described, and the convergence of the simulation algorithm is proved. The convergence of numerical results to exact solutions of the Smoluchowski equation and to finite-difference solutions is tested.     

Numerical differentiation from a viewpoint of regularization theory

In this paper, we discuss the classical ill-posed problem of numerical differentiation, assuming that the smoothness of the function to be differentiated is unknown. Using recent results on adaptive regularization of general ill-posed problems, we propose new rules for the choice of the stepsize in the finite-difference methods, and for the regularization parameter choice in numerical differentiation regularized by the iterated Tikhonov method. These methods are shown to be effective for the differentiation of noisy functions, and the order-optimal convergence results for them are proved.

     

Decay of confined,two-dimensional,spatially periodic arrays of vortices: A numerical investigation

The disarrangement of a perturbed lattice of vortices was studied numerically. The basic state is an exponentially decaying, exact solution of the Navier-Stokes equations. Square arrays of vortices with even numbers of vortex cells along each side were perturbed and their evolution was investigated. Whether the energy in the perturbation grows somewhat before it decays or decays monotonically depends on the initial strength of the vortices of the basic state, the extent of lateral confinement and the structure of the perturbation. The critical condition for temporally local instability, i.e. the critical amplitude of the basic state that must be exceeded to allow energy transfer from the basic state to the perturbation, is discussed. In the strongly confined case of a square lattice of four vortices the appearance of enchancement of global rotation is the result of energy transfer from the basic state to a temporally local unstable mode. Energy is transferred from the basic state to larger-scaled structures (inverse cascade) only if the scales of the larger structures are inherently contained in the initial structure of the perturbation. The initial structure of the double array of vortices is not maintained except for a very special form of perturbation. The facts that large scales decay more slowly than small scales and that, when non-linearities are sufficiently strong, energy is transferred from one scale to another explain the differences in the disarrangement process for different initial strengths of the vortices of the basic state. The stronger vortices, i.e. the vortices perturbed in a manner that increases their strength, tend to dominate the weaker vortices. The pairing and subsequent merging (or capture) of vortices of like sense into larger-scale vortices are described in terms of peaks in the evolution of the square root of the palinstrophy divided by the enstrophy.     

利用时域有限差分法模拟金纳米球及球壳的光学特性

利用时域有限差分方法研究了金纳米球、金纳米球壳及多层球壳的消光特性及电场分布.结果表明:金纳米颗粒的几何参量对消光峰的位置有显著影响.随着SiO2核心半径的增大,金纳米壳的消光峰显著红移.随着金核心半径的增大,gold-silica-gold多层球壳消光谱的低能峰显著红移,而高能峰微弱蓝移.     

Uniform stable conformal convolutional perfectly matched layer for enlarged cell technique conformal finite-difference time-domain method

Based on conformal construction of physical model in a three-dimensional Cartesian grid,an integral-based conformal convolutional perfectly matched layer(CPML) is given for solving the truncation problem of the open port when the enlarged cell technique conformal finite-difference time-domain(ECT-CFDTD) method is used to simulate the wave propagation inside a perfect electric conductor(PEC) waveguide.The algorithm has the same numerical stability as the ECT-CFDTD method.For the long-time propagation problems of an evanescent wave in a waveguide,several numerical simulations are performed to analyze the reflection error by sweeping the constitutive parameters of the integral-based conformal CPML.Our numerical results show that the integral-based conformal CPML can be used to efficiently truncate the open port of the waveguide.