On Multi‐component Ermakov Systems in a Two‐Layer Fluid: Integrable Hamiltonian Structures and Exact Vortex Solutions

By introducing an elliptic vortex ansatz, the 2+1‐dimensional two‐layer fluid system is reduced to a finite‐dimensional nonlinear dynamical system. Time‐modulated variables are then introduced and multicomponent Ermakov systems are isolated. The latter is shown to be also Hamiltonian, thereby admitting general solutions in terms of an elliptic integral representation. In particular, a subclass of vortex solutions is obtained and their behaviors are simulated. Such solutions have recently found applications in oceanic and atmospheric dynamics. Moreover, it is proved that the Hamiltonian system is equivalent to the stationary nonlinear cubic Schrödinger equations coupled with a Steen‐Ermakov‐Pinney equation.     

Symbolic computation of high order compact difference schemes for three dimensional linear elliptic partial differential equations with variable coefficients

We present a symbolic computation procedure for deriving various high order compact difference approximation schemes for certain three dimensional linear elliptic partial differential equations with variable coefficients. Based on the Maple software package, we approximate the leading terms in the truncation error of the Taylor series expansion of the governing equation and obtain a 19 point fourth order compact difference scheme for a general linear elliptic partial differential equation. A test problem is solved numerically to validate the derived fourth order compact difference scheme. This symbolic derivation method is simple and can be easily used to derive high order difference approximation schemes for other similar linear elliptic partial differential equations.     

A variable gravity model of the otolith membrane

This paper describes a finite difference model of the otolith membrane which allows the acceleration due to gravity to vary, thus simulating conditions of gravity on the lunar and planetary surfaces.The differential coefficients of the second-order system of elliptic partial differential equations governing the steady state displacements of points of the membrane are replaced by finite differences. The resulting system of difference equations is seen to be consistent with the system of differential equations and to have a truncation error of order four.A close approximation to the physical boundary of a typical otolith membrane is used and two sets of numerical experiments are carried out which which simulate rotations of the membrane on the Moon and on a number of planets.The displacements at thirty nodes of the membrane are computed by solving the linear system of sixty equations obtained by applying the difference equations to each of the thirty nodes. The numerical results obtained are seen to be in general agreement with experimental results reported in the literature.     

Numerical simulations for the contraction flow using grid generation

We study the incomprssible Navier Stokes equations for the flow inside contraction geometry. The governing equations are expressed in the vorticity-stream function formulations. A rectangular computational domain is arised by elliptic grid generation technique. The numerical solution is based on a technique of automatic numerical generation of acurvilinear coordinate system by transforming the governing equation into computational plane. The transformed equations are approximated using central differences and solved simultaneously by successive over relaxation iteration. The time dependent of the vorticity equation solved by using explicit marching procedure. We will apply the technique on several irregularshapes.     

Error estimations for some meshless boundary interpolation methods

Meshless methods have become quite popular in numerical treatment of partial differential equations because of their simplicity and the fact that they require neither domain nor boundary mesh. In general, however, they convert the original problem to a highly ill‐conditioned linear system of algebraic equations with a dense matrix. Recently, a special technique has been proposed which circumvents this computational difficulty. This method, called Direct Multi‐Elliptic Interpolation Method, is based on a scattered data interpolation which defines the interpolation function as a solution of a higher order multi‐elliptic equation. Here the boundary version of this meshless method which is based on a multi‐elliptic boundary interpolation is considered. Error estimations are derived justifying the interpolation function to be a good approximation of the solution of the original boundary value problem as well. At the same time, the problem of large, dense and ill‐conditioned matrices as well as the mesh generation are completely avoided. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Optimal finite difference grids and rational approximations of the square root I. Elliptic problems

