On the use of high order difference methods for the system of one space second order nonlinear hyperbolic equations with variable coefficients

Implicit difference schemes of O(k⁴ + k²h² + h⁴), where k0, h 0 are grid sizes in time and space coordinates respectively, are developed for the efficient numerical integration of the system of one space second order nonlinear hyperbolic equations with variable coefficients subject to appropriate initial and Dirichlet boundary conditions. The proposed difference method for a scalar equation is applied for the wave equation in cylindrical and spherical symmetry. The numerical examples are given to illustrate the fourth order convergence of the methods.     

The exact solution of a linear integral equation with weakly singular kernel

A space , which is proved to be a reproducing kernel space with simple reproducing kernel, is defined. The expression of its reproducing kernel function is given. Subsequently, a class of linear Volterra integral equation (VIE) with weakly singular kernel is discussed in the new reproducing kernel space. The reproducing kernel method of linear operator equation Au=f, which request the image space of operator A is and operator A is bounded, is improved. Namely, the request for the image space is weakened to be L²[a,b], and the boundedness of operator A is also not required. As a result, the exact solution of the equation is obtained. The numerical experiments show the efficiency of our method.     

Fredholm-Volterra integral equation with singular kernel

The purpose of this paper is to obtain the solution of Fredholm-Volterra integral equation with singular kernel in the space L₂(?1, 1) × C(0,T), 0 ≤t ≤T < ∞, under certain conditions. The numerical method is used to solve the Fredholm integral equation of the second kind with weak singular kernel using the Toeplitz matrices. Also, the error estimate is computed and some numerical examples are computed using the MathCad package.     

Equiconvergence of expansions in eigenfunctions of Sturm-Liouville operators with distributional potentials in Hölder spaces

We establish the equiconvergence of expansions of an arbitrary function in the class L ₂(0, π) in the Fourier series in sines and in the Fourier series in the eigenfunctions of the first boundary value problem for the one-dimensional Schrödinger operator with a nonclassical potential. The equiconvergence is studied in the norm of the Hölder space. The potential is the derivative of a function that belongs to a fractional-order Sobolev space.     

Newton's method for fractal approximation

The problem of fitting a given function in theL ^q norm with a function generated by an iterated function system can be rapidly solved by applying Newton's method on the parameter space of the iterated function system. The key to this is a method for calculating the derivatives of a potential function with respect to the parameters.     

Error estimates for the interpolating moving least-squares method in n-dimensional space

In this paper, the interpolating moving least-squares (IMLS) method is discussed in details. A simpler expression of the approximation function of the IMLS method is obtained. Compared with the moving least-squares (MLS) approximation, the shape function of the IMLS method satisfies the property of Kronecker δ function. Then the meshless method based on the IMLS method can overcome the difficulties of applying the essential boundary conditions. The error estimates of the approximation function and its first and second order derivatives of the IMLS method are presented in n-dimensional space. The theoretical results show that if the weight function is sufficiently smooth and the order of the polynomial basis functions is big enough, the approximation function and its partial derivatives are convergent to the exact values in terms of the maximum radius of the domains of influence of nodes. Then the interpolating element-free Galerkin (IEFG) method based on the IMLS method is presented for potential problems. The advantage of the IEFG method is that the essential boundary conditions can be applied directly and easily. For the purpose of demonstration, some selected numerical examples are given to prove the theories in this paper.     

Transformation methods for finding multiple roots of nonlinear equations

In this paper, to estimate a multiple root p of an equation f(x) = 0, we transform the function f(x) to a hyper tangent function combined with a simple difference formula whose value changes from −1 to 1 as x passes through the root p. Then we apply the so-called numerical integration method to the transformed equation, which may result in a specious approximate root. Furthermore, in order to enhance the accuracy of the approximation we propose a Steffensen-type iterative method, which does not require any derivatives of f(x) nor is quite affected by an initial approximation. It is shown that the convergence order of the proposed method becomes cubic by simultaneous approximation to the root and its multiplicity. Results for some numerical examples show the efficiency of the new method.     

Numerical modelling of aeroelastic behaviour of an airfoil in viscous incompressible flow

