Mechanical quadrature methods and extrapolation for solving nonlinear boundary Helmholtz integral equations

This paper presents mechanical quadrature methods (MQMs) for solving nonlinear boundary Helmholtz integral equations. The methods have high accuracy of order O(h ³) and low computation complexity. Moreover, the mechanical quadrature methods are simple without computing any singular integration. A nonlinear system is constructed by discretizing the nonlinear boundary integral equations. The stability and convergence of the system are proved based on an asymptotical compact theory and the Stepleman theorem. Using the h ³-Richardson extrapolation algorithms (EAs), the accuracy to the order of O(h ⁵) is improved. To slove the nonlinear system, the Newton iteration is discussed extensively by using the Ostrowski fixed point theorem. The efficiency of the algorithms is illustrated by numerical examples.     

Kernel-independent fast multipole method within the framework of regularized Stokeslets

The method of regularized Stokeslets (MRS) uses a radially symmetric blob function of infinite support to smooth point forces and allows for evaluation of the resulting flow field. This is a common method to study swimmers at zero Reynolds number where the Stokeslet is the fundamental solution corresponding to the kernel of the single layer potential. Simulating the collective motion of N micro-swimmers using the MRS results in at least N² pair-wise interactions. Efficient simulation of a large number of swimmers in free space is observed with the implementation of the kernel-independent fast multipole method (FMM) for radial basis functions. We illustrate the complexity of the algorithm on a simple test case where we study regularized point forces, showing that the method is of order N. Additionally, we explore accuracy in time for the MRS where the swimmers are modeled as Kirchhoff rods and the kernel-independent FMM is compared to the direct calculation using the standard MRS. Optimal hydrodynamic efficiency is also explored for different configurations of swimmers.     

2D Burgers equation with large Reynolds number using POD/DEIM and calibration

Model order reduction of the two‐dimensional Burgers equation is investigated. The mathematical formulation of POD/discrete empirical interpolation method (DEIM)‐reduced order model (ROM) is derived based on the Galerkin projection and DEIM from the existing high fidelity‐implicit finite‐difference full model. For validation, we numerically compared the POD ROM, POD/DEIM, and the full model in two cases of Re = 100 and Re = 1000, respectively. We found that the POD/DEIM ROM leads to a speed‐up of CPU time by a factor of O(10). The computational stability of POD/DEIM ROM is maintained by means of a careful selection of POD modes and the DEIM interpolation points. The solution of POD/DEIM in the case of Re = 1000 has an accuracy with error O(10^?3) versus O(10^?4) in the case of Re = 100 when compared with the high fidelity model. For this turbulent flow, a closure model consisting of a Tikhonov regularization is carried out in order to recover the missing information and is developed to account for the small‐scale dissipation effect of the truncated POD modes. It is shown that the computational results of this calibrated ROM exhibit considerable agreement with the high fidelity model, which implies the efficiency of the closure model used. Copyright © 2016 John Wiley & Sons, Ltd.     

APPLICATION OF THE DUAL RECIPROCITY BOUNDARY ELEMENT METHOD TO SOLUTION OF NONLINEAR DIFFERENTIAL EQUATION

In this paper the dual reciprocity boundary element method is employed to solve nonlinear differential equation ∇² u+u+ɛu ³ =b. Results obtained in an example have a good agreement with those by FEM and show the applicability and simplicity of dual reciprocity method(DRM)in solving nonlinear differential equations.     

Differential-equation-based representation of truncation errors for accurate numerical simulation

High-order compact finite difference schemes for two-dimensional convection-diffusion-type differential equations with constant and variable convection coefficients are derived. The governing equations are employed to represent leading truncation terms, including cross-derivatives, making the overall O(h⁴) schemes conform to a 3 × 3 stencil. We show that the two-dimensional constant coefficient scheme collapses to the optimal scheme for the one-dimensional case wherein the finite difference equation yields nodally exact results. The two-dimensional schemes are tested against standard model problems, including a Navier-Stokes application. Results show that the two schemes are generally more accurate, on comparable grids, than O(h²) centred differencing and commonly used O(h) and O(h³) upwinding schemes.     