The main objective of this paper is optimization of second‐order finite difference schemes for elliptic equations, in particular, for equations with singular solutions and exterior problems. A model problem corresponding to the Laplace equation on a semi‐infinite strip is considered. The boundary impedance (Neumann‐to‐Dirichlet map) is computed as the square root of an operator using the standard three‐point finite difference scheme with optimally chosen variable steps. The finite difference approximation of the boundary impedance for data of given smoothness is the problem of rational approximation of the square root on the operator's spectrum. We have implemented Zolotarev's optimal rational approx‐imant obtained in terms of elliptic functions. We have also found that a geometrical progression of the grid steps with optimally chosen parameters is almost as good as the optimal approximant. For bounded operators it increases from second to exponential the convergence order of the finite difference impedance with the convergence rate proportional to the inverse of the logarithm of the condition number. For the case of unbounded operators in Sobolev spaces associated with elliptic equations, the error decays as the exponential of the square root of the mesh dimension. As an example, we numerically compute the Green function on the boundary for the Laplace equation. Some features of the optimal grid obtained for the Laplace equation remain valid for more general elliptic problems with variable coefficients. © 2000 John Wiley & Sons, Inc.     

Multiple vortices for a self-dual CP(1) Maxwell-Chern-Simons model

We prove the existence of at least two doubly periodic vortex solutions for a self-dual CP(1) Maxwell-Chern-Simons model. To this end we analyze a system of two elliptic equations with exponential nonlinearities. Such a system is shown to be equivalent to a fourth-order elliptic equation admitting a variational structure. Tonia Ricciardi: Partially supported by the MIUR National Project Variational Methods and Nonlinear Differential Equations     

Weak oscillations of a gas bubble in a spherical volume of compressible liquid

The following spherically symmetric problem is considered: a single gas bubble at the centre of a spherical flask filled with a compressible liquid is oscillating in response to forced radial excitation of the flask walls. In the long-wave approximation at low Mach numbers, one obtains a system of differential-difference equations generalizing the Rayleigh-Lamb-Plesseth equation. This system takes into account the compressibility of the liquid and is suitable for describing both free and forced oscillations of the bubble. It includes an ordinary differential equation analogous to the Herring-Flinn-Gilmore equation describing the evolution of the bubble radius, and a delay equation relating the pressure at the flask walls to the variation of the bubble radius. The solutions of this system of differential-difference equations are analysed in the linear approximation and numerical analysis is used to study various modes of weak but non-linear oscillations of the bubble, for different laws governing the variation of the pressure or velocity of the liquid at the flask wall. These solutions are compared with numerical solutions of the complete system of partial differential equations for the radial motion of the compressible liquid around the bubble.     

An efficient meshfree point collocation moving least squares method to solve the interface problems with nonhomogeneous jump conditions

We are going to study a simple and effective method for the numerical solution of the closed interface boundary value problem with both discontinuities in the solution and its derivatives. It uses a strong‐form meshfree method based on the moving least squares (MLS) approximation. In this method, for the solution of elliptic equation, the second‐order derivatives of the shape functions are needed in constructing the global stiffness matrix. It is well‐known that the calculation of full derivatives of the MLS approximation, especially in high dimensions, is quite costly. In the current work, we apply the diffuse derivatives using an efficient technique. In this technique, we calculate the higher‐order derivatives using the approximation of lower‐order derivatives, instead of calculating directly derivatives. This technique can improve the accuracy of meshfree point collocation method for interface problems with nonhomogeneous jump conditions and can efficiently estimate diffuse derivatives of second‐ and higher‐orders using only linear basis functions. To introduce the appropriate discontinuous shape functions in the vicinity of interface, we choose the visibility criterion method that modifies the support of weight function in MLS approximation and leads to an efficient computational procedure for the solution of closed interface problems. The proposed method is applied for elliptic and biharmonic interface problems. For the biharmonic equation, we use a mixed scheme, which replaces this equation by a coupled elliptic system. Also the application of the present method to elasticity equation with discontinuities in the coefficients across a closed interface has been provided. Representative numerical examples demonstrate the accuracy and robustness of the proposed methodology for the closed interface problems. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1031–1053, 2015     

半导体问题的交替方向有限元方法

该文用交替方向有限元方法求解半导体问题的Energy Trans port (ET)模型。对模型中椭圆型的电子位势方程采用交替方向迭代法，对流占优扩散的电子浓度和空穴浓度方程采用特征交替方向有限元方法，热传导方程利用Patch逼近采用交替方向有限元方法求解。利用微分方程的先验估计理论和技巧，分别得到了椭圆型方程和抛物型方程的最优H+1和L+2误差估计。     