In this paper, the problem of the numerical approximation of a two-dimensional incompressible viscous fluid flow interacting with a flexible structure is considered. Due to high Reynolds numbers in the range 10⁴ − 10⁶ the turbulent character of the flow is considered and modelled with the aid of Reynolds equations coupled with the k − ω turbulence model. The structure motion is described by a system of ordinary differential equations for three degrees of freedom: vertical displacement, rotation and rotation of the aileron. The problem is discretized in space by the Galerkin Least-Squares stabilized finite element method and the computational domain is treated with the aid of Arbitrary Lagrangian Eulerian method.     

GLOBAL EXISTENCE, UNIFORM DECAY AND EXPONENTIAL GROWTH FOR A CLASS OF SEMI-LINEAR WAVE EQUATION WITH STRONG DAMPING

In this paper, we consider the nonlinearly damped semi-linear wave equation associated with initial and Dirichlet boundary conditions. We prove the existence of a local weak solution and introduce a family of potential wells and discuss the invariants and vacuum isolating behavior of solutions. Furthermore, we prove the global existence of solutions in both cases which are polynomial and exponential decay in the energy space respectively, and the asymptotic behavior of solutions for the cases of potential well family with 0相似文献   

Discrete spectrum of a model Schr?dinger operator on the half-plane with Neumann conditions

We study the eigenpairs of a model Schr?dinger operator with a quadratic potential and Neumann boundary conditions on a half-plane. The potential is degenerate in the sense that it reaches its minimum all along a line that makes the angle ?? with the boundary of the half-plane. We show that the first eigenfunctions satisfy localization properties related to the distance to the minimum line of the potential. We investigate the densification of the eigenvalues below the essential spectrum in the limit ?? ?? 0, and we prove a full asymptotic expansion for these eigenvalues and their associated eigenvectors. We conclude the paper by numerical experiments obtained by a finite element method. The numerical results confirm and enlighten the theoretical approach.     

A numerical comparison of Chebyshev methods for solving fourth order semilinear initial boundary value problems

In solving semilinear initial boundary value problems with prescribed non-periodic boundary conditions using implicit-explicit and implicit time stepping schemes, both the function and derivatives of the function may need to be computed accurately at each time step. To determine the best Chebyshev collocation method to do this, the accuracy of the real space Chebyshev differentiation, spectral space preconditioned Chebyshev tau, real space Chebyshev integration and spectral space Chebyshev integration methods are compared in the L² and W^2,2 norms when solving linear fourth order boundary value problems; and in the L^∞([0,T];L²) and L^∞([0,T];W^2,2) norms when solving initial boundary value problems. We find that the best Chebyshev method to use for high resolution computations of solutions to initial boundary value problems is the spectral space Chebyshev integration method which uses sparse matrix operations and has a computational cost comparable to Fourier spectral discretization.     

Efficient long-time computations of time-domain boundary integrals for 2D and dissipative wave equation

Linear hyperbolic partial differential equations in a homogeneous medium, e.g., the wave equation describing the propagation and scattering of acoustic waves, can be reformulated as time-domain boundary integral equations. We propose an efficient implementation of a numerical discretization of such equations when the strong Huygens’ principle does not hold.For the numerical discretization, we make use of convolution quadrature in time and standard Galerkin boundary element method in space. The quadrature in time results in a discrete convolution of weights W_j with the boundary density evaluated at equally spaced time points. If the strong Huygens’ principle holds, W_j converge to 0 exponentially quickly for large enough j. If the strong Huygens’ principle does not hold, e.g., in even space dimensions or when some damping is present, the weights are never zero, thereby presenting a difficulty for efficient numerical computation.In this paper we prove that the kernels of the convolution weights approximate in a certain sense the time domain fundamental solution and that the same holds if both are differentiated in space. The tails of the fundamental solution being very smooth, this implies that the tails of the weights are smooth and can efficiently be interpolated. Further, we hint on the possibility to apply the fast and oblivious convolution quadrature algorithm of Schädle et al. to further reduce memory requirements for long-time computation. We discuss the efficient implementation of the whole numerical scheme and present numerical experiments.     

A novel finite volume method for the Riesz space distributed-order advection–diffusion equation

