A High-Order Kernel-Free Boundary Integral Method for Incompressible Flow Equations in Two Space Dimensions

This paper presents a fourth-order kernel-free boundary integral method for the time-dependent, incompressible Stokes and Navier-Stokes equations defined on irregular bounded domains. By the stream function-vorticity formulation, the incompressible flow equations are interpreted as vorticity evolution equations. Time discretization methods for the evolution equations lead to a modified Helmholtz equation for the vorticity, or alternatively, a modified biharmonic equation for the stream function with two clamped boundary conditions. The resulting fourth-order elliptic boundary value problem is solved by a fourth-order kernel-free boundary integral method, with which integrals in the reformulated boundary integral equation are evaluated by solving corresponding equivalent interface problems, regardless of the exact expression of the involved Green's function. To solve the unsteady Stokes equations, a four-stage composite backward differential formula of the same order accuracy is employed for time integration. For the Navier-Stokes equations, a three-stage third-order semi-implicit Runge-Kutta method is utilized to guarantee the global numerical solution has at least third-order convergence rate. Numerical results for the unsteady Stokes equations and the Navier-Stokes equations are presented to validate efficiency and accuracy of the proposed method.     

Polarized light transfer above a reflecting surface

The existence and uniqueness are established for the solution of the equation of transfer of polarized light in a homogeneous atmosphere of finite optical thickness, assuming reflection by the planetary surface. A general L_p-space formulation is adopted. The boundary value problem is first written as a vector-valued integral equation. Using monotonicity properties of the spectral radii of the integral operators involved as well as recent half-range completeness results for kinetic equations with reflective boundary conditions, the present results follow as a corollary.     

Adapted kussmaul formulation of acoustic scattering problems: an overview

Problems of exterior acoustic scattering may be conveniently formulated by means of boundary integral equations. The problem seeks to find a wave function which gives velocity potential profile, pressure density profile, etc. of the acoustic wave at points in space. At the background of the formulations are two theories viz. (Helmholtz) Potential theory and the Green's representation formula. Potential theory gives rise to the so-called indirect formulation and the Green's representation formula to the direct formulations. Classical boundary integral formulations fail at the eigenfrequencies of the interior domain. That is, if a solution is sought of the exterior problem by first solving a homogeneous boundary integral equation, one is inevitably led to the conclusion that these homogeneous boundary equations have nontrivial solutions at certain wave-numbers which are the eigenvalues of the corresponding interior problem. At lower wave-numbers, these eigenfrequencies are thinly distributed but the higher the wave-number, the denser it becomes. This is a well-known drawback for both time-harmonic acoustics and elastodynamics. This is not a physical difficulty but arises entirely as a result of a deficiency in the integral equation is representation. Why then use It? The use has many advantages notably in that the meshing region is reduced from the infinite domain exterior to the body to its finite surface. This created the need for some robust formulations. A proof of the Kussmaul [1] formulation is presented. The formulation has a hypersingular kernel in the integral operator, which creates a havoc in computation (e.g., ill conditioning). The hyper-singularity can be avoided [2], as a result a new formulation is proposed. This paper presents a broad overview of the Adapted Kussmaul Formulation (AKF).     

On the variational formulation of a transmission problem for the Biot equations

In this article, a variational formulation for the transmission problem of the fluid–bone interaction is formulated. The formulation is based on a modified Biot system of equations for the cancellous bone together with a boundary integral equation formulation of the pressure in the water. Existence and uniqueness for the weak solution of the interaction problem are established in appropriate Sobolev spaces.     

Scalar problem of plane wave diffraction by a system of nonintersecting screens and inhomogeneous bodies

The scalar problem of plane wave diffraction by a system of bodies and infinitely thin screens is considered in a quasi-classical formulation. The solution is sought in the classical sense but is defined not in the entire space ?³ but rather everywhere except for the screen edges. The original boundary value problem for the Helmholtz equation is reduced to a system of weakly singular integral equations in the regions occupied by the bodies and on the screen surfaces. The equivalence of the integral and differential formulations is proven, and the solvability of the system in the Sobolev spaces is established. The integral equations are approximately solved by the Bubnov-Galerkin method. The convergence of the method is proved, its software implementation is described, and numerical results are presented.     

BOUNDARY ELEMENT ANALYSIS FOR A CLASS OF MULTI-DIMENSIONAL LINEAR PARABOLIC EQUATIONS

In this paper,by me as of beundary element method,we try to deal with the initial -boundary value problem for a class of linear parunolic equations,which is a linear heat conduction equation. We tresent a boundary integral equation for the solution to the problem and its variational formalation The well-posedness of the variational formulation is proved. And the error estimates for the approsutate solutions are provided. The results of this paper are more general than those of[1]     

On the construction and investigation of solutions of boundary value problems of thermoviscoelasticity