Study of the contraction flow using grid generation technique

We study a model procedure to solve the incompressible Navier-Stokes equations on the flow inside contraction geometry. The governing equations are expressed in the primitive variable formulation. A rectangular computational plane is arises by elliptic grid generation technique. The numerical solution is based on a technique of automatic numerical generation of a curvilinear coordinate system. By transformed the governing equation into computational plane. The time dependent momentum equations are solved explicitly for the velocity field using the explicit marching procedure, the continuity equation is applied at each grid point in the solution of pressure equation, while the successive over relaxation (SOR) method is used for the Neumann problem for pressure. We will apply the technique on several irregular-shape.     

A power penalty approach to a mixed quasilinear elliptic complementarity problem

In this paper, a power penalty approximation method is proposed for solving a mixed quasilinear elliptic complementarity problem. The mixed complementarity problem is first reformulated as a double obstacle quasilinear elliptic variational inequality problem. A nonlinear elliptic partial differential equation is then defined to approximate the resulting variational inequality by using a power penalty approach. The existence and uniqueness of the solution to the partial differential penalty equation are proved. It is shown that, under some mild assumptions, the sequence of solutions to the penalty equations converges to the unique solution of the variational inequality problem as the penalty parameter tends to infinity. The error estimates of the convergence of this penalty approach are also derived. At last, numerical experimental results are presented to show that the power penalty approximation method is efficient and robust.

     

A pendulum with an elliptic-type parametric excitation: Stability charts for a damped and undamped system

In this paper, a pendulum parametrically excited by the excitation which has the form of the Jacobi cn elliptic function is considered. Three cases related to the value of the elliptic parameter are distinguished: the case when it is smaller than zero, when it ranges between zero and unity, and when it is higher than unity. First, interpretations of the excitation with such elliptic parameter are given in terms of its period, higher harmonic content and the amplitude. These interpretations enable one to consider the elliptic-type excitation as a type of multi-cosine excitation whose frequency and amplitude are related mutually in a particular way. Stability charts are determined for damped and undamped systems. When the elliptic parameter is equal to zero, the governing equations considered transform to the well-known Mathieu equation. In all other cases, the governing equations considered can be seen as a new generalisation of the Mathieu equation. The influence of an arbitrary real elliptic parameter on the location and shape of the transition curves and instability tongues is investigated, illustrated and discussed in all three cases, which represent new and so far unknown results.     

A semi-analytic approach for the nonlinear dynamic response of circular plates

This paper presents a new semi-analytic perturbation differential quadrature method for geometrically nonlinear vibration analysis of circular plates. The nonlinear governing equations are converted into a linear differential equation system by using Linstedt–Poincaré perturbation method. The solutions of nonlinear dynamic response and the nonlinear free vibration are then sought through the use of differential quadrature approximation in space domain and analytical series expansion in time domain. The present method is validated against analytical results using elliptic function in several examples for both clamped and simply supported circular plates, showing that it has excellent accuracy and convergence. Compared with numerical methods involving iterative time integration, the present method does not suffer from error accumulation and is able to give very accurate results over a long time interval.     

Unified Orbital Description of the Envelope Dynamics in Two‐Dimensional Simple Periodic Lattices

The propagation of wave envelopes in two‐dimensional (2‐D) simple periodic lattices is studied. A discrete approximation, known as the tight‐binding (TB) approximation, is employed to find the equations governing a class of nonlinear discrete envelopes in simple 2‐D periodic lattices. Instead of using Wannier function analysis, the orbital approximation of Bloch modes that has been widely used in the physical literature, is employed. With this approximation the Bloch envelope dynamics associated with both simple and degenerate bands are readily studied. The governing equations are found to be discrete nonlinear Schrödinger (NLS)‐type equations or coupled NLS‐type systems. The coefficients of the linear part of the equations are related to the linear dispersion relation. When the envelopes vary slowly, the continuous limit of the general discrete NLS equations are effective NLS equations in moving frames. These continuous NLS equations (from discrete to continuous) also agree with those derived via a direct multiscale expansion. Rectangular and triangular lattices are examples.     