In this paper, we investigate the finite volume method (FVM) for a distributed-order space-fractional advection–diffusion (AD) equation. The mid-point quadrature rule is used to approximate the distributed-order equation by a multi-term fractional model. Next, the transformed multi-term fractional equation is solved by discretizing in space by the finite volume method and in time using the Crank–Nicolson scheme. We use a novel technique to deal with the convection term, by which the Riesz fractional derivative of order 0 < γ < 1 is transformed into a fractional integral form. An important contribution of our work is the use of nodal basis function to derive the discrete form of our model. The unique solvability of the scheme is also discussed and we prove that the Crank–Nicolson scheme is unconditionally stable and convergent with second-order accuracy. Finally, we give some examples to show the effectiveness of the numerical method.     

Galerkin finite element approximation of symmetric space-fractional partial differential equations

In this paper, symmetric space-fractional partial differential equations (SSFPDE) with the Riesz fractional operator are considered. The SSFPDE is obtained from the standard advection-dispersion equation by replacing the first-order and second-order space derivatives with the Riesz fractional derivatives of order 2β ∈ (0, 1) and 2α ∈ (1, 2], respectively. We prove that the variational solution of the SSFPDE exists and is unique. Using the Galerkin finite element method and a backward difference technique, a fully discrete approximating system is obtained, which has a unique solution according to the Lax-Milgram theorem. The stability and convergence of the fully discrete schemes are derived. Finally, some numerical experiments are given to confirm our theoretical analysis.     

A new off-step high order approximation for the solution of three-space dimensional nonlinear wave equations

In this paper, we propose a new high accuracy numerical method of O(k² + k²h² + h⁴) based on off-step discretization for the solution of 3-space dimensional non-linear wave equation of the form u_tt = A(x,y,z,t)u_xx + B(x,y,z,t)u_yy + C(x,y,z,t)u_zz + g(x,y,z,t,u,u_x,u_y,u_z,u_t), 0 < x,y,z < 1,t > 0 subject to given appropriate initial and Dirichlet boundary conditions, where k > 0 and h > 0 are mesh sizes in time and space directions respectively. We use only seven evaluations of the function g as compared to nine evaluations of the same function discussed in  and . We describe the derivation procedure in details of the algorithm. The proposed numerical algorithm is directly applicable to wave equation in polar coordinates and we do not require any fictitious points to discretize the differential equation. The proposed method when applied to a telegraphic equation is also shown to be unconditionally stable. Comparative numerical results are provided to justify the usefulness of the proposed method.     

On the intrinsic and the spatial numerical range

For a bounded function f from the unit sphere of a closed subspace X of a Banach space Y, we study when the closed convex hull of its spatial numerical range W(f) is equal to its intrinsic numerical range V(f). We show that for every infinite-dimensional Banach space X there is a superspace Y and a bounded linear operator such that . We also show that, up to renormig, for every non-reflexive Banach space Y, one can find a closed subspace X and a bounded linear operator T∈L(X,Y) such that .Finally, we introduce a sufficient condition for the closed convex hull of the spatial numerical range to be equal to the intrinsic numerical range, which we call the Bishop-Phelps-Bollobás property, and which is weaker than the uniform smoothness and the finite-dimensionality. We characterize strong subdifferentiability and uniform smoothness in terms of this property.     

The exponential dichotomy on general approximation scheme

This paper is devoted to the numerical analysis of abstract parabolic problems u′(t) = Au(t), u(0) = u ⁰ with hyperbolic generator A. We develop a general approach to establish a discrete dichotomy in a very general setting in the case of discrete approximation in space and time. It is a well-known fact that the phase space in the neighborhood of the hyperbolic equilibrium can be split in such a way that the original initial value problem is reduced to initial value problems with exponentially decaying solutions in opposite time directions. We use the theory of compact approximation principle and collectively condensing approximation to show that such a decomposition of the flow persists under rather general approximation schemes. The main assumption of our results are naturally satisfied, in particular, for operators with compact resolvents and condensing semigroups and can be verified for the finite element method as well as finite difference methods.     

Compact supported endographs and fuzzy sets

A metric is defined on a space of functions from a locally compact metric space X into the unit interval I in terms of the Hausdorff metric distance between their compact supported endographs in X × I. Convergence in this metric is shown to be equivalent to the conjunction of the Hausdorff metric convergence of supports in X and two conditions involving numerical values of the functions. The space of nonempty compact subsets of X with the Hausdorff metric is imbedded in the above function space by the characteristic function on subsets of X. Applications of these results to fuzzy subsets of X and fuzzy dynamical systems on X are indicated.