Finite element implementation of boundary conditions for the pressure Poisson equation of incompressible flow

In this paper we address the problem of the implementation of boundary conditions for the derived pressure Poisson equation of incompressible flow. It is shown that the direct Galerkin finite element formulation of the pressure Poisson equation automatically satisfies the inhomogeneous Neumann boundary conditions, thus avoiding the difficulty in specifying boundary conditions for pressure. This ensures that only physically meaningful pressure boundary conditions consistent with the Navier-Stokes equations are imposed. Since second derivatives appear in this formulation, the conforming finite element method requires C¹ continuity. However, for many problems of practical interest (i.e. high Reynolds numbers) the second derivatives need not be included, thus allowing the use of more conventional C⁰ elements. Numerical results using this approach for a wall-driven contained flow within a square cavity verify the validity of the approach. Although the results were obtained for a two-dimensional problem using the p-version of the finite element method, the approach presented here is general and remains valid for the conventional h-version as well as three-dimensional problems.     

Gross-Pitaevskii Equation as the Mean Field Limit of Weakly Coupled Bosons

We consider the dynamics of N boson systems interacting through a pair potential N ^?1 V _a (x _i ?x _j ) where V _a (x)=a ^?3 V(x/a). We denote the solution to the N-particle Schrödinger equation by Ψ _{N, t} . Recall that the Gross-Pitaevskii (GP) equation is a nonlinear Schrödinger equation and the GP hierarchy is an infinite BBGKY hierarchy of equations so that if u _t solves the GP equation, then the family of k-particle density matrices solves the GP hierarchy. Under the assumption that a=N ^?? for 0N→∞ the limit points of the k-particle density matrices of Ψ _{N, t} are solutions of the GP hierarchy with the coupling constant in the nonlinear term of the GP equation given by ∫V(x)dx. The uniqueness of the solutions of this hierarchy remains an open question.     

Determination of the Boundaries between the Domains of Stability and Instability for the Hill Equation

The stability problem for the Hill equation containing two parameters is analyzed using the Mathematica computer algebra system. The characteristic constant is found as a series expansion in powers of a small parameter e. It is shown that the domains of instability are located only between the curves a = a(e) on the a-e plane crossing the axis e = 0 at the points a = (2k – 1)² / 4, k = 1, 2, 3, ....The corresponding curves are found as power series in e with accuracy O(e ⁶).     

The asymptotic theory of semilinear perturbed telegraph equation and its application

I.IntroductionConsiderthefollowingsemilinearperturbedtelegraphequationuII-u., P'u==sj(t,x,u,ul,u,,e)(-ooo)(l.l)u(o,x)=u,(x,e)(-ooo,u=u(t,x),fuoandulsatisfycertainsmoo…     

BEM for transient 2D elastodynamics using multiquadric functions

A new boundary element model for transient dynamic analysis of 2D structures is presented. The dual reciprocity method (DRM) is reformulated for the 2D elastodynamics by using the multiquadric radial basis functions (MQ). The required kernels for displacement and traction particular solutions are derived. Some terms of these kernels are found to be singular; therefore, a new smoothing technique is proposed to solve this problem. Hence, the limiting values of relevant kernels are computed. The validity and strength of the proposed formulation are demonstrated throughout several numerical applications. It is proven from the results that the present formulation is more stable than the traditional DRM, which uses the conical (1 + R) function, especially in predicting results in the far time zone.     

The Laplace transform method for Burgers' equation

The Laplace transform method (LTM) is introduced to solve Burgers' equation. Because of the nonlinear term in Burgers' equation, one cannot directly apply the LTM. Increment linearization technique is introduced to deal with the situation. This is a key idea in this paper. The increment linearization technique is the following: In time level t, we divide the solution u(x, t) into two parts: u(x, t^k) and w(x, t), t^k⩽t⩽t^k+1, and obtain a time‐dependent linear partial differential equation (PDE) for w(x, t). For this PDE, the LTM is applied to eliminate time dependency. The subsequent boundary value problem is solved by rational collocation method on transformed Chebyshev points. To face the well‐known computational challenge represented by the numerical inversion of the Laplace transform, Talbot's method is applied, consisting of numerically integrating the Bromwich integral on a special contour by means of trapezoidal or midpoint rules. Numerical experiments illustrate that the present method is effective and competitive. Copyright © 2009 John Wiley & Sons, Ltd.     