The uncoupled mixed boundary value problem of thermoviscoelasticity is considered in a quasistatic formulation. The temperature distribution is assumed nonstationary and inhomogeneous. The influence of the temperature on the viscoelastic properties of the material is taken into account by the introduction of a reduced time. The equations of state of the material are written in differential form as a system of kinetic equations in some tensor-type strain parameters. The system mentioned is equivalent to a Volterra integral equation with kernel in the form of a sum of exponents. The differential approach used is apparently more convenient for numerical realization /1/ (especially in nonuniform problems) and results in a substantially different mathematical formulation as compared with that based on the integral form of writing the equations of state investigated in /2,3/. Precisely for going over to the boundary value problem are the kinetic differential equations converted into an operator differential equation in Hubert space. The existence, uniqueness, and stability of the solution of the problem formulated are established, and conditions for the convergence of the Galerkin approximations and the stability of the difference approximations in time are formulated.     

Integrability and Self-Similarity in Transient Stimulated Raman Scattering

Summary. The phenomenon of stimulated Raman scattering (SRS) can be described by three coupled PDEs which define the pump electric field, the Stokes electric field, and the material excitation as functions of distance and time. In the transient limit these equations are integrable, i.e., they admit a Lax pair formulation. Here we study this transient limit. The relevant physical problem can be formulated as an initial-boundary value (IBV) problem where both independent variables are on a finite domain. A general method for solving IBV problems for integrable equations has been introduced recently. Using this method we show that the solution of the equations describing the transient SRS can be obtained by solving a certain linear integral equation. It is interesting that this equation is identical to the linear integral equation characterizing the solution of an IBV problem of the sine-Gordon equation in light-cone coordinates. This integral equation can be solved uniquely in terms of the values of the pump and Stokes fields at the entry of the Raman cell. The asymptotic analysis of this solution reveals that the long-distance behavior of the system is dominated by the underlying self-similar solution which satisfies a particular case of the third Painlevé transcendent. This result is consistent with both numerical simulations and experimental observations. We also discuss briefly the effect of frequency mismatch between the pump and the Stokes electric fields. Received December 10, 1996; second revision received October 10, 1997; final revision received January 20, 1998     

Contact of shells with a sharp linear stamp of finite length

The action of a plane, absolutely rigid stamp on a transversely isotropic shell is investigated. The use of the equations of shells with finite shear stiffness enables the correct formulation of the problem of the action on a shell by a stamp of fixed length. The problem is reduced to an integral equation. Applying the Fourier transform, the kernel of the integral equation is represented in the form of an expansion with respect to Chebyshev polynomials. By the representation of the solution of the integral equation in the form of a product, of a series of Chebyshev polynomials and a function that takes into account the singularities of the solution at the boundary of the contact zone, the considered problem is reduced to the solving of an infinite system of linear algebraic equations, whose coefficients have been determined by the methods of numerical integration. As an example a problem for a cylindrical shell has been solved.Translated from Teoreticheskaya i Prikladnaya Mekhanika, No. 20, pp. 59–63, 1989.     

String bowing revisited

This work explores an integral equation formulation of a generalised bowing problem incorporating air resistance and bow width. The Keller/Friedlander results are obtained as a special case. Some illustrative results are presented for a lightly damped string bowed at a point. The related integral equations are solved using Filon's method, a clever but relatively unknown algorithm for the numerical evaluation of integrals with rapidly oscillating trignometric kernels.     

Green's function-based integral approaches to linear transient boundary-value problems and their stability characteristics (I)

Two Green's function-based formulations are applied to the governing differential equation which describes unsteady heat or mass transport in an isotropic homogeneous 1-D domain. In this first part of a two series of papers, the linear form of the differential equation is addressed. The first formulation, herein denoted the quasi-steady Green element (QSGE) formulation, uses the Laplace differential operator as auxiliary equation to obtain the singular integral representation of the governing equation, while the second, denoted the transient Green element (TGE), uses the transient heat equation as auxiliary equation. The mathematical simplicity of the Green's function of the first formulation enhances the ease of solution of the integral equations and the resultant discrete equations. From the point of computational convenience, therefore, the first formulation is preferred. The stability characteristics of the two formulations are evaluated by examining how they propagate various Fourier harmonics in speed and amplitude. We found that both formulations correctly reproduce the theoretical speed of the harmonics, but fail to propagate the amplitude of the small harmonics correctly for Courant value of about unity. The QSGE formulation with difference weighting values between 0.67 and 0.75, and the TGE formulation provide optimal performance in numerical stability.     

Transient elastodynamic three-dimensional problems in cracked bodies

The general formulation of the transient elastodynamic second boundary value problem in an isotropic linear elastic body with a crack of arbitrary shape by combining the boundary integral equation method and the Laplace transform with respect to time is presented in this paper. Both finite and infinite elastic bodies are considered. A numerical solution of the transformed boundary integral equations is proposed.     

Smooth boundary integration method in a linear two-dimensional problem of incompressible fluid and solid