A linear,decoupled fractional time‐stepping method for the nonlinear fluid–fluid interaction

In this paper, a linear decoupled fractional time stepping method is proposed and developed for the nonlinear fluid–fluid interaction governed by the two Navier–Stokes equations. Partitioned time stepping method is applied to two‐physics problems with stiffness of the coupling terms being treated explicitly and is also unconditionally stable. As for each fluid, the velocity and pressure are respectively determined by just solving one vector‐valued quasi‐elliptic equation and the Possion equation with homogeneous Neumann boundary condition per time step. Therefore, the cost of the fluid–fluid interaction is dominant to solve four simple linear equations, which greatly reduces the computational cost of the whole system. The method exploits properties of the fluid–fluid system to establish its stability and convergence with the same results as the standard scheme. Finally, numerical experiments are presented to show the performance of the proposed method.     

Traveling wave solutions for fractional partial differential equations arising in mathematical physics by an improved fractional Jacobi elliptic equation method

In this paper, combining with a new generalized ansätz and the fractional Jacobi elliptic equation, an improved fractional Jacobi elliptic equation method is proposed for seeking exact solutions of space‐time fractional partial differential equations. The fractional derivative used here is the modified Riemann‐Liouville derivative. For illustrating the validity of this method, we apply it to solve the space‐time fractional Fokas equation and the the space‐time fractional BBM equation. As a result, some new general exact solutions expressed in various forms including the solitary wave solutions, the periodic wave solutions, and Jacobi elliptic functions solutions for the two equations are found with the aid of mathematical software Maple. Copyright © 2016 John Wiley & Sons, Ltd.     

A symmetric integrated radial basis function method for solving differential equations

In this article, integrated radial basis functions (IRBFs) are used for Hermite interpolation in the solution of differential equations, resulting in a new meshless symmetric RBF method. Both global and local approximation‐based schemes are derived. For the latter, the focus is on the construction of compact approximation stencils, where a sparse system matrix and a high‐order accuracy can be achieved together. Cartesian‐grid‐based stencils are possible for problems defined on nonrectangular domains. Furthermore, the effects of the RBF width on the solution accuracy for a given grid size are fully explored with a reasonable computational cost. The proposed schemes are numerically verified in some elliptic boundary‐value problems governed by the Poisson and convection‐diffusion equations. High levels of the solution accuracy are obtained using relatively coarse discretisations.     

On the hydrodynamic limit of Ginzburg‐Landau wave vortices

In this paper, we study the hydrodynamic limit of the finite Ginzburg‐Landau wave vortices, which was established in [16]. Unlike the classical vortex method for incompressible Euler equations, we prove here that the densities approximated by the vortex blob method associated with the Ginzburg‐Landau wave vortices tend to the solutions of the pressure‐less compressible Euler‐Poisson equations. The convergence of such approximation is proved before the formation of singularities in the limit system as the blob sizes and the grid sizes tend to zero in appropriate rates. © 2002 John Wiley & Sons, Inc.     

Operator adapted spectral element methods I: harmonic and generalized harmonic polynomials

Summary. The paper presents results on the approximation of functions which solve an elliptic differential equation by operator adapted systems of functions. Compared with standard polynomials, these operator adapted systems have superior local approximation properties. First, the case of Laplace's equation and harmonic polynomials as operator adapted functions is analyzed and rates of convergence in a Sobolev space setting are given for the approximation with harmonic polynomials. Special attention is paid to the approximation of singular functions that arise typically in corners. These results for harmonic polynomials are extended to general elliptic equations with analytic coefficients by means of the theory of Bergman and Vekua; the approximation results for Laplace's equation hold true verbatim, if harmonic polynomials are replaced with generalized harmonic polynomials. The Partition of Unity Method is used in a numerical example to construct an operator adapted spectral method for Laplace's equation that is based on approximating with harmonic polynomials locally. Received May 26, 1997 / Revised version received September 21, 1998 / Published online September 7, 1999