An optimized Schwarz method with two‐sided Robin transmission conditions for the Helmholtz equation

Optimized Schwarz methods are working like classical Schwarz methods, but they are exchanging physically more valuable information between subdomains and hence have better convergence behaviour. The new transmission conditions include also derivative information, not just function values, and optimized Schwarz methods can be used without overlap. In this paper, we present a new optimized Schwarz method without overlap in the 2d case, which uses a different Robin condition for neighbouring subdomains at their common interface, and which we call two‐sided Robin condition. We optimize the parameters in the Robin conditions and show that for a fixed frequency an asymptotic convergence factor of 1 – O(h^1/4) in the mesh parameter h can be achieved. If the frequency is related to the mesh parameter h, h = O(1/ω^γ) for γ?1, then the optimized asymptotic convergence factor is 1 – O(ω^(1–2γ)/8). We illustrate our analysis with 2d numerical experiments. Copyright © 2007 John Wiley & Sons, Ltd.     

Nonstationary pandom vibpation analysis of linear elastic structures with finite element method

At present, the finite element method is an efficient method for analyzing structural dynamic problems. When the physical quantities such as displacements and stresses are resolved in the spectra and the dynamic matrices are obtained in spectral resolving form, the relative equations cannot be solved by the vibration mode resolving method as usual. For solving such problems, a general method is put forward in this paper. The excitations considered with respect to nonstationary processes are as follows: P(t)={P_i(t)},P_i(t)=a_i(t)P_i(t), a_i(t) is a time function already known. We make Fourier transformation for the discretized equations obtained by finite element method, and by utilizing the behaviour of orthogonal increment of spectral quantities in random process^[1], some formulas of relations about the spectra of excitation and response are derived. The cross power spectral denisty matrices of responses can be found by these formulas, then the structrual safety analysis can be made. When a_i(t)=l (i= 1,2,…n), the. method stated in this paper will be reduced to that which is used in the special case of stationary process.     

Comparison of three spectral methods for the Benjamin-Ono equation: Fourier pseudospectral, rational Christov functions and Gaussian radial basis functions

The Benjamin-Ono equation is especially challenging for numerical methods because (i) it contains the Hilbert transform, a nonlocal integral operator, and (ii) its solitary waves decay only as O(1/|x|²). We compare three different spectral methods for solving this one-space-dimensional equation. The Fourier pseudospectral method is very fast through use of the Fast Fourier Transform (FFT), but requires domain truncation: replacement of the infinite interval by a large but finite domain. Such truncation is unnecessary for a rational basis, but it is simple to evaluate the Hilbert Transform only when the usual rational Chebyshev functions TB_n(x) are replaced by their cousins, the Christov functions; the FFT still applies. Radial basis functions (RBFs) are slow for a given number of grid points N because of the absence of a summation algorithm as fast as the FFT; because RBFs are meshless, however, very flexible grid adaptation is possible.     

Extended Lagrangian formalism and the corresponding energy relations