This paper presents a new application of a theoretical and computational method of smooth boundary integration which belongs to the methods of boundary integral equations. Smooth integration is not a method of approximation. In its final analytical form, a smooth-kernel integral equation is computerized easily and accurately.

Smooth integration is associated with a “pressure-vorticity” formulation which covers linear problems in elasticity and fluid mechanics. The solution presented herein is essentially the same as that reported in an earlier paper for regular elasticity. The constraint of incompressibility does not cause difficulties in the pressure-vorticity formulation.

The linear fluid mechanics problem formulated and solved in this paper covers Stokes' problem of a slow viscous flow, and has a wider interpretation. The translational inertia forces are incorporated in the linear problem, as in Euler's dynamic theory of inviscid flow. The centrifugal inertia forces are left for the non-linear problem. The linear problem is perceived as a step in solution of the non-linear problems.     

The coupling of boundary integral and finite element methods for the bidimensional exterior steady stokes problem

In this paper, we represent a new numerical method for solving the steady-state Stokes equations in an unbounded plane domain. The technique consists in coupling the boundary integral and the finite element methods. An artificial smooth boundary is introduced separating an interior inhomogeneous region from an exterior one. The solution in the exterior domain is represented by an integral equation over the artificial boundary. This integral equation is incorporated into a velocitypressure formulation for the interior region, and a finite element method is used to approximate the resulting variational problem. This is studied by means of an abstract framework, well adapted to the model problem, in which convergence results and optimal error estimates are derived. Computer results will be discussed in a forthcoming paper.     

The collocation method for mixed boundary value problems on domains with curved polygonal boundaries

Summary. We consider an indirect boundary integral equation formulation for the mixed Dirichlet-Neumann boundary value problem for the Laplace equation on a plane domain with a polygonal boundary. The resulting system of integral equations is solved by a collocation method which uses a mesh grading transformation and a cosine approximating space. The mesh grading transformation method yields fast convergence of the collocation solution by smoothing the singularities of the exact solution. A complete stability and solvability analysis of the transformed integral equations is given by use of a Mellin transform technique, in a setting in which each arc of the polygon has associated with it a periodic Sobolev space. Received April 15, 1995 / Revised version received April 10, 1996     

Diffraction d'une Onde Acoustique par une Paroi Absorbante: Nouvelles Equations Intégrales

We present here a new formulation of the exterior Robin problem of the reduced wave equation, using boundary integrals with the jumps of the function and its normal derivative as variables. Our system of integral equations present an important advantage for any numerical purpose: it can be transformed in a variational problem where the sesqui-linear form is coercive and is exempt of hyper singular integrals.     

On a Class of Coupled Variational Formulations for a Nonlinear Parabolic Equation

51.Introducti0nSince198O)stheoriesandapplicationsofboundaryelementmethods(BEM)orboundaryintegralmethods(BIM)havemadegreatsuccessesfortheparaboliclnit1alboundaryvalueproblems(seeL1-12j),andtheapproachhasbeenappliedtonumericalsolutionsofinitialboundaryva1ueproblemssuccessfully(seeL1-5j'L8j).Thepropertiesofboundaryelementoperatorshavebeenstudiedbyboundaryintegralmethodsbymanyauthors(see.[4j,L6J'[7j'L12J).Theseresultsprovideabasisforconvergencesanderrorestimatesfornumericalapproximationofbou…     

Spectral Approach to D‐bar Problems

We present the first numerical approach to D‐bar problems having spectral convergence for real analytic, rapidly decreasing potentials. The proposed method starts from a formulation of the problem in terms of an integral equation that is numerically solved with Fourier techniques. The singular integrand is regularized analytically. The resulting integral equation is approximated via a discrete system that is solved with Krylov methods. As an example, the D‐bar problem for the Davey‐Stewartson II equations is considered. The result is used to test direct numerical solutions of the PDE.© 2017 Wiley Periodicals, Inc.     

非对称载荷作用的外部圆形裂纹问题

使用边界积分方程方法,研究了三维无限弹性体中受非对称载荷作用的外部圆形裂纹问题。通过使用Fourier级数和超几何函数,将问题的二维边界奇异积分方程简化为Abel型方程,获得了一般非对称载荷作用的外部圆形裂纹问题的应力强度因子精确解,比用Hankel变换法得到的结果更为一般。结果表明:边界积分方程法在解析分析方面还有很大的潜力。     

Second-kind integral formulations of the capacitance problem

The standard approach to calculating electrostatic forces and capacitances involves solving a surface integral equation of the first kind. However, discretizations of this problem lead to ill-conditioned linear systems and second-kind integral equations usually solve for the dipole density, which can not be directly related to electrostatic forces. This paper describes a second-kind equation for the monopole or charge density and investigates different discretization schemes for this integral formulation. Numerical experiments, using multipole accelerated matrix–vector multiplications, demonstrate the efficiency and accuracy of the new approach. This revised version was published online in June 2006 with corrections to the Cover Date.