In this paper an extended Lagrangian formalism for the rheonomic systems with the nonstationary constraints is formulated, with the aim to examine more completely the energy relations for such systems in any generalized coordinates, which in this case always refer to some moving frame of reference. Introducing new quantities, which change according to the law τ_a=φ_a(t), it is demonstrated that these quantities determine the position of this moving reference frame with respect to an immobile one. In the transition to the generalized coordinates qⁱ they are taken as the additional generalized coordinates q^a=τ_a, whose dependence on time is given a priori. In this way the position of the considered mechanical system relative to this immobile frame of reference is determined completely.Based on this and using the corresponding d'Alembert–Lagrange's principle, an extended system of the Lagrangian equations is obtained. It is demonstrated that they give the same equations of motion qⁱ=qⁱ(t) as in the usual Lagrangian formulation, but substantially different energy relations. Namely, in this formulation two different types of the energy change law dE/dt and the corresponding conservation laws are obtained, which are more general than in the usual formulation. So, under certain conditions the energy conservation law has the form E=T+U+P=const, where the last term, so-called rheonomic potential expresses the influence of the nonstationary constraints.Afterwards, a detailed analysis of the obtained results and their connection with the usual formulation of mechanics are given. It is demonstrated that so formulated energy relations are in full accordance with the corresponding ones in the usual vector formulation, when they are expressed in terms of the rheonomic potential. Finally, the obtained results are illustrated by several simple, but characteristic examples.     

Numerical studies of the fingering phenomena for the thin film equation

We present a new interpretation of the fingering phenomena of the thin liquid film layer through numerical investigations. The governing partial differential equation is h_t + (h²?h³)_x = ??·(h³?Δh), which arises in the context of thin liquid films driven by a thermal gradient with a counteracting gravitational force, where h = h(x, y, t) is the liquid film height. A robust and accurate finite difference method is developed for the thin liquid film equation. For the advection part (h²?h³)_x, we use an implicit essentially non‐oscillatory (ENO)‐type scheme and get a good stability property. For the diffusion part ??·(h³?Δh), we use an implicit Euler's method. The resulting nonlinear discrete system is solved by an efficient nonlinear multigrid method. Numerical experiments indicate that higher the film thickness, the faster the film front evolves. The concave front has higher film thickness than the convex front. Therefore, the concave front has higher speed than the convex front and this leads to the fingering phenomena. Copyright © 2010 John Wiley & Sons, Ltd.     

Some particular solutions for penny-shaped crack problem by using hypersingular integral equation or differential-integral equation

Summary A hypersingular integral equation or a differential-integral equation is used to solve the penny-shaped crack problem. It is found that, if a displacement jump (crack opening displacement COD) takes the form of (a ²−x ²−y ²)^1/2 x ^m y ⁿ , where a denotes the radius of the circular region, the relevant traction applied on the crack face can be evaluated in a closed form, and the stress intensity factor can be derived immediately. Finally, some particular solutions of the penny-shaped crack problem are presented in this paper. Received 1 July 1997; accepted for publication 13 October 1997     

Mixed boundary value problem of heat conduction for infinite slab

Summary An axisymmetric steady state heat conduction boundary value problem having mixed boundary conditions on both faces of an infinite slab, is reduced to a pair of Fredholm integral equations of the second kind. For large values of h, the slab thickness, a solution correct to O(h ^–6) is obtained by expanding the kernels in power series.Presently at Imperial College, London.     

Wideband fast multipole boundary element method: Application to acoustic scattering from aerodynamic bodies

A wideband adaptive multi‐level fast multipole method (MLFMM) is used to accelerate the matrix–vector products arising from a boundary element method (BEM) formulation which solves the Burton–Miller boundary integral equation (BIE). The wideband MLFMM presented here applies a plane wave expansion formulation with fast interpolation and filtering for calculations in the high‐frequency regime and a partial wave expansion formulation with rotation‐coaxial translation in the low‐frequency regime. The iterative solvers GMRES, Bi‐CGSTAB and CGS are tested and compared and a block diagonal preconditioner is used to improve the condition number of the BEM matrices and to accelerate the convergence of the iterative solvers. Details on the implementation of the formulations are described, including the treatment of singular integrals. Results for acoustic scattering from a wing plus engine nacelle configuration for a prescribed source in a subsonic uniform flow are presented for a broad range of frequencies in order to assess the implemented capability. Copyright © 2010 John Wiley & Sons, Ltd.     

Dynamics of inner mappings

We study properties of invariant sets of dynamical systems generated by inner mappings. We prove that if x is a nonwandering point of a finitely multiple inner mapping, then not only its positive trajectory O⁺(x) consists of nonwandering points but also the negative trajectory O⁻(x) contains at least one partial semitrajectory consisting of nonwandering points.